Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modelling and simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95

Virtual Airspace Modeling and Simulation (VAMS) Project First Technical Interchange Meeting

A three-day NASA Virtual Airspace and Modeling Project (VAMS) Technical Interchange Meeting (TIM) was held at the NASA Ames Research Center in Mountain View, CA. on May 21 through May 23,2002. The purpose of this meeting was to share initial concept information sponsored by the VAMS Project. An overall goal of the VAMS Project is to develop validated, blended, robust and transition-able air transportation system concepts over the next five years that will achieve NASA's long-term Enterprise Aviation Capacity goals. This document describes the presentations at the TIM, their related questions and answers, and presents the TIM recommendations.

Robert Beard↗

Arctic amplification is caused by sea-ice loss under increasing CO 2

Warming in the Arctic has been much faster than the rest of the world in both observations and model simulations, a phenomenon known as the Arctic amplification (AA) whose cause is still under debate. By analyzing data and model simulations, here we show that large AA occurs only from October to April and only over areas with significant sea-ice loss. AA largely disappears when Arctic sea ice is fixed or melts away. Periods with larger AA are associated with larger sea-ice loss, and models with bigger sea-ice loss produce larger AA. Increased outgoing longwave radiation and heat fluxes from the newly opened waters cause AA, whereas all other processes can only indirectly contribute to AA by melting sea-ice. We conclude that sea-ice loss is necessary for the existence of large AA and that models need to simulate Arctic sea ice realistically in order to correctly simulate Arctic warming under increasing CO 2 .

54 ENVIRONMENTAL SCIENCES↗

ROCKE-3D 2.0: An Updated General Circulation Model for Simulating the Climates of Rocky Planets

We present the second generation of ROCKE-3D (Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics), the generalized 3-dimensional General Circulation Model (GCM) for use in Solar System and exoplanetary simulations of rocky planet climates. ROCKE-3D version 2.0 is a descendant of GISS ModelE2.1, the flagship GCM of the NASA Goddard Institute for Space Studies (GISS) used in the most recent Intergovernmental Panel for Climate Change (IPCC) assessment. ROCKE-3D is a continuous effort to expand the capabilities of GISS ModelE to handle a broader range of atmospheric conditions, including different atmospheric planet sizes, gravities, pressures, rotation rates, more diverse chemistries and atmospheric compositions, diverse ocean and land distributions and topographies, and potential basic biosphere functions. We present in this release many more supported configurations which can serve as starting points to simulate terrestrial planets of interest. Two different radiation schemes are supported, the GISS radiation, valid for Earth-like planets, and SOCRATES, which is more generalized but requires more CPU resources. In addition, three different atmospheric compositions are supported, three ocean configurations (prescribed, q-flux, and dynamic), and two resolutions: the medium resolution (4 x 5 degrees in latitude x longitude) used in ROCKE-3D 1.0, and the fine resolution, which has double the resolution in the atmosphere and 4 times the resolution in the ocean. We quantify how the different component choices affect model results, and discuss strengths and limitations of using each component, together with how one can select which component to use. ROCKE-3D is publicly available and tutorial sessions are available for the community, greatly facilitating its use by any interested group.

ROCKE-3D↗

Improved Models for Simulation of Cryogenic Spray Cooling for Propellant Tank Chilldown

This paper presents advances in spray boiling models for cryogenic fluids that exhibit strong variation in thermodynamic properties; it accounts for the interaction of the spray droplets with walls and the complex heat transfer and boiling processes that play a vital role in chilling solid surfaces. The model has been developed by corroborating and validating it with test data in liquid nitrogen obtained at the University of Connecticut. The test data includes detailed visualization and plate temperature measurements for liquid nitrogen sprays impinging on a warm test plate and subsequent transient chilldown of the plate. Simulations of transient plate chilldown in liquid nitrogen were carried out and temperature histories at various locations on the plate were compared with test measurements. The final paper will present sensitivity studies and additional comparisons with test data.

Computational Fluid Dynamics, Cryogenic, Spray Chi↗

Modeling and Simulation of Electrostatics of Ge$_{\text{1-x}}$Sn$_{\text{x}}$ Layers Grown on Ge Substrates

This work introduces a comprehensive simulation tool that provides a robust 1D Schrödinger – Poisson solver for modeling the electrostatics of heterostructures with an arbitrary number of layers, and non-uniform doping profiles along with the treatment of partial ionization of dopants at low temperatures. The effective masses are derived from the first-principles calculations. The solver is used to characterize three Ge 1-x Sn x /Ge heterostructures with non-uniform doping profiles and determine the subband structure at various temperatures. Here, the simulation results of the sheet carrier densities show excellent agreement with the experimentally extracted data, thus demonstrating the capabilities of the solver.

42 ENGINEERING↗

Simulation and Prediction of Warm Season Drought in North America

This presentation presents our recent work on model simulation and prediction of warm season drought in North America. The emphasis will be on the contribution from the leading modes of subseasonal atmospheric circulation variability, which are often present in the form of stationary Rossby waves. Here we take advantage of the results from observations, reanalyses, and simulations and reforecasts performed using the NASA Goddard Earth Observing System (GEOS-5) atmospheric and coupled General Circulation Model (GCM). Our results show that stationary Rossby waves play a key role in Northern Hemisphere (NH) atmospheric circulation and surface meteorology variability on subseasonal timescales. In particular, such waves have been crucial to the development of recent short-term warm season heat waves and droughts over North America (e.g. the 1988, 1998, and 2012 summer droughts) and northern Eurasia (e.g., the 2003 summer heat wave over Europe and the 2010 summer drought and heat wave over Russia). Through an investigation of the physical processes by which these waves lead to the development of warm season drought in North America, it is further found that these waves can serve as a potential source of drought predictability. In order to properly represent their effect and exploit this source of predictability, a model needs to correctly simulate the Northern Hemisphere (NH) mean jet streams and be able to predict the sources of these waves. Given the NASA GEOS-5 AGCM deficiency in simulating the NH jet streams and tropical convection during boreal summer, an approach has been developed to artificially remove much of model mean biases, which leads to considerable improvement in model simulation and prediction of stationary Rossby waves and drought development in North America. Our study points to the need to identify key model biases that limit model simulation and prediction of regional climate extremes, and diagnose the origin of these biases so as to inform modeling group for model improvement.

Rossby Waves↗

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau-to-River (P2R) Model is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by Central Plateau Cleanup Company (CPCCo) at the Hanford Site in Washington State. Figure 1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the composite analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3 documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Buffer-IPyC separation process in TRISO fuel particles simulated with Bison code

During High Temperature Gas-cooled Reactor (HTGR) operation, tristructural isotropic (TRISO) coated-particle fuel undergoes irradiation-induced changes in morphology and thermomechanical properties. Experimental results from the Advanced Gas Reactor (AGR) Fuel Development and Qualification Program show, among other things, the mechanism of gap formation between the buffer and inner pyrolytic carbon (IPyC) layers, which could be explored further via computational simulations using the Bison code. Two simulation models were developed, the debonding restricted model, where no gap formation between buffer and IPyC layers is permitted, and the debonding enabled model, where the gap between those layers is created. The inputs of the simulated models are based on the irradiation conditions from the AGR-1 experiment. The research included simulations on spherical and aspherical fuel types. Under the specific temperatures and fluences of the AGR-1 irradiation experiment, and the Bison simulations, it was concluded that the most common scenario is a gap formation along the buffer-IPyC interface, while the least possible scenario is the situation where there is no gap formation at the buffer-IPyC junction. The computational results confirmed that the sphericity of the fuel influences the thickness of the gap that occurs at the buffer-IPyC junction, in a way that with increasing aspect ratio the gap thickness increases. The results obtained for spherical and aspherical fuel are nearly identical. Finally, performed simulations match conclusions observed from the AGR-1 experiment, which as such shows that the Bison code is a good computational method for simulating the TRISO fuel. Future simulations will include the validation of performed research and comparison of the results between Bison and PARFUME codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A conservative, interface-resolved, compressible framework for the modeling and simulation of liquid/gas phase change

This paper presents a method for simulating evaporation in a compressible, interface-resolved framework appropriate for modeling problems of engineering interest. In order to achieve robustness and broad applicability, the method has been designed to discretely enforce consistent mass and thermal energy transport at the phase interface, to globally conserve mass, momentum, and energy, and to be capable of modeling compressible and incompressible systems. Verification is performed via the Sod-shock test, one-dimensional heat conduction, evaporation from a planar interface, and evaporation of three-dimensional droplets. Convergence with increasing mesh resolution is demonstrated in all tested configurations, and conservation is maintained near machine precision for a translating droplet. Conservation and accurate phase change rates are preserved at the low numerical resolutions commonly encountered in engineering calculations. Following verification, the method is validated by comparison to an empirical correlation for evaporating droplets in high temperature crossflow, and the presentation concludes with the simulation of an iso-octane spray at conditions representative of gasoline direct injection. In conclusion, successful verification, validation, and demonstrated practical utility suggest the method to be an accurate, efficient, and robust approach for the study of phase change in engineering systems.

97 MATHEMATICS AND COMPUTING↗

Modeling and Simulation of Rotary Sloshing in Launch Vehicles

The nonlinear dynamics of propellant sloshing during orbital ascent are usually neglected in the flight control analysis of large boost vehicles under the assumption that the viscous damping of the fluid is sufficient to suppress nonlinear phenomena and confine the fluid to small, planar free surface displacements. In this case, the sloshing dynamics can be modeled using a spring-mass-damper or linearized pendulum mechanical analog. However, large, smooth-wall tanks without significant internal hardware or ring baffles are still susceptible to nonlinear effects. In particular, rotary sloshing can present a risk to flight control as it involves the formation of a stable limit cycle which can lead to undesirable roll coupling. The underlying phenomena of jump resonance does not manifest in linear models, but can be reproduced using a nonlinear spherical pendulum or the Bauer paraboloid model developed during the Apollo/Saturn program. In this paper, a detailed analysis of the rotary sloshing dynamics of these mechanical analogs is presented, and discussed in the context of flight control stability. High-fidelity simulations of a representative boost vehicle are used to verify the semi-analytical predictions of the nonlinear dynamic response.

Modeling↗

Development and testing of a mouse simulated space flight model

The development and testing of a mouse model for simulating some aspects of weightlessness that occurs during space flight, and the carrying out of immunological experiments on animals undergoing space flight is examined. The mouse model developed was an antiorthostatic, hypokinetic, hypodynamic suspension model similar to one used with rats. The study was divided into two parts. The first involved determination of which immunological parameters should be observed on animals flown during space flight or studied in the suspension model. The second involved suspending mice and determining which of those immunological parameters were altered by the suspension. Rats that were actually flown in Space Shuttle SL-3 were used to test the hypotheses.

Sonnenfeld, Gerald↗

Towards Multifluid Multiphysics Continuum Plasma Simulation for Modeling Magnetically-driven Experiments on Z

Magnetically driven experiments supporting pulsed-power utilize a wide range of configurations, including wire-arrays, gas-puffs, flyer plates, and cylindrical liners. This experimental flexibility is critical to supporting radiation effects, dynamic materials, magneto-inertial-fusion (MIF), and basic high energy density laboratory physics (HEDP) efforts. Ultimately, the rate at which these efforts progress is limited by our understanding of the complex plasma physics of these systems. Our effort has been to begin to develop an advanced algorithmic structure and a R&D code implementation for a plasma physics simulation capability based on the five-moment multi-fluid / full-Maxwell plasma model. This model can be used for inclusion of multiple fluid species (e.g., electrons, multiple charge state ions, and neutrals) and allows for generalized collisional interactions between species, models for ionization/recombination, magnetized Braginskii collisional transport, dissipative effects, and can be readily extended to incorporate radiation transport physics. In the context of pulsed-power simulations this advanced model will help to allow SNL to computationally simulate the dense continuum regions of the physical load (e.g. liner implosions, flyer plates) as well as partial power-flow losses in the final gap region of the inner MITL. In this report we briefly summarize results of applying a preliminary version of this model in the context of verification type problems, and some initial magnetic implosion relevant prototype problems. The MIF relevant prototype problems include results from fully-implicit / implicit-explicit (IMEX) resistive MHD as well as full multifluid EM plasma formulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient simulation and model reformulation of two-dimensional electrochemical thermal behavior of lithium-ion batteries

Lithium-ion batteries are an important technology to facilitate efficient energy storage and enable a shift from petroleum based energy to more environmentally benign sources. Such systems can be utilized most efficiently if good understanding of performance can be achieved for a range of operating conditions. Mathematical models can be useful to predict battery behavior to allow for optimization of design and control. An analytical solution is ideally preferred to solve the equations of a mathematical model, as it eliminates the error that arises when using numerical techniques and is usually computationally cheap. An analytical solution provides insight into the behavior of the system and also explicitly shows the effects of different parameters on the behavior. However, most engineering models, including the majority of battery models, cannot be solved analytically due to non-linearities in the equations and state dependent transport and kinetic parameters. The numerical method used to solve the system of equations describing a battery operation can have a significant impact on the computational cost of the simulation. In this paper, a model reformulation of the porous electrode pseudo three dimensional (P3D) which significantly reduces the computational cost of lithium ion battery simulation, while maintaining high accuracy, is discussed. This reformulation enables the use of the P3D model into applications that would otherwise be too computationally expensive to justify its use, such as online control, optimization, and parameter estimation. Furthermore, the P3D model has proven to be robust enough to allow for the inclusion of additional physical phenomena as understanding improves. In this study, the reformulated model is used to allow for more complicated physical phenomena to be considered for study, including thermal effects.

25 ENERGY STORAGE↗

Analysis of radiant cooling system integrated with cooling tower for composite climatic conditions

Increased demand for cooling leads to consumption of a significant amount of energy by heating, ventilation, and air conditioning (HVAC) systems in buildings. Building envelope acts as a thermal barrier and plays a significant role in improving building energy efficiency. Radiant cooling systems often utilize building structure in the form of thermal storage and provide thermal comfort. The radiant cooling system has the potential of saving peak power in the buildings. In present study experimental and a simulation study is performed for chiller and cooling tower integrated chilled ceiling radiant cooling system for chilled water source and a parallel dedicated outdoor air system (DOAS) is used for providing latent load and fresh air requirements of the zone. Cooling tower based radiant cooling system is compared with a chiller-based radiant cooling system for achieving annual energy saving. Experiments are conducted for the chiller and cooling tower operated radiant cooling system. Based on experimental data, whole building simulation models of both cooling tower and chiller-based systems were calibrated. For both the simulation model of chiller and cooling tower operated radiant cooling system simulation were carried out for round the year. The annual simulation results show that the cooling tower operated radiant cooling system saves 14% energy as compared to the chilled water operated radiating cooling system.

Srivastava, Prateek↗

Lunar Lander Landing Modeling and Simulation in Simscape

The landing gear is a critical subsystem of a lander for a successful landing on the moon or a planet. During a landing event, structural loading, energy absorption and toppling stability are important factors for the landing gear development. These aspects can be studied using simulations. The landing gear comprising the primary struct, secondary structs, footpad, and joints, is easily modeled using Simscape Multibody toolbox. The energy absorption mechanisms in the structs is readily modeled using Simulink and Multibody blocks based on the relative motions of the inner and outer cylinder of the structures. During touchdown, the footpads strike the lunar surface, and the contact forces can be modeled including the soil mechanical properties. Slosh dynamics are included using a mechanical pendulum model. Tip-over controllers are developed using Simulink blocks and toolboxes. This model was intensively tested manually then additionally with Monte Carlo simulation. It offers the opportunity for landing stability study, i.e. the max slope a lander can land without tipping over in 3D realistic landing situations.

Wei Wu↗

Computational Fluid Dynamics Modeling to Simulate a Combined Reforming Process for Syngas and Hydrogen Production

An Oxygen Transport Membrane (OTM) combined reforming technology for producing syngas and hydrogen integrates the advantages of multiple processes—steam methane reforming (SMR), autothermal reforming (ATR), an air separation unit (ASU)—into a single integrated technology. The OTM consists of a primary reforming tube, in which desulfurized natural gas is partially reformed by steam at high pressure in the presence of a metal catalyst. This process is followed in series by a ceramic OTM with a secondary reformer, in which residual methane reforms and O 2 - ions react with a portion of the CO and H 2 fuel to provide the heat to support both primary and secondary reforming. Although the OTM combined reformer technology for syngas and H 2 production has been substantially developed in the last decade, several challenges that affect the overall production efficiency and reliability are yet to be fully understood, addressed, and resolved. Therefore, developing Computational Fluid Dynamics (CFD) models that incorporate fluid dynamics, mass transport, kinetics, heat transport, and structural mechanics is critical to understanding and minimizing the probability of tube failures during the startup and operation. In this report, an exhaustive literature review was performed to survey the current state of technology for producing syngas and H 2 using either conventional or renewable energy sources. The feedstocks reviewed include natural gas and coal for the conventional technologies, whereas biomass, solar, wind, and nuclear energy for the renewable technologies. The existing industry-grade COMSOL multiphysics models of OTM were upgraded for the latest software release. In addition, they were improved to help achieve grid and solver independence and were successfully ported on the ORNL high-performance computing clusters to speed up their run times. A 42% reduction in the simulation run time was achieved. A new higher-fidelity CFD model of an OTM tube was developed in the StarCCM+ simulation platform. This new model was designed to simulate various physics using first principles, e.g., turbulent flow, heat transfer, and chemical reactions while avoiding unnecessary simplifications. The resulting predictions were qualitatively assessed and provided useful insights into the multiphysics complexity of an OTM tube.

08 HYDROGEN↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗