SEARCH · Search NASA

Results for “plasma global models”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2021

A boundary value “reservoir problem” and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma–wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remain an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1987

Characteristics of cometary picked-up ions in a global model of Giacobini-Zinner

Energetic ions observed during the International Cometary Explorer (ICE) spacecraft flyby of comet Giacobini-Zinner provide information about both the constitution of comets and the plasma physical processes associated with their interaction with the solar wind. In this investigation the details of ion 'pickup,' in the limit where small-scale fluctuations in the plasma and magnetic field are neglected, are modeled by following the motion of a large number of initially cold, heavy (mass 18) ions in a global magnetohydrodynamic model of the local plasma and magnetic field. The results indicate how the background or macroscopic velocity and magnetic field structure of the comet can affect the average spatial and spectral characteristics of the observed cometary ions. These effects, which occur by virtue of forces associated with the compression and the curvature of the magnetic field in the presence of the stagnating plasma flow, can explain the double maxima in the time series of the energetic ion flux observed along the ICE trajectory.

Kimmel, C. D.↗

NASA NTRS1984

Nightward ion flow in the Venus ionosphere - Implications of momentum balance

Using global empirical models of Venus ionospheric conditions, the plasma flow field consistent with the horizontal momentum equation in both viscid and inviscid forms is solved for. It is found that plasma viscosity is negligible except at low altitudes and that the observed plasma flows are consistent with the inviscid solution above 300 km but are larger than the calculated flows at lower altitudes. This is probably due to downward momentum advection. The strong vertical shear in the calculated and observed flows just beyond the terminator may produce turbulence there, manifested as the observed transterminator waves in plasma density and magnetic field. The implications of ionospheric superrotation are discussed.

Elphic, R. C.↗

NASA NTRS1992

Helium escape from the Earth's atmosphere - The charge exchange mechanism revisited

We have studied the escape of neutral helium from the terrestrial atmosphere through exothermic charge exchange reactions between He(+) ions and the major atmospheric constituents N2, O2 and O. Elastic collisions with the neutral background particles were treated quantitatively using a recently developed kinetic theory approach. An interhemispheric plasma transport model was employed to provide a global distribution of He(+) ions as a function of altitude, latitude and local solar time and for different levels of solar ionization. Combining these ion densities with neutral densities from an MSIS model and best estimates for the reaction rate coefficients of the charge exchange reactions, we computed the global distribution of the neutral He escape flux. The escape rates show large diurnal and latitudinal variations, while the global average does not vary by more than a factor of three over a solar cycle. We find that this escape mechanism is potentially important for the overall balance of helium in the Earth's atmosphere. However, more accurate values for the reaction rate coefficients of the charge exchange reactions are required to make a definitive assessment of its importance.

Lie-Svendsen, O.↗

NASA NTRS1984

Application of a global solar wind/planetary obstacle interaction computational model: Earth, Venus, Mars, Jupiter and Saturn studies

This work relates to a series of collaborative investigations involving the application of a computational model for the determination of the detailed plasma and magnetic field properties associated with the global interaction of the solar wind with various planetary obstacles throughout the solar system. The theoretical method is based on an established single fluid, steady, dissipationless, magnetohydrodynamic continuum model, and is appropriate for the calculation of supersonic, super-Alfvenic solar wind flow past planetary obstacles. The investigations undertaken relate to studies of various solar wind interaction phenomena with Venus, Earth, Jupiter, and Saturn.

Stahara, S. S.↗

DOE OSTI2022

An E and B gyrokinetic simulation model for kinetic Alfvén waves in tokamak plasmas

The gyrokinetic particle simulation is a powerful tool for studies of transport, nonlinear phenomenon, and energetic particle physics in tokamak plasmas. While most gyrokinetic simulations make use of the scalar and vector potentials, a new model (GK-E&B) has been developed by using the E and B field in a general form and has been implemented in simulating kinetic Alfvén waves in uniform plasma. In our work, the Chen et al. GK-E&B model has been expressed, in general, tokamak geometry using the local orthogonal coordinates and general tokamak coordinates. Its reduction for uniform plasma is verified, and the numerical results show good agreement with the original work. The theoretical dispersion relation and numerical results in the local model in screw pinch geometry are also in excellent agreement. Numerical results show excellent performance in a realistic parameter regime of burning plasmas with high values of β/(M e k$^{2}_{⊥}$ρ$^{2}_{i}$), which is a challenge for traditional methods due to the “cancellation” problem. As one application, the GK-E&B model is implemented with kinetic electrons in the local single flux surface limit. With the matched International Tokamak Physics Activity-Toroidicity-induced Alfvén Eigenmodes parameters adopted, numerical results show the capability of the GK-E&B in treating the parallel electron Landau damping for realistic tokamak plasma parameters. As another application, the global GK-E&B model has been implemented with the dominant electron contribution in the cold electron limit. Its capability in simulating the finite E || due to the finite electron mass is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Self-organized confinement in whole-device modeling of laboratory magnetospheres

Turbulent self-organization driven by global ideal interchange modes in a dipole-confined plasma is explored with self-consistent, whole-plasma simulations using a flux-tube averaged magnetohydrodynamic model in dipole magnetic geometry. We show the existence of robust particle pinch driven by ideal interchange-mode fluctuations, in which the particles are transported up the density gradient. It is found that the plasma profiles in a dipole field spontaneously relax to a marginally stable state as centrally peaked pressure and density are created by global interchange-mode transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Global gyro-kinetic ion temperature gradient and trapped electron mode turbulence modelling in X -point geometry in negative and positive triangularity

Comparative modelling of ion temperature gradient (ITG)/trapped electron mode (TEM) turbulence in negative (NT) and positive (PT) triangularity plasma shapes was done using the nonlinear global full-f gyrokinetic particle code JOREK-GK in the realistic X-point tokamak geometry including the Scrape Off Layer (SOL) for TCV and DIII-D parameters. A comparison of JOREK-GK code with the gyrokinetic codes GS2 and GENE was done using NT/PT triangularity TCV L-modes parameters showing good agreement between codes in linear growth rates and clear beneficial effect of NT as compared to PT. Global non-linear modelling of the ITG/TEM saturated turbulence for realistic DIII-D NT pulses was done and compared with numerically constructed PT equilibrium with the same plasma profiles. Existence of longer correlation length of density fluctuations in PT compared to NT was demonstrated. Stronger and more sheared zonal flows are generated via Reynolds stress in NT compared to PT. These factors are stabilizing for TEM/ITG turbulence in NT and lead to smaller heat fluxes and heat conductivities in NT compared to PT configuration. Weak dependence of plasma confinement on collisionality and plasma rotation was found in modelling of DIII-D NT shots similar to the experiment. The confinement scaling with normalized ion gyro-radius ρ* was estimated both for NT and PT. Bohm-like scaling was obtained in both configurations, however with better confinement for NT compared to PT which could be favourable factor for reactor size machines.

ITG↗

NASA NTRS1995

Visualization of numerical simulation results: An interactive tool for mission oriented theory

With the advent of multispacecraft missions, it became obvious that theoretical models are needed, not only to interpret local spacecraft observations, but also to link single point measurements. Models that address space plasma phenomena on local and global scales, that involve nonlinear processes, and that use self-consistent approaches, were developed. An interactive system for visualizing and analyzing numerical simulation results was also developed. This system is based on visualization modules which use the application visualization system and the NCAR graphic packages. Examples from recent simulations illustrate how these modules can be used for displaying and manipulating simulation results.

Berchem, Jean↗

DOE OSTI2025

Two Mechanisms Limiting the Emitted Electron Current from a Cathode to an Anode

It is known that the current of emitted electrons flowing through a plasma can saturate upon formation of a potential well adjacent to the cathode (the “space charge effect”). Here, we demonstrate another saturation mechanism that will often set a more restrictive limit on the global current. When “backflow saturation” occurs, the cathode sheath weakens to allow emitted electrons that already entered the plasma to backflow to the cathode. This effect could not be captured by studies modelling the cathode sheath by itself because its origin is coupled to processes in the interior plasma and anode sheath. By modeling a full plasma diode, we show that depending on conditions the global current can be limited in four ways; by backflow alone, by space charge alone, by both mechanisms in a stable cooperative form, or by both in a competing oscillatory form. Published by the American Physical Society 2025

Laboratory plasma↗

NASA NTRS2013

Plasma Ion Composition Measurements for Europa

Jupiter magnetospheric interactions and surface composition, both important to subsurface ocean detection for the Galilean icy moons Europa, Ganymede, and Callisto, can be measured using plasma ion mass spectrometry on either an orbiting spacecraft or one designed for multiple flybys of these moons. Detection of emergent oceanic materials at the Europa surface is more likely than at Ganymede and Callisto. A key challenge is to resolve potential intrinsic Europan materials from the space weathering patina of iogenic species implanted onto the sensible surface by magnetospheric interactions. Species resolved measurements of pickup ion currents are also critical to extraction of oceanic induced magnetic fields from magnetospheric interaction background dominated by these currents. In general the chemical astrobiological potential of Europa should be determined through the combination of surface, ionospheric, and pickup ion composition measurements. The requisite Ion Mass Spectrometer (IMS) for these measurements would need to work in the high radiation environment of Jupiter's magnetosphere between the orbits of Europa and Ganymede, and beyond. A 3D hybrid model of the moon magnetosphere interaction is also needed to construct a global model of the electric and magnetic fields, and the plasma environment, around Europa. Europa's ionosphere is probably usually dominated by hot pickup ions with 100-1000 eV temperatures, excursions to a "classical" cold ionosphere likely being infrequent. A field aligned ionospheric wind driven by the electron polarization electric field should arise and be measurable.

Plasma↗

DOE OSTI2025

Core-edge integrated predictive studies of ST40 and NSTX plasmas with the scrape-off layer box model

The ability to model the interplay between the core and edge of tokamak plasmas is crucial to designing both the plasma operating scenario of a fusion pilot plant and the design of the tokamak itself. Scrape-off-layer (SOL) models that are tailored to integrated scenario modeling need to have fast turn-around time and minimal computational burden to enable wide parameter-space coverage for design scoping. The SOL 0-D Box model is a reduced SOL model based on global power and particle balance that captures the essential physics of SOL transport with little computational cost. The usage of the 0-D Box model in core-edge coupled simulations has been demonstrated in both interpretive and predictive modes on a variety of devices. This paper presents a sensitivity study of the 0-D Box model to the input SOL heat-flux width for an ST40 plasma. This study demonstrates that accurate prediction of this width is crucial to predicting global performance parameters of a plasma scenario, such as energy confinement time and flux consumption. We also present an extension of the Box model to 1-D to allow for parallel variation of plasma parameters along the magnetic field lines. The 1-D Box model is then compared with SOLPS-ITER simulations of an NSTX plasma. Advantages and limitations of the Box model are discussed, and future directions are outlined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Modeled F region response to auroral dynamics based upon Dynamics Explorer auroral observations

Auroral images from the Dynamics Explorer 1 (DE 1) scanning auroral imager have been combined with in situ auroral precipitation data from the DE 2 low-altitude plasma instrument, to form a time-dependent global auroral energy flux model. This model has both good time (12 min) and spatial (100 km) resolution compared to that currently available for global-scale ionospheric and thermospheric modeling. The development and comparison of this model with others are discussed. Data from an aurorally active period, November 25, 1981, are presented and used as a case study for this model. Using a global ionospheric model, the effect of the DE auroral model is contrasted with that of a conventional empirical auroral energy flux model. Major differences in the modeled F region ionosphere are predicted from this comparative study.

Sojka, J. J.↗

NASA NTRS1993

Magnetosheath-ionspheric plasma interactions in the cusp/cleft. 2: Mesoscale particle simulations

Ionospheric plasma flowing out from the cusp can be an important source of plasma to the magnetosphere. One source of free energy that can drive this outflow is the injection of magnetosheath plasma into the cusp. Two-dimensional (three velocity) mesoscale particle simulations are used to investigate the particle dynamics in the cusp during southward interplanetary magnetic field. This mesoscale model self-consistently incorporates (1) global influences such as the convection of plasma across the cusp, the action of the mirror force, and the injection of the magnetosheath plasma, and (2) wave-particle interactions which produce the actual coupling between the magnetosheath and ionospheric plasmas. It is shown that, because the thermal speed of the electrons is higher than the bulk motion of the magnetosheath plasma, an upward current is formed on the equatorward edge of the injection region with return currents on either side. However, the poleward return currents are the stronger due to the convection and mirroring of many of the magnetosheath electrons. The electron distribution in this latter region evolves from upward directed streams to single-sided loss cones or possibly electron conics. The ion distribution also shows a variety of distinct features that are produced by spatial and/or temporal effects associated with varying convection patterns and wave-particle interactions. On the equatorward edge the distribution has a downflowing magnetosheath component and an upflowing cold ionospheric component due to continuous convection of ionospheric plasma into the region. In the center of the magnetosheath region, heating from the development of an ion-ion streaming instability causes the suppression of the cold ionospheric component and the formation of downward ionospheric streams. Further poleward there is velocity filtering of ions with low pitch angles, so that the magnetosheath ions develop a ring-beam distribution and the ensuing wave instabilities generate downward ionospheric conics. These downward ionospheric components are eventually turned by the mirror force, leading to the production of upward conics at elevated energies throughout the region.

Winglee, R. M.↗

NASA NTRS1998

Modeling Magnetotail Ion Distributions with Global Magnetohydrodynamic and Ion Trajectory Calculations

On February 9, 1995, the Comprehensive Plasma Instrumentation (CPI) on the Geotail spacecraft observed a complex, structured ion distribution function near the magnetotail midplane at x approximately -30 R(sub E). On this same day the Wind spacecraft observed a quiet solar wind and an interplanetary magnetic field (IMF) that was northward for more than five hours, and an IMF B(sub y) component with a magnitude comparable to that of the RAF B(sub z) component. In this study, we determined the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time, using the time-dependent electric and magnetic fields obtained from a global MHD simulation. The Wind observations were used as input for the MHD model. The ion distribution function observed by Geotail at 1347 UT was found to consist primarily of particles from the dawn side low latitude boundary layer (LLBL) and from the dusk side LLBL; fewer than 2% of the particles originated in the ionosphere.

El-Alaoui, M.↗

DOE OSTI2020

An Alfvenic reconnecting plasmoid thruster

A new concept for the generation of thrust for space propulsion is introduced. Energetic thrust is generated in the form of plasmoids (confined plasma in closed magnetic loops) when magnetic helicity (linked magnetic field lines) is injected into an annular channel. Using a novel configuration of static electric and magnetic fields, the concept utilizes a current-sheet instability to spontaneously and continuously create plasmoids via magnetic reconnection. Here, the generated low-temperature plasma is simulated in a global annular geometry using the extended magnetohydrodynamic model. Because the system-size plasmoid is an Alfvenic outflow from the reconnection site, its thrust is proportional to the square of the magnetic field strength and does not ideally depend on the mass of the ion species of the plasma. Exhaust velocities in the range of 20 to $500\ \textrm {km}\ \textrm {s}^{-1}$ , controllable by the coil currents, are observed in the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Overview of Liquid-Metal PFC R&D at the University of Illinois Urbana-Champaign

The design and implementation of future flowing liquid-lithium plasma-facing components (LLPFCs) will be dependent on several factors. Of course, one of the most important is the need to be able to deal with high heat fluxes incident on the surface of the LLPFCs, but there are also several other important liquid-metal behaviors that have been identified for their critical impact on the feasibility of a LLPFC. One of these is the ability to constantly wet 100% of the plasma-facing component area and the best way to achieve that. Another key point is knowing and understanding the erosion and corrosion of the surfaces subject to a flowing liquid-lithium system and the ability for hydrogen and helium uptake by the system. The Center for Plasma Material Interactions (CPMI) has been tasked with looking at these various issues. The Mock-up Entry module for EAST device was used to investigate wetting and erosion effects and to design a suitable distribution and collection system with a liquid-lithium loop. The vapor shielding effects of lithium on the surface were also modeled and studied. A model coupling CRANE, an open-source global reaction network solver, and Zapdos, a plasma transport solver, is being developed to better understand the dynamics of the vapor cloud. Experiments on the Magnum-PSI at the Dutch Institute for Fundamental Energy Research have been carried out to study the vapor shielding effect and obtain experimental benchmarks to verify the model. Also, initial experiments using the Hybrid Illinois Device for Research and Applications have been performed to understand the pumping effects of lithium on helium. Experiments with a drop of liquid lithium (~100 mg) into a helium plasma have shown the ability of lithium to take out the cold recycling helium gas as well as hydrogen and oxygen impurity gases. The improvement in plasma performance was significant, and further understanding of this effect will have impacts on how future LLPFCs will be designed. Further investigation into the exact mechanism for helium pumping by lithium needs to be performed in the future. Finally, this paper presents a summary of the results obtained at the CPMI.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗