Search NASA⌕ Search

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.

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 OSTI2022

Comparison of Energetic Electron Fluxes Measured by GPS and THEMIS Spacecraft in the Inner Magnetosphere

Key elements of space weather models are energetic electron fluxes in the inner magnetosphere and the outer radiation belt. Flux depletion is driven by various loss processes: scattering into atmosphere, magnetopause shadowing. Flux enhancement is driven by various acceleration processes: local wave-particle interactions, radial transport, plasma sheet injections. Many of these processes operate on ~ hour timescales. Such mesoscale flux variations are not well traced by equatorial spacecraft with much longer orbits. Energetic electron detectors onboard the Global Positioning System (GPS) constellation provide a unique opportunity for probing such ~ hour-scale flux variations. Measurements from up to 23 identically instrumented GPS satellites cover a wide energy and L-shell range with a subhour time resolution. However, their orbits are inclined and thus all measurements at L-shell >4.3 are off-equatorial. Here, in this report, we present a comparison of equatorial THEMIS and nonequatorial GPS measurements of omnidirectional ≤600 keV electron fluxes. Such a comparison allows us to derive coefficients for using off-equatorial GPS fluxes to infer the equatorial values. These coefficients depend on particle energy and L-shell. We demonstrate a new data set derived from GPS measurements and discuss how it can be used to investigate mesoscale dynamics of energetic electron fluxes in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Novel intrinsic helical cores and MHD dynamo flux pumping evidence in DIII-D

Evidence of a MHD dynamo brought about by intrinsic helical core structures in negative triangularity (NT) and ITER baseline scenario (IBS) plasmas demonstrates paths to operating tokamaks without seeds for deleterious global MHD modes. A stationary sawtooth- and edge localized mode-free state maintains high performance for long durations in a strong NT shape in the DIII-D tokamak. Disparity between minimum safety factor values predicted by magnetic flux diffusion theory and derived from experiment suggests anomalous flux diffusion by an MHD dynamo, driven by flows associated with quasi-interchange modes. Stability calculations with the GATO code find NT and IBS experimental equilibria unstable to quasi-interchange modes, in agreement with observed saturated activity. IBS discharges with similar magnetic winding structure occupy robust helical core states following a transition from sawtooth instabilities. These states agree with prior modeling of helical core bifurcation thresholds with the VMEC equilibrium code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Gauged global strings

We investigate the string solutions and cosmological implications of the gauge U(1) Z × global U(1) PQ model. With two hierarchical symmetry-breaking scales, the model exhibits three distinct string solutions: a conventional global string, a global string with a heavy core, and a gauge string as a bound state of the two global strings. This model reveals rich phenomenological implications in cosmology. During the evolution of the universe, these three types of strings can form a Y-junction configuration. Intriguingly, when incorporating this model with the QCD axion framework, the heavy-core global strings emit more axion particles compared to conventional axion cosmic strings due to their higher tension. This radiation significantly enhances the QCD axion dark matter abundance, thereby opening up the QCD axion mass window. Consequently, axions with masses exceeding ~ 10 −5 eV have the potential to constitute the whole dark matter abundance. Furthermore, in contrast to conventional gauge strings, the gauge strings in this model exhibit a distinctive behavior by radiating axions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Soft X-ray measurements of transport and MHD activity in the core and edge NSTX plasma (Final Technical Report)

A comprehensive suite of innovative diagnostic instrumentation and techniques has been developed by the Johns Hopkins collaboration for measurements of transport and MHD activity in the core and edge of the National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory (PPPL). These include ultrafast multi-energy poloidal soft X-ray arrays, high spatial resolution multi-energy tangential soft X-ray arrays, high spectral resolution transmission grating imaging VUV spectrometers and ultrafast bolometric arrays, which measure over a broad range of spatial, spectral, and temporal scales the emission from ohmic, RF, and beam driven NSTX plasmas. The instrumentation was designed to fulfill the specific measurement needs of the Spherical Torus (ST) magnetic fusion concept and is accompanied by a comprehensive set of numerical tools for data analysis and modeling, aimed at determining as many as possible plasma parameters and at extracting as much as possible physics information from NSTX experiments. The Johns Hopkins diagnostics have been and continue to be instrumental for essentially all the NSTX experiments. In addition, our own research and experiments using these diagnostics has led over the years to fundamental conclusions for the ST fusion physics, such as the establishment of near neoclassical impurity transport in the NSTX core and edge plasma, and the critical role of the Global Alfven Eigenmode (GAE) MHD activity for the ST core electron thermal transport. The Johns Hopkins Diagnostics continue to be operational on the upgraded NSTX-U and will provide key plasma information in its upcoming high field, current and power experimental campaigns. The present report summarizes the main technical characteristics and measurement capabilities of the diagnostic instrumentation and methods developed within our research and describes a key physics result achieved during this research. Success in our research was made possible also by the excellent support and cooperation of the NSTX research team and of PPPL researchers, who are warmly acknowledged here.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Erosion and impurity transport for the edge localized mode suppression window in KSTAR

A 3D analysis of plasma wall interactions and global impurity transport for the edge localized mode suppression window in KSTAR during H-Mode discharges has shown that carbon erosion at the divertor plates is a strong function of the resonant magnetic perturbation (n = 1) coil current and relative phasing. The Generalized Perturbation Equilibrium Code was used to determine a realistic initial perturbed plasma equilibrium, and EMC3-EIRENE was used to calculate the resulting scrape-off layer plasma used in this study as a fixed background for the ERO2.0 plasma–material interaction model. The resulting transport leads to deposition of impurities along the targets positioned at the high-field side of the device. An attempt at calculating the resulting effective charge state has demonstrated a similar dependence on the perturbation coil current and has been able to determine a window for the experimentally observed values of Zeff by including contributions of all ionized carbon charge states and deuterium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Global particle buildup simulations with gas puff scan: application to WEST discharge

This paper deals with the distribution of sources, transport, and exhaust of particles in a tokamak. Knowledge and understanding of all the physical phenomena involved in the global particle buildup are necessary to study and predict density regimes and subsequently to develop optimized scenarios for tokamak operation in order to control heat and particle exhaust. Neutral particles and their interactions with plasma are central in this perspective. This paper discusses the impact of varying the intensity of particle fueling in 2D transport simulations of a WEST discharge. Simulations are performed with an updated version of SOLEDGE-HDG that allows a more realistic transport of neutrals using a self-consistent diffusive model based on charge exchange and ionization processes. New code capabilities allow the entire WEST poloidal cross section to be simulated in a realistic configuration for both geometry and the range of control parameters. A gas puff scan illustrates the main features of the sheath-limited, high-recycling, and detached regimes, such as the buildup of the temperature gradient and the pressure drop in the scrape-off layer (SOL), the target temperature falling to 1 eV, and the ionization source moving away from the targets, as well as the particle flux rollover. A crude estimate of wall erosion is also provided, showing the respective role of each plasma wall component in each of these regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Hamiltonian kinetic-Hall magnetohydrodynamics with fluid and kinetic ions in the current and pressure coupling schemes

We present two generalized hybrid kinetic-Hall magnetohydrodynamics (MHD) models describing the interaction of a two-fluid bulk plasma, which consists of thermal ions and electrons, with energetic, suprathermal ion populations described by Vlasov dynamics. The dynamics of the thermal components are governed by standard fluid equations in the Hall MHD limit with the electron momentum equation providing an Ohm's law with Hall and electron pressure terms involving a gyrotropic electron pressure tensor. The coupling of the bulk, low-energy plasma with the energetic particle dynamics is accomplished through the current density (current coupling scheme; CCS) and the ion pressure tensor appearing in the momentum equation (pressure coupling scheme; PCS) in the first and the second model, respectively. The CCS is a generalization of two well-known models, because in the limit of vanishing energetic and thermal ion densities, we recover the standard Hall MHD and the hybrid kinetic-ions/fluid-electron model, respectively. This provides us with the capability to study in a continuous manner, the global impact of the energetic particles in a regime extending from vanishing to dominant energetic particle densities. The noncanonical Hamiltonian structures of the CCS and PCS, which can be exploited to study equilibrium and stability properties through the energy-Casimir variational principle, are identified. As a first application here, we derive a generalized Hall MHD Grad–Shafranov–Bernoulli system for translationally symmetric equilibria with anisotropic electron pressure and kinetic effects owing to the presence of energetic particles using the PCS.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Three-dimensional Dynamics of Strongly Twisted Magnetar Magnetospheres: Kinking Flux Tubes and Global Eruptions

Abstract The origins of the various outbursts of hard X-rays from magnetars (highly magnetized neutron stars) are still unknown. We identify instabilities in relativistic magnetospheres that can explain a range of X-ray flare luminosities. Crustal surface motions can twist the magnetar magnetosphere by shifting the frozen-in footpoints of magnetic field lines in current-carrying flux bundles. Axisymmetric (2D) magnetospheres exhibit strong eruptive dynamics, i.e., catastrophic lateral instabilities triggered by a critical footpoint displacement of ψ crit ≳ π . In contrast, our new three-dimensional (3D) twist models with finite surface extension capture important non-axisymmetric dynamics of twisted force-free flux bundles in dipolar magnetospheres. Besides the well-established global eruption resulting (as in 2D) from lateral instabilities, such 3D structures can develop helical, kink-like dynamics, and dissipate energy locally (confined eruptions). Up to 25% of the induced twist energy is dissipated and available to power X-ray flares in powerful global eruptions, with most of our models showing an energy release in the range of the most common X-ray outbursts, ≲10 43 erg. Such events occur when significant energy builds up while deeply buried in the dipole magnetosphere. Less energetic outbursts likely precede powerful flares, due to intermittent instabilities and confined eruptions of a continuously twisting flux tube. Upon reaching a critical state, global eruptions produce the necessary Poynting-flux-dominated outflows required by models prescribing the fast radio burst production in the magnetar wind—for example, via relativistic magnetic reconnection or shocks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

An explicit scheme to enforce charge conservation in transient Particle-in-Cell simulations with Maxwell-Boltzmann electrons

The Maxwell-Boltzmann electron model is a very popular approximation of the electronic behavior in electrostatic plasmas, adopted in all cases where the electron dynamics is governed only by a balance between electric and pressure forces. In such a model, the electron dynamics is reduced to solving for a nonlinear Poisson problem together with an additional expression enforcing charge conservation. In this work we derive an expression of charge conservation which can be conveniently applied to explicit schemes to update the reference Boltzmann density, particularly useful when simulating plasma sheaths. The scheme can be equally applied to both steady-state and transient problems. Two examples are shown, a steady-state plasma sheath, and a radio-frequency magnetic presheath. Our proposed scheme allows to enforce global charge conservation locally in time, and can thus be applied to the simulation of transient phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Verification of a fully kinetic ion model for electromagnetic simulations of high-frequency waves in toroidal geometry

For the research of high-frequency electromagnetic waves in tokamaks, an electromagnetic simulation model, in which the ion dynamics is described by a six-dimensional Vlasov equation and the electron dynamics is described by a drift kinetic equation, is formulated and implemented in the global gyrokinetic toroidal code (GTC). Analytic dispersion relations are derived in reduced systems and compared with various theories to verify the model. Linear simulations of a generalized ion Bernstein wave and ion cyclotron emission are verified by comparing the GTC simulation results with analytic dispersion relation theory and magnetoacoustic cyclotron instability theory, respectively, in cylindrical geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Computational investigation of incident ion angles and material erosion at rough graphite and silicon carbide divertor surfaces

We present a computational investigation of the dependence of material erosion on the incident ion angle at rough graphite and silicon carbide divertor surfaces. Ion angle distributions (IADs) for D plasmas at NSTX-U and DIII-D divertors were calculated by an equation-of-motion model that traces the ion trajectories in the sheath. Then, the effective sputtering yields and ion shadowed area fractions were calculated by a Monte Carlo micro-patterning and roughness code that applied the calculated IADs to surface topographic data that were obtained from optical confocal microscopy of rough graphite and SiC divertor surfaces from NSTX-U and DIII-D experiments. The calculations found that the effective sputtering yields, the sputtering pattern, and the shadowed area are determined by the detailed surface topology rather than the root mean square roughness R RMS , which represents deviations from a flat surface. The suppression of the effective sputtering yields for rough surfaces compared to the yield for a smooth surface was accounted for by the change of the mean local incident ion angle (LIIA) (θ'). The mean surface inclination angle distribution (SIAD) (δ) was found to be a useful parameter to estimate the LIIA from the calculated IADs. We report global empirical formulas for the mean LIIA and fraction of the area shadowed from the main ions for D plasmas for rough surfaces with B-field incident angles α = 85°–89° as a function of the mean SIAD (δ). We propose the use of the mean LIIA (θ') to estimate the sputtering yield for rough surfaces from the angular dependence of the sputtering yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Bayesian calibration of a hybrid nuclear collision model using $\textit{p}$₋Pb and Pb-Pb data at energies available at the CERN Large Hadron Collider

Here, we posit a unified hydrodynamic and microscopic description of the quark-gluon plasma (QGP) produced in ultrarelativistic p₋Pb and Pb-Pb collisions at $\sqrt{s}_\text{NN}$ = 5.02 TeV and evaluate our assertion using Bayesian inference. Specifically, we model the dynamics of both collision systems using initial conditions with parametric nucleon substructure, a preequilibrium free streaming stage, event-by-event viscous hydrodynamics, and a microscopic hadronic afterburner. Free parameters of the model, which describe the initial state and QGP medium are then simultaneously calibrated to fit charged-particle yields, mean $p_T$, and flow cumulants. We argue that the global agreement of the calibrated model with the experimental data strongly supports the existence of hydrodynamic flow in small collision systems at ultrarelativistic energies, and that the flow produced develops at length scales smaller than a single proton. Posterior estimates for the model's input parameters are obtained, and new insights into the temperature dependence of the QGP transport coefficients and event-by-event structure of the proton are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

The aim of the study is to analyze the stability of the energetic particle modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. Additionally, the analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n = 1 EPM and n = 2 global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis ($\rlap{-} \iota_{0}$). In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given $\rlap{-} \iota_{0}$, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n = 1 EPM, n = 2 toroidal AE (TAE) and n = 3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP β threshold with respect to the co- and no-ECCD simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

On the Effect of Beating during Nonlinear Frequency Chirping

Spectroscopic analyses of energetic particle (EP) driven bursts of MHD fluctuations in magnetically confined plasmas often exhibit chirps that occur simultaneously in groups of two or more. While the superposition of oscillations at multiple frequencies necessarily causes beating in the signal acquired by a localized external probe, self-consistent hybrid simulations of chirping EP modes in a JT-60U tokamak plasma have demonstrated the possibility of global beating, where the mode's electromagnetic field vanishes globally between beats and reappears with opposite phase [Bierwage et al., Nucl. Fusion 57, 016036 (2017)]. This implies that there can be a single coherent field mode that oscillates at multiple frequencies simultaneously when it is resonantly driven by multiple density waves in EP phase space. Conversely, this means that the EP density waves are mutually coupled and interfere with each other via the jointly driven field, a mechanism ignored in some theories of chirping. In this thesis-style treatise, we study the role of field pulsations in general and beating in particular using the Hamiltonian guiding center orbit-following code ORBIT with a reduced wave-particle interaction model in realistic geometry. Beating is found to drive the evolution of EP phase space structures. Additionally, a key mechanism is the pulsation of effective phase space islands combined with the alternation of their effective O- and X-points due to phase jumps between each beat. Observations: (1) Beating causes density wave fronts to advance radially in a pulsed manner and the resulting chirps become staircase-like. (2) The pulsations facilitate convective transfer of material between neighboring layers of phase space density waves. On the one hand, this may inhibit the early detachment of solitary phase space vortices. On the other hand, it facilitates the accumulation of hole and clump fragments into larger structures. (3) Long-range chirping is observed when massive holes or clumps detach and drift away from the turbulent belt around the seed resonance. It is remarkable that the detached vortices remain robust and,on average, maintain their concentric nested layers while being visibly perturbed by the field's continued beating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Meso-Scale Electrodynamic Coupling of the Earth Magnetosphere-Ionosphere System

Within the fully integrated magnetosphere-ionosphere system, many electrodynamic processes interact with each other. We review recent advances in understanding three major meso-scale coupling processes within the system: the transient field-aligned currents (FACs), mid-latitude plasma convection, and auroral particle precipitation. (1) Transient FACs arise due to disturbances from either dayside or nightside magnetosphere. As the interplanetary shocks suddenly compress the dayside magnetosphere, short-lived FACs are induced at high latitudes with their polarity successively changing. Magnetotail dynamics, such as substorm injections, can also disturb the current structures, leading to the formation of substorm current wedges and ring current disruption. (2) The mid-latitude plasma convection is closely associated with electric fields in the system. Recent studies have unraveled some important features and mechanisms of subauroral fast flows. (3) Charged particles, while drifting around the Earth, often experience precipitating loss down to the upper atmosphere, enhancing the auroral conductivity. Recent studies have been devoted to developing more self-consistent geospace circulation models by including a better representation of the auroral conductance. It is expected that including these new advances in geospace circulation models could promisingly strengthen their forecasting capability in space weather applications. The remaining challenges especially in the global modeling of the circulation system are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Impact of various DIII-D diagnostics on the accuracy of neural network surrogates for kinetic EFIT reconstructions

Abstract Kinetic equilibrium reconstructions make use of profile information such as particle density and temperature measurements in addition to magnetics data to compute a self-consistent equilibrium. They are used in a multitude of physics-based modeling. This work develops a multi-layer perceptron (MLP) neural network (NN) model as a surrogate for kinetic Equilibrium Fitting (EFITs) and trains on the 2019 DIII-D discharge campaign database of kinetic equilibrium reconstructions. We investigate the impact of including various diagnostic data and machine actuator controls as input into the NN. When giving various categories of data as input into NN models that have been trained using those same categories of data, the predictions on multiple equilibrium reconstruction solutions (poloidal magnetic flux, global scalars, pressure profile, current profile) are highly accurate. When comparing different models with different diagnostics as input, the magnetics-only model outputs accurate kinetic profiles and the inclusion of additional data does not significantly impact the accuracy. When the NN is tasked with inferring only a single target such as the EFIT pressure profile or EFIT current profile, we see a large increase in the accuracy of the prediction of the kinetic profiles as more data is included. These results indicate that certain MLP NN configurations can be reasonably robust to different burning-plasma-relevant diagnostics depending on the accuracy requirements for equilibrium reconstruction tasks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Dijet Acoplanarity in CUJET3 as a Probe of the Nonperturbative Color Structure of QCD Perfect Fluids

Using the CUJET3=DGLV+VISHNU jet-medium interaction framework, we show that dijet azimuthal acoplanarity in high energy A + A collisions is sensitive to possible non-perturbative enhancement of the jet transport coefficient, $\hat{q}$(T,E), in the QCD crossover temperature T ~ 150–300 MeV range. With jet-medium couplings constrained by global RHIC& LHC χ 2 fits to nuclear modification data on R AA (p T > 20) GeV, we compare predictions of the medium induced dijet transverse momentum squared, $Q_s^2~\langle \hat{q}L\rangle~ΔΦ^2E^2$, in two models of the temperature, T, and jet energy E dependence of the jet medium transport coefficient, $\hat{q}$(T,E)). In one model, wQGP, the chromo degrees of freedom (dof) are approximated by a perturbative dielectric gas of quark and gluons dof. In the second model, sQGMP, we consider a nonperturbative partially confined semi-Quark-Gluon-Monopole-Plasma with emergent color magnetic dof constrained by lattice QCD data. Unlike the slow variation of the scaled jet transport coefficient, $\hat{q}_{wQGP}/T^3$, the sQGMP model $\hat{q}_{sQGMP}/T^3$ features a sharp maximum in the QCD confinement crossover T range. We show that the dijet path averaged medium induced azimuthal acoplanarity, Δφ 2 , in sQGMP is robustly ~ 2 times larger than in perturbative wQGP. even though the radiative energy loss in both models is very similar as needed to fit the same R AA data. Future A+A dijet acoplanarity measurements constrained together with single jet R AA and ν n measurements therefore appears to be a promising strategy to search for possible signatures of critical opalescence like phenomena in the QCD confinement temperature range.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

MHD modeling of shattered pellet injection in JET

Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗