Engineering Papers⌕ Search

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.

At least 703 records · Page 39

Interpretative transport modeling of the WEST boundary plasma: main plasma and light impurities

Understanding impurity transport in tokamak plasmas is crucial to control radiative losses and material migration in future magnetic fusion reactors. In this work we deploy the SolEdge2D-EIRENE code to model the boundary plasma in a WEST discharge, satisfactorily reproducing measurements of both upstream and divertor plasma conditions. Additionally, the spatial distribution of oxygen, studied here as a representative light impurity, is compared to vacuum ultraviolet spectroscopy measurements acquired with an oscillating line of sight. The simulation captures a key feature of the experiment, namely a factor of sime2 higher oxygen brightness in the inner divertor region compared to the outer one. This spatial asymmetry in oxygen concentration is interpreted by analyzing the balance of friction forces and thermal gradient forces that the light impurity exchanges with the main plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temporal features of the outflow of heavy ionospheric ions in response to a high altitude plasma cavity

Using a hydrodynamic model for the plasma, it is demonstrated that the auroral plasma cavity is capable of drawing an appreciably large flux of oxygen ions, which are normally gravitationally bound. This escape mechanism of O(+) does not involve any additional heating or acceleration of the plasma in the ionosphere. The temporal evolution of the outflow shows that it starts near the cavity and penetrates into the ionosphere with a flux front moving down with the ion-acoustic speed. A steady outflow is reached in a few hours for a cavity at a height of one R(e) if the cavity is maintained. This time reduces to less than 1 hour for cavity heights less than 3000 km. During the transient state the outflux can appreciably exceed that in the steady state. In view of the horizontal convection, the transient outflow is of special significance. As the cavity descends, the ionospheric flux of O(+) increases; for the lower edge of the cavity in the altitude range 3000-10,000 km, the steady flux is found to be in the range 10 to the 7th to 10 to the 9th. An analysis for the steady state shows that outflow is controlled by the density scale height in the cavity.

Singh, Nagendra↗

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid computational model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, along with the associated values of net charge, electric field, and electric potential were calculated throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the accelerating grids) to examine the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown. The results from this simulation are compared to existing literature on ion-ion plasma thrusters.

Sankaran, K.↗

Computational Investigation of the Near-Field Plasma Plume in Ion-Ion Propulsion

A two-fluid numerical model of plasma flows was developed to investigate the plume of an ion-ion propulsion system. The densities of positive and negative ions, and the associated values of net charge, electric field, and electric potential were calculated as a function of time throughout the domain. The computational domain was chosen to be large enough (25 thruster diameters downstream of the exit plane) to allow for examining the neutralization of the plume. The resulting plasma electric potential and charge neutrality at the downstream end of the domain are shown and they indicate that it is possible to alternatively accelerate oppositely charged ions without the need for an electron-emitting neutralizer and without facing any electric potential hills that could cause stagnation. However, compared to existing literature on ion-ion plasma thrusters, the results from this simulation predict a longer length-scale for voltage decay.

Sankaran, K.↗

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling an Iodine Hall Thruster Plume in the Iodine Satellite (ISAT)

An iodine-operated 200-W Hall thruster plume has been simulated using a hybrid-PIC model to predict the spacecraft surface-plume interaction for spacecraft integration purposes. For validation of the model, the plasma potential, electron temperature, ion current flux, and ion number density of xenon propellant were compared with available measurement data at the nominal operating condition. To simulate iodine plasma, various collision cross sections were found and used in the model. While time-varying atomic iodine species (i.e., I, I+, I2+) information is provided by HPHall simulation at the discharge channel exit, the molecular iodine species (i.e., I2, I2+) are introduced as Maxwellian particles at the channel exit. Simulation results show that xenon and iodine plasma plumes appear to be very similar under the assumptions of the model. Assuming a sticking coefficient of unity, iodine deposition rate is estimated.

Numerical Simulation↗

Modeling an Iodine Hall Thruster Plume in the Iodine Satellite (ISAT)

An iodine-operated 200-W Hall thruster plume has been simulated using a hybrid-PIC model to predict the spacecraft surface-plume interaction for spacecraft integration purposes. For validation of the model, the plasma potential, electron temperature, ion current flux, and ion number density of xenon propellant were compared with available measurement data at the nominal operating condition. To simulate iodine plasma, various collision cross sections were found and used in the model. While time-varying atomic iodine species (i.e., I, I+, I2+) information is provided by HP Hall simulation at the discharge channel exit, the molecular iodine species (i.e., I2, I2+) are introduced as Maxwellian particles at the channel exit. Simulation results show that xenon and iodine plasma plumes appear to be very similar under the assumptions of the model. Assuming a sticking coefficient of unity, iodine deposition rate is estimated.

Iodine↗

Plasma Beta Above a Solar Active Region: Rethinking the Paradigm

In this paper, we present a model of the plasma beta above an active region and discuss its consequences in terms of coronal magnetic field modeling. The beta-plasma model is representative and derived from a collection of sources. The resulting beta variation with height is used to emphasize the assumption that the magnetic pressure dominates over the plasma pressure must be carefully considered depending on what part of the solar atmosphere is being considered. This paper points out (1) that the paradigm that the coronal magnetic field can be constructed from a force-free magnetic field must be used in the correct context, since the forcefree region is sandwiched between two regions which have beta greater than 1, (2) that the chromospheric MgIICIV magnetic measurements occur near the beta-minimum, and (3) that, moving from the photosphere upwards, beta can return to 1 at relatively low coronal heights, e.g. R approximately 1.2R(sub)s.

Gary, G. Allen↗

Low-Power Operation and Plasma Characterization of a Qualification Model SPT-140 Hall Thruster for NASA Science Missions

Hall thruster systems based on commercial product lines can potentially lead to lower cost electric propulsion (EP) systems for deep space science missions. A 4.5-kW SPT-140 Hall thruster presently under qualification testing by SSL leverages the substantial heritage of the SPT-100 being flown on Russian and US commercial satellites. The Jet Propulsion Laboratory is exploring the use of commercial EP systems, including the SPT-140, for deep space science missions, and initiated a program to evaluate the SPT-140 in the areas of low power operation and thruster operating life. A qualification model SPT-140 designated QM002 was evaluated for operation and plasma properties along channel centerline, from 4.5 kW to 0.8 kW. Additional testing was performed on a development model SPT-140 designated DM4 to evaluate operation with a Moog proportional flow control valve (PFCV). The PFCV was commanded by an SSL engineering model PPU-140 Power Processing Unit (PPU). Performance measurements on QM002 at 0.8 kW discharge power were 50 mN of thrust at a total specific impulse of 1250 s, a total thruster efficiency of 0.38, and discharge current oscillations of under 3% of the mean current. Steady-state operation at 0.8 kW was demonstrated during a 27 h firing. The SPT-140 DM4 was operated in closed-loop control of the discharge current with the PFCV and PPU over discharge power levels of 0.8-4.5 kW. QM002 and DM4 test data indicate that the SPT-140 design is a viable candidate for NASA missions requiring power throttling down to low thruster input power.

Garner, Charles E.↗

Hermes-3

Hermes-3 is a plasma simulation model built on BOUT++, developed as part of research studies of the edge and divertor regions of magnetically confined fusion plasmas such as tokamaks. The main purpose of this model is multi-species simulation of plasmas containing a mixture of deuterium, tritium, helium and other plasma and neutral gas species.

Dudson, BenjaminD↗

Empirical modeling of the quiet time nightside magnetosphere

Empirical modeling of plasma pressure and magnetic field for the quiet time nightside magnetosphere is investigated. Two models are constructed for this study. One model, referred to here as T89R, is basically the magnetic field model of Tsyganenko (1989) but is modified by the addition of an inner eastward ring current at a radial distance of approximately 3 R(sub E) as suggested by observation. The other is a combination of the T89R model and the long version of the magnetic field model of Tsyganenko (1987) such that the former dominates the magnetic field in the inner magnetosphere, whereas the latter prevails in the distant tail. The distribution of plasma pressure, which is required to balance the magnetic force for each of these two field models, is computed along the tail axis in the midnight meridian. The occurrence of pressure anisotropy in the inner magnetospheric region is also taken into account by determining an empirical fit to the observed plasma pressure anisotropy. This effort is the first attempt to obtain the plasma pressure distribution in force equilibrium with magnetic stresses from an empirical field model with the inclusion of pressure anisotropy. The inclusion of pressure anisotropy alters the plasma pressure by as much as a factor of approximately 3 in the inner magnetosphere. The deduced plasma pressure profile along the tail axis is found to be in good agreement with the observed quiet time plasma pressure for geocentric distances between approximately 2 and approximately 35 R(sub E).

Lui, A. T. Y.↗

Empirical modeling of the quiet time nightside magnetosphere

Empirical modeling of plasma pressure and magnetic field for the quiet time nightside magnetosphere is investigated. Two models are constructed for this study. One model, referred to here as T89R, is basically the magnetic field model of Tsyganenko (1989) but is modified by the addition of an inner eastward ring current at a radial distance of approximately 3 RE as suggested by observation. The other is a combination of the T89R model and the long version of the magnetic field model of Tsyganenko (1987) such that the former dominates the magnetic field in the inner magnetosphere while the latter prevails in the distant tail. The distribution of plasma pressure which is required to balance the magnetic force for each of these two field models is computed along the tail axis in the midnight meridian. The occurrence of pressure anisotropy in the inner magnetospheric region is also taken into account by determining an empirical fit to the observed plasma pressure anisotropy. This represents the first effort to obtain the plasma pressure distribution in force equilibrium with magnetic stresses from an empirical field model with the inclusion of pressure anisotropy. The inclusion of pressure anisotropy alters the plasma pressure by as much as a factor of approximately 3 in the inner magnetosphere. The deduced plasma pressure profile along the tail axis is found to be in good agreement with the observed quiet time plasma pressure for geocentric distances between approximately 2 and approximately 35 RE.

Lui, A. T. Y.↗

Towards Predictive Plasma Science and Engineering through Revolutionary Multi-Scale Algorithms and Models (Final Report)

This report describes the high-level accomplishments from the Plasma Science and Engineering Grand Challenge LDRD at Sandia National Laboratories. The Laboratory has a need to demonstrate predictive capabilities to model plasma phenomena in order to rapidly accelerate engineering development in several mission areas. The purpose of this Grand Challenge LDRD was to advance the fundamental models, methods, and algorithms along with supporting electrode science foundation to enable a revolutionary shift towards predictive plasma engineering design principles. This project integrated the SNL knowledge base in computer science, plasma physics, materials science, applied mathematics, and relevant application engineering to establish new cross-laboratory collaborations on these topics. As an initial exemplar, this project focused efforts on improving multi-scale modeling capabilities that are utilized to predict the electrical power delivery on large-scale pulsed power accelerators. Specifically, this LDRD was structured into three primary research thrusts that, when integrated, enable complex simulations of these devices: (1) the exploration of multi-scale models describing the desorption of contaminants from pulsed power electrodes, (2) the development of improved algorithms and code technologies to treat the multi-physics phenomena required to predict device performance, and (3) the creation of a rigorous verification and validation infrastructure to evaluate the codes and models across a range of challenge problems. These components were integrated into initial demonstrations of the largest simulations of multi-level vacuum power flow completed to-date, executed on the leading HPC computing machines available in the NNSA complex today. These preliminary studies indicate relevant pulsed power engineering design simulations can now be completed in (of order) several days, a significant improvement over pre-LDRD levels of performance.

42 ENGINEERING↗

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Particle control via cryopumping and its impact on the edge plasma profiles of Alcator C-Mod

At the high n e proposed for high-field fusion reactors, it is uncertain whether ionization, as opposed to plasma transport, will be most influential in determining n e at the pedestal and separatrix. A database of Alcator C-Mod discharges is analyzed to evaluate the impact of source modification via cryopumping. The database contains similarly-shaped H-modes at fixed I P = 0.8 MA and B t = 5.4 T, spanning a large range in P net and ionization. Measurements from an edge Thomson scattering system are combined with those from a midplane-viewing Ly α camera to evaluate changes to n e and T e in response to changes in ionization rates, S ion ∙ $n^{sep}_e$ and $T^{ped}_e$ are found to be most sensitive to changes to $S^{sep}_{ion}$, as opposed to $n^{ped}_e$ and $T^{sep}_e$. Dimensionless quantities, namely α MHD and v*, are found to regulate attainable pedestal values. Select discharges at different values of P net and in different pumping configurations are analyzed further using SOLPS-ITER. It is determined that changes to plasma transport coefficients are required to self-consistently model both plasma and neutral edge dynamics. Pumping is found to modify the poloidal distribution of atomic neutral density, n 0 , along the separatrix, increasing n 0 at the active X-point. Opaqueness to neutrals from high n e in the divertor is found to play a role in mediating neutral penetration lengths and hence, the poloidal distribution of neutrals along the separatrix. Pumped discharges thus require a larger particle diffusion coefficient than that inferred purely from 1D experimental profiles at the outer midplane.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Models of the plasmaspheric thermal plasma distribution

Current understanding of the thermal plasma in the atmosphere and its coupling to the ionosphere is reviewed. Existing models appear adequate to explain the gross behavior of the cold thermal plasma, but there remain some vexing problems. Notably, (1) why does the density in flux tubes appear to saturate at lower values than are predicted theoretically, (2) what causes the sunset peak in measured Te, and (3) why does the equatorial plasmapause signature differ in latitude from the ionosphere signatures. The more difficult problem of what happens during the early stages of refilling after a magnetic storm, when the high altitude plasma is likely to be supersonic and collisionless, has received much attention, but the results are not definite. A number of papers have dealt with the interaction of supersonic counterstreaming fluxes and there are now models that can handle the transition from supersonic to subsonic flows although the transition from a collisionless to a collision-dominated plasma remains difficult to deal with.

Richards, P. G.↗

Multi‐Scale Model‐Informed Deep Learning for Plasma‐Nanoparticle Interaction

The Overarching Goal of this proposed research is to understand and quantitively determine the interactions between non-thermal plasma (hot electrons, reactive radicals, vibrationally excited species) and surface reactions on influencing the activity and selectivity of the desired reactions via developing multi-scale model informed deep learning algorithm. Investigating non-thermal plasma-surface interaction is feasible due to the low bulk temperature in the discharge region. To investigate the role of plasma-nanoparticle interaction on enhancing the reaction kinetics, we will focus on ammonia cracking to generate clean hydrogen over earth-abundant, non-critical metallic nanoparticles, which is of great significance for decarbonization. We hypothesize that (1) reactive radicals interacting with surface reaction species via Eley–Rideal mechanism will significantly lower the energetics of the potential rate-limiting step of nitrogen formation; (2) the surface will be charged heterogeneously under non-thermal plasma conditions and the charged site will lower the energetics of ammonia cracking through Langmuir– Hinshelwood mechanism; (3) vibrationally excited ammonia will further promote the initial N-H bond cleavage. To access the hypothesis, we will (1) reveal the surface charge effects on tunning the reaction energetics via interpretable, physics-informed deep learning accelerated density functional theory (DFT) calculations; (2) determine the reactive radicals interacting with surface reaction species on tuning the reaction energetics via DFT; (3) reveal the surface charge effects on tunning the reaction energetics via DFT and deep learning models, (4) quantify how vibrationally excited species, reactive radicals, and surface charging effects on enhancing the catalysis via developing DFT-based microkinetic modeling (MKM) and active learning. Deep and active learning of plasma-nanoparticle interactions effects on enhancing ammonia cracking to generate hydrogen represents a new paradigm for designing high performance plasma materials. The fundamental science of how plasma-nanoparticle interactions will change the plasma kinetics and will improve the energy efficiency for decarbonization and sustainability. The interpretable and physics-informed machine learning model will accelerate low temperature plasma chemistry and material discovery with physics rules and model interpretation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗