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 721 records · Page 40

Simulations of vertical displacement oscillatory modes and global Alfvén Eigenmodes in JET geometry

Vertical Displacement Oscillatory Modes (VDOM), with frequency in the Alfvén range, are natural modes of oscillation of magnetically confined laboratory plasmas with elongated cross-section. These axisymmetric modes arise from the interaction between the plasma current, which is in equilibrium with currents flowing in external coils, and perturbed currents induced on a nearby conducting wall. The restoring force exerted by these perturbed currents on the vertical motion of the plasma column leads to its oscillatory behavior. An analytic model for VDOM was proposed based on an idealized 'straight tokamak' equilibrium with uniform equilibrium current density. This article introduces the first numerical simulations of VDOM in a realistic JET tokamak configuration, using the extended-MHD code NIMROD and drawing comparisons with Global Alfvén Eigenmodes (GAE). The results show qualitative agreement with analytic predictions regarding mode frequency and radial structure, supporting the identification of VDOM as a fundamental oscillation mode in tokamak plasmas. VDOM and GAE are modeled in a representative JET discharge, where axisymmetric perturbations with toroidal mode number n = 0 driven unstable by fast ions were observed. The two modes are examined separately using a forced oscillator within the NIMROD code, which enables a comparison of their characteristics and helps identify the experimentally observed mode possibly as a GAE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A model for trapping and re-solution regarding intra-granular bubbles in UO 2 , linked to atomic-scale simulations

In the literature, a clear definition of the irradiation re-solution frequency of gas from bubbles in the UO 2 fuel is absent. Moreover, for intra-granular bubbles, a detailed calculation of the cumulated displaced gas quantities in function of the distance from the radius of the bubble after a re-solution event has never been published. The assessment of these two elements is very useful if we want to increase the adherence of fission gas release codes to our present knowledge of the behavior of fission gases. Hence, we suggest to link the definition of the re-solution frequency to atomic-scale simulations. Furthermore, we present the cumulated displaced gas quantities obtained from Molecular Dynamics calculations, from which we have derived a re-solution profile that can be exploited to better consider the irradiation re-solution phenomenon inside Fission Gas Release codes. On top of that, we have built a new trapping/re-solution model for intra-granular bubbles linked to Molecular Dynamics simulations that can be easily incorporated into Fission Gas Release codes. In conclusion, we also check that the model is properly built through the comparison of the new model against a reference.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SICAS, a new integrated divertor-SOL-pedestal-core framework for self-consistent modeling of ion and impurity transport in fusion devices

We present a new framework for core-edge integration studies named SICAS (SOLPS-ITER coupled to ASTRA-STRAHL) which enables high fidelity simulations of the core, edge, and divertor regions encompassing the transport of ions as well as the impurities through the entire plasma domain. SICAS handles the exchanging of the particle and power fluxes as well as transport coefficients to ensure consistency through the codes. An overlapping region is defined from the inside of the separatrix of the plasma to the SOLPS-ITER core boundary allowing for matching profiles and fluxes between the two codes for a self-consistent approach. The results presented here demonstrated the flexibility of SICAS to simulate different configurations, scenarios, divertor geometries, and plasma species with good agreement with DIII-D experimental data. This tool opens new possibilities in integrated modeling of fusion devices integrating all relevant phenomena in the core and the divertor plasmas. These capabilities are required for the interpretation of current experiments as well as the design of new devices.

core-edge integration↗

Initial steady-state core simulation capability or thermal and pool-type molten salt reactors, coupling reactor physics, thermal-hydraulics, and evolving chemistry

This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.

42 - ENGINEERING↗

Berkeley eXtensible Environment (BXE) v3

The Berkeley eXtensible Environment (BXE) provides a cloud environment for hardware designers and computer architects to design, build, and simulate their custom architectures on an on-premises FPGA cluster. Utilizing the Chipyard, MoSAIC, and FireSim frameworks, users are provided an environment where they can assemble SoC designs from an existing library of components or import their own source code. Once their designs are ready, they can utilize the FireSim framework provided by BXE to deploy and simulate their designs on the FPGA. Users aren't limited to a single FPGA; they can deploy multiple instances across multiple FPGAs, acting like a rack of servers, or partition their large design across multiple FPGAs, ganging multiple FPGAs into a single simulated system.

Fatollahi-Fard, Farzin↗

Observation and Numerical Simulation of Cold Ions Energized by EMIC Waves

This is the first report of significant energization (up to 7,000 eV) of low-energy He + ions, which occurred simultaneously with H-band electromagnetic ion cyclotron (EMIC) wave activity, in a direction mostly perpendicular to the ambient magnetic field. The event was detected by the Arase satellite in the dayside plasmatrough region off the magnetic equator on 15 May 2019. The peak energy of the He + flux enhancements is mostly above 1,000 eV. At some interval, the He + ions are energized up to ~7,000 eV. The H-band waves are excited in a frequency band between the local crossover and helium gyrofrequencies and are close to a linear polarization state with weakly left-handed or right-handed polarization. The normal angle of the waves exhibits significant variation between 0° and 80°, indicating a non-parallel propagation. Here, we run a hybrid code with parameters estimated from the Arase observations to examine the He + energization. The simulations show that cold He + ions are energized up to more than 1,000 eV, similar to the spacecraft observations. From the analysis of the simulated wave fields and cold plasma motions, we found that the ratio of the wave frequency to He + gyrofrequency is a primary factor for transverse energization of cold He + ions. As a consequence of the numerical analysis, we suggest that the significant transverse energization of He + ions observed by Arase is attributed to H-band EMIC waves excited near the local helium gyrofrequency.

79 ASTRONOMY AND ASTROPHYSICS↗

MONTE CARLO CROSS SECTION LOOKUP KERNEL FOR THE CEREBRAS WSE-2 IN CSL

This is a small kernel that was used to collect data for an upcoming paper. We would like to have the code be open source so that the reviewers (and then readers) of the paper can see the whole code, and can reproduce/verify our results. This is not a fully featured application, it cannot produce any useful simulation results, it just executes a small abstracted kernel using synthetic data. The purpose of the kernel is to understand the basic performance characteristics of an HPC kernel on novel AI accelerator architectures. The main kernel is written in the CSL coding language for use with the Cerebras WSE-2 AI accelerator. The kernel represented is the Monte Carlo cross section lookup kernel, which is a small kernel used by the Monte Carlo neutral particle transport algorithm. There is also a baseline kernel written in CUDA that we will include in the repository to form a basis for comparing the WSE-2 to GPU.

Tramm, John↗

Inferring demographic and selective histories from population genomic data using a 2-step approach in species with coding-sparse genomes: an application to human data

Abstract The demographic history of a population, and the distribution of fitness effects (DFE) of newly arising mutations in functional genomic regions, are fundamental factors dictating both genetic variation and evolutionary trajectories. Although both demographic and DFE inference has been performed extensively in humans, these approaches have generally either been limited to simple demographic models involving a single population, or, where a complex population history has been inferred, without accounting for the potentially confounding effects of selection at linked sites. Taking advantage of the coding-sparse nature of the genome, we propose a 2-step approach in which coalescent simulations are first used to infer a complex multi-population demographic model, utilizing large non-functional regions that are likely free from the effects of background selection. We then use forward-in-time simulations to perform DFE inference in functional regions, conditional on the complex demography inferred and utilizing expected background selection effects in the estimation procedure. Throughout, recombination and mutation rate maps were used to account for the underlying empirical rate heterogeneity across the human genome. Importantly, within this framework it is possible to utilize and fit multiple aspects of the data, and this inference scheme represents a generalized approach for such large-scale inference in species with coding-sparse genomes.

Soni, Vivak (ORCID:0000000294969562)↗

Integrated Neutronics Modeling for Inertial Fusion Energy Systems: Development and Application to LD-FIRST

Lawrence Livermore National Laboratory (LLNL) is proposing a new Laser Driven Fusion Integration Research and Science Test Facility (LD-FIRST) with the goal of providing an experimental testbed for future Inertial Fusion Energy (IFE) systems. However, IFE systems require detailed and accurate multiphysics modeling to quantify material damage, thermal loading, and tritium breeding within complex chamber environments. This article presents the first step in an integrated multiphysics framework that couples meshed CAD-based geometry within Monte Carlo neutronic simulations to enable high-fidelity analysis of IFE chamber concepts, with future coupling to external codes. The neutronics workflow utilizes OpenMC and its third-party capability to use CAD-based geometries through DAGMC and tally on unstructured meshes with Libmesh to evaluate neutron transport behavior, geometric fidelity, and material performance under reactor-relevant conditions. The use of tailored tallies on unstructured meshes in this framework allows direct transfer without interpolating to CFD simulation tools. Two IFE chambers were evaluated, both conceived by LLNL: HYLIFE-II and Laser IFE (LIFE). This work produced high-fidelity conformal surface and volumetric meshes of the HYLIFE-II and LIFE chambers with mapped spatial insight into material damage, thermal loading, and tritium breeding. The HYLIFE-II model was built utilizing available resources and used as a test case to verify that the neutronics framework can handle complex geometries. The LIFE chamber CAD was provided by LLNL and was the main focus of this work. This work analyzes multiple ternary alloy breeding materials for the LIFE chamber, across different 6 Li enrichments to produce data relevant to the LD-FIRST project. This work also investigates the level of model fidelity for the LIFE chamber, and results show that inclusion of detailed first wall and coolant structures increased the predicted tritium breeding ratio (TBR) by ~30%, highlighting the sensitivity of tritium breeding and the need for a high-fidelity simulation framework for IFE chambers. These developments provide a scalable toolset for the design and optimization of next-generation IFE chambers, forming a solid foundation for future coupled multiphysics analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Enabling Parallel Performance and Portability of Solid Mechanics Simulations Across CPU and GPU Architectures

Efficiently simulating solid mechanics is vital across various engineering applications. As constitutive models grow more complex and simulations scale up in size, harnessing the capabilities of modern computer architectures has become essential for achieving timely results. This paper presents advancements in running parallel simulations of solid mechanics on multi-core CPUs and GPUs using a single-code implementation. This portability is made possible by the C++ matrix and array (MATAR) library, which interfaces with the C++ Kokkos library, enabling the selection of fine-grained parallelism backends (e.g., CUDA, HIP, OpenMP, pthreads, etc.) at compile time. MATAR simplifies the transition from Fortran to C++ and Kokkos, making it easier to modernize legacy solid mechanics codes. We applied this approach to modernize a suite of constitutive models and to demonstrate substantial performance improvements across different computer architectures. This paper includes comparative performance studies using multi-core CPUs along with AMD and NVIDIA GPUs. Results are presented using a hypoelastic–plastic model, a crystal plasticity model, and the viscoplastic self-consistent generalized material model (VPSC-GMM). The results underscore the potential of using the MATAR library and modern computer architectures to accelerate solid mechanics simulations.

Morgan, Nathaniel (ORCID:0000000276118449)↗

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗

Efficient Simulation of Logical Magic State Preparation Protocols

Developing space- and time-efficient logical magic state preparation (MSP) protocols will likely be an essential step toward building a large-scale fault-tolerant quantum computer. Motivated by this need, we introduce a scalable method for simulating logical MSP protocols under the standard circuit-level noise model. When applied to protocols based on code-switching, magic state cultivation, and magic state distillation, our method yields a complexity polynomial in (i) the number of qubits and (ii) the nonstabilizerness, e.g., stabilizer rank or Pauli rank, of the target encoded magic state. The efficiency of our simulation method is rooted in a curious fact: every circuit-level Pauli error in these protocols propagates to a Clifford error at the end. This property is satisfied by a large family of protocols, including those that repeatedly measure a transversal Clifford that squares to a Pauli. We provide a proof-of-principle numerical simulation that prepares a magic state using such logical Clifford measurements. Our work enables practical simulation of logical MSP protocols without resorting to approximations or resource-intensive state-vector simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Validation of the SCALE/Polaris–PARCS Code Procedure With the ENDF/B-VII.1 AMPX 56-Group Library: Boiling Water Reactor

The SCALE/Polaris–PARCS code procedure has been used in the confirmatory analysis for boiling water reactors by the US Nuclear Regulatory Commission. In this study, the SCALE/Polaris v6.3.0–PARCS v3.4.2 code procedure with the Evaluated Nuclear Data File (ENDF)/B-VII.1 AMPX 56-group library was validated by comparing the simulated results with the measured data for operating boiling water reactors, including Peach Bottom Unit 2 cycles 1–3, Hatch Unit 1 cycles 1–3, and Quad Cities Unit 1 cycles 1–3. The uncertainties and biases of the SCALE/Polaris–PARCS code package for boiling water reactor physics analysis were evaluated in the validation for key nuclear parameters such as reactivity and traversing in-core probe data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparison of interlaminar damage modeling strategies for hybrid composite/aluminum laminates subjected to low-velocity impact

Low-velocity impact of hybrid metal-composite structures was investigated experimentally and computationally. Composite laminates consisting of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) were joined with a 6061-T6 aluminum plate using an epoxy adhesive. Two variations of the structure were studied; one consisting of all plies oriented at 0° and one consisting of all plies oriented at 45°. A drop tower was used to impact structures at a range of energies, including energies above and below the threshold at which the aluminum layer was perforated. Numerical simulations were implemented using Sierra/SM, an in-house transient dynamics finite element code developed at Sandia National Laboratories. A Hosford plasticity model was used to describe the response of the aluminum layer. A newly implemented orthotropic continuum damage mechanics (CDM) constitutive model was used to represent the composite laminate. This 3D-CDM model was compared to a cohesive zone model (2D-CDM/CZM) to investigate efficacy of aluminum perforation energy prediction, delamination prediction, and computational cost. Accuracy of each model was evaluated using the experimental results. Each showed good agreement with the tests for both the force and velocity histories, as well as the observed damage mechanisms. The 2D-CDM/CZM model was marginally more accurate in capturing both the composite and aluminum behavior — this model averaged error percentages of -11.2% and 10.8% for residual velocity and peak force, respectively. Meanwhile, the 3D-CDM model predictions yielded average error percentages of -35.5% (velocity) and 22.6% (force). However, the 3D-CDM model generally resulted in a decreased computational cost; the average run time was 14% shorter than the 2D-CDM/CZM model and 3x as many timesteps per hour were computed using the same computational resources. In conclusion, new experimental data on the impact and perforation resistance of metal-composite laminates is presented in addition to numerical predictions of the impact behavior.

Carbon fiber↗

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Charge-exchange losses of beam ions in MAST Upgrade: experiments and modelling

Dedicated experiments were performed on MAST Upgrade to study beam-ion losses caused by charge exchange (CX) with edge neutrals. The fuelling was switched from the high-field side to the low-field side mid-discharge. Direct measurements suggest a strong increase in the neutral density around the plasma and a decrease in the beam-ion density, which is qualitatively explained by CX losses. Measurements by a resistive bolometer have suggested particle bombardment during neutral beam injection, providing a unique opportunity to separate CX from other loss mechanisms. To verify and quantify CX losses, the orbit-following code ASCOT, which accounts both for CX neutralization and reionization, was used to simulate beam-particle power loads on the bolometer. Simulations reproduce measured bolometer power loads during high-field-side fuelling, verifying CX losses of approximately 10% of the off-axis beam power. Toroidally symmetric simulations overestimate power loads on the bolometer during low-field-side fuelling, which is explained by toroidal asymmetry in the neutral density distribution, as is demonstrated by toroidally asymmetric simulations. Results suggest significantly higher CX losses during low-field-side fuelling, up to about 50% of off-axis beam power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗