On the numerical solution of Poisson's equation
Numerical integration of Poisson equation for computer simulated plasma
COMPUTER SIMULATION↗
SEARCH · Engineering Papers
Results for “plasma simulations”
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.
Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length
A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
3D full wave fast wave modeling with realistic HHFW antenna geometry and SOL plasma in NSTX-U
This paper reports the significant advancement of our ability to model and to understand how RF waves interact with the SOL plasma, by developing for the full torus 3D scrape-off layer (SOL) plasma simulation together with the antenna and core plasma in the NSTX-U device in HHFW frequency regime for a cold plasma model. The present paper extends the previous modeling (Bertelli et al 2020 AIP Conf. Proc. 2254 030001) to a full 3D device geometry including realistic antenna geometry in order to capture a more realistic experimental situation including the fully 3D wave field effects and the antenna plasma interaction in the SOL plasma, and, at the same time, resulting core wave propagation. Here, the central tool of the present work is the Petra-M code, which is a newly developed state-of-the-art generic electromagnetic simulation tool for modeling RF wave propagation based on MFEM (http://mfem.org), an open source scalable C++ finite element method library.
3D RF modeling↗
Procedures for Including Secondary Electron Emission in Numerical Simulations of Plasma-Insulator Interactions
Previous Monte Carlo simulations provide a data base for properties of secondary electron emission (SEE) from insulators and metals. Incident primary electrons are considered at energies up to 1200 eV. The behavior of secondary electrons is characterized by (1) yield vs. primary energy E(sub p), (2) distribution vs. secondary energy E(sub s), and (3) distribution vs. angle of emission theta. Special attention is paid to the low energy range E(sub p) up to 50 eV, where the number and energy of secondary electrons is limited by the finite band gap of the insulator. For primary energies above 50 eV the SEE yield curve can be conveniently parameterized by a Haffner formula. The energy distribution of secondary electrons is described by an empirical formula with average energy about 8.0 eV. The angular distribution of secondaries is slightly more peaked in the forward direction than the customary cos theta distribution. Empirical formulas and parameters are given for all yield and distribution curves. Procedures and algorithms are described for using these results to find the SEE yield, and then to choose the energy and angle of emergence of each secondary electron. These procedures can readily be incorporated into numerical simulations of plasma-solid surface interactions in low earth orbit.
Beyst, Brian↗
Hybrid-gyrokinetic simulations of low- n toroidal Alfvén eigenmodes using gKPSP
Here, we report a benchmark study of toroidal Alfvén eigenmode (TAE) simulation using the hybrid-gyrokinetic code GyroKinetic Plasma Simulation Program (gKPSP). A simulation capability for energetic particles based on the gyrokinetic δf method has been newly implemented in the gKPSP code. Benchmark simulations have been performed in both circular and realistic tokamak geometries. Good agreement has been found with previously reported results, demonstrating the new capability of the gKPSP code. We have investigated the effects of the distribution function on TAE stability by examining both isotropic and anisotropic slowing-down distributions of energetic particles. The slowing-down distribution produces a higher linear growth rate than a Maxwellian distribution, while a growth rate scan with the anisotropy shows an opposite trend. This can be attributed to competition between Landau damping and the linear drive, which are correlated with the fraction of resonant passing particles and their distribution in phase space.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Massively parallel axisymmetric fluid model for streamer discharges
A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Auroral plasma transport processes in the presence of kV potential structures
We have simulated plasma transport processes in the presence of a quasi-two-dimensional current filament, that generated kV potential structure in the auroral region. The simulation consists of a set of one-dimensional flux tube simulations with different imposed time-dependent, field-aligned currents. The model uses the 16 moment system of equations and simultaneously solves coupled continuity and momentum equations and equations describing the transport along the magnetic field lines of parallel and perpendicular thermal energy and heat flows for each species. The lower end of the simulation is at an altitude of 800 km, in the collisional topside ionosphere, while the upper end is at 10 R(sub E) in the magnetosphere. The plasma consists of hot electrons and protons of magnetospheric origin and low-energy electrons, protons, and oxygen ions of ionospheric origin. The dynamical interaction of the individual current filaments with ionospheric and magnetospheric plasma generates a potential structure in the horizontal direction and kilovolt field-aligned potential drops along the field lines. The side-by-side display exhibits the evolution of the implied potential structure in the horizontial direction. In the presence of this potential structure and parallel electric field ionospheric plasma density is depleted and velocity is reduced, while density enhancement and increased velocity is observed in magnetospheric plasma. The ionospheric and magnetospheric electron temperatures increase below 2 R(sub E) due to magnetic mirror force on converging geomagnetic field lines. The primary cross-field motion produced by the horizontal E field (E x B drift) is perpendicular to both of the significant spatial directions and is thus ignorable in this geometry. The effects of other cross-field drift processes are discussed. The simulation thus provides insight into the dynamical evolution of two-dimensional potential structures driven by an imposed finite width, field-aligned current profile.
Ganguli, Supriya B.↗
Test and Analysis Capabilities of the Space Environment Effects Team at Marshall Space Flight Center
Marshall Space Flight Center has developed world-class space environmental effects testing facilities to simulate the space environment. The combined environmental effects test system exposes temperature-controlled samples to simultaneous protons, high- and low-energy electrons, vacuum ultraviolet (VUV) radiation, and near-ultraviolet (NUV) radiation. Separate chambers for studying the effects of NUV and VUV at elevated temperatures are also available. The Atomic Oxygen Beam Facility exposes samples to atomic oxygen of 5 eV energy to simulate low-Earth orbit (LEO). The LEO space plasma simulators are used to study current collection to biased spacecraft surfaces, arcing from insulators and electrical conductivity of materials. Plasma propulsion techniques are analyzed using the Marshall magnetic mirror system. The micro light gas gun simulates micrometeoroid and space debris impacts. Candidate materials and hardware for spacecraft can be evaluated for durability in the space environment with a variety of analytical techniques. Mass, solar absorptance, infrared emittance, transmission, reflectance, bidirectional reflectance distribution function, and surface morphology characterization can be performed. The data from the space environmental effects testing facilities, combined with analytical results from flight experiments, enable the Environmental Effects Group to determine optimum materials for use on spacecraft.
Finckenor, M. M.↗
Simulations of auroral plasma processes - Electric fields, waves and particles
Plasma processes driven by current sheets of finite thicknesses in an ambient magnetized plasma are studied using a 2 1/2 dimensional particle-in-cell code, and similarities are found between simulated plasma processes and those observed in the auroral plasma. Current sheets are shown to be bounded by large perpendicular electric fields occurring near their edges above the conducting boundary. Shaped potential structures form when the current sheets are narrow, and when the current sheets are wide, potential structures develop a significant parallel potential drop such that the electrons are accelerated upwards. Downward parallel electric fields of variable strength are noted in the downward current region, and double layer formation is seen in both narrow and wide current sheets. High frequency oscillations near the electron plasma frequency and its harmonic are seen, and low frequency waves are observed.
Singh, Nagendra↗
Plasma flow measurements in a simulated low earth orbit plasma
The employment of large, higher power solar arrays for space operation has been considered, taking into account a utilization of high operating voltages. In connection with the consideration of such arrays, attention must be given to the fact that the ambient environment of space contains a tenuous low energy plasma which can interact with the high voltage array causing power 'leakage' and arcing. An investigation has been conducted with the aim to simulate the behavior of such an array in low-earth-orbit (LEO). During the experiments, local concentrations of the 'leakage' current were observed when the panel was at a high voltage. These concentrations could overload or damage a small area of cells in a large string. It was hypothesized that this effect was produced by electrostatic focusing of the particles by the sheath fields. To verify this experimentally, an end-effect Langmuir probe was employed. The obtained results are discussed.
Gabriel, S. B.↗
Coupled radiation transport and plasma physics for ICF simulations [Slides]
Describing high-temperature plasmas requires inclusion of several nonlinear coupled physics, including ion kinetics - for each isotope species present in the plasma, electron kinetics, thermal radiative transfer (TRT), electromagnetic description, and nuclear reactions (fusion burn). Simulating plasmas in ICF hohlraums aims to solve kinetic ions (particles), fluid electrons, energetic electrons, electromagnetics, kinetic radiation (MG DP solver), and laser-plasma interactions.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Conservative discontinuous Galerkin schemes for nonlinear Dougherty–Fokker–Planck collision operators
Here, we present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker–Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the ‘Lenard–Bernstein’ or ‘Dougherty’ operator. Several novel algorithmic innovations, based on the concept of weak equality, are reported. These weak equalities are used to define weak operators that compute primitive moments, and are also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integrations by parts are used to construct the discrete weak form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the kinetic-energy term, $v^2$ is not included in the basis set used in the discretization. A series of benchmark problems is presented and shows that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Koopman wavefunctions and Clebsch variables in Vlasov–Maxwell kinetic theory
Motivated by recent discussions on the possible role of quantum computation in plasma simulations, here, we present different approaches to Koopman's Hilbert-space formulation of classical mechanics in the context of Vlasov–Maxwell kinetic theory. The celebrated Koopman–von Neumann construction is provided with two different Hamiltonian structures: one is canonical and recovers the usual Clebsch representation of the Vlasov density, the other is non-canonical and appears to overcome certain issues emerging in the canonical formalism. Furthermore, the canonical structure is restored for a variant of the Koopman–von Neumann construction that carries a different phase dynamics. Going back to van Hove's prequantum theory, the corresponding Koopman–van Hove equation provides an alternative Clebsch representation which is then coupled to the electromagnetic fields. Finally, the role of gauge transformations in the new context is discussed in detail.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Advancing Multiscale Simulation of Plasma-Surface Interfaces
We report the development of an atomistic-informed, surface-state-dependent predictive model for particle exchange in a carbon-tungsten plasma-surface interface. The predictive model uses machine learning (ML) techniques to learn the energy and angular distributions for particle exchange and rate functions for surface state evolution from molecular dynamics simulations of cumulative bombardment of tungsten by energetic carbon ions. Each predictive component is sensitive to the energy and trajectory of incident plasma species and the surface state. The surface state is represented by a set of surface state descriptors, which were derived from the atomistic surface state for each independent carbon bombardment event. These descriptors are representative of the composition and degree of amorphization of the outermost angstrom of surface material and were chosen to optimize predictive performance for particle exchange at the interface. The distributions for particle exchange (reflection/sputtering) are demonstrated to vary with each surface state descriptor, motivating the development of surface-state-dependent particle exchange models for plasma simulations. The performance of various ML methods was compared, including polynomial quantile regression, artificial neural networks, k-nearest neighbors, and random forest algorithms, with polynomial regression performing the best for interpolation and extrapolation of learned relationships. In addition to the particle exchange model, a neutral network was developed and used to identify data sufficiency throughout surface descriptor space, which will enable real-time feedback during future data production to ensure data is produced where it is most needed, and we provide commentary on improvements to the data production workflow for future endeavors.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Artificial Intelligence/Deep Learning FRNN Software for Prediction & Real-Time Control of DIII-D Plasma Control System (PCS)
This collaborative project integrated an improved version of the Artificial Intelligence/Deep Learning FRNN prediction and control software into the real-time DIII-D PCS (plasma control system). A key AI/DL software challenge is to build a modern high-performance computing (HPC) enabled “synthetic plasma simulator” capable of carrying out HPC-driven real-time plasma control applications. This involves development of a deep learning framework to train the surrogate model for a first-principles-based instability analysis simulator (“SGTC”) derived from the global gyrokinetic code GTC. The role of SGTC is to provide accurate and detailed plasma instability information from a real-time AI-based simulator capability to complement the deep learning prediction and control from experimentally-measured signals, such as ECE Imaging, supplemented by synthetic SGTC-ECEI.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
High-Frequency Density Oscillations from a Plasma Source Used for Simulating Low-Earth Orbit Plasma Environment
We present data from ground-based, vacuum-chamber tests demonstrating the ability to modulate the output of a plasma source capable of producing a low-Earth orbit (LEO) type plasma. We obtained plasma oscillations up to 2.5 kHz impingent on stationary test equipment, which corresponds to meter-level ionospheric structures in LEO. This plasma source is, therefore, suitable for developing scientific instruments that measure the LEO plasma environment, in situ, with meter-level spatial resolution. Measurements were made using a fixed-bias collector and an electrometer sampling at 40 kHz. A mechanical aperture was established at the output of the plasma source via two concentric grids. The outer grid was free to rotate in the azimuthal direction with respect to the fixed inner grid. An identical, alternating hole pattern in the two grids resulted in a variable aperture that cycles through 90 open/close cycles per revolution. The frequency of the plasma oscillations is limited by the mechanism used to spin the grids and the bearing assembly on which the grids rotate. Higher frequencies are obtainable by upgrading the drive mechanism, allowing the possibility of centimeter-level spatial resolution.
McTernan, Jesse K.↗
Simulation of space plasma phenomena
The importance of computer modeling in investigating physical phenomena is discussed, noting that modern computers can follow the time evolution of systems containing many millions of degrees of freedom, all of which are simultaneously interacting with each other. The two types of models used in simulating plasma, that is, particle and fluid, are described. Attention is given to particle models and some recent results obtained concerning the diffuse aurora. The simulations are shown to give clear evidence that a type of nonlinear cyclotron resonance process for the cold particles is responsible for the saturation observed in these calculations. It is also responsible for the cold electron heating and it is found that at saturation the rate of heating of the cold electrons essentially balances the rate at which energy is fed into unstable waves.
Dawson, J. M.↗
Description of space simulation experiments utilizing plasma produced by contact ionization
Space simulation experiments utilizing thermally ionized plasmas
Quinn, R. G.↗