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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2010

Modification of the Plasma in the Near-Vicinity of Enceladus by the Enveloping Dust

The plasma near Saturn's equator is quasi-corotating, but those fluid elements entering the near-vicinity of the moon Enceladus become uniquely modified. Besides the solid body, the Moon has a surrounding dust envelop that we show herein to be detected approx.20 Enceladus radii (1 R(sub E) = 252 km) both north and south of the body. Previous reports indicate that co rotating plasma slows down substantially in the near-vicinity of En eel ad us. We show herein that the commencement of this plasma slow down matches closely with Cassini's entry into the dense portions of the enveloping dust in the northern hemisphere above the Moon. We also examine in detail the source of the dust about 400 km above the south polar fissures. We find that a large positive potential must exist between the south pole of the moon and the spacecraft to account for ions streaming away from the pole on connecting magnetic field lines

Farrell, W. M.↗
NASA NTRSDate not supplied

The Circulation of the Plasmasphere Fluid during the Erosion Event on September 8, 2017

A strong solar wind pressure pulse triggered the magnetic storm on September 7, 2017. Near the end of September 7, the z-component of the interplanetary magnetic field (IMF Bz) dropped from 9 to -10 nT in 30 min. The IMF Bz remained at the level of -10 nT for 2 hours and then had another rapid drop to -31 nT in 30 min. The sudden plunge of IMF Bz and the associated strong convection electric field stirred up the storm main phase with Dst falling from 0 to -122 nT from 2200 UT on September 7 to 0200 UT on the 8th. Severe plasmasphere erosion was observed on September 8 by multiple spacecraft, such as the Van Allen Probes and the Arase satellite. In this study, we examine the fate of the eroded plasmasphere particles during this event by model simulation as well as satellite data analysis. The simulation tool we use is the Space Weather Modeling Framework (SWMF)/Block-Adaptive Tree Solarwind Roe-type Upwind Scheme (BATS-R-US) model coupled with the Comprehensive Inner Magnetosphere-Ionosphere (CIMI) model. One of the distinctive capabilities of the SWMF/BATSRUS-CIMI model is that it treats the cold plasmas in the plasmasphere as a separate fluid in the MHD equations. As a result, the transport and circulation of the plasmasphere fluid in the global magnetosphere can be traced and the impacts of this cold fluid on the global magnetosphere can be evaluated. In this paper, we will show how the drainage plume is formed during the storm and how the plasmasphere fluid is transported to the flank and lobe regions and eventually to the plasma sheet and reenters into the plasmasphere. We will validate our simulation by plasmasphere signatures observed in both the inner and outer magnetosphere.

Mei-Ching Fok↗
DOE OSTI2014

The full Zakharov equations in nonextensive q-plasma

The wave-wave interactions among electromagnetic mode, Langmuir and ion-sound modes are theoretically investigated in a nonextensive plasma, with the help of the two-fluid theory and the two-scale method. The full Zakharov equations for the nonextensive distribution are derived and analyzed numerically. The results show that the collapsing of the three waves would be affected by the nonextensive parameter q. When q increasing, the collapsing would become faster and the rate of density perturbation would be greater. And the localized density structure would form a plasma channel much easier.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Equation of state of dense neon and krypton plasmas in the partial ionization regime

The compression behaviors of dense neon and krypton plasmas over a wide pressure-temperature range are investigated by self-consistent fluid variational theory. The ionization degree and equation of state of dense neon and krypton are calculated in the density-temperature range of 0.01–10 g/cm{sup 3} and 4–50 kK. A region of thermodynamic instability is found which is related to the plasma phase transition. The calculated shock adiabat and principal Hugoniot of liquid krypton are in good agreement with available experimental data. The predicted results of shock-compressed liquid neon are presented, which provide a guide for dynamical experiments or numerical first-principle calculations aimed at studying the compression properties of liquid neon in the partial ionization regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is &#x3B3; &#x3C4; A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of &#x3B3; &#x3C4; A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A note on the cylindrical solitary waves in an electron-acoustic plasma with vortex electron distribution

In the present work, we consider the propagation of nonlinear electron-acoustic non-planar waves in a plasma composed of a cold electron fluid, hot electrons obeying a trapped/vortex-like distribution, and stationary ions. The basic nonlinear equations of the above described plasma are re-examined in the cylindrical coordinates through the use reductive perturbation method in the long-wave approximation. The modified cylindrical Korteweg-de Vries equation with fractional power nonlinearity is obtained as the evolution equation. Due to the nature of nonlinearity, which is fractional, this evolution equation cannot be reduced to the conventional Korteweg–de Vries equation. An analytical solution to the evolution equation, by use of the method developed by Demiray [Appl. Math. Comput. 132, 643 (2002); Comput. Math. Appl. 60, 1747 (2010)] and a numerical solution by employing a spectral scheme are presented and the results are depicted in a figure. The numerical results reveal that both solutions are in good agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1992

Drift Kelvin-Helmholtz instabilities in space plasmas

Drift Kelvin-Helmholtz instabilities of a finite-beta plasma in equilibrium electric and magnetic fields which are perpendicular to each other are studied using two fluid equations. Three types of these instabilities are considered including the magnetosonic instability of a finite beta-homogeneous plasma, the electrostatic drift instability of an inhomogeneous low-beta plasma, and the magneto-acoustic instability of a high-beta inhomogeneous isothermal plasma. It is shown that the electric field has either stabilizing or destabilizing effect depending on conditions under consideration.

Sharma, Avadhesh C.↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2016

Wakes and precursor soliton excitations by a moving charged object in a plasma

We study the evolution of nonlinear ion acoustic wave excitations due to a moving charged source in a plasma. Our numerical investigations of the full set of cold fluid equations go beyond the usual weak nonlinearity approximation and show the existence of a rich variety of solutions including wakes, precursor solitons, and “pinned” solitons that travel with the source velocity. These solutions represent a large amplitude generalization of solutions obtained in the past for the forced Korteweg deVries equation and can find useful applications in a variety of situations in the laboratory and in space, wherever there is a large relative velocity between the plasma and a charged object.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

A high-order relativistic two-fluid electrodynamic scheme with consistent reconstruction of electromagnetic fields and a multidimensional Riemann solver for electromagnetism

In various astrophysics settings it is common to have a two-fluid relativistic plasma that interacts with the electromagnetic field. While it is common to ignore the displacement current in the ideal, classical magnetohydrodynamic limit, when the flows become relativistic this approximation is less than absolutely well-justified. In such a situation, it is more natural to consider a positively charged fluid made up of positrons or protons interacting with a negatively charged fluid made up of electrons. The two fluids interact collectively with the full set of Maxwell's equations. As a result, a solution strategy for that coupled system of equations is sought and found here. Our strategy extends to higher orders, providing increasing accuracy. The primary variables in the Maxwell solver are taken to be the facially-collocated components of the electric and magnetic fields. Consistent with such a collocation, three important innovations are reported here. The first two pertain to the Maxwell solver. In our first innovation, the magnetic field within each zone is reconstructed in a divergence-free fashion while the electric field within each zone is reconstructed in a form that is consistent with Gauss' law. In our second innovation, a multidimensionally upwinded strategy is presented which ensures that the magnetic field can be updated via a discrete interpretation of Faraday's law and the electric field can be updated via a discrete interpretation of the generalized Ampere's law. This multidimensional upwinding is achieved via a multidimensional Riemann solver. The multidimensional Riemann solver automatically provides edge-centered electric field components for the Stokes law-based update of the magnetic field. It also provides edge-centered magnetic field components for the Stokes law-based update of the electric field. The update strategy ensures that the electric field is always consistent with Gauss' law and the magnetic field is always divergence-free. This collocation also ensures that electromagnetic radiation that is propagating in a vacuum has both electric and magnetic fields that are exactly divergence-free. Coupled relativistic fluid dynamic equations are solved for the positively and negatively charged fluids. The fluids' numerical fluxes also provide a self-consistent current density for the update of the electric field. Our reconstruction strategy ensures that fluid velocities always remain sub-luminal. Our third innovation consists of an efficient design for several popular IMEX schemes so that they provide strong coupling between the finite-volume-based fluid solver and the electromagnetic fields at high order. This innovation makes it possible to efficiently utilize high order IMEX time update methods for stiff source terms in the update of high order finite-volume methods for hyperbolic conservation laws. We also show that this very general innovation should extend seamlessly to Runge–Kutta discontinuous Galerkin methods. The IMEX schemes enable us to use large CFL numbers even in the presence of stiff source terms. Several accuracy analyses are presented showing that our method meets its design accuracy in the MHD limit as well as in the limit of electromagnetic wave propagation. Several stringent test problems are also presented. We also present a relativistic version of the GEM problem, which shows that our algorithm can successfully adapt to challenging problems in high energy astrophysics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Modifications produced on a large magnetized plasma column by a floating end-plate that is partially emissive: Experiment and theory

An experiment is performed on a large plasma device operated by the Basic Plasma Science Facility at the University of California, Los Angeles, in which an electrically floating structure is placed near the end of the 18-m magnetized plasma column. The structure consists of a flat carbon plate that acts as a mask for a smaller, ring-shaped LaB6 emissive surface whose temperature can be externally controlled. This configuration has been previously used to study electron heat transport and pressure-driven avalanches [B. Van Compernolle and G. J. Morales, Phys. Plasmas 24, 112302 (2017)] by biasing the LaB6 ring-cathode with respect to a distant anode in a cold afterglow plasma. In contrast, the present study is performed during the active portion of the steady-state discharge in which the nominal plasma parameters are determined by the injection of an electron beam from a BaO cathode at the opposite end. It is found that, even without an applied bias on the LaB6 cathode, the self-consistent potential and current profiles are modified near the end plate as the LaB6 temperature is increased, resulting in density increases on the field lines in contact with the ring-cathode. In the absence of enhanced ionization, at the largest cathode temperatures, the ambient density can be doubled. A theoretical model is presented that provides a quantitative explanation for the experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Multifluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. The interstellar pickup protons are treated in the model as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. The simulation assumes that the background magnetic field on the Sun is either a dipole (aligned or tilted with respect to the solar rotation axis) or one that is deduced from solar magnetograms.

Goldstein, Melvyn L.↗
DOE OSTI2011

Traveling Wave Solutions of the Gardner Equation and Motion of Plane Curves Governed by the mKdV Flow

The Gardner equation is well-known in the mathematical literature since the late sixties of 20th century. Initially, it appeared in the context of the construction of local conservation laws admitted by the KdV equation. Later on, the Gardner equation was generalized and found to be applicable in various branches of physics (solid-state and plasma physics, fluid dynamics and quantum field theory). In this paper, we examine the travelling wave solutions of the Gardner equation and derive the full set of solutions to the corresponding reduced equation in terms of Weierstrass and Jacobi elliptic functions. Then, we use the travelling wave solutions of the focusing mKdV equation and obtain in explicit analytic form exact solutions of a special type of plane curve flow, known as the mKdV flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Self-focusing of circularly polarized laser pulse propagating through a magnetized non-Maxwellian plasma

Self-focusing of an intense circularly polarized laser pulse propagating through a magnetized non-Maxwellian plasma is investigated. Based on a relativistic two-fluid model, nonlinear equation describing dynamics of the slowly varying amplitude is obtained. The evolution of laser spot size is studied and effect of non-Maxwellian distribution of charge density on the spot size is considered. It is shown that the existence of super-thermal particles leads to the enhancement of the self-focusing quality of plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

RAYLEIGH-TAYLOR INSTABILITY IN PARTIALLY IONIZED COMPRESSIBLE PLASMAS

We study the modification of the classical criterion for the linear onset and growing rate of the Rayleigh-Taylor instability (RTI) in a partially ionized plasma in the two-fluid description. The plasma is composed of a neutral fluid and an electron-ion fluid, coupled by means of particle collisions. The governing linear equations and appropriate boundary conditions, including gravitational terms, are derived and applied to the case of the RTI in a single interface between two partially ionized plasmas. The limits of collisionless, no gravity, and incompressible fluids are checked before addressing the general case. We find that both compressibility and ion-neutral collisions lower the linear growth rate, but do not affect the critical threshold of the onset of the RTI. The configuration is always unstable when a lighter plasma is below a heavier plasma regardless the value of the magnetic field strength, the ionization degree, and the ion-neutral collision frequency. However, ion-neutral collisions have a strong impact on the RTI growth rate, which can be decreased by an order of magnitude compared to the value in the collisionless case. Ion-neutral collisions are necessary to accurately describe the evolution of the RTI in partially ionized plasmas such as prominences. The timescale for the development of the instability is much longer than in the classical incompressible fully ionized case. This result may explain the existence of prominence fine structures with life times of the order of 30 minutes. The timescales derived from the classical theory are about one order of magnitude shorter and incompatible with the observed life times.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Modulation instability of an intense laser beam in the hot magnetized electron-positron plasma in the quasi-neutral limit

The aim of the present study is to investigate the problem of modulation instability of an intense laser beam in the hot magnetized electron-positron plasma. Propagation of the intense circularly polarized laser beam along the external magnetic field is studied using a relativistic fluid model. A nonlinear equation describing the interaction of the laser pulse with the magnetized hot pair plasma is derived based on the quasi-neutral approximation, which is valid for the hot plasma. Also, the nonlinear dispersion equation for the hot plasma is obtained. The growth rate of the instability is calculated and its dependence on temperature and external magnetic field are considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theory of ignition and burn propagation in inertial fusion implosions

A detailed analytic model is presented here to investigate the physics of burn propagation in inertially confined plasmas. The onset of ignition and burn propagation occurs when alpha heating of the hot spot causes rapid ablation of shell mass into the hot spot. This allows large energy gains to be achieved since most of the fuel mass is located in the shell. Here, we first present a comprehensive review of previous analytic models that have been used to describe the physics of hot-spot evolution and ignition; we then show that a proper description of a propagating burn wave requires a comprehensive model of hot spot and shell evolution that includes proper mass conservation in the shell, fusion reactivity, and fuel depletion. The analytic theory is in good agreement with detailed radiation-hydrodynamic simulations that predict the onset of burn propagation as occurring when the yield enhancement caused by alpha heating is between 15- and 25-fold, $f_α$ ~ 1.4, where $f_α$ = alpha energy deposited/hot-spot energy at bang time, and the hot-spot burnup fraction is approximately 2%. We show that the definition of ignition is not sensitive to the alpha-particle stopping power nor asymmetries provided that the absorbed fraction of alpha particles $θ_α$ is correctly accounted for. Finally, we use the results of 2-D simulations to show that even when $θ_α$ is small and unknown (as is true in hot spots with mid modes that have significant leakage of alpha particles into the surrounding cold bubbles), one can still relate the experimentally measureable parameter $\chi^{53}_α$ to the yield amplification and the burning-plasma parameter $Q^{\text{hs}}_α$ = alpha energy deposited/total input work delivered to the hot spot

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗