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 451 records · Page 25

Two-dimensional hybrid Monte Carlo–fluid modelling of dc glow discharges: Comparison with fluid models, reliability, and accuracy

Two-dimensional hybrid Monte Carlo–fluid numerical code is developed and applied to model the dc glow discharge. The model is based on the separation of electrons into two parts: the low energetic (slow) and high energetic (fast) electron groups. Ions and slow electrons are described within the fluid model using the drift-diffusion approximation for particle fluxes. Fast electrons, represented by suitable number of super particles emitted from the cathode, are responsible for ionization processes in the discharge volume, which are simulated by the Monte Carlo collision method. Electrostatic field is obtained from the solution of Poisson equation. The test calculations were carried out for an argon plasma. Main properties of the glow discharge are considered. Current-voltage curves, electric field reversal phenomenon, and the vortex current formation are developed and discussed. The results are compared to those obtained from the simple and extended fluid models. Contrary to reports in the literature, the analysis does not reveal significant advantages of existing hybrid methods over the extended fluid model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Fluid model of the sheath in front of a floating electrode immersed in a magnetized plasma with oblique magnetic field: Some comments on ion source terms and ion temperature effects

A one-dimensional fluid model of the magnetized plasma-wall transition region in front of a floating electrode immersed in a magnetized plasma with oblique magnetic field is presented. The Boltzmann relation is assumed for the electrons, while the positive ions obey the ion continuity and momentum exchange equation. The ions are assumed to be isothermal. By comparison with a two-fluid model, it is shown that assuming the Boltzmann relation for the electrons implies that there is no creation or annihilation of the electrons. Consequently, there should not be any creation and annihilation of the positive ions either. The models that assume the Boltzmann relation for the electrons and a non-zero ion source term at the same time are therefore inconsistent, but such models have nevertheless been used extensively by many authors. So, in this work, an extensive comparison of the results obtained using the zero source term on one hand and three different non-zero source terms on the other hand is made. Four different ion source terms are considered in total: the zero source term and three different non-zero ion source terms. When the zero source term is used, the model becomes very sensitive to the boundary conditions, and in some cases, the solutions exhibit large amplitude oscillations. If any of the three non-zero ion source terms is used, those problems are eliminated, but also the consistency of the model is broken. The model equations are solved numerically in the entire magnetized plasma-wall transition region. For zero ion temperature, the model can be solved even if a very small ion velocity is selected as a boundary condition. For finite ion temperature, the system of equations becomes stiff, unless the ion velocity at the boundary is increased slightly above the ion thermal velocity. A simple method how to find a solution with a very small ion velocity at the boundary also for finite ion temperature in the entire magnetized plasma-wall transition region is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Final Technical Report on “Center for Nonlinear Simulation of Energetic Particles in Burning Plasmas (CSEP)”

The gyrokinetic turbulence code GEM is used to study reverse shear Alfven eigenmodes in the DIII-D tokamak, cylindrical tearing mode, the nonlinear saturation mechanism of Alfven eigenmodes, and verify other gyrokinetic codes. It is used to study the stability of Alfven waves in the presence of fusion product alpha particles in ITER. New algorithm features and physics capabilities are added to GEM, including a novel discretization method for the Poisson equation, a vorticity approach to the fluid electron model and a low-n solver in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-dimensional modelling of the solar wind-LISM interaction including neutrals: A Boltzmann equation approach

The importance of interstellar neutrals in understanding and modelling the global interaction of the solar wind with the local interstellar medium is becoming increasingly apparent. Unfortunately the self-consistent inclusion of a neutral interstellar component into time-dependent, dynamical models is formidably difficult due to the extremely large mean free paths associated with the neutrals and the creation of essentially different neutral distributions from different interaction regions of the solar wind and LISM. In full generality, one has to address the problem by treating the neutrals kinetically with the appropriate extinction and creation source terms. In this paper, a limited set of simulations will be presented in which the solar wind and interstellar plasma is described as a 2D fully compressible time-dependent fluid while the interstellar neutral distribution is derived by solving the appropriate Boltzmann equation directly.

Zank, G. P.↗

Linearly Coupled Electrostatic and Shear Alfven Waves in Dense Plasma in the Presence of Stationary Dust

Low frequency electrostatic and electromagnetic waves in a dense magnetoplasma are studied. The dispersive contribution of electron quantum effects in an electron-ion plasma in the presence of positively or negatively charged dust particles in the background is emphasized. By employing the quantum hydrodynamic model, a linear dispersion relation is derived which shows coupling of electrostatic and shear Alfven modes which shows influence of electron quantum effects and dust density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ion acoustic solitary waves and double layers in a plasma with two temperature electrons featuring Tsallis distribution

The propagation properties of large amplitude ion acoustic solitary waves (IASWs) are studied in a plasma containing cold fluid ions and multi-temperature electrons (cool and hot electrons) with nonextensive distribution. Employing Sagdeev pseudopotential method, an energy balance equation has been derived and from the expression for Sagdeev potential function, ion acoustic solitary waves and double layers are investigated numerically. The Mach number (lower and upper limits) for the existence of solitary structures is determined. Positive as well as negative polarity solitary structures are observed. Further, conditions for the existence of ion acoustic double layers (IADLs) are also determined numerically in the form of the critical values of q{sub c}, f and the Mach number (M). It is observed that the nonextensivity of electrons (via q{sub c,h}), concentration of electrons (via f) and temperature ratio of cold to hot electrons (via β) significantly influence the characteristics of ion acoustic solitary waves as well as double layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The role of electron equation of state in heating partition of protons in a collisionless plasma

One of the outstanding questions related to the solar wind is the heating of solar wind plasma. Addressing this question requires a self consistent treatment of the kinetic physics of a collisionless plasma. A hybrid code (with particle ions and fluid electrons) is one of the most convenient computational tools, which allows us to explore self consistent ion kinetics, while saving us computational time as compared to the full particle in cell codes. A common assumption used in hybrid codes is that of isothermal electrons. In this paper, we discuss the role that the equation of state for electrons could potentially play in determining the ion kinetics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear electrostatic excitations of charged dust in degenerate ultra-dense quantum dusty plasmas

The linear and nonlinear properties of low-frequency electrostatic excitations of charged dust particles (or defects) in a dense collisionless, unmagnetized Thomas-Fermi plasma are investigated. A fully ionized three-component model plasma consisting of electrons, ions, and negatively charged massive dust grains is considered. Electrons and ions are assumed to be in a degenerate quantum state, obeying the Thomas-Fermi density distribution, whereas the inertial dust component is described by a set of classical fluid equations. Considering large-amplitude stationary profile travelling-waves in a moving reference frame, the fluid evolution equations are reduced to a pseudo-energy-balance equation, involving a Sagdeev-type potential function. The analysis describes the dynamics of supersonic dust-acoustic solitary waves in Thomas-Fermi plasmas, and provides exact predictions for their dynamical characteristics, whose dependence on relevant parameters (namely, the ion-to-electron Fermi temperature ratio, and the dust concentration) is investigated. An alternative route is also adopted, by assuming weakly varying small-amplitude disturbances off equilibrium, and then adopting a multiscale perturbation technique to derive a Korteweg-de Vries equation for the electrostatic potential, and finally solving in terms for electric potential pulses (electrostatic solitons). A critical comparison between the two methods reveals that they agree exactly in the small-amplitude, weakly superacoustic limit. The dust concentration (Havnes) parameter h=Z{sub d0}n{sub d0}/n{sub e0} affects the propagation characteristics by modifying the phase speed, as well as the electron/ion Fermi temperatures. Our results aim at elucidating the characteristics of electrostatic excitations in dust-contaminated dense plasmas, e.g., in metallic electronic devices, and also arguably in supernova environments, where charged dust defects may occur in the quantum plasma regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integro-differential modeling of ICRH wave propagation and damping at arbitrary cyclotron harmonics and wavelengths in tokamaks

Both at low and higher cyclotron harmonics, properly accounting for finite Larmor radius effects is crucial in many ion cyclotron resonance frequency heating scenarios creating high energy tails. The present paper discusses ongoing work to extend the 1D TOMCAT wave equation solver [D. Van Eester and R. Koch, Plasma Phys. Contr. Fusion 40 (1998) 1949] to arbitrary harmonics and arbitrary wavelengths. Rather than adopting the particle position, the guiding center position is used as the independent variable when writing down an expression for the dielectric response. Adopting a philosophy originally due to Kaufman [A.N. Kaufman, Phys. Fluids 15 (1972) 1063], the relevant dielectric response in the Galerkin formalism is written in a form where the electric field and the test function vector appear symmetrically, which yields a power balance equation that guarantees non-negative absorption for any wave type for Maxwellian plasmas. Moreover, this choice of independent variable yields intuitive expressions that can directly be linked to the corresponding expressions in the RF diffusion operator. It also guarantees that a positive definite power transfer from waves to particles is ensured for any of the wave modes in a plasma in which all populations have a Maxwellian distribution, as is expected from first principles. Rather than relying on a truncated Taylor series expansion of the dielectric response, an integro-differential approach that retains all finite Larmor radius effects [D. Van Eester and E. Lerche, Plasma Phys. Control. Fusion 55 (2013) 055008] is proposed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Plasma dark matter direct detection

Dark matter in spiral galaxies like the Milky Way may take the form of a dark plasma. Hidden sector dark matter charged under an unbroken U(1)' gauge interaction provides a simple and well defined particle physics model realising this possibility. The assumed U(1)' neutrality of the Universe then implies (at least) two oppositely charged dark matter components with self-interactions mediated via a massless 'dark photon' (the U(1)' gauge boson). In addition to nuclear recoils such dark matter can give rise to keV electron recoils in direct detection experiments. In this context, the detailed physical properties of the dark matter plasma interacting with the Earth is required. This is a complex system, which is here modelled as a fluid governed by the magnetohydrodynamic equations. These equations are numerically solved for some illustrative examples, and implications for direct detection experiments discussed. In particular, the analysis presented here leaves open the intriguing possibility that the DAMA annual modulation signal is due primarily to electron recoils (or even a combination of electron recoils and nuclear recoils). The importance of diurnal modulation (in addition to annual modulation) as a means of probing this kind of dark matter is also emphasised.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonlinear dust acoustic waves and shocks

We describe experiments on (1) nonlinear dust acoustic waves and (2) dust acoustic shocks performed in a direct current (DC) glow discharge dusty plasma. First, we describe experiments showing nonlinear dust acoustic waves characterized by waveforms of the dust density that are typically sharper in the wave crests and flatter in the wave troughs (compared to sinusoidal waves), indicating the development of wave harmonics. We discuss this behavior in terms of a second-order fluid theory for dust acoustic waves. Second, experimental observations of the propagation and steepening of large-amplitude dust acoustic waves into dust acoustic shock waves are presented. The observed shock wave evolution is compared with numerical calculations based on the Riemann solution of the fully nonlinear fluid equations for dust acoustic waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasi-static magnetic compression of field-reversed configuration plasma: amended scalings and limits from two-dimensional MHD equilibrium

In this work, several key scaling laws of the quasi-static magnetic compression of field reversed configuration (FRC) plasma (Spencer et al 1983 Phys. Fluids 26 1564) are amended from a series of two-dimensional FRC MHD equilibriums numerically obtained using the Grad–Shafranov equation solver NIMEQ. Based on the new scaling for the elongation and the magnetic fields at the separatrix and the wall, the empirically stable limits for the compression ratio, the fusion gain, and the neutron yield are evaluated, which may serve as a more accurate estimate for the upper ceiling of performance from the magnetic compression of FRC plasma as a potential fusion energy as well as neutron source devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On hydromagnetic wave interactions in collisionless, high- β plasmas

We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of ion temperature on ion-acoustic solitary waves in a magnetized plasma in presence of superthermal electrons

Obliquely propagating ion-acoustic soliatry waves are examined in a magnetized plasma composed of kappa distributed electrons and fluid ions with finite temperature. The Sagdeev potential approach is used to study the properties of finite amplitude solitary waves. Using a quasi-neutrality condition, it is possible to reduce the set of equations to a single equation (energy integral equation), which describes the evolution of ion-acoustic solitary waves in magnetized plasmas. The temperature of warm ions affects the speed, amplitude, width, and pulse duration of solitons. Both the critical and the upper Mach numbers are increased by an increase in the ion temperature. The ion-acoustic soliton amplitude increases with the increase in superthermality of electrons. For auroral plasma parameters, the model predicts the soliton speed, amplitude, width, and pulse duration, respectively, to be in the range of (28.7-31.8) km/s, (0.18-20.1) mV/m; (590-167) m, and (20.5-5.25) ms, which are in good agreement with Viking observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Relativistic Modeling Capabilities in Perseus Extended MHD Simulation Code for HED Plasmas

This project is submitted in response to DE-FOA-0001153, and will be carried out by Dr. Charles Seyler and Dr. Nathaniel Hamlin as post doc under the direction of Professor Seyler. The objectives are to enable PERSEUS to self-consistently model relativistic high-energy-density (HED) laboratory and astrophysical phenomena. Electrons and ions in a plasma are generally modeled either as fluids using a two-fluid code, like PERSEUS, or as particles using a Particle-In- Cell (PIC) code. Nonrelativistic PERSEUS solves the two-fluid equations, formulated in terms of a generalized Ohm’s law (GOL), so as to model about nine orders of magnitude in density variation using a local semi-implicit method. Relativistic PERSEUS preserves this structure and therefore retains these advantageous properties, enabling it to model a broader range of relativistic HED phenomena than previous codes. We have also overcome a major technical challenge associated with solving a relativistic system of equations relating conserved quantities (momentum, energy, etc.) to primitive variables (density, velocity, and pressure). We have thus far developed a relativistic version of PERSEUS and demonstrated that it recovers expected nonrelativistic results. We are currently using relativistic PERSEUS to model two HED phenomena, namely laser-plasma interactions and X-pinches, both of which show spectroscopic evidence of relativistic electrons. In an X-pinch, a section of wire carrying sufficiently large current forms a plasma column that pinches down to a point-like source, generating X-rays for use in radiographic imaging. Laser-plasma interactions have been simulated using PIC codes, and we use their results as a benchmark for comparison. However, neither PIC codes nor existing two-fluid schemes can simulate X-pinches, due to the broad range of densities which have thus far made such a simulation computationally prohibitive. In both laser-plasma interaction and Xpinch simulations, we have observed relativistic phenomena that PERSEUS could not previously model. This includes relativistic channeling of a laser into an over-dense deuterium gas (i.e. induced transparency), consistent with PIC simulation results, along with energetic electrons and ions following an X-pinch and likely accompanied by strong X-ray emission, consistent with experimental observations of energetic X-rays. This proposal includes our latest simulation results. We will improve computational efficiency with techniques that allow the use of a less refined numerical grid. This will enable extensions to three-dimensional geometries. We will improve modeling of laser-plasma interactions and X-pinches, including radiative transport processes, in order to improve comparisons with PIC simulations and experimental diagnostics. We will also model power loss through electrode surface plasmas in magnetically insulated transmission lines (MITLs). We will model relativistic astrophysical phenomena in both astrophysical and HED regimes, which could shed light on the validity of laboratory astrophysics for certain problems. These phenomena are currently modeled using PIC codes and existing relativistic two-fluid codes, neither of which can model the range of length scales of which relativistic PERSEUS is capable. There is strong potential for this project to have a broad impact on both the astrophysical and laboratory plasma communities. In the laboratory community, our work is applicable to those concerned with laser-plasma interactions such as laser-driven fusion, and those in the dense Z pinch community concerned with X-pinches and electrode surface plasmas. An improved understanding of the mechanisms of hard (i.e. high-energy) X-ray generation will likely make important contributions toward the use of X-pinches for radiography through the harnessing or mitigation of hard X-rays, depending on the application. The modeling of electrode surface plasmas, including vacuum electron flow, will improve our understanding of power loss to various loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evolution of rogue waves in dusty plasmas

The evolution of rogue waves associated with the dynamics of positively charged dust grains that interact with streaming electrons and ions is investigated. Using a perturbation method, the basic set of fluid equations is reduced to a nonlinear Schrödinger equation (NLSE). The rational solution of the NLSE is presented, which proposed as an effective tool for studying the rogue waves in Jupiter. It is found that the existence region of rogue waves depends on the dust-acoustic speed and the streaming densities of the ions and electrons. Furthermore, the supersonic rogue waves are much taller than the subsonic rogue waves by ∼25 times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗