Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Magnetized 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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Observation and analysis of abrupt changes in the interplanetary plasma velocity and magnetic field.

This paper presents a limited study of the physical nature of abrupt changes in the interplanetary plasma velocity and magnetic field based on 19 day's data from the Pioneer 6 spacecraft. The period was chosen to include a high-velocity solar wind stream and low-velocity wind. Abrupt events were accepted for study if the sum of the energy density in the magnetic field and velocity changes was above a specified minimum. A statistical analysis of the events in the high-velocity solar wind stream shows that Alfvenic changes predominate. This conclusion is independent of whether steady state requirements are imposed on conditions before and after the event. Alfvenic changes do not dominate in the lower-speed wind. This study extends the plasma field evidence for outwardly propagating Alfvenic changes to time scales as small as 1 min (scale lengths on the order of 20,000 km).

Martin, R. N.↗
DOE OSTI2020

Adaptive Hybridizable Discontinuous Galerkin discretization of the Grad–Shafranov equation by extension from polygonal subdomains

In this work, we propose a high-order adaptive numerical solver for the semilinear elliptic boundary value problem modeling magnetic plasma equilibrium in axisymmetric confinement devices. In the fixed boundary case, the equation is posed on curved domains with piecewise smooth curved boundaries that may present corners. The solution method we present is based on the hybridizable discontinuous Galerkin method and sidesteps the need for geometry-conforming triangulations thanks to a transfer technique that allows to approximate the solution using only a polygonal subset as computational domain. Moreover, the solver features automatic mesh refinement driven by a residual-based a posteriori error estimator. As the mesh is locally refined, the computational domain is automatically updated in order to always maintain the distance between the actual boundary and the computational boundary of the order of the local mesh diameter. Numerical evidence is presented of the suitability of the estimator as an approximate error measure for physically relevant equilibria with pressure pedestals, internal transport barriers, and current holes on realistic geometries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

An efficient high-order numerical solver for diffusion equations with strong anisotropy

In this paper, we present an interior penalty discontinuous Galerkin finite element scheme for solving diffusion problems with strong anisotropy arising in magnetized plasmas for fusion applications. Additionally, we demonstrate the accuracy produced by the high-order scheme and develop an efficient preconditioning technique to solve the corresponding linear system, which is robust to the mesh size and anisotropy of the problem. Several numerical tests are provided to validate the accuracy and efficiency of the proposed algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Conservative discontinuous Galerkin interpolation: Sheared boundary conditions

Local studies of accretion disks and laboratory magnetized plasmas employ analytical coordinate mappings that introduce sheared boundary conditions (BCs). Here, we present a discontinuous Galerkin (DG) algorithm to apply such BCs based on projections and quadrature-free integration. The procedure is high-order accurate, preserves moments exactly and works in multiple dimensions. Tests of increasing complexity are provided, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) BCs. Results show that the algorithm is (p + 1)-order accurate in the DG representation and (p + 2) -order accurate in the cell averages, with p being the order of the polynomial basis. Quantification of the algorithm's hyperdiffusion and discussion of aliasing errors are given. This technique enables conservative local simulations of plasma turbulence with DG, not possible until now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

TokaMaker: An open-source time-dependent Grad-Shafranov tool for the design and modeling of axisymmetric fusion devices

Here, in this paper, we present a new static and time-dependent MagnetoHydroDynamic (MHD) equilibrium code, TokaMaker, for axisymmetric configurations of magnetized plasmas, based on the well-known Grad-Shafranov equation. This code utilizes finite element methods on an unstructured triangular grid to enable capturing accurate machine geometry and simple mesh generation from engineering-like descriptions of present and future devices. The new code is designed for ease of use without sacrificing capability and speed through a combination of Python, Fortran, and C/C++ components. A detailed description of the numerical methods of the code, including a novel formulation of the boundary conditions for free-boundary equilibria, and validation of the implementation of those methods using both analytic test cases and cross-code validation is shown. Results show expected convergence across tested polynomial degree for analytic and cross-code test cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Collisional effects on resonant particles in quasilinear theory

A careful examination of the effects of collisions on resonant wave–particle interactions leads to an alternate interpretation and deeper understanding of the quasilinear operator originally formulated by Kennel & Engelmann ( Phys. Fluids , vol. 9, 1966, pp. 2377–2388) for collisionless, magnetized plasmas, and widely used to model radio frequency heating and current drive. The resonant and nearly resonant particles are particularly sensitive to collisions that scatter them out of and into resonance, as for Landau damping as shown by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). As a result, the resonant particle–wave interactions occur in the centre of a narrow collisional boundary when the collision frequency $\unicode[STIX]{x1D708}$ is very small compared to the wave frequency $\unicode[STIX]{x1D714}$ . The diffusive nature of the pitch angle scattering combined with the wave–particle resonance condition enhances the collision frequency by $(\unicode[STIX]{x1D714}/\unicode[STIX]{x1D708})^{2/3}\gg 1$ , resulting in an effective resonant particle collisional interaction time of $\unicode[STIX]{x1D70F}_{\text{int}}\sim (\unicode[STIX]{x1D708}/\unicode[STIX]{x1D714})^{2/3}/\unicode[STIX]{x1D708}\ll 1/\unicode[STIX]{x1D708}$ . A collisional boundary layer analysis generalizes the standard quasilinear operator to a form that is fully consistent with Kennel–Englemann, but allows replacing the delta function appearing in the diffusivity with a simple integral (having the appropriate delta function limit) retaining the new physics associated with the narrow boundary layer, while preserving the entropy production principle. The limitations of the collisional boundary layer treatment are also estimated, and indicate that substantial departures from Maxwellian are not permitted.

Physics↗
DOE OSTI2020

Theory of the tertiary instability and the Dimits shift within a scalar model

The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in a magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry–Horton system has been proposed by St-Onge ( J. Plasma Phys. , vol. 83, 2017, 905830504) as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$ , where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U'' \doteq {\rm d}^2 U/{\rm d} x^2 $ . Then, depending on $U''$ , the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $| {\rm d} U / {\rm d} x |$ . Our analytic predictions are in agreement with direct numerical simulations of the modified Terry–Horton system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Electromagnetic pinned solitons for space debris detection

Electromagnetic “pinned” solitons in the form of stationary nonlinear waves are studied within the framework of an inertial magneto-hydrodynamic model. These structures, that can arise when a charged source moves in a magnetized plasma, have a velocity that is equal to the source velocity and, hence, appear as “pinned” structures that envelope the source. We investigate the excitation of such solitons in the Low Earth Orbit region due to the passage of charged orbital debris objects. The spatial size of these electromagnetic solitons, typically of the order of a few ion skin depths, can be very large in this region. Such solitons can be detected using a variety of ground- or orbit-based radio sounding techniques and may provide a convenient additional means of tracking small sized orbital debris objects that are difficult to spot optically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Neutrino and pair creation in reconnection-powered coronae of accreting black holes

A ubiquitous feature of accreting black hole systems is their hard X-ray emission which is thought to be produced through Comptonization of soft photons by electrons and positrons in the vicinity of the black hole, in a region with optical depth of order unity. The origin and composition of this Comptonizing region, known as the corona, is a matter open for debate. In this paper we investigate the role of relativistic protons accelerated in black-hole magnetospheric current sheets for the pair enrichment and neutrino emission of AGN coronae. Our model has two free parameters, namely the proton plasma magnetization σ p , which controls the peak energy of the neutrino spectrum, and the Eddington ratio λ X,Edd (defined as the ratio between X-ray luminosity L X and Eddington luminosity L Edd ), which controls the amount of energy transferred to secondary particles. For sources with λ X,Edd ≳ λ Edd,crit (where λ Edd,crit ∼ 10 -1 for σ p = 10 5 or ∼ 10 -2 for σ p = 10 7 ), proton-photon interactions and γγ annihilation produce enough secondary pairs to achieve Thomson optical depths τ T ∼ 0.1-10. In the opposite case of λ X,Edd ≲ λ Edd,crit , the coronal pairs cannot originate only from hadronic interactions. Additionally, we find that the neutrino luminosity scales as L 2 X /L Edd for λ X,Edd ≲ λ Edd,crit , while it is proportional to L X for higher λ X,Edd values. We apply our model to four Seyfert galaxies, including NGC 1068, and discuss our results in light of recent IceCube observations.

active galactic nuclei↗
DOE OSTI2024

Assessing the impact of alpha particles on thermal confinement in JET D-T plasmas through global GENE-Tango simulations

The capability of the global, electromagnetic gyrokinetic GENE code interfaced with the transport Tango solver is exploited to address the impact of fusion alpha particles (in their dual role of fast particles and heating source) on plasma profiles and performance at JET in the discharges with the highest quasi-stationary peak fusion power during the DTE2 experimental campaigns. Employing radially global nonlinear electromagnetic GENE-Tango simulations, we compare results with/without alpha particles and alpha heating. Our findings reveal that alpha particles have a negligible impact on turbulent transport, with GENE-Tango converging to similar plasma profiles regardless of their inclusion as a kinetic species in GENE. On the other hand, alpha heating is found to contribute to the peaking of the electron temperature profiles, leading to a 1 keV drop on the on-axis electron temperature when alpha heating is neglected in Tango. The minimal impact of alpha particles on turbulent transport in this JET discharge–despite this being the shot with the highest fusion output–is attributed to the low content of fusion alpha in this discharge. To assess the potential impact of alpha particles on turbulent transport in regimes with higher alpha particle density, as expected in ITER and fusion reactors, we artificially increased the alpha particle concentration to levels expected for ITER. By performing global nonlinear GENE standalone simulations, we found that increasing the alpha particle density beyond five times the nominal value lead to significant overall turbulence destabilization. These results demonstrate that an increased alpha particle concentration can significantly impact transport properties under simulated JET experimental conditions. However, these findings cannot be directly extrapolated to ITER due to the substantial differences in parameters such as plasma size, magnetic field, plasma current, and thermal pressure.

energetic particles↗
DOE OSTI2021

Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds

Gyrokinetic theory is arguably the most important tool for numerical studies of transport physics in magnetized plasmas. However, exact local energy-momentum conservation laws for the electromagnetic gyrokinetic system have not been found despite the continuous effort. Without such local conservation laws, energy and momentum can be instantaneously transported across spacetime, which is unphysical and casts doubt on the validity of numerical simulations based on the gyrokinetic theory. The standard Noether procedure for deriving conservation laws from corresponding symmetries does not apply to gyrokinetic systems because the gyrocenters and electromagnetic fields reside on different manifolds. To overcome this difficulty, we develop a high-order field theory on heterogeneous manifolds for classical particle-field systems and apply it to derive exact, local conservation laws, in particular the energy-momentum conservation laws, for the electromagnetic gyrokinetic system. A weak Euler-Lagrange equation is established to replace the standard Euler-Lagrange equation for the particles. It is discovered that an induced weak Euler-Lagrange current enters the local conservation laws. And it is the new physics captured by the high-order field theory on heterogeneous manifolds. A recently developed gauge-symmetrization method for high-order electromagnetic field theories using the electromagnetic displacement-potential tensor is applied to render the derived energy-momentum conservation laws electromagnetic gauge-invariant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Natural hot-ion modes in a rotating plasma

We report that in steady state, the fuel cycle of a fusion plasma requires inward particle fluxes of fuel ions. These particle flows are also accompanied by heating. In the case of classical transport in a rotating cylindrical plasma, this heating can proceed through several distinct channels depending on the physical mechanisms involved. Some channels directly heat the fuel ions themselves, whereas others heat electrons. Which channel dominates depends, in general, on the details of the temperature, density and rotation profiles of the plasma constituents. However, remarkably, under relatively few assumptions concerning these profiles, if the α particles, the by-products of the fusion reaction, can be removed directly by other means, then a hot-ion mode tends to emerge naturally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Growth and Saturation of the Electron Drift Instability in a Crossed Field Plasma

The linear growth and nonlinear energy transfer of the electron drift instability (EDI) are experimentally measured in the plume of a low-temperature, Hall effect discharge. A frequency-based bispectral analysis technique applied to fast ion density fluctuation measurements shows a growth rate function that is qualitatively similar to predictions from the linear instability dispersion relation, but an order of magnitude smaller. Calculation of the nonlinear transfer function indicates multiple three-wave interactions between high-frequency resonances of the instability in addition to an inverse energy cascade toward lower-frequency modes. Furthermore, these results are discussed in the context of recent theoretical, numerical, and experimental efforts on the EDI in Hall effect discharges and how the EDI may impact anomalous cross field transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Topological modes and spectral flows in inhomogeneous PT-symmetric continuous media

In classical Hermitian continuous media, the spectral-flow index of topological modes is linked to the bulk topology via index theorem. However, the interface between two bulks is usually non-Hermitian due to the inhomogeneities of system parameters. We show that the connection between topological modes and bulk topology still exists despite the non-Hermiticity at the interface if the system is endowed with PT symmetry. The theoretical framework developed is applied to the Hall magnetohydrodynamic model to identify a topological mode called topological Alfvén sound wave in magnetized plasmas. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Measurement of the Alfvén Wave Parametric Decay Instability Growth Rate

Alfvén waves, a fundamental mode of magnetized plasmas, are ubiquitous in space and laboratory plasmas. The nonlinear behavior of these modes is thought to play a key role in important problems in space plasma, such as the heating of the solar corona and solar wind turbulence. In particular, theoretical predictions show that these Alfvén waves may be unstable to various parametric instabilities, but space observations of these processes are limited. We demonstrate the first measurement of the Alfvén wave parametric decay instability (PDI) growth rate. Experiments are conducted on the Large Plasma Device at UCLA in which a high amplitude 𝛿⁢𝐵/𝐵 0 ∼ 0.7% pump Alfvén wave is launched from one end of the device and a smaller seed Alfvén wave is launched from the other side. When the frequency of the seed wave is chosen to match the backward wave expected from PDI, damping of the seed wave is reduced. We compare this reduction in damping to the theoretically expected PDI growth rate while accounting for acoustic mode damping. Results show agreement between measurements and theoretical predictions. As a result, this not only provides critical validation for PDI theories and simulations that could help interpret future space observations but also suggests a new way of studying similar nonlinear wave phenomena.

Alfvén waves↗
DOE OSTI2023

Topological Langmuir-cyclotron wave

A theory is developed to describe the topological Langmuir-cyclotron wave (TLCW), a topological excitation in magnetized plasmas recently identified by numerical simulations. As a topological wave in a continuous medium, the TLCW propagates unidirectionally without scattering in complex boundaries and can be explored as an effective mechanism to energize particles. We show that, because momentum space in continuous media is contractible in general, the topology of the wave bundles is trivial over momentum space that contains no degeneracy points. This is in stark contrast to condensed matters with periodic lattice structures that impose nontrivial topology on momentum space. In continuous media without lattice structures, nontrivial topology of the eigenmode bundles manifests over phase space, and it is the nontrivial topology over phase space that underpins the topological excitations, such as the TLCW. It is shown that the TLCW can be faithfully modeled by a generic tilted Dirac cone in phase space, whose entire spectrum, including the spectral flow, is given.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗