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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electric Field Measurements to Investigate Sheath Formation in a Nanosecond Pulsed Discharge

The objective of the research was to study the temporally and spatially resolved dynamics of sheath formation in a nanosecond pulsed plasma generated in a parallel plate electrode geometry at near atmospheric pressures. We have performed LIF dip measurements at the Sandia Low-Temperature Plasma Research Facility showing the capability to measure electric fields in near atmospheric pressure plasmas although with insufficient spatial resolution near the electrodes. Stark polarization spectroscopy implemented at the University of Minnesota was shown to be able to probe the formation of the sheaths in nanosecond pulsed plasmas. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field. High resolution Stark polarization spectroscopy showed sheaths with dimensions of the order of 50 μm underlining the need of micrometer spatial resolution for sheath characterization in near atmospheric pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Townsend to glow discharge transition for a nanosecond pulse plasma in helium: space charge formation and resulting electric field dynamics

Stark polarization spectroscopy is used to investigate the temporal evolution of the electric field distribution in the cathode region of a nanosecond pulsed discharge in helium at 120 Torr. The measurements are performed on the He I transition at 492.19 nm, during the early stages of the discharge formation. The experimental results are compared with the predictions of a 1D fluid model. Time-resolved ICCD images show that the discharge develops as a diffuse, cathode-directed ionization wave with a Townsend-like feature before transitioning into a glow-like structure. Near anode instabilities characterized by filament formation were observed near the high voltage electrode. Within 30 ns, a reduction of the sheath thickness to about 250 μm is observed, coinciding with a gradual increase of the discharge current and proportional increase in electric field at the cathode. Additionally, the cathode electric field corresponding to this sheath with a thickness of 250 μm is about 40 kV cm -1 . A subsequent steep increase of the discharge current leads to a further reduction of the sheath width. The electric field evolution as obtained by the fluid model is in excellent agreement with the measurements and shows that an enhanced ionization near the cathode is causing the space charge formation responsible for the increase in electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
DOE OSTI2025

Analytical models for the enhancement of fusion reactivity by turbulence

The reactivity of fusion plasma depends not only on its local density and temperature but also, through a recently identified kinetic effect, on the relative velocities of nearby fluid elements. Turbulence on fine spatial scales, therefore, enhances fusion reactivity. The enhancement is quantified here for general subsonic turbulent flows. Leveraging this effect in the design of inertial confinement fusion experiments could enable substantial energy savings.

Fluid dynamics↗
DOE OSTI2023

Verification of gyrokinetic particle simulations of neoclassical tearing modes in fusion plasmas

Abstract The ability to simulate neoclassical tearing modes (NTMs) in the gyrokinetic toroidal code (GTC) has been developed and verified, in which ions are treated with a gyrokinetic model and electrons are treated as a resistive fluid. The simulation results demonstrate that the neoclassical bootstrap current effect can destabilize an otherwise stable classical tearing mode. In the cylindrical geometry, GTC simulations in the magnetohydrodynamic limit show quantitative agreement with the modified Rutherford theory, both in terms of the scaling law in the small island limit and in terms of the saturation level and pressure flattening effect in the large island limit. The toroidal effects are slightly destabilizing for the NTM, while the kinetic effects of thermal ions are stabilizing for the NTM and increase its excitation threshold.

Physics↗
NASA NTRS2000

Laser Lightcraft Performance

The purpose of this study is to establish the technical ground for modeling the physics of laser powered pulse detonation phenomenon. The principle of the laser power propulsion is that when high-powered laser is focused at a small area near the surface of a thruster, the intense energy causes the electrical breakdown of the working fluid (e.g. air) and forming high speed plasma (known as the inverse Bremsstrahlung, IB, effect). The intense heat and high pressure created in the plasma consequently causes the surrounding to heat up and expand until the thrust producing shock waves are formed. This complex process of gas ionization, increase in radiation absorption and the forming of plasma and shock waves will be investigated in the development of the present numerical model. In the first phase of this study, laser light focusing, radiation absorption and shock wave propagation over the entire pulsed cycle are modeled. The model geometry and test conditions of known benchmark experiments such as those in Myrabo's experiment will be employed in the numerical model validation simulations. The calculated performance data will be compared to the test data.

Chen, Yen-Sen↗
DOE OSTI2023

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Thrust for interorbital propulsion - A question of lifetime

The most promising operational window for the use of magnetoplasmadynamic (MPD) thrusters is identified to be at megawatt power levels for orbital maneuvering. For such applications, the operation of a steady-state MPD thruster system imposes stringent requirements on the lifetime of thruster surfaces which will erode due to interactions with the plasma working fluid. Basic erosion mechanisms are presented and the problems associated with measuring the erosion rates are discussed. An experimental approach that would allow the development of a phenomenological model for erosion is proposed.

Gabriel, S. B.↗
NASA NTRS1992

Electron acceleration by Alfven waves in the magnetosphere

The self-consistent electron kinetics of Alfven waves on the electron inertial scale is studied using a 2D hybrid-kinetic description. The ions follow a fluid description for Alfven waves at frequencies below the ion cyclotron frequency. The parallel electron dynamics is treated kinetically using particle-in-cell techniques. In this model, the electron plasma mode is eliminated, and only the physics of the Alfven waves is retained. At sufficiently large amplitudes, it is found that oblique Alfven waves break due to finite electron inertia in a cold plasma. The consequence of wave breaking is the formation of an electron beam which can be unstable to the beam-plasma instability. The electrons supporting the parallel current thermalize into a non-Maxwellian distribution with an energetic tail up to several keV, assuming a reasonable magnetospheric Alfven speed. In hot plasma simulations, electron trapping is the principal mechanism of electron acceleration. It is proposed that wave breaking or electron trapping of oblique Alfven waves at 1 R(E) can result in electron acceleration and may explain some observed auroral phenomena.

Hui, C.-H.↗
NASA NTRS2012

Chemical and Hydrodynamical Models of Cometary Comae

Multi-fluid modelling of the outflowing gases which sublimate from cometary nuclei as they approach the Sun is necessary for understanding the important physical and chemical processes occurring in this complex plasma. Coma chemistry models can be employed to interpret observational data and to ultimately determine chemical composition and structure of the nuclear ices and dust. We describe a combined chemical and hydrodynamical model [1] in which differential equations for the chemical abundances and the energy balance are solved as a function of distance from the cometary nucleus. The presence of negative ions (anions) in cometary comae is known from Giotto mass spectrometry of 1P/Halley. The anions O(-), OH(-), C(-), CH(-) and CN(-) have been detected, as well as unidentified anions with masses 22-65 and 85-110 amu [2]. Organic molecular anions such as C4H(-) and C6H(-) are known to have a significant impact on the charge balance of interstellar clouds and circumstellar envelopes and have been shown to act as catalysts for the gas-phase synthesis of larger hydrocarbon molecules in the ISM, but their importance in cometary comae has not yet been fully explored. We present details of new models for the chemistry of cometary comae that include atomic and molecular anions and calculate the impact of these anions on the coma physics and chemistry af the coma.

Charnley, Steven↗
NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. 1: Theoretical foundations

The theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions (3-D) is discussed. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 solar radii where the flow is supersonic) is approximately described by the nonlinear, single fluid, polytropic (magneto-) hydrodynamic equations. Efficient numerical techniques for solving this complex system of coupled, hyperbolic partial differential equations are outlined. The formulation is inviscid and nonmagnetic, but methods allow for the potential inclusion of both features with only modest modifications. One simple, highly idealized, hydrodynamic model stream is examined to illustrate the fundamental processes involved in the 3-D dynamics of stream evolution. Spatial variations in the rotational stream interaction mechanism were found to produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction.

Pizzo, V. J.↗
DOE OSTI2020

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

TRANSP-TGLF core predictive modeling of the JET DT baseline scenario

In recent years, an intense modeling activity has been focused on preparing and analyzing the second JET Deuterium–Tritium (D–T) experimental campaign DTE2. Among the numerous scientific outcomes of this campaign was the unique opportunity to test and validate the state of the art of modeling tools with fusion-relevant DT plasmas using the full metallic ITER-like wall in different scenarios. This work reports on the core predictive modeling of plasma density, electron and ion temperatures (n e , T e , T i ) performed using TRANSP (Pankin et al 2025 Comput. Phys. Commun. 312 109611) coupled with Trapped Gyro-Landau Fluid (TGLF)-SAT2 (Kinsey et al 2008 Phys. Plasmas 15 055908), (Staebler et al 2021 Nucl. Fusion 61 116007) for the JET D–T baseline scenario (I p = 3.5 MA, q 95 = 3, β N < 2, with pellet pacing) (Garzotti et al 2025 Plasma Phys. Control. Fusion 67 075011). The sensitivity to different input parameters, as the $\vec{E}$ x $\vec{B}$ shear parameterization and the values of the kinetic quantities at the boundary of the prediction domain (ρ = 0.85) has been assessed, identifying the confidence interval of the prediction results. In particular, the dependence of the electron density profile on the particle source parameters has been studied, identifying the ionization source as the main cause for the density gradient under-prediction obtained by TRANSP-TGLF, reported in (Hyun-Tae et al 2023 Nucl. Fusion 63 112004).

JET↗
DOE OSTI2022

Digitally manufactured air plasma-on-water reactor for nitrate production

The sustainable production of food to support the increasing world population is one of humanity’s most pressing challenges. Plasma activated water, produced using renewable energy, can help fulfill plants’ needs in sustainable agriculture approaches. The design, implementation, and characterization of a digitally manufactured air plasma-on-water reactor (POWR) for the synthesis of nitrate as green nitrogen fertilizer is presented. The interaction of air plasma-generated reactive oxygen and nitrogen species with water produces nitrate (NO 3 - ) and related species, which are the main nitrogen-containing nutrients for plants. The mild conditions of the operation of the POWR opens the possibility to use plastics, particularly through digital manufacturing strategies such as 3D-printing, for its fabrication. A pin-to-plate reactor configuration powered by high-voltage alternating power is chosen due to its simplicity and efficacy. A computational thermal-fluid model is used to evaluate the design and attain expected operational characteristics. The experimental characterization of the POWR encompassed design and operation parameters, namely electrode-water spacing, air flow rate, and voltage level. A machine learning approach is implemented to extract and quantify characteristic features of the plasma–water interaction, such plasma volume and plasma–water interface area. Experimental results revealed that the nitrate production rate varies linearly with dimensionless plasma volume. The design, fabrication, and characterization methods presented can be adapted to other POWRs and help enable on-demand nitrogen fertilizer production at low environmental and economic cost.

Physics↗
DOE OSTI2025

Critical fluid dynamics in two and three dimensions

We describe a numerical method for simulating stochastic fluid dynamics near a critical point in the Ising universality class. This theory is known as model H, and is expected to govern the nonequilibrium dynamics of quantum chromodynamics (QCD) near a possible critical endpoint of the phase transition between a hadron liquid and the quark-gluon plasma. The numerical algorithm is based on a Metropolis scheme, and automatically ensures that the distribution function of the hydrodynamic variables in equilibrium is independent of the transport coefficients and only governed by the microscopic free energy. We verify dynamic scaling near the critical point of a two and three-dimensional fluid and extract the associated critical exponent z. Here, we find z≃3 in three dimensions, and z≃2 for a two-dimensional fluid. In a finite system, we observe a crossover between the mean field value z=4 and the true critical exponent z≃3 (z≃2 in d=2). This crossover is governed by the values of the correlation length and the renormalized shear viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modeling of thermal and kinetic processes in non-equilibrium plasma ignition applied to a lean combustion engine

In recent years novel ignition systems have been developed to enable stable and efficient engine operations with lean mixtures. Among them, radio-frequency corona ignition systems create discharges that involve a much wider region compared to traditional spark, and produce non equilibrium plasma with high levels of active radicals and excited species. These devices considerably increase the early flame growth speed and extend stable operating limits. With the aim of expanding the knowledge on high efficiency lean-burn SI engines, this paper investigates and compares the combustion development generated by spark and corona ignitions through computational fluid dynamics, within the Reynolds-Averaged Navier-Stokes framework for turbulence modeling. In order to simultaneously take thermal and chemical effects into account, the Perfectly Stirred Reactor combustion model is used. Experimental data are also collected for validation in an optical access engine, for different mixture levels, from stoichiometric to very lean. Furthermore, the faster burn rate generated by the corona system in the initial stage of the combustion is well predicted by the simulations, in all the relative air-fuel ratio conditions. Remarkably, as the mixture becomes lean, simulations are able to capture the non-linear transition from fast to slow kernel growth, before a self-sustainable flame propagation is established. This correlates very well with the measured engine cyclic variability and the corresponding steep change in the duration of the flame kernel formation. Ultimately, this study highlights the important role of the atomic oxygen, as active radical, in promoting and enhancing the combustion initiated by a corona discharge, in addition to the volumetric ignition effect. By contrast, the validated simulations allow to explain that the high-temperature thermal plasma generated in a traditional spark discharge is insensitive to kinetic aspects.

42 ENGINEERING↗
NASA NTRS2005

Propagation of Interplanetary Disturbances in the Outer Heliosphere

Contents include the following: 1. We have developed a one-dimensional, spherically symmetric, multi-fluid MHD model that includes solar wind protons and electrons, pickup ions, and interstellar neutral hydrogen. This model advances the existing solar wind models for the outer heliosphere in two important ways: one is that it distinguishes solar wind protons from pickup ions, and the other is that it allows for energy transfer from pickup ions to the solar wind protons. Model results compare favorably with the Voyager 2 observations. 2. 2. Solar wind slowdown and interstellar neutral density. The solar wind in the outer heliosphere is fundamentally different from that in the inner heliosphere since the effects of interstellar neutrals become significant. 3. ICME propagation from the inner to outer heliosphere. Large coronal mass ejections (CMEs) have major effects on the structure of the solar wind and the heliosphere. The plasma and magnetic field can be compressed ahead of interplanetary CMEs. 4. During the current solar cycle (Cycle 23), several major CMEs associated with solar flares produced large transient shocks which were observed by widely-separated spacecraft such as Wind at Earth and Voyager 2 beyond 60 AU. Using data from these spacecraft, we use the multi-fluid model to investigate shock propagation and interaction in the heliosphere. Specifically, we studied the Bastille Day 2000, April 2001 and Halloween 2003 events. 5. Statistical properties of the solar wind in the outer heliosphere. In a collaboration with L.F. Burlaga of GSFC, it is shown that the basic statistical properties of the solar wind in the outer heliosphere can be well produced by our model. We studied the large-scale heliospheric magnetic field strength fluctuations as a function of distance from the Sun during the declining phase of a solar cycle, using our numerical model with observations made at 1 AU during 1995 as input. 6. Radial heliospheric magnetic field events. The heliospheric magnetic field (HMF) direction, on average, conforms well to the Parker spiral.

Wang, Chi↗