Search NASA⌕ Search

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.

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Implementation of a 9-point stencil in SOLPS-ITER and implications for Alcator C-Mod divertor plasma simulations

The SOLPS-ITER code suite is used worldwide for plasma edge modeling, the interpretation of experiments, as well as for the design of the ITER divertor. The numerical scheme of the plasma solver of the code, B2.5, is based on the assumption of perfectly field-aligned grids, while in practice grids are often strongly distorted to match divertor target shapes. Neglecting these grid distortion leads to qualitatively and quantitatively incorrect results for fluid neutral simulations, and may affect results in cold (detached) divertors even when using kinetic neutral simulations. In this contribution, we present the first results of a newly implemented 9-point stencil in B2.5 to properly handle misaligned grids. The new scheme is then applied to fluid neutral simulations of a well-diagnosed and previously modeled Alcator C-Mod discharge. Results are compared with the original 5-point scheme neglecting grid distortion effects, as well as with simulations including a full kinetic neutral model. We conclude that the 9-point stencil is essential to correctly model the transport of fluid neutrals on distorted grids, and to capture the effects of divertor closure on the fluid neutral behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients in stellarator plasmas

Abstract MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/ .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Simulations of a molecular plasma in collisional-radiative nonequilibrium

A code for the simulation of nonequilibrium plasmas is being developed, with the capability to couple the plasma fluid-dynamics for a single fluid with a collisional-radiative model, where electronic states are treated as separate species. The model allows for non-Boltzmann distribution of the electronic states. Deviations from the Boltzmann distributions are expected to occur in the rapidly ionizing regime behind a strong shock or in the recombining regime during a fast expansion. This additional step in modeling complexity is expected to yield more accurate predictions of the nonequilibrium state and the radiation spectrum and intensity. An attempt at extending the code to molecular plasma flows is presented. The numerical techniques used, the thermochemical model, and the results of some numerical tests are described.

Cambier, Jean-Luc↗
DOE OSTI2024

An implicit particle code with exact energy and charge conservation for studies of dense plasmas in axisymmetric geometries

A collisional particle code based on implicit energy- and charge-conserving methods in axisymmetric geometries is presented. A new particle pusher for axisymmetric systems is introduced that is compatible with exact energy and charge conservation and yields improved accuracy compared to other methods. How to appropriately treat all aspects of the algorithm near the r = 0 axis of symmetry is described in detail. Here, the axisymmetric model is verified by simulating the free expansion of a plasma sphere in 2D cylindrical and 1D spherical geometries. The algorithm's ability to study the dynamic compression of a dense plasma is illustrated by simulating the dynamic Z-pinch in 1D cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Studies of Plasma Flow Past Jupiter's Galilean Satellites

We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus, using our semi-implicit time-dependent 3D MHD code to model the plasma interactions. We have used the same code to model the plasma interaction at Ganymede.

Linker, Jon A.↗
DOE OSTI2023

Verification methods for drift–diffusion reaction models for plasma simulations

Abstract Compared to other computational physics areas such as codes for general computational fluid dynamics, the documentation of verification methods for plasma fluid codes remains under developed. Current analytical solutions for plasma are often highly limited in terms of testing highly coupled physics, due to the harsh assumptions needed to derive even simple plasma equations. This work highlights these limitations, suggesting the method of manufactured solutions (MMSs) as a potential option for future verification efforts. To demonstrate the flexibility of MMS in verifying these highly coupled systems, the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework was utilized. Thanks to the MOOSE framework’s robustness and modularity, as well as to its physics module capabilities and ecosystem applications (i.e. Zapdos and the chemical reaction network) developed for plasma physics modeling and simulation, this report lays the groundwork for a structured method of conducting plasma fluid code verification.

DeChant, Corey (ORCID:0000000293779732)↗
DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Time-dependent models of magnetized pair plasmas

Attention is given to a numerical code developed to study the time evolution of electron-positron plasmas. The code solves in a self-consistent manner kinetic equations describing the effects of Compton scattering, two-photon pair production, pair annihilation, cooling of pairs via Coulomb scattering, e-e bremsstrahlung, and synchrotron radiation. The kinetic equations are derived under the approximation of homogeneous and isotropic particle distributions on the basis of a study by Coppi and Blandford (1993). Both stationary and time-varying output radiation spectra are computed. Good qualitative agreement with previous calculations is found, except where the differences are attributable to the improved treatment of the microphysics.

Coppi, Paolo S.↗
DOE OSTI2024

Validation of prediction capability of operating space for plasma initiation in MAST-U

DYON is a plasma initiation modelling code that solves the differential equation system of the full circuit equations (plasma current, active coil currents and eddy currents in full passive structures) and 0D global energy and particle balance equations (Kim 2022 Nucl. Fusion 62 126012). In order to test the capability of the full electromagnetic plasma initiation model to predict individual discharges in experiments and thus the operating space in the device, a dedicated experimental database was built in MAST-U by scanning the prefilled gas pressure p 0 and the induced loop voltage V loop . In the experimental operating space of p 0 and V loop the lower and the upper limits of p 0 are determined by the plasma breakdown failure and the plasma burn-through failure, respectively. The lower limit of V loop is determined by the plasma burn-through failure. By directly reading the control room data used in each discharge (i.e. currents in the solenoid, poloidal field coils, and toroidal field coils, p 0 , and gas puffing rate), the full electromagnetic DYON consistently predicted the failed breakdown, failed burn-through, and successful plasma initiation discharges in the experimental database, demonstrating its capability to predict the operating space for inductive plasma initiation. The Paschen curve calculated with the effective connection length in MAST-U indicates a much higher p 0 required for plasma breakdown than the experimental data, indicating that individual field line evaluation is necessary to calculate the quantitative requirements for Townsend breakdown. The demonstration in this paper shows that the full electromagnetic DYON could be a useful simulation tool to assess the feasibility of inductive plasma initiation and to optimise operating scenarios in future devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

A general concurrent algorithm for plasma particle-in-cell simulation codes

The general concurrent particle-in-cell (GCPIC) algorithm has been used to implement an electrostatic particle-in-cell code on a 32-node hypercube parallel computer. The GCPIC algorithm decomposes the PIC code by dividing the particle simulation physical domain into subdomains that are equal in number to the number of processors; all subdomains will accordingly possess approximately equal numbers of particles. The portion of the code which updates particle positions and velocities is nearly 100 percent efficient when the number of particles increases linearly with that of hypercube processors.

Liewer, Paulett C.↗
NASA NTRS1987

A theoretical investigation of laser-sustained plasma thruster

A numerical code has been successfully developed for the investigation of thruster performance using a laser-sustained hydrogen plasma as the propellant. The plasma was sustained using a 10.6-micron CO2 laser beam focused at different positions within the thruster. The physical model assumed that plasma is in thermodynamical equilibrium (LTE), and geometric ray tracing was adopted to describe the laser beam. The steady-state, axisymmetric, Navier-Stokes equations coupled with the laser power absorption process have been solved numerically. A pressure based Navier-Stokes numerical solver using body-fitted coordinates was used to calculate the laser-supported rocket flow which includes both subsonic and supersonic flow regions. From the limited parametric study, which did not try to optimize the rocket performance, it was found that better performance was obtained when the laser beam was focused closer to the rocket throat.

Jeng, San-Mou↗
NASA NTRS2000

Simulation and Experimental Measurements of Inductively Coupled CF4 and CF4/Ar Plasmas

The recently developed code SEMS (semiconductor equipment modeling software)is applied to the simulation of CF4 and CF4/Ar inductively coupled plasmas (ICP). This work builds upon the earlier nitrogen, transformer coupled plasma (TCP) SEMS research by demonstrating its accuracy for more complex reactive mixtures, moving closer to the realization of a virtual plasma reactor. Attention is given to the etching of and/or formation of carbonaceous films on the quartz dielectric window and diagnostic aperatures. The simulations are validated through comparisons with experimental measurements using FTIR (Fourier Transform Infrared) and UV absorption spectroscopy for CFx and SiFx neutral radicals, QMS (quadrupole mass spectrometry) for the ions, and Langmuir probe measurements of electron number density and temperature in an ICP GEC reference cell.

Hash, D. B.↗
DOE OSTI2022

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effect of convective transport in edge and scrape-off-layer plasmas of ADITYA-U tokamak

The 2-D edge plasma fluid transport code, UEDGE has been used to simulate the edge and scrape-off layer regions of circular limiter plasmas of ADITYA-U. The computational mesh defining the limiter geometry of ADITYA-U has been introduced in the UEDGE code through an in-house developed routine. It has successfully been integrated with the UEDGE code and the measured radial profile of electron density (n e ) in the edge region is modelled using the simulations. The plasma discharge for the present study has plasma current I P ~ 130 kA, duration ~ 150 ms, chord-averaged electron density ~ 1.3 x 10 19 m -3 , and central electron temperature ~ 270 eV has been observed. The maximum toroidal magnetic field is ~ 1.2 T. It has been found that an inward convective velocity, v conv ~ 1.5 m/s is required in addition to the constant perpendicular diffusion coefficient, D ⟂ ~ 0.2 m 2 /s to adequately match the measured n e profile in typical discharges of ADITYA-U. The value of D ⟂ ~ 0.2 m 2 /s is found to lie between the estimated neoclassical diffusivity and Bohm diffusivity.

Physics - Plasma physics↗
DOE OSTI2022

RF-transpond: A 1D coupled cold plasma wave and plasma transport model for ponderomotive force driven density modification parallel to B

The RF-Transpond code couples a fluid plasma transport solver with a frequency domain cold plasma RF wave solver in a 1D domain parallel to a strong background magnetic field. A ponderomotive force term proportional to parallel gradients in the electric field strength is included in the transport model in order to describe ponderomotive effects in the scrape-off layer (SOL) of fusion plasmas. The transport and wave codes are verified independently and a coupled case corresponding to experimental parameters from the LArge Plasma Device (LAPD) is presented. The density perturbation ratio R n , calculated to describe ponderomotive force driven modifications, is up to 20% for the simulation inputs used.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗