Search NASA⌕ Search

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.

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

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Initial validation of ALFRED: A Monte Carlo code built on Geant4 for TREAT energy deposition

Predicting the energy deposited in the specimen during an experiment in the Transient Reactor Test (TREAT) Facility is a complex problem due to the nature of the transients occurring in the reactor. In addition, the many particles contributing to energy deposition have different behavior in time and space. ALFRED, a new Geant4 based code, is developed to transport and simulate each particle generated in the core. This code is verified against OpenMC (Open Monte Carlo) on the Godiva benchmark and a simple TREAT model. Next, the energy deposition in TREAT is calculated: 181.05 ± 0.01 MeV for the “instantaneous” energy deposition (which accounts for the energy deposited within 1 s after neutron emission) in fuel and 189.90 ± 0.01 MeV for the total energy deposition in fuel. We discuss these results in this paper with previous calculations and experimental evaluations. This work demonstrates ALFRED’s potential as a high-fidelity tool for computing the spatial and temporal energy deposition in TREAT paving the way for a better understanding of the energy coupling factors in TREAT.

73 - NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A full-scope, high-fidelity simulator-based hardware-in-the-loop testbed for comprehensive nuclear power plant cybersecurity research

Nuclear power plant (NPP) cybersecurity research often relies on hardware-in-the-loop (HIL) testbeds that integrate real hardware components into simulated environments. These testbeds allow researchers to identify vulnerabilities, evaluate attack impacts, and test security measures in a controlled setting. Furthermore, previous HIL testbeds lacked fidelity to accurately represent real nuclear systems, limiting the scope of cybersecurity analysis. This study presents the creation of a HIL testbed, devised upon a full-scope, high-fidelity NPP simulator, to facilitate realistic and comprehensive cybersecurity research. To demonstrate its capabilities, the control logic for the steam generator water level was migrated from the simulator to an external programmable logic controller. As a practical application of the developed testbed, supply chain attack scenarios were simulated by injecting malicious code into the controller logic, and the effects of manipulating sensor inputs and control commands were observed. While this HIL testbed provides more detailed simulations, enhanced realism, and wider applicability compared to other options utilizing a less complex simulator, it is also more intricate and costly. For this reason, we include a detailed comparison with some alternative architectures to aid fellow researchers and practitioners in the selection of a suitable HIL architecture based on specific research objectives.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Relations between Shot Noise, Gain Bandwidth, and Saturation of Instabilities

There are numerous instabilities present in charged particle beams that undergo exponential growth and reach saturation. In various applications, such as free-electron lasers or micro-bunching light sources, achieving saturation is desirable. Conversely, there are applications where these instabilities are utilized as linear broad-band amplifiers for signals embedded in the charged beam. In the latter scenario, the saturation of an instability induces non-linear distortions in the imprinted signal, thereby limiting the useful range of such amplifiers. Accurate evaluation of these instabilities necessitates a complete and comprehensive modeling approach that includes shot noise within the beam. Unfortunately, such modeling is not always feasible or practical. In this paper, we introduce a methodology utilizing the frequency and bandwidth of the instability as key parameters. Through this, we derive an estimation for the range of linear instability growth. Our derivation is conducted in a model-independent manner, making it applicable to a broad spectrum of instabilities. To validate our approach, we employ established and thoroughly benchmarked simulations with a free electron laser (FEL) code as well as self-consistent 3-dimensional simulation of plasma-cascade instability using code SPACE.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Milestone in predicting core plasma turbulence: successful multi-channel validation of the gyrokinetic code GENE

On the basis of several recent breakthroughs in fusion research, many activities have been launched around the world to develop fusion power plants on the fastest possible time scale. In this context, high-fidelity simulations of the plasma behavior on large supercomputers provide one of the main pathways to accelerating progress by guiding crucial design decisions. When it comes to determining the energy confinement time of a magnetic confinement fusion device, which is a key quantity of interest, gyrokinetic turbulence simulations are considered the approach of choice – but the question, whether they are really able to reliably predict the plasma behavior is still open. The present study addresses this important issue by means of careful comparisons between state-of-the-art gyrokinetic turbulence simulations with the GENE code and experimental observations in the ASDEX Upgrade tokamak for an unprecedented number of simultaneous plasma observables.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An experimentally informed design process for future inertial confinement fusion facilities

The achievement of ignition in the laboratory has renewed interest in defining the requirements for a future high-gain inertial confinement fusion (ICF) facility. Our best chance of predicting future ICF performance is with 3-D radiation hydrodynamic simulations that have been benchmarked against experimental data, but their high computational cost is prohibitive for use in practical design studies. We introduce a hierarchical approach where 3-D simulations are tuned to match experimental measurements and used to train 3-D degradation models in 1-D simulations allowing for accurate predictions over the entire OMEGA direct-drive database. A genetic algorithm was used in combination with the trained 1-D simulations to search for optimal direct-drive implosion designs at driver energies ranging from 20 kJ to 10 MJ. As the fidelity of 3-D codes improves, this approach will provide a viable experimentally informed tool for defining the next ICF facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Photochemically Induced Acousto-optics Fluid Simulations

PIAFS is a finite-difference code to solve the compressible Navier-Stokes equations with chemical heating on Cartesian grids. It models chemical reactions of air (oxygen and carbon dioxide) with ozone subject to radiation. It uses a high-order WENO spatial discretization and explicit Runge-Kutta time integration. It is capable of parallel simulations using MPI. The code is written in C/C++.

Oudin, AlbertineN [Lawrence Livermore National Lab↗
DOE OSTI2025

Magnetic geometry effects on turbulent density fluctuations in Wendelstein 7-X

Line-integrated density fluctuation measurements from the core of the optimised stellarator Wendelstein 7-X (W7-X) are studied in the context of magnetic field geometry changes. The amplitude as well as the spectral distribution of fluctuations is affected by varying the magnetic field configuration. Weaker turbulent fluctuations are observed in configurations with higher rotational transform and lower magnetic mirror depth. Simulations with the gyrokinetic codes stella and GENE-3D are used to investigate the role of magnetic geometry for the experimental observations. The stabilising effect of the high rational transform configuration is reproduced by simulations and can be associated with stronger shaping of the magnetic field. This leads to increased finite Larmor radius and local shear effects, which limit the toroidal extent of ion temperature gradient (ITG) modes. For W7-X configurations with higher magnetic mirror depth, gyrokinetic theory and simulations predict weaker fluctuations, since they are closer to having the maximum-J property, which has a stabilising effect on trapped-electron mode and ITG. The opposite trend is observed experimentally, which can therefore not be explained by geometry effects on turbulence alone.

Wendelstein-7X↗
DOE OSTI2025

X-Ray Spectra from General Relativistic Radiation Magnetohydrodynamic Simulations of Thin Disks

We compare X-ray emission from several general relativistic, multifrequency, radiation magnetohydrodynamic simulations of thin black hole accretion disks with different accretion rates and spins. The simulations were performed using the M1 closure scheme, resolved with 12 frequency (energy) bins logarithmically spaced from 5 × 10 −3 to 5 × 10 3 ​​​ keV. We apply a general relativistic Monte Carlo transport code to postprocess the simulation data with greater fidelity in frequency resolution and Compton scattering treatment. Despite the relatively few energy bins and Kompaneets approximation to Compton scattering utilized in the M1 method, we find generally good agreement between the methods. Both produce prominent thermal profiles with peaks around 2–2.5 keV, where agreement is particularly strong and representative of the soft state. Both also find weaker (lower-luminosity) thermally sourced emission extending out to 100 keV due to the hotter innermost regions of the disks. Inverse Compton scattering becomes increasingly effective at hardening spectral outputs with increasing black hole spin, and becomes the dominant mechanism for photons that escape with energies between 10 to several hundred keV. At very high rates of spin, the radiation flux in this upscattered component becomes comparable to the thermal flux, a phenomenon typically associated with intermediate states. Beyond 10 4 keV, we observe faint, free–free emission from hot, optically thin coronal regions developing near the horizon, common to both spinning and nonspinning black holes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The shape of convection in 2D and 3D global simulations of stellar interiors

Context. Theoretical descriptions of convective overshooting in stellar interiors often rely on a basic one-dimensional parameterization of the flow called the filling factor for convection. Several different definitions of the filling factor have been developed for this purpose, based on: (1) the percentage of the volume, (2) the mass flux, and (3) the convective flux that moves through the boundary. Aims. We examine these definitions of the filling factor with the goal of establishing their ability to explain differences between 2D and 3D global simulations of stellar interiors that include fully compressible hydrodynamics and realistic microphysics for stars. Methods. We study convection and overshooting in pairs of identical two-dimensional (2D) and three-dimensional (3D) global simulations of stars produced with MUSIC , a fully compressible, time-implicit hydrodynamics code. We examine pairs of simulations for (1) a 3 M ⊙ red giant star near the first dredge-up point, (2) a 1 M ⊙ pre-main-sequence star with a large convection zone, (3) the current sun, and (4) a 20 M ⊙ main-sequence star with a large convective core. Results. Our calculations of the filling factor based on the volume percentage and the mass flux indicate asymmetrical convection near the surface for each star with an outer convection zone. However, near the convective boundary, convective flows achieve inward-outward symmetry for each star that we study; for 2D and 3D simulations, these filling factors are indistinguishable. A filling factor based on the convective flux is contaminated by boundary-layer-like flows, making a theoretical interpretation difficult. We present two possible new alternatives to these frequently used definitions of a filling factor, which instead compare flows at two different radial points. The first alternative is the penetration parameter of Anders et al. (2022, ApJ, 926, 169). The second alternative is a new statistic that we call the plume interaction parameter. We demonstrate that both of these parameters captures systematic differences between 2D and 3D simulations around the convective boundary.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Radiation asymmetry during shattered pellet and massive gas injection in DIII-D

Infrared thermography of the first wall in DIII-D is used to show the toroidal peaking of plasma radiation during mitigated disruptions with shattered pellet injection (SPI) and with massive gas injection (MGI). During MGI, the radiation peak location is shown to be due to the n = 1 magnetohydrodynamic (MHD) activity and continuously variable in toroidal phase based on externally applied error fields, consistent with previous experimental results. Furthermore, the measured toroidal peaking factor (TPF) is in agreement with the previous estimates based on radiometry, with a total TPF of 1.07 ± 0.05 when integrated over the entire duration of the disruption, and a value of 1.3 ± 0.1 during the thermal quench (TQ). For SPI, the location of the radiation peak is found to be determined predominantly by the particle source location, and thus can be varied on DIII-D by injecting from either of two toroidally separated injectors. 3D non-linear MHD simulations with the NIMROD code support this picture. The localized SPI particle source leads to higher peaking than for MGI, with a total TPF of 1.3 ± 0.1 over the full disruption and a TQ TPF of 1.9 +0.5/−0.3, a value also consistent with the NIMROD simulations. This TQ value is close to or potentially exceeding the allowable limit for the previously planned Be first wall in ITER, pointing to the importance of radiation asymmetries depending on the choice of first wall material.

Physics - Plasma physics↗
DOE OSTI2025

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

APOLLO: a facility-scale differentiable virtual accelerator at Fermilab FAST/IOTA

As the design complexity of modern accelerators grows, there is more interest in using advanced simulations that have fast execution time or yield additional insights like gradients. The FAST/IOTA facility has been working on implementing and experimentally validating an end-to-end digital twin that is both fast and gradient-aware, allowing for rapid prototyping of new software and experiments with minimal beam time costs. Our framework integrates physics and ML codes for linac and ring simulation through a set of generic interfaces between surrogate and physics-based sections. To reproduce device inputs and outputs, system state is exposed as a deterministic event loop in a specialized discrete event simulator architecture. Because Fermilab is undergoing control system transition, several APIs were implemented as final user interfaces - a fully asynchronous EPICS soft IOC, a gRPC-based Data Pool Manager (DPM), and legacy ACNET protocols. We discuss implementation details as well as challenges handling live data assimilation and future plans to extend modelling to main complex proton accelerators like PIPII and Booster.

Kuklev, Nikita [Fermilab]↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗
DOE OSTI2025

Gas-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the pebble-bed high-temperature gas-cooled reactor (PB-HTGR) reference plant models for the United States Nuclear Regulatory Commission. These models serve as the foundation for the future development of detailed design evaluation models based on license applications. The reference plant models have been developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, which is the code suite proposed for non-light-water reactor systems safety analysis. It incorporates various simulation tools developed by the Nuclear Energy Advanced Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of the work that was performed was to assess BlueCRAB’s level of readiness for modeling a PB-HTGR. To do so, we first developed numerical models in BlueCRAB that include the key physics for this technology to ensure an adequate level of fidelity for modeling PB-HTGR core performance and for performing multiphysics simulations for equilibrium core conditions and different accident scenarios. Then we simulated transient scenarios, including depressurized and pressurized loss of forced cooling accidents, over-cooling, and control rod withdrawal events with delayed and prompt supercritical reactivity insertions. The analysis in this report includes comparisons of the 2D thermal fluid porous media models in Pronghorn and SAM, and comparisons of coupled SAM/Griffin/SAM and coupled Pronghorn/Griffin for depressurized and pressurized loss of forced cooling, over-cooling, and control rod withdrawal events. In addition, we compare 3D, 2D, and 0D/PKE neutronic models for the two control rod withdrawal scenarios with coupled Pronghorn/Griffin. The comparisons show that the BlueCRAB models lead to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields, such as neutron flux, power, temperature, density, pressure, and velocity, are within the expected ranges, and their distributions can be explained by the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. This work further highlights the need for flexible tools with various levels of fidelity to cover the breadth and depth of needs that may arise in future technical evaluations of the PB-HTGR. We believe that the BlueCRAB capabilities will be a significant asset for confirmatory analyses in order to resolve important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Nonlinear light-output calibration of the oxygenated xylene scintillators used in OMEGA neutron time-of-flight spectrometers

Neutron time-of-flight (nTOF) spectrometers are essential instruments for measuring and evaluating the performance of inertial confinement fusion implosions. The neutron spectrometers utilized for the OMEGA laser include two liquid-based scintillators, each consisting of a large volume filled with xylene that is coupled to four photomultiplier tubes. Analysis of the signal from these detectors requires detailed knowledge of the scintillator’s light output, which is needed to fit the nTOF spectrum, from which the neutron energy spectrum is informed. The light output is nonlinearly proportional to the neutron energy, which, in turn, affects the interpretation of the neutron energy spectrum from a TOF signal. A recent campaign on OMEGA was performed to calibrate the xylene detectors and infer the shape of the light-output curve. The campaign utilized materials with increasing Z placed in the OMEGA target chamber to initiate scattering events with the 14 MeV fusion neutrons. This process leads to the production of backscatter neutrons of varying energies that appear as peaks in the nTOF data. Simulations using a neutron transport code were combined with the measured deuterium–tritium neutron yields to calculate the expected backscattered neutron yields from the well-known scattering cross sections of each material. Furthermore, the neutron-energy dependent light output of the scintillator inferred from the experiment is compared to the light-output curve simulated with a neutron transport code for the following neutron energies: 1.5, 2.5, 6, and 14 MeV.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Study of the heat pipe vapor flow using the MOOSE application Sockeye and the CFD code Nek5000

Heat pipes are efficient heat transfer devices used in various applications, including nuclear microreactors. Heat pipe-cooled microreactors offer advantages in size and cost, which can significantly accelerate their deployment and adoption. Developing accurate models of heat pipes is crucial for a heat pipe-cooled microreactor design and operation. In this work a verification and validation of the heat pipe application Sockeye is presented. Results are compared against RANS simulations, performed with the code Nek5000, and the SPHERE experiment, conducted at INL.

42 ENGINEERING↗