Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Modeling of injected current stream-induced 3D perturbations in local helicity injection plasmas

Solenoid-free tokamak startup techniques are essential for spherical tokamaks and offer a pathway to cost reduction and design simplification in fusion energy systems. Local helicity injection (LHI) is one such approach, employing compact edge current sources to drive open field line current that initiates and sustains tokamak plasmas. The recently commissioned Pegasus-III Experiment provides a platform for advancing this and other solenoid-free startup methods. This study investigates the effect of LHI on magnetic topology in Pegasus-III plasmas. A helical filament model represents the injected current, and the linear plasma response to its three-dimensional field is calculated with M3D-C1. Poincaré mapping reveals substantial flux surface degradation in all modeled cases. The onset of overlapping magnetic structures and large-scale surface deformation begins near Ψ N ≈ 0.37, indicating a broad region of perturbed topology extending toward the edge. In rotating plasmas, both single-fluid and two-fluid models exhibit partial screening of the n = 1 perturbation, with two-fluid calculations showing stronger suppression near the edge. In contrast, the absence of rotation leads to strong resonant field amplification in the single-fluid case, while the two-fluid case with zero electron rotation mitigates this amplification and preserves edge screening. Magnetic probe measurements indicate that modeling the stream with spatial spreading—representing distributed current and/or oscillatory motion—better reproduces measured magnetic power profiles than a rigid filament model. The results underscore the role of rotation and two-fluid physics in screening stream perturbations and point to plasma flow measurements and refined stream models as key steps toward improving predictive fidelity.

Schaefer, Carolyn E. [Univ. of Wisconsin, Madison,↗
DOE OSTI2026

Effects of plasma squareness on RMP response in MAST-U

The effect of the squareness of the plasma boundary shape on the plasma response to the applied n = 2 (n is the toroidal mode number) resonant magnetic perturbation, and consequently on controlling the edge-localized mode (ELM) in MAST-U in the double-null configuration, is numerically systematically investigated. The modeling is performed utilizing three different plasma response models: the single-fluid model, the two-fluid model and the magnetohydrodynamic-kinetic hybrid model. The effect of the plasma squareness on ELM control is found to depend on the equilbrium kinetic profiles. Fixing the equilibrium radial profiles to that from a recent MAST-U discharge, the experimental plasma shape is found to be close to optimal for ELM control in terms of squareness. The optimal coil phasing does not significantly vary with varying plasma squareness. These findings robustly hold with (i) two figures of merit (FoMs) defined for predicting the ELM control feasibility, and (ii) three plasma response models employed for evaluating these FoMs.

MAST-U↗
NASA NTRS1996

Kinetic Space Weather: Toward a Global Hybrid Model of the Polar Ionosphere-Lower Magnetosphere Plasma Transport

During the indicated period of performance, we had a number of publications concerned with kinetic polar ionosphere-lower magnetosphere plasma transport. For the IUGG 1991-4 Quadrennial Report, we reviewed aspects of U.S. accomplishments concerned with polar plasma transport, among other issues. In another review, we examined the computer simulations of multiple-scale processes in space plasmas, including polar plasma outflow and transport. We also examined specifically multiscale processes in ionospheric outflows. We developed a Generalized Semi-Kinetic(GSK) model for the topside-lower magnetosphere which explored the synergistic action of wave heating and electric potentials in the formation of auroral Ion conics, in particular the "pressure cooker" mechanism. We extended the GSK model all the way down to 120 km and applied this code to illustrate the response of the ionosphere- magnetosphere to soft-electron precipitation and convection-driven frictional ion heating, respectively. Later, the convection-driven heating work was extended to a paper for the Journal of Geophysical Research. In addition to the above full published papers, we also presented the first developments of the coupled fluid-semikinetic model for polar plasma transport during this period. The results from a steady-state treatment were presented, with the second presentation being concerned with the effects of photo-electrons on the polar wind, and the first garnering an outstanding student paper award from the American Geophysical Union. We presented the first results from a time-dependent version of this coupled fluid-semikinetic model.

Horwitz, James L.↗
NASA NTRS1986

Density fluctuations in the interstellar medium: Evidence for anisotropic magnetogasdynamic turbulence. II - Stationary structures

A model of anisotropic, plasma-fluid variations was used to investigate the unknown origin of the power spectra of interstellar electron fluctuations inferred by Armstrong, Cordes, and Rickett (1981). The modeled electron variations are interpreted as density components of an anisotropic stationary mode of nonlinear magnetogasdynamics-tangential pressure balances. It is suggested that the wavenumber spectra of electron variations are identical to the spectra of the convecting velocity fields over a wide range of wavenumbers.

Higdon, J. C.↗
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 OSTI2023

Coherent fluctuations in a partially magnetized small magnetron discharge operating on molecular nitrogen

We report experimental data on rotating plasma disturbances in a quasi-steady direct-current magnetron discharge operating with molecular nitrogen at pressures of 27, 40, 67, and 133 Pa. Experiments with a segmented anode measuring the local current flow reveal up to 12 coherent structures (“spokes”) propagating in the +E×B and −E×B azimuthal directions. The current–voltage characteristics are studied and found to follow a negative resistance profile at the lower discharge currents. The existence of coherent structures and the number of spokes observed depend on the discharge current and operating pressure. Lower pressures are more favorable to the retrograde −E×B spoke rotation, which is observed in the 27, 40, and 67 Pa discharges. The results are interpreted with a gradient-driven drift wave model, and using the dispersion relation, we show that the propagation direction depends on the orientation of the local electric field within the plasma. Using an electron fluid model, we show how experimental conditions can lead to field reversals that generate −E×B spoke rotation in the lower pressure discharges.

Physics↗
NASA NTRS1986

Theory of plasma contractors for electrodynamic tethered satellite systems

Recent data from ground and space experiments indicate that plasma releases from an object dramatically reduce the sheath impedance between the object and the ambient plasma surrounding it. Available data is in qualitative accord with the theory developed to quantify the flow of current in the sheath. Electron transport in the theory is based on a fluid model of a collisionless plasma with an effective collision frequency comparable to frequencies of plasma oscillations. The theory leads to low effective impedances varying inversely with the square root of the injected plasma density. To support such a low impedance mode of operation using an argon plasma source for example requires that only one argon ion be injected for each thirty electrons extracted from the ambient plasma. The required plasma flow rates are quite low; to extract one ampere of electron current requires a mass flow rate of about one gram of argon per day.

Parks, D. E.↗
NASA NTRS1987

Theory of plasma contactors for electrodynamic tethered satellite systems

Recent data from ground and space experiments indicate that plasma releases from an object dramatically reduce the sheath impedance between the object and the ambient plasma surrounding it. Available data is in qualitative accord with the theory developed to quantify the flow of current in the sheath. Electron transport in the theory is based on a fluid model of a collisionless plasma with an effective collision frequency comparable to frequencies of plasma oscillations. The theory leads to low effective impedances varying inversely with the square root of the injected plasma density. To support such a low impedance mode of operation using an argon plasma source, for example, requires that only one argon ion be injected for each thirty electrons extracted from the ambient plasma. The required plasma flow rates are quite low; to extract one ampere of electron current requires a mass flow rate of about one gram of argon per day.

Parks, D. E.↗
NASA NTRS1989

Phase conjugation by four-wave mixing in inhomogeneous plasmas

The effects of density, temperature, and velocity gradients on four-wave mixing (FWM) in a plasma are investigated. A fluid model is used in which the stimulated Brillouin terms are included, but pump depletion is neglected. The steady state phase conjugate reflectivity and signal transmission coefficients are calculated and discussed for both degenerate and resonant FWM. The substantial effects of inhomogeneity on the use of FWM as a plasma diagnostic are discussed.

Williams, Edward A.↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Drift-kinetic effects of tungsten on plasma response to RMP in ITER

Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.

ITER↗
DOE OSTI2022

How sheath properties change with gas pressure: modeling and simulation

Abstract Particle-in-cell simulations are used to study how neutral pressure influences plasma properties at the sheath edge. The high rate of ion–neutral collisions at pressures above several mTorr are found to cause a decrease in the ion velocity at the sheath edge (collisional Bohm criterion), a decrease in the edge-to-center density ratio ( h l factor), and an increase in the sheath width and sheath potential drop. A comparison with existing analytic models generally indicates favorable agreement, but with some distinctions. One is that models for the h l factor need to be made consistent with the collisional Bohm criterion. With this and similar corrections, a comprehensive fluid-based model of the plasma boundary transition is constructed that compares well with the simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

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↗
DOE OSTI2021

Propagation of positive discharges in an air bubble having an embedded water droplet

Low temperature plasmas generated inside gas bubbles immersed in water is an effective method of rapidly transferring plasma generated reactive species to the water for applications in biomedicine, agriculture and environment. Reactive species are generally produced in the gas phase plasma and then solvate into the liquid. The large surface-to-volume ratio (SVR) of the bubble accelerates this process. In generating bubbles in water, aerosols and droplets are also contained within the bubble. These droplets also have a large SVR and so can be rapidly plasma activated. However, the presence of the droplets can also impact the propagation of the plasma in the bubble. In this paper, results are discussed from computational and experimental investigations of the formation and evolution of discharges in an air bubble immersed in water with an embedded water droplet. The computations were performed with a two-dimensional plasma hydrodynamics model. Experiments were performed with a quasi-2D bubble apparatus. In bubbles having a droplet, a plasma filament typically bridges from the powered electrode to the droplet, and then from the droplet to the bubble surface. A surface-hugging streamer also occurs on the inner bubble surface and on the surface of the droplet. Both surface streamers result in part from surface charge accumulation and can dominate the formation of reactive species that transport into the droplet. Increasing droplet conductivity suppresses propagation of the surface discharge and leads to a lower density of aqueous reactive species. Increasing conductivity of the surrounding water does not change the overall structure of the discharge but does slightly elevate the discharge intensity. Finally, the size and shape of the embedded droplet can significantly affect the formation and propagation of the streamer.

2D fluid modelling↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Validation of SOLPS-ITER simulations with kinetic, fluid, and hybrid neutral models for JET-ILW low-confinement mode plasmas

For JET L-mode plasmas in low-recycling conditions (electron temperature at the outer strike point, T e,ot ≳ 30 eV ), SOLPS-ITER simulations agree within the error bars for the experimental profiles at the low-field side (LFS) divertor target. The peak Balmer-α (Dα ) emission in the LFS divertor agrees within the error bars of the KS3 filterscope diagnostic, but is approximately 30% lower than the peak value of the KT1 spectrometer. Simulations have been performed with fluid, kinetic, and hybrid models for the neutrals. The large fluid-kinetic discrepancies of more than a factor 2 are successfully corrected by using a hybrid fluid-kinetic approach, for which kinetic atoms are transferred to the fluid population when the local Knudsen number of the atom becomes smaller than a user-defined transition Knudsen number Kn t . The hybrid-kinetic discrepancies are limited to a few % for Kn t ≤ 100 . When increasing the upstream density to high-recycling conditions, at the onset of detachment ( T e,ot ≈ 5 eV ), the simulations predict more than a factor 2 lower peak ion saturation current to the LFS divertor than the experiments. Also the Dα emission is underpredicted with approximately a factor 2. For these high-recycling conditions, the fluid-kinetic discrepancies are limited to maximum 50%, which are again corrected by using the hybrid approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1973

Turbulent heating of colliding streams in the solar wind.

Turbulent heating of colliding plasma streams has previously been observed in the solar wind. The original data were interpreted in terms of a fluid model. It is contended that a plasma-kinetic description is the more appropriate theoretical approach and is necessary in order to better understand the microscopic physical phenomena that underlie all fluid models. Microscopic solar-wind parameters characteristic of conditions during the observations were used, together with the quasi-linear plasma-kinetic theory, to compute the expected magnetic field and temperature enhancements in the interaction region between two counterstreaming plasma beams. The physical mechanism of excitation is the electromagnetic two-stream instability in which Alfven waves are unstable. A total field in the interaction region of about 8 gamma and a change in temperature of about 100,000 K are obtained.

Goldstein, M. L.↗