Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integrated 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Center for Integrated Simulation of Fusion Relevant RF Actuators

This project was part of the “Center for Integrated Simulation of Fusion Relevant RF Actuators” SciDAC-4 project, led by Dr. Paul Bonoli (MIT). Rather than use an acronym (CISFRRFA), the project will be referred to in this document as the “RF-SciDAC4”. The larger SciDAC-4 project goals were to: 1. Develop an integrated simulation of the antenna + sheath + scrape-off-layer + core plasma system which fully utilizes leadership class computing resources to move towards a quantitative predictive capability for the response to RF power. 2. Work closely with the SciDAC-4 Whole Device Modeling (WDM) community to make both our new code development efforts, as well as the established hierarchy of RF tools, available within their environment, and to utilize WDM technologies to implement the couplings below. 3. Validate this predictive capability on appropriately diagnosed experiments including dedicated RF test stands, linear devices, and existing tokamaks. 4. Use these tools to inform design of robust, impurity-mitigating RF heating and current drive sources for future fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Finite-element boundary-integral simulation of thin wires and inhomogeneous penetrable bodies in subsurface multilayered anisotropic media

With the prevailing presence of drilling wells near the subsurface in mature oil and gas fields, the application of electromagnetic methods can be particularly challenging where the electromagnetic field is affected by the steel casing. In the past decades, borehole-to-surface and crosswell electromagnetic methods have been utilized for monitoring of reservoir and underground CO 2 storage. This paper presents a unified finite-element boundary-integral (FEBI) method capable of simultaneously modeling the complex electromagnetic interactions between thin metallic wires (representing steel casings) with 3D trajectory and arbitrary 3D inhomogeneous penetrable bodies (such as CO 2 plumes or hydrocarbon reservoirs) within anisotropic multilayered subsurface environments. Unlike existing approaches that treat these components separately or require dense discretization, or are limited to vertical wells, our unified formulation preserves flexible electromagnetic coupling while delivering improved computational efficiency. Assuming the background formation is multilayered anisotropic media, the surface integral equation method is applied to model the thin wires and boundaries of the inhomogeneous bodies. Meanwhile, the finite element method is applied to model the volume of inhomogeneous bodies. Here, the performance of the proposed FEBI method is assessed through comparison with reference numerical results and its practical significance is demonstrated through CO 2 plume monitoring scenarios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Integrated simulations of premagnetized and self-magnetizing dynamic screw pinch-driven MagLIF

Magnetically driven implosions such as in magnetized liner inertial fusion (MagLIF) on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that dynamically redistribute liner mass during implosion, limiting fusion fuel compression and confinement, which ultimately degrades performance. Driving the implosion with an initially helical drive field that dynamically shifts the direction of the magnetic field surrounding the liner (i.e., a dynamic screw pinch, DSP) is a method proposed to mitigate MRTI in-flight and improve target performance. In DSPs, the axial drive magnetic field component implodes the liner and diffuses through the shocked, melted liner material into the fuel throughout the implosion. Liners can be designed to enable enough axial magnetic flux to diffuse through the liner material to effectively magnetize the fuel region without the need of an initial axial magnetic field (i.e., from external field coils). We present results from three-dimensional radiation-magnetohydrodynamic simulations of MagLIF implosions employing drive magnetic fields composed of axial and azimuthal components (a helical drive field). These simulated DSP-driven MagLIF targets demonstrate improved fuel conditions and thermonuclear yield compared to a traditional MagLIF target implosion. Synthetic x-ray radiography of the imploding liner material and x-ray emission images of the fuel region at the time of peak neutron yield rate indicate superior implosion morphology for DSP-MagLIF implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FRISM (FReight Integrated Simulation Model) [SWR-24-40]

FRISM simulates day-to-day freight activities including end-consumer shopping, distribution channel, and carrier operation with e-commerce dynamics between passenger and freight travel. It outputs an assignment of different shipments to carriers and tour plans of the carriers' vehicles to transport shipments from their origins to their destinations.

Ugirumurera, Juliette↗
DOE OSTI2024

Integrated Simulation of PIP-II at Fermilab

We describe progress towards a community software ecosystem for efficient modeling of the Fermilab PIP-II complex for design validation and virtual test stands.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
DOE OSTI2025

RE-INTEGRATE EMT Simulation Software: DAE Solvers and Automation

Existing electromagnetic transient (EMT) simulation tools face challenges in accelerating EMT simulations, especially for very large-scale power networks. To tackle this issue, next generation EMT simulation tools such as RE-INTEGRATE EMT are being researched upon. Such tools should be equipped with automation capabilities and advanced numerical differential-algebraic equation (DAE) solvers. In this paper, the DAE solvers incorporated within the RE-INTEGRATE EMT simulation tool are discussed. In particular, a modified ODEINT-based DAE solver and the ARKODE solver from SUN-DIALS are leveraged within RE-INTEGRATE EMT. In addition, the automation implemented within RE-INTEGRATE EMT to automate the DAE generation (replacing the need of manual discretization and assembling DAEs) is discussed. Different use cases were implemented using the RE-INTEGRATE EMT tool and were validated with respect to baseline simulations.

Marthi, Phani Ratna Vanamali [ORNL] (ORCID:0000000↗
DOE OSTI2025

RE-INTEGRATE EMT Simulation Tool: Input Data Processing Layer for Bulk Power System

This paper introduces an advanced input data processing layer for EMT simulations of large-scale bulk power systems. The paper proposes two versions of the RE-INTEGRATE EMT simulation tool, RE-INTEGRATE Gen-0 and RE-INTEGRATE Gen-1, which are developed to enhance simulation generalizability, scalability, and accuracy. The framework leverages a generic class design for components to incorporate linear equations, which are generated by discretizing the Differential-Algebraic Equations (DAEs) that represent the dynamics of the components. In addition, the framework employs a parsing algorithm that parses a power system’s raw and dyr files to generate a connectivity graph which is then traversed to form the overall system’s dynamics. The proposed input data processing layer is used to simulate the IEEE 39-bus test system. The obtained results demonstrate the framework’s capability to achieve simulation scalability and accuracy. Further, the results indicate that EMT simulations performed using the proposed automations can effectively handle complex grid configurations.

Mishra, Rahul [ORNL] (ORCID:0000000328205932)↗
DOE OSTI2025

Machine Learned Empirical Numerical Integrator from Simulated Data

Recently, a number of state-of-the-art surrogate machine learning (ML) models have been designed for global weather and climate prediction, which have been trained using reanalysis data products. Reanalysis data products are constructed using numerical model simulations that combine numerical integration of partial differential equations and parameterization schemes. These products are typically only archived and made available using coarsened spatial and temporal resolutions. This study explores the impact of the numerical generation methods used to produce the training datasets and the temporal resolution of those datasets on machine learning surrogate models. Using the nonlinear vector autoregression (NVAR) machine as an explainable ML technique, simple dynamical systems are emulated with ML models trained on data produced by three classical numerical integration schemes. NVAR is validated as a skillful ML method, capable of producing accurate predictions and, more importantly, reconstructing both the underlying dynamics and the numerical integration scheme used to generate the training data. However, the machine fails to generalize predictions on unseen test data generated by different numerical integration schemes, despite the underlying dynamical system being the same. This result provides a word of caution for the growing field of machine learning emulation of weather and climate dynamics. Furthermore, we illustrate using NVAR that training on temporally coarsened data may increase the required complexity of ML models and potentially introduce new numerical challenges. Finally, we discover that empirical integration schemes with arbitrary time-stepping sizes can be constructed directly from the data, which implies a potential for the development of empirical numerical integration schemes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Comparative performance assessment of air-source and ground-source heat pumps using CO₂ and R-410A with water well integration: A simulation study

This study investigates the performance and feasibility of heat pump systems for residential space heating in cold climates, with a particular focus on ground source heat pumps (GSHPs) with water wells. Four configurations are modeled and compared, a CO 2 air-source heat pump (ASHP), an R-410A ASHP, a CO 2 GSHP with water well integration, and an R-410A ground-source heat pump with water well integration. System simulations are conducted at both equipment and whole-building levels, followed by a nationwide analysis across ten representative cities using EnergyPlus. Results show that water-well-coupled GSHPs maintain approximately 87 % of their efficiency and 85.5 % of heating capacity as ambient temperature drops from 47 °F/8.3 °C to −15 °F/-26.1 °C, whereas ASHPs retain only 55 % efficiency and 44.5 % capacity. R-410A systems achieve higher efficiency, while CO 2 systems provide significant environmental advantages, including reduced risks of groundwater contamination from refrigerant leakage. On average, the CO 2 GSHPs deliver 35–45 % heating energy savings relative to ASHPs and demonstrate strong potential as a replacement for gas furnaces in cold climates. This work represents a systematic comparative assessment of CO 2 and R-410A air-source and ground-source heat pumps with water well integration. In conclusion, the findings highlight the technical viability, environmental benefits and deployment potential of CO 2 GSHPs, offering a pathway toward efficient and sustainable residential heating solutions in diverse U.S. climates.

CO 2 GSHP↗
DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Integrating HPC simulations and physical experiments to characterize the effects of gamma radiation on seismic protective devices

Seismic protective systems, composed of seismic isolators and dampers, can substantially reduce the effects of earthquake shaking on nuclear power plants and components therein. To enable the use of these devices to protect equipment inside a plant and close to a source of radiation, the U.S. Department of Energy (DOE) funded a project at the Idaho National Laboratory (INL) and the University at Buffalo to characterize the effects of absorbed gamma dose on their mechanical properties. An early task in the project was to determine the exposure time required in the INL Foss Therapy Services (FTS) 60 Co gamma irradiator to achieve a target absorbed dose in the materials used to construct isolators and dampers, including fluids, polymers, composites, and metals. This task required the novel integration of high-performance computing (HPC), Monte Carlo N-Particle (MCNP) simulations, and irradiation experiments using Fricke dosimetry. An MCNP model of the FTS irradiator at INL was developed and validated using Fricke dosimetry. Simulations of three experiments in the irradiator, two with Fricke vials only and one with Fricke vials and a large-size isolator, predicted the Fricke-measured absorbed dose rate to within 15% in all three cases, providing high confidence in the calculation of the gamma dose absorbed in the materials comprising the seismic protective devices. The simulations demonstrated that the effects of photon scattering on absorbed dose rate in the FTS irradiator are negligible for test articles installed close to the cobalt sources and near the rear of the irradiator. The validated MCNP model of the FTS irradiator is being used to support ongoing DOE-funded experiments on seismic protective devices and could be applied to future, non-seismic-related experiments. In conclusion, the novel validation process successfully deployed for the FTS irradiator at INL could be applied to other irradiators, requiring new MCNP models and simulations, and irradiation experiments using dosimeters.

42 - ENGINEERING↗
DOE OSTI2025

Hydrologic response of artificially drained agricultural watersheds: insights from high-resolution integrated surface/subsurface simulations

Artificial drainage systems comprising subsurface networks of perforated pipes (tile drains) and engineered surface ditches are widely used to remove excess water from poorly drained agricultural regions. Artificial drainage lowers the water table by design but also has inadvertent effects on the watershed-scale hydrologic response with important implications for flood risk and nutrient exports. We investigated the effects of tile drains on watershed-scale hydrologic response in the Portage River, OH, Watershed using a high-resolution physically based integrated surface/subsurface hydrology model with recently developed capabilities to represent artificial drainage. Tile drains were found to enhance streamflow during times of low flow, generally consistent with previous studies. Streamflow flashiness was found to have a non-monotonic dependence on tile spacing with a minimum at intermediate spacings (∼50 m). Flashiness and the event hydrographs for small tile spacing were similar to the situation with no tiles, but flow paths from farm to stream were very different for those two end member cases, emphasizing the limitations of the stream hydrograph in characterizing hydrologic response. For typical tile spacings, peak flow can either be enhanced or attenuated by the presence of tiles, depending on the size of the event and the antecedent meteorological conditions. Tiles enhance peak flow when the event is below a threshold of ∼25 mm or when events arrive in dry conditions. Peak flow is reduced by tiles when events are large and arrive in conditions that are not overly dry. The dependence on event size and antecedent conditions is explained by differences in available storage and flow paths to the streams. These results provide additional insights into how tile drainage modulates event-scale hydrologic response, an important control on flood generation mechanisms and nutrient exports.

agricultural watershed↗
DOE OSTI2025

Refining seasonal performance metrics for room air-conditioning in emerging markets: Integrating building simulations with real-world equipment performance data

Buildings significantly impact worldwide energy consumption, emphasizing the need to reduce the cooling energy demand, especially in warm climates. Minimum energy performance standards (MEPS) and seasonal performance metrics such as the Cooling Seasonal Performance Factor (CSPF) are crucial for improving room air conditioning (RAC) efficiency. However, challenges remain, particularly in emerging markets like Brazil, where seasonal performance metrics have recently been introduced. This study assesses the factors influencing country-level seasonal efficiency metrics and proposes a framework to refine these calculations by considering local climates and expected RAC usage in real-world households via building simulations. Key considerations include outdoor air temperature binning for different climates, RAC usage patterns (i.e., daytime and nighttime usages), envelope thermal performance of households, and urban heat island (UHI) effects. The results reveal that CSPF values can vary significantly based on climate conditions, with observed CSPF ranging from 4.10 to 11.59 Wh/Wh across 577 Brazilian climates. The inclusion of UHI effects led to a reduction in CSPF values by up to 29% during nighttime operations in hot urban areas. Additionally, building envelope efficiency showed contrasting impacts on RAC performance, with CSPFs reaching up to 15.35 Wh/Wh under specific optimized conditions. These findings highlight the need for transparent policymaking in RAC performance databases, facilitating the application of approaches like those proposed in this study and supporting diverse stakeholders in decision-making.

Bavaresco, Mateus↗
DOE OSTI2025

RE-INTEGRATE EMT Simulation Software: Graph Convolutional Network for Sparse Matrix Pattern Detection

The increasing complexity of power networks, driven by proliferation of inverters, presents analytical challenges that simplified models often fail to capture, necessitating Electromagnetic Transient (EMT) simulations. EMT models are represented as discretized differential-algebraic equations (DAEs), forming a linear system Ax = b that is computationally intensive to solve. Due to inherent sparsity of adjacency matrix A, distinct patterns emerge that, when accurately identified, enable efficient solver selection to minimize computation time. However, identifying ideal pattern is complicated by numerous reordering algorithms and limited structural insights. To address this, we introduce a Graph Convolutional Network (GCN) model for classifying sparse matrix patterns common in power system analysis. The model, achieving 96% test accuracy, is validated using PV plant models of 125 MW capacities connected to New England 39-bus transmission system (TS), and further scaled to a 4,992-bus network with 384 PV plants, yielding 191, 616 × 191, 616 sized A matrix. For all cases, the GCN model accurately identifies the matrix’s intrinsic sparse pattern, demonstrating its potential to enhance solver performance in EMT analysis.

Hossain, Md Rifat [Florida International Universit↗
DOE OSTI2024

A Framework for Integrating Quantum Simulation and High Performance Computing

Scientific applications are starting to explore the viability of quantum computing. This exploration typically begins with quantum simulations that can run on existing classical platforms, albeit without the performance advantages of real quantum resources. In the context of high-performance computing (HPC), the incorporation of simulation software can often take advantage of the powerful resources to help scale-up the simulation size. The configuration, installation and operation of these quantum simulation packages on HPC resources can often be rather daunting and increases friction for experimentation by scientific application developers. We describe a framework to help streamline access to quantum simulation software running on HPC resources. This includes an interface for circuit-based quantum computing tasks, as well as the necessary resource management infrastructure to make effective use of the underlying HPC resources. The primary contributions of this work include a classification of different usage models for quantum simulation in an HPC context, a review of the software architecture for our approach and a detailed description of the prototype implementation to experiment with these ideas using two different simulators (TNQVM & NWQ-Sim). We include initial experimental results running on the Frontier supercomputer at the Oak Ridge Leadership Computing Facility (OLCF) using a synthetic workload generated via the SupermarQ quantum benchmarking framework.

Shehata, Amir [ORNL] (ORCID:0000000224531426)↗