Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “plasma simulation”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Geometric electrostatic particle-in-cell algorithm on unstructured meshes

We present a geometric particle-in-cell (PIC) algorithm on unstructured meshes for studying electrostatic perturbations with frequency lower than electron gyrofrequency in magnetized plasmas. In this method, ions are treated as fully kinetic particles and electrons are described by the adiabatic response. The PIC method is derived from a discrete variational principle on unstructured meshes. To preserve the geometric structure of the system, the discrete variational principle requires that the electric field is interpolated using Whitney 1-forms, the charge is deposited using Whitney 0-forms and the electric field is computed by discrete exterior calculus. The algorithm has been applied to study the ion Bernstein wave (IBW) in two-dimensional magnetized plasmas. The simulated dispersion relations of the IBW in a rectangular region agree well with theoretical results. In a two-dimensional circular region with fixed boundary condition, the spectrum and eigenmode structures of the IBW are obtained from simulations. We compare the energy conservation property of the geometric PIC algorithm derived from the discrete variational principle with that of previous PIC methods on unstructured meshes. The comparison shows that the new PIC algorithm significantly improves the energy conservation property.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of collisionality on turbulence in the edge of tokamak plasma using 3D global simulations

Collisionality is one of the key parameters in determining turbulent transport in the plasma edge, regulating phenomena such as 'shoulder formation', separation of scale lengths in the scrape-off layer (SOL), turbulence damping and zonal flow dynamics. Understanding its role is therefore of primary importance for future reactors like ITER. Getting reliable predictions and a better characterization of plasma flow properties when varying collisionality remains, however, a critical challenge for the simulations. This paper focuses on the impact of varying collisionality in a non-isothermal three-dimensional fluid model of the plasma edge. A high field side limited configuration encompassing open and closed magnetic field lines with parameters typical of a medium-sized tokamak is considered. The present model can consistently account for the variations of collisionality and its impact on both the parallel resistivity η ∥ and the ion and electron parallel thermal conductivities χ ∥e,i . Details on mean flow and turbulence properties are given. Changing collisionality leads to significant changes in the flow properties both on the mean and fluctuating quantities. In particular, lowering collisionality decreases the size of coherent structures, the fluctuation levels of turbulence, and steepens the density and temperature equilibrium profiles around the separatrix leading to a global reduction of the turbulent transport. The SOL width is observed to increase with collisionality, eventually resulting in the disappearance of the scale lengths separation between near and far SOL, consistently with previous experimental observations. At low collisionality, where the presence of narrow feature is well-established, a contribution of heat conduction increases up to compete with heat convection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Nonlinear coupling of whistler waves to oblique electrostatic turbulence enabled by cold plasma

Kinetic simulations and theory demonstrate that whistler waves can excite oblique, short-wavelength fluctuations through secondary drift instabilities if a population of sufficiently cold plasma is present. The excited modes lead to heating of the cold populations and damping of the primary whistler waves. The instability threshold depends on the density and temperature of the cold population and can be relatively small if the temperature of the cold population is sufficiently low. Furthermore, this mechanism may thus play a significant role in controlling amplitude of whistlers in the regions of the Earth's magnetosphere where cold background plasma of sufficient density is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A parallel-kinetic-perpendicular-moment model for magnetised plasmas

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

astrophysical plasmas↗
DOE OSTI2026

Role of ion acoustic instability in magnetic reconnection

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Initial position optimization in molecular dynamics simulations for a Coulomb system

A new algorithm for molecular dynamics (MD) simulations is developed to optimize plasma particle distributions at given initial temperatures. By combining velocity scaling and reassignment, the method effectively eliminates the initial rise and oscillation in temperatures observed with randomly distributed positions. These rises and oscillations are undesired numerical artifacts observed in conventional plasma MD simulations, arising from unoptimized particle positions. The algorithm demonstrates temperature relaxation without initial rises or oscillations, as well as precise flow velocity relaxation, enabling accurate measurement of relaxation times. The code is accelerated using graphics processing units for parallel processing, enhancing the study of plasma dynamics. The proposed method for distributing physically valid particles in MD simulations enables accurate studies of intrinsic collision processes in plasmas, including the dynamics of strongly coupled plasmas, plasma–wave interactions, and transport phenomena in magnetized plasmas. The paper concludes with a discussion of potential applications and future enhancements to the algorithm.

Jo, Jawon (ORCID:0009000924193285)↗
DOE OSTI2021

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ion temperature gradient (ITG), trapped electron mode (TEM) and electron temperature gradient (ETG) modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. Finally, it is also confirmed that thermodiffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

3D full wave fast wave modeling with realistic antenna geometry and SOL plasma

This paper reports the significant advancement of our ability to model and to understand how RF waves interact with the SOL plasma, by simulating the full torus 3D SOL plasma together with the antenna and core plasma. We introduce and use a recently developed and open source code, Petra-M, which was constructed on the scalable MFEM C++ finite element library, and performed 3D full wave simulations in the HHFW regime for both NSTX-U and LAPD plasmas. A first full wave simulation for a full 3D torus including a realistic antenna geometry and SOL plasma region for NSTX-U is presented. A scan of the antenna phasing shows a strong interaction between FWs and the SOL plasma for lower antenna phasing, which is consistent with previous NSTX HHFW observations. The effect of the 3D wave field on the fast ion population from NBI beams in NSTX-U is also discussed by using the 3D field obtained from the Petra-M simulations in the full-orbit following particle SPIRAL code. On LAPD, 3D full wave simulations of a new HHFW 4-strap antenna recently installed by TAE Technologies on LAPD are performed showing a qualitative agreement with experimental data.

Bertelli, Nicola↗
DOE OSTI2021

Analysis of Design Alternatives of Actively Cooled RF Window for MPEX

The Materials Plasma Exposure eXperiment (MPEX) has been designed as a linear plasma divertor simulator in order to address plasma material interaction (PMI) science for next-generation fusion devices. It will have the capability to test neutron irradiated samples with plasma fluxes of greater than 10 24 m -2 s -1 . It is expected to operate steady state for up to 10 6 s to consider PMI affects through reactor end of life. The conceptual design of MPEX was completed in 2019, with preliminary design having begun in 2020. The plasma source for MPEX is a helicon antenna, where the energized helical antenna sits outside of the vacuum in order to minimize impurities in the plasma. It is expected to receive up to 200 kW of continuous power, and so the antenna and the window must be actively cooled. The water-cooled copper antenna has been operated at full power on the Proto-MPEX device (which is a test facility to demonstrate the plasma source and heating systems). The water-cooled window, however, is a novel component that must meet numerous competing requirements. It requires a low dielectric loss to allow the Radio Frequency (RF) power to create the plasma within the vacuum boundary. It must be structurally robust to handle the significant heat flux from the plasma and any heat from dielectric coupling. It must be compatible with the coolant (preferably water). It requires a vacuum seal that minimizes impurities into the plasma and does not compromise the structural integrity of the window. Two window designs have been tested. Results from these tests, where temperatures are measured and heat fluxes inferred from infrared camera data, have been correlated with thermal-structural simulations. When these simulations are extrapolated to the full power steady-state heat fluxes that are expected in MPEX, the designs do not appear to have the necessary structural robustness. This study explores design alternatives for the MPEX helicon antenna window, presents analysis results for several of the alternatives, and shows a viable solution that satisfies the requirements for MPEX operation.

Tipton Jr, Joseph↗
DOE OSTI2022

Eulerian simulations of electrostatic waves in plasmas with a single sign of charge

An Eulerian, numerical simulation is used to model the launching of plasma waves in a non-neutral plasma that is confined in a Penning–Malmberg trap. The waves are launched by applying an oscillating potential to an electrically isolated sector at one end of the conducting cylinder that bounds the confinement region and are received by another electrically isolated sector at the other end of the cylinder. The launching of both Trivelpiece–Gould waves and electron acoustic waves is investigated. Adopting a stratagem, the simulation captures essential features of the finite length plasma, while retaining the numerical advantages of a simulation employing periodic spatial boundary conditions. As a benchmark test of the simulation, the results for launched Trivelpiece–Gould waves of small amplitude are successfully compared to a linearized analytic solution for these fluctuations.

Physics↗
DOE OSTI2025

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

SciDAC ISEP: Integrated Simulation of Energetic Particles in Burning Plasmas

The objective of the SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP) is to improve physics understanding of energetic particle (EP) confinement and EP interactions with burning thermal plasmas through large-scale simulations. The ISEP center will develop a multiscale and multiphysics ISEP framework for a predictive capability of EP physics and deliver an EP module incorporating both first-principles simulations and high fidelity reduced transport models to the fusion whole device modeling (WDM) project. The ISEP framework will enable us to perform long time, global kinetic simulations of EP physics in burning plasmas, by utilizing the full power of the next generation supercomputers. Our research and development activities will build on fruitful collaborations with computer scientists and applied mathematicians to offer enabling technologies for performance scalability, portability, solvers, coupling for integration with the fusion WDM project, and long-term preservation of data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Electron-Only Reconnection in Plasma Turbulence

Hybrid-Vlasov–Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast to the usual picture, where a full hierarchy of structure sizes is generally observed. These current sheets are shown to be the sites of electron-only reconnection (e-rec), in which the usual electron exhausts are unaccompanied by ion outflows and which are in qualitative agreement with those recently observed by MMS in the Earth's turbulent magnetosheath, downstream of the bow shock. Some features of the e-rec phenomenology are shown to be consistent with an electron magnetohydrodynamic description. Simulations suggest that this regime of collisionless reconnection may be found in turbulent systems where plasma processes, such as micro-instabilities and/or shocks, overpower the more customary turbulent cascade by directly injecting energy close to the ion-kinetic scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Power-transfer and fixed-point analysis of sawtooth simulations of a current-carrying stellarator

Power-transfer and fixed-point analysis of previous NIMROD simulations (Roberds et al. , Phys. Plasmas , vol. 23, issue 9, 2016, 092513) improved the understanding of the effect of 3D (non-axisymmetric equilibrium) magnetic fields on sawtooth oscillations in the Compact Toroidal Hybrid (CTH) experiment. Computing the locations of order-1 fixed points, their Greene's residues, and local values for the rotational transform results in a description of CTH sawteeth consistent with Kadomtsev's model. A power-transfer analysis quantifies the distribution of energy among toroidal Fourier modes and their nonlinear interactions. The Lorentz power transfer drives sawtooth growth, and it is unambiguously interpreted as the flow of energy from toroidal mode $n'$ to mode $n$ , catalysed by $\boldsymbol {B}_{n-n'}$ . It has been reported previously that the CTH sawtooth frequency increases with the 3D field strength. This is attributed to an increased growth rate of the internal kink that drives sawtooth oscillations. Here, 3D fields remove energy from the kink, eliminating the possibility that these fields are an additional energy source that drives growth. Instead, 3D fields catalyse energy transfer from large-to-small scales, where magnetic reconnection is stronger. It is proposed that this energy transfer increases the reconnection rate at small scales, which is consistent with the increased growth rate observed at higher 3D field strengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An implicit particle code with exact energy and charge conservation for electromagnetic studies of dense plasmas

A collisional particle code based on implicit energy- and charge-conserving methods is presented. A modified version of the particle-suppressed Jacobian-Free Newton-Krylov method that can enhance the solver efficiency is introduced. Mathematically, it is shown that this new approach can be viewed as a fixed-point iteration method for the particle positions. The model can exactly conserve global energy and local charge and can efficiently use time steps larger than the plasma period. In conclusion, the algorithm's ability to simulate dense plasmas accurately and efficiently is quantified by simulating the dynamic compression of a plasma slab via a magnetic piston in 1D planar geometry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Analysis and design of fast flow liquid Li divertor for fusion nuclear science facility (FNSF) using coupled plasma boundary and LM MHD/heat transfer codes *

The SOLPS-ITER code is utilized to analyze the boundary plasma associated with a fast-flow lithium (Li) divertor configuration in the fusion nuclear science facility (FNSF) tokamak and identify operational regimes with acceptable divertor and core conditions. Plasma transport from the SOLPS-ITER code has been coupled with a liquid metal (LM) MHD/heat transfer code to model a Li open-surface divertor design and assess its impact on the scrape-off-layer (SOL) and core plasma performance. Simulations with only Neon (Ne) impurity seeding have been conducted to evaluate its impact on meeting FNSF design demands for the divertor and upstream plasma parameters. Simulation results indicate that Ne seeding significantly mitigates divertor heat flux but potentially reduces both upstream electron and main ion density due to fuel dilution. The combined application of Ne seeding and deuterium (D 2 ) puffing is required to satisfy the FNSF design requirements on upstream density ($n^{OMP}_{e,sep}$~1× 10 20 m -3 ) and divertor energy flux ($q^{Odiv}_{\bot,max}$<10 MW m -2 ). D 2 puffing plays a role in counteracting upstream density drops and augmenting energy and momentum losses through atomic and molecular processes. The inlet Li flow velocity is systematically varied across a wide range to identify acceptable flows and corresponding LM surface temperatures. This comprehensive analysis identifies the acceptable Li flow parameters, LM surface temperature, and emitted Li fluxes necessary to meet the major design constraints. The emitted Li fluxes exhibit minimal impact on the main plasma at surface temperatures up to approximately ~525 °C, corresponding emitted Li fluxes of up to φ Li ~2X10 23 atoms s -1 . Uncertainties in the Li emission processes from the surface are also investigated, primarily influencing Li loss in the lower surface temperature range (<525 °C), with simulation results indicating a minor impact on the divertor and upstream plasma. Conversely, evaporation predominantly drives the Li loss processes at higher surface temperature ranges (>525 °C), contaminating both the divertor and upstream plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas

Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗