Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary 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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Polarization-induced 2D hole gases in pseudomorphic undoped GaN/AlN heterostructures on single-crystal AlN substrates

A high-conductivity two-dimensional (2D) hole gas is the enabler of wide-bandgap p-channel transistors. Compared to commonly used AlN template substrates with high dislocation densities, the recently available single-crystal AlN substrates are promising to boost the speed and power handling capability of p-channel transistors based on GaN/AlN 2D hole gases (2DHGs) thanks to the much lower dislocation densities and the absence of thermal boundary resistance. Using plasma-assisted molecular beam epitaxy, we report the observation of polarization-induced high-density 2DHGs in undoped pseudomorphic GaN/AlN heterostructures on the single-crystal AlN substrates with high structural quality and atomic steps on the surface. The high-density 2DHG persists down to cryogenic temperatures with a record high mobility exceeding 280 cm 2 /V s and a density of 2.2 × 10 13 /cm 2 at 10 K. These results shed light on aspects of improving 2D hole mobilities and indicate significant potential of GaN/AlN 2DHG grown on bulk AlN substrates for future high performance wide-bandgap p-channel transistors.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sheath collapse at critical shallow angle due to kinetic effects

The Debye sheath is known to vanish completely in magnetised plasmas for a sufficiently small electron gyroradius and small angle between the magnetic field and the wall. This angle depends on the current onto the wall. When the Debye sheath vanishes, there is still a potential drop between the wall and the plasma across the magnetic presheath. The magnetic field angle corresponding to the predicted sheath collapse is shown to be much smaller than previous estimates, scaling with the electron-ion mass ratio and not with the square root of the mass ratio. This is shown to be a consequence of the kinetic electron and finite ion orbit width effects, which are not captured by fluid models. The wall potential with respect to the bulk plasma at which the Debye sheath vanishes is calculated. Above this wall potential, it is possible that the Debye sheath will invert.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Characterizing impurity sourcing and transport in the high temperature boundary of DIII-D wide pedestal QH-mode plasmas

Wide Pedestal QH-mode (WPQH) plasmas in the DIII-D tokamak show a sheath limited SOL, where electron density and temperature remain nearly constant along field lines from the midplane to the divertor. Consequently, parallel gradients are weak. The first Langmuir probe measurements in this regime point to a high sheath temperature with target electron temperature T e up to ~150 eV, leading to high carbon self-sputtering. Midplane carbon densities from SOLPS-ITER modeling of these plasmas in a double-null configuration fall far below experimental measurements unless full drifts are activated. In the drift-dependent SOLPS-ITER modeling, anomalous poloidally uniform, radially varying transport is adjusted to match measured radial electron temperature and density profiles in the pedestal and SOL region. The resulting carbon density (C6+) matches measured carbon densities just inside the separatrix. The C2+ density near the outer strike point is also consistent with spectroscopic imaging. With the ion B×∇B drift towards the X-point, carbon in the lower divertors is redistributed from the private flux region to the high field side (HFS) and pushed upstream in the SOL by poloidal E×B drifts. The B × ∇B drift dominates the radial flow of carbon as it moves upstream. Simulations with reversed toroidal field (ion B × ∇B drift away from the X-point) show a radically different behavior, where carbon accumulates on the low field side in lower divertors, and is pushed towards the HFS in the upper divertors, indicating a strong effect of particle drifts on impurity distribution. In double-null configurations which are usually used in WPQH plasmas, these drift effects tend to counterbalance, so that the carbon density and Z eff are reduced at the outer mid-plane only modestly, by around 7% in the modeling. Further predictive modeling indicates that the carbon density can be significantly reduced by an order of magnitude by implementing a single-null shape with reversed toroidal field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Validation of an integrated modeling framework for investigating 3D plasma responses in tokamak plasmas

As contemporary experimental tokamaks are pushed toward reactor-relevant operation, they provide essential testbeds for demonstrating and exploring ELM control strategies for deployment in future fusion pilot plants. Accurate predictions of full plasma responses are essential to guide and optimize these demonstrations. This paper introduces and validates an integrated modeling framework over a historical range of DIII-D operational space. The integrated modeling framework uses only scalar plasma parameters and optional reference boundary to self-consistently and flexibly scan through tokamak operational space and estimate the corresponding ELM-suppression relevant plasma response over 3D coil phase space. The framework generates tightly converged equilibria that satisfy a target set of plasma parameters (I p , β N , l i ), with kinetic profiles constrained by an EPED(NN)-computed pedestal and empirical core model. The plasma response of these modeled equilibria is calculated with GPEC. Validated against 55 distinct DIII-D equilibria, the framework consistently reproduces experimentally constrained equilibrium pressure, q, and other representative profiles using the time-varying information from only evolving scalar plasma parameters; plasma response validation is performed for the n=3 perturbation. In particular, the plasma response of a DIII-D discharge scanning q 95 in search of ELM suppression windows is quantitatively reproduced. By accurately reproducing experimental equilibria and plasma responses across wide parameter variation, the framework supports potential for synthetic parameter scans in key operational, stability, and plasma response dimensions to investigate RMP ELM-suppression experiments and inform predictive RMP scenario optimization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Collisionless cooling of perpendicular electron temperature in the thermal quench of a magnetized plasma

Thermal quench of a nearly collisionless plasma against an isolated cooling boundary or region is an undesirable off-normal event in magnetic fusion experiments, but an ubiquitous process of cosmological importance in astrophysical plasmas. Parallel transport theory of ambipolar-constrained tail electron loss is known to predict rapid cooling of the parallel electron temperature $T_{e\Vert}$ although $T_{e\Vert}$ is difficult to diagnose in actual experiments. Instead direct experimental measurements can readily track the perpendicular electron temperature $T_{e\bot}$ via electron cyclotron emission. The physics underlying the observed fast drop in $T_{e\bot}$ requires a resolution. Here two collisionless mechanisms, dilutional cooling by infalling cold electrons and wave-particle interaction by two families of whistler instabilities, are shown to enable fast $T_{e\bot}$ cooling that closely tracks the mostly collisionless crash of $T_{e\Vert}$. These findings motivate both experimental validation and reexamination of a broad class of plasma cooling problems in laboratory, space, and astrophysical settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Topanga: A kinetic ion plasma code for large-scale ionospheric simulations on magnetohydrodynamic timescales

Topanga is a kinetic ion code developed for simulating large-scale plasma phenomena in the Earth's ionosphere on magnetohydrodynamic timescales. It is a domain-decomposed parallel code that runs on high-performance computing platforms. Features of Topanga include spherical geometry for simplified boundary conditions and computational efficiency; a hybrid plasma model with inertia-less fluid electrons, kinetic ions, and an electric field specified via an Ohm's law; a Maxwell-FDTD (finite difference time domain) plasma model which retains the displacement current in Maxwell's equations and models electron currents in the ionosphere with a tensor conductivity; sponge-layer boundary conditions for absorption of electromagnetic and plasma waves incident on the domain boundaries; and a novel mixed-implicit algorithm for evolving the EM fields inside the Maxwell-FDTD region that is stable over many orders of magnitude in the electron–ion collision frequency. We verify the numerical methods used in Topanga on a pair of test problems. The first test involves modeling a three-dimensional collisionless shock using the hybrid set of equations. The second test involves modeling a spherical TEM mode in vacuum using the Maxwell-FDTD set of equations. Finally, we demonstrate how using the combined set of hybrid and Maxwell-FDTD equations to model the Starfish Prime high-altitude nuclear test recovers a “missing” EM signal on the ground that is not present when using only the hybrid set of equations. The magnitude of this signal in the simulation containing the Maxwell-FDTD region agrees well with the E3a portion of the magnetohydrodynamic electromagnetic pulse from Starfish Prime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anticathode effect on electron kinetics in electron beam generated E × B plasma

Electron beam (e-beam) generated plasmas with applied crossed electric and magnetic (E x B) fields are promising for low damage processing of materials with applications to microelectronics and quantum information systems. In cylindrical e-beam E x B plasmas, radial confinement of electrons and ions is achieved by an axial magnetic field and radial electric field, respectively. To control the axial confinement of electrons, such e-beam generated plasma sources may incorporate a conducting boundary known as an anticathode, which is placed on the axially opposite side of the plasma from the cathode. In this work, it is shown that varying the anticathode voltage bias can control the degree to which the anticathode collects or repels incident electrons, allowing control of warm electron (electron energies in 10–30 eV range) and beam electron population confinement. It is suggested that the effect of the anticathode bias on the formation of these distinct electron populations is also associated with the transition between weak turbulence and strong Langmuir turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Liquid lithium divertor analysis using coupled plasma material interaction model

A liquid lithium divertor can improve performance of future fusion devices by creating efficient power exhaust and improving the energy confinement via pumping of the hydrogen isotopes. In addition, significantly higher heat fluxes can be handled if controlled vapor shielding is used to redistribute the divertor heat flux over a wider area. Design and optimization of such a system calls for an analysis model which includes a strong two-way coupling between the plasma and divertor material. The incoming plasma heat and particle flux will affect the divertor surface temperature, which is a defining factor of the lithium evaporative and sputtered flux going into the plasma. Results of the coupled model based on the plasma edge code SOLPS-ITER and the computational fluid dynamics (CFD) code ANSYS-CFX will be presented for different configurations. An analytical slab flow model is used as a heat transfer boundary condition for SOLPS, defining particle flux from the wall via calculation of the surface temperature. At the final step, results of the SOLPS analysis are verified using a 3D CFD magnetohydrodynamics (MHD) analysis which uses heat and particle flux from SOLPS as a boundary condition. In addition to plasma heat flux, both analytical and CFD temperature models include several plasma material interaction effects, such as lithium evaporation, condensation and sputtering based on deuterium target flux. New adatom sputtering model based on the available experimental data is presented. Analytical model is expanded to include free surface axisymmetric configurations. Results of parametric studies of the divertor configurations with different lithium inlet temperature and velocity will be presented leading to the optimal design resulting in the lowest possible lithium contamination in the core.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Free-boundary MRxMHD equilibrium calculations using the stepped-pressure equilibrium code

Here the stepped-pressure equilibrium code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 , 112 502) is extended to enable free-boundary multi-region relaxed magnetohydrodynamic (MRxMHD) equilibrium calculations. The vacuum field surrounding the plasma inside an arbitrary 'computational boundary', , is computed, and the virtual-casing principle is used iteratively to compute the normal field on so that the equilibrium is consistent with an externally produced magnetic field. Recent modifications to SPEC are described, such as the use of Chebyshev polynomials to describe the radial dependence of the magnetic vector potential, and a variety of free-boundary verification calculations are presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Parametric Instability of Alfvén Waves and Wave Packets in Periodic and Open Systems

The parametric decay instability of Alfvén waves has been widely studied, but few investigations have examined wave packets of finite size and the effect of different boundary conditions on the growth rate. In this paper, we perform a linear analysis of circular and arc-polarized wave trains and wave packets in periodic and open boundary systems in a low- β plasma. We find that both types of wave are 3–5 times more stable in open boundary conditions compared to periodic. Additionally, once the wave packet width ℓ becomes smaller than the system size L , the growth rate decreases nearly with a power law γ ∝ ℓ / L . This study demonstrates that the stability of a pump wave cannot be separated from the laboratory settings, and that the growth rate of daughter waves depends on the conditions downstream and upstream of the pump wave and on the fraction of volume it fills. Our results can explain simulations and experiments of localized Alfvén waves. They also suggest that Alfvénic fluctuations in the solar wind, including sharp impulses known as switchbacks, can be more stable than traditional theory suggests depending on wind conditions.

Alfven waves↗
DOE OSTI2021

Direct evaluation of the phase diagrams of dense multicomponent plasmas by integration of the Clapeyron equations

Accurate phase diagrams of multicomponent plasmas are required for the modeling of dense stellar plasmas, such as those found in the cores of white dwarf stars and the crusts of neutron stars. Those phase diagrams have been computed using a variety of standard techniques, which suffer from physical and computational limitations. Here we present an efficient and accurate method that overcomes the drawbacks of previously used approaches. In particular, finite-size effects are avoided as each phase is calculated separately; the plasma electrons and volume changes are explicitly taken into account; and arbitrary analytic fits to simulation data as well as particle insertions are avoided. Furthermore, no simulations at “uninteresting” state conditions, i.e., away from the phase coexistence curves, are required, which improves the efficiency of the technique. The method consists of an adaptation of the so-called Gibbs-Duhem integration approach to electron-ion plasmas, where the coexistence curve is determined by direct numerical integration of its underlying Clapeyron equation. The thermodynamics properties of the coexisting phases are evaluated separately using Monte Carlo simulations in the isobaric semigrand canonical ensemble (NPT Δ μ ). We describe this Monte Carlo-based Clapeyron integration method, including its basic physical and numerical principles, our extension to electron-ion plasmas, and our numerical implementation. We illustrate its applicability and benefits with the calculation of the melting curve of dense carbon-oxygen plasmas under conditions relevant for the cores of white dwarf stars and provide analytic fits to implement this new melting curve in white dwarf models. While this work focuses on the liquid-solid phase boundary of dense two-component plasmas, a wider range of physical systems and phase boundaries are within the scope of the Clapeyron integration method, which had until now only been applied to simple model systems of neutral particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reduced-model scrape-off layer turbulence (nSOLT) simulations comparing three fueling scenarios

The 2D scrape-off-layer turbulence code (nSOLT) includes 1D Boltzmann neutral–plasma interactions, a model of divertor recycling (introduced here), and a fixed source of plasma concentrated at the core-side boundary. Three fueling methods are considered herein: (1) neutral injection in the far-SOL is accomplished by specifying the density of Franck–Condon distributed neutrals streaming in from the boundary. (2) Divertor recycling is modeled by injecting a fraction of the particle parallel flux in the scrape-off layer (SOL) back into the edge region as a source of plasma. (3) A constant source fuels the edge plasma from the core-side boundary to model pellet injection. For machine parameters (B, R m , and L // ) illustrative of the MAST-U device, and for a deuterium plasma, turbulent equilibria are obtained that share the same plasma fueling rate for each of the three fueling methods, with only one of the sources on in each case. In the presence of self-consistent turbulence, quasi-steady plasma and neutral (deuterium) profiles, fueling efficiencies, SOL transparencies, and heat flux widths are compared. Characteristics of the turbulent fluctuations, including skewness, cross-phases, and power spectra, are described. The calculated fueling efficiencies, SOL transparencies to neutral penetration, and many of the turbulent properties are remarkably similar for all three fueling methods despite significant differences in the plasma profiles. The nonlinear states of the three cases are dominated by separatrix-spanning vortex cells that control particle and heat losses into the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Drift kinetic theory of the NTM magnetic islands in a finite beta general geometry tokamak plasma

In (Imada et al 2019 Nucl. Fusion 59 046016 and references therein) a new 4D drift kinetic nonlinear theory, valid in the limit of a low beta, small inverse aspect ratio, circular cross section, toroidal geometry, to describe the plasma response to the neoclassical tearing mode (NTM) magnetic perturbation is derived. In (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) this theory is reduced in a low collisionality limit, which allows a dimensionality reduction to a 3D problem to efficiently resolve the collisional dissipation layer in the vicinity of the trapped-passing boundary. (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) adopts an improved model for the magnetic drift frequency, which reduces the threshold magnetic island half-width from $8.73\rho_{b i}$, where $\rho_{b i}$ is the trapped ion banana orbit width, to $1.46 \rho_{b i}$, making it in closer agreement with experimental observations for the large aspect ratio tokamak equilibrium. In the present paper, the theory is extended to a high beta, arbitrary tokamak geometry to capture the plasma shaping effects on the NTM threshold physics with the focus on the non-zero triangularity discharges that are known to have a strong impact on the plasma MHD stability. First, it is found that the higher triangularity plasma is more prone to NTMs which is in agreement with the 2/1 tearing mode onset relative frequency measurements in DIII-D. Second, the NTM threshold dependence on the tokamak inverse aspect ratio obtained in (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) is refined and extended to a finite aspect ratio limit. Third, the NTM threshold dependence on poloidal beta is obtained and benchmarked against the EAST threshold island width measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Relativistic Shear Boundary Layer and the Gamma-Ray Emission of GW170817

We present an analysis of the radiation characteristics of kinetic shear boundary layers created by relativistic plasma jets. Using a model of electromagnetic field data based on particle-in-cell simulations of an electron–ion plasma, we solve the motion of individual test electrons and compute their instantaneous radiated power and peak frequency. By analyzing a large number of test electrons in this manner, we find two distinct electron populations present around the shear boundary layer. The most highly radiative electrons execute looping motion due to crossed electric and magnetic fields as they are accelerated along the bulk flow of the jet, and eventually cross the shear boundary interface at steep angles. Electrons that never cross the shear boundary interface radiate much less energy as a group. Summing over all of the highly radiative electrons, we compute the distribution of the total radiated energy as a function of the angle relative to the bulk flow. Furthermore, this result has important potential implications for the observed radiation output of short gamma-ray bursts viewed at large angles from the jet axis, such as the neutron star merger event GW170817/GRB 170817A.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗
DOE OSTI2020

Dissipation of electron-beam-driven plasma wakes

Metre-scale plasma wakefield accelerators have imparted energy gain approaching 10 gigaelectronvolts to single nano-Coulomb electron bunches. To reach useful average currents, however, the enormous energy density that the driver deposits into the wake must be removed efficiently between shots. Yet mechanisms by which wakes dissipate their energy into surrounding plasma remain poorly understood. Here, we report picosecond-time-resolved, grazing-angle optical shadowgraphic measurements and large-scale particle-in-cell simulations of ion channels emerging from broken wakes that electron bunches from the SLAC linac generate in tenuous lithium plasma. Measurements show the channel boundary expands radially at 1 million metres-per-second for over a nanosecond. Simulations show that ions and electrons that the original wake propels outward, carrying 90 percent of its energy, drive this expansion by impact-ionizing surrounding neutral lithium. The results provide a basis for understanding global thermodynamics of multi-GeV plasma accelerators, which underlie their viability for applications demanding high average beam current.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗