Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2024

Quasilinear theory: the lost ponderomotive effects and why they matter

Quasilinear theory (QLT) has been used for modeling wave–plasma interactions for decades but remains largely heuristic. Plasma inhomogeneity, ponderomotive effects, microscopic fluctuations, and collisions are not easily accommodated from first principles in QLT, and typically are ignored entirely, due to the limitations of the standard Fourier–Laplace global-mode approach. This results in inconsistencies, for example, violation of the action conservation for nonresonant waves. However, these issues can be avoided, and the theory can be substantially generalized and corrected, if QLT is formulated using more suitable analytical tools, particularly, the Weyl symbol calculus. Here, an attempt is made to deliver an accessible review of this modern formulation, provide intuitive calculations for special cases, and elaborate on the connection with the ‘oscillation-center QLT’ originally proposed by Dewar (Phys Fluids 16:1102, 1973). A Fokker–Planck equation for a ‘dressed’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, ponderomotive forces, and interactions with background fields for a generic Hamiltonian, so many known formulations of QLT for specific plasma models become corollaries of a single unifying theory. Also, waves are allowed to be off-shell (not constrained by a dispersion relation), which allows them to accommodate microscopic fluctuations. This leads to a collision integral of the Balescu–Lenard type that has all the usual properties but is not restricted to any specific plasma model. For on-shell waves, a generalized version of the classic oscillation-center QLT is obtained. Finally, combined with the wave-kinetic equation, this formulation not only conserves particles, momentum, and energy, like the classic QLT but also reinstates the action conservation for nonresonant waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1980

Magnetoacoustic heating by ion Landau damping

The Vlasov-fluid model of Freidberg (1972) is used to study the resonance heating of a sharp-boundary screw pinch. The analysis provides the first treatment of the magnetoacoustic heating of a cylindrical plasma by means of ion Landau damping, which was identified as a viable dissipative mechanism for the conversion of magnetoacoustic wave energy into ion thermal energy. In addition, local and global energy conservation are considered, and formulae and numerical results for the thermal energy doubling time and the associated induced rf electric fields are presented. It is shown that collisionless absorption can provide a heating mechanism when an equilibrium plasma column is pumped by oscillations of the confining magnetic field at a frequency near the oblique magnetoacoustic frequency.

Turner, L.↗
NASA NTRS2004

Plasma Sail Concept Fundamentals

The mini-magnetospheric plasma propulsion (M2P2) device, originally proposed by Winglee et al., predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to unprecedented speeds of 50 C80 km/s after an acceleration period of 3 mo. Such velocities will enable travel out of the solar system in period of 7 yr almost an order of magnitude improvement over present chemical-based propulsion systems. However, for the parameters of the simulation of Winglee et al., a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the magnetohydrodynamic (MHD) fluid model, normally applied to planetary magnetospheres, that the characteristic scale size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or comparable to the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. A two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail has been adopted. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the forces delivered to the innermost regions of the plasma sail are considerably ( 10 times) smaller than the MHD counterpart, are dominated by the magnetic field pressure gradient, and are directed primarily in the transverse direction.

Khazanov, G. V.↗
DOE OSTI2021

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Striations in moderate pressure dc driven nitrogen glow discharge

Plasma stratification has been studied for more than a century. Despite the many experimental studies reported on this topic, theoretical analyses and numerical modeling of this phenomenon have been mostly limited to rare gases. In this work, a one-dimensional fluid model with detailed kinetics of electrons and vibrationally excited molecules is employed to simulate moderate-pressure (i.e. a few Torrs) dc discharge in nitrogen in a 15.5 cm long tube of radius 0.55 cm. The model also considers ambipolar diffusion to account for the radial loss of ions and electrons to the wall. The proposed model predicts self-excited standing striations in nitrogen for a range of discharge currents. The impact of electron transport parameters and reaction rates obtained from a solution of local two-term and a multi-term Boltzmann equation on the predictions are assessed. In-depth kinetic analysis indicates that the striations result from the undulations in electron temperature caused due to the interaction between ionization and vibrational reactions. Furthermore, the vibrationally excited molecules associated with the lower energy levels are found to influence nitrogen plasma stratification and the striation pattern strongly. A balance between ionization processes and electron energy transport allows the formation of the observed standing striations. Simulations were conducted for a range of discharge current densities from ~0.018 to 0.080 mA cm –2 , for an operating pressure of 0.7 Torr. Parametric studies show that the striation length decreases with increasing discharge current. Furthermore, the predictions from the model are compared against experimental measurements and are found to agree favorably.

42 ENGINEERING↗
NASA NTRS1986

Thermal catastrophe in the plasma sheet boundary layer

This letter presents a first step towards a substorm model including particle heating and transport in the plasma sheet boundary layer (PSBL). The heating mechanism discussed is resonant absorption of Alfven waves. For some assumed MHD perturbation incident from the tail lobes onto the plasma sheet, the local heating rate in the PSBL has the form of a resonance function of the one-fluid plasma temperature. Balancing the local heating by convective transport of the heated plasma toward the central plasma sheet, an 'equation of state" is found for the steady-state PSBL whose solution has the form of a mathematical catastrophe: at a critical value of a parameter containing the incident power flux, the local density, and the convection velocity, the equilibrium temperature jumps discontinuously. Associating this temperature increase with the abrupt onset of the substorm expansion phase, the catastrophe model indicates at least three ways in which the onset may be triggered. Several other consequences related to substorm dynamics are suggested by the simple catastrophe model.

Smith, Robert A.↗
NASA NTRS1978

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

The paper is concerned with the development of the theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 Rs where the flow is supersonic) is approximately described by the nonlinear, single-fluid, polytropic magnetohydrodynamic or hydrodynamic equations. Efficient numerical techniques are outlined for solving this complex system of coupled, hyperbolic partial differential equations. The present formulation is inviscid and nonmagnetic, but the methods used allow for the potential inclusion of both features with only modest modifications. A simple, highly idealized hydrodynamic model stream is examined to illustrate the fundamental processes involved in the three-dimensional dynamics of stream evolution. It is found that spatial variations in the rotational stream interaction mechanism 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. Comparison with simpler models demonstrates the essential nonlinear, multidimensional nature of the interplanetary dynamics.

Pizzo, V.↗
DOE OSTI2022

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗
NASA NTRS1982

Fluid shifts and erythropoiesis - Relevance to the 'anemia' of space flight

To model the fluid shifts thought to occur in man during space flight, cephalic fluid shifts have been induced in man subjected to horizontal or headdown bedrest, in squirrel monkeys exposed to lower body positive pressure, and in rats subjected to antiorthostatic hypokinesia. The influence on erythropoiesis of such fluid redistribution has been studied. Only in man did a cephalic fluid shift consistently and significantly lead to a plasma volume reduction and an increased hematocrit. Although there was evidence for erythrosuppression and the subjects were 'anemic' at the end of the study, serum erythropoietin titers remained normal throughout bedrest. The erythrosuppression probably did not arise due to the increased hematocrit but may have been related to P50 shifts or the loss of body weight. Each model appeared to reproduce different parts of man's physiological response to weightlessness and promises to be useful in unraveling the etiology of the 'anemia' of space flight.

Dunn, C. D. R.↗
DOE OSTI2023

Non-ideal MHD growth of current interchange tearing modes at plasma edge and response to externally-imposed flow

Current interchange tearing modes (CITMs) have been proposed to explain an intermittent eruption of electric current toward a divertor at an edge region of a tokamak device. Extended magnetohydrodynamic (MHD) simulations combined with the Scrape-Off-Layer diffusion model, which has been developed by the authors [Miura et al., Phys. Plasmas 24, 092111 (2017)], are carried out in order to study the growth of a CITM under influences of two-fluid and gyro-viscous effects, or non-ideal MHD effects. Numerical simulations both with and without an externally imposed azimuthal flow show that the growth of a CITM is observed for a relatively small flow velocity, and that a CITM can be a candidate mechanism for the intermittent eruption of the current in a tokamak, whether the growth is under influences of non-ideal MHD effects or not. Here, the growth can be suppressed when an externally imposed azimuthal flow is sufficiently large to cause a finite radial displacement because such a displacement prevents a transition from an interchange to a tearing mode that is the key process of a CITM. Furthermore, it is also found that a stripe pattern similar to so-called streamers is formed in the course of the CITM growth. The basic nature of this streamer-like structure is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Advances in 3D transient plasma dynamics and control through MHD and hybrid fluid-kinetic simulations with JOREK

Transient phenomena and their control are of high relevance in magnetic confinement fusion plasmas to guarantee a stable and safe plasma operation. Interpretative simulations can maximize the insights gained from experiments on present machines and predictive simulations can help in the preparation of design, mitigation techniques and operational scenarios for future devices. In this article, we provide an overview of recent advances and novel scientific results obtained with the 3D non-linear hybrid fluid-kinetic code JOREK, covering physics of plasma transients from the core to the scrape-off layer (SOL) both for tokamak and stellarator devices. Substantial progress was made in the physics understanding, model validation with experiments and experiment interpretation, thus, giving confidence for predictions to devices like DTT, ITER and DEMO. The topics addressed comprise a wide range: the edge physics of new operation scenarios and edge localized mode suppression; major disruptions with a focus on runaway electrons and vertical displacement events as well as disruption mitigation by shattered pellet injection; the physics mechanisms and operational limits of the flux pumping regime for sawtooth control; MHD limits of stellarators and work towards incorporating advanced edge/SOL/exhaust dynamics; continuing improvements of the code for more efficient hybrid simulations on conventional and accelerated high performance computing architectures.

disruptions↗
DOE OSTI2024

Λ polarization from vortex rings as the medium response for jet thermalization

We performed a systematic study on the formation of vorticity rings as the process for jet thermalization in the medium created in high-energy nuclear collisions. In this work, we expanded our previous analysis to a more realistic framework by considering noncentral events and fluctuations in the initial condition. We simulate the formation and evolution of the flow vortex structure in a relativistic viscous hydrodynamic model and study the sensitivity of the proposed “ring observable” (ℛ$^{𝑡}_{Λ}$) that can be measured experimentally through the polarization of Λ hyperons. We show that this observable is robust with respect to fluctuating initial conditions to capture the jet-induced vortex flow signal and further study its dependence on different model parameters, such as the jet's velocity, position, the fluid's shear viscosity, and the collision centrality. The proposed observable is associated with the formation of vorticity in a quark-gluon plasma, showing that the measurement of particle polarization can be a powerful tool to probe different properties of jet-medium interactions and to understand better the polarization induced by the transverse and longitudinal expansions of the medium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Time dependent simulation of cosmic-ray shocks including Alfven transport

Time evolution of plane, cosmic-ray modified shocks was simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven wave energy transport equations. Nonlinear feedback from the cosmic-rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the time scale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic-rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS1994

Time-dependent simulation of cosmic-ray shocks, including Alfven transport

Time evolution of plane, cosmic-ray modified shocks has been simulated numerically for the case with parallel magnetic fields. Computations were done in a 'three-fluid' dynamical model incorporating cosmic-ray and Alfven-wave energy transport equations. Nonlinear feedback from the cosmic rays and Alfven waves is included in the equation of motion for the underlying plasma, as is the finite propagation speed and energy dissipation of the Alfven waves. Exploratory results confirm earlier, steady state analyses that found these Alfven transport effects to be potentially important when the upstream Alfven speed and gas sound speeds are comparable. As noted earlier, Alfven transport effects tend to reduce the transfer of energy through a shock from gas to energetic particles. These studies show as well that the timescale for modification of the shock is altered in nonlinear ways. It is clear, however, that the consequences of Alfven transport are strongly model dependent and that both advection of cosmic rays by the waves and dissipation of wave energy in the plasma will be important to model correctly when quantitative results are needed. Comparison is made between simulations based on a constant diffusion coefficient and more realistic diffusion models allowing the diffusion coefficient to vary in response to changes in Alfven wave intensity. No really substantive differences were found between them.

Jones, T. W.↗
NASA NTRS2012

Extending Magnetohydrodynamics to the Slow Dynamics of Collisionless Plasmas

A fluid approach aimed to provide a consistent description of the slow dynamics of a collisionless plasma, is presented. In this regime, both Landau damping and finite Larmor radius effects cannot be ignored. Two models are discussed; one retains the dynamics at sub-ionic scales, while the other is restricted to scales larger than the ion gyroscale. Special attention is paid to the capability of these approaches to accurately reproduce the properties of linear waves that are known to play an important role, for example, in the small-scale dynamics of solar wind turbulence.

Passot, T.↗
NASA NTRS2003

Fundamentals of the Plasma Sail Concept: MHD and Kinetic Studies

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000] predicts that a 15-km standoff distance (or 20-km cross-sectional dimension) of the magnetic bubble will provide for sufficient momentum transfer from the solar wind to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Such velocities will enable travel out of the solar system in period of about seven years-almost an order of magnitude improvement over present chemical based propulsion systems. However, for the parameters of the simulation of Winglee et al. [2000], a fluid model for the interaction of M2P2 with the solar wind is not valid. It is assumed in the MHD fluid model, normally applied to planetary magnetospheres, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the solar wind. In the case of M2P2, the size of the magnetic bubble is actually less than or, comparable to, the scale of these characteristic parameters. Therefore, a kinetic approach, which addresses the small-scale physical mechanisms, must be used. We have adopted a two-component approach to determining a preliminary estimate of the momentum transfer to the plasma sail. The first component is a self-consistent MHD simulation of the small-scale expansion phase of the magnetic bubble. The fluid treatment is valid to roughly 5 km from the source and the steady-state MHD solution at the 5 km boundary was then used as initial conditions for the hybrid simulation. The hybrid simulations showed that the momentum transfer to the innermost regions of the plasma sail is negligible.

Khazanov, G.↗