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

Drekar v.2.0

Drekar is an application code that solves partial differential equations for fluids that can be optionally coupled to electromagnetics. Drekar solves low-mach compressible and incompressible computational fluid dynamics (CFD), compressible and incompressible resistive magnetohydrodynamics (MHD), and multiple species plasmas interacting with electromagnetic fields. Drekar discretization technology includes continuous and discontinuous finite element formulations, stabilized finite element formulations, mixed integration finite element bases (nodal, edge, face, volume) and an initial arbitrary Lagrangian Eulerian (ALE) capability. Drekar contains the implementation of the discretized physics and leverages the open source Trilinos project for both parallel solver capabilities and general finite element discretization tools. The code will be released open source under a BSD license. The code is used for fundamental research for simulation of fluids and plasmas on high performance computing environments.

Seefeldt, Ben↗
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.↗
NASA NTRS1985

Hybrid simulation codes with application to shocks and upstream waves

Hybrid codes in which part of the plasma is represented as particles and the rest as a fluid are discussed. In the past few years such codes with particle ions and massless, fluid electrons have been applied to space plasmas, especially to collisionless shocks. All of these simulation codes are one-dimensional and similar in structure, except for how the field equations are solved. The various approaches that are used (resistive Ohm's law, predictor-corrector, Hamiltonian) are described in detail and results from the various codes are compared with examples taken from collisionless shocks and low frequency wave phenomena upstream of shocks.

Winske, D.↗
DOE OSTI2014

From Boltzmann equations to steady wall velocities

By means of a relativistic microscopic approach we calculate the expansion velocity of bubbles generated during a first-order electroweak phase transition. In particular, we use the gradient expansion of the Kadanoff-Baym equations to set up the fluid system. This turns out to be equivalent to the one found in the semi-classical approach in the non-relativistic limit. Finally, by including hydrodynamic deflagration effects and solving the Higgs equations of motion in the fluid, we determine velocity and thickness of the bubble walls. Our findings are compared with phenomenological models of wall velocities. As illustrative examples, we apply these results to three theories providing first-order phase transitions with a particle content in the thermal plasma that resembles the Standard Model.

79 ASTRONOMY AND ASTROPHYSICS↗
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↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Study of transport properties with relativistic ponderomotive effect in two-electron temperature plasma

In the present paper we make an analytical investigation to study transport properties with relativistic ponderomotive effect in two-electron temperature plasma. Using fluid model the two-electron temperature are introduced through relativistic ponderomotive force for the transportation of two species of electrons. Applying WKB and paraxial ray approximation the nonlinear dielectric constant and self-focusing equation is evaluated and analyzed with experimental relevance. Numerical calculations are made for different concentration of electron density (10{sup 19}−10{sup 21} per cm{sup 3}) at arbitrary values of laser intensity in the range 10{sup 18}−10{sup 21} W/cm{sup 2}. For a minimum radius depending on the initial conditions it is oscillating between a minimum and maximum value. The hot electrons leading to the increase of the on-axis transportation and favorable effect on relativistic self-focusing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Viscous quark-gluon plasma model through fluid QCD approach

A Lagrangian density for viscous quark-gluon plasma has been constructed within the fluid-like QCD framework. Gauge symmetry is preserved for all terms inside the Lagrangian, except for the viscous term. The transition mechanism from point particle field to fluid field, and vice versa, are discussed. The energy momentum tensor that is relevant to the gluonic plasma having the nature of fluid bulk of gluon sea is derived within the model. By imposing conservation law in the energy momentum tensor, shear viscosity appears as extractable from the equation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Pressure tensor in the presence of velocity shear: Stationary solutions and self-consistent equilibria

Observations and numerical simulations of laboratory and space plasmas in almost collisionless regimes reveal anisotropic and non-gyrotropic particle distribution functions. We investigate how such states can persist in the presence of a sheared flow. We focus our attention on the pressure tensor equation in a magnetized plasma and derive analytical self-consistent plasma equilibria which exhibit a novel asymmetry with respect to the magnetic field direction. These results are relevant for investigating, within fluid models that retain the full pressure tensor dynamics, plasma configurations where a background shear flow is present.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Linear and nonlinear heavy ion-acoustic waves in a strongly coupled plasma

A theoretical study on the propagation of linear and nonlinear heavy ion-acoustic (HIA) waves in an unmagnetized, collisionless, strongly coupled plasma system has been carried out. The plasma system is assumed to contain adiabatic positively charged inertial heavy ion fluids, nonextensive distributed electrons, and Maxwellian light ions. The normal mode analysis is used to study the linear behaviour. On the other hand, the well-known reductive perturbation technique is used to derive the nonlinear dynamical equations, namely, Burgers equation and Korteweg-de Vries (K-dV) equation. They are also numerically analyzed in order to investigate the basic features of shock and solitary waves. The adiabatic effects on the HIA shock and solitary waves propagating in such a strongly coupled plasma are taken into account. It has been observed that the roles of the adiabatic positively charged heavy ions, nonextensivity of electrons, and other plasma parameters arised in this investigation have significantly modified the basic features (viz., polarity, amplitude, width, etc.) of the HIA solitary/shock waves. The findings of our results obtained from this theoretical investigation may be useful in understanding the linear as well as nonlinear phenomena associated with the HIA waves both in space and laboratory plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A self-consistent two-fluid model of a magnetized plasma-wall transition

A self-consistent one-dimensional two-fluid model of the magnetized plasma-wall transition is presented. The model includes magnetic field, elastic collisions between ions and electrons, and creation/annihilation of charged particles. Two systems of differential equations are derived. The first system describes the whole magnetized plasma-wall transition region, which consists of the pre-sheath, the magnetized pre-sheath (Chodura layer), and the sheath, which is not neutral, but contains a positive space charge. The second system of equations describes only the neutral part of the plasma-wall transition region—this means only the pre-sheath and the Chodura layer, but not also the sheath. Both systems are solved numerically. The first system of equations has two singularities. The first occurs when ion velocity in the direction perpendicularly to the wall drops below the ion thermal velocity. The second occurs when the electron velocity in the direction perpendicularly to the wall exceeds the electron thermal velocity. The second system of differential equations only has one singularity, which has also been derived analytically. For finite electron to ion mass ratio, the integration of the second system always breaks down before the Bohm criterion is fulfilled. Some properties of the first system of equations are examined. It is shown that the increased collision frequency demagnetizes the plasma. On the other hand, if the magnetic field is so strong that the ion Larmor radius and the Debye length are comparable, the electron velocity in the direction perpendicularly to the wall reaches the electron thermal velocity before the ion velocity in the direction perpendicularly to the wall reaches the ion sound velocity. In this case, the integration of the model equations breaks down before the Bohm criterion is fulfilled and the sheath is formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1970

Dispersion discontinuities of strong collisionless shocks

Linear fluid equations are used to estimate wave train properties of strong collisionless shocks. Fast shocks exhibit several dispersion changes with increasing Mach number. For perpendicular propagation into a finite-beta plasma, an ion cyclotron radius trailing wave train exists only for (M sub F)2 is smaller than 2. Oblique fast shocks have a leading ion inertia wave train if M sub A is smaller than root of M(+)/M(-) cos theta/2 and a trailing electron inertia train if M sub A is greater than root of M(+)/M(-) cos theta/2. If the downstream sound speed exceeds the flow speed, linear wave theory predicts a trailing ion acoustic structure which probably resides within the magnetic shock. For a turbulent shock model in which an effective electron-ion collision frequency exceeds the lower hybrid frequency, ions decouple from the magnetic field; the shock wave train now trails with electron inertia and electron gyroradius lengths. Comparisons of this turbulent model and observations on the earth's bow shock are made.

Coroniti, F. V.↗
NASA NTRS1982

Collisionless electrostatic interchange instabilities

The linear Vlasov dispersion equation for electrostatic plasma instabilities driven by gravity and weak density gradients perpendicular to a uniform magnetic field is derived and solved numerically. Two interchange instabilities emerge: the well-known fluid mode at long wavelengths and a kinetic model at wavelengths short compared with the ion gyroradius. The properties of both instabilities are studied, as well as the effects of gravity on the universal and lower-hybrid density drift instabilities. The results show that the kinetic interchange generally has a larger growth rate than the fluid interchange instability, indicating that, whenever the latter is present in a collisionless plasma, the former may also be found.

Gray, S. P.↗
NASA NTRS1989

The electron beam instability revisited: Growth above and below f sub p

The growth of electrostatic waves near the plasma frequency (fp) due to an unstable electron beam is investigated by solving the unmagnetized electrostatic dispersion equation numerically. The numerical solutions are compared with analytic theories for reactive (or fluid-type) and kinetic versions of the beam instability, and for the O'Neil/Malmberg connection of the beam and Langmuir modes. Conditions for growth significantly above or below fp are given. Three general results are found: (1) The unstable waves do not grow on a mode with Langmuir dispersion except in the limit of a very dilute beam with growth on O'Neil/Malmberg's connected mode; (2) The properties of the unstable mode depend strongly on beam parameters such as beam density, speed and temperature; and (3) The frequency of maximum growth frequently lies significantly above or below fp, and differs significantly from that predicted by the Langmuir dispersion relation. Results imply important consequences for theories of strong turbulence and nonlinear wave-wave processes, and observational identification of fp from observed wave frequencies.

Cairns, I. H.↗
DOE OSTI2019

Modification of the Alfvén wave spectrum by pellet injection

Alfvén eigenmodes driven by energetic particles are routinely observed in tokamak plasmas. These modes consist of poloidal harmonics of shear Alfvén waves coupled by inhomogeneity in the magnetic field. In this work, further coupling is introduced by 3D inhomogeneities in the ion density during the assimilation of injected pellets. This additional coupling modifies the Alfvén continuum and discrete eigenmode spectrum. We find that the frequencies of Alfvén eigenmodes drop dramatically when a pellet is injected in JET. From these observations, information about the changes in the ion density caused by a pellet can be inferred. To use Alfvén eigenmodes for MHD spectroscopy of pellet injected plasmas, the 3D MHD codes Stellgap and AE3D were generalised to incorporate 3D density profiles. Additionally, a model for the expansion of the ionised pellet plasmoid along a magnetic field line was derived from the fluid equations. Thereby, the time evolution of the Alfvén eigenfrequency is reproduced. By comparing the numerical frequency drop of a toroidal Alfvén eigenmode (TAE) to experimental observations, the initial ion density of a cigar-shaped ablation region of length 4 cm is estimated to be n * =6.8×10 22 m -3 at the TAE location (r/a≈0.75). The frequency sweeping of an Alfvén eigenmode ends when the ion density homogenises poloidally. Modelling suggests that the time for poloidal homogenisation of the ion density at the TAE position is τ h =18±4 ms for inboard pellet injection, and τ h =26±2 ms for outboard pellet injection. By reproducing the frequency evolution of the elliptical Alfvén eigenmode (EAE), the initial ion density at the EAE location (r/a≈0.9) can be estimated to be n * =4.8×10 22 m -3 . Poloidal homogenisation of the ion density takes 2.7 times longer at the EAE location than at the TAE location for both inboard and outboard pellet injection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A multispecies, multifluid model for laser–induced counterstreaming plasma simulations

The interpenetration of counterstreaming plasmas is a critical phenomenon in several application areas, such as astrophysical flows, design of controlled fusion devices, and laser–induced plasma experiments. Multispecies “single-fluid” codes are unable to model this phenomenon due to the single velocity representation for all the species/fluids. Kinetic codes, though capable of modeling interpenetration, are computationally prohibitive for at-scale simulations. Here, we introduce a multifluid model that solves the fluid equations for each ion fluid or stream. This allows distinct flows that interact with each other through electrostatic and collisional forces. We observe and describe our code, EUCLID, that uses a conservative finite-difference formulation to discretize the governing equations in space. The 5th-order Monotonicity-Preserving WENO scheme is used for the upwind approximation of the hyperbolic flux, and the explicit 4th-order Runge–Kutta scheme is used for time integration. The code is verified for several benchmark cases and manufactured solutions. We simulate one- and two-dimensional interactions of counterstreaming plasmas in vacuum as well as in the presence of gas fill, where the setups are representative of laser-induced plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Physics of the intermediate layer between a plasma and a collisionless sheath and mathematical meaning of the Bohm criterion

A transformation of the ion momentum equation simplifies a mathematical description of the transition layer between a quasi-neutral plasma and a collisionless sheath and clearly reveals the physics involved. Balance of forces acting on the ion fluid is delicate in the vicinity of the sonic point and weak effects come into play. For this reason, the passage of the ion fluid through the sonic point, which occurs in the transition layer, is governed not only by inertia and electrostatic force but also by space charge and ion-atom collisions and/or ionization. Occurrence of different scenarios of asymptotic matching in the plasma-sheath transition is analyzed by means of simple mathematical examples, asymptotic estimates, and numerical calculations. In the case of a collisionless sheath, the ion speed distribution plotted on the logarithmic scale reveals a plateau in the intermediate region between the sheath and the presheath. The value corresponding to this plateau has the meaning of speed with which ions leave the presheath and enter the sheath; the Bohm speed. The plateau is pronounced reasonably well provided that the ratio of the Debye length to the ion mean free path is of the order of 10{sup -3} or smaller. There is no such plateau if the sheath is collisional and hence no sense in talking of a speed with which ions enter the sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Fluid equations in the presence of electron cyclotron current drive

Two-fluid equations, which include the physics imparted by an externally applied radiofrequency source near electron cyclotron resonance, are derived in their extended magnetohydrodynamic forms using the formalism of Hegna and Callen [Phys. Plasmas 16, 112501 (2009)]. The equations are compatible with the closed fluid/drift-kinetic model developed by Ramos [Phys. Plasmas 17, 082502 (2010); 18, 102506 (2011)] for fusion-relevant regimes with low collisionality and slow dynamics, and they facilitate the development of advanced computational models for electron cyclotron current drive-induced suppression of neoclassical tearing modes.

42 ENGINEERING↗