SEARCH · Search NASA
Results for “Plasma 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.
Causal explicit algorithm for heat conduction in a plasma
Hyperbolic heat conduction extends standard Spitzer-Harm heat conduction by including a term proportional to the time derivative of the heat flux. The new term arises from a kinetic derivation of the heat flux that includes higher order corrections. Here we present a causal explicit numerical algorithm for solving the nonlinear hyperbolic heat conduction equation in an unmagnetized plasma. The maximum stable timestep for the causal explicit algorithm scales linearly with the cell size, owing to the hyperbolic nature of the problem. This is in contrast to the quadratic scaling of the maximum stable timestep with the cell size for the parabolic forward time centered space algorithm. The favorable scaling of the timestep with the cell size enables a practical explicit implementation of heat conduction in high-performance massively parallel plasma codes. In particular, we have implemented the causal explicit algorithm in the laser plasma interaction code pF3D. We verify the CE algorithm and analyze its convergence rate by simulating a harmonic mode, which has an analytic solution within the context of the HHC model. We also compare simulations using the CE algorithm to those using the forward time centered space algorithm on a pair of test problems: evolution in time of a Gaussian temperature perturbation in a uniform plasma and heat transport in the presence of inverse bremsstrahlung heating by a Gaussian laser speckle.
Demonstration of object location, classification, and characterization by developed deep learning dust ablation trail analysis code package using plasma jets
Based on deep learning, a Dust Ablation Trail Analysis (DATA) code package is developed to detect dust ablation trails in tokamaks, which is intended to analyze a large amount data of tokamak dusts. To validate and benchmark the DATA code package, 2440 plasma jet images are exploited for the training and test of the deep learning DATA code package, since plasma jets resemble the shape and size of dust ablation clouds in tokamaks. After being trained by 1920 plasma jet images, the DATA code package is able to locate 100% plasma jets, classify plasma jets with the accuracy of >99.9%, and output image skeleton information for classified plasma jets. The DATA code package trained by the plasma jet images is also used to analyze the dust ablation trails captured in the Experimental Advanced Superconducting (EAST) tokamak with the satisfactory performance, further verifying its applicability in the fusion dust ablation investigation. Based on its excellent performance presented here, it is demonstrated that our DATA code package is able to automatically identify and analyze dust ablation trails in tokamaks, which can be used for further detailed investigations, such as the three-dimensional reconstruction of dusts and their ablation trails.
Title: Benchmarking of High-performance Low-temperature Plasma Particle-in-cell Codes for Modeling Magnetized Plasmas
Understandinrelationships between magnetic field and plasma conditions
Time-dependent modeling of coupled plasma-wall dynamics
Effects of dynamic wall outgassing on divertor plasma during an ELM-like heat pulse are modelled using newly coupled edge plasma code UEDGE and all reaction-diffusion code FACE. Different divertor regimes are simulated depending on the degree of the plasma detachment. It is shown that hydrogen outgassing from the divertor target plate induced by ELM pulse in the semi-detached regime can trigger deep plasma detachment. The induced detachment does not spontaneously reverse after the pulse due to ineffective implantation of low energy plasma hydrogen in the target material. Furthermore, we demonstrate that plasma neutrals play essential role in heat transport in a detached divertor. The modelled hydrogen dynamics in the wall material demonstrates that the heat pulse leads to enhanced hydrogen transport toward the wall bulk.
Verification and benchmarking relativistic electron beam transport through a background gas
It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. Additionally, we also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.
Modelling of edge plasma dynamics with active wall boundary conditions
A self-consistent 2D model is presented for transport in boundary plasma and plasma-facing material walls. Plasma dynamics in the domain is represented by a 2D collisional plasma fluid model in the edge-plasma code UEDGE, and transport of hydrogen and heat in the wall is represented by a system of reaction–diffusion equations in the 1D wall code FACE. To account for variation of parameters along the wall, in the coupled model multiple instances of the FACE code run in parallel. Here, the coupled model provides a tool for investigating a range of dynamic plasma–material interactions phenomena in 2D. For demonstration of its capability, one application of particular interest is the role of active wall in tokamak strike point sweeping proposed for mitigation of divertor heat loads. In the present study, the coupled calculations are applied to investigation of the impact of heat and hydrogen transport in the material wall on the divertor plasma and target heat load during sweeping of the target strike point for parameters of a high-power tokamak.
Characterization of sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements with 3D hybrid particle-in-cell code
Hard x-rays produced by intense laser-produced fast electrons interacting with solids are a vital source for producing radiographs of high-density objects and implosion cores for inertial confinement fusion. Accurate calculation of hard x-ray sources requires a three-dimensional (3D) simulation geometry that fully models the electron transport dynamics, including electron recirculation and the generation of absolute photon yields. To date, 3D simulations of laser-produced bremsstrahlung photons over tens of picoseconds and code benchmarking have not been performed definitively. In this study, we characterize sub-picosecond laser-produced fast electrons by modeling angularly resolved bremsstrahlung measurements for refluxing and non-refluxing targets using the 3D hybrid particle-in-cell (PIC), Large Scale Plasma code. Bremsstrahlung radiation and escaped electron data were obtained by focusing a 50-TW Leopard laser (15 J, 0.35 ps, 2 × 10 19 W/cm 2 ) on a 100-μm-thick Cu foil and a Cu with a large plastic backing (Cu–CH target). Data for both the Cu and Cu–CH targets were reproduced for simulations with a given set of electron parameters. Comparison of the simulations revealed that the hard x-ray emission from the Cu target was significantly longer in duration than that from the Cu–CH target. Here, the benchmarked hybrid PIC code could prove to be a powerful tool in the design and optimization of time- and angular-dependent bremsstrahlung sources for flash x-ray and gamma-ray radiography.
Formation of edge pressure pedestal and reversed magnetic shear due to toroidal rotation in a tokamak equilibrium
Toroidal rotation is well known to play significant roles in the edge transport and L-H transition dynamics of tokamaks. Our recent calculation finds that a sufficiently strong localized toroidal rotation can directly bring out the formation of edge pressure pedestal with reversed magnetic shear that is reminiscent of an H-mode plasma, purely through the effects of toroidal rotation on the tokamak MHD equilibrium itself. In particular, the enhanced edge toroidal rotation enables a substantial peaking of the parallel current profile near edge in higher β regimes, which leads to the flattening or reversal of the local q (safety factor) profile. Here the formation of pressure pedestal along with the reversed magnetic shear region is shown to be the natural outcome of the MHD tokamak equilibrium in a self-consistent response to the presence of a localized toroidal rotation typically observed in H-mode or QH-mode.
The magnetic structure of long-wavelength magnetohydrodynamic modes in current-carrying stellarator plasmas
Magnetic field fluctuations are observed in current-carrying stellarator plasmas when the rotational transform is close to a rational value at the edge of the plasma. At low plasma pressure, these fluctuations are associated with perturbed currents parallel to the equilibrium magnetic field lines. A model for these magnetohydrodynamic modes in a low-β, three-dimensional stellarator equilibria has been developed. A set of helical current filaments are constrained to mimic the structure of magnetic field lines on rational surfaces derived from three-dimensional (3D) equilibrium reconstructions. Transformation to straight field line coordinates then allows fitting of the poloidal magnetic sensor data to a single harmonic function, which fixes the modeled toroidal mode structure via the field line flow geometry. The developed procedure accurately captures phase and amplitude variation for m/n = 3/2, 3/1, and 4/1 modes in the 3D equilibria of the compact toroidal hybrid experiment.
Platform for 100 s Mbar equation of state measurements on the National Ignition Facility
Equation of state (EOS) measurements in the 100 s Mbar range are needed to underwrite models employed in the simulation of high energy density plasmas. To this end, a platform has been developed for fielding on the National Ignition Facility, capable of producing high-quality impedance match EOS data, wherein a planar, high-pressure, steady shock is driven into a sample package, and sample and reference standard shock velocities are measured. This platform, dubbed planar high pressure, or PHP, was fielded with an initial proof-of-concept shot in January 2023. The first PHP shot, aiming to study gold, demonstrated a pressure close to 400 Mbar, two orders of magnitude higher than previously reported gold EOS data.
Collisional-Radiative Nonequilibrium and Precursor Effects in a Nitrogen Shock Wave
Improvements to a plasma code with a Collisional-Radiative (CR) non-equilibrium model are made, allowing for a more accurate description of the physical processes. The code allows for non-Boltzmann distributions of the electronic excited states by convecting separately each excited state, as a pseudo-specie. Each molecular state has also its own vibrational temperature, while a global rotational temperature is assumed. The free electron temperature is different from those of the excited states, and the electron heat conduction is also included. The CR model also uses a unique coupling between chemistry and vibrational energy (C-V coupling), which is fully coherent, and has the property of establishing thermal equilibrium as well as chemical equilibrium, on its own. We have also included a coupling between electronic excitations and vibrational energy (X-V coupling), which can have a strong influence on the vibrational temperature of some states. The recent improvements include the multi- temperature dependence of the chemical rates for associative ionization, as well as the estimation of the internal energies transferred during this process. Additionally, the distribution of energy into different translational modes (electron and heavy particles) is now correctly modeled. This provides a very rapid heating mechanism for the free electrons, since it is found that the electrons are generated with an average thermal energy of the same order as the heavy particle translational energy. This effect was observed by Gorelov et al in a recent paper, and lead to pronounced peaks in electron temperature immediately behind the shock. We will attempt ro reproduce this phenomenon. The last modification concerns the inclusion of the radiative terms into the calculations, thus enabling us to observe the effect of radiative losses and radiation transport. Preliminary tests have shown that the radiative losses are not negligible, i.e. the shock velocity drops when the radiative emission is included. In addition, we have observed some precursor excitation and ionization phenomena. These may be important at high shock velocities, and this will be investigated in detail. Precursor phenomena were also observed (in air) by Gorelov et al, and our results (for pure nitrogen) will be qualitatively compared with theirs.
XMM-Newton and Chandra Observations of the Galaxy Group NGC 5044: Evidence for Limited Multiphase Hot Gas - 1
Using new XMM and Chandra observations, we present an analysis of the temperature structure of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot gas is favored over single-phase or cooling flow (M = 4.5 +/- 0.2 solar mass/yr) models within the central approx.30 kpc. Alternatively, the data can be fitted equally well if the temperature within each spherical shell varies continuously from approx.T(sub h) to T(sub c) approx. T(sub h)/2, but no lower. The high spatial resolution of the Chandra data allows us to determine that the temperature excursion T(sub h) approaches T(sub c) required in each shell exceeds the temperature range between the boundaries of the same shell in the best-fitting single-phase model. This is strong evidence for a multiphase gas having a limited temperature range. We do not find any evidence that azimuthal temperature variations within each annulus on the sky can account for the range in temperatures within each shell. We provide a detailed investigation of the systematic errors on the derived spectral models considering the effects of calibration, plasma codes, bandwidth, variable NH, and background rate. We find that the RGS gratings and the EPIC and ACIS CCDs give fully consistent results when the same models are fitted over the same energy ranges for each instrument. The cooler component of the 2T model has a temperature (T(sub c) approx. 0.7 keV) similar to the kinetic temperature of the stars. The hot phase has a temperature (T(sub h) approx. 1.4 keV) characteristic of the virial temperature of the solar mass halo expected in the NGC 5044 group. However, in view of the morphological disturbances and X-ray holes visible in the Chandra image within R approx. equals 10 kpc, bubbles of gas heated to approx.T(sub h) in this region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be associated with the evolution of the cold front near R approx. 50 kpc, as suggested by the sharp edge in the EPIC images. Subject headings: cooling flows - galaxies: formation - galaxies: halos - intergalactic medium - X-rays: galaxies: clusters
Doubly-Shadowed Regions in Lunar Polar Craters: Hydrogen Accumulation in the Presence of Recursive Plasma Wakes
Permanently shadowed regions (PSRs) of the Moon have been identified as unique environments of extreme cold and comprise a natural cold trap for sequestering volatiles [Paige et al. 2010]. The diverse chemical composition of the LCROSS impact plume provided evidence for a volatile-rich and chemically-complex PSR environment [Cola prete et al. 2010, Schultz et al. 2010]. Additionally, the polar electrostatic environment is highly complex, with the possibility of strong, localized electric fields that divert solar wind ions directly into polar cold traps [Farrell et al. 2010, Zimmerman et al. 2011]. Thus, regional plasma physics processes couple directly with volatile sequestration. In the present work, kinetic simulations show that recursive plasma wake structure arises in the presence of step-like topographic features (Le. doubly-shadowed craters). Combining the plasma code with a numerical sputtering model demonstrates that solar wind protons can be either a hydrogen source via implantation or a volatile loss mechanism via sputtering, depending on properties of the regolith and solar wind. The present model provides a novel theoretical pathway toward understanding the lunar surface/solar wind physical and chemical interactions for complex topography near the poles.
A New Method to Constrain Supernova Fractions Using X-ray Observations of Clusters of Galaxies
Supernova (SN) explosions enrich the intracluster medium (ICM) both by creating and dispersing metals. We introduce a method to measure the number of SNe and relative contribution of Type Ia supernovae (SNe Ia) and core-collapse supernovae (SNe cc) by directly fitting X-ray spectral observations. The method has been implemented as an XSPEC model called snapec. snapec utilizes a single-temperature thermal plasma code (apec) to model the spectral emission based on metal abundances calculated using the latest SN yields from SN Ia and SN cc explosion models. This approach provides a self-consistent single set of uncertainties on the total number of SN explosions and relative fraction of SN types in the ICM over the cluster lifetime by directly allowing these parameters to be determined by SN yields provided by simulations. We apply our approach to XMM-Newton European Photon Imaging Camera (EPIC), Reflection Grating Spectrometer (RGS), and 200 ks simulated Astro-H observations of a cooling flow cluster, A3112.We find that various sets of SN yields present in the literature produce an acceptable fit to the EPIC and RGS spectra of A3112. We infer that 30.3% plus or minus 5.4% to 37.1% plus or minus 7.1% of the total SN explosions are SNe Ia, and the total number of SN explosions required to create the observed metals is in the range of (1.06 plus or minus 0.34) x 10(exp 9), to (1.28 plus or minus 0.43) x 10(exp 9), fromsnapec fits to RGS spectra. These values may be compared to the enrichment expected based on well-established empirically measured SN rates per star formed. The proportions of SNe Ia and SNe cc inferred to have enriched the ICM in the inner 52 kiloparsecs of A3112 is consistent with these specific rates, if one applies a correction for the metals locked up in stars. At the same time, the inferred level of SN enrichment corresponds to a star-to-gas mass ratio that is several times greater than the 10% estimated globally for clusters in the A3112 mass range.
Performance-Portable GPU Acceleration of the EFIT Tokamak Plasma Equilibrium Reconstruction Code
This paper presents the steps followed to GPU-offload parts of the core solver of EFIT-AI, an equilibrium reconstruction code suitable for tokamak experiments and burning plasmas. For this work, we will focus on the fitting procedure that consists of a Grad–Shafranov (GS) equation inverse solver that calculates equilibrium reconstructions on a grid. We will show profiling results of the original code (CPU-baseline), as well as the directives used to GPU-offload the most time-consuming function, initially to compare OpenACC and OpenMP on NVIDIA and AMD GPUs and later on to assess OpenMP performance portability on NVIDIA, AMD and Intel GPUs. We will make a performance comparison for different spatial grid sizes and show the speedup achieved on NVIDIA A100 (Perlmutter-NERSC), AMD MI250X (Frontier-OLCF) and Intel PVC GPUs (Sunspot-ALCF). Finally, we will draw some conclusions and recommendations to achieve high-performance portability for an equilibrium reconstruction code on the new HPC architectures
Applications of numerical codes to space plasma problems
Solar wind, earth's bowshock, and magnetospheric convection and substorms were investigated. Topics discussed include computational physics, multifluid codes, ionospheric irregularities, and modeling laser plasmas.
Simulation of divertor heat flux width on EAST by BOUT++ transport code
The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.