Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “toroidal rotation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Linear simulation of magnetohydrodynamic plasma response to three-dimensional magnetic perturbations in high-β P plasmas

In this work, we report the numerical analyses of linear magnetohydrodynamics (MHD) plasma response to applied three-dimensional magnetic perturbations (MPs) in a joint DIII-D/EAST collaboration on high-β P (poloidal beta) plasmas, utilizing the extended-MHD code M3D-C1, with the purpose of realizing a better understanding of the existing experiment in which the n=3 MPs were applied to such high-β P plasmas attempting to control large amplitude type-I ELMs. Such high-β P plasmas obtained at the DIII-D tokamak feature an upper-biased double null configuration, a high edge safety factor q 95 ~6.4, and a stable internal transport barrier (ITB) leading to relatively high core pressures. Single-fluid simulations show that the plasma response to n=3 MPs, including both non-resonant/kinking and resonant components, is significantly weaker than that to n=1 or 2 MPs. To survey the impact of q 95 on plasma response to applied MPs, the SEGWAY (Self-consistent Equilibrium Generating Workflow for AnalYsis) module, developed in the OMFIT integrated modelling framework, is employed to generate a series of equilibria with a wide range of q 95 while other key parameters including the normalized beta, electron density at pedestal top, and plasma shape are kept fixed. Compared to the vacuum response, single-fluid M3D-C1 simulations predict a much more significant decrease of resonant plasma response to the applied n=3 MPs at the maximum penetration radii as q 95 increases. In contrast to single-fluid simulation results showing resonant penetration occurs only near the pedestal top where the E×B toroidal rotation frequency is zero, two-fluid simulations show two comparable resonant penetrations locating near the pedestal top and the ITB foot, where the perpendicular electron rotation frequency is zero. Such resonant field penetration near the ITB foot may be responsible for the observed formation of a staircase structure in both electron density and temperature profiles and thereby a considerable deterioration of global plasma performance when MPs are applied in high-β P plasmas. Motivated by this numerical work, we provide some ideas for the future research, with the purpose of realizing effective ELM control in such high-β P plasmas on the DIII-D and EAST devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Integrated control of individual plasma scalars with simultaneous neoclassical tearing-mode suppression

A novel integrated-control architecture has been tested in nonlinear, one-dimensional simulations using the control-oriented transport simulator (COTSIM©) and in DIII-D experiments. Integrated architectures that can perform continuous-mission control while also handling off-normal events will be vital in future reactor-grade tokamaks. Continuous-mission controllers for individual magnetic and kinetic scalars (thermal stored-energy (W), volume-average toroidal rotation (Ω Φ ), and safety factor profile (q) at different spatial locations) have been integrated in this work with event-triggered neoclassical tearing-mode (NTM) suppression controllers by combining them into an architecture augmented by a supervisory and exception handling (S&EH) system and an actuator management (AM) system. Here, the AM system, which enables the integration of competing controllers, solves in real time a nonlinear optimization problem that takes into account the high-level control priorities dictated by the S&EH system. The resulting architecture offers a high level of integration and some of the functionalities that will be required to fulfill the advanced-control requirements anticipated for ITER. Initial simulations using COTSIM suggest that the plasma performance and its MHD stability may be improved under integrated feedback control. In addition, the integrated-control architecture has been implemented in the DIII-D plasma control system and tested experimentally for the first time ever in DIII-D in a high-q min scenario, which is a candidate for steady-state operation in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Overview of physics results from the ADITYA-U tokamak and future experiments

The ADITYA upgrade (ADITYA-U), a medium-sized (R 0 = 75 cm, a = 25 cm) conventional tokamak facility in India, has been consistently producing experiments findings by usingcircular and shaped-plasmas. Recognizing the plasma parameters aligning closely with the design parameters of circular limited plasmas, ADITYA-U shifted its focus toward exploring the operational regime for experimentation on saw-tooth and MHD phenomena. Moreover, ADITYA-U has made consistent advancements toward conducting preliminary plasma shaping experiments through the activation of top and bottom divertor coils utilizing hydrogen as well as deuterium fuels. Confinement is improved by a factor of ~1.5 in D 2 plasmas when compared to H 2 plasmas of ADITYA-U. Further, ADITYA-U operations emphasize preventing disruptions and runaway electrons (REs) to ensure safe operations for future fusion devices. Significant suppression of REs has been achieved in ADITYA-U with the application of pulsed localized vertical magnetic field (LVF) perturbation, thereby establishing the technique's independence from the tokamak device. The successful RE mitigation requires a critical threshold of LVF pulse magnitude, which is approximately 1% of the toroidal magnetic field, and a minimum duration of ~5 ms. Apart from this, several novel findings have been achieved in the ADITYA-U experiments, including the modification of sawtooth duration through gas-puff, the emergence of MHD-induced geodesic acoustic mode-like oscillations, the propagation of fast heat pulses induced by MHD activity, the control of RE dynamics through Gas-puffs, the propagation of pinch-driven cold-pulses, the transport and core accumulations of argon impurities, the mass dependency of plasma toroidal rotation and the detection of 'RICE' scaling, as well as the characterization of edge plasma using wall conditioning methods, such as glow discharge cleaning using a combination of Ar-H 2 mixture, localized wall cleaning by electron cyclotron resonant plasma, and the development of machine learning-based disruption predictions, will be discussed in this paper.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Nonlinear gyrokinetic modelling of high confinement negative triangularity plasmas

Abstract Nonlinear gyrokinetic simulations correctly predict particle as well as ion and electron energy fluxes of high confinement plasmas with a negative triangularity cross sectional shape, showing that core transport in these plasmas is well described by standard gyrokinetic models. Experimentally inferred power balance fluxes are mostly reproduced within one standard deviation across a wide portion of the minor radius. Experimental conditions are reproduced by ion scale simulations, without the need to include density and temperature profile curvature effects. The experimental case is used as baseline to predict that the non-dimensional confinement scaling in negative triangularity plasmas increases strongly with plasma current while slightly degrading at increasing normalized pressure and decreasing collisionality. Recent experiments showed that low toroidal rotation negatively impacts confinement; consistent with the experiment, simulations predict that low rotational shear significantly affects confinement unless the plasma effective charge is maintained above a minimum level. Core confinement is predicted to significantly degrade in low aspect ratio devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The Madison plasma dynamo experiment: A facility for studying laboratory plasma astrophysics

The Madison plasma dynamo experiment (MPDX) is a novel, versatile, basic plasma research device designed to investigate flow driven magnetohydrodynamic instabilities and other high-β phenomena with astrophysically relevant parameters. A 3 m diameter vacuum vessel is lined with 36 rings of alternately oriented 4000 G samarium cobalt magnets, which create an axisymmetric multicusp that contains ∼14 m3of nearly magnetic field free plasma that is well confined and highly ionized (>50%). At present, 8 lanthanum hexaboride (LaB6) cathodes and 10 molybdenum anodes are inserted into the vessel and biased up to 500 V, drawing 40 A each cathode, ionizing a low pressure Ar or He fill gas and heating it. Up to 100 kW of electron cyclotron heating power is planned for additional electron heating. The LaB6cathodes are positioned in the magnetized edge to drive toroidal rotation through J × B torques that propagate into the unmagnetized core plasma. Dynamo studies on MPDX require a high magnetic Reynolds number Rm > 1000, and an adjustable fluid Reynolds number 10 1). Initial results from MPDX are presented along with a 0-dimensional power and particle balance model to predict the viscosity and resistivity to achieve dynamo action.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Overview of the TCV tokamak program: Scientific progress and facility upgrades

The TCV tokamak is augmenting its unique historical capabilities (strong shaping, strong electron heating) with ion heating, additional electron heating compatible with high densities, and variable divertor geometry, in a multifaceted upgrade program designed to broaden its operational range without sacrificing its fundamental flexibility. The TCV program is rooted in a three-pronged approach aimed at ITER support, explorations towards DEMO, and fundamental research. A 1 MW, tangential neutral beam injector (NBI) was recently installed and promptly extended the TCV parameter range, with record ion temperatures and toroidal rotation velocities and measurable neutral-beam current drive. ITER-relevant scenario development has received particular attention, with strategies aimed at maximizing performance through optimized discharge trajectories to avoid MHD instabilities, such as peeling-ballooning and neoclassical tearing modes. Experiments on exhaust physics have focused particularly on detachment, a necessary step to a DEMO reactor, in a comprehensive set of conventional and advanced divertor concepts. The specific theoretical prediction of an enhanced radiation region between the two X-points in the low-field-side snowflake-minus configuration was experimentally confirmed. Fundamental investigations of the power decay length in the scrape-off layer (SOL) are progressing rapidly, again in widely varying configurations and in both D and He plasmas; in particular, the double decay length in L-mode limited plasmas was found to be replaced by a single length at high SOL resistivity. Experiments on disruption mitigation by massive gas injection and electron-cyclotron resonance heating (ECRH) have begun in earnest, in parallel with studies of runaway electron generation and control, in both stable and disruptive conditions; a quiescent runaway beam carrying the entire electrical current appears to develop in some cases. Developments in plasma control have benefited from progress in individual controller design and have evolved steadily towards controller integration, mostly within an environment supervised by a tokamak profile control simulator. TCV has demonstrated effective wall conditioning with ECRH in He in support of the preparations for JT-60SA operation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Perturbed Ion Temperature and Toroidal Flow Profile Measurements in Rotating Neoclassical Tearing Mode Magnetic Islands

The perturbed ion temperature and toroidal flow were measured in rotating Neoclassical Tearing Modes (NTM) in a tokamak for the first time. These toroidally and radially resolved profiles were obtained by impurity ion spectroscopy in a 2,1 NTM in DIII-D. In agreement with drift-kinetic simulations, the electron temperature profile is flat, while the ion temperature gradient is restored across the magnetic island O-point in the presence of fast ions; the perturbed flow has minima in the O-points and maxima at the X-points. Furthermore, these measurements provide the first confirmation of the theoretically expected ion temperature and flow response to a magnetic island needed to predict the NTM on set threshold scaling for ITER and other future tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Toroidal polar topology in strained ferroelectric polymer

Polar topological texture has become an emerging research field for exotic phenomena and potential applications in reconfigurable electronic devices. Here, we report toroidal topological texture self-organized in a ferroelectric polymer, poly(vinylidene fluoride-ran-trifluoroethylene) [P(VDF-TrFE)], that exhibits concentric topology with anticoupled chiral domains. The interplay among the elastic, electric, and gradient energies results in continuous rotation and toroidal assembly of the polarization perpendicular to polymer chains, whereas relaxor behavior is induced along polymer chains. Such toroidal polar topology gives rise to periodic absorption of polarized far-infrared (FIR) waves, enabling the manipulation of the terahertz wave on a mesoscopic scale. Our observations should inform design principles for flexible ferroic materials toward complex topologies and provide opportunities for multistimuli conversions in flexible electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analysis of MHD stability and active mode control on KSTAR for high confinement, disruption-free plasma

Long-pulse plasma operation at high normalized beta, $β_N$, above the $\textit{n}$= 1 ideal MHD no-wall stability limit in KSTAR is presently limited by tearing instabilities rather than resistive wall modes. H-mode plasma operation during the recent KSTAR device campaign produced discharges having strong $\textit{m/n}$= 2/1 tearing instabilities at $β_N$ lower than the ideal MHD no-wall beta limit. The unstable tearing mode consequently reduced plasma confinement and toroidal plasma rotation significantly. We report the experiment confirmed that an extended duration of electron cyclotron heating (ECH) at the initial phase of the discharge plays a critical role in mode destabilization. To study destabilizing mechanisms that affect the mode growth, the stability of the observed tearing modes from plasmas with significantly different $β_N$ is computed by using the resistive DCON code and the M3D-C 1 code. Equilibrium reconstructions that include constraints from internal profile diagnostics, and computed fast particle pressure are used as input for reliable computation of stability. The classical tearing stability index, Δ', from resistive DCON is computed to be unstable when the island is fully saturated with large amplitude, while the unstable mode is not computed by M3D-C 1 . The modified Rutherford equation (MRE) describing the evolution of the island width has been constructed for KSTAR plasmas by using plasma parameters computed by the TRANSP code. The MRE model estimates a saturated island width corresponding to ~10% of the plasma minor radius for equilibrium at high $β_N$ 3 having a stable Δ' from the resistive DCON. In preparation for long-pulse plasma operation at higher beta utilizing increased plasma heating power, a resistive wall mode (RWM) active feedback control algorithm has been completed and enabled on KSTAR. To accurately determine the $\textit{n}$= 1component produced by RWMs, an algorithm has been developed that includes magnetic sensor compensation of the prompt applied field and the field from the induced current on the passive conductors. Use of multiple toroidal sensor arrays is enabled by modifying the sensor toroidal angles assumed in mode decomposition to include the effect of varied mode helicities in the outboard region where the mode measurement is made. This analysis on stability, transport, and control provides the required foundation for disruption prediction and avoidance research on KSTAR.

resistive wall mode (RWM)↗
DOE OSTI2023

Experimental study of the edge radial electric field in different drift configurations and its role in the access to H-mode at ASDEX Upgrade

The formation of the equilibrium radial electric field (Er) has been studied experimentally at ASDEX Upgrade (AUG) in L-modes of “favorable” (ion ∇ B-drift toward primary X-point) and “unfavorable” (ion ∇ B-drift away from primary X-point) drift configurations, in view of its impact on H-mode access, which changes with drift configurations. Edge electron and ion kinetic profiles and impurity velocity and mean-field Er profiles across the separatrix are investigated, employing new and improved measurement techniques. The experimental results are compared to local neoclassical theory as well as to a simple 1D scrape-off layer (SOL) model. It is found that in L-modes of matched heating power and plasma density, the upstream SOL Er and the main ion pressure gradient in the plasma edge are the same for either drift configurations, whereas the Er well in the confined plasma is shallower in unfavorable compared to the favorable drift configuration. The contributions of toroidal and poloidal main ion flows to Er, which are inferred from local neoclassical theory and the experiment, cannot account for these observed differences. Furthermore, it is found that in the L-mode, the intrinsic toroidal edge rotation decreases with increasing collisionality and it is co-current in the banana-plateau regime for all different drift configurations at AUG. This gives rise to a possible interaction of parallel Pfirsch–Schlüter flows in the SOL with the confined plasma. Thus, the different H-mode power threshold for the two drift configurations cannot be explained in the same way at AUG as suggested by LaBombard et al. [Phys. Plasmas 12, 056111 (2005)] for Alcator C-Mod. Finally, comparisons of Er profiles in favorable and unfavorable drift configurations at the respective confinement transitions show that also the Er gradients are all different, which indirectly indicates a different type or strength of the characteristic edge turbulence in the two drift configurations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Investigating nonlinear magnetohydrodynamics in an optimized, reactor-scale quasi-axisymmetric stellarator

We use the new simulation capabilities of the extended-magnetohydrodynamic (MHD) code, M3D-C1, to investigate the nonlinear MHD properties of a reactor-scale quasisymmetric stellarator equilibrium. Our model captures the self-consistent evolution of the magnetic field, temperature, density, and flow profiles without imposing restrictions on the structure of the first. We include the effects of resistivity using a realistic temperature-dependent Spitzer model, along with a model for heat transport that captures the key physical characteristic, namely, strongly anisotropic diffusion in directions perpendicular and parallel to the magnetic field. We consider a quasi-axisymmetric, finite-pressure equilibrium that was optimized for self-consistent bootstrap current, quasi-symmetry, and energetic particle confinement. Our assessment finds that the equilibrium is highly unstable to interchange-like pressure-driven instabilities near the plasma edge. The initially unstable modes rapidly destabilize other modes in the direction of the N-fold rotational symmetry (toroidal, in this case). For this equilibrium, N = 2, meaning destabilization of a large number of even-numbered toroidal Fourier modes. Thus, field-periodicity is likely to be an important factor in the nonlinear MHD stability characteristics of optimized stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Progress in extending high poloidal beta scenarios on DIII-D towards a steady-state fusion reactor and impact of energetic particles

To prepare for steady-state operation of future fusion reactors (e.g. the International Thermonuclear Experimental Reactor and China Fusion Engineering Test Reactor (CFETR)), experiments on DIII-D have extended the high poloidal beta ( β P ) scenario to reactor-relevant edge safety factor q 95 ~ 6.0, while maintaining a large-radius internal transport barrier (ITB) using negative magnetic shear. Excellent energy confinement quality ( H 98y2 > 1.5) is sustained at high normalized beta ( β N ~ 3.5). This high-performance ITB state with Greenwald density fraction near 100% and q min ≥ 3 is achieved with toroidal plasma rotation V tor ~ 0 at ρ ≥ 0.6. This is a key result for reactors expected to have low V tor . At high β P (≥1.9), large Shafranov shift can stabilize turbulence leading to a high confinement state with a low pedestal and an ITB. At lower β P (<1.9), negative magnetic shear in the plasma core contributes to turbulence suppression and can compensate for reduced Shafranov shift to continue to access a large-radius ITB and excellent confinement with low V tor , consistent with the results of gyrofluid transport simulations. These high- β P cases are characterized by weak/no Alfvén eigenmodes (a.e.) and classical fast-ion transport. At high density, the fast-ion deceleration time decreases and Δ β fast is lower; these reduce a.e. drive. The reverse-shear Alfvén eigenmodes are weaker or stable because the negative magnetic shear region is located at higher radius, away from the peaked fast-ion profile. Resistive wall modes can be a limitation at simultaneous high β N , low internal inductance, and low rotation. Analysis suggests that additional off-axis external current drive could provide a more stable path at reduced q 95 . Based on a DIII-D high- β P plasma with large-radius ITB, two scenarios are proposed for CFETR Q = 5 steady-state operation with ~1 GW fusion power: a lower-$l_i$($l_i$ ~ 0.66) and a higher-$l_i$($l_i$ ~ 0.75) case. Using a Landau closure model, multiple energetic particle (EP) effects on the a.e. stability are analyzed modifying the growth rate of the a.e.s triggered by the neutral-beam-injection EPs and alpha particles, although the stabilizing/destabilizing effect is weak for the cases analyzed. The stabilizing effects of the combined EP species β , energy, and density profile in CFETR need further investigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Interplay between beam-driven chirping modes and plasma confinement transitions in spherical tokamak ST40

Experiments on the high field spherical tokamak ST40 have led to the recent observation of interplay between beam-driven modes of sweeping frequency (chirping modes) and transitions to the enhanced global confinement regime (H-mode) and back to the low confinement regime (L-mode). The H-modes of plasma confinement are identified from decreased intensity of D α signal and from clear distinctions in the edge gradients of the visible plasma boundary (observed as a sharp plasma edge in camera images). The beam-driven chirping modes, identified as ideal magnetohydrodynamics beta-induced Alfvén acoustic eigenmodes modes, are observed in Mirnov coil signals, interferometry, and soft x-ray diagnostics. A moderate amplitude ‘primer’ chirping mode usually precedes an H–L transition. This is followed by a ‘dominant’ chirping mode with higher amplitude during the L-mode. The L–H transition back to the improved confinement occurs on a longer time scale of tens of ms, consistent with the slowing down time scale of fast beam ions. A dramatic decrease in toroidal plasma rotation is systematically observed associated with chirping modes sweeping down to zero frequency. Resonance maps built for the beam-driven chirping modes with the ASCOT (accelerated simulation of charged particle orbits in torodoial devices) code show that the resonant beam ions have orbits near the trapped-passing boundary. The ASCOT modelling assesses how losses of the resonant fast ions caused by the chirping modes with high enough amplitude modify the torque, potentially affecting the plasma rotation.

magnetic confinement fusion↗
DOE OSTI2022

Momentum transfer and plasma rotation caused by destabilized eigenmodes in tokamaks

In this work, the influence of magnetohydrodynamic eigenmodes destabilized by energetic ions on the momentum of these ions and concomitant sheared plasma rotation are studied. Two mechanisms affecting rotation are revealed: (i) spatial channelling (SC) – radially separated emission and absorption of the momentum; (ii) mode induced redistribution (MIR) across the magnetic field of the momentum of energetic ions. Forces arising during SC and MIR produce both toroidal and poloidal rotations. In addition, the momentum emission during SC leads to a radial flux of fast ions and generation of a radial electric field. Using the developed theory, estimates were made for the ITER (International Thermonuclear Experimental Reactor) 15 MA baseline scenario. They show that a global toroidicity-induced Alfvén eigenmode destabilized by alpha particles and neutral beam injection can result in significant radial electric field and forces applied to plasma. However, available data are not sufficient for a reliable prediction of the effects of SC and MIR in ITER. In general, one can expect that sheared rotation arising after destabilization of Alfvén modes and fast magnetoacoustic modes by energetic ions will tend to suppress the turbulence and improve plasma performance. The importance of plasma rotation is supported, in particular, by the fact that during the JET DTE1 experimental campaign the best parameters were achieved in a deuterium–tritium discharge where the rotation frequency was largest.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Weakly quasisymmetric near-axis solutions to all orders

Here, we show that the equations satisfied by weakly quasisymmetric magnetic fields can be solved to arbitrarily high order in powers of the distance from the magnetic axis. This demonstration does not consider force balance. The existence of solutions requires an appropriate choice of parameters, most notably the toroidal current or rotational transform profiles. We do not prove that the expansion converges (it is likely divergent but asymptotic), and thus, the demonstration here should not be taken as definitive proof of the existence of global solutions. Instead, we provide a systematic construction of solutions to an arbitrarily high order.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Runaway electron transport in stochastic toroidal magnetic fields

Here, we study the transport and confinement properties of runaway electrons (RE) in the presence of magnetic fields with perturbations producing different levels of stochasticity. We use Kinetic Orbit Runaway Electron Code (KORC) [Carbajal et al. , Phys. Plasmas 24, 042512 (2017) and del-Castillo-Negrete et al. , Phys. Plasmas 25, 056104 (2018)] for simulating the full-orbit (FO) and guiding-center (GC) dynamics of RE in perturbed magnetic fields that exhibit magnetic islands. We extend previous works on this problem [Wingen et al. , Nucl. Fusion 46, 941 (2006); Izzo et al. , Nucl. Fusion 51, 063032 (2011); Papp et al. , Nucl. Fusion 51, 043004 (2011); V. Izzo and P. Parks, Phys. Plasmas 24, 060705 (2017); and Sommariva et al. , Nucl. Fusion 58, 016043 (2018)] by studying in detail full-orbit effects on the RE dynamics. We quantify FO effects on RE transport by performing one-to-one comparisons between FO and GC simulations. It is found that, for the magnetic field configurations considered, GC simulations predict twice the RE losses of FO simulations for 1 MeV and four times the RE losses of FO simulations for 25 MeV. Similarly, we show how different GC and FO dynamics of RE moving around magnetic islands can be, especially in the scenario where the RE Larmor radius is on the order of the size of the magnetic island. We also study the role of rotation of the magnetic islands on RE confinement, and we find that low-frequency toroidal rotation has no observable effect on RE transport in the cases considered. These results shed some light into the potential of avoidance or mitigation mechanisms based on magnetic perturbations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗