Search NASA⌕ Search

SEARCH · Search NASA

Results for “boundary plasma transport”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

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↗
NASA NTRS1992

A two satellite study of nightside flux transfer events in the plasma sheet

The short-term (about 1 min) events during two consecutive substorms on March 23 1979 have been studied using high time resolution measurements of magnetic and electric fields and of the kiloelectronvolt electron flux at X about -21 R(e). Data obtained make it possible to describe the localized magnetic structures transported or propagated within the plasma sheet and to identify these events as plasma sheet nightside magnetic flux transfer events (NFTEs). The best examples of the NFTEs found so far are all observed in a thin portion of plasma sheet where the contraction of the reconnected field lines dominates over the oppositely directed pressure gradient force. Systematic differences in the magnetic variations observed at ISEE-1 and 2 indicate current concentration at the outer plasma sheet boundary during the passage of a NFTE.

Sergeev, V. A.↗
NASA NTRS1998

Sources and Transport of Plasma Sheet Ions During Magnetospheric Substorms

This study investigates the sources and transport of ions observed in the near-Earth plasma sheet during the growth and expansion phases of a magnetospheric substorm that took place on November 24, 1996. The sources and acceleration mechanisms of ions observed at Geotail were determined by calculating the trajectories of thousands of ions backward in time. We found that during the growth phase of the substorm, most of the ions reaching Geotail had origins in the low latitude boundary layer (LLBL) and were already in the magnetosphere when the growth phase began. Late in the growth phase and in the expansion phase more plasma mantle ions reached the Geotail location. Indeed, during the expansion phase more than 90% of the ions were from the mantle. The ions were accelerated enroute to the spacecraft; however, most of the energy gained was achieved by non-adiabatic acceleration during the ions' crossing of the equatorial current sheet just prior to the detection of the ions.

Ashour-Abdalla, M.↗
NASA NTRS1986

DE-1 observations of solar wind-magnetosphere coupling processes in the polar cusp

Measurements of hot plasma distribution functions in the midaltitude polar cusp by the Dynamics Explorer-1 spacecraft have provided new information on the processes of magnetosheath plasma injection and transport and the dependence of these phenomena on solar-wind parameters. Results on magnetosheath plasma injection have included: (1) localization of the region of plasma injection to the high-latitude cusp-magnetosheath interface, (2) movement of this region toward lower latitudes as B(z) becomes more southward, (3) maintenance of quasi-neutrality of suprathermal plasmas within the cusp, and (4) support for merging-type injection at the high-latitude cusp boundary for northward IMF. Plasma flow measurements from Dynamics Explorer-1 have led to a global IMF-dependent convection model which contains both merging and viscous driving forces.

Burch, J. L.↗
NASA NTRS1987

Core plasma in the magnetosphere

The paper presents a review of new findings related to core plasma (which includes ions and electrons with energies less than 50 eV) identified in studies of magnetosphere during 1983-1986. These findings include the cleft ion fountain; the plasmaspheric heavy ion torus; core molecular ions; high core-plasma densities in the plasma-sheet boundary layer; intense transverse heating of core plasmas in the equatorial regions; the supersonic polar wind (H/+/ and He/+/); toroidal or ring heavy-ion distribution functions in the auroral region; conic-to-field-aligned ion signatures of inverted-V auroral events; and the large-scale electron density structure of the magnetosphere. The progress made on modeling core plasma includes simulation of thermal helium heating at the equator; kinetic models of the plasma transport from localized ionospheric regions; time-dependent hydrodynamic models of the polar wind and plasmasphere-ionosphere coupling; and kinetic models of plasmasphere refilling.

Horwitz, J. L.↗
DOE OSTI2020

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
NASA NTRS1993

On nonsinusoidal waves at the Earth's magnetopause

We have examined International Sun Earth Explorer (ISEE) 1 and 2 data during an interval of multiple magnetopause crossings along the flanks of the tail on January 22, 1978. During the event the interplanetary magnetic field was mainly northward. We found that the surface waves on the magnetopause were nonsinusoidal with steepened sunward facing surfaces. This result is consistent with observations reported for an analogous event by Chen et al. (1993). The dawn-dusk electric field E in the low-latitude boundary layer was approx. 10 kV/R(sub E), which is much greater than the approx. 1 kV/R(sub E) dawn-dusk E field typical of the sunward convecting plasma sheet plasma. The ratio of the dawn-dusk E field in the LLBL to the dusk-dawn E in the magnetosheath was approx. 1/3. We propose that the magnetic field in the magnetosheath modulates the wave form in a way that may result in an anomalous transport of momentum across the magnetopause.

Chen, S.-H.↗
DOE OSTI2020

Improved core-edge compatibility using impurity seeding in the small angle slot (SAS) divertor at DIII-D

Impurity seeding studies in the slot divertor (SAS) at DIII-D have revealed a strong relationship between the detachment onset and pedestal characteristics with both target geometry and impurity species. N 2 seeding in the slot has led to the first simultaneous observation of detachment on the entire suite of boundary diagnostics viewing the SAS without degradation of core confinement. SOLPS-ITER simulations with D+C+N, full cross-field drifts and n-n collisions activated are performed for the first time in DIII-D to interpret the behavior. This illustrates a strong effect of divertor configuration and plasma drifts on the recycling source distribution with significant consequences on plasma flows. Flow reversal is found for both main ions and impurities affecting strongly the impurity transport and providing an explanation for the observed dependence on strike point location of the detachment onset and impurity leakage found in the experiments. Matched discharges with either nitrogen or neon injection show that while nitrogen does not significantly affect the pedestal, neon leads to increased pedestal pressure gradient and improved pedestal stability. Little nitrogen penetrates in the core, but a great amount of neon is found in the pedestal consistent with the different ionization potentials of the two impurities. This research demonstrates that neutral and impurity distributions in the divertor can be controlled through variations in strike point locations in a fixed baffle structure. Divertor geometry combined with impurity seeding enables mitigated divertor heat flux balancing core contamination and thus leading to enhanced divertor dissipation and improved core-edge compatibility which are essential for ITER and for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Validation of IMEP on Alcator C-Mod and JET-ILW ELMy H-mode plasmas

Abstract The recently developed integrated model based on engineering parameters (IMEP) (Luda et al 2020 Nucl. Fusion 61 126048; Luda et al 2021 Nucl. Fusion 60 036023), so far validated on ASDEX Upgrade, has been tested on a database of 3 Alcator C-Mod and 55 JET-ILW ELMy (type I) H-mode stationary phases. The empirical pedestal transport model included in IMEP, consisting now of imposing a fixed value of R < ∇ T e > / T e , t o p = − 82.5 , allows an accurate prediction of the pedestal top temperature (when the pedestal top density is fixed to the experimental measurements) across these three machines with different sizes, when the pedestal is peeling–ballooning (PB) limited. Cases far from the ideal PB boundary, corresponding to high edge Spitzer resistivity, are instead strongly overpredicted by IMEP. A comparison between the predictions of Europed and IMEP for a subset of JET-ILW cases shows that IMEP can more accurately reproduce the experimental pedestal width. This allows IMEP to better capture profile effects on the pedestal stability, and therefore to correctly describe the negative effect of fueling on the pedestal pressure for PB limited cases. A strong correlation between the separatrix density and the fueling rate has been identified for a subset of JET-ILW cases, when taking into account different divertor configurations. Overall, these promising results encourage further developments of integrated models to obtain reliable predictions of pedestal and global confinement using only engineering parameters for present and future machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Dawn-dusk asymmetries in the low-latitude boundary layer arising from the Kelvin-Helmholtz instability: A particle simulation

Along the low-latitude boundary layer (LLBL), the Kelvin-Helmholtz (K-H) instability can provide a means for injection of solar wind plasma across closed field lines of the magnetosphere. A fully electromagnetic, two-dimensional (three-velocity) particle code is used to investigate dawn-dusk asymmetries that can arise from velocity differences due to gradient drifts and electric field gradients, and ion acceleration across the narrow field transition layers in the dawn and dusk flanks. The model includes the dawn and dusk sides of the LLBL simultaneously in a slab geometry, incorporating seperate populations of ions and electrons for the magnetosphere and the magnetosheath plasmas. We report several effects: (1) asymmetries in the observed morphology of the turbulent structures, with familiar fluidlike vortex formation on the duskside, but tongues of magnetosheath plasma penetrating into the magnetosphere on the dawnside; (2) the formation of discrete current layers, characterized by strong currents and sharp gradients in the magnetic field, and discrete charge layers, having net charges and constant, weaker currents; (3) increasing asymmetry in dawn/dusk behavior with a decrease of initial currents; (3) increasing asymmetry in dawn/dusk behavior with a decrease of initial boundary layer width; and (4) enhancement of the dawn-to-dusk electric field as magnetosheath particles and momenta are transported across the magnetopause.

Wilber, M.↗
NASA NTRS1999

Origins and Transport of Ions during Magnetospheric Substorms

We investigate the origins and the transport of ions observed in the near-Earth plasma sheet during the growth and expansion phases of a magnetospheric substorm that occurred on November 24, 1996. Ions observed at Geotail were traced backward in time in time-dependent magnetic and electric fields to determine their origins and the acceleration mechanisms responsible for their energization. Results from this investigation indicate that, during the growth phase of the substorm, most of the ions reaching Geotail had origins in the low latitude boundary layer (LLBL) and had alread@, entered the magnetosphere when the growth phase began. Late in the growth phase and in the expansion phase a higher proportion of the ions reaching Geotail had their origin in the plasma mantle. Indeed, during the expansion phase more than 90% of the ions seen by Geotail were from the mantle. The ions were accelerated enroute to the spacecraft; however, most of the ions' energy gain was achieved by non-adiabatic acceleration while crossing the equatorial current sheet just prior to their detection by Geotail. In general, the plasma mantle from both southern and northern hemispheres supplied non-adiabatic ions to Geotail, whereas the LLBL supplied mostly adiabatic ions to the distributions measured by the spacecraft.

Ashour-Abdalla, Maha↗
NASA NTRS1999

The Origin of the Near-Earth Plasma Population During a Substorm on November 24, 1996

We investigate the origins and the transport of ions observed in the near-Earth plasma sheet during the growth and expansion phases of a magnetospheric substorm that occurred on November 24, 1996. Ions observed at Geotail were traced backward in time in time-dependent magnetic and electric fields to determine their origins and the acceleration mechanisms responsible for their energization. Results from this investigation indicate that, during the growth phase of the substorm, most of the ions reaching Geotail had origins in the low latitude boundary layer (LLBL) and had already entered the magnetosphere when the growth phase began. Late in the growth phase and in the expansion phase a higher proportion of the ions reaching Geotail had their origin in the plasma mantle. Indeed, during the expansion phase more than 90% of the ions seen by Geotail were from the mantle. The ions were accelerated enroute to the spacecraft; however, most of the ions' energy gain was achieved by non-adiabatic acceleration while crossing the equatorial current sheet just prior to their detection by Geotail. In general, the plasma mantle from both southern and northern hemispheres supplied non-adiabatic ions to Geotail, whereas the LLBL supplied mostly adiabatic ions to the distributions measured by the spacecraft. Distribution functions computed at the ion sources indicate that ionospheric ions reaching Geotail during the expansion phase were significantly heated. Plasma mantle source distributions indicated the presence of a high-latitude reconnection region that allowed ion entry into the magnetosphere when the IMF was northward. These ions reached Geotail during the expansion phase. Ions from the traditional plasma mantle had access to the spacecraft throughout the substorm.

Ashour-Abdalla, M.↗
NASA NTRS2007

Virtual Energetic Particle Observatory for the Heliospheric Data Environment

The heliosphere is pervaded by interplanetary energetic particles, traditionally also called cosmic rays, from solar, internal heliospheric, and galactic sources. The particles species of interest to heliophysics extend from plasma energies to the GeV energies of galactic cosmic rays still measurably affected by heliospheric modulation and the still higher energies contributing to atmospheric ionization. The NASA and international Heliospheric Network of operational and legacy spacecraft measures interplanetary fluxes of these particles. Spatial coverage extends from the inner heliosphere and geospace to the heliosheath boundary region now being traversed by Voyager 1 and soon by Voyager 2. Science objectives include investigation of solar flare and coronal mass ejection events, acceleration and transport of interplanetary particles within the inner heliosphere, cosmic ray interactions with planetary surfaces and atmospheres, sources of suprathermal and anomalous cosmic ray ions in the outer heliosphere, and solar cycle modulation of galactic cosmic rays. The Virtual Energetic Particle Observatory (VEPO) will improve access and usability of selected spacecraft and sub-orbital NASA heliospheric energetic particle data sets as a newly approved effort within the evolving heliophysics virtual observatory environment. In this presentation, we will describe current VEPO science requirements, our initial priorities and an overview of our strategy to implement VEPO rapidly and at minimal cost by working within the high-level framework of the Virtual Heliospheric Observatory (VHO). VEPO will also leverage existing data services of NASA's Space Physics Data Facility and other existing capabilities of the U.S. and international heliospheric research communities.

Cooper, J. F.↗
NASA NTRS2011

Modeling the Superstorm in November 2003

The superstorm on 20.21 November 2003 was the largest geomagnetic storm in solar cycle 23 as measured by Dst, which attained a minimum value of .422 nT. We have simulated this storm to understand how particles originating from the solar wind and ionosphere get access to the magnetosphere and how the subsequent transport and energization processes contribute to the buildup of the ring current. The global electromagnetic configuration and the solar wind H+ distribution are specified by the Lyon-Fedder-Mobarry (LFM) magnetohydrodynamics model. The outflow of H+ and O+ ions from the ionosphere are also considered. Their trajectories in the magnetosphere are followed by a test-particle code. The particle distributions at the inner plasma sheet established by the LFM model and test-particle calculations are then used as boundary conditions for a ring current model. Our simulations reproduce the rapid decrease of Dst during the storm main phase and the fast initial phase of recovery. Shielding in the inner magnetosphere is established at early main phase. This shielding field lasts several hours and then breaks down at late main phase. At the peak of the storm, strong penetration of ions earthward to L shell of 1.5 is revealed in the simulation. It is surprising that O+ is significant but not the dominant species in the ring current in our calculation for this major storm. It is very likely that substorm effects are not well represented in the models and O+ energization is underestimated. Ring current simulation with O+ energy density at the boundary set comparable to Geotail observations produces excellent agreement with the observed symH. As expected in superstorms, ring current O+ is the dominant species over H+ during the main to mid-recovery phase of the storm.

Fok, Mei-Ching↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2012

The Impact of Geometrical Constraints on Collisionless Magnetic Reconnection

One of the most often cited features associated with collisionless magnetic reconnection is a Hall-type magnetic field, which leads, in antiparallel geometries, to a quadrupolar magnetic field signature. The combination of this out of plane magnetic field with the reconnection in-plane magnetic field leads to angling of magnetic flux tubes out of the plane defined by the incoming magnetic flux. Because it is propagated by Whistler waves, the quadrupolar field can extend over large distances in relatively short amounts of time - in fact, it will extend to the boundary of any modeling domain. In reality, however, the surrounding plasma and magnetic field geometry, defined, for example, by the overall solar wind flow, will in practice limit the extend over which a flux tube can be angled out of the main plain. This poses the question to what extent geometric constraints limit or control the reconnection process and this is the question investigated in this presentation. The investigation will involve a comparison of calculations, where open boundary conditions are set up to mimic either free or constrained geometries. We will compare momentum transport, the geometry of the reconnection regions, and the acceleration if ions and electrons to provide the current sheet in the outflow jet.

Hesse, Michael↗
DOE OSTI2022

Staircase formation by resonant and non-resonant transport of potential vorticity

The E × B staircase is a quasi-periodic pattern of pressure profile corrugations. In this work, we present a new mechanism for E × B staircase formation that involves resonant transport versus non-resonant transport. We start from a potential vorticity evolution system and use quasi-linear theory, a model dispersion relation, and a bi-Lorentzian spectrum approximation, to construct the relation between the fluxes and the profiles. With these fluxes, we close the profile evolution equations and the extended turbulence intensity evolution equation, which together constitute a turbulence-profile evolution system. In this system, the Doppler effect from the E × B mean flow can cause resonance between trapped ion precession motion and the trapped ion mode, which drives a resonant transport contribution to the fluxes. The profiles will be flattened where the resonant transport is switched on. In contrast, for the regions of non-resonant transport, profiles are steeper. A quasi-periodic pattern of profile corrugation (the E × B staircase) spontaneously emerges in this system, which is the two states mentioned above, arranged as alternating layers in space. The feedback processes during the staircase pattern formation are identified. An estimate of the critical value of the boundary heat flux is obtained, above which the staircase formation will be triggered. An estimate scaling of the step size in the staircase pattern is obtained. The resonant turbulent transport is also a mechanism for collisionless saturation of zonal flow. Furthermore, this work is related to internal transport barrier formation and suggests some new scenarios, such as an enhanced confined L mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗