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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

GDB: A global 3D two-fluid model of plasma turbulence and transport in the tokamak edge

Here, the global drift-ballooning (GDB) code is developed to study tokamak edge low frequency turbulence and transport, and their relationship to global profile evolution. The code employs a 3D electromagnetic fluid model that does not discriminate between equilibrium and perturbative contributions, capturing arbitrary amplitude fluctuations. Primitive plasma variables, including the flow profiles, are evolved self-consistently in both closed-flux surfaces and the scrape-off-layer (SOL). A suite of numerical techniques described here handle the linear and non-linear components of the model efficiently so as to support realistic discharge parameters (such as realistic deuterium mass ratio √m i /m e ≈ 60) and yield good scaling on high performance computing (HPC) systems. GDB resolves turbulence slower than the ion gyrofrequency in simulations that capture the millisecond-scale evolution of global plasma profiles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Characterization of rarefaction waves in van der Waals fluids

Here, we calculate the isentropic evolution of an instantaneously heated foil, assuming a van der Waals equation of state with the Maxwell construction. The analysis by Yuen and Barnard [Phys. Rev. E 92, 033019 (2015)] is extended for the particular case of three degrees of freedom. We assume heating to temperatures in the vicinity of the critical point. The self-similar profiles of the rarefaction waves describing the evolution of the foil display plateaus in density and temperature due to a phase transition from the single-phase to the two-phase regime. The hydrodynamic equations are expressed in a dimensionless form and the solutions form a set of universal curves, depending on a single parameter: the dimensionless initial entropy. We characterize the rarefaction waves by calculating how the plateau length, density, pressure, temperature, velocity, internal energy, and sound speed vary with dimensionless initial entropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Direct construction of optimized stellarator shapes. Part 2. Numerical quasisymmetric solutions

Quasisymmetric stellarators are appealing intellectually and as fusion reactor candidates since the guiding-centre particle trajectories and neoclassical transport are isomorphic to those in a tokamak, implying good confinement. Previously, quasisymmetric magnetic fields have been identified by applying black-box optimization algorithms to minimize symmetry-breaking Fourier modes of the field strength $B$ . Here, instead, we directly construct magnetic fields in cylindrical coordinates that are quasisymmetric to leading order in the distance from the magnetic axis, without using optimization. The approach involves solution of a one-dimensional nonlinear ordinary differential equation, originally derived by Garren & Boozer ( Phys. Fluids B, vol. 3, 1991, p. 2805). We illustrate the usefulness and accuracy of this optimization-free approach by providing the results of this construction as input to the codes VMEC and BOOZ_XFORM, confirming the purity and scaling of the magnetic spectrum. The space of magnetic fields that are quasisymmetric to this order is parameterized by the magnetic axis shape along with three other real numbers, one of which reflects the on-axis toroidal current density, and another one of which is zero for stellarator symmetry. The method here could be used to generate good initial conditions for conventional optimization, and its speed enables exhaustive searches of parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Coupling extended magnetohydrodynamic fluid codes with radiofrequency ray tracing codes for fusion modeling

Neoclassical tearing modes are macroscopic ($L$ ~ 1 m) instabilities in magnetic fusion experiments; if unchecked, these modes degrade plasma performance and may catastrophically destroy plasma confinement by inducing a disruption. Fortunately, the use of properly tuned and directed radiofrequency waves ($λ$ ~ 1 mm) can eliminate these modes. Numerical modeling of this difficult multiscale problem requires the integration of separate mathematical models for each length and time scale [T. G. Jenkins, S. E. Kruger, Phys. Plasmas 19 (2012) 122508]; the extended MHD model captures macroscopic plasma evolution while the RF model tracks the flow and deposition of injected RF power through the evolving plasma profiles. The scale separation enables use of the eikonal (ray-tracing) approximation to model the RF wave propagation. In this work we demonstrate a technique, based on methods of computational geometry, for mapping the ensuing RF data (associated with discrete ray trajectories) onto the finite-element/pseudospectral grid that is used to model the extended MHD physics. In the new representation, the RF data can then be used to construct source terms in the equations of the extended MHD model, enabling quantitative modeling of RF-induced tearing mode stabilization. Though our specific implementation uses the NIMROD extended MHD [C. R. Sovinec et al., J. Comp. Phys. 195 (2004) 355] and GENRAY RF [A. P. Smirnov, R. W. Harvey, K. Kupfer, Bull. Amer. Phys. Soc. 39 (1994) 1626] codes, the approach presented can be applied more generally to any code coupling requiring the mapping of ray tracing data onto Eulerian grids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A comparison of explosively driven shock wave radius versus time scaling approaches

Abstract Explosively driven shock wave radius versus time profiles are frequently used to document energy release and relative explosive performance. Recently, two universal shock wave radius versus time profiles have been presented in the literature, which demonstrate the ability to represent explosively driven shock wave profiles for all explosive sources in any fluid environment. These two universal shock wave profiles are examined here relative to each other and relative to a commonly used nonlinear shock wave profile, which is fit to experimental data for individual explosive materials. The nonlinear profile, originally developed by Dewey, is examined here, and a universal non-dimensional form of the equation is proposed. The universal shock wave profiles are all found to be relatively similar, but with slight variations in a transition region of non-dimensional radii $$0.15\lesssim R^*\lesssim 2$$ 0.15 ≲ R ∗ ≲ 2 . The variations in this region result in different estimations of energy release or blast strength between the curve fits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Numerical study of He/CF{sub 3}I pulsed discharge used to produce iodine atom in chemical oxygen-iodine laser

The pulsed discharge for producing iodine atoms from the alkyl and perfluoroalky iodides (CH{sub 3}I, CF{sub 3}I, etc.) is the most efficient method for achieving the pulse operating mode of a chemical oxygen-iodine laser. In this paper, a one-dimensional fluid model is developed to study the characteristics of pulsed discharge in CF{sub 3}I-He mixture. By solving continuity equation, momentum equation, Poisson equation, Boltzmann equation, and an electric circuit equation, the temporal evolution of discharge current density and various discharge products, especially the atomic iodine, are investigated. The dependence of iodine atom density on discharge parameters is also studied. The results show that iodine atom density increases with the pulsed width and pulsed voltage amplitude. The mixture ratio of CF{sub 3}I and helium plays a more significant role in iodine atom production. For a constant voltage amplitude, there exists an optimal mixture ratio under which the maximum iodine atom concentration is achieved. The bigger the applied voltage amplitude is, the higher partial pressure of CF{sub 3}I is needed to obtain the maximum iodine atom concentration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Simulation of edge localized mode heat pulse using drift-kinetic ions and Boltzmann electrons

This article reports on 1D+2V heat pulse propagation studies using the COGENT guiding center kinetic code. The model uses magnetized kinetic ions and a simple Boltzmann electron model. Results agree with previous kinetic and fluid modeling benchmark studies that correspond to the parameters of edge localized modes (ELMs) observed on the JET tokamak. The plasma parameters for the edge pedestal and ensuing ELM dynamics are in the low collisionality regime. Hence, the dominant balance between the assumed Maxwellian ELM source and collisionless parallel advection causes the ion PDF to develop a significantly anisotropic velocity distribution. Adding nonlinear Coulomb ion-ion collisions to the model acts to smooth the sharp features of the ion distribution function, but the anisotropy remains robust due to the low collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A general metriplectic framework with application to dissipative extended magnetohydrodynamics

General equations for conservative yet dissipative (entropy producing) extended magnetohydrodynamics are derived from two-fluid theory. Keeping all terms generates unusual cross-effects, such as thermophoresis and a current viscosity that mixes with the usual velocity viscosity. While the Poisson bracket of the ideal version of this model has already been discovered, we determine its metriplectic counterpart that describes the dissipation. This is done using a new and general thermodynamic point of view to derive dissipative brackets, a means of derivation that is natural for understanding and creating dissipative dynamics without appealing to underlying kinetic theory orderings. Finally, the formalism is used to study dissipation in the Lagrangian variable picture where, in the context of extended magnetohydrodynamics, non-local dissipative brackets naturally emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8

Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A constrained-transport embedded boundary method for compressible resistive magnetohydrodynamics

Motivated by the increased interest in pulsed-power magneto-inertial fusion devices in recent years, we present a method for implementing an arbitrarily shaped embedded boundary on a Cartesian mesh while solving the equations of compressible resistive magnetohydrodynamics. The method is built around a finite volume formulation of the equations in which a Riemann solver is used to compute fluxes on the faces between grid cells, and a face-centered constrained transport formulation of the induction equation. The small time step problem associated with the cut cells is avoided by always computing fluxes on the faces and edges of the Cartesian mesh. We extend the method to model a moving interface between two materials with different properties using a ghost-fluid approach, and show some preliminary results including shock-wave-driven and magnetically-driven dynamical compressions of magnetohydrostatic equilibria. In conclusion, we present a thorough verification of the method and show that it converges at second order in the absence of discontinuities, and at first order with a discontinuity in material properties.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

The role of the electron convection term for the parallel electric field and electron acceleration in MHD simulations

There has been a great concern about the origin of the parallel electric field in the frame of fluid equations in the auroral acceleration region. This paper proposes a new method to simulate magnetohydrodynamic (MHD) equations that include the electron convection term and shows its efficiency with simulation results in one dimension. We apply a third-order semi-discrete central scheme to investigate the characteristics of the electron convection term including its nonlinearity. At a steady state discontinuity, the sum of the ion and electron convection terms balances with the ion pressure gradient. We find that the electron convection term works like the gradient of the negative pressure and reduces the ion sound speed or amplifies the sound mode when parallel current flows. The electron convection term enables us to describe a situation in which a parallel electric field and parallel electron acceleration coexist, which is impossible for ideal or resistive MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Fast elliptic solvers in cylindrical coordinates and the Coulomb collision operator

Highlights: {yields} We describe new fast solvers for elliptic partial differential equations in free space. {yields} We combine integral equation methods with Fourier methods to achieve high order accuracy. {yields} We apply these solvers to the evaluation of the Coulomb collision operator in plasma physics. - Abstract: In this paper, we describe a new class of fast solvers for separable elliptic partial differential equations in cylindrical coordinates (r, {theta}, z) with free-space radiation conditions. By combining integral equation methods in the radial variable r with Fourier methods in {theta} and z, we show that high-order accuracy can be achieved in both the governing potential and its derivatives. A weak singularity arises in the Fourier transform with respect to z that is handled with special purpose quadratures. We show how these solvers can be applied to the evaluation of the Coulomb collision operator in kinetic models of ionized gases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Importance of Ambipolar Electric Field in Driving Ion Loss From Mars: Results From a Multifluid MHD Model With the Electron Pressure Equation Included

The multi-fluid (MF) magnetohydrodynamic model of Mars is improved by solving an additional electron pressure equation. Through the electron pressure equation, the electron temperature is calculated based on the effects from various electron-related heating and cooling processes (e.g., photo-electron heating, electron-neutral collision, and electron-ion collision), and thus, the improved model can calculate the electron temperature and the electron pressure force terms self-consistently. Model results of a typical case using the MF with electron pressure equation included model are compared in detail to identical cases using the MF and multi-species models to identify the effect of the improved physics. Here, we find that when the electron pressure equation is included, the general interaction patterns are similar to those with no electron pressure equation. However, the MF with electron pressure equation included model predicts that the electron temperature is much larger than the ion temperature in the ionosphere, consistent with both Viking and Mars Atmosphere and Volatile EvolutioN (MAVEN) observations. Using our numerical model, we also examined in detail the relative importance of different forces in the plasma interaction region. All three models are also applied to a MAVEN event study using identical input conditions; overall, the improved model matches best with MAVEN observations. All of the simulation cases are examined in terms of the total ion loss, and the results show that the inclusion of the electron pressure equation increases the escape rates by 50–110% in total mass, depending on solar condition and strong crustal field orientation, clearly demonstrating the importance of the ambipolar electric field in facilitating ion escape.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Experimental verification of the Boltzmann relation at the wall in inductively coupled plasmas

The radial profiles of the electron density at the bottom wall were measured from the center to the edge at argon 3 mTorr with an RF power of 200 W. Since the Boltzmann relation in terms of the fluid picture only describes the behavior of the confined particles, the escaping electrons, which are not confined in plasmas, were not expected to follow the Boltzmann relation. However, unexpected results were observed. The spatial profiles of the density of escaping electrons and the floating potential were well described by the Boltzmann relation with the electron temperature. Moreover, it was concluded that the behavior of the escaping electrons is also described by the generalized Boltzmann relation in terms of the kinetic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Structure and computation of two-dimensional incompressible extended MHD

A comprehensive study of the extended magnetohydrodynamic model obtained from the two-fluid theory for electrons and ions with the enforcement of quasineutrality is given. Starting from the Hamiltonian structure of the fully three-dimensional theory, a Hamiltonian two-dimensional incompressible four-field model is derived. In this way, the energy conservation along with four families of Casimir invariants is naturally obtained. The construction facilitates various limits leading to the Hamiltonian forms of Hall, inertial, and ideal MHD, with their conserved energies and Casimir invariants. Basic linear theory of the four-field model is treated, and the growth rate for collisionless reconnection is obtained. Finally, results from nonlinear simulations of collisionless tearing are presented and interpreted using, in particular, normal fields, a product of the Hamiltonian theory that gives rise to simplified equations of motion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Two-point spectral model for variable-density homogeneous turbulence

We present a study of buoyancy-driven variable-density homogeneous turbulence, using a two-point spectral closure model. We compute the time-evolution of the spectral distribution in wave number $k$ of the correlation of density and specific volume $b$($k$), the velocity associated with the turbulent mass flux $a$($k$) , and the turbulent kinetic energy $E$($k$) , using a set of coupled equations. Under the modeling assumptions, each dynamical variable has two coefficients governing spectral transfer among modes. In addition, the velocity $a$($k$) has two coefficients governing the drag between the two fluids. Using a prescribed initial condition for $b$($k$) and starting from a quiescent flow, we first evaluate the relative importance of the different coefficients used to model this system and their impact on the statistical quantities. We next assess the accuracy of the model, relative to direct numerical simulation of the complete hydrodynamical equations, using $b$, $a$, and $E$ as metrics. We show that the model is able to capture the spectral distribution and global means of all three statistical quantities at both low and high Atwood number for a set of optimized coefficients. The optimization procedure also permits us to discern a minimal set of four coefficients which are sufficient to yield reasonable results while pointing to the mechanisms that dominate the mixing process in this problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The kinetic analog of the pressure–strain interaction

Energy transport in weakly collisional plasma systems is often studied with fluid models and diagnostics. However, the applicability of fluid models is limited when collisions are weak or absent, and using a fluid approach can obscure kinetic processes that provide key insights into the physics of energy transport. Kinetic diagnostics retain all of the information in 3D-3V phase space and thereby reach beyond the insights of fluid models to elucidate the mechanisms responsible for collisionless energy transport. In this work, we derive the Kinetic Pressure–Strain (KPS): a kinetic analog of the pressure–strain interaction, which is the channel between flow energy density and internal energy density in fluid models. Through two case studies of electron Landau damping, we demonstrate that the KPS diagnostic can elucidate kinetic mechanisms that are responsible for energy transport in this channel, just as the related field–particle correlation is known to identify kinetic mechanisms of transport between electromagnetic field energy density and kinetic energy density in particle flows. In addition, we show that resonant electrons play a major role in transferring energy between fluid flows and internal energy during the process of Landau damping.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗