On Rostoker's test particle theorem
Rostoker particle theorem for multispecies plasma in framework of Klimontovich formalism
MATHEMATICAL MODEL↗
SEARCH · Engineering Papers
Results for “plasma models”
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.
Quantum Information for Fusion Energy Sciences (Final Technical Report)
The simulation of plasma dynamics is a critical area of Fusion Energy Sciences (FES) due to it’s usefulness in predicting, controlling, and confining plasmas in the context of potential fusion reactors. The simulation of plasmas is a computationally difficult problem in both classical and quantum physics, motivating investigation into the potential of quantum computers to simulate these systems. This project took several concrete steps towards this goal by developing tools for improving the control, characterization, and calibration of quantum gates on a superconducting quantum computer, developing error suppression and mitigation tools to reduce errors on the quantum computer, and utilizing these advancements to simulate reduced models of plasma dynamics on the quantum computer. In order to efficiently simulate plasma physics, an optimal control method which synthesizes, directly at the pulse level, any quantum gate on qubit and qutrit systems was developed. Using four superconducting transmon quantum processors at Rigetti and LLNL, it was demonstrated that any arbitrary quantum gate on qubits and qutrits could be implemented with high fidelity, leading to a significantly reduced length of a gate sequence. A problem of interest in FES is the nonlinear optical process of laser pulse compression within a plasma. Since quantum physics is linear, simulating nonlinear operations is not naturally feasible on a quantum computer, however it is possible to simulated a quantized version of the nonlinear process. A quantization approach to convert nonlinear wave-wave interaction problems to Hamiltonian simulation problems was developed and demonstrated using two qubits on a Rigetti device. In this experiment, a number of error suppression and mitigation techniques were investigated to determine how best to utilize the finite quantum resources. This study provides an example of how plasma problems may be solved on near-term, noisy quantum computing platforms and identified a promising set of techniques. Building on the insights of these experiments, the investigation turned to linear electron-plasma wave physics. A connection was identified between a local one-dimensional lattice spin model and linear wave phenomena, allowing a plasma physics problem to be efficiently mapped to the quantum computer. In this framework, reflection and transmission of plasma waves at a sharp boundary was studied, as well as the propagation of waves through an inhomogeneous plasma medium. In addition to the suite of error suppression and mitigation techniques developed, this experiment introduced the use of a digital-analog gate scheme designed to efficiently simulate the plasma Hamiltonian. With hardware available at the conclusion of the project, simulation at the scale of 9 qubits and 15 timesteps (60 entangling layers) was achieved.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Dynamics of laser-generated magnetic fields using long laser pulses
In this work, we report on the experimental investigation of magnetic field generation with a half-loop gold sheet coil driven by long-duration (10 ns) and high-power (0.5 TW) laser pulses. The amplitude of the magnetic field was characterized experimentally using proton deflectometry. The field rises rapidly in the first 1 ns of laser irradiation, and then increases slowly and continuously up to 10 ns during further laser irradiation. The transient dynamics of current shape were investigated with a two-dimensional (2D) numerical simulation that included Ohmic heating of the coil and the resultant change of electrical resistivity determined by the coil material temperature. The numerical simulations show rapid heating at the coil edges by current initially localized at the edges. This current density then diffuses to the central part of the sheet coil in a way that depends both on normal current diffusion as well as temporal changes of the coil resistance induced by the Ohmic heating. The measured temporal evolution of the magnetic field is compared with a model that determines a solution to the coil current and voltage that is consistent with a plasma diode model of the drive region and a 2D simulation of current diffusion and dynamic resistance due to Ohmic heating in the laser coil.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids
We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Monte Carlo Simulation of Nonlinear Radiation Induced Plasmas
A Monte Carlo simulation model for radiation induced plasmas with nonlinear properties due to recombination was, employing a piecewise linearized predict-correct iterative technique. Several important variance reduction techniques were developed and incorporated into the model, including an antithetic variates technique. This approach is especially efficient for plasma systems with inhomogeneous media, multidimensions, and irregular boundaries. The Monte Carlo code developed has been applied to the determination of the electron energy distribution function and related parameters for a noble gas plasma created by alpha-particle irradiation. The characteristics of the radiation induced plasma involved are given.
Wang, B. S.↗
Uncovering turbulent plasma dynamics via deep learning from partial observations
One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a novel multi-network physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. This technique presents a novel paradigm for the advanced design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids
This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.
graphics processing units↗
Modeling X-Ray Photoionized Plasmas: Ion Storage Ring Measurements of Low Temperature Dielectronic Recombination Rate Coefficients for L-Shell Iron
Iron L-shell ions (Fe XVII to Fe XXIV) play an important role in determining the line emission and thermal and ionization structures of photoionized gases. Existing uncertainties in the theoretical low temperature dielectronic recombination (DR) rate coefficients for these ions significantly affects our ability to model and interpret observations of photoionized plasmas. To help address this issue, we have initiated a laboratory program to produce reliable low temperature DR rates. Here, we present some of our recent results and discuss some of their astrophysical implications.
Savin, D. W.↗
Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas
Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge
Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Microtearing instabilities and electron thermal transport in low and high collisionality NSTX discharges
Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. Here, the dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate. In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Uncovering turbulent plasma dynamics via deep learning from partial observations
One of the most intensely studied aspects of magnetic confinement fusion is edge plasma turbulence which is critical to reactor performance and operation. Drift-reduced Braginskii two-fluid theory has for decades been widely applied to model boundary plasmas with varying success. Towards better understanding edge turbulence in both theory and experiment, we demonstrate that a novel multi-network physics-informed deep learning framework constrained by partial differential equations can accurately learn turbulent fields consistent with the two-fluid theory from partial observations of electron pressure which is not otherwise possible using conventional equilibrium models. Furthermore, this technique presents a novel paradigm for the advanced design of plasma diagnostics and validation of magnetized plasma turbulence theories in challenging thermonuclear environments.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas
We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Pharmacokinetic Modeling of Intranasal Scopolamine in Plasma Saliva and Urine
An intranasal gel formulation of scopolamine (INSCOP) was developed for the treatment of Space Motion Sickness. The bioavailability and pharmacokinetics (PK) were evaluated under the Food and Drug Administration guidelines for clinical trials with an Investigative New Drug (IND). The aim of this project was to develop a PK model that can predict the relationship between plasma, saliva and urinary scopolamine concentrations using data collected from the IND clinical trial with INSCOP.
Wu, L.↗
A survey of the effects of magnetic fields, resistivity, viscosity and thermal conduction on the Rayleigh–Taylor instability
The Rayleigh–Taylor instability (RTI) is ubiquitous in nature and has been investigated in great detail in a number of theoretical, numerical, and experimental studies relevant to laboratory and space plasmas. Here, the RTI can significantly be altered by the presence of magnetic fields, resistivity, viscosity, and thermal conduction. This paper presents a brief survey and summarizes initial findings of the effect of a plasma RTI modified by these effects.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Density profiles in stellarators: an overview of particle transport, fuelling and profile shaping studies at TJ-II
We provide an overview of activities carried out at the TJ-II stellarator for improving our understanding of- and developing plasma physics models for particle density profiles in stellarators. Namely, we report on recent progress in turbulent particle transport simulation, validation of pellet deposition models, density profile shaping for performance control and new experimental techniques for edge turbulence and plasma-neutral interaction.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Radiative dominated cooling of the flare corona and transition region
Models in which radiation dominates cooling flare loops are investigated. The radiative models are found to predict a differential emission measure (Q) proportional to T to the (l+1) power, where l measures the dependence of the radiative loss coefficient on temperature, lamda (T) approximately T to the (-l) power. It is concluded that the radiative models are incapable of explaining the observed temperature dependence of Q for flare coronal and transitional plasma. The models suggest that large mass motions (velocities of the order of the sound speed) may be required.
Antiochos, S. K.↗
Magnetic Lens For Plasma Engine
Low-field electromagnet coils placed downstream of plasma engine, polarized oppositely to higher-field but smaller radius coil in nozzle of engine, reduces divergence of plasma jet, thereby increasing efficiency of engine. Concept tested by computer simulation based on simplified mathematical model of plasma, engine, and coils.
Sercel, Joel C.↗