Search NASA⌕ Search

SEARCH · Search NASA

Results for “radial diffusion”

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

Rapid Enhancements of the Seed Populations in the Heart of the Earth's Outer Radiation Belt: A Multicase Study

To better understand rapid enhancements of the seed populations (100s keV electrons) in the heart of the Earth’s outer radiation belt (L*~3.5–5.0) during different geomagnetic activities, we investigate three enhancement events measured by Van Allen Probes in detail. Observations of the fluxes and the pitch angle distributions of energetic electrons are analyzed to determine rapid enhancements of the seed populations. Our study shows that three specified processes associated with substorm electron injections can lead to rapid enhancements of the seed populations, and the electron energy increases up to 342 keV. In the first process, substorm electron injections accompanied by the transient and intense substorm electric fields can directly lead to rapid enhancements of the seed populations in the heart of the outer radiation belt. In the second process, the substorm injected electrons are first trapped in the outer radiation belt, and subsequently transported into L*<4.5 by the convection electric field. In the third process, the lower-energy electrons are first injected at L*~5.3, and then undergo drift resonance with ultralow frequency (ULF) waves. These accelerated electrons by ULF waves are further transported into L*<4.5 due to the convection electric field. This process is consistent with the radial diffusion. Our results suggest that these specified processes are important for understanding the dynamics of the seed populations in the heart of the outer radiation belt.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

On the role of last closed drift shell dynamics in driving fast losses and Van Allen radiation belt extinction

We present observations of very fast radiation belt loss as resolved using high time resolution electron flux data from the constellation of Global Positioning System (GPS) satellites. The time scale of these losses is revealed to be as short as ~0.5–2 hr during intense magnetic storms, with some storms demonstrating almost total loss on these time scales and which we characterize as radiation belt extinction. The intense March 2013 and March 2015 storms both show such fast extinction, with a rapid recovery, while the September 2014 storm shows fast extinction but no recovery for around 2 weeks. By contrast, the moderate September 2012 storm which generated a three radiation belt morphology shows more gradual loss. Here, we compute the last closed drift shell (LCDS) for each of these four storms and show a very strong correspondence between the LCDS and the loss patterns of trapped electrons in each storm. Most significantly, the location of the LCDS closely mirrors the high time resolution losses observed in GPS flux. The fast losses occur on a time scale shorter than the Van Allen Probes orbital period, are explained by proximity to the LCDS, and progress inward, consistent with outward transport to the LCDS by fast ultralow frequency wave radial diffusion. Expressing the location of the LCDS in L*, and not model magnetopause standoff distance in units of RE, clearly reveals magnetopause shadowing as the cause of the fast loss observed by the GPS satellites.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Observations and Simulations of Dropout Events and Flux Decays in October 2013: Comparing MEO Equatorial With LEO Polar Orbit

We compare ESA PROBA-V observations of electron flux at LEO with those from the NASA Van Allen Probes mostly at MEO for October 2013. Dropouts are visible at all energy during four storms from both satellites. Equatorially trapped electron fluxes are higher than at LEO by 10 2 (<1 MeV) to 10 5 (>2.5 MeV). We observe a quite isotropic structure of the outer belt during quiet times, contrary to the inner belt, and pitch angle dependence of high energy injection. We find a very good overlap of the outer belt at MEO and LEO at ~0.5 MeV. We use test-particle simulations of the energetic electrons trapped in the terrestrial magnetic field to study the outer radiation belt electron flux changes during geomagnetic storms. We show that the Dst (Disturbance storm time) effect during the main phase of a geomagnetic storm results in a betatron mechanism causing outward radial drift and a deceleration of the electrons. This outward drift motion is energy independent, pitch angle-dependent, and represents a significant distance (~1 L-shell at L = 5 for moderate storms). At fixed L-shell, this causes a decay of the LEO precipitating flux (adiabatic outward motion), followed by a return to the normal state (adiabatic inward motion) during main and recovery phases. Dst effect, associated with magnetopause shadowing and radial diffusion can explain the main characteristics of outer radiation belt electron dropouts in October 2013. We also use Fokker-Planck simulations with event-driven diffusion coefficients at high temporal resolution, to distinguish instantaneous loss from the gradual scattering that depopulates the slot region and the outer belt after storms. Simulations reproduce the slot formation and the gradual loss in the outer belt. The typical energy dependence of these losses leads to the absence of scattering for relativistic and ultra-relativistic electrons in the outer belt, oppositely to dropouts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Modeling the Effects of Drift Orbit Bifurcation on Radiation Belt Electrons

Abstract Drift orbit bifurcation (DOB) has been suggested to play a major role in the loss and transport of radiation belt electrons since it violates the second adiabatic invariants of particles and makes the third invariant undefined. Results from our guiding‐center test particle simulations using the Tsyganenko‐1989c magnetic field model show that the DOB could affect a broad region of the outer radiation belt, which can penetrate inside the geosynchronous orbit at Kp ≥ 3, and its effects are more significant further away from the Earth, at a higher Kp and for higher electron energies. Specifically, the short‐term simulation results after one electron drift show both traditional and nontraditional DOB transport of electrons, with the nontraditional DOB, caused by a third minimum of the magnetic field strength near the equator, reported here for the first time. Moreover, our results show large ballistic jumps in the second invariant and radial distance for electrons at high equatorial pitch angles after one drift. In addition, long‐term DOB transport coefficients of electrons over many drifts are calculated based on our simulation results. We find that the pitch angle and radial diffusion coefficients of electrons due to DOB could be comparable to or even larger than those caused by electron interactions with chorus and Ultra‐Low‐Frequency waves, respectively. In sum, our results demonstrate that DOB could cause effective loss and transport of radiation belt electrons even in the absence of waves.

Astronomy & Astrophysics↗
DOE OSTI2023

Modeling Ring Current Proton Fluxes Using Artificial Neural Network and Van Allen Probe Measurements

Abstract Terrestrial ring current dynamics are a critical part of the near‐space environment, in that they directly drive geomagnetic field variations that control particle drifts, and define geomagnetic storms. The present study aims to specify a global and time‐varying distribution of ring current proton using geomagnetic indices and solar wind parameters with their history as input. We train an artificial neural network (ANN) model to reproduce proton fluxes measured by the Radiation Belt Storm Probes Ion Composition Experiment instrument onboard Van Allen Probes. By choosing optimal feature parameters and their history length, the model results show a high correlation and a small error between model specifications and satellite measurements. The modeled results well capture energy‐dependent proton dynamics in association with geomagnetic storms, including inward radial diffusion, acceleration and decay. Our ANN model produces proton fluxes with their corresponding 3D spatiotemporal variations, capturing the latitudinal distribution and local time asymmetry that are consistent with observations and that can further inform theory.

Li, Jinxing↗
DOE OSTI2023

Modeling the Effects of Drift Orbit Bifurcation on the Magnetopause Shadowing Loss of Radiation Belt Electrons

Abstract Last closed drift shell (LCDS) has been identified as a crucial parameter for investigating the magnetopause shadowing loss of radiation belt electrons. However, drift orbit bifurcation (DOB) effects have not been physically incorporated into the LCDS calculation. Here we calculate event‐specific LCDS using different approaches to dealing with the DOB effects, that is, tracing field lines ignoring DOB, tracing test particles rejecting field lines with DOB, and tracing particles including field lines with DOB, and then incorporate them into a radial diffusion model to simulate the fast electron dropout observed by Van Allen Probes in May 2017. The model effectively captures the fast dropout at high L* and exhibits the best agreement with data when LCDS is calculated by tracing test particles with DOB more physically included. This study represents the first quantitative modeling of the DOB effects on radiation belt magnetopause shadowing loss via a more physical specification of LCDS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Modeling the Simultaneous Dropout of Energetic Electrons and Protons by Magnetopause Shadowing

Abstract Magnetopause shadowing (MPS) effect could drive a concurrent dropout of radiation belt electrons and ring current protons. However, its relative role in the dropout of both plasma populations has not been well quantified. In this work, we study the simultaneous dropout of MeV electrons and 100s keV protons during an intense geomagnetic storm in May 2017. A radial diffusion model with an event‐specific last closed drift shell is used to simulate the MPS loss of both populations. The model well captures the fast shadowing loss of both populations at L * > 4.6, while the loss at L * < 4.6, possibly due to the electromagnetic ion cyclotron wave scattering, is not captured. The observed butterfly pitch angle distributions of electron fluxes in the initial loss phase are well reproduced by the model. The initial proton losses at low pitch angles are underestimated, potentially also contributed by other mechanisms such as field line curvature scattering.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Simulating Long–Term Dynamics of Radiation Belt Electrons Using DREAM3D Model

Here, we compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Electrical resistance of the current collector controls lithium morphology

Abstract The electrodeposition of low surface area lithium is critical to successful adoption of lithium metal batteries. Here, we discover the dependence of lithium metal morphology on electrical resistance of substrates, enabling us to design an alternative strategy for controlling lithium morphology and improving electrochemical performance. By modifying the current collector with atomic layer deposited conductive (ZnO, SnO 2 ) and resistive (Al 2 O 3 ) nanofilms, we show that conductive films promote the formation of high surface area lithium deposits, whereas highly resistive films promote the formation of lithium clusters of low surface area. We reveal an electrodeposition mechanism in which radial diffusion of electroactive species is promoted on resistive substrates, resulting in lateral growth of large (150 µm in diameter) planar lithium deposits. Using resistive substrates, similar lithium morphologies are formed in three distinct classes of electrolytes, resulting in up to ten-fold improvement in battery performance. Ultimately, we report anode-free pouch cells using the Al 2 O 3 -modified copper that maintain 60 % of their initial discharge capacity after 100 cycles, displaying the benefits of resistive substrates for controlling lithium electrodeposition.

25 ENERGY STORAGE↗
DOE OSTI2019

A reduced model of neutral-plasma interactions in the edge and scrape-off-layer: Verification comparisons with kinetic Monte Carlo simulations

The 2D scrape-off-layer turbulence code (SOLT) is extended to include neutral-plasma interactions. A Boltzmann equation is derived for the evolution of the bi-normally-averaged neutral distribution function, G( x ,vx,t), in the radial dimension, and this evolution is included in the new code (nSOLT). Neutral-plasma interactions are mediated by charge-exchange (CX) and ionization rates based on poloidally-averaged plasma density and temperature. Good agreement is obtained between asymptotically stationary neutral density profiles from nSOLT simulations and those previously obtained from the Monte Carlo neutral transport code DEGAS 2, for time-averaged NSTX H-mode plasma profiles. The sensitivity of the nSOLT neutral profiles to atomic physics parameters, with and without CX physics, is included in the comparison. In addition, nSOLT simulations that evolve the plasma in 1D, using radial diffusion as a proxy for turbulent (blob) transport, illustrate the convergence to a self-consistent neutral-plasma equilibrium sustained by a neutral source at the far-SOL boundary and plasma heating in the core; equilibria consistent with typical NSTX Ohmic L-mode plasmas are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Towards fast, accurate predictions of RF simulations via data-driven modeling: Forward and lateral models

Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. A single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time to complete, which is acceptable for many purposes, but too slow for integrated modeling and real-time control applications. More accurate simulations with fast electron diffusion are even slower, requiring multiple hours of run time with parallel processing. The machine learning models use a database of 16,000+ GEN-RAY/CQL3D simulations for training, validation, and testing. Latin hypercube sampling methods implemented in πScope ensure that the database covers the range of 9 input parameters (n e0 , T e0 , I p , B t , R 0 , n ∥︀ , Z e f f , V loop , P LHCD ) with sufficient density in all regions of parameter space. The surrogate models reduce the computation time from minutes-hours to ms with high accuracy across the input parameter space. Data-driven surrogate models also allow for solving inverse and “lateral” problems. A surrogate model for the inverse problem maps from a desired current drive or power deposition profile to a set of input parameters that would result in such a profile, while a surrogate model for the lateral problem maps from a measured experimental quantity such as hard x-ray emission to a current drive or power deposition profile. In conclusion, the πScope database creation workflow is flexible and applicable to other RF simulation codes such as TORIC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation of plasma transport in MPS-LD linear plasma device by using BOUT++

Here a linear plasma device (LPD) module based on the 2D transport module under the BOUT++ framework is developed in this paper to simulate plasma transport in the LPD. The LPD module includes three parts, i.e. magnetic field calculation, simulation mesh generation and plasma transport setup. The magnetic field is calculated with the circular current loop using the location and current information of each coil. The mesh generation code can produce a simulation mesh for LPD by employing the magnetic field. The plasma transport model is based on the reduced Braginskii equations, which consist of the continuity equation, momentum equation, and energy equation. The fluid neutral model is applied for neutral particle simulation. Deuterium (D) and helium (He) ions and atoms can be simulated by the model. The first attempt to simulate plasma transport in the new LPD, called multiple plasma simulation linear plasma device (MPS-LD), is presented using the developed model. The D plasma transport in the MPS-LD is simulated and benchmarked against the two-point model, showing the validation of the BOUT++ simulation. The effects of radial transport on the heat load to the target are studied, which illustrates the significant impact of the radial diffusivities $D_⟂$ and $\mathcal{X}_{⟂i,e}$ on the plasma. Moreover, the He impurity injection process during the discharge is studied with emphasis on plasma–impurity interactions. The simulation results show that He injection can reduce the plasma energy load to the target significantly, and the efficiency depends on the He source density and injection velocity. The present work provides an alternative and flexible simulation tool for plasma transport in LPDs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Modeling turbulent impurity transport in the SOL of DIII-D with a reduced model

A novel impurity transport model that approximates SOL turbulence as a fluctuating poloidal electric field is shown to be an acceptable replacement for the traditional approach of assigning an arbitrary radial diffusion coefficient to the impurity ions. The model is implemented in the DIVIMP impurity transport code and applied to an L-Mode tungsten divertor experiment on DIII-D. The poloidal electric field is represented as fluctuating between ±1000 V m –1 based on previous measurements. The resulting intermittent v r = E θ × B T transport causes ions to transport both into the core as well as into the far-SOL. Simultaneous agreement with estimates of the W density just inside the separatrix as well as in the far-SOL is obtained (n W ~ 10 14 m –3 and n W ~ 10 12 m –3 , respectively). Prompt re-deposition of the W ions was necessary to obtain agreement (f redep ~ 99%). Here, we conclude that simulating impurity transport using a physics-based approximation for turbulence in the SOL, versus arbitrarily assigning diffusion coefficients, may enable better reactor scale predictions of core impurity contamination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Heating and ablation of high- Z cryogenic pellets in high temperature plasmas

A kinetic model is developed here to determine the power deposition from energetic electrons into the neutral gas shield of an ablating high-Z pellet. For high-Z, the velocity distribution of the hot electrons is nearly isotropic, and we use this feature to develop solutions to the kinetic equation. In contrast to pre-existing models, we consider the effect of gyro-motion as well as elastic scattering of the hot electrons. Two limits are considered: when the gyro-frequency is much greater than the elastic collision frequency, the hot electrons diffuse longitudinally along the field lines as they slow down; but if the gyro-frequency is much less than the elastic collision frequency, the hot electrons diffuse radially. In both limits, it is possible to express the kinetic equation describing the hot electrons independently of the gas density profile, and therefore, the power deposition model is universal in this respect. The emphasis on elastic scattering yields an ablation rate which scales as Z -7/6 which is different than Z -2/3 shown in Sergeev et al (2006 Plasma Phys. Rep. 32 5). It is also demonstrated that the sheath potential required to maintain ambipolarity in the cloud scales as Z -1/3 . Fluid simulations yield ablation rates that scale with the four-thirds power of the pellet radius in agreement with [2].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The role of B T -dependent flows on W accumulation at the edge of the confined plasma

Abstract Near-separatrix impurity accumulation between the crown and the outer midplane of tokamaks is a common feature in results from codes such as SOLPS-ITER and DIVIMP; however, experimental evidence of accumulation has only recently been obtained and is reported here. The codes find that the poloidal distribution of impurity ions in the scrape-off layer (SOL) depends primarily on toroidal field ( B T )-dependent parallel flow patterns of the background plasma and the parallel ion temperature gradient (∇ ‖ T ion ) force. Experimentally, Mach probes used in L-mode plasmas with favorable (for H-mode access) B T measure fast ( M ∼ 0.3–0.5) inner-target-directed (ITD) background plasma flows at the crown of single-null discharges. This study reports a set of DIVIMP simulations for two similar H-mode discharges from the DIII-D W metal rings campaign differing primarily in B T -direction to assess the effect that fast ITD flows have on the distribution of W ions in the SOL. It is found that for imposed ITD flows of M = 0.3, W ions that otherwise accumulate due to the ∇ ‖ T ion -force are largely flushed out. It is also found that doubling the radial diffusion coefficient from 0.3 to 0.6 m 2 s −1 prevents accumulation due to rapid cross-field transport into the far-SOL, where background plasma flows drain W ions to the divertors. Far-SOL W distributions from DIVIMP are then used to specify input to the impurity transport code 3DLIM, which is used to interpretively model collector probe (CP) deposition patterns measured in the ‘wall-SOL’. It is demonstrated that the deposition patterns are consistent with the DIVIMP predictions of near-SOL accumulation for the unfavorable- B T direction, and little/no accumulation for the favorable- B T direction. The wall-SOL CPs have thus provided the first experimental evidence, albeit indirect, of near-SOL W accumulation—finding it occurs for the unfavorable- B T direction only. For the favorable- B T direction, fast flows can largely prevent accumulation from occurring.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Fusion RF Modeling Machine Learning (FusionML_RF) v1.0

FusionML_RF consists of multiple codes and trained machine learning (ML) models that perform low-cost output modeling from the Genray-CQL3D. Three machine learning techniques (multilayer perceptron, random forest, and Gaussian process) provide fast surrogate models for lower hybrid current drive (LHCD) simulations. For example, completing a single GENRAY/CQL3D simulation without radial diffusion of fast electrons requires several minutes of wall-clock time. On the other hand, these ML models achieve ~ms of inference time with high accuracy across the input parameter space. This software collection consists of multiple components. (1) codes that use ML methods and precomputed Genray-CQL3D simulation output to build regression models that enable approximate computations of Genray-CLQ3D outputs from arbitrary but physically meaningful input parameters (surrogate modeling); (2) three trained models created by the team, using a database of 16,000+ GENRAY/CQL3D simulations, to study the performance of ML models for surrogate modeling; (3) codes that load the trained models and simulation data, and then compute mean squared error between the models' predictions and the ground truth of simulation output data. This collection is being made available in conjunction with a scientific publication about the work to promote reusability and provide an artifact of the scientific work.

Bai, Zhe↗
DOE OSTI2017

Los Alamos Space Weather Summer School: Institutional Computing 2016

During the summer school, students carry out independent research projects on a range of topics related to space weather. In 2016, one student used the LANL Institutional Computing resources. Results of this project were the first to demonstrate that the magnitude of radial diffusion is found to agree well with the early observations of the Earth's radiation belts, indicating this effect should be included in community models of the radiation belts.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗