Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling 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,729 records · Page 96

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↗

Development of Models to Simulate Tracer Tests for Characterization of Enhanced Geothermal Systems

A recent report found that power and heat produced from enhanced (or engineered) geothermal systems (EGSs) could have a major impact on the U.S energy production capability while having a minimal impact on the environment. EGS resources differ from high-grade hydrothermal resources in that they lack sufficient temperature distribution, permeability/porosity, fluid saturation, or recharge of reservoir fluids. Therefore, quantitative characterization of temperature distributions and the surface area available for heat transfer in EGS is necessary for the design and commercial development of the geothermal energy of a potential EGS site. The goal of this project is to provide integrated tracer and tracer interpretation tools to facilitate this characterization. This project was initially focused on tracer development with the application of perfluorinated tracer (PFT) compounds, non-reactive tracers used in numerous applications from atmospheric transport to underground leak detection, to geothermal systems, and evaluation of encapsulated PFTs that would release tracers at targeted reservoir temperatures. After the 2011 midyear review and subsequent discussions with the U.S. Department of Energy Geothermal Technology Program (GTP), emphasis was shifted to interpretive tool development, testing, and validation. Subsurface modeling capabilities are an important component of this project for both the design of suitable tracers and the interpretation of data from in situ tracer tests, be they single- or multi-well tests. The purpose of this report is to describe the results of the tracer and model development for simulating and conducting tracer tests for characterizing EGS parameters.

EGS↗

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↗

One- and Two-Equation Models to Simulate Ion Transport in Charged Porous Electrodes

Energy storage in porous capacitor materials, capacitive deionization (CDI) for water desalination, capacitive energy generation, geophysical applications, and removal of heavy ions from wastewater streams are some examples of processes where understanding of ionic transport processes in charged porous media is very important. In this work, one- and two-equation models are derived to simulate ionic transport processes in heterogeneous porous media comprising two different pore sizes. It is based on a theory for capacitive charging by ideally polarizable porous electrodes without Faradaic reactions or specific adsorption of ions. A two-step volume averaging technique is used to derive the averaged transport equations for multi-ionic systems without any further assumptions, such as thin electrical double layers or Donnan equilibrium. A comparison between both models is presented. The effective transport parameters for isotropic porous media are calculated by solving the corresponding closure problems. An approximate analytical procedure is proposed to solve the closure problems. Numerical and theoretical calculations show that the approximate analytical procedure yields adequate solutions. Lastly, a theoretical analysis shows that the value of interphase pseudo-transport coefficients determines which model to use.

25 ENERGY STORAGE↗

Evaluation of Aqua MODIS Thermal Emissive Bands Stability through Radiative Transfer Modeling

Moderate Resolution Imaging Spectroradiometer (MODIS) on Aqua has been inoperation providing continuous global observations for science research and applications since2002. The long-term stability of thermal emissive bands (TEBs) of Aqua MODIS was monitoredthrough inter-comparisons with measurements by hyperspectral or multi-spectral infraredsensors or through vicarious monitoring over cold targets such as Dome-C and deep convectiveclouds. The radiative transfer modeling (RTM)-based simulation models are developed to per-form the long-term monitoring of the stability of Aqua MODIS TEBs. Long-term EuropeanCentre for Medium-Range Weather Forecasts global atmospheric reanalysis data are used asinputs to the RTM simulation. By confining ocean in low to middle latitudes as the area ofinterest, the long-term stabilities of Aqua MODIS TEBs are monitored through observationminus background (O-B) brightness temperature (BT) bias between MODIS measurements andsimulations from the RTM. In general, the RTM-based O-B analysis shows that the long-termstability of Aqua MODIS TEBs are maintained well. It is shown that Aqua MODIS surfacebands are all radiometrically stable with yearly BT bias drifts <0.005 K∕yearfor B20, B22,and B23 and∼0.01 K∕yearfor B31 and B32. The carbon dioxide (CO2) absorption channelsof Aqua MODIS, e.g., B33 to B36, are stable with a BT bias drift <0.005 K∕year. It is also foundthat after accounting for the long-term growth of greenhouse gasN2O, the O-B bias trends forB24 and B25 are quite stable with yearly BT bias drifts around 0.0002 and0.0055 K∕year,respectively. The stabilities of the moisture-sensitive channels B27 to B29 and the ozone channelB30 of Aqua MODIS are also evaluated and the remnant small variations in the O-B bias trend-ing of these channels are further discussed. The RTM-based inter-comparison of MODIS TEBmeasurements with global atmospheric and climate re-analysis data provides validation of thelong-term radiometric stability of Aqua MODIS TEBs. The analysis in this paper also demon-strates that the comparison of well-calibrated MODIS TEB measurements with RTM simulationhelps quantify the long-term impacts of global greenhouse gas concentration growth on the BTdecreases in the MODIS greenhouse gas-sensitive channels.

MODIS↗

Refined Transistor Model For Simulation Of SEU

Equivalent base resistance added. Theoretical study develops equations for parameters of Gummel-Poon model of bipolar junction transistor: includes saturation current, amplification factors, charging times, knee currents, capacitances, and resistances. Portion of study concerned with base region goes beyond Gummel-Poon analysis to provide more complete understanding of transistor behavior. Extended theory useful in simulation of single-event upset (SEU) caused in logic circuits by cosmic rays or other ionizing radiation.

Zoutendyk, John A.↗

Discrete Fracture Network Modeling and Simulation of Subsurface Transport for the Topopah Springs and Lava Flow Aquifers at Pahute Mesa, FY 15 Progress Report

This progress report for fiscal year 2015 (FY15) describes the development of discrete fracture network (DFN) models for Pahute Mesa. DFN models will be used to upscale parameters for simulations of subsurface flow and transport in fractured media in Pahute Mesa. The research focuses on modeling of groundwater flow and contaminant transport using DFNs generated according to fracture characteristics observed in the Topopah Spring Aquifer (TSA) and the Lava Flow Aquifer (LFA). This work will improve the representation of radionuclide transport processes in large-scale, regulatory-focused models with a view to reduce pessimistic bounding approximations and provide more realistic contaminant boundary calculations that can be used to describe the future extent of contaminated groundwater. Our goal is to refine a modeling approach that can translate parameters to larger-scale models that account for local-scale flow and transport processes, which tend to attenuate migration.

54 ENVIRONMENTAL SCIENCES↗

Automated model generation and parameter estimation of building energy models using an ontology-based framework

This study presents a methodology for automated model generation and parameter estimation of building energy models using semantic modeling and Bayesian estimation. Semantic modeling techniques are used to represent the system components and their interactions, facilitating the automatic generation of a simulation model from dynamic component models. The proposed approach is applied to a case study of a ventilation system where a simulation model is generated, calibrated, and assessed through different performance metrics. These metrics demonstrate the accuracy and reliability of both model point estimates and probabilistic prediction intervals across all model outputs. Overall, the proposed methodology offers a systematic and automated approach to model development and calibration in building energy systems, with potential applications in building performance analysis, monitoring, and optimization.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a three-dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of solar system and exoplanetary terrestrial planets. Its parent model, known as ModelE2, is used to simulate modern Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions, including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (from slower to more rapid than modern Earth’s, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the solar system such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn’s moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents, we can then further expand its capabilities to those exoplanetary rocky worlds that have been discovered in the past, as well as those to be discovered in the future. We also discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) 1.0: A General Circulation Model for Simulating the Climates of Rocky Planets

Resolving Orbital and Climate Keys of Earth and Extraterrestrial Environments with Dynamics (ROCKE-3D) is a three-dimensional General Circulation Model (GCM) developed at the NASA Goddard Institute for Space Studies for the modeling of atmospheres of solar system and exoplanetary terrestrial planets. Its parent model, known as ModelE2, is used to simulate modern Earth and near-term paleo-Earth climates. ROCKE-3D is an ongoing effort to expand the capabilities of ModelE2 to handle a broader range of atmospheric conditions, including higher and lower atmospheric pressures, more diverse chemistries and compositions, larger and smaller planet radii and gravity, different rotation rates (from slower to more rapid than modern Earth's, including synchronous rotation), diverse ocean and land distributions and topographies, and potential basic biosphere functions. The first aim of ROCKE-3D is to model planetary atmospheres on terrestrial worlds within the solar system such as paleo-Earth, modern and paleo-Mars, paleo-Venus, and Saturn's moon Titan. By validating the model for a broad range of temperatures, pressures, and atmospheric constituents, we can then further expand its capabilities to those exoplanetary rocky worlds that have been discovered in the past, as well as those to be discovered in the future. We also discuss the current and near-future capabilities of ROCKE-3D as a community model for studying planetary and exoplanetary atmospheres.

Astrobiology↗