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

Pseudo-Spectral Methods for the Laplace-Beltrami Equation and the Hodge Decomposition on Surfaces of Genus One

The inversion of the Laplace-Beltrami operator and the computation of the Hodge decomposition of a tangential vector field on smooth surfaces arise as computational tasks in many areas of science, from computer graphics to machine learning to computational physics. In this work, we present a high-order accurate pseudo-spectral approach, applicable to closed surfaces of genus one in three-dimensional space, with a view toward applications in plasma physics and fluid dynamics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Radial energy transport by magnetospheric ULF waves: Effects of magnetic curvature and plasma pressure

The 'radial' transport of energy by internal ULF waves, stimulated by dayside magnetospheric boundary oscillations, is analyzed in the framework of one-fluid magnetohydrodynamics. (the term radial is used here to denote the direction orthogonal to geomagnetic flux surfaces.) The model for the inhomogeneous magnetospheric plasma and background magnetic field is axisymmetric and includes radial and parallel variations in the magnetic field, magnetic curvature, plasma density, and low but finite plasma pressure. The radial mode structure of the coupled fast and intermediate MHD waves is determined by numerical solution of the inhomogeneous wave equation; the parallel mode structure is characterized by a Wentzel-Kramer-Brillouin (WKB) approximation. Ionospheric dissipation is modeled by allowing the parallel wave number to be complex. For boudnary oscillations with frequencies in the range from 10 to 48 mHz, and using a dipole model for the background magnetic field, the combined effects of magnetic curvature and finite plasma pressure are shown to (1) enhance the amplitude of field line resonances by as much as a factor of 2 relative to values obtained in a cold plasma or box-model approximation for the dayside magnetosphere; (2) increase the energy flux delivered to a given resonance by a factor of 2-4; and (3) broaden the spectral width of the resonance by a factor of 2-3. The effects are attributed to the existence of an 'Alfven buoyancy oscillation,' which approaches the usual shear mode Alfven wave at resonance, but unlike the shear Alfven mode, it is dispersive at short perpendicular wavelengths. The form of dispersion is analogous to that of an internal atmospheric gravity wave, with the magnetic tension of the curved background field providing the restoring force and allowing radial propagation of the mode. For nominal dayside parameters, the propagation band of the Alfven buoyancy wave occurs between the location of its (field line) resonance and that of the fast mode cutoff that exists at larger radial distances.

Kouznetsov, Igor↗
DOE OSTI2012

Properties of a warm plasma collisional sheath in an oblique magnetic field

The properties of a warm plasma collisional sheath in an oblique magnetic field and the associated sheath criterion are investigated with a two-fluid model. In the fluid framework, a sheath criterion including effects of the magnetic field and collision is established theoretically for a wide range of ion temperature. With the sheath criterion as the plasma-sheath boundary condition, different plasma parameters including potential, electron and ion densities, and ion velocity are calculated for various ion temperatures and ion thermal motions. It is shown that the properties of the sheath depend not only on the plasma balance equations but also on the sheath boundary conditions. In addition, effects of the directions and magnitudes of the magnetic field on the plasma sheath are also discussed under different ion temperatures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Three-dimensional simulation of microwave-induced helium plasma under atmospheric pressure

A three-dimensional model is presented to investigate helium plasma generated by microwave under atmospheric pressure in this paper, which includes the physical processes of electromagnetic wave propagation, electron and heavy species transport, gas flow, and heat transfer. The model is based on the fluid approximation calculation and local thermodynamic equilibrium assumption. The simulation results demonstrate that the maxima of the electron density and gas temperature are 4.79 × 10{sup 17 }m{sup −3} and 1667 K, respectively, for the operating conditions with microwave power of 500 W, gas flow rate of 20 l/min, and initial gas temperature of 500 K. The electromagnetic field distribution in the plasma source is obtained by solving Helmholtz equation. Electric field strength of 2.97 × 10{sup 4 }V/m is obtained. There is a broad variation on microwave power, gas flow rate, and initial gas temperature to obtain deeper information about the changes of the electron density and gas temperature.

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↗
DOE OSTI2020

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Charged particle motion in spherically symmetric distributions of magnetic monopoles

The classical equations of motion of a charged particle in a spherically symmetric distribution of magnetic monopoles can be transformed into a system of linear equations, thereby providing a type of integrability. In the case of a single monopole, the solution was given long ago by Poincaré. In the case of a uniform distribution of monopoles, the solution can be expressed in terms of parabolic cylinder functions (essentially the eigenfunctions of an inverted harmonic oscillator). Further, this solution is relevant to recent studies of nonassociative star products, symplectic lifts of twisted Poisson structures, and fluids and plasmas of electric and magnetic charges.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Effects of electromagnetic wiggler and ion channel guiding on equilibrium orbits and waves propagation in a free electron laser

In this paper, an analysis of equilibrium orbits for electrons by a simultaneous solution of the equation of motion and the dispersion relation for electromagnetic wave wiggler in a free-electron laser (FEL) with ion-channel guiding has been presented. A fluid model has been used to investigate interactions among all possible waves. The dispersion relation has been derived for electrostatic and electromagnetic waves with all relativistic effects included. This dispersion relation has been solved numerically. For group I and II orbits, when the transverse velocity is small, only the FEL instability is found. In group I and II orbits with relatively large transverse velocity, new couplings between other modes are found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

On stability criteria for kinetic magnetohydrodynamics

The existence of a potential energy functional in the zero-Larmor-radius collisionless plasma theory of Kruskal & Oberman (Phys. Fluids, vol. 1, 1958 p. 275), Rosenbluth & Rostoker (Phys. Fluids, vol. 2, 1959, p. 23) allows us to derive easily sufficient conditions for linear stability. However, this kinetic magnetohydrodynamics (KMHD) theory does not have a self-adjointness property, making it difficult to derive necessary conditions. In particular, the standard methods to prove that an instability follows if some trial perturbation makes the incremental potential energy negative, which rely on the self-adjointness of the force operator or on the existence of a complete basis of normal modes, are not applicable to KMHD. This paper investigates KMHD linear stability criteria based on the time evolution of initial-value solutions, without recourse to the classic bounds or comparison theorems of Kruskal–Oberman and Rosenbluth–Rostoker for the KMHD potential energy. The adopted approach does not solve the kinetic equations by integration along characteristics and does not require that the particle orbits be periodic or nearly periodic. Most importantly, the investigation of a necessary condition for stability does not require the self-adjointness of the force operator or the existence of a complete basis of normal modes. Thus, it is thereby shown that stability in isothermal ideal-MHD is a sufficient condition for stability in KMHD and that, with a proviso on the long-time behaviour of oscillations about stable equilibria, stability in the double-adiabatic fluid theory, including the variation of the parallel fluid displacement, would be a necessary condition for stability in KMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulation of a strong steady-state plasma shock in a warm dense matter regime

The structure of collisional plasma shocks has been subject to an extensive, multi-decadal investigation—in the hydrodynamic, hybrid kinetic ion/electron fluid, and fully kinetic ion/electron limits. Despite this thoroughness, all of these studies apply exclusively to classical, weakly coupled plasmas. Here, we report the first results for a planar hydrodynamic simulation of a strong, steady-state shock in a subspace of the warm dense matter (WDM) regime. Specifically, we consider a plasma of fully degenerate electrons with moderate-to-strongly coupled ions. Since the WDM ion and electron transport coefficients and equation of state differ markedly from their non-degenerate, weak-coupling equivalents, we find that the structure of a WDM plasma shock notably deviates from the ideal plasma picture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Magnetohydrodynamic modeling of the solar wind in the outer heliosphere

We present initial results from a solar wind model that accounts for transport of turbulence and treats pickup protons as a separate fluid. The model is based on a numerical solution of the coupled set of mean-field Reynolds-averaged solar wind equations and small-scale turbulence transport equations in the region from 0.3-100 AU. The pickup protons are assumed to be comoving with the solar wind flow and described by separate mass and energy equations. The equations include the terms for energy transfer from pickup protons to solar wind protons and for the plasma heating by turbulent dissipation. The momentum equation contains a term that describes the loss of momentumby the solar wind flow due to the charge exchangewith the interstellar neutral hydrogen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerical study on a canonized Hamiltonian system representing reduced magnetohydrodynamics and its comparison with two-dimensional Euler system

Introducing a Clebsch-like parameterization, we have formulated a canonical Hamiltonian system on a symplectic leaf of reduced magnetohydrodynamics. An interesting structure of the equations is in that the Lorentz-force, which is a quadratic nonlinear term in the conventional formulation, appears as a linear term −ΔQ, just representing the current density (Q is a Clebsch variable, and Δ is the two-dimensional Laplacian); omitting this term reduces the system into the two-dimensional Euler vorticity equation of a neutral fluid. A heuristic estimate shows that current sheets grow exponentially (even in a fully nonlinear regime) together with the action variable P that is conjugate to Q. By numerical simulation, the predicted behavior of the canonical variables, yielding exponential growth of current sheets, has been demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Body Fluids Monitor

Method and apparatus are described for determining volumes of body fluids in a subject using bioelectrical response spectroscopy. The human body is represented using an electrical circuit. Intra-cellular water is represented by a resistor in series with a capacitor; extra-cellular water is represented by a resistor in series with two parallel inductors. The parallel inductors represent the resistance due to vascular fluids. An alternating, low amperage, multifrequency signal is applied to determine a subject's impedance and resistance. From these data, statistical regression is used to determine a 1% impedance where the subject's impedance changes by no more than 1% over a 25 kHz interval. Circuit component, of the human body circuit are determined based on the 1% impedance. Equations for calculating total body water, extra-cellular water, total blood volume, and plasma volume are developed based on the circuit components.

Siconolfi, Steven F.↗
DOE OSTI2017

Investigation of magnetic flux transport and shock formation in a staged Z-pinch

Target preheating is an integral component of magnetized inertial fusion in reducing convergence ratio. In the staged Z-pinch concept, it is achieved via one or more shocks. Previous work [Narkis et al., Phys. Plasmas 23, 122706 (2016)] found that shock formation in the target occurred earlier in higher-Z liners due to faster flux transport to the target/liner interface. However, a corresponding increase in magnitude of magnetic pressure was not observed, and target implosion velocity (and therefore shock strength) remained unchanged. To investigate other means of increasing the magnitude of transported flux, a Korteweg-de Vries-Burgers equation from the 1-D single-fluid, resistive magnetohydrodynamic equations is obtained. Solutions to the nondispersive (i.e., Burgers) equation depend on nondimensional coefficients, whose dependence on liner density, temperature, etc., suggests an increase in target implosion velocity, and therefore shock strength, can be obtained by tailoring the mass of a single-liner gas puff to a double-liner configuration. In the selected test cases of 1-D simulated implosions of krypton on deuterium, the peak Mach number increased from ∼5 to ∼8. While a notable increase was seen, Mach numbers exceeding 10 (implosion velocities exceeding ∼25 cm/μs) are necessary for adequate shock preheating.

Narkis, J. (ORCID:0000000154368933)↗
DOE OSTI2019

Kelvin-Helmholtz instability is the result of parity-time symmetry breaking

Parity-Time (PT)-symmetry is being actively investigated as a fundamental property of observables in quantum physics. We show that the governing equations of the classical two-fluid interaction and the incompressible fluid system are PT-symmetric, and the well-known Kelvin-Helmholtz instability is the result of spontaneous PT-symmetry breaking. It is expected that all classical conservative systems governed by Newton's law admit PT-symmetry, and the spontaneous breaking thereof is a generic mechanism for classical instabilities. Here discovering the PT-symmetry of systems in fluid dynamics and plasma physics and identifying the PT-symmetry breaking responsible for instabilities enable new techniques to classical physics and enrich the physics of PT-symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

The Role of Extracellular Fluid in Biokinetic Modeling

Here, the pharmacokinetic equations of Pierson et al. describing the behavior of bromide in rat provide a general approach to the modeling of extracellular fluid (ECF). The movement of material into ECF spaces is rapid and is completely characterized by tissue volumes and vascular flow rates to and from a tissue, the volumes of the tissue, and the ECF associated with the tissue. Early-time measurements are needed to characterize ECF. Measurements of DTPA disappearance from plasma by Wedeking et al. are discussed as an example of such measurements. In any biokinetic model, the fastest transfer rates are not determinable with the usual datasets, and if determined empirically, these rates will have very large and highly correlated uncertainties, so particular values of these rates, even though the model fits the available data, are not significant. A pharmacokinetic front-end provides values for these fast rates. An example of such a front-end for a 200–g rat is given.

61 RADIATION PROTECTION AND DOSIMETRY↗