Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma modeling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Plasma contactor modeling with NASCAP/LEO - Extending laboratory results to space systems

In the laboratory, hollow cathode-based plasma contactors have been observed to both emit and collect ampere-level electron currents with low impedance. The laboratory behavior of hollow cathode-based plasma contactors and the limited space experience with hollow cathodes suggest that, for many applications, a hollow cathode-based plasma contactor is the ideal device to provide electrical connection with the space plasma. In order to confidently extend the laboratory experience to the low-earth-orbit environment, a series of plasma contactor computer models has been developed. Calculations show that a hollow cathode plasma contactor that collects 0.5 A in the laboratory will only collect 2.4 mA in space. The simplest way to boost the collected current is to increase the gas flow. A mole of gas is enough to collect ampere level currents for 5-1/2 hours.

Davis, V. A.↗
NASA NTRS1988

Constraints on the neutral hydrogen corona at Uranus from its interaction with magnetospheric plasma

Models of collisional reactions between the hot (several keV) protons and neutral atomic hydrogen corona observed at Uranus by Voyager 2 are compared with observations in order to deduce constraints on the coronal density. An upper limit is found for n sub H (5R sub U) of several tens of H/cu cm on the sunward side of Uranus and several H/cu cm on the dark side. The possibility that the hot plasma population might arise from photoionization and electron impact ionization of coronal and interstellar medium hydrogen in the magnetotail, followed by pickup and energization while convecting toward Uranus, is examined and rejected.

Herbert, Floyd↗
NASA NTRS2000

IMAGE EUV Observations and Modeling of the Plasmaspheric Density Trough Associated with the 24 May 2000 Geomagnetic Storm

The IMAGE EUV imager observed a plasmaspheric density trough in association with a geomagnetically active period on 24 May 2000. At EUV wavelengths, this density trough appeared as an Archimedes spiral extending from Earth's shadow to approximately 1800 MLT. We present an analysis of this density trough using simulated EUV images. Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) [Gallagher et al., 20001 and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.

Adrian, M.L.↗
NASA NTRS1990

A two-dimensional model of plasma expansion in the ionosphere

A systematic parameter study is conducted of the motion of artificial plasma clouds across the geomagnetic field at ionospheric altitudes. The study is based on a two-dimensional numerical model. Taken into consideration are the effects solar cycle variations on the cloud evolution, as well as the effects of seasonal and geomagnetic activity conditions, and variations in the ionospheric release conditions. The effects of magnetospheric electric fields and thermospheric winds are also considered. The plasma cloud model and related assumptions are discussed, and numerical results are presented.

Ma, T.-Z.↗
NASA NTRS1988

The kinetic approach in magnetospheric plasma transport modeling

The need for a kinetic approach in magnetospheric plasma transport problems is reviewed, as are the trends in its recent applications. The need for kinetic modeling is particularly obvious when confronted with the astonishing variety of magnetospheric particle measurements that display compelling energy and pitch angle-related spatial and/or temporal dispersion, and various types of highly non-Maxwellian features in the distribution functions. Global problems in which the kinetic approach has recently been applied include solar wind plasma injection and dispersion over the cusp, substorm particle injection near synchronous orbit, synergistic energization of ionospheric ions into ring current populations by waves and induced electric field-driven convection, and ionospheric outflow from restricted source regions into the magnetosphere. Kinetic modeling can include efforts ranging from test-particle techniques to particle-in-cell studies, and this range is considered here. There are some areas where fluid and kinetic approaches have been combined or patched together, and these will be briefly discussed.

Horwitz, J. L.↗
NASA NTRS2010

Recent Simulation Results on Ring Current Dynamics Using the Comprehensive Ring Current Model

Plasma sheet conditions and electromagnetic field configurations are both crucial in determining ring current evolution and connection to the ionosphere. In this presentation, we investigate how different conditions of plasma sheet distribution affect ring current properties. Results include comparative studies in 1) varying the radial distance of the plasma sheet boundary; 2) varying local time distribution of the source population; 3) varying the source spectra. Our results show that a source located farther away leads to a stronger ring current than a source that is closer to the Earth. Local time distribution of the source plays an important role in determining both the radial and azimuthal (local time) location of the ring current peak pressure. We found that post-midnight source locations generally lead to a stronger ring current. This finding is in agreement with Lavraud et al.. However, our results do not exhibit any simple dependence of the local time distribution of the peak ring current (within the lower energy range) on the local time distribution of the source, as suggested by Lavraud et al. [2008]. In addition, we will show how different specifications of the magnetic field in the simulation domain affect ring current dynamics in reference to the 20 November 2007 storm, which include initial results on coupling the CRCM with a three-dimensional (3-D) plasma force balance code to achieve self-consistency in the magnetic field.

Zheng, Yihua↗
DOE OSTI2026

Divertor detachment characterization in negative triangularity discharges in DIII-D via 2D edge-plasma transport modeling

Edge fluid modeling of the first divertor-plasma detachment experiments in negative triangularity (NT) discharges in the DIII-D tokamak is presented using the 2D multifluid edge transport code UEDGE, including cross-field particle drifts. Experiments were performed where the lower single-null magnetic equilibrium had a strong NT (δ≈−0.5), that is, where the magnetic X-point is at a larger major radius than the core magnetic axis. Divertor-plasma detachment was induced by increasing the core plasma density in DIII-D via intrinsic gas puffing. Here density scans are performed with UEDGE to reach a detached plasma and to quantitatively recover the experimental rollover of the ion saturation current on the outer divertor target plate. The simulations cover experiments with both signs of the toroidal magnetic field, B T , where the ion magnetic Grad-B drifts are directed into (forward B T ) and out of (reverse B T ) the divertor region. Consistent with experiments with neutral beam power injection, the NT simulations reproduce: 40% higher density is needed to reach detachment onset with forward B T compared with reverse B T , and the absence of deep detachment is found with reverse B T . Similarly, comparison between Ohmic discharges in NT and positive triangularity (PT) shaping confirms that a substantially higher density is needed to achieve detachment in NT than in PT, with NT requiring an line-average density of at least the Greenwald density or higher. Simulation results suggest that higher densities are needed to reach detachment in negative compared to PT because these discharges have a shorter midplane-to-target distance along the total magnetic field B, a shorter outer divertor poloidal leg length (0.06 m vs 0.2 m), and reduced radial transport [near-scrape-off layer (SOL) D ⊥ /χ ⊥ =0.3/0.5 vs D ⊥ /χ ⊥ =1.0/1.0, all in (m 2 /s)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On Ohm’s law in reduced plasma fluid models

Drift-reduced MHD models are widely used to study magnetised plasma phenomena, in particular for magnetically confined fusion applications, as well as in solar and astrophysical research. This letter discusses the choice of Ohm's law in these models, the resulting dispersion relations for the dynamics parallel to the magnetic field, and the implications for numerical simulations. We find that if electron pressure is included in Ohm's law, then both electromagnetic and finite electron mass effects must also be included in order to obtain physical dispersion relations. A simple modification to the plasma vorticity is also found which improves handling of low density regions, of particular relevance to the simulation of the boundary region of magnetised plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma diagnostics and modeling of lithium-containing plasmas

Thin-film deposition from chemically reactive multi-component plasmas is complex, and the lack of electron collision cross-sections for even the most common metalorganic precursors and their fragments complicates their modeling based on fundamental plasma physics. This study focuses on understanding the plasma physics and chemistry in argon (Ar) plasmas containing lithium bis (trimethylsilyl) amide used to deposit Li x Si y thin films. These films are emerging as potential solid electrolytes for lithium-ion batteries, and the Li-to-Si ratio is a crucial parameter to enhance their ionic conductivity. We deposited Li x Si y films in an axial flow-through plasma reactor and studied the factors that determine the variation of the Li-to-Si ratio in films deposited at various points on a substrate spanning the entire reactor axis. While the Li-to-Si ratio is 1:2 in the precursor, the Li-to-Si ratio is as high as 3:1 in films deposited near the plasma entrance and decreases to 1:1 for films deposited downstream. Optical emission from the plasma is dominated by Li emission near the entrance, but Li emission disappears downstream, which we attribute to the complete consumption of the precursor. We hypothesized that the axially decreasing precursor concentration affects the electron energy distribution function in a way that causes different dissociation efficiencies for the production of Li and Si. We used Li line intensities to estimate the local precursor concentration and Ar line ratios to estimate the local reduced electric field to test this hypothesis. This analysis suggests that the mean electron energy increases along the reactor axis with decreasing precursor concentration. As a result, the decreasing Li-to-Si ratio with axially decreasing precursor concentration may be explained by Li release from the precursor having lower threshold energy than Si release.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1981

Fluid model of plasma outside a hollow cathode neutralizer

The present study analyzes the capability of a fluid model of electron transport to explain observed properties of the external plasma of a hollow cathode neutralizer used to neutralize beams emerging from ion thrusters. Calculations reported here show that when the effective collision frequency in such a model is near the plasma frequency, the resulting electric potential and electron temperature variations are in qualitative agreement with values measured in the plume mode of the hollow cathode. Both theory and experiment show strong variations of temperature and potential within a few centimeters of the cathode orifice.

Parks, D. E.↗
DOE OSTI2024

1D modeling of plasma streamers at ammonia-air flame conditions

Abstract Self-consistent 1D modeling of streamers in ammonia-oxygen-nitrogen-water mixtures has been performed in this work. A fluid model that includes species transport, electrostatic potential, and detailed chemistry was developed and verified. This model is then used to simulate the avalanche, streamer formation and propagation phases, driven by a nanosecond voltage pulse, at different thermochemical conditions derived from a 1D laminar premixed ammonia-air flame. The applicability of the Meek’s criterion in predicting the streamer inception location was successfully confirmed. Streamer formation and propagation duration were found to vary significantly with different thermochemical conditions, due to the difference in ionization rates. The thermochemical state also affected the breakdown characteristics which was tested by maintaining the background reduced electric field constant. Detailed kinetic analyses revealed the importance of O ( 1 D ) in the production of key radicals, such as O, OH, and NH 2 . Furthermore, the contributions of the dissociative electronic excitation of NH 3 towards the production of H and NH 2 radicals have also been reported. Spatial and temporal evolution of the electron energy loss fractions for various inelastic collision processes at different thermochemical states uncovered the input plasma energy spent of fuel dissociation and the large variability in the dominant processes during the avalanche and streamer propagation phases. The methodology and analyses reported in this work are key towards developing effective strategies for controlled nanosecond-pulsed non-equilibrium plasma sources used for ammonia ignition and flame stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Full wave modeling of radio-frequency beams in tokamaks in the electron cyclotron frequency range

Simulation of full wave, without paraxial approximation, high-resolution solution of wave equations in frequency domain in the electron cyclotron resonance (ECR) frequency range for realistic Tokamak plasma parameters became possible by using recently formulated hybrid iterative algorithm [Svidzinski et al., Phys. Plasmas 25, 082509 (2018)] for numerically solving discretized wave equations. This approach combines time evolution and iterative relaxation techniques into iteration cycles. This algorithm is implemented in 2D code FullWave, solving wave equations in Tokamaks in cold and hot plasma models, and it has been tested in 3D full wave iterative RF beams simulation tool, which is presently being developed to model 3D ECRH RF beams in fusion devices using dynamic grid adaptation. The results of 2D full wave modeling, assuming specified toroidal mode number, of ECRH RF beams in DIII-D plasma, performed in the cold and hot plasma models for outboard and top launch scenarios using FullWave are presented. Nonlocal hot plasma response model, based on accurate numerical solution of linearized Vlasov equation, is used to model beam propagation and absorption in the 2nd electron cyclotron harmonic region. Demonstration of capability of the hybrid iterative algorithm to model ECRH RF beams in 3D is made by simulating a substantial part of realistic beam in DIII-D, launched from outboard side of the machine. All relevant physics of RF beam propagation, most of which is not captured in paraxial approximation, such as beam's divergence, interference between the X and O modes in the beam, X-O mode conversion, beam splitting into the X and O mode beams, transformation of beam's cross section, and absorption at the 2nd electron cyclotron harmonic, is captured in the simulations. A numerical technique to find an optimal beam polarization at the launcher to launch a nearly pure X or O mode beam in plasma is developed and tested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Screening of resonant magnetic perturbation fields assuming various plasma flow models

A recently updated version of the MARS-F code [Y. Q. Liu et al., Phys. Plasmas 7, 3681 (2000); L. Li et al., Phys. Plasmas 25, 082512 (2018); and G. L. Xia et al., Nucl. Fusion 59, 126035 (2019)] is utilized to numerically investigate the plasma screening effect on the applied resonant magnetic perturbation (RMP) field, assuming various equilibrium flow models, including the toroidal flow, the parallel flow and their combinations, and poloidal and toroidal projections of the parallel flow. A parallel equilibrium flow with a uniform radial profile is found to have no effect on plasma screening of the RMP field. A sheared parallel flow, however, does change plasma screening. The poloidal projection of the parallel flow weakens plasma screening in the resistive-inertial regime. The effect on the favorable average curvature regime is found, however, to be non-monotonic. With the increasing flow speed, the poloidal projection first weakens Glasser-Green-Johnson (GGJ)-screening. Further increase in the flow speed results in enhanced GGJ-screening again. This non-monotonic behavior is related to the perturbed parallel shielding current, which appears also off the mode rational surface at fast flow due to additional resonances between the RMP perturbation and the sound wave continuum. These results indicate that flow induced plasma screening to the RMP field can have complicated characteristics, which, in turn, can have implications on the RMP field penetration into the plasma in experiments for controlling the edge localized modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

The Gamma-Ray Pulsar Phenomenology in View of 3D Kinetic Global Magnetosphere Models

We develop kinetic plasma models of pulsar magnetospheres with magnetic-field-line-dependent plasma injection that reveal the importance of various magnetosphere regions in regulating the gamma-ray emission. We set different particle injection rates for the so-called open, closed, and separatrix zones. Moderate particle injection rates in open and closed zones ensure a global field structure close to the force-free one, while the dissipation occurs mainly in and around the equatorial current sheet. The particles injected in the separatrix zone affect the particle populations that enter the equatorial current sheet region and, therefore, the corresponding accelerating electric fields, particle energies, the spectral cutoff energy, and gamma-ray efficiency. The separatrix zone models reproduce the recently discovered fundamental plane of gamma-ray pulsars consistent with curvature radiation emission, the gamma-ray light-curve shapes, and the radio-lag vs. peak-separation correlation reported in the Fermi second pulsar catalog. The model beaming factors indicate that the pulsar total gamma-ray luminosities listed in the Fermi catalogs are overestimations of the actual ones. We find that the radiation reaction limited regime starts ceasing to govern the high-energy emission for spin-down powers less than 1034 erg/s. Our results also indicate that toward high magnetic inclination angles, the "Y point" around the rotational equator migrates well inside the light cylinder sparking additional peaks in the gamma-ray pulse profiles. We find that an equivalent enhanced particle injection beyond the Y point strengthens these features making the model gamma-ray light curves inconsistent with those observed.

Neutron stars↗
DOE OSTI2021

Data Needs for Modeling Low-Temperature Non-Equilibrium Plasmas: The LXCat Project, History, Perspectives and a Tutorial

Technologies based on non-equilibrium, low-temperature plasmas are ubiquitous in today’s society. Plasma modeling plays an essential role in their understanding, development and optimization. An accurate description of electron and ion collisions with neutrals and their transport is required to correctly describe plasma properties as a function of external parameters. LXCat is an open-access, web-based platform for storing, exchanging and manipulating data needed for modeling the electron and ion components of non-equilibrium, low-temperature plasmas. The data types supported by LXCat are electron- and ion-scattering cross-sections with neutrals (total and differential), interaction potentials, oscillator strengths, and electron- and ion-swarm/transport parameters. Online tools allow users to identify and compare the data through plotting routines, and use the data to generate swarm parameters and reaction rates with the integrated electron Boltzmann solver. In this review, the historical evolution of the project and some perspectives on its future are discussed together with a tutorial review for using data from LXCat.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗