Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA FLUID EQUATIONS”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Simulations of drift resistive ballooning L-mode turbulence in the edge plasma of the DIII-D tokamak

Results from simulations of electromagnetic drift-resistive ballooning turbulence for tokamak edge turbulence in realistic single-null geometry are reported. The calculations are undertaken with the BOUT three-dimensional fluid code that solves Braginskii-based fluid equations [X. Q. Xu and R. H. Cohen, Contrib. Plasma Phys. 36, 158 (1998)]. The simulation setup models L-mode edge plasma parameters in the actual magnetic geometry of the DIII-D tokamak [J. L. Luxon et al., Fusion Sci. Technol. 48, 807 (2002)]. The computations track the development of drift-resistive ballooning turbulence in the edge region to saturation. Fluctuation amplitudes, fluctuation spectra, and particle and thermal fluxes are compared to experimental data near the outer midplane from Langmuir probe and beam-emission-spectroscopy for a few well-characterized L-mode discharges in DIII-D. The simulations are comprised of a suite of runs in which the physics model is varied to include more fluid fields and physics terms. The simulations yield results for fluctuation amplitudes, correlation lengths, particle and energy fluxes, and diffusivities that agree with measurements within an order of magnitude and within factors of 2 or better for some of the data. The agreement of the simulations with the experimental measurements varies with respect to including more physics in the model equations within the suite of models investigated. The simulations show stabilizing effects of sheared E × B poloidal rotation (imposed zonal flow) and of lower edge electron temperature and density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modelling of nanometer scale dust grains in tokamak

Dust poses a serious threat to tokamak operation and safety. It is important to study the behaviour of dust grains under tokamak's discharge conditions, which depends heavily on their size and charge. Existing simulations mainly address issues on dust grains with radii larger than 1 μm, in which case, the drift effect due to electromagnetic fields can be safely ignored. For nanometer scale dust grains, however, the drift effect becomes significant and a new model based on guiding-centre system needs to be established. In this work, the NDS has been done under BOUT++ framework. The simulation contains two parts. Part one, NDS evaluates the charging and ablation processes of the dust grains. In the second part, the guiding-centre orbits of dust particles are tracked in tokamak plasmas, whose parameters are obtained from BOUT++, a highly desirable C++ code package for performing parallel plasma fluid simulations with an arbitrary number of equations in 3D curvilinear coordinates. The orbit of nanodust dynamics is described by guiding centre equations for simplicity, and these equations are numerically solved by conventional fourth-order Runge Kutta method. Simulations provide results such as trajectories and evolutions of dust particles with different sizes and velocities for different tokamak geometries. Results show tungsten dust grains with a radius of a few nanometers launched from outer midplane will oscillate before totally ablated in C-Mod. The oscillation in this case is driven by the ion drag force. Larger Nanodust with a radius of 100 nm, on the contrary, cannot be completely constrained by the electromagnetic field. In conclusion, the high plasma temperature and density in the seperatrix region causes severe dust ablation, resulting in total ablation within several ms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Moment-Fourier approach to ion parallel fluid closures and transport for a toroidally confined plasma

A general method of solving the drift kinetic equation is developed for an axisymmetric magnetic field. Expanding a distribution function in general moments, a set of ordinary differential equations is obtained. Successively expanding the moments and magnetic-field involved quantities in Fourier series, a set of linear algebraic equations is obtained. The set of full (Maxwellian and non-Maxwellian) moment equations is solved to express the first-order density, temperature, and flow velocity in terms of radial gradients of the zeroth-order pressure and temperature. Closure relations that connect parallel heat flux density and viscosity to the radial gradients and parallel gradients of temperature and flow velocity are also obtained by solving the non-Maxwellian moment equations. The closure relations combined with the linearized fluid equations reproduce the same solution obtained directly from the full moment equations. Furthermore, the method can be generalized to derive closures and transport for an electron-ion plasma and a multi-ion plasma in a general magnetic field.

neoclassical transport↗
NASA NTRS2000

Ion Dynamics Model for Collisionless Radio Frequency Sheaths

Full scale reactor model based on fluid equations is widely used to analyze high density plasma reactors. It is well known that the submillimeter scale sheath in front of a biased electrode supporting the wafer is difficult to resolve in numerical simulations, and the common practice is to use results for electric field from some form of analytical sheath model as boundary conditions for full scale reactor simulation. There are several sheath models in the literature ranging from Child's law to a recent unified sheath model [P. A. Miller and M. E. Riley, J. Appl. Phys. 82, 3689 (1997)l. In the present work, the cold ion fluid equations in the radio frequency sheath are solved numerically to show that the spatiotemporal variation of ion flux inside the sheath, commonly ignored in analytical models, is important in determining the electric field and ion energy at the electrode. Consequently, a semianalytical model that includes the spatiotemporal variation of ion flux is developed for use as boundary condition in reactor simulations. This semianalytical model is shown to yield results for sheath properties in close agreement with numerical solutions.

Bose, Deepak↗
NASA NTRS1986

Rail gun performance and plasma characteristics due to wall ablation

The experiment of Bauer, et al. (1982) is analyzed by considering wall ablation and viscous drag in the plasma. Plasma characteristics are evaluated through a simple fluid-mechanical analysis considering only wall ablation. By equating the energy dissipated in the plasma with the radiation heat loss, the average properties of the plasma are determined as a function of time.

Ray, P. K.↗
NASA NTRS1987

Arc-driven rail accelerator research

Arc-driven rail accelerator research is analyzed by considering wall ablation and viscous drag in the plasma. Plasma characteristics are evaluated through a simple fluid-mechanical analysis considering only wall ablation. By equating the energy dissipated in the plasma with the radiation heat loss, the average properties of the plasma are determined as a function of time and rate of ablation. Locations of two simultaneously accelerating arcs were determined by optical and magnetic probes and fron streak camera photographs. All three measurements provide consistent results.

Ray, Pradosh K.↗
DOE OSTI2017

Nonlinear fluid simulation study of stimulated Raman and Brillouin scatterings in shock ignition

Here, we developed a new nonlinear fluid laser-plasma-instability code ( FLAME ) using a multi-fluid plasma model combined with full electromagnetic wave equations. The completed one-dimensional version of FLAME was used to study laser-plasma instabilities in shock ignition. The simulations results showed that absolute stimulated Raman scattering (SRS) modes growing near the quarter-critical surface were saturated by Langmuir-wave Decay Instabilities (LDI) and pump depletion. The ion-acoustic waves from LDI acted as seeds of Stimulated Brillouin Scattering, which displayed a bursting pattern and caused strong pump depletion. Re-scattering of SRS at the 1/16th-critical surface was also observed in a high temperature case. These results largely agreed with the corresponding Particle-in-Cell simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Generation of excited species in a streamer discharge

At or near atmospheric pressure, most transient discharges, particularly in molecular gases or gas mixture containing molecular gases result in a space-charge dominated transport called a streamer discharge. The excited species generation in such discharges forms the basis for plasma chemistry in most technological applications. In this paper we simulate the propagation of streamers in atmospheric pressure N 2 to understand the energy partitioning in the formation of various excited species and compare the results to a uniform Townsend discharge. The model is fully two-dimensional with azimuthal symmetry. Results show a significantly larger fraction of the energy goes into vibrational excitation of the N 2 ground state in a streamer-type discharge in comparison to a Townsend discharge. For lower applied voltages, the anode-directed (negative) steamer is slightly more efficient in channeling energy for excited species production in comparison to a cathode-directed (positive) streamer. Near 70% overvoltage both types of streamers show very similar energy partitioning but quite different from a Townsend discharge.

42 ENGINEERING↗
DOE OSTI2018

Turbulence in Magnetized Pair Plasmas

Alfvénic-type turbulence in strongly magnetized, low-beta pair plasmas is investigated. A coupled set of equations for the evolution of the magnetic and flow potentials are derived, covering both fluid and kinetic scales. In the fluid (magnetohydrodynamic) range those equations are the same as for electron–ion plasmas, so turbulence at those scales is expected to be of the Alfvénic nature, exhibiting critical balance, dynamic alignment, and transition to a tearing-mediated regime at small scales. The critical scale at which a transition to a tearing-mediated range occurs is derived, and the spectral slope in that range is predicted to be ${k}_{\perp }^{-8/3}$ (or ${k}_{\perp }^{-3}$, depending on details of the reconnecting configuration assumed). At scales below the electron (and positron) skin depth, it is argued that turbulence is dictated by a cascade of the inertial Alfvén wave, which we show to result in the magnetic energy spectrum $\propto {k}_{\perp }^{-11/3}$.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Potential around a dust grain in collisional plasma

The ion neutral collision can lead to interesting phenomena in dust charging, totally different from the expectations based on the traditional orbit motion limited theory. The potential around a dust grain is investigated for the collisional plasma considering the presence of ion neutral collisions. Fluid equations are solved for the one dimensional radial coordinate. It is observed that with the gradual increase in ion neutral collision, the potential structure around the dust grain changes its shape and is different from the usual Debye-Hückel potential. The shift however starts from a certain value of ion neutral collision and the electron-ion density varies accordingly. The potential variation is interesting and reconfirms the fact that there exists a region of attraction for negative charges. The collision modeling is done for the full range of plasma, i.e., considering the bulk and the sheath jointly. The potential variation with collision is also shown explicitly and the variation is found to cope up with the earlier observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Mixing of equations of state for xenon-deuterium using density functional theory

We report on a theoretical study of equation of state (EOS) properties of fluid and dense plasma mixtures of xenon and deuterium to explore and illustrate the basic physics of the mixing of a light element with a heavy element. Accurate EOS models are crucial to achieve high-fidelity hydrodynamics simulations of many high-energy-density phenomena, for example inertial confinement fusion and strong shock waves. While the EOS is often tabulated for separate species, the equation of state for arbitrary mixtures is generally not available, requiring properties of the mixture to be approximated by combining physical properties of the pure systems. Density functional theory (DFT) at elevated-temperature is used to assess the thermodynamics of the xenon-deuterium mixture at different mass ratios. The DFT simulations are unbiased as to elemental species and therefore provide comparable accuracy when describing total energies, pressures, and other physical properties of mixtures as they do for pure systems. The study focuses on addressing the accuracy of different mixing rules in the temperature range 1000-40 000 K for pressures between 100 and 600 GPa (1-6 Mbar), thus, including the challenging warm dense matter regime of the phase diagram. We find that a mix rule taking into account pressure equilibration between the two species performs very well over the investigated range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A scaling law of the neutral penetration length and Balmer- α wing shape in high-temperature plasmas

Hydrogen atoms penetrating deep inside high-temperature magnetically confined plasmas by repetitive charge-exchange collisions result in a particle source of the plasma, which affects the plasma transport significantly. In this paper, we present an approximate solution of the fluid equations for neutral transport and an analytical representation of the neutral penetration length, in a simplified plasma geometry. This analysis predicts a power-law decay in the Balmer-α line wings which reflects the velocity distribution of the neutral atoms, with the power-law exponent analytically represented as well. These scaling laws are compared with a simple Monte–Carlo simulation and spectroscopic observations of Large Helical Device plasmas. Since the Balmer-α line wings are experimentally accessible, our formulation opens the possibility to quickly estimate the neutral penetration length from spectroscopic observations.

neutral opacity↗
NASA NTRS1972

Numerical study of strong plasma shock waves produced in an electromagnetic shock tube.

A set of two-fluid Navier-Stokes equations with classical physical transport coefficients is used to compute the evolution and structure of collisional plasma shock waves in an electromagnetic shock tube. Shock speeds up to 200 cm/microsec and shocked-plasma temperatures of the order of kilovolts are studied. A strong transverse bias magnetic field is employed, which significantly alters the size and shape of the shock profiles, when compared with the zero bias field case. The wave structure is different from that of a two-fluid gasdynamic shock without transverse magnetic field. Especially significant is the effect of the small ion Larmor radius in reducing the gasdynamic shock thickness by at least an order of magnitude in the transverse shock case, which permits collisional shocks to have thicknesses much smaller than the post-shock mean free path. These collisional shock waves produce very hot ions, the ion temperature increasing with shock speed, although the formation time and distances also increase substantially with shock speed.

Schneider, S. H.↗
DOE OSTI2017

Color instabilities in the quark–gluon plasma

When the quark–gluon plasma (QGP) – a system of deconfined quarks and gluons – is in a nonequilibrium state, it is usually unstable with respect to color collective modes. The instabilities, which are expected to strongly influence dynamics of the QGP produced in relativistic heavy-ion collisions, are extensively discussed under the assumption that the plasma is weakly coupled. Here, we begin by presenting the theoretical approaches to study the QGP, which include: field theory methods based on the Keldysh–Schwinger formalism, classical and quantum kinetic theories, and fluid techniques. The dispersion equations, which give the spectrum of plasma collective excitations, are analyzed in detail. We pay particular attention to a momentum distribution of plasma constituents which is obtained by deforming an isotropic momentum distribution. Mechanisms of chromoelectric and chromomagnetic instabilities are explained in terms of elementary physics. The Nyquist analysis, which allows one to determine the number of solutions of a dispersion equation without explicitly solving it, and stability criteria are also discussed. We then review various numerical approaches – purely classical or quantum – to simulate the temporal evolution of an unstable quark–gluon plasma. The dynamical role of instabilities in the processes of plasma equilibration is analyzed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

Longitudinal singular response of dusty plasma medium in weak and strong coupling limits

The longitudinal response of a dusty plasma medium in both weak and strong coupling limits has been investigated in detail using analytic as well as numerical techniques. In particular, studies on singular response of the medium have been specifically investigated here. A proper Galilean invariant form of the generalized hydrodynamic fluid model has been adopted for the description of the dusty plasma medium. For weak non-linear response, analytic reductive perturbative approach has been adopted. It is well known that in the weak coupling regime for the dusty plasma medium, such an analysis leads to the Korteweg-de Vries equation (KdV) equation and predicts the existence of localized smooth soliton solutions. We show that the strongly coupled dust fluid with the correct Galilean invariant form does not follow the KdV paradigm. Instead, it reduces to the form of Hunter-Saxton equation, which does not permit soliton solutions. The system in this case displays singular response with both conservative as well as dissipative attributes. At arbitrary high amplitudes, the existence and spontaneous formation of sharply peaked cusp structures in both weak and strong coupling regimes has been demonstrated numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Symmetries of a reduced fluid-gyrokinetic system

Symmetries of a fluid-gyrokinetic model are investigated using Lie group techniques. Specifically, the nonlinear system constructed, which combines nonlinear fluid equations with a drift-kinetic description of parallel electron dynamics, is studied. Significantly, this model is fully gyrokinetic, allowing for arbitrary $k_{\bot }\unicode[STIX]{x1D70C}_{i}$ , where $k_{\bot }$ is the perpendicular wave vector of the fluctuations and $\unicode[STIX]{x1D70C}_{i}$ the ion gyroradius. The model includes integral operators corresponding to gyroaveraging as well as the moment equations relating fluid variables to the kinetic distribution function. Here, a large variety of exact symmetries is uncovered, some of which have unexpected form. Using these results, new nonlinear solutions are constructed, including a helical generalization of the Chapman–Kendall solution for a collapsing current sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fueling limits in a cylindrical viscosity-limited reactor

Recently, a method to achieve a “natural hot-ion mode” was suggested by utilizing ion viscous heating in a rotating plasma with a fixed boundary. We explore the steady-state solution to the Braginskii equations and find the parameter regime in which a significant temperature difference between ions and electrons can be sustained in a driven steady state. Here, the threshold for this effect occurs at ρ i ≳ 0.1R. An analytic, leading order low flow solution is obtained, and a numerical, moderate Mach number M ≲ 2 is investigated. The limitation is found to be at moderate Mach numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the merit of hot ion mode for tearing mode stabilization

The stabilization of tearing modes with rf driven current benefits from the cooperative feedback loop between rf power deposition and electron temperature within the island. This effect, termed rf current condensation, can greatly enhance and localize the current driven within magnetic islands. It has previously been shown that the condensation effect opens the possibility of passive stabilization with broad rf profiles, as would be typical of LHCD for steady state operation. Here, we show that this self-healing effect can be dramatically amplified by operation in a hot ion mode, due to the additional electron heat source provided by the hotter ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗