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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

Instabilities in uranium plasma and the gas-core nuclear rocket engine

The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P approximates 0.00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.

Tidman, D. A.↗
DOE OSTI2023

Novel features of attractors and transseries in nonconformal Bjorken flows

In this work we investigate the impact of conformal symmetry breaking on hydrodynamization of a far-from-equilibrium fluid. We find a new kind of transseries solutions for the nonconformal hydrodynamic equations of a longitudinal boost invariant expanding plasma. The new transseries solutions unveil a rich physical structure which arises due to the interplay of different physical scales. In the perfect fluid case the nonconformal speed of sound slows down the cooling of the temperature due to the emergence of logarithmic corrections that depends on the mass of the particle. Further, these terms propagate into the perturbative and nonperturbative sectors of the transseries once viscous corrections are included. The logarithmic mass contributions increase the asymptotic value of the Knudsen number while decreasing the damping rate of the transient nonhydrodynamic modes and thus, yielding to an extremely slow hydrodynamization process where flow lines merge to their forward attractor at extremely late times. The early time free streaming expansion is modified and receives logarithmic mass corrections induced by the shear-bulk couplings. The global flow structure and numerical analyses carried out in our work demonstrate the existence of local attraction at early and late times for the shear viscous tensor and bulk viscous pressure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1971

Instabilities in uranium plasma.

The nonlinear evolution of unstable sound waves in a uranium plasma has been calculated using a multiple time-scale asymptotic expansion scheme. The fluid equations used include the fission power density, radiation diffusion, and the effects of the changing degree of ionization of the uranium atoms. The nonlinear growth of unstable waves is shown to be limited by mode coupling to shorter wavelength waves which are damped by radiation diffusion. This mechanism limits the wave pressure fluctuations to values of order delta P/P equal to about .00001 in the plasma of a typical gas-core nuclear rocket engine. The instability is thus not expected to present a control problem for this engine.

Tidman, D. A.↗
DOE OSTI2014

Theoretical and numerical modelling of chaotic electrostatic ion cyclotron (EIC) oscillations by Jerk equation

In the last few years, third order explicit autonomous differential equations, known as jerk equations, have generated great interest as they show features of regular and chaotic motion. In this paper, we have modelled chaotic electrostatic ion cyclotron oscillations using a third order nonlinear ordinary differential equation (ODE) and investigated its nonlinear dynamical properties. The nonlinear ODE has been derived for a plasma system from a two fluid model in the presence of a source term, under the influence of an external magnetic field, which is perpendicular to the direction of the wave vector. It is seen that the equation does not require an external forcing term to obtain chaotic behaviour. The stability of the solutions of the equation has been investigated analytically as well as numerically, and the bifurcation diagram obtained shows a number of interesting phenomena for various regimes of parameters. The coexisting attractors as well as their corresponding basins are shown and the phase space portraits at different conditions are obtained numerically and shown here. The results obtained here are in agreement with preliminary experiments conducted for a similar configuration of a plasma system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Effect of trapped electron on the dust ion acoustic waves in dusty plasma using time fractional modified Korteweg-de Vries equation

The time fractional modified Korteweg-de Vries (TFMKdV) equation is solved to study the nonlinear propagation of small but finite amplitude dust ion-acoustic (DIA) solitary waves in un-magnetized dusty plasma with trapped electrons. The plasma is composed of a cold ion fluid, stationary dust grains, and hot electrons obeying a trapped electron distribution. The TFMKdV equation is derived by using the semi-inverse and Agrawal's methods and then solved by the Laplace Adomian decomposition method. Our results show that the amplitude of the DIA solitary waves increases with the increase of time fractional order β, the wave velocity v{sub 0}, and the population of the background free electrons λ. However, it is vice-versa for the deviation from isothermality parameter b, which is in agreement with the result obtained previously.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Nonlinear ion-acoustic waves in a degenerate plasma with nuclei of heavy elements

The ion-acoustic (IA) solitary waves propagating in a fully relativistic degenerate dense plasma (containing relativistic degenerate electron and ion fluids, and immobile nuclei of heavy elements) have been theoretically investigated. The relativistic hydrodynamic model is used to derive the Korteweg-de Vries (K-dV) equation by the reductive perturbation method. The stationary solitary wave solution of this K-dV equation is obtained to characterize the basic features of the IA solitary structures that are found to exist in such a degenerate plasma. It is found that the effects of electron dynamics, relativistic degeneracy of the plasma fluids, stationary nuclei of heavy elements, etc., significantly modify the basic properties of the IA solitary structures. The implications of this results in astrophysical compact objects like white dwarfs are briefly discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Magnetogenesis through convection in barotropic fluids

It is shown that an unmagnetized plasma with non-uniform bulk velocity can generate magnetic fields through consideration of the non-relativistic isentropic two-fluid equations, even when the initial conditions contain with no fields or currents, uniform densities and pressures, and a divergence-free bulk velocity. This effect does not depend on the baroclinicity of the plasma and is therefore relevant even in barotropic flows, where the Biermann battery is absent. It also does not rely on kinetic effects or shear discontinuities. Instead, our magnetogenesis effect arises from convection terms proportional to the electron mass in the generalized Ohm's law. The resulting magnetic fields are typically weak but may still serve as seed fields for dynamo mechanisms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Magnetic island evolution in hot ion plasmas

Effects of finite ion temperature on magnetic island evolution are studied by means of numerical simulations of a reduced set of two-fluid equations which include ion as well as electron diamagnetism in slab geometry. The polarization current is found to be almost an order of magnitude larger in hot than in cold ion plasmas, due to the strong shear of ion velocity around the separatrix of the magnetic islands. As a function of the island width, the propagation speed decreases from the electron drift velocity (for islands thinner than the Larmor radius) to values close to the guiding-center velocity (for islands of order 10 times the Larmor radius). In the latter regime, the polarization current is destabilizing (i.e., it drives magnetic island growth). This is in contrast to cold ion plasmas, where the polarization current is generally found to have a healing effect on freely propagating magnetic island.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Nonlinear theory of magnetohydrodynamic flows of a compressible fluid in the shallow water approximation

Shallow water magnetohydrodynamic (MHD) theory describing incompressible flows of plasma is generalized to the case of compressible flows. A system of MHD equations is obtained that describes the flow of a thin layer of compressible rotating plasma in a gravitational field in the shallow water approximation. The system of quasilinear hyperbolic equations obtained admits a complete simple wave analysis and a solution to the initial discontinuity decay problem in the simplest version of nonrotating flows. In the new equations, sound waves are filtered out, and the dependence of density on pressure on large scales is taken into account that describes static compressibility phenomena. In the equations obtained, the mass conservation law is formulated for a variable that nontrivially depends on the shape of the lower boundary, the characteristic vertical scale of the flow, and the scale of heights at which the variation of density becomes significant. A simple wave theory is developed for the system of equations obtained. All self-similar discontinuous solutions and all continuous centered self-similar solutions of the system are obtained. The initial discontinuity decay problem is solved explicitly for compressible MHD equations in the shallow water approximation. It is shown that there exist five different configurations that provide a solution to the initial discontinuity decay problem. For each configuration, conditions are found that are necessary and sufficient for its implementation. Differences between incompressible and compressible cases are analyzed. In spite of the formal similarity between the solutions in the classical case of MHD flows of an incompressible and compressible fluids, the nonlinear dynamics described by the solutions are essentially different due to the difference in the expressions for the squared propagation velocity of weak perturbations. In addition, the solutions obtained describe new physical phenomena related to the dependence of the height of the free boundary on the density of the fluid. Self-similar continuous and discontinuous solutions are obtained for a system on a slope, and a solution is found to the initial discontinuity decay problem in this case.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Two-fluid description of wave-particle interactions in strong Buneman turbulence

To understand the nature of anomalous resistivity in magnetic reconnection, we investigate turbulence-induced momentum transport and energy dissipation while a plasma is unstable to the Buneman instability in force-free current sheets. Using 3D particle-in-cell simulations, we find that the macroscopic effects generated by wave-particle interactions in Buneman instability can be approximately described by a set of electron fluid equations. We show that both energy dissipation and momentum transport along electric current in the current layer are locally quasi-static, but globally dynamic and irreversible. Turbulent drag dissipates both the streaming energy of the current sheet and the associated magnetic energy. The net loss of streaming energy is converted into the electron component heat conduction parallel to the magnetic field and increases the electron Boltzmann entropy. The growth of self-sustained Buneman waves satisfies a Bernoulli-like equation that relates the turbulence-induced convective momentum transport and thermal momentum transport. Electron trapping and de-trapping drive local momentum transports, while phase mixing converts convective momentum into thermal momentum. The drag acts like a micro-macro link in the anomalous heating processes. The decrease of magnetic field maintains an inductive electric field that re-accelerates electrons, but most of the magnetic energy is dissipated and converted into the component heat of electrons perpendicular to the magnetic field. This heating process is decoupled from the heating of Buneman instability in the current sheets. Ion heating is weak but ions play an important role in assisting energy exchanges between waves and electrons. Cold ion fluid equations together with our electron fluid equations form a complete set of equations that describes the occurrence, growth, saturation and decay of the Buneman instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Magnetohydrodynamic motion of a two-fluid plasma

Here, the two-fluid Maxwell system couples frictionless electron and ion fluids via Maxwell’s equations. When the frequencies of light waves, Langmuir waves, and single-particle cyclotron motion are scaled to be asymptotically large, the two-fluid Maxwell system becomes a fast-slow dynamical system. This fast-slow system admits a formally-exact single-fluid closure that may be computed systematically with any desired order of accuracy through the use of a functional partial differential equation. In the leading order approximation, the closure reproduces magnetohydrodynamics (MHD). Higher order truncations of the closure give an infinite hierarchy of extended MHD models that allow for arbitrary mass ratio, as well as perturbative deviations from charge neutrality. The closure is interpreted geometrically as an invariant slow manifold in the infinite-dimensional two-fluid phase space, on which two-fluid motions are free of high-frequency oscillations. This perspective shows that the full closure inherits a Hamiltonian structure from two-fluid theory. By employing infinite-dimensional Lie transforms, the Poisson bracket for the all-orders closure may be obtained in closed form. Thus, conservative truncations of the single-fluid closure may be obtained by simply truncating the single-fluid Hamiltonian. Moreover, the closed-form expression for the all-orders bracket gives explicit expressions for a number of the full closure’s conservation laws. Notably, the full closure, as well as any of its Hamiltonian truncations, admits a pair of independent circulation invariants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Pressure–Strain Interaction as the Energy Dissipation Estimate in Collisionless Plasma

The dissipative mechanism in weakly collisional plasma is a topic that pervades decades of studies without a consensus solution. We compare several energy dissipation estimates based on energy transfer processes in plasma turbulence and provide justification for the pressure–strain interaction as a direct estimate of the energy dissipation rate. The global and scale-by-scale energy balances are examined in 2.5D and 3D kinetic simulations. We show that the global internal energy increase and the temperature enhancement of each species are directly tracked by the pressure–strain interaction. The incompressive part of the pressure–strain interaction dominates over its compressive part in all simulations considered. The scale-by-scale energy balance is quantified by scale filtered Vlasov–Maxwell equations, a kinetic plasma approach, and the lag dependent von Kármán–Howarth equation, an approach based on fluid models. We find that the energy balance is exactly satisfied across all scales, but the lack of a well-defined inertial range influences the distribution of the energy budget among different terms in the inertial range. Therefore, the widespread use of the Yaglom relation in estimating the dissipation rate is questionable in some cases, especially when the scale separation in the system is not clearly defined. In contrast, the pressure–strain interaction balances exactly the dissipation rate at kinetic scales regardless of the scale separation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Massively parallel axisymmetric fluid model for streamer discharges

A highly parallelizable fluid plasma simulation tool based upon the first-order drift-diffusion equations is discussed. Atmospheric pressure plasmas have densities and gradients that require small element sizes in order to accurately simulate the plasm resulting in computational meshes on the order of millions to tens of millions of elements for realistic size plasma reactors. To enable simulations of this nature, parallel computing is required and must be optimized for the particular problem. Here, a finite-volume, electrostatic drift-diffusion implementation for low-temperature plasma is discussed. The implementation is built upon the Message Passing Interface (MPI) library in C++ using Object Oriented Programming. The underlying numerical method is outlined in detail and benchmarked against simple streamer formation from other streamer codes. Electron densities, electric field, and propagation speeds are compared with the reference case and show good agreement. Convergence studies are also performed showing a minimal space step of approximately 4 μm required to reduce relative error to below 1% during early streamer simulation times and even finer space steps are required for longer times. Additionally, strong and weak scaling of the implementation are studied and demonstrate the excellent performance behavior of the implementation up to 100 million elements on 1024 processors. Lastly, different advection schemes are compared for the simple streamer problem to analyze the influence of numerical diffusion on the resulting quantities of interest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Coupled modes in magnetized dense plasma with relativistic-degenerate electrons

Low frequency electrostatic and electromagnetic waves are investigated in ultra-dense quantum magnetoplasma with relativistic-degenerate electron and non-degenerate ion fluids. The dispersion relation is derived for mobile as well as immobile ions by employing hydrodynamic equations for such plasma under the influence of electromagnetic forces and pressure gradient of relativistic-degenerate Fermi gas of electrons. The result shows the coexistence of shear Alfven and ion modes with relativistically modified dispersive properties. The relevance of results to the dense degenerate plasmas of astrophysical origin (for instance, white dwarf stars) is pointed out with brief discussion on ultra-relativistic and non-relativistic limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Numerical analysis of a mixture of Ar/NH{sub 3} microwave plasma chemical vapor deposition reactor

A two-dimensional fluid model has been used to investigate the properties of plasma in Ar/NH{sub 3} microwave electron cyclotron resonance discharge at low pressure. The electromagnetic field model solved by the three-dimensional Simpson method is coupled to a fluid plasma model. The finite difference method was employed to discrete the governing equations. 40 species (neutrals, radicals, ions, and electrons) are consisted in the model. In total, 75 electron-neutral, 43 electron-ion, 167 neutral-neutral, 129 ion-neutral, 28 ion-ion, and 90 3-body reactions are used in the model. According to the simulation, the distribution of the densities of the considered plasma species has been showed and the mechanisms of their variations have been discussed. It is found that the main neutrals (Ar*, Ar**, NH{sub 3}{sup *}, NH, H{sub 2}, NH{sub 2}, H, and N{sub 2}) are present at high densities in Ar/NH{sub 3} microwave electron cyclotron resonance discharge when the mixing ratio of Ar/NH{sub 3} is 1:1 at 20 Pa. The density of NH is more than that of NH{sub 2} atom. And NH{sub 3}{sup +} are the most important ammonia ions. But the uniformity of the space distribution of NH{sub 3}{sup +} is lower than the other ammonia ions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

ON THE SPATIAL SCALES OF WAVE HEATING IN THE SOLAR CHROMOSPHERE

Dissipation of magnetohydrodynamic (MHD) wave energy has been proposed as a viable heating mechanism in the solar chromospheric plasma. Here, we use a simplified one-dimensional model of the chromosphere to theoretically investigate the physical processes and spatial scales that are required for the efficient dissipation of Alfvén waves and slow magnetoacoustic waves. We consider the governing equations for a partially ionized hydrogen-helium plasma in the single-fluid MHD approximation and include realistic wave damping mechanisms that may operate in the chromosphere, namely, Ohmic and ambipolar magnetic diffusion, viscosity, thermal conduction, and radiative losses. We perform an analytic local study in the limit of small amplitudes to approximately derive the lengthscales for critical damping and efficient dissipation of MHD wave energy. We find that the critical dissipation lengthscale for Alfvén waves depends strongly on the magnetic field strength and ranges from 10 m to 1 km for realistic field strengths. The damping of Alfvén waves is dominated by Ohmic diffusion for weak magnetic field and low heights in the chromosphere, and by ambipolar diffusion for strong magnetic field and medium/large heights in the chromosphere. Conversely, the damping of slow magnetoacoustic waves is less efficient, and spatial scales shorter than 10 m are required for critical damping. Thermal conduction and viscosity govern the damping of slow magnetoacoustic waves and play an equally important role at all heights. These results indicate that the spatial scales at which strong wave heating may work in the chromosphere are currently unresolved by observations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1978

Radiative cooling in shock-heated hydrogen-helium plasmas

Axial and off-axis radiative cooling of cylindrical shock-heated hydrogen-helium plasmas is investigated theoretically and experimentally. The coupled fluid dynamic-radiative transfer equations are solved by a combination of approximation techniques aimed at simplifying the computation of the flux divergence term, namely, the quasi-isothermal approximation and the exponential approximation developed for the solid angle integration. The accuracy of the approximation schemes has been assessed and found acceptable for applying the methods to the rapid computation of the radiatively coupled flow problem. Radiative cooling experiments were conducted in a 6-inch annular arc accelerator shock tube (ANAA) for an initial pressure of 1 torr and shock speeds from 35 to 45 Km/sec. The results indicate that the lateral cooling is small compared with the axial cooling, and that better agreement is achieved between the data and the theoretical results by inclusion of the lateral temperature gradient.

Poon, P. T. Y.↗