Search NASA⌕ Search

SEARCH · Search NASA

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Simulations and models for the Richtmyer– Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P∝ k m , where k is the wavenumber. The hydrodynamic condition is the same as in the $\theta$-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The $\theta$-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ~1.86 shock and gamma-law = 5/3 fluids with Atwood number A = 0.5. The bubble amplitude h B is found to grow in two stages. Initially, h B undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, h B grows as a power law $~t^{\theta}$ when $\langle{k}\rangle{|h_B|}$ > O(1), where $\langle{k}\rangle$ is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

GeoMicro3D: A novel numerical model for simulating the reaction process and microstructure formation of alkali-activated slag

Highlights: • The GeoMicro3D model was innovatively developed for simulating the reaction process and microstructure formation of AAS. • In GeoMicro3D, the dissolution of slag and reactions of aqueous ions are fundamentally described in the local lattice cell. • A novel approach is proposed to speed up the nucleation simulations. • GeoMicro3D was implemented and verified with the relevant experimental data and thermodynamic calculation results. For the first time, this study developed a novel model, named GeoMicro3D, to simulate the reaction process and microstructure formation of alkali-activated slag. The GeoMicro3D model consists of four modules that are designed to simulate, respectively: (i) the initial spatial distribution of real-shape slag particles in alkaline activator, (ii) the dissolution of slag and diffusion of ions via the transition state theory and lattice Boltzmann method, respectively, (iii) the spatial distribution of reaction products using a nucleation probability theory, and (iv) the chemical reactions with thermodynamic modelling. Afterwards the GeoMicro3D model was implemented and verified. The simulation results were discussed and compared with the relevant experimental data and thermodynamic calculation results using GEMS. A good agreement was found in the comparisons, showing the strong simulation capability of GeoMicro3D.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Anthropogenic emissions of NOx over China: Reconciling the difference of inverse modeling results using GOME-2 and OMI measurements

Inverse modeling using satellite observations of nitrogen dioxide (NO2) columns has been extensively used to estimate nitrogen oxides (NOx) emissions in China. Recently, the Global Ozone Monitoring Experiment-2 (GOME-2) and Ozone Monitoring Instrument (OMI) provide independent global NO2 column measurements on a nearly daily basis at around 9:30 and 13:30 local time across the equator, respectively. Anthropogenic NOx emission estimates by applying previously developed monthly inversion (MI) or daily inversion (DI) methods to these two sets of measurements show substantial differences. We improve the DI method by conducting model simulation, satellite retrieval, and inverse modeling sequentially on a daily basis. After each inversion, we update anthropogenic NOx emissions in the model simulation with the newly obtained a posteriori results. Consequently, the inversion-optimized emissions are used to compute the a priori NO2 profiles for satellite retrievals. As such, the a priori profiles used in satellite retrievals are now coupled to inverse modeling results. The improved procedure was applied to GOME-2 and OMI NO2 measurements in 2011. The new daily retrieval-inversion (DRI) method estimates an average NOx emission of 6.9 Tg N/yr over China, and the difference between using GOME-2 and OMI measurements is 0.4 Tg N/yr, which is significantly smaller than the difference of 1.3 Tg N/yr using the previous DI method. Using the more consistent DRI inversion results, we find that anthropogenic NOx emissions tend to be higher in winter and summer than spring (and possibly fall) and the weekday-to-weekend emission ratio tends to increase with NOx emission in China.

Gu, Dasa↗
DOE OSTI2016

Nonlinear electromagnetic formulation for particle-in-cell simulation of lower hybrid waves in toroidal geometry

An electromagnetic particle simulation model has been formulated and verified for nonlinear processes of lower hybrid (LH) waves in fusion plasmas. Electron dynamics are described by the drift kinetic equation using either kinetic momentum or canonical momentum. Ion dynamics are treated as the fluid system or by the Vlasov equation. Compressible magnetic perturbation is retained to simulate both the fast and slow LH waves. Numerical properties are greatly improved by using the electron continuity equation to enforce the consistency between electrostatic potential and vector potential, and by using the importance sampling scheme. The simulation model has been implemented in the gyrokinetic toroidal code (GTC), and verified for the dispersion relation and nonlinear particle trapping of the electromagnetic LH waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

RACORO continental boundary layer cloud investigations. 3. Separation of parameterization biases in single-column model CAM5 simulations of shallow cumulus

Climatically important low-level clouds are commonly misrepresented in climate models. The FAst-physics System TEstbed and Research (FASTER) project has constructed case studies from the Atmospheric Radiation Measurement (ARM) Climate Research Facility's Southern Great Plain site during the RACORO aircraft campaign to facilitate research on model representation of boundary-layer clouds. This paper focuses on using the single-column Community Atmosphere Model version 5 (SCAM5) simulations of a multi-day continental shallow cumulus case to identify specific parameterization causes of low-cloud biases. Consistent model biases among the simulations driven by a set of alternative forcings suggest that uncertainty in the forcing plays only a relatively minor role. In-depth analysis reveals that the model's shallow cumulus convection scheme tends to significantly under-produce clouds during the times when shallow cumuli exist in the observations, while the deep convective and stratiform cloud schemes significantly over-produce low-level clouds throughout the day. The links between model biases and the underlying assumptions of the shallow cumulus scheme are further diagnosed with the aid of large-eddy simulations and aircraft measurements, and by suppressing the triggering of the deep convection scheme. It is found that the weak boundary layer turbulence simulated is directly responsible for the weak cumulus activity and the simulated boundary layer stratiform clouds. Increased vertical and temporal resolutions are shown to lead to stronger boundary layer turbulence and reduction of low-cloud biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Internally cooled membrane-based absorber for dehumidification and water heating: Validated model and simulation study

This report presents a comprehensive model for advanced membrane-based absorber components accounting for three separate fluid streams, one of them providing internal cooling, with complex flow patterns. The model was implemented with various liquid desiccants and their properties (including LiCl, CaCl 2 , and [emim][OMS] ionic liquid [IL]). The model has been validated using data from two laboratory prototype absorbers, one with CaCl 2 as the working fluid and one with an IL as the working fluid. Entropy analysis was further carried out to verify the model and understand effect of internal cooling on the entropy variation of three fluid streams. The developed model and codes are expected to enable detailed configuration optimization and provide in-depth understanding of three-fluid heat and mass exchanger (HMX) performance. The paper also describes parametric studies of the HMX components in utilizing the latent heat removed in space cooling to heat water and numerically explores the unique benefit of its application in a semi-open absorption heat pump water heater. The results show that a single three-fluid HMX has the potential to achieve simultaneous dehumidification and water heating efficiently at cost effectiveness, particularly in hot and humid climates.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

Speeding up low-mass planetary microlensing simulations and modeling: The caustic region of influence

Extensive simulations of planetary microlensing are necessary both before and after a survey is conducted: before to design and optimize the survey and after to understand its detection efficiency. The major bottleneck in such computations is the computation of light curves. However, for low-mass planets, most of these computations are wasteful, as most light curves do not contain detectable planetary signatures. In this paper, I develop a parameterization of the binary microlens that is conducive to avoiding light curve computations. I empirically find analytic expressions describing the limits of the parameter space that contain the vast majority of low-mass planet detections. Through a large-scale simulation, I measure the (in)completeness of the parameterization and the speed-up it is possible to achieve. For Earth-mass planets in a wide range of orbits, it is possible to speed up simulations by a factor of ∼30-125 (depending on the survey's annual duty-cycle) at the cost of missing ∼1% of detections (which is actually a smaller loss than for the arbitrary parameter limits typically applied in microlensing simulations). The benefits of the parameterization probably outweigh the costs for planets below 100 M{sub ⊕}. For planets at the sensitivity limit of AFTA-WFIRST, simulation speed-ups of a factor ∼1000 or more are possible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Advection, moistening, and shallow‐to‐deep convection transitions during the initiation and propagation of M adden‐ J ulian O scillation

Abstract Using observations from the 2011 AMIE/DYNAMO field campaign over the Indian Ocean and a high‐resolution regional model simulation, the processes that lead to the rapid shallow‐to‐deep convection transitions associated with the initiation and eastward propagation of the Madden‐Julian Oscillation (MJO) are examined. By tracking the evolution of the depth of several thousand individual model simulated precipitation features, the role of and the processes that control the observed midtropospheric moisture buildup ahead of the detection of deep convection are quantified at large and convection scales. The frequency of shallow‐to‐deep convection transitions is found to be sensitive to this midlevel moisture and large‐scale uplift. This uplift along with the decline of large‐scale drying by equator‐ward advection causes the moisture buildup leading to the initiation of the MJO. Convection scale moisture variability and uplift, and large‐scale zonal advection play secondary roles.

Hagos, Samson↗
DOE OSTI2021

Simulation and Modeling of Time-Resolved X-Ray Detector for the Saturn Accelerator

In this work, we present the technology-aided computer design (TCAD) device simulation and modeling of a silicon p-i-n diode for detecting time-dependent X-ray radiation. We show that the simulated forward and reverse breakdown current–voltage characteristics agree well with the measured data under nonradiation environment by only calibrating carrier lifetimes for the forward bias case and avalanche model critical fields for the reverse bias condition. Using the calibrated parameters and other nominal material properties, we simulated the radiation responses of the p-i-n diode and compared with experimental data when the diode was exposed to X-ray radiation at Sandia’s Saturn facility and the Idaho State University (ISU) TriMeV facility. For Saturn’s Gaussian dose-rate pulses, we show three findings from TCAD simulations. First, the simulated photocurrents are in excellent agreement with the measured data for two dose-rate pulses with peak values of 1.16×10 10 and 1.88×10 10 rad(Si)/s. Second, the simulation results of high dose-rate pulses predict increased delayed photocurrents with longer time tails in the diode electrical responses due to excess carrier generation. Third, simulated peak values of diode radiation responses versus peak dose rates at different bias conditions provide useful guidance to determine the dose-rate range that the p-i-n diode can reliably detect in experiment. For TriMeV’s non-Gaussian dose-rate pulse, our simulated diode response is in decent agreement with the measured data without further calibration. We also studied the effects of device geometry, recombination process, and dose-rate enhancement via TCAD simulations to understand the higher measured response in the time after the peak dose-rate radiation for the p-i-n diode exposed to TriMeV irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Modelica-based modeling and simulation of district cooling systems: A case study

While equation-based object-oriented modeling language Modelica can evaluate practical energy improvements for district cooling systems, few have adopted Modelica for this type of large-scale thermo-fluid system. Further, to our best knowledge, district cooling modeling studies have yet to include hydraulics in piping networks alongside plant models featuring realistic mechanical systems and controls. These are critical details to include when looking to make energy and control improvements in many physical system installations. To fill these gaps, this study released new open-source district cooling models at the Modelica Buildings Library and leveraged these models for a real-world case study at the University of Colorado Boulder. Here, the site includes six buildings connected to a central chiller plant featuring a waterside economizer. Several energy saving strategies are pursued based on the validated model, including control setpoint optimization, equipment modification, and pump setpoint adjustments. Results indicate that a combination of the studied measures can save the campus annually 84.6 MWh of energy, 8.9% of electricity costs, 58.0 metric tons of carbon dioxide emissions, while the waterside economizer cuts down chillers’ run times by 201 days/year, reducing maintenance costs and extending chiller life.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗
DOE OSTI2015

Modeling and Simulations for the High Flux Isotope Reactor Cycle 400

A concerted effort over the past few years has been focused on enhancing the core model for the High Flux Isotope Reactor (HFIR), as part of a comprehensive study for HFIR conversion from high-enriched uranium (HEU) to low-enriched uranium (LEU) fuel. At this time, the core model used to perform analyses in support of HFIR operation is an MCNP model for the beginning of Cycle 400, which was documented in detail in a 2005 technical report. A HFIR core depletion model that is based on current state-of-the-art methods and nuclear data was needed to serve as reference for the design of an LEU fuel for HFIR. The recent enhancements in modeling and simulations for HFIR that are discussed in the present report include: (1) revision of the 2005 MCNP model for the beginning of Cycle 400 to improve the modeling data and assumptions as necessary based on appropriate primary reference sources HFIR drawings and reports; (2) improvement of the fuel region model, including an explicit representation for the involute fuel plate geometry that is characteristic to HFIR fuel; and (3) revision of the Monte Carlo-based depletion model for HFIR in use since 2009 but never documented in detail, with the development of a new depletion model for the HFIR explicit fuel plate representation. The new HFIR models for Cycle 400 are used to determine various metrics of relevance to reactor performance and safety assessments. The calculated metrics are compared, where possible, with measurement data from preconstruction critical experiments at HFIR, data included in the current HFIR safety analysis report, and/or data from previous calculations performed with different methods or codes. The results of the analyses show that the models presented in this report provide a robust and reliable basis for HFIR analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

A coupled DEM-IMB-LBM model for simulating methane hydrate exploitation involving particle dissolution

The coupled discrete element and lattice Boltzmann method using an immersed moving boundary scheme was extended to simulate methane hydrate exploitation involving mass transport and particle dissolution. In this coupled DEM-IMB-LBM model, a new Dirichlet-type thermal boundary condition is extended to simulate moving curved boundaries with constant concentration. A novel periodic boundary including an efficient searching algorithm for particle contact is proposed to reduce the computational cost and boundary effect. So, this model is validated by two numerical examples: a circular particle with concentration convection-diffusion moving in a horizontal channel and mass transport from a cylinder particle in a simple shear flow. The numerical results obtained from the proposed model agree well with previous studies. To further demonstrate the capacity of the proposed model, simulations of methane hydrate exploitation including two formations in marine sediments are carried out. The numerical results indicate that the coupled DEM-IMB-LBM is not only capable of simulating the dissolution of hydrate particles at the grain level, but also recover the sand erosion and migration process in a fundamental perspective during the methane hydrate exploitation process.

42 ENGINEERING↗
DOE OSTI2020

EGS Collab Modeling and Simulation Working Group Teleconference Series (1-98)

This submission contains the presentation slides and recordings from the first 98 EGS Collab Modeling and Simulation Working Group teleconferences. These teleconferences served three objectives for the project: 1) share simulation results, 2) communicate field activities and results to the simulation teams, and 3) hold open scientific discussions on EGS topics.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Status of Research Activities on LVDTs: Modeling and Simulation Development and Performance Characterization of Commercial LVDTs

FY-22 research focused on two major areas. The first area was modeling and simulation with linked testing that will advance the Linear Variable Differential Transformers (LVDT) technology for use by stakeholders requiring LVDTs in upcoming irradiation tests. The second area focused on testing LVDTs purchased from U.S. suppliers that were investigated in FY-21. These cross-cutting development activities will ensure stakeholders have the current state-of-the-art LVDT-based technologies for deployment in future irradiation tests. Two stand-alone reports are included in the appendixes. The purpose of this report is to provide a convenient body of work to access the reports and to evaluate the data and conclusions provided. Appendix A, “Testing of the RDP Translation Transducer,” contains a complete evaluation of an RDP Electrosense linear transducer that has the capability of replicating the form, fit, and function of LVDTs currently being used at Idaho National Laboratory (INL). Appendix B is an extract from the report, “Advanced Sensors and Instrumentation Program at Boise State University (BSU), August 2022, Linear Variable Differential Transformers.” This report contains modeling and testing results conducted at BSU and INL.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗