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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Photon polarizability and its effect on the dispersion of plasma waves

High-frequency photons travelling in plasma exhibit a linear polarizability that can influence the dispersion of linear plasma waves. We present a detailed calculation of this effect for Langmuir waves as a characteristic example. Here, two alternative formulations are given. In the first formulation, we calculate the modified dispersion of Langmuir waves by solving the governing equations for the electron fluid, where the photon contribution enters as a ponderomotive force. In the second formulation, we provide a derivation based on the photon polarizability. Then, the calculation of ponderomotive forces is not needed, and the result is more general.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The effect of a dust size distribution on electrostatic sheaths in unmagnetized dusty plasmas

In this work, the structure of plasma sheaths in presence of dust particles with different sizes is investigated numerically in a multifluid framework, where the dust size distribution is modeled by Gauss' law. For this, we have established a 1D, stationary, unmagnetized, and weakly collisional electronegative dusty plasma sheath model. The electrons and negative ions are considered in a local thermodynamic equilibrium, therefore, described by a Boltzmann distribution. On the other hand, positive ions and dust grains are described by fluid equations. The charging process is described by the orbit motion limited model. It is shown that taking into account dust grains with different sizes reduces considerably the sheath thickness. The behavior of dust surface potential is not affected, but the dust charge number is reduced, as well as the electrostatic force. It results in a decrease of layered structure. The presence of negative ions makes the structure of the electrostatic potential more oscillatory. The other physical parameters are also analyzed and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Role of the edge stochastic layer in density pump-out by resonant magnetic perturbations

The role of the edge stochastic layer on particle transport is addressed in DIII-D plasmas with applied resonant magnetic perturbations (RMPs) causing density pump-out (a term used to denote the density reduction at the top of the pedestal due to the applied RMPs). Using an analytical model that adds the stochastic parallel transport of electrons (Harvey et al 1981 Phys. Rev. Lett. 47 102) to the fluid equations, the ambipolar radial electric field and particle flux are calculated simultaneously. In this model the nonambipolar electron particle flux, driven by the stochastic magnetic field, is predominantly balanced by the nonambipolar perpendicular ion flux, driven by toroidal viscosity, across a narrow stochastic layer of order two percent of the plasma radius (0.98 < ψn < 1). The model reproduces the level of density pump-out and its dependence on RMP amplitude near the separatrix for three plasma discharges for which the pedestal foot is at medium to high collisionality $(2\leqslant {\nu }_{\text{sep}}^{\ast }\leqslant 30)$. In this work, the experimentally observed and numerically tested inverse density dependence of density pump-out (Hu et al 2020 Nucl. Fusion 60 076001) is accurately captured by the model: an increase in collisionality with density in the pedestal foot results in a decrease in stochastic diffusivity, and hence a decrease in the level of pump-out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Major disruptions, inverse cascades, and the Strauss equations

Current-carrying plasmas in a strong d.c. magnetic field are subject to violent disruptions above certain thresholds. At present difficult to verify, explanations are typically sought in terms of 'tearing modes'. An alternative explanation is in terms of inverse magnetic helicity cascades, generated from a variety of possible sources of small-scale MHD turbulence. Strongly anisotropic MHD plasmas may be described by the Strauss equations. Indications of turbulent inverse cascade behavior for the Strauss equations are sought, in parallel with earlier examples from MHD and fluid mechanics.

Montgomery, D.↗
DOE OSTI2019

Study of up–down poloidal density asymmetry of high- $\mathrm{Z}$ impurities with the new impurity version of XGCa

Addition of multispecies impurity ions to the total-f gyrokinetic particle-in-cell code XGCa is reported, including a cross-verification of neoclassical physics against the NEO code. This new version of the neoclassical gyrokinetic code XGCa is used to benchmark and confirm the previous reduced-equation-based prediction that high- $\mathrm{Z}$ impurity particles in the Pfirsch–Schlüter regime can exhibit a significant level of up–down poloidal asymmetry, through the large parallel friction force, and thus influence the radial plasma transport significantly. The study is performed in a plasma with weak toroidal rotation. In comparison, when the impurity particles are in the plateau regime, the up–down poloidal asymmetry becomes weak, with the parallel friction force becoming weaker than the parallel viscous force. It is also found that the linearization of the perturbed distribution function, based on the small poloidal asymmetry assumption, can become invalid. Finally, using the numerical data from XGCa, each term in the parallel fluid force-balance equation have been analysed to find that both the main ions and the electrons respond to the poloidal potential variation adiabatically when the high- $\mathrm{Z}$ tungsten possesses large poloidal variation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Energy and momentum flow in electromagnetic fields and plasma

The energy momentum tensor for a perfect fluid in a magnetic field is used to predict the momentum density, energy density, momentum flow, and energy flow of the fluid and the electromagnetic field. It is shown that taking the momentum flow from the energy momentum tensor, rather than starting with differential magnetohydrodynamic equations, can produce more accurate results on the basis of magnetic field data. It is suggested that the use of the energy momentum tensor has the potential for application to analysis of data from the more dynamic regions of the solar system, such as the plasma boundaries of Venus, the Jovian ionosphere, and the terrestrial magnetopause.

Parish, J. L.↗
DOE OSTI2015

Extension of the flow-rate-of-strain tensor formulation of plasma rotation theory to non-axisymmetric tokamaks

A systematic formalism for the calculation of rotation in non-axisymmetric tokamaks with 3D magnetic fields is described. The Braginskii Ωτ-ordered viscous stress tensor formalism, generalized to accommodate non-axisymmetric 3D magnetic fields in general toroidal flux surface geometry, and the resulting fluid moment equations provide a systematic formalism for the calculation of toroidal and poloidal rotation and radial ion flow in tokamaks in the presence of various non-axisymmetric “neoclassical toroidal viscosity” mechanisms. The relation among rotation velocities, radial ion particle flux, ion orbit loss, and radial electric field is discussed, and the possibility of controlling these quantities by producing externally controllable toroidal and/or poloidal currents in the edge plasma for this purpose is suggested for future investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Mixing of the Interstellar and Solar Plasmas at the Heliospheric Interface

From the ideal MHD perspective, the heliopause is a tangential discontinuity that separates the solar wind plasma from the local interstellar medium plasma. There are physical processes, however, that make the heliopause permeable. They can be subdivided into kinetic and MHD categories. Kinetic processes occur on small length and time scales, and cannot be resolved with MHD equations. On the other hand, MHD instabilities of the heliopause have much larger scales and can be easily observed by spacecraft. The heliopause may also be a subject of magnetic reconnection. In this paper, we discuss mechanisms of plasma mixing at the heliopause in the context of Voyager 1 observations. Numerical results are obtained with a Multi-Scale Fluid-Kinetic Simulation Suite (MS-FLUKSS), which is a package of numerical codes capable of performing adaptive mesh refinement simulations of complex plasma flows in the presence of discontinuities and charge exchange between ions and neutral atoms. The flow of the ionized component is described with the ideal MHD equations, while the transport of atoms is governed either by the Boltzmann equation or multiple Euler gas dynamics equations. The code can also treat nonthermal ions and turbulence produced by them.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2-D, two-temperature, single fluid magnetohydrodynamic code which incorporates classical plasma transport coefficients and Hall effects was developed to predict steady-state, self-field magnetoplasmadynamic (MPD) thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects, the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries, assuming a fully ionized argon propellant. A maximum predicted specific impulse of 1680 s is obtained for a thruster with an anode radius of 2.5 cm, a cathode radius of 0.5 cm, and equal electrode lengths of 2.5 cm. A scaling relation is developed to predict, within limits, the onset of cylindrical, self-field thruster instability as a function of geometry and operating condition.

Lapointe, M.↗
NASA NTRS1992

Numerical simulation of geometric scale effects in cylindrical self-field MPD thrusters

A 2D, two-temperature, single fluid MHD code which incorporates classical plasma transport coefficients and Hall effects has been developed to predict steady-state, self-field MPD thruster performance. The governing equations and numerical methods of solution are outlined and discussed. Experimental comparisons are used to validate model predictions. The model accurately predicts thrust and reproduces trends in the discharge voltage for discharge currents below experimentally measured onset values. However, because the model does not include electrode effects the calculated voltage drops are significantly lower than experimentally measured values. Predictions of thrust and flow efficiency are made for a matrix of fifteen cylindrical thruster geometries assuming a fully ionized argon propellant.

Lapointe, M.↗
DOE OSTI2011

Gyrokinetic simulations of the tearing instability

Linear gyrokinetic simulations covering the collisional-collisionless transitional regime of the tearing instability are performed. It is shown that the growth rate scaling with collisionality agrees well with that predicted by a two-fluid theory for a low plasma beta case in which ion kinetic dynamics are negligible. Electron wave-particle interactions (Landau damping), finite Larmor radius, and other kinetic effects invalidate the fluid theory in the collisionless regime, in which a general non-polytropic equation of state for pressure (temperature) perturbations should be considered. We also vary the ratio of the background ion to electron temperatures and show that the scalings expected from existing calculations can be recovered, but only in the limit of very low beta.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Nonlinear propagation of coherent electromagnetic waves in a dense magnetized plasma

We present an investigation of the nonlinear propagation of high-frequency coherent electromagnetic waves in a uniform quantum magnetoplasma. Specifically, we consider nonlinear couplings of right-hand circularly polarized electromagnetic-electron-cyclotron (CPEM-EC) waves with dispersive shear Alfven (DSA) and dispersive compressional Alfven (DCA) perturbations in plasmas composed of degenerate electron fluids and non-degenerate ion fluids. Such interactions lead to amplitude modulation of the CPEM-EC wave packets, the dynamics of which is governed by a three-dimensional nonlinear Schroedinger equation (NLSE) with the frequency shift arising from the relativistic electron mass increase in the CPEM-EC fields and density perturbations associated with the DSA and DCA perturbations. Accounting for the electromagnetic and quantum forces, we derive the evolution equation for the DSA and DCA waves in the presence of the magnetic field-aligned ponderomotive force of the CPEM-EC waves. The NLSE and the driven DSA and DCA equations are then used to investigate the modulational instability. The relevance of our investigation to laser-plasma interaction experiments and the cores of white dwarf stars is pointed out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)], in which low-density plasmas ( 10 19 to 10 20 cm – 3 ) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, in this work we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Fokker Planck and Krook theory of energetic electron transport in a laser produced plasma

Various laser plasma instabilities, such as the two plasma decay instability and the stimulated Raman scatter instability, produce large quantities of energetic electrons. How these electrons are transported and heat the plasma are crucial questions for laser fusion. This paper works out a Fokker Planck and Krook theory for such transport and heating. The result is a set of equations, for which one can find a simple asymptotic approximation for the solution, for the Fokker Planck case, and an exact solution for the Krook case. These solutions are evaluated and compared with one another. They give rise to expressions for the spatially dependent heating of the background plasma, as a function of the instantaneous laser and plasma parameters, in either planar or spherical geometry. These formulas are simple, universal (depending weakly only on the single parameter Z, the charge state), and can be easily be incorporated into a fluid simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Viscous effects on the Rayleigh-Taylor instability with background temperature gradient

Here we studied the growth rate of the compressible Rayleigh-Taylor instability in the presence of a background temperature gradient, Θ, using a normal mode analysis. The effect of Θ variation is examined for three interface types corresponding to the combinations of the viscous properties of the fluids (inviscid-inviscid, viscous-viscous, and viscous-inviscid) at different Atwood numbers, At, and when at least one of the fluids' viscosity is non-zero, as a function of the Grashof number. For the general case, the resulting ordinary differential equations are solved numerically; however, dispersion relations for the growth rate are presented for several limiting cases. An analytical solution is found for the inviscid-inviscid interface and the corresponding dispersion equation for the growth rate is obtained in the limit of large Θ. For the viscous-inviscid case, a dispersion relation is derived in the incompressible limit and Θ=0. Compared to Θ=0 case, the role of Θ<0 (hotter light fluid) is destabilizing and becomes stabilizing when Θ>0 (colder light fluid). The most pronounced effect of Θ ≠ 0 is found at low At and/or at large perturbation wavelengths relative to the domain size for all interface types. On the other hand, at small perturbation wavelengths relative to the domain size, the growth rate for the Θ<0 case exceeds the infinite domain incompressible constant density result. The results are applied to two practical examples, using sets of parameters relevant to Inertial Confinement Fusion coasting stage and solar corona plumes. The role of viscosity on the growth rate reduction is discussed together with highlighting the range of wavenumbers most affected by viscosity. The viscous effects further increase in the presence of background temperature gradient, when the viscosity is temperature dependent.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗