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 865 records · Page 48

Pickup Ion Phase Space Distributions at Titan in a Three Dimensional Exosphere

The composition and structure of neutral exospheres imbedded in moving plasmas can be determined by measurements of the velocity distributions of their pickup ion progeny. In turn, the velocity distributions are dependent on the spatial structure of the neutral source gases. Since Titan's neutral exosphere extends into the Saturn's magnetosphere (or solar wind) and well above its ionopause, it serves as a good place to analyze such characteristics. They are analyzed using pickup ion measurements made by the Cassini Plasma Spectrometer (CAPS) at Titan [e.g., Hartle et al., 2006] and an ion kinetic model. An early version of the model [Hartle and Sittler, 2007] is an expression describing the phase space density of pickup ions, which is derived from the Vlasov equation with an ion source that explicitly accounts for the velocity and spatial variation of the exosphere source gases. The current version used here includes exosphere source gases in three dimensions. A fundamental parameter of the phase space densities is the ratio of the gyroradius to the neutral scale height alpha, = r(sub g)/H. Titan's exosphere structure yields pickup ions whose phase space distributions are beam-like when alpha >> 1 and fluid-like when alpha << 1. Downstream from the source peak, the light pickup ions, with alpha << 1, are easily observed because their phase space densities are almost uniform over the orbit phases. On the other hand, the phase space distributions of the heavier ions, with alpha >> 1, peak over narrow velocity and spatial ranges. This beam-like nature makes it considerably more difficult to observe heavy ions because their downstream positions and viewing directions are narrowly constrained. Examples of these extremes will be discussed.

Hartle, Richard↗

Time-dependent relaxed magnetohydrodynamics: Inclusion of cross helicity constraint using phase-space action

In this work, a phase-space version of the ideal magnetohydrodynamic (MHD) Lagrangian is derived from first principles and shown to give a relabeling transformation when a cross-helicity constraint is added in Hamilton's Action Principle. A new formulation of time-dependent “relaxed” magnetohydrodynamics is derived using microscopic conservation of mass and macroscopic constraints on total magnetic helicity, cross helicity, and entropy under variations of density, pressure, fluid velocity, and magnetic vector potential. This gives Euler–Lagrange equations consistent with previous work on both ideal and relaxed MHD equilibria with flow, but generalizes the relaxation concept from statics to dynamics. The application of the new dynamical formalism is illustrated for short-wavelength linear waves, and the interface connection conditions for Multiregion Relaxed MHD (MRxMHD) are derived. The issue of whether E+u×B=0 should be a constraint is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma expansion and evolution of density perturbations in the polar wind - Comparison of semikinetic and transport models

Comparisons are made between transport and semikinetic models in a study of the time evolution of plasma density perturbations in the polar wind. The situations modeled include plasma expansion into a low-density region and time evolution of localized density enhancements and cavities. The results show that the semikinetic model generally yields smoother profiles in density, drift velocity, and ion temperature than the transport model, principally because of ion velocity dispersion. While shocks frequently develop in the results of the transport model, they do not occur in the semikinetic results. In addition, in the semikinetic results, two ion streams, or double-humped distributions, frequently develop. In the transport model results the bulk parameters, at a given time, often have a one-to-one correspondence in the locations of their local minima or maxima. This is a consequence of the coupling of the fluid equations. There is, however, no such relationship among the moments produced by the semikinetic model where the local moment maxima and minima are often shifted in altitude. In general, incorporation of enhanced heat fluxes in the transport model leads to somewhat improved agreement with the semikinetic results.

Ho, C. W.↗

Magnetic-field generation and its effect on ablative Rayleigh–Taylor instability in diffusive ablation fronts

The effects of self-generated magnetic fields on the ablative Rayleigh–Taylor (RT) instability are investigated in the linear regime. The main governing parameters are the Froude number (Fr), which stands for the ratio between ablative convection and acceleration of the target, and the Mach number at the ablation front (Ma), assumed to be small (isobaricity). During the development of the RT instability, magnetic fields are generated due to misalignment between pressure and density gradients (Biermann-battery effect). They accumulate at the section of the ablation front where the Nernst and the plasma velocities cancel each other. The magnetic field modifies the dynamics of the instability through the Righi–Leduc term, which acts as a heat source in the energy equation. It is found that the B fields affect perturbations with short wavelengths up to the most unstable wave in the spectrum. Here, the B field plays a destabilizing role for moderate Froude numbers and becomes stabilizing for large Froude numbers. For plastic ablators, the Fr threshold is found to be Fr = 5.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hermes-3: Multi-component plasma simulations with BOUT++

A new open source tool for fluid simulation of multi-component plasmas is presented, based on a flexible software design that is applicable to scientific simulations in a wide range of fields. Hermes-3 is built on plasma simulation framework BOUT++, consolidating earlier SD1D and Hermes models into a single code that can be configured at run-time to solve plasma models in 1D, 2D or 3D, either for transport (steady-state) or turbulent (time-evolving) problems, with an arbitrary number of ion and neutral species. Here, we describe the improved numerical algorithms and software design that have been implemented in Hermes-3. To demonstrate the capabilities of this tool, applications relevant to the boundary of tokamak plasmas are presented: 1D simulations of diveror plasmas evolving equations for all charge states of neon and deuterium; 2D transport simulations of tokamak equilibria in single-null X-point geometry with plasma ion and neutral atom species; and simulations of the time-dependent propagation of plasma filaments (blobs). Hermes-3 is publicly available on Github under the GPL-3 open source license. The repository includes documentation and a suite of unit, integrated and convergence tests.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A computational model for exploring particle acceleration during reconnection in macroscale systems

A new computational model suitable for exploring the self-consistent production of energetic electrons during magnetic reconnection in macroscale systems is presented. The equations are based on the recent discovery that parallel electric fields are ineffective drivers of energetic particles during reconnection so that the kinetic scales which control the development of such fields can be ordered out of the equations. The resulting equations consist of a magnetohydrodynamic (MHD) backbone with the energetic component represented by macro-particles described by the guiding center equations. Crucially, the energetic component feeds back on the MHD equations so that the total energy of the MHD fluid and the energetic particles is conserved. The equations correctly describe the firehose instability, whose dynamics plays a key role in throttling reconnection and in controlling the spectra of energetic particles. The results of early tests of the model, including the propagation of Alfvén waves in a system with pressure anisotropy and the growth of firehose modes, establish that the basic algorithm is stable and produces reliable physics results in preparation for further benchmarking with particle-in-cell models of reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of toroidally localized harmonics of global Alfvén eigenmodes

Abstract The spherical tokamak NSTX and the upgraded, higher toroidal field, NSTX-U, routinely create plasmas with large un-thermalized populations of super-Alfvénic fast ions. This population of non-thermal energetic ions excites a broad spectrum of Alfvénic waves, from the lower frequency (≈30 kHz–≈200 kHz) toroidal Alfvén eigenmodes (Cheng and Chance 1986 Phys. Fluids 29 3695), to the higher frequency (400 kHz–3 MHz) global Alfvén eigenmodes (GAEs) and compressional Alfvén eigenmodes (Goedbloed 1975 Phys. Fluids 18 1258). In this paper we present evidence that the GAE non-linearly excite modes, presumably GAE, at frequencies consistent with non-linear or 3-wave coupling. The observation of the excitation of 2 nd harmonic GAE through the intrinsic non-linearity of Global Alfvén modes demonstrates that the non-linear terms can act as an exciter-antenna inside the plasma, broadcasting at harmonics of the mode frequency and with concomitantly shorter wavelengths. As with experiments using external antenna to excite otherwise weakly stable Alfvénic modes, this data can provide information about the stability of modes at harmonics of the GAE. It may also provide information on the nature of the non-linearities in the wave dispersion equation or a direct measure of the mode amplitude (Smith et al 2006 Phys. Plasmas 13 042504). We also report that the short wavelength 2 nd harmonic GAE can be strongly toroidally localized. The observation that the shorter wavelength ( n ≈ 20) 2 nd harmonic modes can be toroidally localized potentially has implications for the impact of fast-ion driven instabilities on fast-ion confinement in ITER and future fusion reactors (Gorelenkov et al 2014 Nucl. Fusion 54 125001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linearized Boltzmann-Langevin model for heavy quark transport in hot and dense QCD matter

In relativistic heavy-ion collisions, the production of heavy quarks at large transverse momenta is strongly suppressed compared to proton-proton collisions. In addition, an unexpectedly large azimuthal anisotropy was observed for the emission of charmed hadrons in noncentral collisions. Both observations pose challenges to the theoretical understanding of the coupling between heavy quarks and the quark-gluon plasma produced in these collisions. Transport models for the evolution of heavy quarks in a QCD medium offer the opportunity to study these effects; two of the most successful approaches are based on the linearized Boltzmann transport equation and the Langevin equation. In this work, we develop a hybrid transport model that combines the strengths of both of these approaches: Heavy quarks scatter with medium partons using matrix-elements calculated in perturbative QCD, while between these discrete hard scatterings they evolve using a Langevin equation with empirical transport coefficients to capture the nonperturbative soft part of the interaction. With the hybrid transport model coupled to a state-of-the-art event-by-event bulk evolution model based on 2+1D relativistic viscous fluid dynamics, we study the azimuthal anisotropy and nuclear modification factor of heavy quarks in Pb+Pb collisions at √s = 5.02 TeV. The parameters related to heavy-flavor transport are calibrated using a Bayesian analysis comparing them to available D -meson and B -meson data at the Large Hadron Collider. Lastly, using the calibrated model, we study the implications on heavy-flavor transport properties and predict observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement of the sound speed in dense fluid deuterium along the cryogenic liquid Hugoniot

Many experimental studies, spanning several decades of research and using various dynamic compression schemes, have been conducted to investigate cryogenic liquid deuterium under strong shock compression. The consensus emerging from these studies is that of a progressive dissociation of the D 2 molecules into an electrically conducting, atomic plasma, when subjected to shock pressures exceeding ~50 GPa. While state-of-the-art numerical simulations based on density-functional-theory or quantum Monte-Carlo techniques capture this behavior quite well, subtle differences subsist between these simulations and the available experimental data regarding the pressure-density compressibility. Here, leveraging a recently developed analysis method for high-resolution Doppler interferometric velocity data, we present Eulerian sound speed measurements in compressed deuterium to shock pressures between 50 and 200 GPa. These results, extracted from laser-driven shockwave experiments, are found to agree with several of the most accurate equation of state models for deuterium at those conditions up to ~150 GPa. However, the data indicate that these models fail to reproduce the experimentally observed sound speed at higher pressures, approaching 200 GPa. In particular, we unveil a discrepancy between the experimental results and the equation of state model that is most commonly used in inertial confinement fusion at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-Linear Dynamics and Emergence in Laboratory Fusion Plasmas

Turbulent behaviour of laboratory fusion plasma system is modelled using extended Hasegawa-Wakatani equations. The model is solved numerically using finite difference techniques. We discuss non-linear effects in such a system in the presence of the micro-instabilities, specifically a drift wave instability. We explore particle dynamics in different range of parameters and show that the transport changes from diffusive to non-diffusive when large directional flows are developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Excitation of flow instabilities due to nonlinear scale invariance

A novel route to instabilities and turbulence in fluid and plasma flows is presented in kinetic Vlasov-Maxwell model. New kind of flow instabilities is shown to arise due to the availability of new kinetic energy sources which are absent in conventional treatments. The present approach is based on a scale invariant nonlinear analytic formalism developed to address irregular motions on a chaotic attractor or in turbulence in a more coherent manner. We have studied two specific applications of this turbulence generating mechanism. The warm plasma Langmuir wave dispersion relation is shown to become unstable in the presence of these multifractal measures. In the second application, these multifractal measures are shown to induce naturally non-Gaussian, i.e., a stretched, Gaussian distribution and anomalous transport for tracer particles from the turbulent advection-diffusion transport equation in a Vlasov plasma flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiphysics simulation of corona discharge induced ionic wind

Ionic wind devices or electrostatic fluid accelerators are becoming of increasing interest as tools for thermal management, in particular for semiconductor devices. In this work, we present a numerical model for predicting the performance of such devices; its main benefit is the ability to accurately predict the amount of charge injected from the corona electrode. Our multiphysics numerical model consists of a highly nonlinear, strongly coupled set of partial differential equations including the Navier-Stokes equations for fluid flow, Poisson's equation for electrostatic potential, charge continuity, and heat transfer equations. To solve this system we employ a staggered solution algorithm that generalizes Gummel's algorithm for charge transport in semiconductors. Predictions of our simulations are verified and validated by comparison with experimental measurements of integral physical quantities, which are shown to closely match.

36 MATERIALS SCIENCE↗

Direct construction of optimized stellarator shapes. Part 1. Theory in cylindrical coordinates

The confinement of the guiding-centre trajectories in a stellarator is determined by the variation of the magnetic field strength $B$ in Boozer coordinates $(r,\unicode[STIX]{x1D703},\unicode[STIX]{x1D711})$ , but $B(r,\unicode[STIX]{x1D703},\unicode[STIX]{x1D711})$ depends on the flux surface shape in a complicated way. Here we derive equations relating $B(r,\unicode[STIX]{x1D703},\unicode[STIX]{x1D711})$ in Boozer coordinates and the rotational transform to the shape of flux surfaces in cylindrical coordinates, using an expansion in distance from the magnetic axis. A related expansion was done by Garren and Boozer ( Phys. Fluids B, vol. 3, 1991 a , 2805) based on the Frenet–Serret frame, which can be discontinuous anywhere the magnetic axis is straight, a situation that occurs in the interesting case of omnigenity with poloidally closed $B$ contours. Our calculation in contrast does not use the Frenet–Serret frame. The transformation between the Garren–Boozer approach and cylindrical coordinates is derived, and the two approaches are shown to be equivalent if the axis curvature does not vanish. The expressions derived here help enable optimized plasma shapes to be constructed that can be provided as input to VMEC and other stellarator codes, or to generate initial configurations for conventional stellarator optimization.

Physics↗

Multi-fluid Modeling of Magnetosonic Wave Propagation in the Solar Chromosphere: Effects of Impact Ionization and Radiative Recombination

In order to study chromospheric magnetosonic wave propagation including, for the first time, the effects of ion–neutral interactions in the partially ionized solar chromosphere, we have developed a new multi-fluid computational model accounting for ionization and recombination reactions in gravitationally stratified magnetized collisional media. The two-fluid model used in our 2D numerical simulations treats neutrals as a separate fluid and considers charged species (electrons and ions) within the resistive MHD approach with Coulomb collisions and anisotropic heat flux determined by Braginskiis transport coefficients. The electromagnetic fields are evolved according to the full Maxwell equations and the solenoidality of the magnetic field is enforced with a hyperbolic divergence-cleaning scheme. The initial density and temperature profiles are similar to VAL III chromospheric model in which dynamical, thermal, and chemical equilibrium are considered to ensure comparison to existing MHD models and avoid artificial numerical heating. In this initial setup we include simple homogeneous flux tube magnetic field configuration and an external photospheric velocity driver to simulate the propagation of MHD waves in the partially ionized reactive chromosphere. In particular, we investigate the loss of chemical equilibrium and the plasma heating related to the steepening of fast magnetosonic wave fronts in the gravitationally stratified medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Self-consistent interaction of linear gravitational and electromagnetic waves in non-magnetized plasma

Here, this paper explores the hybridization of linear metric perturbations with linear electromagnetic (EM) perturbations in non-magnetized plasma for a general background metric. The local wave properties are derived from first principles for inhomogeneous plasma, without assuming any symmetries of the background metric. First, we derive the effective (“oscillation-center”) Hamiltonian that governs the average dynamics of plasma particles in a prescribed quasimonochromatic wave that involves metric perturbations and EM fields simultaneously. Then, using this Hamiltonian, we derive the backreaction of plasma particles on the wave itself and obtain gauge-invariant equations that describe the resulting self-consistent gravito-electromagnetic (GEM) waves in a plasma. The transverse tensor modes of gravitational waves are found to have no interaction with the plasma and the EM modes in the geometrical-optics limit. However, for $longitudinal$ GEM modes with large values of the refraction index, the interplay between gravitational and EM interactions in plasma can have a strong effect. In particular, the dispersion relation of the Jeans mode is significantly affected by electrostatic interactions. As a spin-off, our calculation also provides an alternative resolution of the so-called Jeans swindle.

79 ASTRONOMY AND ASTROPHYSICS↗

Transient Relativistic Fluid Dynamics in a General Hydrodynamic Frame

In this work, we propose a new theory of second-order viscous relativistic hydrodynamics which does not impose any frame conditions on the choice of the hydrodynamic variables. It differs from Mueller-Israel-Stewart theory by including additional transient degrees of freedom, and its first-order truncation reduces to Bemfica-Disconzi-Noronha-Kovtun theory. Conditions for causality and stability are explicitly given in the conformal regime. As an illustrative example, we consider Bjorken flow solutions to our equations and identify variables which make a hydrodynamic attractor manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The limits of the Bohm criterion in collisional plasmas

The sheath formation within a low-pressure collisional plasma is analysed by means of a two-fluid model. The Bohm criterion takes into account the effects of the electric field and the inertia of the ions. Numerical results yield that these effects contribute to the space charge formation, only, if the collisionality is lower than a relatively small threshold. It follows that a lower and an upper limit of the drift speed of the ions exist where the effects treated by Bohm can form a sheath. This interval becomes narrower as the collisionality increases and vanishes at the mentioned threshold. Above the threshold, the sheath is mainly created by collisions and the ionisation. Under these conditions, the sheath formation cannot be described by means of Bohm like criteria. In a few references, a so-called upper limit of the Bohm criterion is stated for collisional plasmas where the momentum equation of the ions is taken into account, only. However, the present paper shows that this limit results in an unrealistically steep increase of the space charge density towards the wall, and, therefore, it yields no useful limit of the Bohm velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗