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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Kinetic stability of Chapman–Enskog plasmas

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$ , and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$ , the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma ( viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$ , the electron inertial scale $d_e$ and the plasma $\beta$ . This work is of general consequence in emphasising the fact that high- $\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$ . In this situation, the determination of transport coefficients via the standard CE approach is not valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Theory of spatially non-symmetric kinetic equilibria for collisionless plasmas

The problem posed by the possible existence/non-existence of spatially non-symmetric kinetic equilibria has remained unsolved in plasma theory. For collisionless magnetized plasmas, this involves the construction of stationary solutions of the Vlasov-Maxwell equations. In this paper, the issue is addressed for non-relativistic plasmas both in astrophysical and laboratory contexts. The treatment is based on a Lagrangian variational description of single-particle dynamics. Starting point is a non-perturbative formulation of gyrokinetic theory, which allows one to construct 'a posteriori' with prescribed order of accuracy an asymptotic representation for the magnetic moment. In terms of the relevant particle adiabatic invariants generalized bi-Maxwellian equilibria are proved to exist. These are shown to recover, under suitable assumptions, a Chapman-Enskog form which permits an analytical treatment of the corresponding fluid moments. In particular, the constrained posed by the Poisson and the Ampere equations are analyzed, both for quasi-neutral and non-neutral plasmas. The conditions of existence of the corresponding non-symmetric kinetic equilibria are investigated. As a notable feature, both astrophysical and laboratory plasmas are shown to exhibit, under suitable conditions, a kinetic dynamo, whereby the equilibrium magnetic field can be self-generated by the equilibrium plasma currents.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Effect of polarization force on the Mach cones in a complex plasma

We report the modifications of compressional Mach cone propagation characteristics due to the polarization force acting on micron size dust particles embedded in a non-uniform plasma. We solve the hydrodynamic fluid equations for highly charged dust particles to investigate the Mach cone by incorporating the polarization force in the momentum equation and observe the structural change on lateral wakes at different polarization force for a given Mach number and Epstein drag force. We also notice that the maximum amplitude of normalized dust density perturbation decreases with the increase of polarization interaction when the other parameters remain constant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Foundations of magnetohydrodynamics

In this tutorial, a derivation of magnetohydrodynamics (MHD) valid beyond the usual ideal gas approximation is presented. Non-equilibrium thermodynamics is used to obtain conservation equations and linear constitutive relations. When coupled with Maxwell's equations, this provides closed fluid equations in terms of material properties of the plasma, described by the equation of state and transport coefficients. These properties are connected to microscopic dynamics using the Irving–Kirkwood procedure and Green–Kubo relations. Symmetry arguments and the Onsager–Casimir relations allow one to vastly simplify the number of independent coefficients. Importantly, expressions for current density, heat flux, and stress (conventionally Ohm's law, Fourier's law, and Newton's law) take different forms in systems with a non-ideal equation of state. The traditional form of the MHD equations, which is usually obtained from a Chapman–Enskog solution of the Boltzmann equation, corresponds to the ideal gas limit of the general equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids

This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.

graphics processing units↗
DOE OSTI2013

On the fully nonlinear acoustic waves in a plasma with positrons beam impact and superthermal electrons

Arbitrary amplitude ion-acoustic waves in an unmagnetized plasma consisting of cold positive ions, superthermal electrons, and positrons beam are reported. The basic set of fluid equations is reduced to an energy-balance like equation. The latter is numerically analyzed to examine the existence regions for solitary and shock waves. It is found that only solitary waves can propagate, however, the model cannot support shocks. The effects of superthermality and beam parameters (via, positrons concentration and streaming velocity) on the existence region, as well as solitary wave profile have been discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A STUDY OF ALFVÉN WAVE PROPAGATION AND HEATING THE CHROMOSPHERE

Alfvén wave propagation, reflection, and heating of the chromosphere are studied for a one-dimensional solar atmosphere by self-consistently solving plasma, neutral fluid, and Maxwell's equations with incorporation of the Hall effect and strong electron-neutral, electron-ion, and ion-neutral collisions. We have developed a numerical model based on an implicit backward difference formula of second-order accuracy both in time and space to solve stiff governing equations resulting from strong inter-species collisions. A non-reflecting boundary condition is applied to the top boundary so that the wave reflection within the simulation domain can be unambiguously determined. It is shown that due to the density gradient the Alfvén waves are partially reflected throughout the chromosphere and more strongly at higher altitudes with the strongest reflection at the transition region. The waves are damped in the lower chromosphere dominantly through Joule dissipation, producing heating strong enough to balance the radiative loss for the quiet chromosphere without invoking anomalous processes or turbulences. The heating rates are larger for weaker background magnetic fields below ∼500 km with higher-frequency waves subject to heavier damping. There is an upper cutoff frequency, depending on the background magnetic field, above which the waves are completely damped. At the frequencies below which the waves are not strongly damped, the interaction of reflected waves with the upward propagating waves produces power at their double frequencies, which leads to more damping. The wave energy flux transmitted to the corona is one order of magnitude smaller than that of the driving source.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Improved analytic modeling of the linear Rayleigh–Taylor instability with plasma transport

The Rayleigh–Taylor (R–T) instability with transport effects has long been studied analytically in the neutral fluid regime. How this picture may differ in plasmas has been given less consideration, however. Leveraging previous numerical simulations with plasma viscosity and diffusion of R–T at a binary plasma interface, we show how plasma transport and kinetic effects alter the linear R–T dispersion relation. Our insight informs a much improved analytical dispersion relation, which better matches the simulation data than previously established models. Additionally, we demonstrate the importance of employing asymptotically correct viscosity coefficients in fluid codes for accurately capturing plasma R–T growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Inhibition of stimulated Raman scattering due to the excitation of stimulated Brillouin scattering

The nonlinear coupling between stimulated Raman scattering (SRS) and stimulated Brillouin scattering (SBS) of intense laser in underdense plasma is studied theoretically and numerically. Based upon the fluid model, their coupling equations are derived, and a threshold condition of plasma density perturbations due to SBS for the inhibition of SRS is given. Particle-in-cell simulations show that this condition can be achieved easily by SBS in the so-called fluid regime with kLλD<0.15, where kL is the Langmuir wave number and λD is the Debye length [Kline et al., Phys. Plasmas 13, 055906 (2006)]. SBS can reduce the saturation level of SRS and the temperature of electrons in both homogeneous and inhomogeneous plasma. Numerical simulations also show that this reduced SRS saturation is retained even if the fluid regime condition mentioned above is violated at a later time due to plasma heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Nodal soliton solutions for generalized quasilinear Schrödinger equations

This paper is concerned with constructing nodal radial solutions for generalized quasilinear Schrödinger equations in R{sup N} which arise from plasma physics, fluid mechanics, as well as high-power ultashort laser in matter. For any given integer k ⩾ 0, by using a change of variables and minimization argument, we obtain a sign-changing minimizer with k nodes of a minimization problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Modelling general-relativistic plasmas with collisionless moments and dissipative two-fluid magnetohydrodynamics

ABSTRACT Relativistic plasmas are central to the study of black hole accretion, jet physics, neutron star mergers, and compact object magnetospheres. Despite the need to accurately capture the dynamics of these plasmas and the implications for relativistic transients, their fluid modelling is typically done using a number of (overly) simplifying assumptions, which do not hold in general. This is especially true when the mean free path in the plasma is large compared to the system size, and kinetic effects start to become important. Going beyond common approaches used in the literature, we describe a fully relativistic covariant 14-moment based two-fluid system appropriate for the study of electron–ion or electron–positron plasmas. This generalized Israel-Stewart-like system of equations of motion is obtained directly from the relativistic Boltzmann–Vlasov equation. This new formulation can account for non-ideal effects, such as anisotropic pressures and heat fluxes, not present in previous formulations of two-fluid magnetohydrodynamics. We show that a relativistic two-fluid plasma can be recast as a single fluid coupled to electromagnetic fields with (potentially large) out-of-equilibrium corrections. We keep all electron degrees of freedom, which provide self-consistent evolution equations for electron temperature and momentum. The out-of-equilibrium corrections take the form of a collisional 14-moment closure previously described in the context of viscous single fluids. The equations outlined in this paper are able to capture the full two-fluid character of collisionless plasmas found in black hole accretion and flaring processes around compact objects, as well Braginskii-like two-fluid magnetohydrodynamics applicable to weakly collisional plasmas inside accretion discs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

INDUCED MAXIMUM MAGNETIC FIELD IN COSMIC OUTFLOW SYSTEM BY A RELATIVISTIC CURRENT FILAMENTATION INSTABILITY: EXACT ANALYTICAL MODEL

We present an analytical study of current filamentation instability (CFI) in a fully relativistic cold plasma system, including arbitrary currents. For our purposes, we employ the cold fluid equations, together with Maxwell's equations as well as the plasma shell concept and boost frame method, to obtain an exact solution of the instability growth rate. A simple relation is found for the maximum growth rate of the CFI (for any arbitrary current system), which remarkably is used to calculate the large magnitude of an induced magnetic field in astrophysical environments such as active galactic nuclei (AGNs), microquasars, supernova remnants (SNRs), and stellar winds. We find that the magnetic field is amplified in the SNR up to the level required to justify the recent discovery of the year-scale variability in the X-ray emission of SNRs. Also, the maximum magnetic field of two and three orders higher (or two orders lower) than that of the SNR has been derived for microquasars and AGNs (or stellar winds), respectively. Moreover, making use of the exact analytical solution of the CFI, it is shown that the maximum magnetic field up to around 10{sup 8} G can be detected from a classical cold counterstreaming system after a saturation time.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Ionospheric plasma flow about a system of electrically biased flat plates

The steady state interaction of two electrically biased parallel plates immersed in a flowing plasma characteristic of low earth orbit is studied numerically. Fluid equations are developed to describe the motion of the cold positively charged plasma ions, and are solved using finite-differences in two dimensions on a Cartesian grid. The behavior of the plasma electrons is assumed to be described by the Maxwell-Boltzmann distribution. Results are compared to an analytical and a particle simulation technique for a simplified flow geometry consisting of a single semi-infinite negatively biased plate. Comparison of the extent of the electrical disturbance into the flowing plasma and the magnitude of the current collected by the plate is very good. The interaction of two equally biased parallel plates is studied as a function of applied potential. The separation distance at which the current collected by either plate decreases by five and twenty percent is determined as a function of applied potential. The percent decreases were based on a non-interacting case. The decrease in overall current is caused by a decrease in ionic density in the region between the plates. As the separation between the plates decreases, the plates collect the ions at a faster rate than they are supplied to the middle region by the oncoming plasma flow. The docking of spacecraft in orbit is simulated by moving two plates of unequal potential toward one another in a quasi-static manner. One plate is held at a large negative potential while the other floats electrically in the resulting potential field. It is found that the floating plate does not charge continuously negative as it approaches the other more negatively biased plate. Instead, it charges more and then less negative as ionic current decreases and then increases respectively upon approach. When the two plates come into contact, it is expected that the electrically floating plate will charge rapidly negative to a potential near that of the other plate.

Herr, Joel L.↗
DOE OSTI2022

Inviscid damping of an elliptical vortex subject to an external strain flow

Inviscid spatial Landau damping is studied experimentally for the case of oscillatory motion of a two-dimensional vortex about its elliptical equilibrium in the presence of an applied strain flow. Here, the experiments are performed using electron plasmas in a Penning–Malmberg trap. They exploit the isomorphism between the two-dimensional Euler equations for an ideal fluid and the drift-Poisson equations for the plasma, where plasma density is the analog of vorticity. Perturbed elliptical vortex states are created using E x B strain flows, which are generated by applying voltages to electrodes surrounding the plasma. Measurements of spatial Landau damping (also called critical-layer damping) are in agreement with previous studies in the absence of an applied strain, where the damping is due to a resonance between the local fluid motion and the vortex oscillations. Interestingly, the damping rate does not change significantly over a wide range of applied strain rates. This can be accurately predicted from the initial vorticity profile, even though the resonant frequency is reduced substantially due to the applied strain. For higher amplitude perturbations, nonlinear trapping oscillations also exhibit behavior similar to the strain-free case. In principle, higher-order effects of the applied strain, such as separatrix crossing of peripheral vorticity and interactions with harmonics of the fundamental resonance, are expected to change the damping rate. However, this occurs only for conditions that are not realized in the experiments described here. Vortex-in-cell simulations are used to investigate the possible roles of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Repeated-cascade theory of strong turbulence in a magnetized plasma

A two-dimensional Navier-Stokes equation of vorticity in fluid turbulence is used to model drift turbulence in a plasma with a strong constant magnetic field and a constant mean density gradient. The nonlinear eddy diffusivity is described by a time-integrated Lagrangian correlation of velocities, and the repeated-cascade method is employed to choose the rank accounting for nearest-neighbor interactions, to calculate the Lagrangian correlation, and to close the correlation hierarchy. As a result, the diffusivity becomes dependent on the plasma's induced diffusion and is represented by a memory chain that is cut off by similarity and inertial randomization. Spectral laws relating the kinetic-energy spectrum to the -5, -5/2, -3, and -11 powers of wavenumber are derived for the velocity subranges of production, approach to inertia, inertia, and dissipation, respectively. It is found that the diffusivity is proportional to some inverse power of the magnetic field, that power being 1, 2/3, 5/6, and 2, respectively, for the four velocity subranges.

Tchen, C. M.↗