Search NASA⌕ Search

SEARCH · Search NASA

Results for “plasma fluid modeling”

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 217 records · Page 12

Characterizing magnetized plasmas with dynamic mode decomposition

Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas

Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computer simulation of the effect of dDAVP with saline loading on fluid balance after 24-hour head-down tilt

Fluid loading (FL) before Shuttle reentry is a countermeasure currently in use by NASA to improve the orthostatic tolerance of astronauts during reentry and postflight. The fluid load consists of water and salt tablets equivalent to 32 oz (946 ml) of isotonic saline. However, the effectiveness of this countermeasure has been observed to decrease with the duration of spaceflight. The countermeasure's effectiveness may be improved by enhancing fluid retention using analogs of vasopressin such as lypressin (LVP) and desmopressin (dDAVP). In a computer simulation study reported previously, we attempted to assess the improvement in fluid retention obtained by the use of LVP administered before FL. The present study is concerned with the use of dDAVP. In a recent 24-hour, 6 degree head-down tilt (HDT) study involving seven men, dDAVP was found to improve orthostatic tolerance as assessed by both lower body negative pressure (LBNP) and stand tests. The treatment restored Luft's cumulative stress index (cumulative product of magnitude and duration of LBNP) to nearly pre-bedrest level. The heart rate was lower and stroke volume was marginally higher at the same LBNP levels with administration of dDAVP compared to placebo. Lower heart rates were also observed with dDAVP during stand test, despite the lower level of cardiovascular stress. These improvements were seen with only a small but significant increase in plasma volume of approximately 3 percent. This paper presents a computer simulation analysis of some of the results of this HDT study.

NASA Discipline Cardiopulmonary↗

Global Λ-hyperon polarization in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV

Global hyperon polarization, $\overline{P}$ H , in Au + Au collisions over a large range of collision energy, $\sqrt{S_{NN}}$, was recently measured and successfully reproduced by hydrodynamic and transport models with intense fluid vorticity of the quark-gluon plasma. While na¨ıve extrapolation of data trends suggests a large $\overline{P}$ H as the collision energy is reduced, the behavior of $\overline{P}$ H at small $\sqrt{S_{NN}}$ < 7.7 GeV is unknown. Operating the STAR experiment in fixed-target mode, we measured the polarization of Λ hyperons along the direction of global angular momentum in Au + Au collisions at $\sqrt{S_{NN}}$ = 3 GeV. As a result, the observation of substantial polarization of 4.91 ± 0.81 (stat.) ± 0.15 (syst.) % in these collisions may require a reexamination of the viscosity of any fluid created in the collision, of the thermalization timescale of rotational modes, and of hadronic mechanisms to produce global polarization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Theory of the low-latitude boundary layer

A one-dimensional steady state fluid mechanical model is developed of the low-latitude plasma boundary layer inside the dawn and dusk magnetopause. Momentum transfer in the layer is produced by viscosity and/or mass diffusion. Coupling to the ionosphere is achieved via field-aligned currents, the magnitude of which is limited by parallel potential drops. These currents flow into and out of the ionosphere in the manner described by Iijima and Potemra. The higher-latitude (region 1) currents are associated with the boundary layer proper, while the lower-latitude (region 2) ones are associated with a region of sunward return flow adjacent to the boundary layer. The parallel potential drops have a magnitude of typically 2-3 kV and a north-south extent of 100-200 km. The calculated potential profile corresponds reasonably well to observed inverted V precipitation events.

Sonnerup, B. U. OE.↗

Comments on the photospheric dynamo model of Henoux and Somov

A detailed model for a photospheric dynamo has been presented by Henoux and Somov (1987), who used the three-fluid model to treat the properties of the weakly ionized plasma. Only the equations for the two ionized components were solved. The equation for the neutral component is considered, and it is argued that the model is unacceptable becaused of an implied impossibly large unbalanced stress on the neutral gas. It is argued more generally that all existing photospheric dynamo models are untenable.

Melrose, D. B.↗

Multifluid Simulations of the Global Solar Wind Including Pickup Ions and Turbulence Modeling

I will describe a three-dimensional magnetohydrodynamic model of the solar wind that takes into account turbulent heating of the wind by velocity and magnetic fluctuations as well as a variety of effects produced by interstellar pickup protons. The interstellar pickup protons are treated in the model as one fluid and the protons and electrons are treated together as a second fluid. The model equations include a Reynolds decomposition of the plasma velocity and magnetic field into mean and fluctuating quantities, as well as energy transfer from interstellar pickup protons to solar wind protons that results in the deceleration of the solar wind. The model is used to simulate the global steady-state structure of the solar wind in the region from 0.3 to 100 AU. The simulation assumes that the background magnetic field on the Sun is either a dipole (aligned or tilted with respect to the solar rotation axis) or one that is deduced from solar magnetograms.

Goldstein, Melvyn L.↗

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A two-fluid solar-wind model with intermittent Alfvénic turbulence

In one of the leading theories for the origin of the solar wind, photospheric motions launch Alfvén waves (AWs) that propagate along open magnetic-field lines through the solar atmosphere and into the solar wind. The radial variation in the Alfvén speed causes some of the AWs to reflect, and counter-propagating AWs subsequently interact to produce Alfveńic turbulence, in which AW energy cascades from long wavelengths to short wavelengths and dissipates, heating the plasma. In this paper we develop a one-dimensional two-fluid solar-wind model that includes Alfvénic turbulence, proton temperature anisotropy and a novel method for apportioning the turbulent heating rate between parallel proton heating, perpendicular proton heating and electron heating. We employ a turbulence model that accounts for recent observations from NASA’s Parker Solar Probe, which find that AW fluctuations in the near-Sun solar wind are intermittent and less anisotropic than in previous models of anisotropic magnetohydrodynamic turbulence. Our solar-wind model reproduces a wide range of remote observations of the corona and in-situ measurements of the solar wind, and our turbulent heating model consists of analytic equations that could be usefully incorporated into other solar-wind models and numerical models of more distant astrophysical plasmas.

astrophysical plasmas↗

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Nanosecond-Pulsed Spark Discharge and Flame Kernel Evolution

Dilute combustion, either using exhaust gas recirculation or with excess air, is considered a promising strategy to improve the thermal efficiency of internal combustion engines. However, the dilute air-fuel mixture, especially under intensified turbulence and high-pressure conditions, poses significant challenges for ignitability and combustion stability, which may limit the attainable efficiency benefits. In-depth knowledge of the flame kernel evolution to stabilize ignition and combustion in a challenging environment is crucial for effective engine development and optimization. To date, a comprehensive understanding of ignition processes that result in the development of fully predictive ignition models usable by the automotive industry does not yet exist. Spark-ignition consists of a wide range of physics that includes electrical discharge, plasma evolution, joule-heating of gas, and flame kernel initiation and growth into a self-sustainable flame. In this study, an advanced approach is proposed to model spark-ignition energy deposition and flame kernel growth. To decouple the flame kernel growth from the electrical discharge, a nanosecond-pulsed high-voltage discharge is used to trigger spark-ignition in an optically accessible small ignition test vessel with a quiescent mixture of air and methane. Initial conditions for the flame kernel, including its thermodynamic state and species composition, are derived from a plasma-chemical equilibrium calculation. The geometric shape and dimension of the kernel are characterized using a multi-dimensional thermal plasma solver. Here, the proposed modeling approach is evaluated using a high-fidelity computational fluid dynamics procedure to compare the simulated flame kernel evolution against flame boundaries from companion Schlieren images.

42 ENGINEERING↗

Development, Implementation, and Verification of Partially-Ionized Collisional Multifluid Plasma Models in Drekar

We consider the development of multifluid models for partially ionized multispecies plasmas. The models are composed of a standard set of five-moment fluid equations for each species plus a description of electromagnetics. The most general model considered utilizes a full set of fluid equations for each charge state of each atomic species, plus a set of fluid equations for electrons. The fluid equations are coupled through source terms describing electromagnetic coupling, ionization, recombination, charge exchange, and elastic scattering collisions in the low-density coronal limit. The form of each of these source terms is described in detail, and references for required rate coefficients are identified for a diverse range of atomic species. Initial efforts have been made to extend these models to incorporate some higher-density collisional effects, including ionization potential depression and three- body recombination. Some reductions of the general multifluid model are considered. First, a reduced multifluid model is derived which averages over all of the charge states (including neutrals) of each atomic species in the general multifluid model. The resulting model maintains full consistency with the general multifluid model from which it is derived by leveraging a quasi-steady-state collisional ionization equilibrium assumption to recover the ionization fractions required to make use of the general collision models. Further reductions are briefly considered to derive certain components of a single-fluid magnetohydrodynamics (MHD) model. In this case, a generalized Ohm's law is obtained, and the standard MHD resistivity is expressed in terms of the collisional models used in the general multifluid model. A number of numerical considerations required to obtain robust implementations of these multifluid models are discussed. First, an algebraic flux correction (AFC) stabilization approach for a continuous Galerkin finite element discretization of the multifluid system is described in which the characteristic speeds used in the stabilization of the fluid systems are synchronized across all species in the model. It is demonstrated that this synchronization is crucial in order to obtain a robust discretization of the multifluid system. Additionally, several different formulations are considered for describing the electromagnetics portion of the multifluid system using nodal continuous Galerkin finite element discretizations. The formulations considered include a parabolic divergence cleaning method and an implicit projection method for the traditional curl formulation of Maxwell's equations, a purely- hyperbolic potential-based formulation of Maxwell's equations, and a mixed hyperbolic-elliptic potential-based formulation of Maxwell's equations. Some advantages and disadvantages of each formulation are explored to compare solution robustness and the ease of use of each formulation. Numerical results are presented to demonstrate the accuracy and robustness of various components of our implementation. Analytic solutions for a spatially homogeneous damped plasma oscillation are derived in order to verify the implementation of the source terms for electromagnetic coupling and elastic collisions between fluid species. Ionization balance as a function of electron temperature is evaluated for several atomic species of interest by comparing to steady-state calculations using various sets of ionization and recombination rate coefficients. Several test problems in one and two spatial dimensions are used to demonstrate the accuracy and robustness of the discretization and stabilization approach for the fluid components of the multifluid system. This includes standard test problems for electrostatic and electromagnetic shock tubes in the two-fluid and ideal shock-MHD limits, a cylindrical diocotron instability, and the GEM challenge magnetic reconnection problem. A one-dimensional simplified prototype of an argon gas puff configuration as deployed on Sandia's Z-machine is used as a demonstration to exercise the full range of capabilities associated with the general multifluid model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A comparison of two plasma models

The time dependent behavior of a plasma which surrounds a highly biased conducting sphere is considered. This plasma is treated as either a cold two component fluid or as a warm plasma whose time development can be found by solving the Vlasov equation. Both models demonstrate oscillatory behavior, but the electric fields predicted by the models are quantitatively quite different in regions close to the surface of the sphere and very similar otherwise. A broadening of the electron distribution function with time is observed indicating local heating of the plasma near the surface of the sphere.

Cottam, Russell↗

Plasma Sheet Circulation Pathways

Global simulations of Earth's magnetosphere in the solar wind compute the pathways of plasma circulation through the plasma sheet. We address the pathways that supply and drain the plasma sheet, by coupling single fluid simulations with Global Ion Kinetic simulations of the outer magnetosphere and the Comprehensive Ring Current Model of the inner magnetosphere, including plasmaspheric plasmas. We find that the plasma sheet is supplied with solar wind plasmas via the magnetospheric flanks, and that this supply is most effective for northward IMF. For southward IMF, the innermost plasma sheet and ring current region are directly supplied from the flanks, with an asymmetry of single particle entry favoring the dawn flank. The central plasma sheet (near midnight) is supplied, as expected, from the lobes and polar cusps, but the near-Earth supply consists mainly of slowly moving ionospheric outflows for typical conditions. Work with the recently developed multi-fluid LFM simulation shows transport via plasma "fingers" extending Earthward from the flanks, suggestive of an interchange instability. We investigate this with solar wind ion trajectories, seeking to understand the fingering mechanisms and effects on transport rates.

Moore, Thomas E.↗

Nonlinear modeling of the scaling law for the m/n = 3/2 error field penetration threshold

The scaling law for the m/n=3/2 error field (EF) penetration threshold is predicted numerically based on nonlinear single-fluid and two-fluid modeling using the TM1 code. The simulated penetration threshold of radial magnetic field b r at the plasma edge is scaled to the electron density n e , temperature T e , viscous time τ μ , toroidal field B t and the natural frequency ω in the form of b r /B t ∝n e αn T e αT τ μ αμ B t αB ω αω by scanning these parameters separately. Here, α n , α T , α μ , α B and α ω are the scaling coefficients on n e , T e , τ μ , B t and ω, respectively. Single-fluid modeling shows that the 3/2 EF threshold scales as b r /B t ∝n e 0.56 T e 0.6 τ μ -0.59 B t -1.15 ω, which is similar with the analytical scaling law in both the Rutherford and visco-resistive regimes. Yet, two-fluid modeling shows that the scaling law differs significantly in particular regarding the dependence on plasma rotation. In detail, the scaling coefficient α n on density decreases from 0.67 to 0.56 and α T on temperature decreases from 0.67 to 0.32, while α μ on viscous time is around -0.45 and α B on toroidal field decreases slightly from -1.15 to -1, when the ratio |ω E /ω *e | between plasma rotation frequency ω E and diamagnetic drift frequency ω *e varies from 0 to 10. Scans of the plasma rotation reveals that the penetration threshold linearly depends on the perpendicular rotation frequency (or natural frequency) ω ⊥e =ω E +ω *e , and there is a minimum in the required field amplitude when ω ⊥e 0. In addition, the enduring mystery of non-zero penetration threshold at zero plasma natural frequency in EF experiments is resolved by two-fluid simulations. We report that the very small island and smooth bifurcation in EF penetration near zero frequency is hard to detect in the experiment, leading to a finite penetration threshold within the capability of the experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solar wind flow past Venus - Theory and comparisons

Advanced computational procedures are applied to an improved model of solar wind flow past Venus to calculate the locations of the ionopause and bow wave and the properties of the flowing ionosheath plasma in the intervening region. The theoretical method is based on a single-fluid, steady, dissipationless, magneto-hydrodynamic continuum model and is appropriate for the calculation of axisymmetric supersonic, super-Alfvenic solar wind flow past a nonmagnetic planet possessing a sufficiently dense ionosphere to stand off the flowing plasma above the subsolar point and elsewhere. Determination of time histories of plasma and magnetic field properties along an arbitrary spacecraft trajectory and provision for an arbitrary oncoming direction of the interplanetary solar wind have been incorporated in the model. An outline is provided of the underlying theory and computational procedures, and sample comparisons of the results are presented with observations from the Pioneer Venus orbiter.

Spreiter, J. R.↗

Combined Influence of Rotation and Scrape-Off Layer Drifts on Recycling Asymmetries in Tokamak Plasmas

Coupled 2D fluid-kinetic simulations of a DIII-D high confinement tokamak plasma show that plasma rotation coupled with drift effects near the plasma edge play a significant role in the creation of the observed poloidal distribution of neutrals. It is observed that including either drift or rotation effects enhances particle flux at the inner target in the case of ion 𝐵×∇𝐵 drift toward the 𝑋-point. However, the particle flux asymmetry is significantly higher with the combination of drifts and rotation than either effect alone. The heightened particle flux asymmetry allows for improved simulation of the strong in-out asymmetry of the Lyman-𝛼 brightness profiles measured in the experiment. Enhancement of radial transport of parallel momentum changes the upstream scrape-off layer flow pattern, increasing the fraction of deuterium flux that reaches the inboard divertor entrance while lowering that which arrives at the outboard. In conclusion, this Letter indicates that by combining drifts, rotation, and viscous coupling, existing boundary plasma models can achieve a satisfactory agreement with experimentally measured neutral asymmetries.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗