Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA CONFINEMENT”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Comparison of different interpolation operators including nonlinear subdivision schemes in the simulation of particle trajectories

In this work, we compare different interpolation operators in the context of particle tracking with an emphasis on situations involving velocity field with steep gradients. Since, in this case, most classical methods give rise to the Gibbs phenomenon (generation of oscillations near discontinuities), we present new methods for particle tracking based on subdivision schemes and especially on the Piecewise Parabolic Harmonic (PPH) scheme which has shown its advantage in image processing in presence of strong contrasts. First an analytic univariate case with a discontinuous velocity field is considered in order to highlight the effect of the Gibbs phenomenon on trajectory calculation. Theoretical results are provided. Then, we show, regardless of the interpolation method, the need to use a conservative approach when integrating a conservative problem with a velocity field deriving from a potential. Finally, the PPH scheme is applied in a more realistic case of a time-dependent potential encountered in the edge turbulence of magnetically confined plasmas, to compare the propagation of density structures (turbulence bursts) with the dynamics of test particles. This study highlights the difference between particle transport and density transport in turbulent fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Progress in reducing ICRF-specific impurity release in ASDEX upgrade and JET

Use of new 3-strap ICRF antennas with all-tungsten (W) limiters in ASDEX Upgrade results in a reduction of the W sources at the antenna limiters and of the W content in the confined plasma by at least a factor of 2 compared to the W-limiter 2-strap antennas used in the past. The reduction is observed with a broad range of plasma shapes. In multiple locations of antenna frame, the limiter W source has a minimum when RF image currents are decreased by cancellation of the RF current contributions of the central and the outer straps. In JET with ITER-like wall, ITER-like antenna produces about 20% less of main chamber radiation and of W content compared to the old A2 antennas. However the effect of the A2 antennas on W content is scattered depending on which antennas are powered. Experiments in JET with trace nitrogen (N 2 ) injection show that a presence of active ICRF antenna close to the midplane injection valve has little effect on the core N content, both in dipole and in -90° phasing. This indicates that the effect of ICRF on impurity transport across the scape-off-layer is small in JET compared to the dominant effect on impurity sources leading to increased impurity levels during ICRF operation.

1CRF↗
DOE OSTI2018

Experimental study of the stability and dynamics of a two-dimensional ideal vortex under external strain

The dynamics of two-dimensional (2-D) ideal fluid vortices is studied experimentally in the presence of an irrotational strain flow. Laboratory experiments are conducted using strongly magnetized pure electron plasmas, a technique which is made possible by the isomorphism between the drift–Poisson equations describing plasma dynamics transverse to the field and the 2-D Euler equations describing an ideal fluid. The electron plasma system provides an excellent opportunity to study the dynamics of a 2-D Euler fluid due to weak dissipation and weak 3-D effects, simple diagnosis and precise control. The plasma confinement apparatus used here was designed specifically to study vortex dynamics under the influence of external flow by applying boundary conditions in two dimensions. Additionally, vortex-in-cell simulations are carried out to complement the experimental results and to extend the parameter range of the studies. It is shown that the global dynamics of a quasi-flat vorticity profile is in good quantitative agreement with the theory of a piecewise-constant elliptical patch of vorticity, including the equilibria, dynamical orbits and stability properties. Deviations from the elliptical patch theory are observed for non-flat vorticity profiles; they include inviscid damping of the orbits and modified stability limits. The dependence of these phenomena on the flatness of the initial profile is discussed. The relationship of these results to other theoretical, numerical and experimental studies is also discussed.

Mechanics↗
DOE OSTI2017

A drift-kinetic analytical model for scrape-off layer plasma dynamics at arbitrary collisionality

We report a drift-kinetic model to describe the plasma dynamics in the scrape-off layer region of tokamak devices at arbitrary collisionality is derived. Our formulation is based on a gyroaveraged Lagrangian description of the charged particle motion, and the corresponding drift-kinetic Boltzmann equation that includes a full Coulomb collision operator. Using a Hermite–Laguerre velocity space decomposition of the gyroaveraged distribution function, a set of equations to evolve the coefficients of the expansion is presented. By evaluating explicitly the moments of the Coulomb collision operator, distribution functions arbitrarily far from equilibrium can be studied at arbitrary collisionalities. A fluid closure in the high-collisionality limit is presented, and the corresponding fluid equations are compared with previously derived fluid models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Finite Larmor radius effects on E x B weak turbulence transport

Transport of test particles in two-dimensional weak turbulence with waves propagating along the poloidal direction is studied using a reduced model. Finite Larmor radius (FLR) effects are included by gyroaveraging over one particle orbit. For low wave amplitudes the motion is mostly regular with particles trapped in the potential wells. As the amplitude increases the trajectories become chaotic and the Larmor radius modifies the orbits. For a thermal distribution of Finite Larmor radii the particle distribution function (PDF) is Gaussian for small $\unicode[STIX]{x1D70C}_{th}$ (thermal gyroradius) but becomes non-Gaussian for large $\unicode[STIX]{x1D70C}_{th}$ . Furthermore, the time scaling of transport is diffusive, as characterized by a linear dependence of the variance of the PDF with time. An explanation for this behaviour is presented that provides an expression for an effective diffusion coefficient and reproduces the numerical results for large wave amplitudes which implies generalized chaos. When a shear flow is added in the direction of wave propagation, a modified model is obtained that produces free-streaming particle trajectories in addition to trapped ones; these contribute to ballistic transport for low wave amplitude but produce super-ballistic transport in the chaotic regime. As in the previous case, the PDF is Gaussian for low $\unicode[STIX]{x1D70C}_{th}$ becoming non-Gaussian as it increases. The perpendicular transport presents the same behaviour as in the case with no flow but the diffusion is faster in the presence of the flow.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Practical gyrokinetics

We report thousands of gyrokinetic papers have been published since the introduction of the gyrokinetic change of variables. The intent here is not to review the field, but rather the goal is to present a tutorial providing insight into why gyrokinetics is an appropriate description of turbulent transport. The focus is on turbulent transport in axisymmetric tokamaks, but many of the ideas and techniques are applicable to stellarators and other magnetic fusion devices with nested surfaces of constant pressure. Besides the origins of gyrokinetics and recent insights, gyrokinetic orderings and gyrokinetic variables are summarized. Then a compact, but careful, derivation of the simplest electrostatic gyrokinetic equation for tokamaks is presented, along with a brief mention of gyrokinetics for stellarators. The advantages of assuming scale separation between the finer spatial scales and faster time variation of the turbulence and the global behaviour and slow temporal evolution of the background are stressed. Moreover, the procedure for removing the adiabatic or Maxwell–Boltzmann response is emphasized. Scale separation allows the near Maxwellian behaviour of the background and the rapid variation of the turbulence to be described by separate, local gyrokinetic equations. The turbulent fluctuations are found by solving the nonlinear gyrokinetic equation for the non-adiabatic portion of the fluctuating distribution function. Also, generalizations of the gyrokinetic equation so that it retains electromagnetic and flow modifications are considered in detail along with some symmetry properties. In addition, ambipolarity, particle transport and heat transport are addressed, along with a discussion of the complications associated with describing momentum transport and profile evolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Tokamak parameter regimes with low toroidal momentum diffusivity

The strong, sharp flow structures that are seen frequently in tokamak cores, and large amplitude spontaneous global toroidal rotation are both surprising in light of current theories where toroidal flow evolution is dominantly diffusive. Mechanisms for spontaneously generating strong poloidal shear flows have been extensively investigated, but these processes were thought not to apply to toroidal flows. We confirm, however, that there is a regime with near-zero toroidal momentum diffusivity, where toroidal flow structures are spontaneously generated, as shown in earlier global gyrokinetic simulations. This allows strong rotation with negligible applied external torque, and the regime where this occurs is also favourable for strong turbulence stabilisation by rotation. In conclusion, the transition to low momentum diffusivity occurs for tight aspect ratio, low safety factor and at low levels of heat flux, in agreement with a simple zero-dimensional model for the flow–turbulence interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Direct construction of optimized stellarator shapes. Part 3. Omnigenity near the magnetic axis

>The condition of omnigenity is investigated, and applied to the near-axis expansion of Garren & Boozer ( Phys. Fluids B, vol. 3 (10), 1991 a , pp. 2805–2821). Due in part to the particular analyticity requirements of the near-axis expansion, we find that, excluding quasi-symmetric solutions, only one type of omnigenity, namely quasi-isodynamicity, can be satisfied at first order in the distance from the magnetic axis. Our construction provides a parameterization of the space of such solutions, and the cylindrical reformulation and numerical method of Landreman & Sengupta ( J. Plasma Phys. , vol. 84 (6), 2018, 905840616); Landreman et al. ( J. Plasma Phys. , vol. 85 (1), 2019, 905850103), enables their efficient numerical construction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with ion cyclotron range-of-frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach $\textit{Q}$ ≈ 11, which is well above its mission objective of $\textit{Q}$ > 2. To build confidence that SPARC will be successful, physics-based integrated modelling has also been performed. The TRANSP code coupled with the theory-based trapped gyro-Landau fluid (TGLF) turbulence model and EPED predictions for pedestal stability find that $\textit{Q}$ ≈ 9 is attainable in standard H-mode operation and confirms $\textit{Q}$ > 2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modelled with the Monte–Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fast-ion physics in SPARC

Potential loss of energetic ions including alphas and radio-frequency tail ions due to classical orbit effects and magnetohydrodynamic instabilities (MHD) are central physics issues in the design and experimental physics programme of the SPARC tokamak. The expected loss of fusion alpha power due to ripple-induced transport is computed for the SPARC tokamak design by the ASCOT and SPIRAL orbit-simulation codes, to assess the expected surface heating of plasma-facing components. We find good agreement between the ASCOT and SPIRAL simulation results not only in integrated quantities (fraction of alpha power loss) but also in the spatial, temporal and pitch-angle dependence of the losses. If the toroidal field (TF) coils are well-aligned, the SPARC edge ripple is small (0.15–0.30 %), the computed ripple-induced alpha power loss is small ( ~0.25% ) and the corresponding peak surface power density is acceptable ( 244 kW m -2 ). However, the ripple and ripple-induced losses increase strongly if the TF coils are assumed to suffer increasing magnitudes of misalignment. Surface heat loads may become problematic if the TF coil misalignment approaches the centimetre level. Ripple-induced losses of the energetic ion tail driven by ion cyclotron range of frequency (ICRF) heating are not expected to generate significant wall or limiter heating in the nominal SPARC plasma scenario. Because the expected classical fast-ion losses are small, SPARC will be able to observe and study fast-ion redistribution due to MHD including sawteeth and Alfvén eigenmodes (AEs). SPARC's parameter space for AE physics even at moderate Q is shown to reasonably overlap that of the demonstration power plant ARC (Sorbom et al., Fusion Engng Des., vol. 100, 2015, p. 378), and thus measurements of AE mode amplitude, spectrum and associated fast-ion transport in SPARC would provide relevant guidance about AE behaviour expected in ARC.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An Alfvenic reconnecting plasmoid thruster

A new concept for the generation of thrust for space propulsion is introduced. Energetic thrust is generated in the form of plasmoids (confined plasma in closed magnetic loops) when magnetic helicity (linked magnetic field lines) is injected into an annular channel. Using a novel configuration of static electric and magnetic fields, the concept utilizes a current-sheet instability to spontaneously and continuously create plasmoids via magnetic reconnection. Here, the generated low-temperature plasma is simulated in a global annular geometry using the extended magnetohydrodynamic model. Because the system-size plasmoid is an Alfvenic outflow from the reconnection site, its thrust is proportional to the square of the magnetic field strength and does not ideally depend on the mass of the ion species of the plasma. Exhaust velocities in the range of 20 to $500\ \textrm {km}\ \textrm {s}^{-1}$ , controllable by the coil currents, are observed in the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ion temperature gradient (ITG), trapped electron mode (TEM) and electron temperature gradient (ETG) modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. Finally, it is also confirmed that thermodiffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Adjoint methods for quasi-symmetry of vacuum fields on a surface

Adjoint methods can speed up stellarator optimisation by providing gradient information more efficiently compared with finite-difference evaluations. Adjoint methods are herein applied to vacuum magnetic fields, with objective functions targeting quasi-symmetry and a rotational transform value on a surface. Here to measure quasi-symmetry, a novel way of evaluating approximate flux coordinates on a single flux surface without the assumption of a neighbourhood of flux surfaces is proposed. The shape gradients obtained from the adjoint formalism are evaluated numerically and verified against finite-difference evaluations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Modelling the nonlinear plasma response to externally applied three-dimensional fields with the Stepped Pressure Equilibrium Code

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma image classification using cosine similarity constrained convolutional neural network

Plasma jets are widely investigated both in the laboratory and in nature. Astrophysical objects such as black holes, active galactic nuclei and young stellar objects commonly emit plasma jets in various forms. With the availability of data from plasma jet experiments resembling astrophysical plasma jets, classification of such data would potentially aid in not only investigating the underlying physics of the experiments but also the study of astrophysical jets. In this work we use deep learning to process all of the laboratory plasma images from the Caltech Spheromak Experiment spanning two decades. We found that cosine similarity can aid in feature selection, classify images through comparison of feature vector direction and be used as a loss function for the training of AlexNet for plasma image classification. We also develop a simple vector direction comparison algorithm for binary and multi-class classification. Using our algorithm we demonstrate 93 % accurate binary classification to distinguish unstable columns from stable columns and 92 % accurate five-way classification of a small, labelled data set which includes three classes corresponding to varying levels of kink instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The DESC stellarator code suite Part 3: Quasi-symmetry optimization

The DESC stellarator optimization code takes advantage of advanced numerical methods to search the full parameter space much faster than conventional tools. Only a single equilibrium solution is needed at each optimization step thanks to automatic differentiation, which efficiently provides exact derivative information. A Gauss–Newton trust-region optimization method uses second-order derivative information to take large steps in parameter space and converges rapidly. With just-in-time compilation and GPU portability, high-dimensional stellarator optimization runs take orders of magnitude less computation time with DESC compared to other approaches. This paper presents the theory of the DESC fixed-boundary local optimization algorithm along with demonstrations of how to easily implement it in the code. Example quasi-symmetry optimizations are shown and compared to results from conventional tools. Three different forms of quasi-symmetry objectives are available in DESC, and their relative advantages are discussed in detail. In the examples presented, the triple product formulation yields the best optimization results in terms of minimized computation time and particle transport. This paper concludes with an explanation of how the modular code suite can be extended to accommodate other types of optimization problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The DESC stellarator code suite. Part 1. Quick and accurate equilibria computations

Three-dimensional equilibrium codes are vital for stellarator design and operation, and high-accuracy equilibria are also necessary for stability studies. This paper details comparisons of two three-dimensional equilibrium codes: VMEC, which uses a steepest-descent algorithm to reach a minimum-energy plasma state, and DESC, which minimizes the magnetohydrodynamic (MHD) force error in real space directly. Accuracy as measured by satisfaction of MHD force balance is presented for each code, along with the computation time. It is shown that DESC is able to achieve more accurate solutions, especially near axis. The importance of higher-accuracy equilibria is shown in DESC's better agreement of stability metrics with asymptotic formulae. DESC's global Fourier–Zernike basis also yields solutions with analytic derivatives explicitly everywhere in the plasma volume, provides improved accuracy in the radial direction versus conventional finite differences and allows for exponential convergence. Further, DESC can compute a solution with the same accuracy as VMEC in order-of-magnitude less time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗