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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Thermalization & hydrodynamics in Bjorken & Gubser flows

The dynamical scaling behavior of hydrodynamic and non-hydrodynamic moments of the distribution function is studied using third-order Chapman-Enskog hydrodynamics and anisotropic hydrodynamics for systems undergoing Bjorken and Gubser expansions, where exact solutions of the Boltzmann equation in Relaxation Time Approximation (RTA) are available for comparison. While dimensionless quantities like normalized shear, pressure anisotropy and normalized entropy show at late times universal scaling relations with small (large) Knudsen number for Bjorken (Gubser) flows, dimensionful quantities like the entropy per unit rapidity do not. Although the two hydrodynamic approximation schemes describe the exact attractors for normalized shear with high accuracy, their description for the normalized entropy is less precise. We attribute this to non-negligible couplings to non-hydrodynamic modes in the entropy evolution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Ion kinetic effects on the ignition and burn of inertial confinement fusion targets: A multi-scale approach

In this article, we study the hydrodynamics and burn of the thermonuclear fuel in inertial confinement fusion pellets at the ion kinetic level. The analysis is based on a two-velocity-scale Vlasov-Fokker-Planck kinetic model that is specially tailored to treat fusion products (suprathermal α-particles) in a self-consistent manner with the thermal bulk. The model assumes spherical symmetry in configuration space and axial symmetry in velocity space around the mean flow velocity. A typical hot-spot ignition design is considered. Compared with fluid simulations where a multi-group diffusion scheme is applied to model α transport, the full ion-kinetic approach reveals significant non-local effects on the transport of energetic α-particles. This has a direct impact on hydrodynamic spatial profiles during combustion: the hot spot reactivity is reduced, while the inner dense fuel layers are pre-heated by the escaping α-suprathermal particles, which are transported farther out of the hot spot. We show how the kinetic transport enhancement of fusion products leads to a significant reduction of the fusion yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Equation-of-state, sound speed, and reshock of shock-compressed fluid carbon dioxide

Mechanical equation-of-state data of initially liquid and solid CO2 shock-compressed to terapascal conditions are reported. Diamond-sapphire anvil cells were used to vary the initial density and state of CO2 samples that were then further compressed with laser-driven shock waves, resulting in a data set from which precise derivative quantities, including Grüneisen parameter and sound speed, are determined. Reshock states are measured to 800 GPa and map the same pressure-density conditions as the single shock using different thermodynamic paths. The compressibility data reported here do not support current density-functional-theory calculations, but are better represented by tabular equation-of-state models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2001

Properties and Radial Trends of Coronal Mass Ejecta and Their Associated Shocks Observed by Ulysses in the Ecliptic Plane

In this paper, magnetic and plasma measurements are used to analyze 17 interplanetary coronal mass ejections (CMEs) identified by Ulysses during its in-ecliptic passage to Jupiter. We focus on the expansion characteristics of these CMEs (as inferred from the time rate of change of the velocity profiles through the CMEs) and the properties of 14 forward shocks unambiguously associated with these CMEs. We highlight radial trends from 1 to 5.4 AU. Our results indicate that the CMEs are generally expanding at all heliocentric distances. With regard to the shocks preceding these ejecta, we note the following: (1) There is a clear tendency for the shock speed (in the upstream frame of reference) to decrease with increasing heliocentric distance as the CMEs transfer momentum to the ambient solar wind and slow down; (2) 86% of the shock fronts are oriented in the ecliptic plane such that their normals point westward (i.e., in the direction of planetary motion about the Sun), (3) 86% of the shocks are propagating toward the heliographic equator; and (4) no clear trend was found in the strength of the shocks versus heliocentric distance. These results are interpreted using simple dynamical arguments and are supported by fluid and magnetohydrodynamic (MHD) simulations.

Riley, Pete↗
DOE OSTI2025

Thermodynamic properties of “near-perfect” gas

When applied thermodynamics requires simple representations of thermodynamic quantities, one often finds quantities such as pressure and the thermal energy density expressed as a monomial linear product of density and temperature power laws. This is a simple generalization of the product of density and temperature used for the perfect gas that admits a broader range of thermodynamic behavior into analytic fluid calculations or provides an analytic form that is readily fit to tabulated equation-of-state data for an arbitrary material. This paper reviews the thermodynamic properties of this generalized perfect-gas model, treating it as a class of “near-perfect-gas” models that are defined and unified here, for the first time, through a shared Helmholtz free energy. This generalization from perfect to near-perfect-gas models preserves some perfect-gas properties (e.g., constant adiabatic exponents) but not others (e.g., constant specific heats) and is constrained by thermodynamic consistency. It is important to be aware of the properties of this class of near-perfect-gas models when using them as surrogates for real nonperfect materials that may have properties that cannot be captured by near-perfect-gas models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A crude model to study radio frequency induced density modification close to launchers

The interplay between radio frequency (RF) waves and the density is discussed by adopting the general framework of a 2-time-scale multi-fluid treatment, allowing to separate the dynamics on the RF time scale from that on the time scale on which macroscopic density and flows vary as a result of the presence of electromagnetic and/or electrostatic fields. The focus is on regions close to launchers where charge neutrality is incomplete and waves are commonly evanescent. The fast time scale dynamics influences the slow time scale behavior via quasilinear terms (the Ponderomotive force for the case of the equation of motion). Electrons and ions are treated on the same footing. Also, both fast and slow waves are retained in the wave description. Although this work is meant as a subtopic of a large study—the wave induced “convective cell” physics at hand is of a 2- or 3-dimensional nature while this paper limits itself to a single dimension—a few tentative examples are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Transfer Request from Rutgers University to University of Illinois Urbana Champaign

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons will be exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges must be developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients will be calculated across beam energies. New techniques will be developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project will provide essential guidance to the upcoming Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Dynamical Aspects of the Quark Gluon Plasma

Shortly after the Big Bang the entire universe was filled with a nearly perfect fluid known as the Quark Gluon Plasma. Relativistic heavy ion collisions can now reproduce this fluid in the laboratory where the Quark Gluon Plasma exhibits a rapid but smooth cross-over phase transition into hadrons at vanishing net-baryon densities. Recent experiments plan to explore finite baryon densities where a critical point is expected. If found, this would mark the first discovery of a critical point in a relativistic system described by a fundamental theory of nature, which would have far-reaching consequences for high-energy nuclear physics and nuclear astrophysics (such as in neutron star mergers). Characteristic temperatures of equilibrium (e.g. the inflection point of the entropy density) and transport coefficients (e.g. minimum of the shear viscosity over entropy density) vary widely at a cross-over phase transition but converge at a critical point, and extracting the behavior of these characteristic temperatures is a major focal point of this research. Specifically, the interplay between strange and light hadrons is exploited to study the flavor hierarchy in the cross-over region. To investigate this, a viscous relativistic hydrodynamics framework with two conserved charges is being developed into a new open-source code along with initial conditions that contain baryon number and strangeness. Flow observables sensitive to the equation of state and transport coefficients are calculated across beam energies. New techniques are being developed to study this flavor hierarchy from first principles and to extract the characteristic temperatures from experimental data. Through this new dynamical framework, this project provides essential guidance to the Beam Energy Scan II runs at Relativistic Heavy-Ion Collider and the future Facility for Antiproton and Ion Research facility in the search for the Quantum Chromodynamic critical point and subsequent investigation of the baryon-rich Quark Gluon Plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Microphysics Studies of Hohlraum Dynamics (LDRD 17-ERD-060 Final Report)

The goal of this project was demonstrate the feasibility and utility of modeling the effects of kinetic physics in gas-filled hohlraums using high-fidelity simulation tools. We have successfully achieved this goal. In our work, we have found that hybrid particle-in-cell simulations (Thoma et al. 2017, Higginson et al. 2019a) utilizing quasi-neutral, fluid electrons and kinetic ions are powerful tools to understand the kinetic physics of interpenetration and mix relevant to low-gas-fill (near-vacuum) hohlraums. We have benchmarked the code, Chicago, against OMEGA laser experiments designed to mimic conditions found in hohlraums; the code was able to reproduce data from x-ray imaging and Thomson Spectra resolved over both time and space. The experiments highlighted kinetic-ion effects such increased penetration depth that were only able to be resolved with a kinetic model. We have run preliminary simulations with Chicago of full-scale hohlraum experiments that include laser raytracing, equation-of-state and ionization, radiation transport, and accurate experimental geometries. Continuation of this work will predictive hohlraum models at low-gas-fill and identify the better regimes of hohlraum design with the potential for improved coupling and capsule symmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

As a Matter of State: The Role of Thermodynamics in Magnetohydrodynamic Turbulence

Turbulence simulations play a key role in advancing the general understanding of the physical properties of turbulence and in interpreting astrophysical observations of turbulent plasmas. For the sake of simplicity, however, turbulence simulations are often conducted in the isothermal limit. Given that the majority of astrophysical systems are not governed by isothermal dynamics, we aim to quantify the impact of thermodynamics on the physics of turbulence, through varying adiabatic index, γ, combined with a range of optically thin cooling functions. Here, we present a suite of ideal magnetohydrodynamics simulations of thermally balanced stationary turbulence in the subsonic, super-Alfvénic, high ${\beta }_{{\rm{p}}}$ (ratio of thermal to magnetic pressure) regime, where turbulent dissipation is balanced by two idealized cooling functions (approximating linear cooling and free–free emission) and examine the impact of the equation of state by considering cases that correspond to isothermal, monatomic, and diatomic gases. We find a strong anticorrelation between thermal and magnetic pressure independent of thermodynamics, whereas the strong anticorrelation between density and magnetic field found in the isothermal case weakens with increasing γ. Similarly, the linear relation between variations in density and thermal pressure with sonic Mach number becomes steeper with increasing γ. This suggests that there exists a degeneracy in these relations with respect to thermodynamics and Mach number in this regime, which is dominated by slow magnetosonic modes. These results have implications for attempts to infer (e.g.,) Mach numbers from (e.g.,) Faraday rotation measurements, without additional information regarding the thermodynamics of the plasma. However, our results suggest that this degeneracy can be broken by utilizing higher-order moments of observable distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

Symmetries of resistive two-fluid magnetohydrodynamics under reversals of toroidal field, current, and rotation

The discrete symmetries of time-dependent, non-axisymmetric, two-fluid magnetohydrodynamics (MHD) equations are considered. Solutions of the time-independent, axisymmetric, ideal-MHD equations remain solutions under reversals of the toroidal field, current density, or rotation. Introducing non-axisymmetry, resistivity, and two-fluid effects into the equations each break different symmetries. Symmetry groups for solutions possessing up-down spatial symmetry, and the effect of flipping the magnetic geometry across the horizontal midplane, are also considered. It is shown that poloidal velocity may or may not be reversed under a simultaneous reversal of the toroidal field and a vertical reflection of the magnetic geometry, depending on the symmetry of the velocity. Because the symmetry groups of ideal, resistive, and two-fluid MHD are distinct, it may be possible to ascertain the dominant physical mechanism of various phenomena through empirical observations of symmetry-breaking. These results, which hold for nonlinear solutions in arbitrary geometry, should also be of use in testing numerical codes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Electromagnetic effect on geodesic acoustic mode with adiabatic electrons

The geodesic acoustic mode (GAM) is analytically investigated by taking into account the finite-orbit-width (FOW) resonance effect to the second order and the finite β effect. The general dispersion relation is derived from the gyro-kinetic equations in the presence of nonzero δA ∥ , the parallel component of the perturbed magnetic vector potential. Transparent and concise expressions for the GAM frequency and Landau damping rate in the presence of the second order FOW effect and finite β effect are first presented. It is clearly shown that the m = ±2 harmonics dominant δA ∥ and the kinetic expression of δA ∥ have the same form as the fluid one. For the real frequency, the electromagnetic effect introduces a term on the order of q 2 β, which is comparable to the second order electrostatic terms, namely, the terms introduced by the second order FOW resonance effect. While for the collisionless damping rate, δA ∥ does not directly introduce β–dependent terms, but affects the damping rate via modifying the real frequency. Besides, our analytical result shows good agreement with the numerical examinations.

58 GEOSCIENCES↗
DOE OSTI2018

An adaptive, implicit, conservative, 1D-2V multi-species Vlasov–Fokker–Planck multi-scale solver in planar geometry

We consider a 1D-2V Vlasov–Fokker–Planck multi-species ionic description coupled to fluid electrons. We address temporal stiffness with implicit time stepping, suitably preconditioned. To address temperature disparity in time and space, we extend the conservative adaptive velocity-space discretization scheme proposed in [Taitano et al., J. Comput. Phys., 318, 391–420, (2016)] to a spatially inhomogeneous system. Here in this approach, we normalize the velocity-space coordinate to a temporally and spatially varying local characteristic speed per species. We explicitly consider the resulting inertial terms in the Vlasov equation, and derive a discrete formulation that conserves mass, momentum, and energy up to a prescribed nonlinear tolerance upon convergence. Our conservation strategy employs nonlinear constraints to enforce these properties discretely for both the Vlasov operator and the Fokker–Planck collision operator. Numerical examples of varying degrees of complexity, including shock-wave propagation, demonstrate the favorable efficiency and accuracy properties of the scheme.

1D2V↗
DOE OSTI2021

Cross-Code verification and sensitivity analysis to effectively model the electrothermal instability

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. Here, a sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2011

Modelling of Ion Cyclotron Wall Conditioning plasmas

Ion Cyclotron Wall Conditioning (ICWC) is envisioned in ITER to clean the wall from impurities, to control the wall isotopic ratio and the hydrogen recycling in the presence of the toroidal magnetic field. Various experiments and modelling are advancing to consolidate this technique. In this contribution the modeling of ICWC is presented, which can be divided in two parts: plasma description and plasma wall interaction. Firstly a 0D plasma model, based on a set of energy and particle balance equations for Maxwellian Hydrogen and Helium species, is presented. The model takes into account elementary collision processes, coupled RF power, particle confinement, wall recycling, and active gas injection and pumping. The RF plasma production process is based mainly on electron collisional ionization. The dependency of the plasma parameters, the Hydrogen and Helium partial pressures and neutral or ionic fluxes on pressure and RF power are quantitatively in good agreement with those obtained experimentally on TORE SUPRA. Secondly an extension of the 0D model including the description of the wall interaction is presented and compared to TORE SUPRA multi-pulse ICWC discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

An Equation of State for Liquid Pb83Li17

Thermal hydraulic and accident analysis codes such as RELAP5-3D and MELCOR for Fusion rely on an equation of state to specify all the thermodynamic properties of fusion-relevant working fluids such as PbLi. The existing liquid metal fluid properties in both RELAP5-3D and MELCOR are based on a five parameter "soft sphere" equation of state for which parameter sets that approximately reproduce experiment data were available in the literature. The PbLi parameters were based on a linear, mass-weighted average of the parameters for pure Pb and pure Li, which does not always result in adequate agreement between the equation and experimental data for PbLi, much of which was not available at the time the original equations were implemented. In order to address this shortcoming, we describe here a modification and non-linear least squares fitting of the parameters to closely match available experiment data for the density, sound speed, and specific heat capacity of liquid PbLi over a range of temperatures. The new equation has been implemented in both RELAP5-3D and MELCOR, and we make comparisons with both codes to experiment data to verify the implementation. Transport properties such as the thermal conductivity and viscosity, though not derived from the equation of state, have also been updated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

BHR equations re-derived with immiscible particle effects

Compressible and variable density turbulent flows with dispersed phase effects are found in many applications ranging from combustion to cloud formation. These types of flows are among the most challenging to simulate. While the exact equations governing a system of particles and fluid are known, computational resources limit the scale and detail that can be simulated in this type of problem. Therefore, a common method is to simulate averaged versions of the flow equations, which still capture salient physics and is relatively less computationally expensive. Besnard et al. [1] developed such a model for variable density miscible turbulence, where ensemble-averaging was applied to the flow equations to yield a set of filtered equations. Besnard et al. further derived transport equations for the Reynolds stresses, the turbulent mass flux, and the density-specific volume covariance, to help close the filtered momentum and continuity equations. We re-derive the exact BHR closure equations which include integral terms owing to immiscible effects. Physical interpretations of the additional terms are proposed along with simple models. The goal of this work is to extend the BHR model to allow for the simulation of turbulent flows where an immiscible dispersed phase is non-trivially coupled with the carrier phase.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2016

Molecular shear heating and vortex dynamics in thermostatted two dimensional Yukawa liquids

It is well known that two-dimensional macroscale shear flows are susceptible to instabilities leading to macroscale vortical structures. The linear and nonlinear fate of such a macroscale flow in a strongly coupled medium is a fundamental problem. A popular example of a strongly coupled medium is a dusty plasma, often modelled as a Yukawa liquid. Recently, laboratory experiments and molecular dynamics (MD) studies of shear flows in strongly coupled Yukawa liquids indicated the occurrence of strong molecular shear heating, which is found to reduce the coupling strength exponentially leading to the destruction of macroscale vorticity. To understand the vortex dynamics of strongly coupled molecular fluids undergoing macroscale shear flows and molecular shear heating, MD simulation has been performed, which allows the macroscopic vortex dynamics to evolve, while at the same time “removes” the microscopically generated heat without using the velocity degrees of freedom. We demonstrate that by using a configurational thermostat in a novel way, the microscale heat generated by shear flow can be thermostatted out efficiently without compromising the large scale vortex dynamics. In the present work, using MD simulations, a comparative study of shear flow evolution in Yukawa liquids in the presence and absence of molecular or microscopic heating is presented for a prototype shear flow, namely, Kolmogorov flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗