Search NASA⌕ Search

SEARCH · Search NASA

Results for “non-equilibrium plasma”

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 55 records · Page 3

Termination of the solar wind in the hot, partially ionized interstellar medium

Theoretical foundations for understanding the problem of the termination of the solar wind are reexamined in the light of most recent findings concerning the states of the solar wind and the local interstellar medium. The investigation suggests that a simple extention of Parker's (1961) analytical model provides a useful approximate description of the combined solar wind, interstellar wind plasma flowfield under conditions presently thought to occur. A linear perturbation solution exhibiting both the effects of photoionization and charge exchange is obtained for the supersonic solar wind. A numerical algorithm is described for computing moments of the non-equilibrium hydrogen distribution function and associated source terms for the MHD equations. Computed using the algorithm in conjunction with the extended Parker solution to approximate the plasma flowfield, profiles of hydrogen number density are given in the solar wind along the upstream and downstream axes of flow with respect to the direction of the interstellar wind. Predictions of solar Lyman-alpha backscatter intensities to be observed at 1 a.u. have been computed, in turn, from a set of such hydrogen number density profiles varied over assumed conditions of the interstellar wind.

Lombard, C. K.↗

Computer Simulation of Motion of Heliospheric Termination Shock

A time-dependent computer model is used to study a non-equilibrium structure of the heliosphere, resulting from the interaction between the solar wind and the interstellar plasma. An interaction of the heliospheric termination shock with various disturbances such as density or velocity jumps or interplanetary shocks in upstream solar wind as well as pressure jumps in local interstellar medium is investigated. Consequences of latitudinal variation in dynamic pressure are also discussed. The initial boundary conditions of unsteady calculations are given by the solar wind parameters as functions of time and/or latitude on an inner boundary, and/or by the local interstellar pressure as a function of time (roughly simulating the effect of the local interstellar medium) on an outer boundary.

Ratkiewicz, R. E.↗

Non-equilibrium ionization in Puppis A

High resolution X-ray spectroscopy of the brightest knot of emission in the Puppis A supernova remnant shows that it is made up of ionizing plasma, far from equilibrium. Flux measurements in several X-ray lines make it possible to determine the nonequilibrium conditions: electron temperature, ion populations, and time since the knot was heated by the supernova shock. Imaging and spectroscopic data from the Einstein Observatory together suggest that this knot is a cloud of density about 10 per cu cm which has recently been shocked to a temperature 7,000,000 K. Radio and optical data on the region appear consistent with this picture.

Winkler, P. F.↗

Coronal and Prominence Plasmas

Various aspects of solar prominences and the solar corona are discussed. The formation of prominences, prominence diagnostics and structure, prominence dissappearance, large scale coronal structure, coronal diagnostics, small scale coronal structure, and non-equilibrium/coronal heating are among the topics covered.

Poland, Arthur I.↗

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A study of physical mechanisms for filament eruption and coronal mass ejection via numerical simulation

It is well-known from both observation and theory that photospheric shear motion has played a key role in physical processes of the energy build-up and release for the solar flare. In order to further our understanding of the occurrence of solar flares it is necessary to investigate the triggering mechanism. One popular scenario for the onset of eruptive solar flares is that in response to photospheric shear motions the magnetic field evolves slowly through a series of magnetohydrodynamic-equilibria until a threshold is reached where magnetohydrodynamic (MHD) non-equilibrium sets in. Thus, a magnetic eruption occurs, causing the solar flare. To substantiate this claim we have employed our newly developed three-dimensional time-dependent MHD code with gravity to simulate the evolution of the coronal field. We use plasma beta = 0.1 to closely approximate the condition in the actual corona. Some preliminary results are presented.

Wu, S. T.↗

Non-equilibrium of high-temperature nitrogen and oxygen

A stationary colloidal-radiative model for determination of production of electrons, excited atoms, and spectral line intensities in high-temperature atomic N and O is presented. The populations of several atomic levels and partition functions, and the intensities of several lines, are calculated in a wide range of conditions. Transport of radiation is included by coupling the rate equations for production of the electrons and excited atoms with the Holstein escape factors which are not constant, but depend on plasma conditions.

Kunc, J. A.↗

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisional-radiative non-equilibrium in high-temperature helium

A stationary, nonlinear collisional-radiative model for determination of production of electrons, excited atoms and atomic line intensities for high-temperature helium is presented. The populations of several atomic levels and the intensities of spectral and continuum lines are calculated in a wide range of conditions (Te = 12,000-22,000 K, Nt = 10 to the 12 - 10 to the 20th per cu cm. Transport of radiation is included by coupling the rate equations for production of the electrons and excited atoms with the radiation escape factors which are not constant but depend on plasma conditions.

Kunc, J. A.↗

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

stagnation heat flux↗

Turbulence statistical analysis of the L-H transition and RMPs in KSTAR

Here, we investigate the turbulence statistics associated with low-to-high confinement (L-H) transitions and externally applied resonant magnetic perturbations (RMPs) in KSTAR. Time-series fluctuations of electron density n e , electron temperature T e , and the time derivative of the poloidal magnetic field dB θ /dt (Mirnov coils) are analysed using information-geometric measures (information rate Γ and information length $\mathcal{L}$ = ∫ Γ dt), together with kurtosis κ and variance σ 2 . In low-density upper single-null plasmas (n e ~ 1.2 x 10 19 m -3 ), a ~80 kHz magnetic mode coupling n e , T e , dB θ /dt and emerges prior to the L-H transition and persists into the edge-localised modes H-mode. Edge-localised RMPs (ERMPs) suppress this coherent mode but enhance intermittency, producing frequent bursts that abruptly reshape the time-dependent probability density functions (PDFs) and generate large spikes in Γ (with smaller changes in κ), signalling ERMP-driven departures from quasi-stationarity. The impact of ERMPs on background fluctuation levels depends on density, radial location, and the fluctuating variable itself ($\tilde{n}$, $\tilde{T}$, $\dot{B}$ θ ), whereas $\mathcal{L}$ provides a robust, regime-agnostic measure of cumulative statistical reorganisation and spatial decorrelation. In particular, at low density we observe weaker coupling between $\tilde{n}$ and $\tilde{T}$, along with a tendency toward decreased radial correlation-most clearly for $\tilde{T}$-under ERMPs. Overall, information geometry cleanly captures intermittent events, quantifies non-equilibrium PDF evolution, and offers a compact, cross-diagnostic metric for assessing resonant magnetic perturbation effects on edge transport and correlation across densities, radial locations, and confinement states.

Kim, Eun-jin [Coventry Univ. (United Kingdom); Seo↗

Turbulent mixing layers in the interstellar medium of galaxies

We propose that turbulent mixing layers are common in the interstellar medium (ISM). Injection of kinetic energy into the ISM by supernovae and stellar winds, in combination with density and temperature inhomogeneities, results in shear flows. Such flows will become turbulent due to the high Reynolds number (low viscosity) of the ISM plasma. These turbulent boundary layers will be particularly interesting where the shear flow occurs at boundaries of hot (approximately 10(exp 6) K) and cold or warm (10(exp 2) - 10(exp 4) K) gas. Mixing will occur in such layers producing intermediate-temperature gas at T is approximately equal to 10(exp 5.0) - 10(exp 5.5) that radiates strongly in the optical, ultraviolet, and EUV. We have modeled these layers under the assumptions of rapid mixing down to the atomic level and steady flow. By including the effects of non-equilibrium ionization and self-photoionization of the gas as it cools after mixing, we predict the intensities of numerous optical, infrared, and ultraviolet emission lines, as well as absorption column densities of C 4, N 5, Si 4, and O 6.

Slavin, J. D.↗

Large-scale density structures in the outer heliosphere

The Plasma Science experiment on the Voyager 2 spacecraft has measured the solar wind density from 1 to 38 AU. Over this distance, the solar wind density decreases as the inverse square of the heliocentric distance. However, there are large variations in this density at a given radius. Such changes in density are the dominant cause of changes in the solar wind ram pressure in the outer heliosphere and can cause large perturbations in the location of the termination shock of the solar wind. Following a simple model suggested by Suess, we study the non-equilibrium, dynamic location of the termination shock as it responds to these pressure changes. The results of this study suggest that the termination shock is rarely if ever at its equilibrium distance and may depart from that distance by as much as 50 AU at times.

Belcher, J. W.↗

Thermochemically-Closed Sonic-Flow Inversion for Enthalpy and Temperature in Multispecies Arc-Jet Flows

A thermochemically-closed sonic-flow inversion framework (TSIF) is developed to infer bulk enthalpy and total temperature upstream of a choked nozzle in arc-jet flows. The formulation recasts a pressure-rise total enthalpy quantification technique as an inverse problem in characteristic-velocity c * space using measured mass flow rate, upstream total pressure, gas composition, and nozzle throat geometry as inputs. Unlike calorimetric energy-balance approaches or optical diagnostics, the method relies primarily on routinely measured facility quantities combined with explicit thermochemical closure. Thermochemical states are obtained using NASA’s open-source Chemical Equilibrium with Applications (CEA) code, enabling construction of a chemistry-consistent relation between characteristic velocity, total enthalpy, and total temperature under equilibrium or frozen assumptions. A discharge coefficient is self-calibrated using cold-flow (arc-off) operation data and applied to hot-flow (arc-on) measurements, enabling upstream losses to be accounted for without empirical correlations. The framework is applied to air, N 2 , and CO 2 –N 2 arc-jet flows and demonstrates expected trends for the inferred thermochemical states as function of arc power, specific energy input, mass-flow, heater configuration, and test gas. In the air limit, under equilibrium assumptions, the method recovers the classical high-enthalpy asymptotic correlation of Winovich with a mean residual of 4.4%, demonstrating compatibility with established sonic-flow scaling, while extending applicability to arbitrary multi-species mixtures and non-equilibrium chemistry. The framework provides a mixture-flexible methodology for determining bulk thermochemical states in modern arc-jet environments using routine facility pressure, mass-flow, gas-composition, and nozzle-geometry information together with a cold-flow calibration.

inviscid theory↗

In situ studies on heavy ion irradiation and partial oxidation induced stacking faults and nanograins in tungsten nanolaminates

Understanding the microstructural evolution of tungsten (W) under extreme irradiation environments is critical for its application as a plasma-facing material in future fusion reactors. Heavy ion irradiation study offers expedited damage accumulation and thus allows irradiation tolerance property prediction in a short period of time. In this study we explore in situ Kr ion irradiation of W nanolaminates at 800 °C and uncover a complex interplay of irradiation-induced transformations, including the emergence of stacking faults, nanograin formation, and generation of thickened grain boundaries defined as GB regions with the thickness of up to 15 nm. High-resolution transmission electron microscopy studies reveal the formation of extended planar faults and a metastable hexagonal close-packed (hcp) phase within the body-centered cubic (bcc) matrix. These structural transitions are facilitated by irradiation-induced shear. They are further stabilized by the presence of oxygen, which promotes stacking fault formation and vacancy trapping. Elevated temperatures enhance defect mobility, enabling dynamic recrystallization into ultrafine grains and the thickening of GBs due to defect absorption and impurity segregation. These findings unravel non-equilibrium phase transformation pathways in heavy ion irradiated W and highlight the critical role of impurity-mediated defect dynamics in governing its radiation tolerance properties.

Wazeer, Adil [Purdue University]↗

Spectral deconvolution without the deconvolution: Extracting temperature from x-ray Thomson scattering spectra without the source-and-instrument function

X-ray Thomson scattering (XRTS) probes the dynamic structure factor of the system, but the measured spectrum is broadened by the combined source-and-instrument function (SIF) of the setup. In order to extract properties such as temperature from an XRTS spectrum, the broadening by the SIF needs to be removed. Recent work [Dornheim et al. Nat. Commun. 13 , 7911 (2022)] has suggested that the SIF may be deconvolved using the two-sided Laplace transform. However, the extracted information can depend strongly on the shape of the input SIF, and the SIF is in practice challenging to measure accurately. Here, we propose an alternative approach: we demonstrate that considering ratios of Laplace-transformed XRTS spectra collected at different scattering angles is equivalent to performing the deconvolution, but without the need for explicit knowledge of the SIF. From these ratios, it is possible to directly extract the temperature from the scattering spectra, when the system is in thermal equilibrium. We find the method to be generally robust to spectral noise and physical differences between the spectrometers, and we explore situations in which the method breaks down. Furthermore, the fact that consistent temperatures can be extracted for systems in thermal equilibrium indicates that non-equilibrium effects could be identified by inconsistent temperatures of a few eV between the ratios of three or more scattering angles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electrified vapour deposition at ultrahigh temperature and atmospheric pressure for nanomaterials synthesis

Vapour-phase synthesis methods have shown promise for the scalable synthesis of nanomaterials and coatings. However, the vaporization of different precursors for the synthesis of a broad nanomaterial space, particularly at atmospheric pressure, while maintaining compositional and structural control of the final product is challenging. Here we report the generation of an ultrahigh-temperature atomic vapour at atmospheric pressure based on electrified heating, for the growth of multi-elemental nanomaterials and thin films. This process relies on a reactor design whereby solid-state precursors are vaporized within a semi-confined space beneath an electrified heater that can reach ~3,000 K. The proximity of the heater rapidly breaks down the bonds of metal salt precursors and decomposes them into an atomic vapour that expands into a high-temperature (>2,000 K), highly reactive and high-flux vapour (10 21 –10 22 atoms per cm 2 per second) that travels upwards in a directional flow. When mixed with entrained ambient gases, the highly reactive atomic species rapidly nucleate and grow into the desired final products, including alloys, oxides, sulfides and thin films, which can be deposited on a low-temperature substrate. This EVD approach can synthesize a broad range of functional nanomaterials at atmospheric pressure, including single-phase multi-elemental nanomaterials formed under thermodynamically non-equilibrium conditions.

36 MATERIALS SCIENCE↗

Classical-quantum simulation of non-equilibrium Marshak waves

In the radiation hydrodynamic simulations used to design inertial confinement fusion (ICF) and pulsed power experiments, nonlinear radiation diffusion tends to dominate CPU time. This raises the interesting question of whether a quantum algorithm can be found for nonlinear radiation diffusion which provides a quantum speedup. Recently, such a quantum algorithm was introduced based on a quantum algorithm for solving systems of nonlinear partial differential equations (PDEs) which provides a quadratic quantum speedup. Here, we apply this quantum PDE (QPDE) algorithm to the problem of a non-equilibrium Marshak wave propagating through a cold, semi-infinite, optically thick target, where the radiation and matter fields are not assumed to be in local thermodynamic equilibrium. The dynamics is governed by a coupled pair of nonlinear PDEs which are solved using the QPDE algorithm, as well as two standard PDE solvers: (i) Python's py-pde solver; and (ii) the KULL ICF simulation code developed at Lawrence-Livermore National Laboratory. We compare the simulation results obtained using the QPDE algorithm and the standard PDE solvers and find excellent agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗