Search NASA⌕ Search

SEARCH · Search NASA

Results for “Advanced 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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

CTF Theory Manual: Version 4.3

Coolant-Boiling in Rod Arrays—Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and three-dimensional Cartesian forms of its governing equations are available for the solution. The code was originally developed by Pacific Northwest Laboratory in 1980, and had been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently rebranded as CTF. CTF was later adopted in the early 2010s by Oak Ridge National Laboratory (ORNL) for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, which led to a significant advancement of the code software quality, modeling accuracy, testing systems, and capabilities for improved support of modeling of common LWR nominal and transient behavior. As part of the improvement process, it was necessary to generate sufficient documentation for the public domain code which had lacked such material upon being adopted by RDFMG. This document serves as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code, as well as the numerical solution process utilized. Additional documents available in the CTF documentation suite include the user manual and verification and validation manual.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Guest Editorial Special Section on Advanced Medium-Voltage Power Electronics for Grid Interactive Applications

Medium-voltage power electronics (MVPE) plays essential roles in power grid modernization and links the MV distribution grid with low-voltage consumers and prosumers. Various MVPE devices, such as solid-state transformers or circuit breakers, inverter-based resources, power flow controllers, etc., bring the benefits of voltage conversion and power regulation in small footprint, power quality and efficiency improvements, and enhancements of grid controllability, flexibility, stability, and resilience. The MVPE also makes it possible for sustainable energy systems, such as solar/wind farms and energy storage generating facilities, to directly access to MV grids without multistage conversions. With their intrinsic intelligence and communications, MVPE enables many new smart grid functions and applications, e.g., dc interconnections and electric vehicle charging, which were not envisioned by traditional power grids otherwise. In addition, the integration of physical power processing units with cyber components forms a cyber-physical system, which is essential for long-term sustainability, development, and environmental preservation. Nonetheless, technical challenges on MVPE device reliability, scalable and efficient converter topologies, control stability, large-scale modeling and simulation, to name a few, need to be addressed and advanced to the next level. In conclusion, this Special Section on Advanced MV Power Electronics for Grid Interactive Applications in IEEE Transactions on Power Electronics (TPEL) provides an insight on some of the recent advances in MVPE and emerging challenges and potential solutions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗
DOE OSTI2017

Evaluation of Stimulation by Shear Slip in Fractured Rock Using a 3D Coupled Thermo-Poro-Mechanical FEM

This paper presents a 3D finite element model for simulating shear slip on fractures in rock in response to coupled thermo-poro-mechanical processes. The modeling is part of a broader effort to investigate the role of pore pressure and cooling by cold water injection on shear slip and permeability increase in granitic rock using laboratory shearing experiments under triaxial conditions. In particular, 3D thermo-poro-mechanical finite element modeling and analysis of injection experiment in fractured rock has been carried out to analyze the role of pore pressure, temperature, fracture deformation and their interactions. 4-noded tetrahedral elements are employed for intact rock deformation and transport processes within the matrix. To represent the mechanical response of a fracture surface, zero thickness contact interface element is developed based on recently developed element partitioning algorithm and procedures for applying hydraulic pressure on the fracture surfaces. A Mohr-Coulomb type criterion is adapted to capture the slip on the fracture and to quantify its reversible and irreversible deformation. The numerical model has been calibrated using results from well controlled, advanced laboratory experiments. Excellent agreement between modeling and experimental observations is achieved. Simulation results illustrate that pore pressure increase and rock matrix/fracture surface cooling cause the fracture system to deform and slip. Fracture slip is promoted due to its normal stress reduction associated with cooling effect of cold fluid injection. The numerical model provides a physics-based understanding of the role of coupled processes on shear stimulation phenomenon and the resulting permeability enhancement.

58 GEOSCIENCES↗
DOE OSTI2018

Two-Step Procedure for Liquid-Salt-Cooled-Reactor Analysis

Here, this paper presents the development of a lattice physics–to–core simulator two-step procedure for the rapid analysis of the Advanced High Temperature Reactor (AHTR). Lattice physics, reflector, and control blade models were developed from which cross-section libraries could be generated for a nodal core simulator. Few-group structures for the core simulator were also generated to account for the neutronic characteristics of AHTR. After developing the AHTR two-step procedure, cross-section libraries were generated using the SERPENT continuous-energy Monte Carlo code. These libraries were then used in the core simulator NESTLE to perform full-core calculations, which were in turn benchmarked against reference SERPENT full-core models. Benchmarking results showed reasonable accuracy of the developed two-step procedure but revealed an inherent inadequacy in the one-dimensional radial reflector model and showed a likely need for a greater number of energy groups than were used in this study.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Model fusion with physics-guided machine learning: Projection-based reduced-order modeling

The unprecedented amount of data generated from experiments, field observations, and large-scale numerical simulations at a wide range of spatiotemporal scales has enabled the rapid advancement of data-driven and especially deep learning models in the field of fluid mechanics. Although these methods are proven successful for many applications, there is a grand challenge of improving their generalizability. This is particularly essential when data-driven models are employed within outer-loop applications like optimization. In this work, we put forth a physics-guided machine learning (PGML) framework that leverages the interpretable physics-based model with a deep learning model. Leveraging a concatenated neural network design from multi-modal data sources, the PGML framework is capable of enhancing the generalizability of data-driven models and effectively protects against or inform about the inaccurate predictions resulting from extrapolation. We apply the PGML framework as a novel model fusion approach combining the physics-based Galerkin projection model and long- to short-term memory (LSTM) network for parametric model order reduction of fluid flows. We demonstrate the improved generalizability of the PGML framework against a purely data-driven approach through the injection of physics features into intermediate LSTM layers. Our quantitative analysis shows that the overall model uncertainty can be reduced through the PGML approach, especially for test data coming from a distribution different than the training data. Moreover, we demonstrate that our approach can be used as an inverse diagnostic tool providing a confidence score associated with models and observations. The proposed framework also allows for multi-fidelity computing by making use of low-fidelity models in the online deployment of quantified data-driven models.

42 ENGINEERING↗
DOE OSTI2025

Review of recent WEC-Sim (v6.1) advanced features

WEC-Sim (Wave Energy Converter SIMulator) is an opensource software for modeling the motions, loads and power generation of wave energy converters. WEC-Sim performs simulations in the time domain using hydrodynamic coefficients calculated by boundary element method (BEM) frequency-domain potential flow solvers such as WAMIT, NEMOH, Capytaine, or Ansys AQWA. WEC-Sim development is ongoing, including various features, applications, and example cases used to demonstrate potential use-cases to meet the needs of the growing marine energy industry. Through input from a broad user base and an extensive team of developers and collaborators, new features of WEC-Sim are developed to expand the software’s use cases and improve overall functionality. Three new WEC-Sim features highlighted in this paper include updating WEC-Sim to be compatible with MoorDyn Version 2, incorporation of second order excitation loads (quadratic transfer functions) and allowing for dynamically changing hydrodynamics.

13 HYDRO ENERGY↗
DOE OSTI2022

A Multicomponent Reactive Transport Model for Integrated Surface‐Subsurface Hydrology Problems

Abstract Despite the widespread use of integrated hydrology models in a variety of applications, consideration of multicomponent reactive transport is still not common. The implementation of these processes requires coupling transport at the surface‐subsurface interface and efficient solution of the non‐linear geochemical model that is consistent with the integrated hydrology solution. The Advanced Terrestrial Simulator provides a flexible multiphysics framework that facilitated this process. In this work, the integrated reactive transport process kernel (PK) was weakly coupled to the integrated hydrology PK. In turn, integrated transport and reactions were coupled using an operator splitting approach. This splitting enabled an explicit solution of the integrated transport problem, including a novel algorithm to calculate exchange fluxes across the surface‐subsurface interface and a point‐by‐point solution of the geochemical problem. Geochemical capabilities were added using well‐established external codes, but rather than using a custom interface to each, a generic interface was used that clearly specifies the variables and operations used by the chemistry PK. The implementation is demonstrated with two example simulations: transport of a tracer in a soil column as it saturates over time and water ponds on the surface and reactive transport in a hillslope driven by successive wet‐dry cycles that result in infiltration, runoff and exfiltration processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model: Modeling Archive

This archive contains the data and Python scripts required to reproduce the analyses and figures in the study: Gomez-Velez, J. D., Rathore, S. S., Cohen, M. J., & Painter, S. L. (2025). Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model. Submitted to Water Resources Research. The analysis utilizes the subgrid model Advection Dispersion Equation with Lagrangian Subgrids (ADELS) implemented in the Advanced Terrestrial Simulator (ATS; https://amanzi.github.io/ats/stable/). In this case, the ATS and Amanzi versions are (1) ATS version 1.5.1_f5ba18f8 and (2) Amanzi version 1.6-dev_53444cca4. The repository includes a Jupyter Notebook and the necessary data (Pandas DataFrames stored as pickle files) to generate the figures for the manuscript. Additionally, it contains Python scripts to create ATS input files, run the ATS simulations, and post-process the results. Finally, it provides routines for parameter estimation using the Single-Station Metabolism (SSM) model with the Differential Evolution Adaptive Metropolis (DREAM) Markov Chain Monte Carlo (MCMC) algorithm with ZS enhancements (DREAM-ZS).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Role of Neutrals Versus Transport in Determining the Pedestal Density Structure: Final Technical Report

In fusion devices the plasma density plays a crucial role in determining the fusion reaction rate and has a direct impact on the fusion gain of a given device. This density is in general regulated by the particle sources and transport near the plasma edge, which give rise to an edge density pedestal. When predicting the performance of future devices, this density pedestal is often prescribed, rather than predicted, due to a lack of models which allow confident extrapolation. This project aims to advance these models through the focused validation of theoretical models related to the transport of fueling neutral particles, and through interpretive transport modeling in present day fusion plasmas, in which the penetration of neutrals is altered to better simulate future reactor-like conditions. Achievements in theory and model validation under this project have advanced our understanding of how much of the edge density profile is set by transport versus direct ionization, enabling interesting projections to future burning plasma devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Proceedings of the 1998 International Computational Accelerator Physics Conference

The 1998 International Computational Accelerator Physics Conference, ICAP’98, was held September 14-18, 1998, at the Monterey Convention Center in Monterey, California. The meeting continued the tradition of four earlier Computational Accelerator Physics conferences, held in La Jolla, California (1988), Los Alamos, New Mexico (1990), Pleasanton, California (1993), and Williamsburg, Virginia (1996). ICAP’98 was organized by the Stanford Linear Accelerator Center (SLAC), Los Alamos National Laboratory (LANL), and the National Energy Research Scientific Computing Center (NERSC), and received support from the U.S. Department of Energy. The aim of this conference series is to discuss recent advances in Computational Accelerator Physics and their applications to solving important problems related to present and future accelerators. The ICAP’98 program included plenary sessions on modeling challenges in next-generation accelerators, state-of-the art methods in computational accelerator physics, and a view of what high-end computing would be like in the 21st century. There were also focused sessions on accelerator design and control, modeling high intensity beam phenomena, electromagnetic modeling, advanced algorithms, software environments and libraries, visualization and graphic interfaces, modeling advanced accelerator concepts, modeling the beam-beam interaction, computing wakefields, and combined particle/field simulations. In addition to these oral sessions there was a poster session and a vendor display. Prior to the start of the technical program, a half-day tutorial on High Performance Computing was hosted by NERSC. The ICAP series recognizes that computing and, in particular, advanced computing, is playing an increasingly important role in the successful design and development of major accelerator facilities, and in research and development leading to advances in accelerator technology. We expect that this new emphasis on high-end computing will continue in future ICAP conferences.

Ko, Kwok↗
DOE OSTI2025

Evaluating E3SM Global Storm‐Resolving Model Simulations of Deep Convection: Insights From DP‐SCREAM During TRACER

Global Storm-Resolving Models (GSRMs) are becoming increasingly vital for advancing climate modeling and improving the prediction of extreme weather events. Houston, a coastal region frequently affected by deep convective storms, offers an ideal setting to evaluate the ability of GSRMs to simulate deep convection. This study assesses the performance of the Doubly Periodic Simple Cloud-Resolving E3SM (Energy Exascale Earth System Model) Atmosphere Model (DP-SCREAM) using observations from the TRacking Aerosol Convection interactions ExpeRiment (TRACER) campaign. DP-SCREAM effectively reproduces the diurnal cycles of clouds and precipitation, demonstrating much greater skill than the E3SM single column model. The DP-SCREAM is demonstrated to be applicable to coastal regions, partially due to the forcing data sets already capturing the influence of breezes. DP-SCREAM also replicates biases persistent in the global version of SCREAM: the underrepresentation of boundary layer shallow clouds, a lack of mid-level congestus clouds, and the popcorn convection, characterized by small and disorganized convective cells generating the strongest precipitation. To investigate these issues, two sensitivity experiments were conducted: increasing the mixing length and scaling up the buoyancy flux within the Simplified Higher Order Closure scheme. Increasing the mixing length improved mid-level congestus representation and reduced unrealistic early morning fog occurrence. Enhancing buoyancy flux only marginally improved the bias of underproduced big convective cells. In conclusion, an additional resolution sensitivity test at 0.5 km grid spacing demonstrated that a refined horizontal resolution alone is insufficient to resolve these biases.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

CTF Theory Manual: Version 4.4

Coolant-Boiling in Rod Arrays – Two Fluids (COBRA-TF) is a thermal/hydraulic (T/H) simulation code designed for light-water reactor (LWR) vessel analysis. It uses a two-fluid, three-field (i.e., fluid film, fluid drops, and vapor) modeling approach. Both subchannel and 3D Cartesian forms of its governing equations are available for the solution. The code was originally developed by Pacific Northwest Laboratory in 1980 and has been used and modified by several institutions over the last few decades. COBRA-TF also found use at the Pennsylvania State University (PSU) by the Reactor Dynamics and Fuel Modeling Group (RDFMG) and has been improved, updated, and subsequently rebranded as CTF. CTF was later adopted in the early 2010s by Oak Ridge National Laboratory (ORNL) for use in the Consortium for Advanced Simulation of Light Water Reactors (CASL) program, which led to a significant advancement of the code’s software quality, modeling accuracy, testing systems, and capabilities for improved support in modeling common LWR nominal and transient behavior. As part of the improvement process, it was necessary to generate sufficient documentation for the public domain code which had lacked such material upon being adopted by RDFMG. This document serves as a theory manual for CTF, detailing the many two-phase heat transfer, drag, and important accident scenario models contained in the code, as well as the numerical solution process utilized. Additional documents available in the CTF documentation suite include the user manual and the verification and validation manual.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Airborne measurements of isoprene and monoterpene emissions from southeastern U.S. forests

Isoprene and monoterpene emission rates are essential inputs for atmospheric chemistry models that simulate atmospheric oxidant and particle distributions. Process studies of the biochemical and physiological mechanisms controlling these emissions are advancing our understanding and the accuracy of model predictions but efforts to quantify regional emissions have been limited by a lack of constraints on regional distributions of ecosystem emission capacities. We used an airborne wavelet-based eddy covariance measurement technique to characterize isoprene and monoterpene fluxes with high spatial resolution during the 2013 SAS (Southeast Atmosphere Study) in the southeastern United States. The fluxes measured by direct eddy covariance were comparable to emissions independently estimated using an indirect inverse modeling approach. Isoprene emission factors based on the aircraft wavelet flux estimates for high isoprene chemotypes (e.g., oaks) were similar to the MEGAN2.1 biogenic emission model estimates for landscapes dominated by oaks. Aircraft flux measurement estimates for landscapes with fewer isoprene emitting trees (e.g., pine plantations), were about a factor of two lower than MEGAN2.1 model estimates. The tendency for high isoprene emitters in these landscapes to occur in the shaded understory, where light dependent isoprene emissions are diminished, may explain the lower than expected emissions. This result demonstrates the importance of accurately representing the vertical profile of isoprene emitting biomass in biogenic emission models. Airborne measurement-based emission factors for high monoterpene chemotypes agreed with MEGAN2.1 in landscapes dominated by pine (high monoterpene chemotype) trees but were more than a factor of three higher than model estimates for landscapes dominated by oak (relatively low monoterpene emitting) trees. This results suggests that unaccounted processes, such as floral emissions or light dependent monoterpene emissions, or vegetation other than high monoterpene emitting trees may be an important source of monoterpene emissions in those landscapes and should be identified and included in biogenic emission models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Steepest-descent algorithm for simulating plasma-wave caustics via metaplectic geometrical optics

The design and optimization of radiofrequency-wave systems for fusion applications is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at wave cutoffs and caustics. To accurately model the wave behavior in these regions, more advanced and computationally expensive “full-wave” simulations are typically used, but this is not strictly necessary. A new generalized formulation called metaplectic geometrical optics (MGO) has been proposed that reinstates GO near caustics. The MGO framework yields an integral representation of the wavefield that must be evaluated numerically in general. We present an algorithm for computing these integrals using Gauss-Freud quadrature along the steepest-descent contours. Benchmarking is performed on the standard Airy problem, for which the exact solution is known analytically. Furthermore, the numerical MGO solution provided by the new algorithm agrees remarkably well with the exact solution and significantly improves on previously derived analytical approximations of the MGO integral.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electromagnetics and Advanced Manufacturing Simulation Development within the MOOSE Ecosystem

The MOOSE electromagnetics module (EMM) was developed to provide a general simulation tool for modeling a wide range of problems encountered in computational electromagnetics and enable MOOSE and MOOSE-based application coupling to Maxwell's equations to solve a wide class of new problems. Development of the electromagnetics module has also enabled development of new advanced manufacturing models in the MOOSE ecosystem, including that of the electric field-assisted sintering (EFAS) process. Recent results from this effort show good agreement with experimental results, and these two capabilities will be used in future projects to develop models for novel materials development and manufacturing for extreme environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Electromagnetic Transient Simulation of Large-Scale Inverter-Based Resources With High-Granularity

The power grid is undergoing a significant transformation with the rapid increase in inverter-based resources (IBRs), including large-scale photovoltaic (PV) plants. Ensuring reliable and resilient grid operation in this new paradigm necessitates high-granularity electromagnetic transient (EMT) modeling that accurately captures the behavior of individual inverters and their interactions within IBR plants. Central to this approach is the detailed representation of both the IBR plant’s collector system and the dynamics of individual inverters. To achieve this, a high-granularity EMT model of a large-scale PV plant has been developed using advanced simulation algorithms, including matrix splitting and the Schur complement. These proposed techniques significantly enhance simulation speed, numerical stability, and accuracy while improving the modularity and efficiency of the collector system’s representation. The effectiveness of the proposed methods is validated through simulations of a representative large-scale PV plant consisting of 125 individual PV inverters, 25 IBR unit transformers, and a 52-bus collector system.

Choi, Jongchan [Oak Ridge National Laboratory (ORN↗
DOE OSTI2025

First JSNS 2 measurement of the electron neutrino flux using the 12 C⁡(𝜈 𝑒 ,𝑒 − )⁢ 12 N g.s. reaction

JSNS2 (J-PARC Sterile Neutrino Search at J-PARC Spallation Neutron Source) is an experiment searching for sterile neutrinos through the observation of ν¯μ→ν¯e appearance oscillations, using neutrinos produced by muon decay-at-rest. A key aspect of the experiment involves accurately understanding the neutrino flux and the quantities of pions and muons, which are progenitors of (anti)neutrinos, given that their production rates have yet to be measured. We present the first electron-neutrino flux measurement using C12(νe,e−)12Ng.s. reaction in JSNS2, yielding a flux of (6.7±1.6(stat)±1.7(syst))×10−9 cm−2 proton−1 at the JSNS2 detector location, located at 24 meters distance from the mercury target. This flux measurement is consistent with predictions from simulations based on hadron models.

Particle decays↗