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 325 records · Page 18

Observation of a magneto-Rayleigh-Taylor instability in magnetically collimated plasma jets

We present the direct experimental observation of the formation of a diamagnetic cavity and magneto-Rayleigh-Taylor (MRT) instability in a β ≈ 1 high energy density plasma. Proton radiography is used to measure the two dimensional path-integrated magnetic field in a laser-produced plasma propagating parallel to a preimposed magnetic field. Flutelike structures, associated with the MRT instability, are observed to grow at the surface of the cavity, with a measured wavelength of 1.2 mm and growth time of 4 ns. These measurements are in good agreement with predictions of three dimensional magnetohydrodynamic simulations using the GORGON code. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrating machine learning interatomic potentials with hybrid reverse Monte Carlo structure refinements in RMCProfile

Structure refinement with reverse Monte Carlo (RMC) is a powerful tool for interpreting experimental diffraction data. To ensure that the under-constrained RMC algorithm yields reasonable results, the hybrid RMC approach applies interatomic potentials to obtain solutions that are both physically sensible and in agreement with experiment. To expand the range of materials that can be studied with hybrid RMC, we have implemented a new interatomic potential constraint in RMCProfile that grants flexibility to apply potentials supported by the Large-scale Atomic/Molecular Massively Parallel Simulator ( LAMMPS ) molecular dynamics code. This includes machine learning interatomic potentials, which provide a pathway to applying hybrid RMC to materials without currently available interatomic potentials. To this end, we present a methodology to use RMC to train machine learning interatomic potentials for hybrid RMC applications.

Cuillier, Paul↗

A high transmission tender X-ray monochromator employing a matched pair of multilayer grating and mirror

It is now well established that multilayer coated gratings may offer high diffraction efficiencies over the tender X-ray range, from about 1 keV to 5 keV, covering the gap between single layer coated grating monochromators and crystal monochromators. Nevertheless, few beamlines in the world are using such gratings in their monochromator. The successful implementation of a multilayer grating monochromator requires producing a matched pair of a multilayer grating and a multilayer mirror, and this matching is not straightforward as it must account for different incidence angles and refraction effects on each element. Here we review the realization of the multilayer grating monochromator of the SIRIUS beamline which has been successfully in service for several years. We show how, by alternating computer simulation with our diffraction code and measurements, we could optimize the monochromator transmission on a very wide energy range. After the grating was coated, it was found that the angle of optimal efficiency versus photon energy was significantly different from what was predicted by a simple conformal model of binary layers. Layer interdiffusion and profile smoothing during the deposition process must be added to the multilayer model to reproduce the measured data. The critical adjustment of the mirror multilayer period is achieved by the lateral translation of the mirror, which was given a small transverse period gradient. The monochromator is thus providing high transmission efficiency in the 1 to 5 keV energy range, more than 30% over 2.5 keV and up to 46% at 4.6 keV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Virtual Time III, Part 3: Throttling and Message Cancellation

This is Part 3 of a trio of papers that unify in a natural way the two historically distinct parallel discrete event synchronization paradigms, optimistic and conservative, combining the best properties of both into a single framework called Unified Virtual Time (UVT). In this part, we survey the synchronization effects that can be achieved by restricting to corner cases the relationships permitted among the control variables, GVT, CVT, TVT, and LVT, which were defined in Part 1. Here we also survey various throttling policies from the literature and describe how they can be implemented in UVT by controlling the value of TVT, including policies that can take advantage of rollback in addition to LP blocking. A significant result is a new category of efficient and higher precision throttling algorithms for optimistic execution that are based on optimistic lookahead, defined in a way that is symmetric to what we now call the conservative lookahead information that is traditionally used for conservative synchronization. Finally, we present a novel algorithm allowing the choice between lazy and aggressive cancellation to be made on a message-by-message basis using either external logic expressed in the model code, or policy code internal to the simulator, or a mixture of both.

throttling↗

Gyrokinetic PIC Study on RMP Affected Neoclassical Transport in Toroidal Plasmas

In magnetically confined fusion plasmas, the breaking of ‘magnetic flux-surfaces’ due to resonant magnetic perturbations (RMPs) can generate magnetic islands and alter field topology to significantly impact plasma confinement and transport. Here, this work investigates the effect of magnetic islands on neoclassical radial energy transport within the core plasma of an analytic circular tokamak using the XGC-S global gyrokinetic particle-in-cell code. Findings from our simulations revealed substantial enhancements in electron neoclassical radial energy diffusivity in and around the islands, in addition to a newly observed two-peak structure at the O/X-points and outer island boundary in the electron diffusivity profile.

RMP↗

Data driven methods to recognize patterns in EIC weak-strong simulation

Beam-Beam simulations are currently being studied in preparation for future EIC experiments to study beam-beam effects and, in turn, maximize luminosity. Weak-strong methods are studied for single-particle dynamics during collision. 1 million macro-particles for 1 million turns are typically tracked, corresponding to only 10 seconds in the EIC. The goal of this study is to predict beam properties over the scale of hours. A potential solution focuses on using data-driven methods such as machine learning methods to analyze and extend the insights of the beam properties such as long-term nonlinear effects. This would aid in long-term predictions where results would be more efficiently acquired than a typical tracking simulation. Some limitations such as inaccurate predictions and spatial complexity are also discussed. These methods can then be applied to strong-strong simulations in the future studies.

Accelerator Physics↗

FEM Analysis of Hybrid LTS/HTS Cos-Theta Dipole Magnet

After the recent results on Bi-2212 superconductive magnets realized and tested in canted cosine-theta and solenoid designs, respectively, at Lawrence Berkeley National Laboratory (LBNL) and NHMFL, the first Bi-2212 stress-managed cosine-theta insert magnet is in the assembly phase at Fermilab. This insert will be part of the first hybrid cosine-theta magnet made of Nb3Sn outer layers within the US-MDP effort to reach 20 T bore field. This paper reports the analytical analysis of the cosine-theta Nb3Sn/Bi-2212 hybrid magnet. As an introduction, it shows the validation of the finalized HTS coils' design optimized using the ROXIE code. Subsequently, it reports the parameters, logic, and implementation method of the 2D electromagnetic and mechanical FEM analysis of the LTS/HTS hybrid magnet. Results from the ANSYS detailed model made with sub-modeling geometry are compared with the homogeneous model implemented in the past. Moreover, a Python code was implemented to simulate the current degradation due to stresses in the detail-modeled conductor areas. The current degradation is introduced in the simulation dynamics for both conductors as an iteration process during the energization load step only after applying pre-stress and thermal loads. The numerical and graphic results of the 2D cosine-theta LTS/HTS magnet cross-section will be described and analyzed.

D'Agliano, A.↗

Geant4 Event Biasing and Fast Simulation

Geant4 offers advanced event biasing techniques to significantly accelerate simulations involving rare events. Various biasing methods, such as leading particle selection, cross-section biasing, radioactive decay enhancement, and bremsstrahlung splitting, enable efficient event sampling, though they require careful handling. Additionally, Geant4 provides a Fast Simulation Interface, allowing the replacement of standard processes in specific region and for selected particles, enabling faster execution or external code integration. Applications of fast simulation include electromagnetic shower modeling in calorimeters, machine learning inference, and offloading tasks to specialized hardware like GPUs, making Geant4 a powerful tool for computationally demanding simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Solenoidal Fields from Permanent Magnets for a Polarized Electron Preinjector

This technical note describes the design of a compact focusing solenoid constructed from permanent bar magnets, intended for the low-energy section of a spin-polarized electron injector. Typically, quadrupole triplets provide appropriate focusing for the low-energy spin rotator section; however, their large footprint precludes the insertion of sufficient beam diagnostics. The primary motivation of this study is to provide beam focusing while preserving the vertical spin orientation of the electron beam. We propose a solenoidal field composed of multiple permanent magnets arranged to naturally yield a zero longitudinal magnetic field integral ($\int B_zdz = 0$). This design serves as a compact, passive, and non-coupling alternative to traditional quadrupole triplets. We present the magnet arrangement, field calculation methodology, spin tracking simulations using space charge codes, and calculations demonstrating that the required field strength is achievable.

43 PARTICLE ACCELERATORS↗

Experimental Study of Alfvén Wave Reflection from an Alfvén-speed Gradient Relevant to the Solar Coronal Holes

Abstract We report the first experimental detection of a reflected Alfvén wave from an Alfvén-speed gradient under conditions similar to those in coronal holes. The experiments were conducted in the Large Plasma Device at the University of California, Los Angeles. We present the experimentally measured dependence of the coefficient of reflection versus the wave inhomogeneity parameter, i.e., the ratio of the wavelength of the incident wave to the length scale of the gradient. Two-fluid simulations using the Gkeyll code qualitatively agree with and support the experimental findings. Our experimental results support models of wave heating that rely on wave reflection at low heights from a smooth Alfvén-speed gradient to drive turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. II. General Relativistic Module

Black hole (BH) environments often host plasmas that are fully collisionless or contain intrinsically collisionless regions, including relativistic jets and coronae, where particle energization is ubiquitous. Capturing the physics of these systems requires numerical methods capable of modeling relativistic, magnetized, collisionless plasmas in strong gravitational fields. In this work, we introduce the general relativistic module for Entity—the first open-source, coordinate-agnostic, performance-portable, particle-in-cell code. The code enables fast axisymmetric simulations of collisionless plasmas around BHs on any modern high-performance computing architecture (both GPUs and CPUs).

Galishnikova, Alisa [Flatiron Institute, New York,↗

Thermal-Mechanic Modeling of Fusion Components Using The MOOSE Framework

With fusion energy rapidly developing, there is a great need for an open source tool that can rapidly model different fusion tokamak designs. The open-source codes Multiphysics Object Oriented Simulation Environment (MOOSE), Tritium Migration Analysis Program Version 8 (TMAP8) , and Fusion ENergy Integrated Multiphysi-X (FENIX) are being used to demonstrate thermal-mechanic modeling of fusion components. In this study, an analysis is done on the fusion nuclear science facility (FNSF) blanket, and ITER toroidal magnet and diverter mono block. This study will lay the groundwork for fusion component analysis in the MOOSE framework and FENIX program.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗

First operation and validation of simulations for the divertor cryo-vacuum pump in Wendelstein 7-X

Ten cryo-vacuum pumps (CVPs) were installed in the subdivertor region of each island divertor in the stellarator Wendelstein 7-X (W7-X) and operated for the first time during the recently completed plasma campaign OP2.1. A pumping speed of 70 ± 1 $\frac{m^3}{s}$ was measured during dedicated tests with known hydrogen gas injection. Based on a conductance model, the estimated pumping speed ranges from 86-93 $\frac{m^3}{s}$ for different sticking coefficients between 0.6 and 0.8. After completion of the initial tests the CVPs were operated successfully throughout the campaign, with regeneration performed once a week. Neutral gas pressures in the subdivertor in the range of 10 −4 mbar are well within the molecular flow regime and limit the particle exhaust capabilities of the CVPs. Simulations of the neutral gas pressure in the three-dimensional complex geometry of the subdivertor were performed using the DIVGAS code based on the direct simulation Monte Carlo method and a model implemented in the steady-state thermal package in ANSYS, which are in agreement with the measured values during plasma operation.

Cryo-vacuum pumping↗

Hydrodynamic expansion and near-infrared absorption of x-ray heated aluminum plasmas

We use x-ray pulses from dense argon plasmas at the Z Machine (Sandia National Laboratories) to generate hypersonic aluminum plasmas akin to material ejecta during proposed planetary defense missions, fusion reactor wall excursions, and other high-energy density processes. Near-infrared absorption is used to diagnose the controlled expansion of the plasmas through cylindrical cavities following their generation from x-ray heating of solid aluminum 7075 alloy. The data are compared to multidimensional radiation hydrodynamics simulations utilizing the ALEGRA multiphysics code, accounting for the dynamics of radiation scattering, material phase change, plasma expansion, thermal re-irradiation, and interactions with the cavity and with the infrared beams. To allow for accurate simulation, density functional theory is used to apply the Hagen–Rubens relation for the far-infrared and is adjoined with a detailed configuration accounting model using the Propaceos code, producing opacities spanning 10 −1 –10 4 eV photon energy for aluminum 7075 alloy, and in comparison with pure aluminum. The model is found to agree with experimental data in the higher-fluence regime when the Hagen–Rubens relation is applied. The ejected material, which is observed to travel up to 55 km/s, is comprised of a strongly ionized, non-LTE plasma front at ∼10 eV temperature followed by a weakly ionized LTE gas at higher density. The present findings lend some confidence to the broad-range equation of state and infrared opacity models for weakly ionized aluminum plasmas while demonstrating an approach to their future refinement, with potential application to astrophysical plasmas and other extreme processes.

Adiabatic process↗

Consist v0.1.0

A Python library for provenance tracking, intelligent caching, and data virtualization in scientific simulation workflows. It automatically records code, configuration, and input data to skip redundant computations and enables querying results across many runs without manual bookkeeping. Designed to support multi-model simulation workflows like the BEAM CORE toolset at LBL, but designed to be extensible to a wide range of research workflows. Combines lineage tracking features as provided by OpenLineage with deterministic hashing like SnakeMake, and adds powerful analysis tools on model outputs.

Needell, Zachary [Lawrence Berkeley National Labor↗

Acceleration of the particle-in-cell code Osiris with graphics processing units

Fully relativistic particle-in-cell (PIC) simulations are crucial for advancing our knowledge of plasma physics. Modern supercomputers based on graphics processing units (GPUs) offer the potential to perform PIC simulations of unprecedented scale, but require robust and feature-rich codes that can fully leverage their computational resources. In this work, this demand is addressed by adding GPU acceleration to the PIC code Osiris. An overview of the algorithm, which features a CUDA extension to the underlying Fortran architecture, is given. Detailed performance benchmarks for thermal plasmas are presented, which demonstrate excellent weak scaling on NERSC's Perlmutter supercomputer and high levels of absolute performance. The robustness of the code to model a variety of physical systems is demonstrated via simulations of Weibel filamentation and laser-wakefield acceleration run with dynamic load balancing. Finally, measurements and analysis of energy consumption are provided that indicate that the GPU algorithm is up to ~14 times faster and ~7 times more energy efficient than the optimized CPU algorithm on a node-to-node basis. The described development addresses the PIC simulation community's computational demands both by contributing a robust and performant GPU-accelerated PIC code and by providing insight into efficient use of GPU hardware.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗