Data-driven modeling of Landau damping by physics-informed neural networks
Not Available
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
SEARCH · Engineering Papers
Results for “space plasma”
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.
The SpacePy space science package at 12 years
For over a decade, the SpacePy project has contributed open-source solutions for the production and analysis of heliophysics data and simulation results. Here we introduce SpacePy’s functionality for the scientific user and present relevant design principles. We examine recent advances and the future of SpacePy in the broader scientific Python ecosystem, concluding with some of the work that has used SpacePy.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Modeling Electron Acceleration and Transport in the Early Impulsive Phase of the 2017 September 10th Solar Flare
The X8.2-class limb flare on 2017 September 10 is among the best studied solar flare events owing to its great similarity to the standard flare model and the broad coverage by multiple spacecraft and ground-based observations. These multiwavelength observations indicate that electron acceleration and transport are efficient in the reconnection and flare looptop regions. However, there lacks a comprehensive model for explaining and interpreting the multi-faceted observations. In this work, we model the electron acceleration and transport in the early impulsive phase of this flare. We solve the Parker transport equation that includes the primary acceleration mechanism during magnetic reconnection in the large-scale flare region modeled by MHD simulations. We find that electrons are accelerated up to several MeV and fill a large volume of the reconnection region, similar to the observations shown in microwaves. The electron spatial distribution and spectral shape in the looptop region agree well with those derived from the microwave and hard X-ray emissions before magnetic islands grow large and dominate the acceleration. Future emission modelings using the electron maps will enable direct comparison with microwave and hard X-ray observations. These results shed new light on the electron acceleration and transport in a broad region of solar flares within a data-constrained realistic flare geometry.
79 ASTRONOMY AND ASTROPHYSICS↗
Anisotropic Energy Transfer and Conversion in Magnetized Compressible Turbulence
We present a spatial filtering (or coarse-graining) analysis on 3D magnetized magnetohydrodynamic (MHD) turbulence simulations. The filtered compressible MHD formulae show transfer of kinetic and magnetic energies from large to small scales, as well as energy conversion between kinetic, magnetic, and thermal energies. The anisotropic filtering enables separate analyses of the energy flows perpendicular and parallel to the global mean magnetic field. Anisotropy in energy cascade is demonstrated by the larger perpendicular energy cascade rate and also the larger perpendicular wavenumbers associated with the peak energy transfer rate. We also find that the “inertial range” along the parallel (perpendicular) direction in the anisotropic energy cascade formulation is no longer strictly dissipation-free, because it includes the dissipation in the perpendicular (parallel) direction. A change in the driving force (kinetic only versus kinetic and magnetic) affects the energy conversion between kinetic and magnetic energies. While the compressibility of the driving force changes the partition of different channels of energy transfer and conversion, and also increases the total energy transfer rate, the global energy flow remains unaffected by compressibility qualitatively. Our analysis can be applied to multispacecraft observations of turbulence in the solar wind or a planetary magnetosphere.
79 ASTRONOMY AND ASTROPHYSICS↗
A Study of Slow-mode Shocks in the Near-Earth Magnetotail with MMS Observations and Hybrid Simulations
The structure of the magnetic reconnection boundary, particularly the presence of slow-mode shocks in the near-Earth magnetotail was studied by using magnetospheric multiscale (MMS) observations and 2.5D hybrid simulations. A total of 51 crossings of MMS from 2017 to 2021 were analyzed. We found that the detection percentage of slow-mode shocks in the near-Earth magnetotail is 41%–55%. Previous studies have only reported one slow-mode shock event in the near-Earth magnetotail and a slow-mode shock detection percentage of 10% or lower in the mid-to-distant magnetotail. It was observed that if the high-energy beam region data is removed from the slow-mode shock downstream observations then the detection of slow-mode shocks reduces, implying that the kinetic effects play an important role in the detection of slow-mode shocks. For the crossings where the interface was not identified as a slow-mode shock, it was found that the turbulence in those crossings can change the mass flux values and disrupt the detection of slow-mode shock. However, the macroscopic slow-mode shock-like structure stably exists around the magnetic reconnection interface, as most of the conditions for slow-mode shocks were satisfied. This result suggests that slow-mode shocks are a general feature of magnetic reconnection geometry. We find that the lack of detection of slow-mode shocks in previous observations and simulations can be explained by taking into account the kinetic structure of slow-mode shocks and the presence of turbulence.
79 ASTRONOMY AND ASTROPHYSICS↗
A Numerical Model for the Dynamics of Pickup Ions Outside the Heliopause and IBEX “Ribbon” Observation
The leading mechanism for the origin of the Interstellar Boundary Explorer (IBEX) “ribbon” of enhanced energetic neutral atoms (ENAs) from the outer heliosphere is the so-called secondary ENA process. In this mechanism, neutralized ions escape the heliosphere, and after several times of charge exchange collision,s some propagate back toward Earth in directions nearly perpendicular to the local interstellar magnetic field. However, the physical processes governing the distribution of the pickup ions (PUIs) outside the heliopause are still unclear. In this study, we build a new global model of the IBEX ribbon, where a key component is to calculate the dynamics of the PUIs outside the heliopause by solving the gyrophase-averaged focused transport equations on top of the background magnetohydrodynamic-kinetic model for the heliosphere-interstellar medium interaction. We discuss how the properties of the simulated ribbon change with different scattering parameters and show simulation results of the ENA sky maps, ribbon centers, ribbon profiles at several azimuthal slices, the ribbon ENA source region, and the velocity distribution of PUIs in the ENA source region for different scattering frequencies. Our results show that the model can provide reasonable ribbon sky maps comparable to the observed ribbon when the PUI mean free path is on the order of several thousand astronomical units.
79 ASTRONOMY AND ASTROPHYSICS↗
Energy Conversion and Electron Acceleration and Transport in 3D Simulations of Solar Flares
Recent observations and simulations indicate that solar flares undergo extremely complex 3D evolution, making 3D particle transport models essential for understanding electron acceleration and interpreting flare emissions. In this study, we investigate this problem by solving Parker’s transport equation with 3D MHD simulations of solar flares. By examining energy conversion in the 3D system, we evaluate the roles of different acceleration mechanisms, including reconnection current sheet (CS), termination shock (TS), and supra-arcade downflows (SADs). We find that large-amplitude turbulent fluctuations are generated and sustained in the 3D system. The model results demonstrate that a significant number of electrons are accelerated to hundreds of keV and even a few MeV, forming power-law energy spectra. These energetic particles are widely distributed, with concentrations at the TS and in the flare looptop region, consistent with results derived from recent hard X-ray (HXR) and microwave (MW) observations. By selectively turning particle acceleration on or off in specific regions, we find that the CS and SADs effectively accelerate electrons to several hundred keV, while the TS enables further acceleration to MeV. However, no single mechanism can independently account for the significant number of energetic electrons observed. Instead, the mechanisms work synergistically to produce a large population of accelerated electrons. Our model provides spatially and temporally resolved electron distributions in the whole flare region and at the flare footpoints, enabling synthetic HXR and MW emission modeling for comparison with observations. These results offer important insights into electron acceleration and transport in 3D solar flare regions.
79 ASTRONOMY AND ASTROPHYSICS↗
Thermodynamics of the Inner Heliosheath
We derive annual sky maps of the proton temperature in the inner heliosheath (IHS), and track their temporal evolution over the years 2009–2016 of Interstellar Boundary Explorer observations. Other associated thermodynamic parameters also determined are the density, kappa (the parameter that characterizes kappa distributions), temperature rate, polytropic index, and entropy. We exploit the theory of kappa distributions and their connection with polytropes, to (i) express a new polytropic quantity Π that remains invariant along streamlines where temperature and density may vary, (ii) parameterize the proton flux in terms of the Π invariant and kappa, and (iii) derive the temperature and density, respectively, from the slope and intercept of the linear relationship between kappa and logarithm of Π. We find the following thermodynamic characteristics: (1) temperature sky maps and histograms shifted to their lowest values in 2012 and their highest in 2015; (2) temperature negatively correlated with density, reflecting the subisothermal polytropic behavior; (3) temperature positively correlated with kappa, revealing characteristics of the mechanism responsible for generating kappa distributions; (4) processes in IHS are subisothermal tending toward isobaric, consistent with previously published results; (5) linear relationship between kappa and polytropic indices, revealing characteristics of the particle potential energy; and (6) entropy positively correlated with polytropic index, aligned with the underlying theory that entropy increases toward the isothermal state where the kappa distribution reduces to the Maxwell–Boltzmann description.
79 ASTRONOMY AND ASTROPHYSICS↗
Turbulent Acceleration of Interstellar Pickup Ions at the Heliospheric Termination Shock Forms the Global ENA Spectrum
The heliospheric energetic neutral atom spectrum observed by the Interstellar Boundary Explorer (IBEX) reveals that the heliosheath proton distribution is consistent with a power law. The origin of the spectrum is likely from interstellar pickup ions (PUIs) accelerated at the heliospheric termination shock (HTS). We present an explanation of the proton spectrum origin using a test particle simulation of PUIs accelerated at the HTS. PUIs experience preferential heating by the motional electric field in the shock foot, but do not develop a power-law tail without the presence of turbulence at wavenumbers (k) close to the PUI gyroradius scale (R g ). Voyager 2 observations of the magnetic field downstream of the HTS indicate a moderate amount of turbulence at kR g ≅ 1, (δB /B 0 ) 2 @ 0.01, which we find to be sufficient for producing a downstream suprathermal PUI tail but not at intensities observed by IBEX. Within the shock ramp, however, Voyager observed the turbulence power at much smaller scales to be nearly 100 times stronger, suggesting the possibility of strong turbulence at the PUI gyroradius scale. We show that a proton distribution can develop a power law downstream of the HTS consistent with IBEX observations if (δB /B 0 ) 2 ≳ 0.1 at kR g ≅ 1 in the shock foot. Shock drift acceleration of PUIs by the motional electric field is aided by interactions with turbulence upstream of the shock overshoot. Steepening of the IBEX proton spectrum in directions farther from the heliospheric nose suggests the HTS compression ratio and/or turbulence power weakens near the heliotail.
79 ASTRONOMY AND ASTROPHYSICS↗
Collisionless Larmor coupling and blob formation in a laser-plasma expanding into a magnetized ambient plasma
Collisionless Larmor coupling is a fundamental process in space and astrophysical plasmas that enables momentum transfer between an expanding plasma and a magnetized ambient medium. Here, in this paper, we report on the laboratory experimental study of Larmor coupling leading to the formation of a plasma blob associated with a laser-driven, super-Alfvénic plasma flow on the Large Plasma Device at the University of California, Los Angeles. The high-repetition rate enables systematic spatial and temporal scans of the plasma evolution using Doppler spectroscopy, as well as measurements of the magnetic field, electrostatic field, and self-emission of both debris and ambient ions using filtered imaging. We observe the self-focusing of the laser-produced plasma and the formation of a secondary diamagnetic cavity associated with a blob composed of background ions. Doppler spectroscopy reveals the transverse velocity distribution of the background ions, providing direct evidence of ion energization via Larmor coupling. The systematic spatial and temporal scans enabled by the high-repetition-rate experiment allow for a detailed characterization of the ion dynamics. These experimental observations are supported by numerical simulations that provide more insight into the kinetic-scale physics associated with blob formation as well as the role of the ambient plasma density.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Arbitrary-velocity laser pulses in plasma waveguides
Space-time structured laser pulses feature an intensity peak that can travel at an arbitrary velocity while maintaining a near-constant profile. These pulses can propagate in uniform media, where their frequencies are correlated with continuous transverse wave vectors, or in structured media, such as a waveguide, where their frequencies are correlated with discrete mode numbers. Here, we demonstrate the formation and propagation of arbitrary-velocity laser pulses in a plasma waveguide where the intensity can be orders of magnitude higher than in a solid-state waveguide. The flexibility to control the velocity of the peak intensity in a plasma waveguide enables new configurations for plasma-based sources of radiation and energetic particles, including THz generation, laser wakefield acceleration, and direct laser acceleration.
42 ENGINEERING↗
A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes
We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.
97 MATHEMATICS AND COMPUTING↗
Collisionless relativistic magnetic reconnection driven by electron vortices in laser-plasma interaction
Magnetic reconnection (MR) is a fundamental process in space and laboratory plasmas. The appearance of high power lasers opens a new way to investigate MR under the relativistic condition. In this paper, relativistic collisionless MR driven by two ultra-intense lasers and a pair of asymmetric targets is studied numerically via the kinetic simulations. The static magnetic fields produced by the electron vortex structures with opposite magnetic polarities approach each other driven by the magnetic pressure and the density gradient. The antiparallel magnetic fields annihilate accompanied with the topological variation and the corresponding magnetic field energy is being dissipated to the kinetic energy of the nonthermal charged particles. Besides the outflows along the current sheet, a fast particle bunch is accelerated perpendicularly contributed by the displacement current.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Mechanisms of energetic-particle transport in magnetically confined plasmas
Super-thermal ions and electrons occur in both space and fusion plasmas. Because these energetic particles (EP) have large velocities, EP orbits necessarily deviate substantially from magnetic surfaces. Orbits are described by conserved constants of motion that define topological boundaries for different orbit types. Electric and magnetic field perturbations produced by instabilities can disrupt particle orbits, causing the constants of motion to change. The statistics of the “kicks” associated with these perturbations determines the resulting cross field transport. A unifying theme of this tutorial is the importance of the perturbation’s phase at the particle’s position Θ = k • r – ωt, where k and ω are the wavevector and frequency of the perturbation, r is the EP position, and t is the time. A distinction is made between field perturbations that resonate with an aspect of the orbital motion and those that do not. Resonance occurs when the wave phase returns to its initial value in an integer multiple of an orbital period. Convective transport occurs when resonant particles experience an unvarying wave phase. Alternatively, multiple wave-particle resonances usually decorrelate the phase, resulting in diffusive transport. Large orbits increase the number of important resonances and can cause chaotic orbits even for relatively small amplitude waves. In contrast, in the case of non-resonant perturbations, orbital phase averaging reduces transport. Here, large field perturbations introduce additional effects, including nonlinear resonances at fractional values of the orbital motion. In summary, large orbits are a blessing and a curse: For non-resonant modes, orbit-averaging reduces transport but, for resonant transport, large orbits facilitate jumps across topological boundaries and enhance the number of important resonances.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Amplification of turbulence through multiple planar shocks
We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.
Astrophysical fluid dynamics↗
Latent space mapping: Revolutionizing predictive models for divertor plasma detachment control
The inherent complexity of boundary plasma, characterized by multi-scale and multi-physics challenges, has historically restricted high-fidelity simulations to scientific research due to their intensive computational demands. Consequently, routine applications such as discharge control and scenario development have relied on faster but less accurate empirical methods. This work introduces DivControlNN, a novel machine-learning-based surrogate model designed to address these limitations by enabling quasi-real-time predictions (i.e., ~ 0.2 ms) of boundary and divertor plasma behavior. Trained on over 70,000 2D UEDGE simulations from KSTAR tokamak equilibria, DivControlNN employs latent space mapping to efficiently represent complex divertor plasma states, achieving a computational speed-up of over 10 8 compared to traditional simulations while maintaining a relative error below 20% for key plasma property predictions. During the 2024 KSTAR experimental campaign, a prototype detachment control system powered by DivControlNN successfully demonstrated detachment control on its first attempt, even for a new tungsten divertor configuration and without any fine-tuning. These results highlight the transformative potential of DivControlNN in overcoming diagnostic challenges in future fusion reactors by providing fast, robust, and reliable predictions for advanced integrated control systems.
Artificial neural networks↗
Parameter space for magnetization effects in high-energy-density plasmas
Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Velocity-space distribution function of fast ions in a sawtoothing plasma
This study explores the influence of sawtooth oscillations on the velocity space distribution of fast ions in tokamak plasma discharges. The relevant Fokker–Planck equation for fast ions is solved analytically. Two distinct effects arising from the temperature drop associated with a sawtooth crash and their impact on the distribution function of fast ions are considered. The first effect involves the modulation of the fusion alpha particle source on the timescale of the sawtooth period, linked to the drop in fusion yield resulting from the sawtooth temperature relaxations. The second effect is tied to the increase of the slowing-down time during the sawtooth ramp, causing particles born later in the sawtooth cycle to experience reduced slowing down compared to those born right after the crash, creating an accumulation-like mechanism at higher energies. In regimes where the sawtooth period is shorter than the fast ion slowing-down time, the combined influence of these effects gives rise to fast ion distribution functions that transiently exhibit positive slopes in velocity space.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗