Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Space plasma physics”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

The evolution of curvature in planar, photoionization-driven heat fronts

Photoionized plasmas are common in astrophysics and cosmology, especially in space near compact objects, and there are effects from photoionization in high-energy-density plasmas due to the large radiation fields present. Photoionized plasmas are an active area of laboratory research and there are currently experiments to study photoionization-supported heat fronts. These photoionization fronts differ from the physics of diffusive radiation waves, commonly called Marshak waves, that are also an active area of research. This work uses a geometric argument to describe the expected evolution of the photoionization front curvature, in a planar geometry. It then compares this curvature to that of a Marshak wave as a method of diagnosing a heat front experiment. It is found that while the curvature of a planar Marshak wave increases in time, it decreases for a photoionization front. A comparison of radiation energy and electron heat fluxes through the container for the heat front propagating medium demonstrates that the geometric argument for the photoionization front curvature is sufficient. Furthermore, this comparison also demonstrates that wall losses are not significant in a photoionization front because the post-front region is very optically thin. A discussion of the implication this work has on material choice in the targets for an experiment follows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Towards fast and accurate predictions of radio frequency power deposition and current profile via data-driven modelling: applications to lower hybrid current drive

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 modelling and real-time control applications. The machine learning models use a database of more than 16 000 GENRAY/CQL3D simulations for training, validation and testing. Latin hypercube sampling methods ensure that the database covers the range of nine input parameters ( $n_{e0}$ , $T_{e0}$ , $I_p$ , $B_t$ , $R_0$ , $n_{\|}$ , $Z_{{\rm eff}}$ , $V_{{\rm loop}}$ and $P_{{\rm LHCD}}$ ) with sufficient density in all regions of parameter space. The surrogate models reduce the inference time from minutes to $\sim$ ms with high accuracy across the input parameter space.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven surrogate modeling of hPIC ion energy-angle distributions for high-dimensional sensitivity analysis of plasma parameters' uncertainty

In this work, we present a data-driven strategy for effective construction of a surrogate model in high-dimensional parameter space for the ion energy-angle distribution (IEAD) output of hPIC simulations of plasma-surface interactions. The methodology is based on a bin-by-bin least-squares fitting of the IEAD in the parameter space. The fitting is performed in a transformed coordinate system to normalize the IEAD, and it employs sparse grids for sampling the parameter space to overcome sampling challenges in high dimensions. The surrogate model is significantly cheaper computationally than direct hPIC simulations yet maintains high fidelity to them, providing a fast emulator for hPIC simulations. Sensitivity analysis based on the surrogate model is utilized to characterize the dependence of the ion impact angle and energy moments on the physical parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Operational space for lower hybrid heating scenarios in the full tungsten environment of WEST

In tungsten—W—Environment in Steady-state Tokamak (WEST), the lower hybrid current drive (LHCD) system is key for achieving long pulse operation by providing most of the non-inductive plasma current, as well as a crucial source of electron heating. Therefore, determining the operational space for its application is fundamental. In the present study, the LHCD operational space is deeply analyzed for 0.5 MA pulses. This space is bounded by three limits: (i) the ratio of the LHCD power over density must be above a threshold to compensate tungsten radiation with enough core heating, (ii) the line-averaged density must be high enough to allow good coupling of the hybrid wave with the plasma, and (iii) fast electron ripple losses must be below a limit to avoid reaching a thermal threshold on plasma-facing components. If the tungsten radiation peak or burn-through phase is not safely overcome, a maximum electron temperature of 1.5 keV is obtained, confinement is degraded, and magnetohydrodynamic activity is frequently triggered, potentially causing a disruption. From experimental measurements and interpretative simulations, we highlight the main mechanisms that prevent the plasma from heating up during LHCD power ramp-up. Three parameters play a major role: plasma density, tungsten concentration and LHCD power deposition. A strategy to overcome this limitation is found: a precise density ramp-up performed simultaneously with the increase in LHCD power. Additionally, we show that boronization greatly facilitates the burn-through of tungsten by lowering its content during the heating phase. Finally, taking into account the three constraints given above, the LHCD operational space is determined at power ramp-up and during constant heating phases.

lower hybrid heating and current drive↗
DOE OSTI2024

Use of current-potential patches to obtain fundamental improvements to the coils of magnetic fusion devices

A central issue in the design of tokamaks or stellarators is the coils that produce the external magnetic fields. The freedom that remains unstudied in the design of coils is enormous. This freedom could be quickly studied computationally at low cost with high reliability. In particular, the space between toroidal field and modular coils that blocks access to the plasma chamber could be increased by a large factor. The concept of current-potential patches, which was developed in Todd Elder's thesis, provides a method for separating the study of the feasibility of coils with attractive physics properties from the engineering design of specific coils. In addition to enhanced accessibility, coils can be designed for increased plasma-coil separation, insensitivity to coil position errors, and plasma control.

Boozer, Allen H. [Columbia Univ., New York, NY (Un↗
DOE OSTI2022

First-principles study of L -shell iron and chromium opacity at stellar interior temperatures

Recently developed free-energy density functional theory (DFT)-based methodology for optical property calculations of warm dense matter has been applied for studying L-shell opacity of iron and chromium at T = 182 eV. We use Mermin–Kohn–Sham density functional theory with a ground-state and a fully temperature-dependent generalized gradient approximation exchange-correlation (XC) functionals. It is demonstrated that the role of XC at such a high-T is negligible due to the total free-energy of interacting system being dominated by the noninteracting free-energy term in agreement with estimations for the homogeneous electron gas. Furthermore, our DFT predictions are compared to the radiative emissivity and opacity of dense plasmas model, to the real-space Greens function method, and to experimental measurements. Good agreement is found between all three theoretical methods, and in the bound–continuum region for Cr when compared to the experiment, while the discrepancy between direct DFT calculations and the experiment for Fe remains essentially the same as for plasma-physics models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Spatial Intermittency of Particle Distribution in Relativistic Plasma Turbulence

Abstract Relativistic magnetically dominated turbulence is an efficient engine for particle acceleration in a collisionless plasma. Ultrarelativistic particles accelerated by interactions with turbulent fluctuations form nonthermal power-law distribution functions in the momentum (or energy) space, f ( γ ) d γ ∝ γ − α d γ , where γ is the Lorenz factor. We argue that in addition to exhibiting non-Gaussian distributions over energies, particles energized by relativistic turbulence also become highly intermittent in space. Based on particle-in-cell numerical simulations and phenomenological modeling, we propose that the bulk plasma density has lognormal statistics, while the density of the accelerated particles, n , has a power-law distribution function, P ( n ) dn ∝ n − β dn . We argue that the scaling exponents are related as β ≈ α + 1, which is broadly consistent with numerical simulations. Non-space-filling, intermittent distributions of plasma density and energy fluctuations may have implications for plasma heating and for radiation produced by relativistic turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Validation of theory-based models for the control of plasma currents in W7-X divertor plasmas

A theory-based model for the control of plasma currents for steady-state operation in W7-X is proposed and intended for model-based plasma control. The conceptual outline implies the strength of physics-based models: it offer approaches applicable to future conditions of fusion devices or next-step machines. The application at extrapolated settings is related to the validity range of the theory model. Therefore, the predictive power of theory-based control models could be larger than for data-driven approaches and limitations can be predicted from the validity range for the prediction of bootstrap currents in W7-X. The model predicts the L/R response when density or heating power is changed. The model is based on neoclassical bootstrap current calculations and validated for different discharge conditions. While the model was found to be broadly applicable for conducted electron-cyclotron-heated discharges in W7-X, limits were found for cases when the polarization of the electron cyclotron heating was changed from X2 to O2-heating. The validity assessment attempts to quantify the potential of the derived model for model-based control in the operational space (density, heating power) of W7-X.

neoclassical modelling↗
DOE OSTI2021

Filamentary surface plasma discharge flow length and time scales

Nanosecond surface dielectric barrier discharges (ns-SDBDs) are a class of plasma actuators that utilize a high-voltage pulse of nanosecond duration between two surface-mounted electrodes to create an electrical breakdown of air, along with rapid heating. These actuators usually produce multiple filaments when operated at high pulse frequencies, and the rapid heating leads to the formation of shock waves and complex flow fields. In this work we replicate a single filament of the ns-SDBDs and characterize the induced flow using velocity measurements from particle image velocimetry and density measurements from background-oriented schlieren. The discharge is produced by a high voltage electrical pulse between two copper electrodes on an acrylic base. A hot gas kernel characterizes the flow field formed close to the electrodes that expands and cools over time and a vortex ring that propagates away from the surface while entraining cold ambient fluid. The gas density deficit inside the kernel displays a power-law decay over time. Based on the observations, we develop a simplified theoretical model based on vortex-driven cooling and perform a scaling analysis to obtain the induced flow length and time scales. The results show that the cooling process's time scales correspond to a circulation-based time scale of the vortex ring, and the length scale of the kernel corresponds to the vortex ring radius. Furthermore, these findings can guide the choice of optimal filament spacing and pulse frequencies in the design, deployment, and operation of ns-SDBDs for flow control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

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↗
DOE OSTI2025

Analytical and numerical studies of dark current in radiofrequency structures for short-pulse high-gradient acceleration

High-gradient acceleration is a key research area that could enable compact linear accelerators for future colliders, light sources, and other applications. In the pursuit of high-gradient operation, rf breakdown limits the attainable accelerating gradient in normal-conducting rf structures. Recent experiments at the Argonne Wakefield Accelerator suggest a promising approach: using short rf pulses with durations of a few nanoseconds. Experimental studies show that these 𝒪⁡(1⁢ ⁢ns) rf pulses can mitigate breakdown limitations, resulting in higher gradients. For example, an electric field of nearly 400 MV/m was achieved in an 𝑋-band photoemission gun driven by 6-ns-long rf pulses, with rapid rf conditioning and low dark current observed. Despite these promising results, the short-pulse regime remains an underexplored parameter space, and rf breakdown physics under nanosecond-long pulses requires further investigation. In this paper, we present analytical and numerical simulations of dark current dynamics in accelerating cavities operating in the short-pulse regime. We study breakdown-associated processes spanning different time scales, including field emission, multipacting, and plasma formation, using simulations of the 𝑋-band photogun cavities. The results reveal the advantages of using short rf pulses to reduce dark current and mitigate rf breakdown, offering a path toward a new class of compact accelerators with enhanced performance and reduced susceptibility to breakdown.

Linear accelerators↗
DOE OSTI2021

VLBI Data Processing on Coronal Radio-sounding Experiments of Mars Express

The ESA’s Mars Express solar corona experiments were performed at two solar conjunctions in the years 2015 and 2017 by a number of radio telescopes in the European VLBI Network. This paper presents the methods to measure the frequency and phase fluctuations of the spacecraft radio signal, and the applications to study the characteristics of the plasma turbulence effects on the signal at a single station and at multiple stations via cross correlation. The power spectra of the frequency fluctuations observed between 4.9 and 76.3 R {sub s} have a power-law shape close to a Kolmogorov spectrum over the frequency interval ν {sub lo} < ν < ν {sub up}, where the nominal value of ν {sub lo} is set to 3 mHz and ν {sub up} is in the range of 0.03–0.15 Hz. The rms of the frequency fluctuations is presented as a function of the heliocentric distance. Furthermore, we analyze the variations of the electron column density fluctuations at solar offsets 4.9 R {sub s} and 9.9 R {sub s} and the cross-correlation products between the VLBI stations. The power density of the differential fluctuations between different stations decreases at ν < 0.01 Hz. Finally, the fast flow speeds of solar wind >700 km s{sup -1} are derived from the cross correlation of frequency fluctuations at ν < 0.01 Hz. The fast flow speeds of solar wind correspond to the high heliolatitude of the coronal region that the radio rays passed. The VLBI observations and analysis methods can be used to study the electron column density fluctuations and the turbulence at multiple spatial points in the inner solar wind by providing multiple lines of sight between the Earth and the spacecraft.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

The development of a Sustained High Power Density (SHPD) facility

A Sustained High Power Density (SHPD) facility is being planned within the U.S. as a non DT device which can be a bridge to a Compact Fusion Pilot Plant (CFPP) DT operated facility with the future capacity to produce a levelised cost of electricity at a rate competitive with current power generating systems. To accomplish this a review of ITER cost data and recent next step physics/engineering studies that promoted technology development, physics enhancements and innovated configuration improvements that increase operating availability have been injected into a proposed PPPL SHPD design. Much of the physics and engineering design basis for this effort centers on a series of ST studies carried out in recent years attempting to maximize the low aspect ratio (AR) ST performance within realistic engineering constraints of this compact fusion device. One issue found in the low AR design is the lack of volt-seconds generated by a small central solenoid that resides in the available space left after sizing the TF structure to support the high current density/HTS TF coil. Although not required in a larger pilot plant device, a solenoid wrapped around the plasma side of the SC TF coil was considered for the small low aspect ratio SHPD device. Adding a solenoid to the plasma side of the TF imposes a major change to the device configuration – a change which has been found to have cost and assembly issues when compared with an OH solenoid located at the machine center. Finally, the inclusion of a liquid metal first wall divertor system, the integration of local pairs of outboard DCLL blanket segments and a summary of physics performance conditions will be presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anomalous Retrograde Drifts in Obstructed Magnetron Microdischarges: a Consequence of a Field Reversal in the Anode Sheath?

The goal of this project was to carry out experiments at the Plasma Collaborative Research Facility (Princeton Plasma Physics Laboratory) aimed at confirming or refuting the hypothesis that an electric field reversal occurs in the near-anode region of a strongly obstructed planar magnetron microdischarges. Magnetron discharges are magnetized plasmas with E x B field topologies, and variants of such discharges are widely used in materials processing. The topologies are also similar to that of Hall thrusters used in space propulsion. Understanding the structure of these discharges will aid in developing reliable simulations that are intended to facilitate the further development of such technologies, as the field reversals, presumed to be due to strong axial-driven electron transport, provides further insight into transport properties. The possible existence of such a field reversal is hinted at because of the retrograde-propagating (i.e., negative E x B) spokes or gradient-driven drift instabilities seen in the high speed videos and segmented anode probing of the plasma. The presence of this field reversal should be reflected in the ion dynamics, particularly, in an ion drift towards the anode. The proposed experiments to be conducted at the PCRF are to characterize this anode-directed ion drift using state-of-the-art laser-induced fluorescence velocimetry. In particular, we seek to see if the near-anode ion velocity reverses direction when the external circuit is varied to alter the electron current demand imposed at the anode. Such an ion velocity reversal would confirm the mechanism, which is that the highly obstructed (small gapped) discharge produces strong axial gradients in plasma density, driving more electrons towards the anode then that demanded by the external circuit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Fast physics-based launcher optimization for electron cyclotron current drive

With the increased urgency to design fusion pilot plants, fast optimization of electron cyclotron current drive (ECCD) launchers is paramount. Traditionally, this is done by coarsely sampling the 4D parameter space of possible launch conditions consisting of (1) the launch location (constrained to lie along the reactor vessel), (2) the launch frequency, (3) the toroidal launch angle, and (4) the poloidal launch angle. For each initial condition, a ray-tracing simulation is performed to evaluate the ECCD efficiency. Unfortunately, this approach often requires a large number of simulations (sometimes millions in extreme cases) to build up a dataset that adequately covers the plasma volume, which must then be repeated every time the design point changes. Here we adopt a different approach. Rather than launching rays from the plasma periphery and hoping for the best, we instead directly reconstruct the optimal ray for driving current at a given flux surface using a reduced physics model coupled with a commercial ray-tracing code. Repeating this throughout the plasma volume requires only hundreds of simulations, constituting a significant speedup. The new method is validated on two separate example tokamak profiles, and is shown to reliably drive localized current at the specified flux surface with the same optimal efficiency as obtained from the traditional approach.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗