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 361 records · Page 20

The role of magnetohydrodynamics in heliospheric space plasma physics research

Magnetohydrodynamics (MHD) is a fairly recent extension of the field of fluid mechanics. While much remains to be done, it has successfully been applied to the contemporary field of heliospheric space plasma research to evaluate the 'macroscopic picture' of some vital topics via the use of conducting fluid equations and numerical modeling and simulations. Some representative examples from solar and interplanetary physics are described to demonstrate that the continuum approach to global problems (while keeping in mind the assumptions and limitations therein) can be very successful in providing insight and large scale interpretations of otherwise intractable problems in space physics.

Dryer, Murray↗

Science Enhancements by the MAVEN Participating Scientists

NASA implemented a Participating Scientist Program and released a solicitation for the Mars Atmosphere and Volatile EvolutioN mission (MAVEN) proposals on February 14, 2013. After a NASA peer review panel evaluated the proposals, NASA Headquarters selected nine on June 12, 2013. The program's intent is to enhance the science return from the mission by including new investigations that broaden and/or complement the baseline investigations, while still addressing key science goals. The selections cover a broad range of science investigations. Included are: a patching of a 3D exosphere model to an improved global ionosphere-thermosphere model to study the generation of the exosphere and calculate the escape rates; the addition of a focused study of upper atmosphere variability and waves; improvement of a multi-fluid magnetohydrodynamic model that will be adjusted according to MAVEN observations to enhance the understanding of the solar-wind plasma interaction; a global study of the state of the ionosphere; folding MAVEN measurements into the Mars International Reference Ionosphere under development; quantification of atmospheric loss by pick-up using ion cyclotron wave observations; the reconciliation of remote and in situ observations of the upper atmosphere; the application of precise orbit determination of the spacecraft to measure upper atmospheric density and in conjunction with other Mars missions improve the static gravity field model of Mars; and an integrated ion/neutral study of ionospheric flows and resultant heavy ion escape. Descriptions of each of these investigations are given showing how each adds to and fits seamlessly into MAVEN mission science design.

Thermosphere↗

Science Enhancements by the MAVEN Participating Scientists

NASA implemented a Participating Scientist Program and released a solicitation for the Mars Atmosphere and Volatile EvolutioN mission (MAVEN) proposals on February 14, 2013. After a NASA peer review panel evaluated the proposals, NASA Headquarters selected nine on June 12, 2013. The program’s intent is to enhance the science return from the mission by including new investigations that broaden and/or complement the base-line investigations, while still addressing key science goals. The selections cover a broad range of science investigations. Included are: a patching of a 3D exosphere model to an improved global ionosphere-thermosphere model to study the generation of the exosphere and calculate the escape rates; the addition of a focused study of upper atmosphere variability and waves; improvement of a multi-fluid magnetohydrodynamic model that will be adjusted according to MAVEN observations to enhance the understanding of the solar-wind plasma interaction; a global study of the state of the ionosphere; folding MAVEN measurements into the Mars International Reference Ionosphere under development; quantification of atmospheric loss by pick-up using ion cyclotron wave observations; the reconciliation of remote and in situ observations of the upper atmosphere; the application of precise orbit determination of the spacecraft to measure upper atmospheric density and in conjunction with other Mars missions improve the static gravity field model of Mars; and an integrated ion/neutral study of ionospheric flows and resultant heavy ion escape. Descriptions of each of these investigations are given showing how each adds to and fits seamlessly into MAVEN mission science design.

Mars↗

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Final Technical Report - Center for Simulation of Fusion Relevant RF Actuators

We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-Dimensional Magnetohydrodynamic Modeling of the Solar Wind Including Pickup Protons and Turbulence Transport

To study the effects of interstellar pickup protons and turbulence on the structure and dynamics of the solar wind, we have developed a fully three-dimensional magnetohydrodynamic solar wind model that treats interstellar pickup protons as a separate fluid and incorporates the transport of turbulence and turbulent heating. The governing system of equations combines the mean-field equations for the solar wind plasma, magnetic field, and pickup protons and the turbulence transport equations for the turbulent energy, normalized cross-helicity, and correlation length. The model equations account for photoionization of interstellar hydrogen atoms and their charge exchange with solar wind protons, energy transfer from pickup protons to solar wind protons, and plasma heating by turbulent dissipation. Separate mass and energy equations are used for the solar wind and pickup protons, though a single momentum equation is employed under the assumption that the pickup protons are comoving with the solar wind protons.We compute the global structure of the solar wind plasma, magnetic field, and turbulence in the region from 0.3 to 100 AU for a source magnetic dipole on the Sun tilted by 0 deg - .90 deg and compare our results with Voyager 2 observations. The results computed with and without pickup protons are superposed to evaluate quantitatively the deceleration and heating effects of pickup protons, the overall compression of the magnetic field in the outer heliosphere caused by deceleration, and the weakening of corotating interaction regions by the thermal pressure of pickup protons.

Usmanov, Arcadi V.↗

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of high-speed (orbital) releases of electron attachment materials in the ionosphere

The dynamics of ionospheric plasma irregularities produced by the release of electron attachment materials at orbital velocities across the geomagnetic field is studied. A two-dimensional electrostatic fluid model which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species transport is developed. Numerical simulations are performed to study the behavior at early and at late times after the release. At early times, of the order of or less than the attachment material neutral collision time, the negative ion cloud produced by the release may structure owing to the shear in the E x B velocity within the cloud. At high altitudes the cloud may bifurcate and form vortices on the back. At lower altitudes where ion-neutral collisional effects dominate, this structuring is suppressed. At late times, after a plasma depletion has formed due to neutralization chemistry, the cloud structures by the E x B interchange instability. Depending on the release altitude, the depletion structures by the collisional or inertial limit of this instability.

Scales, W. A.↗

Multidimensional modeling of non-equilibrium plasma generated by a radio-frequency corona discharge

Low-temperature plasma (LTP) ignition concepts rely on the production of radical and charged species to speed up the onset of combustion in spark-ignition engines. These features are responsible for the superior performance of LTP igniters under extremely dilute combustion operation that is not achievable by conventional spark igniters. Additionally, LTP discharges extend the lifetime of the igniters, due to the avoidance of spark processes. For these reasons, the engine research community and the automotive industry have shown growing interest in this technology in the recent years. As of today, computational fluid-dynamics (CFD) codes typically used by the multi-dimensional engine modeling community do not have reliable models to describe LTP ignition processes. One key missing piece of information is the physical and chemical properties of the plasma and their effect on combustion ignition. Most non-equilibrium plasma simulations reported in literature are based on simplified, canonical geometries, with simple discharge excitation schemes. Here we conduct multi-dimensional modeling of the non-equilibrium plasma generated by an application-relevant radio-frequency (RF) corona discharge in air. Three test cases are simulated, characterized by different environmental pressure levels and peak electrode voltage values at room temperature. Streamer penetration, electron number density, atomic oxygen production, and bulk gas temperature distribution in the first 10 sinusoidal pulses are presented and discussed. This model can be used as a key tool for an in-depth understanding of RF-corona discharge for automotive applications and provides the basis for future implementations of dedicated LTP ignition models in CFD codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An improved ten-moment closure for reconnection and instabilities

The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Turbulence Transition in Magnetically Confined Hydrogen and Deuterium Plasmas

In this study, we discovered a turbulence transition in a large helical device. The turbulence level and turbulence-driven energy transport decrease to a specific transition density and increase above it. The ruling turbulences below and above the transition density were ion-temperature gradient (ITG) and resistive-interchange (RI) turbulences, consistent with the predictions of gyrokinetic theory and two-fluid MHD model, respectively. Isotope experiments on hydrogen (H) and deuterium (D) clarified the role of transitions. In the ITG regime, turbulence levels and energy transport were comparable in the H and D plasmas. In contrast, in the RI regime, they were clearly suppressed in the D plasma. The results provide crucial knowledge for understanding isotope effects and future optimization of stellarator and heliotron devices. Published by the American Physical Society 2024

Physics↗

EUTERPE: A global gyrokinetic code for stellarator geometry

The current state of the EUTERPE code is described with emphasis on the implemented models and their numerical implementation. The code solves the multi-species electromagnetic gyrokinetic equations in the full volume of a three-dimensional domain. Noise reduction of the particle-in-cell method is achieved by using a δf-method and Fourier filters. The field equations are discretized with B-splines and the resulting system of equations is solved iteratively. For linear simulations a phase-factor transformation is applied in order to strongly reduce the necessary grid resolution. Apart from the full gyrokinetic model, other numerically less expensive hybrid models are also implemented. They are mainly tailored for comparison with fluid theory and for studying the interaction of the bulk plasma with fast particles. The code is parallelized for CPUs by particle and domain decomposition. Good scalability up to several thousand nodes is demonstrated.

97 MATHEMATICS AND COMPUTING↗

Thermally conductive magnetohydrodynamic flows in helmet-streamer coronal structures

The behavior of thermally conductive plasma flows in helmet-streamer coronal structures is investigated within the framework of the axisymmetric nonrotating one-fluid MHD model. Continuous subsonic-supersonic solutions satisfying observed boundary conditions at the sun as well as the vanishing of the temperature at infinity are obtained and presented. Special attention is paid to the combined effects of conductive flow (and corresponding thermal force) and rapidly diverging magnetic field on the critical points. In this, the heliocentric distance of the neutral point determining the separation between closed and open field lines (cusp) is treated as a free parameter. These thermally conductive solutions are contrasted with those provided by corresponding isothermal models.

Cuperman, S.↗

Interpretation of high-speed flows in the plasma sheet

Pursuing an idea suggested by Pontius and Wolf (1990), we propose that the `bursty bulk flows' observed by Baumjohann et al. (1990) and Angelopoulos et al. (1992) are `bubbles' in the Earth's plasma sheet. Specifically, they are flux tubes that have lower values of pV(exp 5/3) than their neighbors, where p is the thermal pressure of the particles and V is the volume of a tube containing one unit of magnetic flux. Whether they are created by reconnection or some other mechanism, the bubbles are propelled earthward by a magnetic buoyancy force, which is related to the interchange instability. Most of the major observed characteristics of the bursty bulk flows can be interpreted naturally in terms of the bubble picture. We propose a new `stratified fluid' picture of the plasma sheet, based on the idea that bubbles constitute the crucial transport mechanism. Results from simple mathematical models of plasma sheet transport support the idea that bubbles can resolve the pressure balance inconsistency, particularly in cases where plasma sheet ions are lost by gradient/curvature drift out the sides of the tail or bubbles are generated by reconnection in the middle of plasma sheet.

Chen, C. X.↗

Linear simulations of a wide pedestal quiescent H-mode plasma with the extended-MHD code NIMROD

We present linear simulations of a Wide Pedestal QH (WPQH)-mode DIII-D plasma with the NIMROD code to address the role of ion multispecies collisionality on WPQH-modes. We use a full extended-MHD model that captures Alfvénic and magnetoacoustic waves, including two-fluid, ion gyroviscosity and cross heat flux. We show that two-fluid and ion gyroviscous effects destabilize electron-directed peeling–ballooning modes in the pedestal. The linear growth rates of the modes decrease when carbon (C) is added to a deuterium plasma because the multispecies collisional effects increase the plasma resistivity. When replacing C with tungsten (W), the multispecies collisional effects increase the resistivity further, as well as the stabilizing effect on the instabilities. The resistivity impact on mode stability is reversed when both two-fluid and ion gyroviscous effects are ignored, which is consistent with previous works, proving that multispecies collisional effects, together with two-fluid and ion gyroviscosity, should be considered to address the stability of plasmas with significant impurity content.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Heating and Acceleration of Solar Wind Ions by Turbulent Wave Spectrum in Inhomogeneous Expanding Plasma

Near the Sun (< 10Rs) the acceleration, heating, and propagation of the solar wind are likely affected by the background inhomogeneities of the magnetized plasma. The heating and the acceleration of the solar wind ions by turbulent wave spectrum in inhomogeneous plasma is studied using a 2.5D hybrid model. The hybrid model describes the kinetics of the ions, while the electrons are modeled as massless neutralizing fluid in an expanding box approach. Turbulent magnetic fluctuations dominated by power-law frequency spectra, which are evident from in-situ as well as remote sensing measurements, are used in our models. The effects of background density inhomogeneity across the magnetic field on the resonant ion heating are studied. The effect of super- Alfvenic ion drift on the ion heating is investigated. It is found that the turbulent wave spectrum of initially parallel propagating waves cascades to oblique modes, and leads to enhanced resonant ion heating due to the inhomogeneity. The acceleration of the solar wind ions is achieved by the parametric instability of large amplitude waves in the spectrum, and is also affected by the inhomogeneity. The results of the study provide the ion temperature anisotropy and drift velocity temporal evolution due to relaxation of the instability. The non-Maxwellian velocity distribution functions (VDFs) of the ions are modeled in the inhomogeneous solar wind plasma in the acceleration region close to the Sun.

Ofman, Leon↗

Dynamics of millicharged dark matter in supernova remnants

Dark matter could have a small electromagnetic charge, provided the charge-to-mass ratio is much less than that of electrons or protons. This candidate, commonly known as millicharged dark matter (mDM), would form a plasma and interact with the interstellar medium and electromagnetic fields within galaxies. In general, understanding the dynamics of mDM requires consideration of collective plasma effects. It has been proposed that mDM can be accelerated in supernova remnants, forming a dark cosmic ray population that would leave distinct experimental signatures. In this work, we study a microphysical model where mDM is shocked by a supernova remnant and isotropized in the frame of the expanding fluid. We find that for | q χ / m χ | ≳ 10 – 13 e / MeV , the isotropization length for electromagnetic plasma instabilities is much shorter than the size of the supernova remnant. This is a necessary, though not sufficient, first step for formation of a Fermi-accelerated mDM component, and determining the size of this component requires further study. Finally, we discuss additional implications of mDM interactions in supernova remnants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗