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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Experimental study of the stability and dynamics of a two-dimensional ideal vortex under external strain

The dynamics of two-dimensional (2-D) ideal fluid vortices is studied experimentally in the presence of an irrotational strain flow. Laboratory experiments are conducted using strongly magnetized pure electron plasmas, a technique which is made possible by the isomorphism between the drift–Poisson equations describing plasma dynamics transverse to the field and the 2-D Euler equations describing an ideal fluid. The electron plasma system provides an excellent opportunity to study the dynamics of a 2-D Euler fluid due to weak dissipation and weak 3-D effects, simple diagnosis and precise control. The plasma confinement apparatus used here was designed specifically to study vortex dynamics under the influence of external flow by applying boundary conditions in two dimensions. Additionally, vortex-in-cell simulations are carried out to complement the experimental results and to extend the parameter range of the studies. It is shown that the global dynamics of a quasi-flat vorticity profile is in good quantitative agreement with the theory of a piecewise-constant elliptical patch of vorticity, including the equilibria, dynamical orbits and stability properties. Deviations from the elliptical patch theory are observed for non-flat vorticity profiles; they include inviscid damping of the orbits and modified stability limits. The dependence of these phenomena on the flatness of the initial profile is discussed. The relationship of these results to other theoretical, numerical and experimental studies is also discussed.

Mechanics↗
DOE OSTI2016

Solitary fast magnetosonic waves propagating obliquely to the magnetic field in cold collisionless plasma

Solutions describing solitary fast magnetosonic (FMS) waves (FMS solitons) in cold magnetized plasma are obtained by numerically solving two-fluid hydrodynamic equations. The parameter domain within which steady-state solitary waves can propagate is determined. It is established that the Mach number for rarefaction FMS solitons is always less than unity. The restriction on the propagation velocity leads to the limitation on the amplitudes of the magnetic field components of rarefaction solitons. It is shown that, as the soliton propagates in plasma, the transverse component of its magnetic field rotates and makes a complete turn around the axis along which the soliton propagates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Experimental evidence of intermittent chaos in a glow discharge plasma without external forcing and its numerical modelling

Intermittent chaos was observed in a glow discharge plasma as the system evolved from regular type of relaxation oscillations (of larger amplitude) to an irregular type of oscillations (of smaller amplitude) as the discharge voltage was increased. Floating potential fluctuations were analyzed by different statistical and spectral methods. Features like a gradual change in the normal variance of the interpeak time intervals, a dip in the skewness, and a hump in the kurtosis with variation in the control parameter have been seen, which are strongly indicative of intermittent behavior in the system. Detailed analysis also suggests that the intrinsic noise level in the experiment increases with the increasing discharge voltage. An attempt has been made to model the experimental observations by a second order nonlinear ordinary differential equation derived from the fluid equations for an unmagnetized plasma. Though the experiment had no external forcing, it was conjectured that the intrinsic noise in the experiment could be playing a vital role in the dynamics of the system. Hence, a constant bias and noise as forcing terms were included in the model. Results from the theoretical model are in close qualitative agreement with the experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Traveling waves in Hall-magnetohydrodynamics and the ion-acoustic shock structure

Hall-magnetohydrodynamics (HMHD) is a mixed hyperbolic-parabolic partial differential equation that describes the dynamics of an ideal two fluid plasma with massless electrons. We study the only shock wave family that exists in this system (the other discontinuities being contact discontinuities and not shocks). We study planar traveling wave solutions and we find solutions with discontinuities in the hydrodynamic variables, which arise due to the presence of real characteristics in Hall-MHD. We introduce a small viscosity into the equations and use the method of matched asymptotic expansions to show that solutions with a discontinuity satisfying the Rankine-Hugoniot conditions and also an entropy condition have continuous shock structures. The lowest order inner equations reduce to the compressible Navier-Stokes equations, plus an equation which implies the constancy of the magnetic field inside the shock structure. We are able to show that the current is discontinuous across the shock, even as the magnetic field is continuous, and that the lowest order outer equations, which are the equations for traveling waves in inviscid Hall-MHD, are exactly integrable. We show that the inner and outer solutions match, which allows us to construct a family of uniformly valid continuous composite solutions that become discontinuous when the diffusivity vanishes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Ballooning theory of the second kind-two dimensional tokamak modes

The 2-D ballooning transform, devised to study local high toroidal number (n) fluctuations in axisymmetric toroidal system (like tokamaks), yields a well-defined partial differential equation for the linear eigenmodes. In this paper, such a ballooning equation of the second kind is set up for ion temperature gradient driven modes pertinent to a 2-D non-dissipative fluid plasma; the resulting partial differential equation is numerically solved, to calculate the global eigenvalues, and the 2-D mode structure is presented graphically along with analytical companions. The radial localization of the mode results from translational symmetry breaking for growing modes and is a vivid manifestation of spontaneous symmetry breaking in tokamak physics. The eigenmode, poloidally ballooned at {theta}={+-}{pi}/2, is radially shifted from associated rational surface. The global eigenvalue is found to be very close to the value obtained in 1-D parameterized ({lambda}={+-}{pi}/2) case. The 2-D eigenmode theory is applied to estimate the toroidal seed Reynolds stress [Y. Z. Zhang, Nucl. Fusion Plasma Phys. 30, 193 (2010)]. The solution obtained from the relatively simplified ballooning theory is compared to the solution of the basic equation in original coordinate system (evaluated via FFTs); the agreement is rather good.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Nonlinear electron acoustic waves in presence of shear magnetic field

Nonlinear electron acoustic waves are studied in a quasineutral plasma in the presence of a variable magnetic field. The fluid model is used to describe the dynamics of two temperature electron species in a stationary positively charged ion background. Linear analysis of the governing equations manifests dispersion relation of electron magneto sonic wave. Whereas, nonlinear wave dynamics is being investigated by introducing Lagrangian variable method in long wavelength limit. It is shown from finite amplitude analysis that the nonlinear wave characteristics are well depicted by KdV equation. The wave dispersion arising in quasineutral plasma is induced by transverse magnetic field component. The results are discussed in the context of plasma of Earth's magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Fluid description of multi-component solar partially ionized plasma

We derive self-consistent formalism for the description of multi-component partially ionized solar plasma, by means of the coupled equations for the charged and neutral components for an arbitrary number of chemical species, and the radiation field. All approximations and assumptions are carefully considered. Generalized Ohm's law is derived for the single-fluid and two-fluid formalism. Our approach is analytical with some order-of-magnitude support calculations. After general equations are developed, we particularize to some frequently considered cases as for the interaction of matter and radiation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Solar wind for a magnetized plasma with tensor plasma pressure.

A solar wind model for a magnetized solar wind is presented using one-fluid hydromagnetic equations with generalized polytrope equations of state for the two tensor components of the plasma pressure. Fluid and magnetic field variables are calculated at the Earth using certain boundary conditions at the Sun. The azimuthal velocity agrees with observed values and explains the considerable loss of angular momentum from the Sun by the solar wind. The results are good for most variables but suggest that a two-fluid model with electrons at higher temperatures and smaller temperature anisotropy ratios than the ions would give improved agreement for some quantities.

Tan, M.↗
DOE OSTI2013

UBIQUITOUS NON-THERMALS IN ASTROPHYSICAL PLASMAS: RESTATING THE DIFFICULTY OF MAINTAINING MAXWELLIANS

This paper outlines the rather narrow conditions on a radiatively decoupled plasma where a Maxwell-Boltzmann (MB) distribution can be assumed with confidence. The complementary non-thermal distribution with non-perturbative kurtosis is argued to have a much broader purview than has previously been accepted. These conditions are expressed in terms of the electron Knudsen number, K{sub e} , the ratio of the electron mean free path to the scale length of electron pressure. Rather generally, f(v < v{sub 2}(K{sub e} )) will be Gaussian, so that MB atomic or wave particle effects controlled by speeds v < v{sub 2} {identical_to} w(15/8K{sub e} ){sup 1/4} will remain defensible, where w is the most probable speed. The sufficient condition for Spitzer-Braginskii plasma fluid closure at the energy equation requires globally K{sub e} (s) {<=} 0.01; this global condition pertains to the maximum value of K{sub e} along the arc length s of the magnetic field (to its extremities) provided that contiguous plasma remains uncoupled from the radiation field. The non-thermal regime K{sub e} > 0.01 is common in all main-sequence stellar atmospheres above approximately 0.05 stellar radii from the surface. The entire solar corona and wind are included in this regime where non-thermal distributions with kurtosis are shown to be ubiquitous, heat flux is not well modeled by Spitzer-Braginskii closure, and fluid modeling is qualitative at best.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Closure models for the feedback of energetic particles on plasma turbulence

Energetic particles interact with the plasma surrounding them, resonating with certain types of plasma waves to stabilize them while destabilizing others, and changing the character of the background turbulence in ways that have not been fully quantified or understood. Interaction with the turbulent background plasma is key to the acceleration of many types of energetic particles including high-energy cosmic rays, solar energetic particles, and pick-up ions. The acceleration of particles is a process that would ideally be described by a kinetic model, a type of model that follows a probability distribution function (PDF) for all particles in 7-dimensional (x, y, z, v x , v y , v z , t) space. Because of the high dimensionality of a kinetic model, simulations that solve kinetic equations use the largest computational resources currently available, and are yet unable to simulate a realistic number of particles, reach the large scales necessary for astrophysical problems, and use high-precision numerical methods. Two available alternatives to kinetic plasma models have been explored for this problem, with limited success. One is a multi-fluid model produced by a cumulant discarding closure, which evolves coupled equations for the velocity, magnetic field, and internal energy for both the background plasma and the fluid of energetic particles. However, simulations that solve multi-fluid magnetohydrodynamic (MHD) equations are able to include the interaction with energetic particles only in crude ways, typically as an add-on pressure term. The second alternative is to use a hybrid method to couple a fluid description of the background plasma to a kinetic model or a Fokker–Planck model for the energetic particles. These methods are hampered by the physical modeling of the coupling. In this work, we develop a new model, which follows the PDF for all particles; this can be viewed as a step toward physical realism above a multi-fluid MHD model, while also being more computationally efficient than a kinetic model. The equations we develop model both the background plasma and the energetic particles self-consistently. Over the last decade, similar PDF methods have been developed to a high level of sophistication to model reactive flows and turbulent combustion for engineering applications. For treatment of the feedback of the energetic particles on a background plasma, a PDF closure approach should evaluate the mean characteristics, including the density, with better statistical quality than will particle-sampling procedures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Comparisons of time explicit hybrid kinetic-fluid code Architect for Plasma Wakefield Acceleration with a full PIC code

Architect, a time explicit hybrid code designed to perform quick simulations for electron driven plasma wakefield acceleration, is described. In order to obtain beam quality acceptable for applications, control of the beam-plasma-dynamics is necessary. Particle in Cell (PIC) codes represent the state-of-the-art technique to investigate the underlying physics and possible experimental scenarios; however PIC codes demand the necessity of heavy computational resources. Architect code substantially reduces the need for computational resources by using a hybrid approach: relativistic electron bunches are treated kinetically as in a PIC code and the background plasma as a fluid. Cylindrical symmetry is assumed for the solution of the electromagnetic fields and fluid equations. In this paper both the underlying algorithms as well as a comparison with a fully three dimensional particle in cell code are reported. The comparison highlights the good agreement between the two models up to the weakly non-linear regimes. In highly non-linear regimes the two models only disagree in a localized region, where the plasma electrons expelled by the bunch close up at the end of the first plasma oscillation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

On the space-charge formation in a collisional magnetized electronegative plasma

The plasma sheath formation in the vicinity of a surrounding wall of magnetized plasma is studied in the presence of the electronegative ions and the positive ion-neutral background collisions. Fluid equations are used to treat the plasma particles species. By using the Sagdeev potential, the influence of the collisions and the magnetic field on the Bohm criterion are investigated. The space-charge profiles are obtained in the presence of a magnetic field in different collision frequencies as well as electronegative ions concentration. It is shown that the collision and the magnetic field raise a space-charge peak, while the presence of the electronegative ions results in damping the peaks. Moreover, it is observed that in the case of high magnetic field, some fluctuations emerge in the space-charge profiles. The influences of the magnetic field and electronegative ion concentration as well as negative ion temperature on the positive ion kinetic energy reaching the plasma surrounding wall and positive ion velocity perpendicular to the sheath axis are investigated. Finally, the net current through the sheath region is obtained for different collisionality and magnetic field values in both electropositive and electronegative plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Plasma ion stratification by weak planar shocks

We derive fluid equations for describing steady-state planar shocks of a moderate strength (0 < M - 1 ≲ 1 with M the shock Mach number) propagating through an unmagnetized quasineutral collisional plasma comprising two separate ion species. In addition to the standard fluid shock quantities, such as the total mass density, mass-flow velocity, and electron and average ion temperatures, the equations describe shock stratification in terms of variations in the relative concentrations and temperatures of the two ion species along the shock propagation direction. We have solved these equations analytically for weak shocks (0 < M - 1 << 1), with the results depending on M, ratios of the ion masses and charges, and the upstream mass fraction of one of the ion species. These analytical results are instrumental for gaining understanding in the behavior of weak shocks, and they have been used to verify kinetic simulations of shocks in multi-ion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Eigenvalue Solver for Fluid and Kinetic Plasma Models in Arbitrary Magnetic Topology

ArbiTER (Arbitrary Topology Equation Reader) is a new code for solving linear eigenvalue problems arising from a broad range of physics and geometry models. The primary application area envisioned is boundary plasma physics in magnetic confinement devices; however ArbiTER should be applicable to other science and engineering fields as well. The code permits a variable numbers of dimensions, making possible application to both fluid and kinetic models. Additionally, the use of specialized equation and topology parsers permits a high degree of flexibility in specifying the physics and geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

A fast non-Fourier method for Landau-fluid operators

An efficient and versatile non-Fourier method for the computation of Landau-fluid (LF) closure operators [Hammett and Perkins, Phys. Rev. Lett. 64, 3019 (1990)] is presented, based on an approximation by a sum of modified-Helmholtz-equation solves (SMHS) in configuration space. This method can yield fast-Fourier-like scaling of the computational time requirements and also provides a very compact data representation of these operators, even for plasmas with large spatial nonuniformity. As a result, the method can give significant savings compared with direct application of “delocalization kernels” [e.g., Schurtz et al., Phys. Plasmas 7, 4238 (2000)], both in terms of computational cost and memory requirements. The method is of interest for the implementation of Landau-fluid models in situations where the spatial nonuniformity, particular geometry, or boundary conditions render a Fourier implementation difficult or impossible. Systematic procedures have been developed to optimize the resulting operators for accuracy and computational cost. The four-moment Landau-fluid model of Hammett and Perkins has been implemented in the BOUT++ code using the SMHS method for LF closure. Excellent agreement has been obtained for the one-dimensional plasma density response function between driven initial-value calculations using this BOUT++ implementation and matrix eigenvalue calculations using both Fourier and SMHS non-Fourier implementations of the LF closures. The SMHS method also forms the basis for the implementation, which has been carried out in the BOUT++ code, of the parallel and toroidal drift-resonance LF closures. The method is a key enabling tool for the extension of gyro-Landau-fluid models [e.g., Beer and Hammett, Phys. Plasmas 3, 4046 (1996)] to codes that treat regions with strong profile variation, such as the tokamak edge and scrapeoff-layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

IMEX and exact sequence discretization of the multi-fluid plasma model

Multi-fluid plasma models, where an electron fluid is modeled in addition to multiple ion and neutral species as well as the full set of Maxwell's equations, are useful for representing physics beyond the scope of classic MHD. This advantage presents challenges in appropriately dealing with electron dynamics and electromagnetic behavior characterized by the plasma and cyclotron frequencies and the speed of light. For physical systems, such as those near the MHD asymptotic regime, this requirement drastically increases runtimes for explicit time integration even though resolving fast dynamics may not be critical for accuracy. Implicit time integration methods, with efficient solvers, can help to step over fast time-scales that constrain stability, but do not strongly influence accuracy. As an extension, Implicit-explicit (IMEX) schemes provide an additional mechanism to choose which dynamics are evolved using an expensive implicit solve or resolved using a fast explicit solve. In this study, in addition to IMEX methods we also consider a physics compatible exact sequence spatial discretization. Here, this combines nodal bases (H-Grad) for fluid dynamics with a set of vector bases (H-Curl and H-Div) for Maxwell's equations. This discretization allows for multi-fluid plasma modeling without violating Gauss' laws for the electric and magnetic fields. This initial study presents a discussion of the major elements of this formulation and focuses on demonstrating accuracy in the linear wave regime and in the MHD limit for both a visco-resistive and a dispersive ideal MHD problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗