SEARCH · Search NASA

Results for “simulation codes”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Dataset of simulated vibrational density of states and X-ray diffraction profiles of mechanically deformed and disordered atomic structures in Gold, Iron, Magnesium, and Silicon

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

DOE OSTI2024

CFD unified approach under Eulerian–Lagrangian framework for methanol and gasoline direct injection sprays in evaporative and flash boiling conditions

Innovative synthetic fuels for advanced propulsion systems, such as methanol and ammonia, and synthetic blended fuels (E00, E10, and E30), known for their high volatility, are often injected directly into combustion chambers. It follows that Eulerian–Lagrangian spray models need to accurately capture the spray collapse as a consequence of flash boiling onset and be capable of proficiently handling the preferential evaporation of multi-component fuels in evaporative scenarios. So, we performed the assessment of an Eulerian–Lagrangian CFD code for simulating methanol and E00 gasoline blend sprays in both early and late injection conditions involving flash boiling conditions and preferential evaporation. The adoption of an effervescent breakup model and of a non-equilibrium phase transition model for the discrete phase allows the adoption of a setup that is almost completely free from specific constant tuning, especially for what concerns the breakup model. We validated the simulations using experimental PLV maps of methanol and E00 sprays issued from the ECN Spray M injector. The results highlight a significantly different morphology of the methanol spray compared to the E00 one under late injection conditions. Under stratified combustion, low-volatile fuels are likely to be ignited first, and the flame propagates toward the high-volatile fuels. In conclusion, the spray collapse was also correctly reproduced, inducing the presence of a low-pressure zone and modifying the spray morphology.

E00↗

DOE OSTI2024

Dependence of high- Z redeposition on the field-to-surface pitch angle and other sheath parameters in tokamaks

Accurately predicting redeposition is vital for high-Z plasma-facing component (PFC) survivability in magnetic confinement fusion. In this study, we categorize high-Z redeposition into three mechanisms: geometric-driven (prompt), sheath-driven (local), and scrape-off-layer-driven (far) redeposition. To investigate these mechanisms, we employ Monte Carlo transport codes to simulate azimuthally symmetric tungsten source erosion and redeposition in a tokamak. By iteratively analyzing critical parameters, we evaluate redeposition scaling for each mechanism. Specifically, we investigate the impact of magnetic-field-to-PFC pitch angle assumptions on PFC losses into the scrape-off layer. Our findings reveal significant pitch angle sensitivity due to an asymmetric prompt vs local redeposition trade-off. These results enhance our understanding of redeposition phenomena in fusion plasma environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Numerical challenges in modeling gravothermal collapse in Self-Interacting Dark Matter halos

When dark matter has a large cross section for self scattering, halos can undergo a process known as gravothermal core collapse, where the inner core rapidly increases in density and temperature. To date, several methods have been used to implement Self-Interacting Dark Matter (SIDM) in N-body codes, but there has been no systematic study of these different methods or their accuracy in the core-collapse phase. In this paper, we compare three different numerical implementations of SIDM, including the standard methods from the GIZMO and Arepo codes, by simulating idealized dwarf halos undergoing significant dark matter self interactions (σ/m = 50 cm 2 /g). When simulating these halos, we also vary the mass resolution, time-stepping criteria, and gravitational force-softening scheme. The various SIDM methods lead to distinct differences in a halo's evolution during the core-collapse phase, as each results in spurious scattering rate differences and energy gains/losses. The use of adaptive force softening for gravity can lead to numerical heating that artificially accelerates core collapse, while an insufficiently small simulation time step can cause core evolution to stall or completely reverse. Additionally, particle numbers must be large enough to ensure that the simulated halos are not sensitive to noise in the initial conditions. Even for the highest-resolution simulations tested in this study (10 6 particles per halo), we find that variations of order 10% in collapse time are still present. Here, the results of this work underscore the sensitivity of SIDM modeling on the choice of numerical implementation and motivate a careful study of how these results generalize to halos in a cosmological context.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Relativistic gas accretion onto supermassive black hole binaries from inspiral through merger

Accreting supermassive black hole binaries are powerful multimessenger sources emitting both gravitational and electromagnetic (EM) radiation. Understanding the accretion dynamics of these systems and predicting their distinctive EM signals is crucial to informing and guiding upcoming efforts aimed at detecting gravitational waves produced by these binaries. To this end, accurate numerical modeling is required to describe both the spacetime and the magnetized gas around the black holes. In this paper, we present two key advances in this field of research. First, we have developed a novel 3D general relativistic magnetohydrodynamics (GRMHD) framework that combines multiple numerical codes to simulate the inspiral and merger of supermassive black hole binaries starting from realistic initial data and running all the way through merger. Throughout the evolution, we adopt a simple but functional prescription to account for gas cooling through photon emission. Next, we have applied our new computational method to follow the time evolution of a circular, equal-mass, nonspinning black hole binary for ∼200 orbits, starting from a separation of 20⁢𝑟 𝑔 and reaching the postmerger evolutionary stage of the system. We have shown how mass continues to flow toward the binary even after the binary “decouples” from its surrounding disk, but the accretion rate onto the black holes diminishes. We have identified how the minidisks orbiting each black hole are slowly drained and eventually dissolve as the binary compresses. We confirm previous findings that the system’s luminosity decreases by a factor of a few during inspiral; however, we observe an abrupt increase by ∼50% in this quantity at the time of merger, likely accompanied by an equally abrupt change in spectrum. Lastly, we have demonstrated that during the inspiral, fluid ram pressure regulates the fraction of the magnetic flux transported to the binary that attaches to the black holes’ horizons.

Accretion disk & black-hole plasma↗

DOE OSTI2024

CHARM-SYCL & IRIS: A Tool Chain for Performance Portability on Extremely Heterogeneous Systems

Performance portability is becoming crucial as high-performance computing systems become increasingly heterogeneous. We have many options for CPUs and accelerators (e.g., GPUs) but also for non-Von Neumann architectures such as field-programmable gate arrays. This paper presents the CHARM-SYCL unified programming environment for multiple accelerator types as a performance-portable programming environment. It uses the IRIS library developed at Oak Ridge National Laboratory as the back end accelerator runtime. IRIS has a high-performance scheduler to distribute tasks across accelerators. This design allows us to run an application from the same source on multiple systems with multiple configurations. We provide three types of portability with CHARM-SYCL: Portable Workflow, Compiler and Runtime Portability, and Application and Performance Portability. We implement a Monte Carlo simulation benchmark code on the CHARM-SYCL execution environment and demonstrate that our programming environment can accommodate extremely heterogeneous systems.

Fujita, Norihisa↗

DOE OSTI2024

Bayesian calibration of stochastic agent based model via random forest

SAND2024-11403O The Bayesian Calibration of Stochastic Agent-based Model via Random Forest is a code that was developed in concurrence with an article by the same name that was written for a journal. The code reproduces simulation results and plots from the article. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

DOE OSTI2025

Impact of increased smoothing by spectral dispersion bandwidth on stimulated Brillouin scattering in laser driven Hohlraums

Experiments were conducted at the National Ignition Facility (NIF) to investigate the impact of increased smoothing by spectral dispersion (SSD) bandwidth on the production of stimulated Brillouin scattering (SBS) within an indirect-drive inertial confinement fusion (ICF) Hohlraum. This was done in a subscale gold Hohlraum driven by 192 laser beams depositing 1.1 MJ of energy. The laser bandwidth was increased from 45 to 118 GHz (before frequency tripling) on the 30° and 50° cones, where backscatter could be measured using the Full Aperture Backscatter Station (FABS). It was expected that this 2.6-fold increase in bandwidth would suppress SBS generated within the Hohlraum plasma and lower the backscattered SBS energy on the 50° cones by over a factor of four. Experimental results, however, show that this SSD change only reduced the 50° cone SBS during the main capsule drive by -18±31% and -4.5±7.8% over the entire pulse. This is small compared to expected shot-to-shot SBS reproducibility (∼30%), such that the result can be considered within normal performance fluctuations. New 3D parallel paraxial code (pF3D) simulations, accounting for beam refractive intensification reproduce this result, suggesting that closer to 300 GHz of bandwidth would have been required to mitigate SBS to the expected level. Delivering such a high bandwidth is not feasible at NIF when operating at high peak power and would potentially prevent NIF’s ability to use cross-beam energy transfer (CBET) for implosion symmetry tuning.

Physics↗

DOE OSTI2024

Conceptual Design of Neutron Sensor Qualification Device

The Advanced Sensors and Instrumentation Program at Idaho National Laboratory has been formulating strategies to qualify sensors for use in nuclear environments, particularly in irradiation experiments and advanced reactors. When qualifying neutron sensors for use in high-temperature environments, the wide range of neutron flux levels and representative energy spectra presents significant challenges. This paper discusses the development of the Neutron Sensor Qualification Device (NQD), which is designed to test neutron sensors in high temperature controlled environments with known neutron spectra, addressing the spatial and spectral complexities of neutron fluxes in reactor cores. The proposed NQD will be situated in the exposure room at the Armed Forces Radiobiology Research Institute, thus affording a unique capability to expose sensors to high neutron and gamma fluxes. To achieve thermal control, the device will utilize a radiation-hardened tube furnace, accommodating multiple sensors and neutron activation dosimetry wires. Titanium, iron, and cobalt dosimeter wires are chosen from the American Society for Testing and Materials and International Reactor Dosimetry and Fusion File libraries as references for providing energy-dependent fluence measurements. The design ensures precise sensor positioning to minimize mutual shielding and flux perturbation, which are evaluated via Monte Carlo N particle Transport Code (MCNP) simulations. These simulations have informed the development of guidelines on sensor placement within the NQD. The NQD is essential to the qualification of neutron sensors for advanced reactor technologies. It enables controlled testing of a statistically significant number of sensors, thereby supporting assessments of sensor performance across various neutron flux levels and temperatures. This paper highlights the detailed planning for the NQD prototype, along with its inaugural irradiation (scheduled for fiscal year [FY] 2025). The results from this initial testing will be fundamental in evaluating the device’s performance and establishing measurement uncertainty for in-pile neutron sensor measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

OpenEdge: A collaborative, open-source, multi-purpose direct simulation Monte Carlo for plasma simulation in magnetic fusion environments

OpenEdge is a collaborative, open-source, object-oriented Direct Simulation Monte Carlo (DSMC) code, designed specifically for plasma simulations in magnetic fusion environments. Here, the code features include advanced structures, robust capabilities, and an effective parallelization strategy, all of which significantly enhance performance. It includes specialized modules for managing complex particle interactions, including collisions, ionization/recombination, and reflection/sputtering. Benchmarks and performance analyses have confirmed its efficiency and scalability. Versatile and adaptable, OpenEdge is applied across a broad spectrum of plasma-material interaction studies and charged particle transport in various fusion research settings.

Boundary plasma↗

DOE OSTI2026

SCALE HTR-PROTEUS Benchmark Model

This dataset contains input and result files of computational simulations of HTR-PROTEUS benchmark with the latest version of SCALE code system. The simulations cover criticality control rod worth calculations as well as sensitivity analysis and uncertainty quantification. Users wanting to reproduce results from this dataset are required to obtain a license to the SCALE code system for which details on the distribution can be found here: https://www.ornl.gov/scale/releases

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2025

A GPU ‐Accelerated 3D Unstructured Mesh Based Particle Tracking Code for Multi‐Species Impurity Transport Simulation in Fusion Tokamaks

ABSTRACT This paper presents the multi‐species global impurity transport capability developed in a GPU‐accelerated fully 3D unstructured mesh‐based code, GITRm, to simultaneously track multiple impurity species and handle interactions of these impurities with mixed‐material surfaces. Different computational approaches to model particle‐surface interaction or surface response have been developed and compared. Sheath electric field is taken into account by employing a fast distance‐to‐boundary calculation, which is carried out in parallel on distributed or partitioned meshes on multiple GPUs without the need for any inter‐process communication during the simulation. Several example cases, including two for the DIII‐D tokamak, that is, one with the SAS‐V divertor and the other with the collector probes, are used to demonstrate the utility of the current multi‐species capability. For the DIII‐D probe case, the capability of GITRm to resolve the spatial distribution of particles in localized regions, such as diagnostic probes, within non‐axisymmetric tokamak geometries is demonstrated. These simulations involve up to 320 million particles and utilize up to 48 GPUs.

Nath, Dhyanjyoti D. [Scientific Computation Resear↗

DOE OSTI2026

A Simulator for Neyer Tests of Explosives

Explosives and explosive devices such as detonators are typically tested by applying a range of stimuli such as voltage or mechanical shock, and recording binary “detonated/did not detonate” responses. These are analyzed using maximum likelihood or generalized linear models to provide estimates of quantities such as the all-fire and no-fire points. Given that the true threshold for detonation is unknown a priori , sequential design methods are typically used to optimize the set of test points. One popular method, implemented in commercial software, is Neyer’s algorithm. To support simulation and experimental design, we have developed code in the R programming language to duplicate the functions of the Neyer software. We provide code for the simulator along with a description and examples of usage.

42 ENGINEERING↗

DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Self-consistent equilibrium and transport simulations for NSTX-U plasmas enhanced via machine learning surrogate models

The Control-Oriented Transport SIMulator (COTSIM) is an advanced equilibrium and transport code designed for simulating tokamak discharges at computational speeds suitable for control applications. COTSIM’s modular framework enables users to select models that balance accuracy with speed according to specific needs, allowing the code to operate from fast to faster-than-real-time performance levels. This work presents recent enhancements to COTSIM’s predictive accuracy for NSTX-U scenarios, achieved by integrating neural-network-based surrogate models and self-consistent equilibrium calculations. To improve source deposition predictions, a surrogate model for NUBEAM has been incorporated. Additionally, a surrogate model for the Multi-Mode Module (MMM) now supports predictions of anomalous thermal, momentum, and particle diffusivities—key factors for modeling the evolution of temperature and rotation. Each surrogate model was specifically trained for the NSTX-U operational regime to enhance COTSIM’s accuracy while maintaining computational efficiency. Moreover, COTSIM now couples fixed-boundary equilibrium solvers with its transport solvers, enabling self-consistent predictions of plasma profiles and equilibrium evolution over the discharge. Simulation results demonstrate strong agreement between COTSIM and TRANSP predictions for NSTX-U discharges. These substantial advancements expand COTSIM’s utility in model-based control applications for NSTX-U. Potential applications include simultaneous optimization of equilibrium and transport scenarios, integration into digital twins, real-time profile estimation (e.g., temperature and rotation) from limited or noisy measurements, and advanced feedback-based scenario control.

Equilibrium and transport modeling↗

DOE OSTI2025

Application of a Physics-Informed Convolutional Neural Network for Monitoring the Temperature Fields in High-Temperature Gas Reactors

Here, this work presents current advances in applying a physics-informed convolutional neural network (CNN) to evaluate temperature distributions in advanced reactors. Our goal is to demonstrate that the CNN can reconstruct temperature fields within the solid region of a prismatic fuel assembly in a high-temperature gas reactor (HTGR) with sensor data available in only a few cooling channels. Before that, we showcase the superior performance of the physics-informed CNN in comparison to a purely data-driven multilayer perceptron (MLP), considering a canonical heated channel setup. This analysis shows the advantages of our approach and justifies its choice. The datasets employed here are obtained upon numerical simulations performed with codes under the Nuclear Energy Advanced Modeling and Simulation program. This work is important, as industry experience indicates that the assembly material in HTGR concepts is prone to large thermal-mechanical loads nearing operational limits. This makes it crucial to characterize peak temperatures and their distributions near hot spots. Modern thermocouples are unreliable in these types of harsh environments because of the high neutron fluxes and elevated temperatures involved. The CNN-based field reconstruction represents an attractive solution, enabling sensor arrays in less aggressive locations and augmenting indirect predictions for less accessible regions. The results show that the CNN reduces prediction errors by orders of magnitude in comparison to the MLP, considering the simple yet well-representative heated channel case. In the case of the HTGR fuel assembly, the CNN can successfully reconstruct temperature fields over various cooling regimes. Furthermore, we also explore the algorithm’s ability to detect abnormalities. Interestingly, the CNN proves it has the capacity to detect blockage in one of the noninstrumented cooling channels.

Machine learning↗

DOE OSTI2024

Present and future of $\mathcal{C}$osmo $\mathcal{L}$attice

We discuss the present state and planned updates of Cosmo Lattice, a cutting-edge code for lattice simulations of non-linear dynamics of scalar-gauge field theories in an expanding background. We first review the current capabilities of the code, including the simulation of interacting singlet scalars and of Abelian and non-Abelian scalar-gauge theories. We also comment on new features recently implemented, such as the simulation of gravitational waves from scalar and gauge fields. Secondly, we discuss new extensions of Cosmo Lattice that we plan to release publicly. We comment on new physics modules, which include axion-gauge interactions $\phi F\tilde{F}$, non-minimal gravitational couplings $\phi ^2R$, creation and evolution of cosmic-defect networks, and magnetohydrodynamics. Finally, we also discuss new technical features, including evolvers for non-canonical interactions, arbitrary initial conditions, simulations in 2+1 dimensions, and higher-accuracy spatial derivatives.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗