SEARCH · Search NASA
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.
Contents include the following: Radiation belt and plasma model environment. Environment hazards for systems and humans. Need for new models. How models are used. Model requirements. How can space weather community help?
Thermal relaxation rate in one and two dimensional electron plasma models with positive background, using particle in cell simulation
Here, we present a one-dimensional (1-D) model composed of aligned, electrostatically interacting charged disks, conceived to address in a computable model the validity of the Bogoliubov assumption on the decay of particle correlations in the Born–Bogoliubov–Green–Kirkwood–Yvon hierarchy. This assumption is a basic premise of plasma kinetic theory. The disk model exhibits spatially 1-D features at short distances, but retains 3-D features at large distances. Here the collective dynamics of this model plasma is investigated by solving the corresponding Vlasov equation. In addition, the implementation of the model for the numerical validation of the Bogoliubov assumption is formulated.
We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.
Propagation in magneto-ionic media described by plasma models
This year, following the study of liquid-wall options, our edge-plasma modeling has returned to analyzing issues associated with a solid tungsten-wall, Fusion Nuclear Science Facility (FNSF)-size device, primarily using the 2D UEDGE plasma/neutral transport code. The focus is to predict plasma/neutral heat and particle fluxes to various wall components and determine the concentration of injected and wall-sputtered impurities that reach the core plasma region. Our models have been advanced in a number of areas compared to the initial modeling done in the 2017 timeframe: this year we include each individual charge-state of impurity ions instead of assuming a fixed concentration of impurities relative to the deuterium-tritium (DT) fuel plasma, including the interaction between inner and outer divertor legs. Here edge fueling of the impurity and pumping in the private-flux region is incorporated. Increasingly detailed atomic-physics rate files for impurity radiation have been compared, as well as the role of molecules. The impact of spatial variation of the turbulence-driven radial plasma transport is considered to contribute to the understanding of the heat-flux-limit effects of choosing a double-null or single-null divertor configuration. Detailed plasma/neutral wall fluxes, densities, and temperatures have been provided to the ERO group for their more detailed studies of wall sputtering and near-surface impurity transport.
Transport properties and distribution function dynamics in one dimensional plasma models
Over 40 years of ground and spacecraft plasmaspheric measurements have resulted in many statistical descriptions of plasmaspheric properties. In some cases, these properties have been represented as analytical descriptions that are valid for specific regions or conditions. For the most part, what has not been done is to extend regional empirical descriptions or models to the plasmasphere as a whole. In contrast, many related investigations depend on the use of representative plasmaspheric conditions throughout the inner magnetosphere. Wave propagation, involving the transport of energy through the magnetosphere, is strongly affected by thermal plasma density and its composition. Ring current collisional and wave particle losses also strongly depend on these quantities. Plasmaspheric also plays a secondary role in influencing radio signals from the Global Positioning System satellites. The Global Core Plasma Model (GCPM) is an attempt to assimilate previous empirical evidence and regional models for plasmaspheric density into a continuous, smooth model of thermal plasma density in the inner magnetosphere. In that spirit, the International Reference Ionosphere is currently used to complete the low altitude description of density and composition in the model. The models and measurements on which the GCPM is currently based and its relationship to IRI will be discussed.
Etching of semiconductor materials is reliant on plasma properties. Quantities such as ion and neutral fluxes, both in magnitude and in direction, are often determined by reactor geometry (height, radius, position of the coils, etc.) In order to obtain accurate etching profiles, one must also model the plasma as a whole to obtain local fluxes and distributions. We have developed a set of three models that simulates C12 plasmas for etching of silicon, ion and neutral trajectories in the plasma, and feature profile evolution. We have found that the location of the peak in the ion densities in the reactor plays a major role in determining etching uniformity across the wafer. For a stove top coil inductively coupled plasma (ICP), the ion density is peaked at the top of the reactor. This leads to nearly uniform neutral and ion fluxes across the wafer. A side coil configuration causes the ion density to peak near the sidewalls. Ion fluxes are thus greater toward the wall's and decrease toward the center. In addition, the ions bombard the wafer at a slight angle. This angle is sufficient to cause slanted profiles, which is highly undesirable.
We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.
The paper studies charge particle impact as a mechanism for the production of Io's ionosphere. Pioneer 10 thermal plasma measurements and magnetospheric plasma models which explain the observed spatial distribution of neutral hydrogen and sodium atoms in the vicinity of Jupiter's satellite Io imply electron fluxes of about 10 to the 10th/sq cm/s. The fluxes and the temperature (about 100 eV) of this plasma suggest that electron impact ionization is the dominant process in forming the ionosphere of Io. It is found that the surface number density of the neutral species required to match the observed electron density profiles is about 10 to the 9th/cu cm or less. This value is two orders of magnitude lower than previous estimates.
The solution of filtered Vlasov-Maxwell equations, rather than the Vlasov-Maxwell equations themselves, is shown to ameliorate the velocity space filamentation problem in collisionless plasma models. Exact field solutions and filtered velocity distribution functions are obtained without introducing errors. Proper selection of the filter width is demonstrated to inhibit development of velocity space filamentation and, it is conjectured, position space filamentation. The results of sample calculations of both filtered and nonfiltered field solutions illustrate a high degree of agreement between both solutions, with significant savings in computational time and memory requirements with the filtered solutions.
Discrete particle processes in the presence of a strong external magnetic field were investigated. These processes include equations of state and other equilibrium thermodynamic relations, thermal relaxation phenomena, transport properties, and microscopic statistical fluctuations in such quantities as the electric field and the charge density. Results from the equilibrium statistical mechanics of two-dimensional plasmas are discussed, along with nonequilibrium statistical mechanics of the electrostatic guiding-center plasma (a two-dimensional plasma model).
A number of previously published studies apply extended-MHD computations to model macroscopic dynamics of tokamak vertical displacement events (VDEs). The characteristic plasma-wall contact and resistive-wall diffusion imply sensitivity to boundary conditions in a general sense, but which conditions have significant influence depends on how the plasma is modeled. This work examines these dependencies by comparing results of axisymmetric extended-MHD computations with different sets of boundary conditions, plasma thermal-conduction models, and wall resistivity values. The geometry and plasma parameters of the computations represent a modest-sized tokamak. The forced-displacement scenario and computational setup are similar to those used in previous 3D computations. The results show that for a given wall resistivity, the VDE time-scale is most sensitive to variations in the boundary and thermal-conduction parameters that affect electron thermal transport. The electrical conductivity depends on electron temperature, and the dependence on thermal transport stems from its influence on the electrical circuit that includes the open-field halo current. Lastly, conditions that lead to hotter, broader halo regions slow the evolution. Significant sensitivity to the boundary condition on plasma flow-velocity exists when electron thermal conduction is restricted and electron energy loss is convective, which is expected for conditions at the entrance of the magnetic presheath.
Models that treat events associated with quasi-stellar objects as plasma phenomena
The time dependent behavior of a plasma which surrounds a highly biased conducting sphere is considered. This plasma is treated as either a cold two component fluid or as a warm plasma whose time development can be found by solving the Vlasov equation. Both models demonstrate oscillatory behavior, but the electric fields predicted by the models are quantitatively quite different in regions close to the surface of the sphere and very similar otherwise. A broadening of the electron distribution function with time is observed indicating local heating of the plasma near the surface of the sphere.
The Saturn magnetosphere model of Richardson and Sittler (1990) is extended to include the outer magnetosphere. The inner magnetospheric portion of this model is updated based on a recent reanalysis of the plasma data near the Voyager 2 ring plane crossing. The result is an axially symmetric model of the plasma parameters which is designed to provide accurate input for models needing either in situ or line-of-sight data and to be a useful tool for Cassini planning.
Inconsistency of plasma metal interface, discussing /warm/ plasma-filled waveguide