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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Long‐Term Large‐Scale Atmospheric Forcing Data From Three‐Dimensional Constrained Variational Analysis for the ARM SGP Site

Here, this study presents a long‐term three‐dimensional large‐scale forcing data set (VARANAL3D) derived from the three‐dimensional constrained variational analysis (3DCVA) method at the Atmospheric Radiation Measurement (ARM) program Southern Great Plains (SGP) site from 2004 to 2018. Building on the same input data sets as the conventional continuous forcing data set (VARANAL), VARANAL3D maintains overall consistency in domain‐averaged fields while introducing spatial variability, offering critical insights into the influence of mesoscale synoptic systems on cloud‐related processes. Evaluations are conducted across four cloud and precipitation regimes: Clear‐sky, Shallow‐clouds, Afternoon‐precipitation, and Nocturnal‐precipitation, presenting high consistency of the domain‐mean forcing data sets while emphasizing the role of subdomain forcing variability particularly in precipitating regimes. Single column model (SCM) simulations demonstrate that subdomain VARANAL3D forcing improves cloud and precipitation representation, with the ensemble outperforming domain‐mean forcing in three cloudy and precipitating regimes. Overall, these results highlight VARANAL3D's value for investigating the impacts of spatial variability of large‐scale forcing on atmospheric processes. The VARANAL3D data set provides new opportunities for evaluating model physics, advancing the development of scale‐aware parameterizations and deepening our understanding of cloud and precipitation dynamics.

Environmental sciences↗
DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗
DOE OSTI2026

Impact of high-fidelity temperature feedback modeling on fusion blanket neutronics

Accurate modeling and simulation of neutron and gamma transport and interactions within fusion tokamaks represents a critical area of research for fusion energy development. This paper presents a method to couple computational fluid dynamics (CFD) results to high-fidelity neutronics analyses, to explore the fusion neutronics’ feedback to blanket temperature and density. This method is demonstrated on a sector model of a tokamak with immersion molten salt blanket, and the relative differences in tritium breeding ratio (TBR) with spatially-resolved temperature modeling are explored. Additionally, the influence of higher fidelity neutronics results on the subsequent CFD analysis is explored. Results show that global TBR changes $<$0.1% due to temperature spatial distribution modeling. However, local TBR in areas with high temperature variability can have up to ∼4% difference.

blanket fusion↗
DOE OSTI2025

CO2 Plume Imaging with Accelerated Deep Learning-based Data Assimilation Considering Multiple Realizations: Application to the Illinois Basin-Decatur Carbon Sequestration Project

We propose a fast and efficient deep learning workflow for near real-time data assimilation, forecasting and visualization of CO2 plume evolution in saline aquifer and demonstrate its application at a field site. Unlike the previous work, this study incorporates the impact of spatial heterogeneity using multiple realizations. In the proposed workflow, a neural network model utilizes available monitoring data such as downhole pressure measurements as input and predicts the propagating pressure ‘front’ using the diffusive time of flight (DTOF) map which is considered as representative reservoir image of the flow field. The DTOF is the arrival time of pressure front propagation, which can be computed by the Fast Marching Method rapidly without flow simulations. Reservoir model calibration can be implemented by selecting the training data samples that describe the predicted DTOF map based on observed data. The power and efficacy of our workflow is demonstrated by application to the Illinois Basin-Decatur Project.

CO2 plume imaging↗
DOE OSTI2025

Studies of the influence of heat pipe model and coupling strategy in heat pipe microreactor simulations using Sockeye

This work demonstrates various modeling approaches for thermally coupling heat pipes to a graphite monolithic block in a microreactor context, accounting for the presence of a gap. We show that heat pipe model selection has a significant impact on the transient behavior of a three-dimensional coupled microreactor assembly problem by comparing results obtained using different heat pipe models from the heat pipe code Sockeye. Additionally, we show that the choice of coupling strategy has significant consequences for both accuracy and performance. Furthermore, we found that our heat flux transfer strategy exhibited greater accuracy than our temperature transfer strategy, even when comparing a loose coupling of the heat flux transfer strategy to a tight coupling of the temperature transfer strategy. Although this research uses only theoretical test problems, it provides important insights in modeling thermal fluids phenomena in heat pipe microreactors.

Direwolf↗
DOE OSTI2025

The performance of missing transverse momentum reconstruction and its significance with the ATLAS detector using 140 $\hbox {fb}^{-1}$ of $\sqrt{s}=13$ TeV pp collisions

This paper presents the reconstruction of missing transverse momentum ($p_{\text {T}}^{\text {miss}}$ ) in proton–proton collisions, at a center-of-mass energy of 13 TeV. This is a challenging task involving many detector inputs, combining fully calibrated electrons, muons, photons, hadronically decaying $\tau$ -leptons, hadronic jets, and soft activity from remaining tracks. Possible double counting of momentum is avoided by applying a signal ambiguity resolution procedure which rejects detector inputs that have already been used. Several $p_{\text {T}}^{\text {miss}}$ ‘working points’ are defined with varying stringency of selections, the tightest improving the resolution at high pile-up by up to 39% compared to the loosest. The $p_{\text {T}}^{\text {miss}}$ performance is evaluated using data and Monte Carlo simulation, with an emphasis on understanding the impact of pile-up, primarily using events consistent with leptonic Z decays. The studies use $140~\text {fb}^{-1}$ of data, collected by the ATLAS experiment at the Large Hadron Collider between 2015 and 2018. The results demonstrate that $p_{\text {T}}^{\text {miss}}$ reconstruction, and its associated significance, are well understood and reliably modelled by simulation. Finally, the systematic uncertainties on the soft $p_{\text {T}}^{\text {miss}}$ component are calculated. After various improvements the scale and resolution uncertainties are reduced by up to 76% and 51%, respectively, compared to the previous calculation at a lower luminosity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Nuclear Data Impact on Key Metrics for a Representative Molten Chloride Fast Reactor Model

Nuclear data are an essential component of the foundation on which all modeling and simulation methods and tools are relying upon, from the front end to the back end of the nuclear fuel cycle. In this study, the impact of uncertainties in nuclear data is investigated for a representative molten chloride fast reactor, for several important metrics, including eigenvalue, reactivity differences, and nuclide inventories in fuel at 5-yr irradiation. Uncertainty of keff for a full core model was found to be similar between the fresh fuel and the irradiated fuel states (1.7-1.8%), with its primary driver being the uncertainty in the 235U (n,γ) cross section. The results obtained for the reactivity differences show large uncertainties, of over 100%, in elastic scattering sensitivities of several nuclides, which led to large uncertainties of temperature reactivity differences for cladding and reflector. These results provide evidence that the currently applied methods may not be sufficiently adequate for ensuring the reliable determination of such metrics.

Procop, Germina [ORNL] (ORCID:0000000342226393)↗
DOE OSTI2025

Modeling Co2 Flow Through Faulted/Fractured Reservoirs Using Tedfm in Corner-Point Grids

Interest in underground CO2 storage has increased significantly over the last decade, driven by growing concern about global warming and rising levels of greenhouse gases in the atmosphere. Given that CO2 accounts for 80% of these greenhouse gases, carbon capture, utilization, and storage (CCUS) is considered one of the most direct approaches to achieving the net-zero carbon target. Although CO2 storage in deep saline aquifers and depleted gas reservoirs has been studied extensively, most studies use commercial simulators that model faults/fractures by simply modifying the transmissibility in the direction perpendicular to the fault surfaces. This work shows that this simplistic approach ignores the accelerated flow in the directions parallel to the fault plane, leading to significantly higher leakage along the fault surface. To accurately model CO2 flow in faulted reservoirs, we present the first transient embedded discrete-fracture model for corner-point grids (tEDFM-CPG). By comparing the tEDFM-CPG results with high-resolution reference solutions, we show that this approach is accurate and efficient at predicting CO2 flow in faulted/fractured reservoirs. In conclusion, this work presents the use of mixed reality (MR) to efficiently observe CO2 gas migration in the interior of these corner-point grid systems.

02 PETROLEUM↗
DOE OSTI2025

Griffin Capability Improvements in Support of Ex-core Deep-Penetration Problems

Advanced reactor designs, especially portable reactors that are designed to be located closer to humans and operate autonomously, require the ability to accurately compute the ex-core neutron and gamma flux solutions in terms of shielding design optimization to reduce dose rates at the vessel boundary and detector signal prediction to drive the reactor control system. The Nuclear Energy Advanced Modeling and Simulation program has prioritized improvements to the Griffin discrete ordinates (SN) solver for deep-penetration problems in fiscal year 2025. Significant advancements have been made to the Griffin methodologies for solving ex-core deep-penetration problems for steady-state, fixed-source and transient calculations. This work presents the methodology improvements as well as a comprehensive demonstration with a Transient Test Reactor model and measurements.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Modeling CO 2 flow through faulted/fractured reservoirs using tEDFM in corner-point grids

The interest in underground CO 2 storage has increased significantly over the last decade because of the rising concern about global warming due to the growing levels of greenhouse gases in the atmosphere. Considering that CO 2 accounts for 80% of these greenhouse gases, carbon capture, utilization, and storage (CCUS) is regarded as one of the most direct approaches to achieving the net zero carbon target. Although CO 2 storage in deep saline aquifers and depleted gas reservoirs has been studied extensively, most studies use commercial simulators that model faults/fractures by simply modifying the transmissibility in the direction perpendicular to the fault surfaces. Here, this work shows that this simplistic approach ignores the accelerated flow in the directions parallel to the fault plane, leading to significantly higher leakage along the fault surface. To accurately model the flow of CO 2 in faulted reservoirs, we present the first transient embedded discrete fracture model for corner-point grids (tEDFM-CPG). By comparing the results of the tEDFM-CPG to high-resolution reference solutions, we show that this approach is accurate and efficient at predicting CO 2 flow in faulted/fractured reservoirs. Finally, this work presents the use of mixed reality (MR) to efficiently observe CO 2 gas migration in the interior of these corner-point grid systems.

25 ENERGY STORAGE↗
DOE OSTI2025

Coupled Hydrological-Thermal-Biogeochemical Modeling for Predicting Arctic Carbon Emissions (CH4PACE)

The Arctic, with nearly 1,700 billion metric tons of carbon stored in its soils, is experiencing significant warming trends that raise concerns about the potential release of greenhouse gases, particularly methane. This report investigates the conversion of carbon stored in Arctic soils and permafrost into methane and its subsequent atmospheric release over the coming decades. Recent data indicate a rise in global methane concentrations, highlighting the importance of understanding the mechanisms of methane generation under anaerobic conditions prevalent in thawing permafrost. Our study employs a multi-faceted approach that integrates fieldwork, benchtop experiments, and numerical simulations to explore the subsurface hydrological-thermal-biogeochemical systems influencing methane emissions. We focus on the relationship between microbial activity in decomposing organic matter under oxygen-limited conditions and the transport of methane out of these saturated wetland soils. Key hydrological parameters, including permeability, porosity, and thermal conductivity, were measured on soil samples from Fairbanks, Alaska, and used to parameterize PFLOTRAN, a coupled hydrological-thermal-reactive transport simulator. Our model enables prediction of the evolution of the active layer and its impact on methane production and release. Our findings underscore the critical need for a mechanistic understanding of methane dynamics in permafrost environments to improve predictions of methane emissions in Earth systems models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Quantitative Modeling of High-Energy Electron Scattering in Thick Samples Using Monte Carlo Techniques

Cryo-electron microscopy (cryo-EM) is a powerful tool for imaging biological samples but is typically limited by sample thickness, which is restricted to a few hundred nanometers depending on the electron energy. However, there is a growing need for imaging techniques capable of studying biological samples up to 10 µm in thickness while maintaining nanoscale resolution. This need motivates the use of mega-electron-volt scanning transmission electron microscopy (MeV-STEM), which leverages the high penetration power of MeV electrons to generate high-resolution images of thicker samples. In this study, we employ Monte Carlo simulations to model electron–sample interactions and explore the signal decay of imaging electrons through thick specimens. By incorporating material properties, interaction cross-sections for energy loss, and experimental parameters, we investigate the relationship between the incident and transmitted beam intensities. Key factors such as detector collection angle, convergence semi-angle, and the material properties of samples were analyzed. Our results demonstrate that the relationship between incident and transmitted beam intensities follows the Beer–Lambert law over thicknesses ranging from a few microns to several tens of microns, depending on material composition, electron energy, and collection angles. The linear depth of silicon dioxide reaches 3.9 µm at 3 MeV, about 6 times higher than that at 300 keV. Meanwhile, the linear depth of amorphous ice reaches 17.9 µm at 3 MeV, approximately 11.5 times higher than that at 300 keV. These findings are crucial for advancing the study of thick biological and semiconductor samples using MeV-STEM.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Co-Simulation of PSS/E, OpenDSS, and PSCAD for Power Systems Stability Analysis With Inverter-Based Resources: Preprint

The increasing penetration of inverter-based resources (IBRs) is reshaping the dynamic behavior of power systems. IEEE standard 1547-2018 suggests that distributed energy resources (DERs) should provide grid services such as voltage and frequency supports. On the other hand, dynamic events caused by IBRs such as sub-synchronous oscillation have been reported. These developments necessitate improvement in the current modeling capabilities to better understand the interdependencies within power systems. These include interactions between transmission and distribution systems, among IBRs themselves, and between IBRs and conventional resources. In this paper, we present a co-simulation model integrating PSS/E, OpenDSS, and PSCAD to analyze IBR impacts on the stability of transmission and distribution systems. A key challenge in developing a co-simulation model is ensuring interoperability among different simulators (interfacing and data flow) while maintaining accurate results. Using the developed model, we simulate the impact of IBRs on power systems in two test cases: 1) fault ride through (FRT) capability during a generation trip contingency; 2) IBR-induced sub-synchronous oscillation. The results show the effectiveness of the co-simulation model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Gene-Metabolite Association Prediction with Interactive Knowledge Transfer Enhanced Graph for Metabolite Production

Identifying gene targets for enhancing metabolite production in metabolic engineering is challenging due to the vast research literature and the approximation in genome-scale metabolic model (GEM) simulations. Here, to address this, we propose the Gene-Metabolite Association Prediction task, which automates gene discovery for given metabolite-gene pairs, accompanied by a benchmark dataset of 2474 metabolites and 1947 genes for Saccharomyces cerevisiae (SC) and Issatchenkia orientalis (IO). This task is complicated by incomplete metabolic graphs and metabolic heterogeneity. We introduce an Interactive Knowledge Transfer mechanism based on Metabolism Graphs (IKT4Meta) to enhance prediction accuracy by integrating cross-metabolism knowledge. Using Pretrained Language Models (PLMs) to generate inter-graph links mitigates heterogeneity issues, while intra-graph links are propagated via these anchors. Gene-metabolite predictions are then performed on the enriched graphs integrating multiple microorganisms’ knowledge. Experiments show that IKT4Meta outperforms baselines by up to 12.3% in link prediction.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Effects of Dark Matter Self-interactions on Sagittarius and its Stream

This work explores how assumptions regarding the particle-physics nature of dark matter can alter the evolution of the Sagittarius (Sgr) dwarf spheroidal galaxy and its expansive stellar stream. We run a large suite of N-body simulations to model the infall of a spherically symmetric Sgr-like dwarf, exploring how the presence of dark matter self-interactions impacts its evolution. For a scattering cross section of σ/m χ = 30 cm 2 g −1 (at orbital velocity scales), these interactions result in significantly less stellar mass and little to no dark matter bound to the progenitor at the present day. To isolate the cause of this mass loss, we introduce a novel technique for controlling which pairs of dark matter simulation particles can interact. This enables us to identify ram-pressure evaporation—the scattering of satellite and host dark matter particles—as the primary source of the enhanced mass loss. The rapid disintegration of the Sgr progenitor when self-interactions are allowed alters some key properties of the resulting stellar stream, most dramatically suppressing the presence of a “spur” on the apocenter of the trailing stream arm that correlates with the mass of the satellite at last pericenter. We demonstrate how the effects on the Sgr system scale with the particular choice of self-interaction cross section, which affects the degree of ram-pressure evaporation. These findings generalize beyond the Sgr system, underscoring that dwarf stellar streams and dwarf galaxies with close passages may serve as sensitive probes for dark matter self-interactions.

Milky Way dark matter halo↗
DOE OSTI2025

A systematic decision-making methodology to formalize the selection of degree of realism in screening analysis of probabilistic risk assessment

In the nuclear power domain, Probabilistic Risk Assessment (PRA) is used to inform decision-making for Nuclear Power Plants (NPPs). Recently, there has been an increase in the utilization of modeling and simulation (M&S) to support the estimation of PRA inputs. Risk analysts should carefully select the PRA items that require M&S and their degree of realism (DoR) with consideration of the required resources. To support this selection, this article formulates a systematic decision-making approach for the DoR selection. The DoR selection is made based on two predictive decision-making attributes: the predicted differences in safety risk estimate (ΔSaRi) and the cost of analysis (ΔCAN). This research also develops and quantifies causal models to estimate ΔSaRi and ΔCAN. The causal model-based prediction of ΔSaRi and ΔCAN helps reduce the trial-and-error nature of the DoR selection in the PRA screening analysis and provides insights for DoR selection and the gradual refinements of PRA realism. This approach is demonstrated for a case study on fire PRA of NPPs, where an adequate DoR is selected from two fire models: an engineering correlation and a zone model.

Alkhatib, Sari [Department of Nuclear, Plasma, and↗