Techniques for testing in a simulated ionospheric plasma.
Test environment and technique permitting ionospheric plasma simulation for spacecraft instruments performance evaluation
Burrowbridge, D. R.↗
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.
Unsupervised discovery of nonlinear plasma physics using differentiable kinetic simulations
Plasma supports collective modes and particle–wave interactions that lead to complex behaviour in, for example, inertial fusion energy applications. While plasma can sometimes be modelled as a charged fluid, a kinetic description is often crucial for studying nonlinear effects in the higher-dimensional momentum–position phase space that describes the full complexity of the plasma dynamics. We create a differentiable solver for the three-dimensional partial-differential equation describing the plasma kinetics and introduce a domain-specific objective function. Using this framework, we perform gradient-based optimization of neural networks that provide forcing function parameters to the differentiable solver given a set of initial conditions. We apply this to an inertial-fusion-relevant configuration and find that the optimization process exploits a novel physical effect.
plasma nonlinear phenomena↗
Space plasma laboratory - Experiment in simulated ionospheric plasma
Several experiments with simulated ionospheric plasma were performed in the Johnson Space Center's large (17 m diameter x 27 m height) vacuum chamber. Two of these, the RF discharge and the electron-beam plasma discharge, are described here. The RF experiment involved measurement of the loading on a 5 m dipole and a 5 m loop antenna caused by the presence of a plasma. The results indicate that when the RF frequency is near the electron gyro frequency a self-sustained plasma discharge is set up and the real part of the complex antenna impedance increases. The beam-plasma experiment was carried out to study a discharge produced by an electron beam with an energy of about 1 keV when the electron beam was injected parallel to the ambient magnetic field. As the beam current was raised to a critical level, a spontaneous transition occurred whereby a luminous halo was formed around the electron beam indicating a sustained discharge.
Konradi, A.↗
Simulations of divertor plasmas with inverse sheaths
The effect of strong electron emission from material surfaces has been proposed to form an “inverse sheath”: a region with a positive potential relative to the near-wall plasma which prevents the flow of ions to the wall [M. D. Campanell, “Negative plasma potential relative to electronemitting surfaces,” Phys. Rev. E. 88, 033103 (2013); M. D. Campanell and M. V. Umansky, “Strongly emitting surfaces unable to float below plasma potential,” Phys. Rev. Lett. 116, 1–5 (2016); M. D. Campanell and G. R. Johnson, “Thermionic cooling of the target plasma to a sub-ev temperature,” Phys. Rev. Lett. 122, 1–5 (2019)]. We assess the viability of this regime in a tokamak device using the 2D edge plasma transport code UEDGE [T. Rognlien et al ., “A fully implicit, time dependent 2-D fluid code for modeling tokamak edge plasmas,” J. Nucl. Mater. 196–198, 347–351 (1992)]. Since the UEDGE code does not consider the sheath region directly, we apply boundary conditions at the divertor targets which emulate the physics of both “standard” and “inverse” sheath regimes [R. Masline et al ., “Influence of the inverse sheath on divertor plasma performance in tokamak edge plasma simulations,” Contrib. Plasma Phys. 60, e201900097 (2020)]. Using these boundary conditions, we perform scoping studies to assess plasma parameters near the target by varying the density at the core-edge interface. We observe a smooth transition in the resultant profiles of plasma parameters for the standard sheath, and a bifurcation across the simulation set for plasmas with an inverse sheath. The cause of this bifurcation is assessed by performing the parameter scan both with and without impurity radiation; we observe that the bifurcation persists in both cases, indicating that this bifurcation is caused by plasma recombination.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Testing of spacecraft systems in a simulated ionospheric plasma
Laboratory simulation of interactions between ionospheric plasma and diverse spacecraft systems
Burrowbridge, D. R., Jr.↗
Testing of spacecraft systems in a simulated ionospheric plasma
Laboratory simulation of interactions between ionospheric plasma and diverse spacecraft systems
Burrowbridge, D. R., Jr.↗
Kinetic Simulation of Inhomogeneous Plasma with a Variable Sized Grid System
Space plasmas are usually inhomogeneous and irregular. Particle in cell (PIC) plasma simulations, however, often consider the plasma to be homogeneous and periodic. Here non-periodic PIC simulations with inhomogeneities in the density, temperature and background magnetic field are considered. Boundary conditions and a variable sized grid system are discussed, along with an application to a magnetospheric plasma physics problem in the auroral zone.
Schriver, David↗
SOLEDGE3X full vessel plasma boundary simulations of ITER non-active phase plasmas
The onset of detachment in the ITER machine is analyzed in this work through the help of 2D-axisymmetric boundary plasma simulations with the SOLEDGE3X-EIRENE code, which features a numerical domain for the plasma solver extending up to the first wall. The plasma boundary is computed in scenarios from the first non-active phase of ITER, in pure H and at 20 MW. This set of simulations is used in two aspects: first, to study the plasma detachment in the divertor, and second, the plasma conditions, fluxes, and beryllium erosion at the first wall. Here, the code results are also compared to those obtained with the well-established SOLPS-ITER code, which includes a plasma numerical domain only covering the main SOL. Results show an increase in the SOL width λ q with increasing density, and a detailed analysis is carried out, for the first time, on each of the different plasma-neutral interactions in the code’s physics model in EIRENE. The gross beryllium erosion rates of first wall panels are estimated from 2D simulations, with the aim of assessing their sensitivity to two parameters: the divertor density regime, and the presence of density shoulders in the far-SOL formed by enhanced perpendicular transport at this location. The erosion contributions from neutrals and ions are considered in each case, and the charge-exchange atoms fluxes and energy distributions are provided, highlighting the two atom populations (cold and charge-exchange).
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Two dimensional PIC simulations of plasma heating by the dissipation of Alfven waves
Two dimensional plasma particle simulations of the evolution of large amplitude circularly polarized Alfven waves propagating parallel to the magnetic field show that the waves decay via both one- and two- dimensional parametric decay instabilities. For parameters studied, one-dimensional processes dominate the simulations, but two-dimensional decay processes, including the recently predicted filamentation instability are also observed. The daughter waves generated by the parametric decay are primarily damped by the ions, leading to ion heating. The parametric decay processes efficiently convert the ordered fluid ion motion in the Alfven wave into ion thermal energy. These processes may be important for the dissipation of Alfven waves in the solar wind, the corona and other space plasma environments. The computations were performed on the Intel Touchstone parallel supercomputer.
Liewer, P. C.↗
Particle simulation of plasmas on the massively parallel processor
Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.
Gledhill, I. M. A.↗
A BOUT++ extension for full annular tokamak edge MHD and turbulence simulations
For tokamak edge plasma simulation, a plasma simulation framework BOUT++ employs a dual coordinate system to simulate moderate-n and high-n plasma instability with reasonable computational cost, where n is the toroidal mode number. This coordinate system however limits the computational domain to the toroidal wedge (full torus divided into N parts in the toroidal direction) for computational efficiency and the use of flute-ordering approximation in the field solver calculating the flow potential from the vorticity which may not be valid for low-n modes. Improving numerical treatment of low-n modes is however indispensable to address simulations of low-n current-driven edge localized mode (ELM), ELM control by resonant magnetic perturbations (RMPs), edge turbulence with RMPs and so on. In this work, BOUT++ is extended to simulate the interplay between $n=0$, low-n and high-n plasma components in a full annular tokamak edge domain through hybrid modeling of the flow potential and the vorticity. Low-n modes of flow potential are calculated in an orthogonal flux surface coordinate and high-n modes in the dual coordinate system separately in Fourier space. Finally, the proposed scheme can capture an interplay between $n=1$ global modes and high-n turbulence during pedestal collapse in a full annular torus domain with a circular cross section.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
3D Electromagnetic Plasma Particle Simulations on a MIMD Parellel Computer
This paper describes computer particle simulation in plasma physics research; also discussed are the results of plasma wave radiation by a moving current source in a magnetized plasma.
N/A↗
Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++
Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Solar Array in Simulated LEO Plasma Environment
Six different types of solar arrays have been tested in large vacuum chambers. The low earth orbit plasma environment was simulated in plasma vacuum chambers, where the parameters could be controlled precisely. Diagnostic equipment included spherical Langmuir probes, mass spectrometer, low-noise CCD camera with optical spectrometer, video camera, very sensitive current probe to measure arc current, and a voltage probe to register variations in a conductor potential. All data (except video) were obtained in digital form that allowed us to study the correlation between external parameters (plasma density, additional capacitance, bias voltage, etc) and arc characteristics (arc rate, arc current pulse width and amplitude, gas species partial pressures, and intensities of spectral lines). Arc inception voltages, arc rates, and current selections are measured for samples with different coverglass materials and thickness, interconnect designs, and cell sizes. It is shown that the array with wrapthrough interconnects have the highest arc threshold and the lowest current collection. Coverglass design with overhang results in decrease of current collection and increase of arc threshold. Doubling coverglass thickness cases the increase in arc inception voltage. Both arc inception voltage and current collection increase significantly with increasing a sample temperature to 80 C. Sustained discharges are initiated between adjacent cells with potential differences of 40 V for the sample with 300 micron coverglass thickness and 60 V for the sample with 150 micron coverglass thickness. Installation of cryogenic pump in large vacuum chamber provided the possibility of considerable outgassing of array surfaces which resulted in significant decrease of arc rate. Arc sites were determined by employing a video-camera, and it is shown that the most probable sites for arc inception are triple-junctions, even though some arcs were initiated in gaps between cells. It is also shown that the arc rate increases with increasing of ion collection current. The analysis of optical spectra (240-800 nm) reveals intensive narrow atomic lines (Ag, H) and wide molecular bands (OH, CH, SiH, SiN) that confirms a complicated mechanism of arc plasma generation. The results obtained seem to be important for the understanding of the arc inception mechanism, which is absolutely essential for progress toward the design of high-voltage solar array for space application.
Vayner, Boris↗
Solar Array in Simulated LEO Plasma Environment
Six different types of solar arrays have been tested in large vacuum chambers. The low Earth orbit plasma environment was simulated in plasma vacuum chambers, where the parameters could be controlled precisely. Diagnostic equipment included spherical Langmuir probes, mass spectrometer, low-noise CCD camera with optical spectrometer, video camera, very sensitive current probe to measure arc current, and a voltage probe to register variations in a conductor potential. All data (except video) were obtained in digital form that allowed us to study the correlation between external parameters (plasma density, additional capacitance, bias voltage, etc) and arc characteristics (arc rate, arc current pulse width and amplitude, gas species partial pressures, and intensities of spectral lines). Arc inception voltages, arc rates, and current collections are measured for samples with different coverglass materials and thickness, interconnect designs, and cell sizes. It is shown that the array with wrapthrough interconnects have the highest arc threshold and the lowest current collection. Coverglass design with overhang results in decrease of current collection and increase of arc threshold. Doubling coverglass thickness causes the increase in arc inception voltage. Both arc inception voltage and current collection increase significantly with increasing a sample temperature to 80 C. Sustained discharges are initiated between adjacent cells with potential differences of 40 V for the sample with 300 m coverglass thickness and 60 V for the sample with 150 m coverglass thickness. Installation of cryogenic pump in large vacuum chamber provided the possibility of considerable outgassing of array surfaces which resulted in significant decrease of arc rate. Arc sites were determined by employing a video-camera, and it is shown that the most probable sites for arc inception are triple-junctions, even though some arcs were initiated in gaps between cells. It is also shown that the arc rate increases with increasing of ion collection current. The analysis of optical spectra (240-800 nm) reveals intensive narrow atomic lines (Ag, H) and wide molecular bands (OH, CH, SiH, SiN) that confirms a complicated mechanism of arc plasma generation. The results obtained seem to be important for the understanding of the arc inception mechanism, which is absolutely essential for progress toward the design of high-voltage solar array for space application.
Vayner, Boris↗
(abstract) 3D Electromagnetic Plasma Particle Simulations
A 3D electromagnetic plasma particle-in-cell code has been developed using the General Concurrent PIC algorithm. The GCPIC algorithm uses a domain decomposition to divide the computation among the processors. Particles must be exchanged between processors as they move. The efficiencies for 1-, 2-, and 3-dimensional partitions of the three dimensional domain are compared, and the algorithm is found to be very efficient even when a large fraction (e.g., 30%) of the particles must be exchanged at every time step. This PIC code has been used to perform simulations of a variety of space plasma physics problems. Results of three applications will be discussed: 1) plasma disturbances induced by moving conducting bodies in a magnetized plasma; 2) plasma plume interactions; and 3) solar wind termination shock.
electromagnetic plasma particles simulations 3D al↗
Numerical simulation of plasma double layers
Numerical simulation results are presented for a plasma double layer, the computer model being a finite one-dimensional particle-in-cell plasma with specified potential difference across the system. A single pulse is formed which crosses the system with constant velocity; this is followed by the formation of a potential drop across a limited region of the plasma. An approximate expression relating the spatial extent of the double layer and the potential drop is presented. Electron and ion beams are generated which tend to lead to instabilities in the upstream and downstream regions.
Joyce, G.↗
Efficacy of the radial pair potential approximation for molecular dynamics simulations of dense plasmas
Macroscopic simulations of dense plasmas rely on detailed microscopic information that can be computationally expensive and is difficult to verify experimentally. In this work, we delineate the accuracy boundary between microscale simulation methods by comparing Kohn–Sham density functional theory molecular dynamics (KS-MD) and radial pair potential molecular dynamics (RPP-MD) for a range of elements, temperature, and density. By extracting the optimal RPP from KS-MD data using force matching, we constrain its functional form and dismiss classes of potentials that assume a constant power law for small interparticle distances. Our results show excellent agreement between RPP-MD and KS-MD for multiple metrics of accuracy at temperatures of only a few electron volts. The use of RPPs offers orders of magnitude decrease in computational cost and indicates that three-body potentials are not required beyond temperatures of a few eV. Due to its efficiency, the validated RPP-MD provides an avenue for reducing errors due to finite-size effects that can be on the order of ∼20%.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗