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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Longitudinal oscillations in a non-uniform spatially dispersive plasma

Longitudinal oscillations of the electron fluid in the hydrodynamic model of a metal are examined with pressure effects taken into account. It is well-known that this entails spatial dispersion. The equilibrium electron number density is taken to be non-uniform and a non-self-adjoint fourth order differential equation obeyed by the electric potential is derived. A velocity potential necessary for the description of sound waves is introduced in the standard fashion and the generalized version of Bloch orthogonality appropriate to a non-uniform background is deduced. We observe a duality between electric and velocity potentials in the sense that the respective differential operators are adjoint to each other. The spectrum is calculated in the special case of an exponential profile for the equilibrium electron number density. The surface plasmons are connected with the analytic properties of the scattering amplitude in the complex plane. The phase shift at threshold is expressed in terms of the number of surface plasmon modes via an expression reminiscent of Levinson’s statement in quantum mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Relativistic modeling capabilities in PERSEUS extended MHD simulation code for HED plasmas

We discuss the incorporation of relativistic modeling capabilities into the PERSEUS extended MHD simulation code for high-energy-density (HED) plasmas, and present the latest hybrid X-pinch simulation results. The use of fully relativistic equations enables the model to remain self-consistent in simulations of such relativistic phenomena as X-pinches and laser-plasma interactions. By suitable formulation of the relativistic generalized Ohm’s law as an evolution equation, we have reduced the recovery of primitive variables, a major technical challenge in relativistic codes, to a straightforward algebraic computation. Our code recovers expected results in the non-relativistic limit, and reveals new physics in the modeling of electron beam acceleration following an X-pinch. Through the use of a relaxation scheme, relativistic PERSEUS is able to handle nine orders of magnitude in density variation, making it the first fluid code, to our knowledge, that can simulate relativistic HED plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Three-fluid, three-dimensional magnetohydrodynamic solar wind model with eddy viscosity and turbulent resistivity

We have developed a three-fluid, three-dimensional magnetohydrodynamic solar wind model that incorporates turbulence transport, eddy viscosity, turbulent resistivity, and turbulent heating. The solar wind plasma is described as a system of co-moving solar wind protons, electrons, and interstellar pickup protons, with separate energy equations for each species. Numerical steady-state solutions of Reynolds-averaged solar wind equations coupled with turbulence transport equations for turbulence energy, cross helicity, and correlation length are obtained by the time relaxation method in the corotating with the Sun frame of reference in the region from 0.3 to 100 AU (but still inside the termination shock). The model equations include the effects of electron heat conduction, Coulomb collisions, photoionization of interstellar hydrogen atoms and their charge exchange with the solar wind protons, turbulence energy generation by pickup protons, and turbulent heating of solar wind protons and electrons. The turbulence transport model is based on the Reynolds decomposition and turbulence phenomenologies that describe the conversion of fluctuation energy into heat due to a turbulent cascade. In addition to using separate energy equations for the solar wind protons and electrons, a significant improvement over our previous work is that the turbulence model now uses an eddy viscosity approximation for the Reynolds stress tensor and the mean turbulent electric field. The approximation allows the turbulence model to account for driving of turbulence by large-scale velocity gradients. Using either a dipole approximation for the solar magnetic field or synoptic solar magnetograms from the Wilcox Solar Observatory for assigning boundary conditions at the coronal base, we apply the model to study the global structure of the solar wind and its three-dimensional properties, including embedded turbulence, heating, and acceleration throughout the heliosphere. The model results are compared with plasma and magnetic field observations on WIND, Ulysses, and Voyager 2 spacecraft.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Linear stability analysis of the Vlasov-Poisson equations in high density plasmas in the presence of crossed fields and density gradients

The equations for the single-particle orbits in a nonneutral high density plasma in the presence of inhomogeneous crossed fields are obtained. Using these orbits, the linearized Vlasov equation is solved as an expansion in the orbital radii in the presence of inhomogeneities and density gradients. A model distribution function is introduced whose cold-fluid limit is exactly the same as that used in many previous studies of the cold-fluid equations. This model function is used to reduce the linearized Vlasov-Poisson equations to a second-order ordinary differential equation for the linearized electrostatic potential whose eigenvalue is the perturbation frequency.

Kaup, D. J.↗
DOE OSTI2022

Modelling of electrostatic ion-scale turbulence in divertor tokamaks with the gyrokinetic code $\mathrm{COGENT}$

Here, continuum gyrokinetic simulations of electrostatic ion scale turbulence are presented for the case of a diverted (single-null) tokamak geometry. The simulation model, implemented in the finite-volume code COGENT, solves the long-wavelength limit of the full-F gyrokinetic equation for ion species coupled to a vorticity equation for electrostatic potential variations, where a fluid model is used for an electron response. The model describes the ion scale ion temperature gradient (ITG) and resistive drift modes as well as neoclassical ion physics effects. Different turbulence regimes are observed depending on the plasma profiles, and the roles of a self-consistent background electric field and an X-point geometry are explored. In particular, increasing the pedestal density gradient and the corresponding radial electric field is demonstrated to suppress the ITG turbulence, whereas the same edge plasma background can still be destabilized by the resistive modes when the plasma resistivity is increased. The effects of X-point geometry are assessed by comparing cross-separatrix simulations with counterpart calculations performed for a toroidal annulus geometry. For the simulation parameters considered, similar global behaviour is observed in both cases, whereas strong local suppression of turbulence fluctuations is demonstrated near the X-point for the case of a single-null geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Transition between inverse and direct energy cascades in multiscale optical turbulence

Transition between inverse and direct energy cascades in multiscale optical turbulence. Multiscale turbulence naturally develops and plays an important role in many fluid, gas, and plasma phenomena. Statistical models of multiscale turbulence usually employ Kolmogorov hypotheses of spectral locality of interactions (meaning that interactions primarily occur between pulsations of comparable scales) and scale-invariance of turbulent pulsations. However, optical turbulence described by the nonlinear Schrodinger equation exhibits breaking of both the Kolmogorov locality and scale-invariance. A weaker form of spectral locality that holds for multi-scale optical turbulence enables a derivation of simplified evolution equations that reduce the problem to a single scale modeling. We present the derivation of these equations for Kerr media with random inhomogeneities. Then, we find the analytical solution that exhibits a transition between inverse and direct energy cascades in optical turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1975

Multifluids description of dynamics of upper atmosphere

A multifluids model to investigate ionospheric dynamics was established on kinetic theory. Its resultant equations are used to examine the following dynamic problems in the gamma region of 80-2000 Km of the ionosphere: (1) propagation of acoustic modes in the 500-2,000 Km of the ionosphere (two fluid model); (2) the relation between the cross field plasma drift instabilities and type I and type II ionospheric irregularities; and (3) time dependent neutral wind structure and horizontal pressure gradient.

Wu, S. T.↗
NASA NTRS2010

Multi-Fluid Block-Adaptive-Tree Solar Wind Roe-Type Upwind Scheme: Magnetospheric Composition and Dynamics During Geomagnetic Storms, Initial Results

The magnetosphere contains a significant amount of ionospheric O{+}, particularly during geomagnetically active times. The presence of ionospheric plasma in the magnetosphere has a notable impact on magnetospheric composition and processes. We present a new multifluid MHD version of the BATS-R-US model of the magnetosphere to track the fate and consequences of ionospheric outflow. The multi-fluid MHD equations are presented as are the novel techniques for overcoming the formidable challenges associated with solving them. Our new model is then applied to the May 4, 1998 and March 31, 2001 geomagnetic storms. The results are juxtaposed with traditional single- fluid MHD and multispecies MHD simulations from a previous study, thereby allowing us to assess the benefits of using a more complex model with additional physics. We find that our multi-fluid MHD model (with outflow) gives comparable results to the multi-species MHD model (with outflow), including a more strongly negative Dst, reduced CPCP, and a drastically improved magnetic field at geosynchronous orbit, as compared to single-fluid MHD with no outflow. Significant differences in composition and magnetic field are found between the multi-species and multi-fluid approach further away from the Earth. We further demonstrate the ability to explore pressure and bulk velocity differences between H{+} and O(+}, which is not possible when utilizing the other techniques considered.

Gkocer, A.↗
DOE OSTI2014

Manufactured solutions for the three-dimensional Euler equations with relevance to Inertial Confinement Fusion

We present a set of manufactured solutions for the three-dimensional (3D) Euler equations. The purpose of these solutions is to allow for code verification against true 3D flows with physical relevance, as opposed to 3D simulations of lower-dimensional problems or manufactured solutions that lack physical relevance. Of particular interest are solutions with relevance to Inertial Confinement Fusion (ICF) capsules. While ICF capsules are designed for spherical symmetry, they are hypothesized to become highly 3D at late time due to phenomena such as Rayleigh–Taylor instability, drive asymmetry, and vortex decay. ICF capsules also involve highly nonlinear coupling between the fluid dynamics and other physics, such as radiation transport and thermonuclear fusion. The manufactured solutions we present are specifically designed to test the terms and couplings in the Euler equations that are relevant to these phenomena. Example numerical results generated with a 3D Finite Element hydrodynamics code are presented, including mesh convergence studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Hamiltonian fluid closures of the Vlasov-Ampère equations: From water-bags to N moment models

Moment closures of the Vlasov-Ampère system, whereby higher moments are represented as functions of lower moments with the constraint that the resulting fluid system remains Hamiltonian, are investigated by using water-bag theory. The link between the water-bag formalism and fluid models that involve density, fluid velocity, pressure and higher moments is established by introducing suitable thermodynamic variables. The cases of one, two, and three water-bags are treated and their Hamiltonian structures are provided. In each case, we give the associated fluid closures and we discuss their Casimir invariants. We show how the method can be extended to an arbitrary number of fields, i.e., an arbitrary number of water-bags and associated moments. The thermodynamic interpretation of the resulting models is discussed. Finally, a general procedure to derive Hamiltonian N-field fluid models is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The fastVFP code for solution of the Vlasov–Fokker–Planck equation

Abstract We describe the fast VFP code for solution of the Vlasov–Fokker–Planck equation for non-local electron transport and the generation of magnetic field, especially for application to laser-produced plasmas. We describe the essential features of the code that make it fast and robust and suitable for inclusion as a transport package in a fluid simulation. We present a few sample results that demonstrate the abilities of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Cylindrical Damped Solitary Propagation in Superthermal Plasmas

Wave properties of damped solitons in a collisional unmagnetized four-components dusty fluid plasma system contains superthermal distributed electrons, mobile ions, and negative-positive dusty grains have been examined. To study dissipative DIA mode properties, a reductive perturbation (RP) analysis is used under convenient geometrical coordinate transformation, three-dimensional damped Kadomtsev–Petviashvili (3D-CDKP) equation in cylindrical coordinates is obtained. Effects of collisional parameters on damped soliton pulse structures are studied. More specifically, the impact of axial, radial, and polar coordinates with the time on solitary propagation are examined. This investigation may be viable in plasma of the Earth’s mesosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

Physics of Reconnection and MMS Mission

Reconnection is the most important process driving the Earth's magnetosphere. Key to the success of the MMS science plan is the coupling of theory and observation. Determining the kinetic processes occurring in the diffusion region and physical parameters that control the rate of magnetic reconnection are among primary objectives of the MMS mission. Analysis of the role played by particle inertial effects in the diffusion region where the plasma is unmagnetized will be presented. The reconnection electric field in he diffusion region is supported primarily by particle non-gyrotropic effects. At the quasi-steady stage the reconnection electric field serves to accelerate and heat the incoming plasma population to maintain the current flow in the diffusion region the pressure balance. The primary mechanism controlling the dissipation in the vicinity of the reconnection site is incorporated into the fluid description in terms of non-gyrotropic corrections to the. induction and energy equations. The results of kinetic and fluid simulations illustrating the physics of magnetic reconnection will be presented. We will dem~:tistrate that kinetic nongyrotropic effects can significantly alter the global magnetosphere evolution and location of reconnection sites.

Kuznetsova, M. M.↗
DOE OSTI2013

Modeling of plasma chemistry in an atmospheric pressure Ar/NH{sub 3} cylindrical dielectric barrier discharge described using the one-dimensional fluid model

The keynote of our research is to study the gas phase chemistry in an atmospheric pressure Ar/NH{sub 3} cylindrical dielectric barrier discharge, which is very important to produce the iron-nitride magnetic fluid. For this purpose, a home-made one dimensional fluid model with the Scharfetter-Gummel method has been developed. The equations solved are the particle balances, assuming a drift-diffusion approximation for the fluxes, and the electron energy equation. The self-consistent electric field is obtained by the simultaneous solution of Poisson's equation. The simulations were carried out for the different ammonia concentrations (2%, 3.5%, and 7%), at a voltage of 1 kV, and a driving frequency of 20 kHz. It concluded that the major ion products of Ar are Ar{sup +} and Ar{sub 2}{sup +}. Ar{sup +} is the most important positive ions, followed by Ar{sub 2}{sup +}. It is shown that the NH{sup +} density is smaller than that of the other ammonia ions. The density of NH{sub 4}{sup +} is more than that of the other ammonia ions when the ammonia concentration increased. The diffuse mode can be established after the discharge was ignited, and the mode changes to filamentary mode with an increase in ammonia concentration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transient behavior of drift and ionization in atmospheric pressure nitrogen discharge

The fluid models are frequently used to describe a non-thermal plasma such as a streamer discharge. The required electron transport data and rate coefficients for the fluid model are parametrized using the local field approximation (LFA) in first order models and the local-mean-energy approximation (LMEA) in second order models. We performed Monte Carlo simulations in Nitrogen gas with step changes in the E/N (reduced electric field) to study the behavior of the transport properties in the transient phase. During the transient phase of the simulation, we extract the instantaneous electron mean energy, which is different from the steady state mean electron energy, and the corresponding transport parameters and rate coefficients. Our results indicate that the mean electron energy is not a suitable parameter for mobility/drift of electrons due to big difference in momentum relaxation and energy relaxation. Furthermore, the high energy threshold rates such as ionization show a strong correlation to mean electron energy. In second order models where the energy-balance equation is solved, we suggest that it would rather be appropriate to use the local electric field to find electron drift velocity in gases such as Nitrogen and the local mean electron energy to determine the ionization and excitation rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Axial magnetic field and toroidally streaming fast ions in the dense plasma focus are natural consequences of conservation laws in the curved axisymmetric geometry of the current sheath

Direct measurement of axial magnetic field in the PF-1000 dense plasma focus (DPF), and its reported correlation with neutron emission, call for a fresh look at previous reports of existence of axial magnetic field component in the DPF from other laboratories, and associated data suggesting toroidal directionality of fast ions participating in fusion reactions, with a view to understand the underlying physics. In this context, recent work dealing with application of the hyperbolic conservation law formalism to the DPF is extended in this paper to a curvilinear coordinate system, which reflects the shape of the DPF current sheath. Locally unidirectional shock propagation in this coordinate system enables construction of a system of 7 one-dimensional hyperbolic conservation law equations with geometric source terms, taking into account all the components of magnetic field and flow velocity. Rankine-Hugoniot jump conditions for this system lead to expressions for the axial magnetic field and three components of fluid velocity having high ion kinetic energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗