SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Machine learning-based real-time kinetic profile reconstruction in DIII-D

Abstract Kinetic equilibrium reconstruction plays a vital role in the physical analysis of plasma stability and control in fusion tokamaks. However, the traditional approach is subjective and prone to human biases. To address this, the consistent automatic kinetic equilibrium reconstruction (CAKE) method was introduced, providing objective results. Nonetheless, its offline nature limits its application in real-time plasma control systems (PCSs). To address this limitation, we present RTCAKENN, a machine learning model that approximates 7 CAKE-level output profiles, namely pressure, inverse q , toroidal current density, electron temperature and density, carbon ion impurity temperature and rotation profiles, using real-time available inputs. The deep neural network consists of an encoder layer, where the scalars and interdependent inputs such as plasma boundary coordinates and motional Stark effect data are encoded using multi-layer perceptrons (MLPs), while profile inputs are encoded by 1D convolutional layers. The encoded data is passed through a MLP for latent feature extraction, before being decoded in the decoding layers, which consist of upsampling and convolutional layers. RTCAKENN has been implemented in the DIII-D PCS and our model achieves accuracy comparable to CAKE and surpasses existing real-time alternatives. Through clever dropout training, RTCAKENN exhibits robustness and can operate even in the absence of Thomson scattering data or charge exchange recombination data. It executes in under 8 ms in the real-time environment, enabling future application in real-time control and analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Numerical implications of including drifts in SOLPS-ITER simulations of EAST

The inclusion of drifts in plasma edge codes like SOLPS-ITER is required to match simulation data with experimental profiles. However, this remains numerically challenging. In this paper, the effect of some numerical factors on the final plasma solution is investigated. This study is performed on three EAST simulations in the upper single null configuration: an attached purely deuterium case, an attached case with limited Ne-seeding, and a detached Ne-seeded case. The effects of the anomalous conductivity and anomalous thermo-electric coefficient on the plasma potential are investigated. Next, the effect of the employed grids is shown. In order to investigate these effects, accurate drift simulations are needed. Therefore, the employed time step and numerical parameters are discussed for the three studied simulations. For all presented simulations, it is verified that the restriction of the grid to the first flux surface tangent to the main chamber wall is sufficient within the context of non-extended simulations. This means that the main power dissipation takes place inside the simulated domain, and only a small fraction of the power is leaving the B2.5 grid through the grid boundary closest to the first wall. Finally, the effect of drifts on the asymmetry between the inner and outer target for EAST simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Neural network based prediction of no-wall β N limits due to ideal external kink instabilities

Multi-layer perceptron based neural networks (NNs) are trained to predict the Troyon no-wall beta limits due to the onset of low-n (n = 1 – 3, n is the toroidal mode number) ideal external kink instabilities in tokamak plasmas. It is revealed that a well-trained NN can predict the n = 1 no-wall beta limit within 10% relative error. The NN performance is somewhat worse for the n = 2 and 3 no-wall beta limits, but still a relative error within 20% is achievable. The trained NNs well reproduce the known dependences of the no-wall beta limits on the plasma pressure (pressure peaking factor) and current (plasma internal inductance) profiles. Other scalings are further easily established with NNs, for parametric dependences such as on the aspect ratio, the elongation and triangularity of the plasma boundary shape. These NN-based Troyon beta limit predictors can be incorporated into integrated modeling platforms, or directly implemented as a real time stability estimator for the purpose of disruption avoidance or mitigation during experiments.

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↗

NASA NTRS1985

Principles of magnetohydrodynamic simulation in space plasmas

Attention is given to the philosophical as well as physical principles that are essential to the establishment of MHD simulation studies for solar plasma research, assuming the capabilities of state-of-the-art computers and emphasizing the importance of 'local' MHD simulation. Solar-terrestrial plasma space is divided into several elementary regions where a macroscopic elementary energy conversion process could conceivably occur; the local MHD simulation is defined as self-contained in each of the regions. The importance of, and the difficulties associated with, the boundary condition are discussed in detail. The roles of diagnostics and of the finite difference method are noted.

Sato, T.↗

DOE OSTI2020

Topological Gaseous Plasmon Polariton in Realistic Plasma

Nontrivial topology in bulk matter has been linked with the existence of topologically protected interfacial states. In this work, we show that a gaseous plasmon polariton (GPP), an electromagnetic surface wave existing at the boundary of magnetized plasma and vacuum, has a topological origin that arises from the nontrivial topology of magnetized plasma. Because a gaseous plasma cannot sustain a sharp interface with discontinuous density, one must consider a gradual density falloff with scale length comparable to or longer than the wavelength of the wave. We show that the GPP may be found within a gapped spectrum in present-day laboratory devices, suggesting that platforms are currently available for experimental investigation of topological wave physics in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Towards the optimization of direct laser acceleration

Experimental measurements using the OMEGA EP laser facility demonstrated direct laser acceleration (DLA) of electron beams to (505 ± 75) MeV with (140 ± 30) nC of charge from a low-density plasma target using a 400 J, picosecond duration pulse. Similar trends of electron energy with target density are also observed in self-consistent two-dimensional particle-in-cell simulations. The intensity of the laser pulse is sufficiently large that the electrons are rapidly expelled from along the laser pulse propagation axis to form a channel. The dominant acceleration mechanism is confirmed to be DLA and the effect of quasi-static channel fields on energetic electron dynamics is examined. A strong channel magnetic field, self-generated by the accelerated electrons, is found to play a comparable role to the transverse electric channel field in defining the boundary of electron motion.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Highest fusion performance without harmful edge energy bursts in tokamak

The path of tokamak fusion and International thermonuclear experimental reactor (ITER) is maintaining high-performance plasma to produce sufficient fusion power. This effort is hindered by the transient energy burst arising from the instabilities at the boundary of plasmas. Conventional 3D magnetic perturbations used to suppress these instabilities often degrade fusion performance and increase the risk of other instabilities. This study presents an innovative 3D field optimization approach that leverages machine learning and real-time adaptability to overcome these challenges. Implemented in the DIII-D and KSTAR tokamaks, this method has consistently achieved reactor-relevant core confinement and the highest fusion performance without triggering damaging bursts. This is enabled by advances in the physics understanding of self-organized transport in the plasma edge and machine learning techniques to optimize the 3D field spectrum. The success of automated, real-time adaptive control of such complex systems paves the way for maximizing fusion efficiency in ITER and beyond while minimizing damage to device components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1991

Extension of convection modeling into the high-latitude ionosphere - Some theoretical difficulties

The Rice Convection Model (RCM) is extended and merged with empirical models so as to cover the entire high-latitude ionosphere with the aim of providing precipitation and electric field inputs for ionosphere and thermosphere modelers and producing a model in which the boundaries of the precipitation and electric field patterns maintain physically consistent relationships to each other. The computed auroral electron energy flux, plotted as a function of latitude, exhibited an exaggerated two-peak structure. When no floor was placed under the precipitation rate, the minimum between the two peaks was much too deep to be consistent with typical observations. The regions of excessively weak precipitation map to equatorial distances of 15-35 RE and thus to the regions of the plasma sheet that were not included in previous self-consistent convection calculations.

Wolf, R. A.↗

DOE OSTI2022

Validation of the strike point position estimation with the local expansion method for MAST upgrade on the DIII-D tokamak

A local expansion method has been proposed for estimating the strike point position for the advanced divertor configuration on the MAST Upgrade tokamak. The paper discusses the application and assesses the performance of the technique on a long-legged divertor plasma configuration on an operating device - the DIII-D tokamak. A comparison of the spatial location of the outer strike point estimated with the local expansion method against the plasma boundary reconstruction and divertor diagnostics on the DIII-D tokamak is reported. A good agreement with the equilibrium reconstruction and diagnostic data is achieved with respect to estimation of the spatial location of the outer strike point for the long-legged divertor plasma discharge.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Narrow operator models of stellarator equilibria in Fourier Zernike basis

Numerical computation of the ideal magnetohydrodynamic (MHD) equilibrium magnetic field is at the base of stellarator optimisation and provides the starting point for solving more sophisticated partial differential equations like transport or turbulence models. Conventional approaches solve for a single stationary point of the ideal MHD equations, which is fully defined by three invariants and the numerical scheme employed by the solver. We present the first numerical approach that can solve for a continuous distribution of equilibria with fixed boundary and rotational transform, varying only the pressure invariant. This approach minimises the force residual by optimising parameters of multilayer perceptrons that map from a scalar pressure multiplier to the Fourier Zernike basis as implemented in the modern stellarator equilibrium solver DESC.

fusion plasma↗

DOE OSTI2022

Energetic particle-induced geodesic acoustic modes on DIII-D

Various properties of the energetic particle-induced geodesic acoustic mode (EGAM) are explored in this large database analysis of DIII-D experimental data. EGAMs are n = 0 modes with m = 0 electrostatic potential fluctuations (where n/m = toroidal/poloidal mode number), m = 1 density fluctuations, and m = 2 magnetic fluctuations. The fundamental frequency (~20–40 kHz) of the mode is typically below that of the traditional geodesic acoustic mode frequency. EGAMs are most easily destabilized by beams in the counter plasma current (counter-I p ) direction as compared to co-Ip and off-axis beams. During counter beam injection, the mode frequency is found to have the strongest linear dependence (correlation coefficient r = –0.71) with the safety factor (q). Here, the stability of the mode in the space of q and poloidal beta (β p ) shows a clear boundary for the mode stability. The stability of the mode depends more strongly on damping rate than on fast-ion drive for a given injection geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Time-Dependent Boundary Modeling to Inform Design of SPARC Diagnostic and Actuators

The purpose of this work was to inform the design of diagnostic and actuator systems for a new experimental facility for Commonwealth Fusion Systems (CFS), called SPARC [Creely2020], using time-dependent plasma boundary simulations with the Scrape Off Layer Plasma Simulator (SOLPS) code [Wiesen2015]. The SPARC tokamak is a device that aims to demonstrate fusion energy production and high-field plasma operating scenarios. The developed technologies will be incorporated into a net-electricity pilot that will demonstrate approximately 200MW of electric power. A major difference from other private industry fusion companies is the use of a standard aspect ratio, but at a high magnetic field strength, which provides a well-established physics basis and a modest extrapolation to reactor operation. Time-dependent simulations were used to inform the design of main-chamber plasma facing components, which are protected by feedback control of the heat and particle loads during high power operation. The simulations provide the timescales and magnitudes of plasma and neutral particle response useful for determining the position and operation of diagnostics and actuators that will enable control. The plasma response and phase-space diagrams also provide input into advanced model-based control schemes, which are envisioned for later SPARC operation. Model-based control can reduce the risk of component failure, which would require expensive in-situ repairs and the associated delays, by predicting the system response to actuators over a short time-horizon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

A new theory of sources of Birkeland currents

An approach to collisionless plasma shows the existence of current orthogonal to B along the low latitude boundary layer of the magnetosphere driven by electric field which is orthogonal to both B and the layer. In this case relationship P perpendicular + (sq B/8 pi) - (sq epsilon E perpendicular/8 pi) = constant holds on a line orthogonal to B and the layer, where epsilon is the dielectric constant of the plasma for electric fields orthogonal to B. Across the geomagnetic tail there flows a current in the direction of the dawn-dusk electric field, and in this case a relationship P perpendicular + (sq B/8 pi) + (sq epsilon E perpendicular/8 pi = constant, holds aong a line orthogonal to E and B. Divergence of both these currents is shown to be a source of Birkeland currents. Some of the boundary layer current is continuous with current across the tail. Electric currents of physically similar origin flow in interplanetary space, and when the magnetosphere interrupts them, additional Birkeland currents are driven.

Cole, K. D.↗

NASA NTRS1984

A new theory of sources of Birkeland currents

An approach to collisionless plasma shows the existence of current orthogonal to B along the low latitude boundary layer of the magnetosphere driven by electric field which is orthogonal to both B and the layer. In this case relationship P perpendicular + (sq B/8 pi) - (sq epsilon E perpendicular/8 pi) = constant holds on a line orthogonal to B and the layer, where epsilon is the dielectric constant of the plasma for electric fields orthogonal to B. Across the geomagnetic tail there flows a current in the direction of the dawn-dusk electric field, and in this case a relationship P perpendicular + (sq B/8 pi) + (sq epsilon E perpendicular/8 pi) = constant, holds along a line orthogonal to E and B. Divergence of both these currents is shown to be a source of Birkeland currents. Some of the boundary layer current is continuous with current across the tail. Electric currents of physically similar origin flow in interplanetary space, and when the magnetosphere interrupts them, additional Birkeland currents are driven.

Cole, K. D.↗

DOE OSTI2022

Eulerian simulations of electrostatic waves in plasmas with a single sign of charge

An Eulerian, numerical simulation is used to model the launching of plasma waves in a non-neutral plasma that is confined in a Penning–Malmberg trap. The waves are launched by applying an oscillating potential to an electrically isolated sector at one end of the conducting cylinder that bounds the confinement region and are received by another electrically isolated sector at the other end of the cylinder. The launching of both Trivelpiece–Gould waves and electron acoustic waves is investigated. Adopting a stratagem, the simulation captures essential features of the finite length plasma, while retaining the numerical advantages of a simulation employing periodic spatial boundary conditions. As a benchmark test of the simulation, the results for launched Trivelpiece–Gould waves of small amplitude are successfully compared to a linearized analytic solution for these fluctuations.

Physics↗

NASA NTRS2013

Analytic Theory of Titans Schumann Resonance: Constraints on Ionospheric Conductivity and Buried Water Ocean

This study presents an approximate model for the atypical Schumann resonance in Titan's atmosphere that accounts for the observations of electromagnetic waves and the measurements of atmospheric conductivity performed with the Huygens Atmospheric Structure and Permittivity, Wave and Altimetry (HASI-PWA) instrumentation during the descent of the Huygens Probe through Titan's atmosphere in January 2005. After many years of thorough analyses of the collected data, several arguments enable us to claim that the Extremely Low Frequency (ELF) wave observed at around 36 Hz displays all the characteristics of the second harmonic of a Schumann resonance. On Earth, this phenomenon is well known to be triggered by lightning activity. Given the lack of evidence of any thunderstorm activity on Titan, we proposed in early works a model based on an alternative powering mechanism involving the electric current sheets induced in Titan's ionosphere by the Saturn's magnetospheric plasma flow. The present study is a further step in improving the initial model and corroborating our preliminary assessments. We first develop an analytic theory of the guided modes that appear to be the most suitable for sustaining Schumann resonances in Titan's atmosphere. We then introduce the characteristics of the Huygens electric field measurements in the equations, in order to constrain the physical parameters of the resonating cavity. The latter is assumed to be made of different structures distributed between an upper boundary, presumably made of a succession of thin ionized layers of stratospheric aerosols spread up to 150 km and a lower quasi-perfect conductive surface hidden beneath the non-conductive ground. The inner reflecting boundary is proposed to be a buried water-ammonia ocean lying at a likely depth of 55-80 km below a dielectric icy crust. Such estimate is found to comply with models suggesting that the internal heat could be transferred upwards by thermal conduction of the crust, while convective processes cannot be ruled out.

TITAN↗

NASA NTRS2014

A Vortical Dawn Flank Boundary Layer for Near-Radial IMF: Wind Observations on 24 October 2001

We present an example of a boundary layer tailward of the dawn terminator which is entirely populated by rolled-up flow vortices. Observations were made by Wind on 24 October 2001 as the spacecraft moved across the region at the X plane approximately equal to −13 Earth radii. Interplanetary conditions were steady with a near-radial interplanetary magnetic field (IMF). Approximately 15 vortices were observed over the 1.5 hours duration of Wind's crossing, each lasting approximately 5 min. The rolling up is inferred from the presence of a hot tenuous plasma being accelerated to speeds higher than in the adjoining magnetosheath, a circumstance which has been shown to be a reliable signature of this in single-spacecraft observations. A blob of cold dense plasma was entrained in each vortex, at whose leading edge abrupt polarity changes of field and velocity components at current sheets were regularly observed. In the frame of the average boundary layer velocity, the dense blobs were moving predominantly sunward and their scale size along the X plane was approximately 7.4 Earth radii. Inquiring into the generation mechanism of the vortices, we analyze the stability of the boundary layer to sheared flows using compressible magnetohydrodynamic Kelvin-Helmholtz theory with continuous profiles for the physical quantities. We input parameters from (i) the exact theory of magnetosheath flow under aligned solar wind field and flow vectors near the terminator and (ii) the Wind data. It is shown that the configuration is indeed Kelvin-Helmholtz (KH) unstable. This is the first reported example of KH-unstable waves at the magnetopause under a radial IMF.

Magnetosheath↗