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 937 records · Page 52

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↗

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↗

Isotropic model of fractional transport in two-dimensional bounded domains

A two-dimensional fractional Laplacian operator is derived and used to model nonlocal, nondiffusive transport. This integro-differential operator appears in the long-wavelength, fluid description of quantities undergoing non-Brownian random walks without characteristic length scale. To study bounded domains, a mask function is introduced that modifies the kernel in the fractional Laplacian and removes singularities at the boundary. Green's function solutions to the fractional diffusion equation are presented for the unbounded domain and compared to the one-dimensional Cartesian approximations. A time-implicit numerical integration scheme is presented to study fractional diffusion in a circular disk with azimuthal symmetry. Numerical studies of steady-state reveal temperature profiles in which the heat flux and temperature gradient are in the same direction, i.e., uphill transport. The response to off-axis heating, scaling of confinement time with system size, and propagation of cold pulses are investigated.

ADVECTION-DISPERSION EQUATIONS↗

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗

Hugoniot and release measurements in diamond shocked up to 26 Mbar [Hugoniot and release measurements in diamond shocked up to 25 Mbar]

The equation of state (EOS) of carbon in its high-pressure solid and liquid phases is of interest to planetary astrophysics and inertial confinement fusion. Of particular interest are the high-pressure shock and release responses of diamond as these provide rigorous constraints on important paths through the EOS. This paper presents experimental Hugoniot and release data for both single-crystal diamond (SCD) and nanocrystalline diamond (NCD), which is comprised of nanometer-scale diamond grains and is ~5% less dense than SCD. We find that NCD has a stiffer Hugoniot than SCD that can be attributed to porosity. As a result, a Grüneisen parameter of ~1 was derived from the data, which suggests increased coordination in the high-pressure fluid carbon compared to ambient diamond.

36 MATERIALS SCIENCE↗

Chiral hydrodynamics in strong external magnetic fields

We construct the general hydrodynamic description of (3+1)-dimensional chiral charged (quantum) fluids subject to a strong external magnetic field with effective field theory methods. We determine the constitutive equations for the energy-momentum tensor and the axial charge current, in part from a generating functional. Furthermore, we derive the Kubo formulas which relate two-point functions of the energy-momentum tensor and charge current to 27 transport coefficients: 8 independent thermodynamic, 4 independent non-dissipative hydrodynamic, and 10 independent dissipative hydrodynamic transport coefficients. Five Onsager relations render 5 more transport coefficients dependent. We uncover four novel transport effects, which are encoded in what we call the shear-induced conductivity, the two expansion-induced longitudinal conductivities and the shear-induced Hall conductivity. Remarkably, the shear-induced Hall conductivity constitutes a novel non-dissipative transport effect. As a demonstration, we compute all transport coefficients explicitly in a strongly coupled quantum fluid via holography.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Hybrid (Vlasov-Fluid) simulation of ion-acoustic solitons chain formation including trapped electrons

Disintegration of a Gaussian profile into ion-acoustic solitons in the presence of trapped electrons [H. Hakimi Pajouh and H. Abbasi, Phys. Plasmas 15, 082105 (2008)] is revisited. Through a hybrid (Vlasov-Fluid) model, the restrictions associated with the simple modified Korteweg de-Vries (mKdV) model are studied. For instance, the lack of vital information in the phase space associated with the evolution of electron velocity distribution, the perturbative nature of mKdV model which limits it to the weak nonlinear cases, and the special spatio-temporal scaling based on which the mKdV is derived. Remarkable differences between the results of the two models lead us to conclude that the mKdV model can only monitor the general aspects of the dynamics, and the precise picture including the correct spatio-temporal scales and the properties of solitons should be studied within the framework of hybrid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Elliptic flow in Pb+Pb collisions at {radical}(s{sub NN}) = 2.76 TeV: Hybrid model assessment of the first data

We analyze the elliptic flow parameter v{sub 2} in Pb+Pb collisions at {radical}(s{sub NN}) = 2.76 TeV and in Au+Au collisions at {radical}(s{sub NN})=200 GeV using a hybrid model in which the evolution of the quark gluon plasma is described by ideal hydrodynamics with a state-of-the-art lattice QCD equation of state, and the subsequent hadronic stage by a hadron cascade model. For initial conditions, we employ Monte Carlo versions of the Glauber and the Kharzeev-Levin-Nardi models and compare results with each other. We demonstrate that the differential elliptic flow v{sub 2}(p{sub T}) hardly changes when the collision energy increases, whereas the integrated v{sub 2} increases due to the enhancement of mean transverse momentum. The amount of increase of both v{sub 2} and mean p{sub T} depends significantly on the model of initialization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Shock-driven hydrodynamic instability of a sinusoidally perturbed, high-Atwood number, oblique interface

A shock incident on an interface will cause any initial perturbations on that interface to grow. When the shock front is parallel to the interface, the perturbations grow due to the Richtmyer-Meshkov (RM) process. When there is some tilt between the shock front and the interface, shear flow will result across the postshock interface. Recent experiments on the OMEGA EP laser have studied the hydrodynamic instability growth which results from a supported shock interacting with a sinusoidally perturbed, oblique interface. The observed instability growth was dominated by Richtmyer-Meshkov at early times but became Kelvin-Helmholtz (KH)-like at late times. Previously, this instability growth was described using an analytic model for the deposition of baroclinic vorticity on the interface by a shock combined with a discrete vortex model. Here, we utilize the same baroclinic vorticity deposition model in conjunction with a desingularized, periodic Birkhoff-Rott equation to model instability evolution. The Birkhoff-Rott equation takes into account the vorticity distribution along the interface, whereas the discrete vortex model assumed that all vorticity over each wavelength of the perturbation is confined to a point. We compare the new model to xRAGE simulations and experiments. In conclusion, the model is found to overpredict both the instability growth and shear across the interface by about a factor of two, but correctly predicts that the growth is RM-like at early times and KH-like at late times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Δ, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. In this work, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Δ. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Δ variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfvén wave resonances as rotation slows down. These are important features that can alter the non-axisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multigrid treatment of implicit continuum diffusion

Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A generalized two-fluid picture of non-driven collisionless reconnection and its relation to whistler waves

A generalized, intuitive two-fluid picture of 2D non-driven collisionless magnetic reconnection is described using results from a full-3D numerical simulation. The relevant two-fluid equations simplify to the condition that the flux associated with canonical circulation Q=m e ∇×u e +q e B is perfectly frozen into the electron fluid. In the reconnection geometry, flux tubes defined by Q are convected with the central electron current, effectively stretching the tubes and increasing the magnitude of Q exponentially. This, coupled with the fact that Q is a sum of two quantities, explains how the magnetic fields in the reconnection region reconnect and give rise to strong electron acceleration. The Q motion provides an interpretation for other phenomena as well, such as spiked central electron current filaments. The simulated reconnection rate was found to agree with a previous analytical calculation having the same geometry. Energy analysis shows that the magnetic energy is converted and propagated mainly in the form of the Poynting flux, and helicity analysis shows that the canonical helicity ∫P·Q dV as a whole must be considered when analyzing reconnection. A mechanism for whistler wave generation and propagation is also described, with comparisons to recent spacecraft observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum Ornstein-Zernike theory for two-temperature two-component plasmas

Laboratory plasma production almost always preferentially heats either the ions or electrons, leading to a two-temperature state. In this state, density functional theory molecular dynamic simulation is the state of the art for modeling bulk material properties. We construct a statistical mechanics model for the two temperature limit that is theoretically consistent with the molecular dynamics method. We proceed to derive the electron-ion multi-temperature quantum Ornstein-Zernike equations for the first time. This allows the construction of a two-temperature two-component plasma model using the average atom from which we can compute bulk material properties at a fraction of the computation time of the two-temperature density functional theory simulation. The accuracy of the model is benchmarked against ion pair correlation and self-diffusion results from ab initio simulation. Here, we proceed to compute the viscosity and ion thermal conductivity as a function of both ion and electron temperature.

Ab initio molecular dynamics↗

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3

In order to provide a comprehensive theoretical description of MgSiO 3 at extreme conditions, we combine results from path integral Monte Carlo and density functional molecular dynamics simulations and generate a consistent equation of state for this material. We consider a wide range of temperature and density conditions from 10 4 to 10 8 K and from 0.321 to 64.2 g/cm -3 (0.1- to 20-fold the ambient density). We study how the L and K shell electrons are ionized with increasing temperature and pressure. We derive the shock Hugoniot curve and compare with experimental results. Our Hugoniot curve is in good agreement with the experiments, and we predict a broad compression maximum that is dominated by the K shell ionization of all three nuclei while the peak compression ratio of 4.70 is obtained when the Si and Mg nuclei are ionized. To conclude, we analyze the heat capacity and structural properties of the liquid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A theory of self-organized zonal flow with fine radial structure in tokamak

The (low frequency) zonal flow-ion temperature gradient (ITG) wave system, constructed on Braginskii's fluid model in tokamak, is shown here to be a reaction-diffusion-advection system; it is derived by making use of a multiple spatiotemporal scale technique and two-dimensional (2D) ballooning theory. For real regular group velocities of ITG waves, two distinct temporal processes, sharing a very similar meso-scale radial structure, are identified in the nonlinear self-organized stage. The stationary and quasi-stationary structures reflect a particular feature of the poloidal group velocity. The equation set posed to be an initial value problem is numerically solved for JET low mode parameters; the results are presented in several figures and two movies that show the spatiotemporal evolutions as well as the spectrum analysis—frequency-wave number spectrum, auto power spectrum, and Lissajous diagram. This approach reveals that the zonal flow in tokamak is a local traveling wave. For the quasi-stationary process, the cycle of ITG wave energy is composed of two consecutive phases in distinct spatiotemporal structures: a pair of Cavitons growing and breathing slowly without long range propagation, followed by a sudden decay into many Instantons that carry negative wave energy rapidly into infinity. A spotlight onto the motion of Instantons for a given radial position reproduces a Blob-Hole temporal structure; the occurrence as well as the rapid decay of Caviton into Instantons is triggered by zero-crossing of radial group velocity. A sample of the radial profile of zonal flow contributed from 31 nonlinearly coupled rational surfaces near plasma edge is found to be very similar to that observed in the JET Ohmic phase [J. C. Hillesheim et al., Phys. Rev. Lett. 116, 165002 (2016)]. The theory predicts an interior asymmetric dipole structure associated with the zonal flow that is driven by the gradients of ITG turbulence intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fast, purely growing collisionless reconnection as an eigenfunction problem related to but not involving linear whistler waves

If either finite electron inertia or finite resistivity is included in 2D magnetic reconnection, the two-fluid equations become a pair of second-order differential equations coupling the out-of-plane magnetic field and vector potential to each other to form a fourth-order system. The coupling at an X-point is such that out-of-plane even-parity electric and odd-parity magnetic fields feed off each other to produce instability if the scale length on which the equilibrium magnetic field changes is less than the ion skin depth. The instability growth rate is given by an eigenvalue of the fourth-order system determined by boundary and symmetry conditions. The instability is a purely growing mode, not a wave, and has growth rate of the order of the whistler frequency. The spatial profile of both the out-of-plane electric and magnetic eigenfunctions consists of an inner concave region having extent of the order of the electron skin depth, an intermediate convex region having extent of the order of the equilibrium magnetic field scale length, and a concave outer exponentially decaying region. If finite electron inertia and resistivity are not included, the inner concave region does not exist and the coupled pair of equations reduces to a second-order differential equation having non-physical solutions at an X-point.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗