Showing 1–29 of 29 results for author: Paz-Soldan, C

Power handling in a highly-radiative negative triangularity pilot plant

Authors: M. A. Miller, D. Arnold, M. Wigram, A. O. Nelson, J. Witham, G. Rutherford, H. Choudhury, C. Cummings, C. Paz-Soldan, D. G. Whyte

Abstract: This work explores power handling solutions for high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept \cite{rutherford_manta_2024}. The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only… ▽ More This work explores power handling solutions for high-field, highly-radiative negative triangularity (NT) reactors based around the MANTA concept \cite{rutherford_manta_2024}. The divertor design is kept as simple as possible, opting for a standard divertor with standard leg length. FreeGS is used to create an equilibrium for the boundary region, prioritizing a short outer leg length of only $\sim$50 cm ($\sim$40\% of the minor radius). The UEDGE code package is used for the boundary plasma solution, to track plasma temperatures and fluxes to the divertor targets. It is found that for $P_\mathrm{SOL}$ = 25 MW and $n_\mathrm{sep} = 0.96 \times 10^{20}$ m$^{-3}$, conditions consistent with initial core transport modeling, little additional power mitigation is necessary. For external impurity injection of just 0.13\% Ne, the peak heat flux density at the more heavily loaded outer targets falls to 7.8 MW/m$^{2}$, while the electron temperature $T_\mathrm{e}$ remains just under 5 eV. Scans around the parameter space reveal that even at densities lower than in the primary operating scenario, $P_\mathrm{SOL}$ can be increased up to 50 MW, so long as a slightly higher fraction of extrinsic radiator is used. With less than 1\% neon (Ne) impurity content, the divertor still experiences less than 10 MW/m$^{2}$ at the outer target. Design of the plasma-facing components includes a close-fitting vacuum vessel with a tungsten inner surface as well as FLiBe-carrying cooling channels fashioned into the VV wall directly behind the divertor targets. For the seeded heat flux profile, Ansys Fluent heat transfer simulations estimate that the outer target temperature remains at just below 1550\degree C. Initial scoping of advanced divertor designs shows that for an X-divertor, detachment of the outer target becomes much simpler, and plasma fluxes to the targets drop considerably with only 0.01\% Ne content. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Error Field Predictability and Consequences for ITER

Authors: Matthew Pharr, Nikolas Logan, Carlos Paz-Soldan, Jong-Kyu Park, Chistopher Hansen

Abstract: ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of $n=1$ error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils… ▽ More ITER coil tolerances are re-evaluated using the modern understanding of coupling to least-stable plasma modes and an updated center-line-traced model of ITER's coil windings. This reassessment finds the tolerances to be conservative through a statistical, linear study of $n=1$ error fields (EFs) due to tilted, shifted misplacements and nominal windings of central solenoid and poloidal field coils within tolerance. We also show that a model-based correction scheme remains effective even when metrology quality is sub-optimal, and compare this to projected empirical correction schemes. We begin with an analysis of the necessity of error field correction (EFC) for daily operation in ITER using scaling laws for the EF penetration threshold. We then consider the predictability of EF dominant mode overlap across early planned ITER scenarios and, as measuring EFs in high power scenarios can pose risks to the device, the potential for extrapolation to the ITER Baseline Scenario (IBS). We find that carefully designing a scenario matching currents proportionally to those of the IBS is far more important than plasma shape or profiles in accurately measuring an optimal correction current set. △ Less

Submitted 9 June, 2024; v1 submitted 3 June, 2024; originally announced June 2024.

Comments: 11 pages, 13 Figures. Submitted to Nuclear Fusion

MANTA: A Negative-Triangularity NASEM-Compliant Fusion Pilot Plant

Authors: MANTA Collaboration, G. Rutherford, H. S. Wilson, A. Saltzman, D. Arnold, J. L. Ball, S. Benjamin, R. Bielajew, N. de Boucaud, M. Calvo-Carrera, R. Chandra, H. Choudhury, C. Cummings, L. Corsaro, N. DaSilva, R. Diab, A. R. Devitre, S. Ferry, S. J. Frank, C. J. Hansen, J. Jerkins, J. D. Johnson, P. Lunia, J. van de Lindt, S. Mackie , et al. (16 additional authors not shown)

Abstract: The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicin… ▽ More The MANTA (Modular Adjustable Negative Triangularity ARC-class) design study investigated how negative-triangularity (NT) may be leveraged in a compact, fusion pilot plant (FPP) to take a ``power-handling first" approach. The result is a pulsed, radiative, ELM-free tokamak that satisfies and exceeds the FPP requirements described in the 2021 National Academies of Sciences, Engineering, and Medicine report ``Bringing Fusion to the U.S. Grid". A self-consistent integrated modeling workflow predicts a fusion power of 450 MW and a plasma gain of 11.5 with only 23.5 MW of power to the scrape-off layer (SOL). This low $P_\text{SOL}$ together with impurity seeding and high density at the separatrix results in a peak heat flux of just 2.8 MW/m$^{2}$. MANTA's high aspect ratio provides space for a large central solenoid (CS), resulting in ${\sim}$15 minute inductive pulses. In spite of the high B fields on the CS and the other REBCO-based magnets, the electromagnetic stresses remain below structural and critical current density limits. Iterative optimization of neutron shielding and tritium breeding blanket yield tritium self-sufficiency with a breeding ratio of 1.15, a blanket power multiplication factor of 1.11, toroidal field coil lifetimes of $3100 \pm 400$ MW-yr, and poloidal field coil lifetimes of at least $890 \pm 40$ MW-yr. Following balance of plant modeling, MANTA is projected to generate 90 MW of net electricity at an electricity gain factor of ${\sim}2.4$. Systems-level economic analysis estimates an overnight cost of US\$3.4 billion, meeting the NASEM FPP requirement that this first-of-a-kind be less than US\$5 billion. The toroidal field coil cost and replacement time are the most critical upfront and lifetime cost drivers, respectively. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Characterization of the ELM-free Negative Triangularity Edge on DIII-D

Authors: A. O. Nelson, L. Schmitz, T. Cote, J. F. Parisi, S. Stewart, C. Paz-Soldan, K. E. Thome, M. E. Austin, F. Scotti, J. L. Barr, A. Hyatt, N. Leuthold, A. Marinoni, T. Neiser, T. Osborne, N. Richner, A. S. Welander, W. P. Wehner, R. Wilcox, T. M. Wilks, J. Yang

Abstract: Tokamak plasmas with strong negative triangularity (NT) shaping typically exhibit fundamentally different edge behavior than conventional L-mode or H-mode plasmas. Over the entire DIII-D database, plasmas with sufficiently negative triangularity are found to be inherently free of edge localized modes (ELMs), even at injected powers well above the predicted L-H power threshold. A critical triangula… ▽ More Tokamak plasmas with strong negative triangularity (NT) shaping typically exhibit fundamentally different edge behavior than conventional L-mode or H-mode plasmas. Over the entire DIII-D database, plasmas with sufficiently negative triangularity are found to be inherently free of edge localized modes (ELMs), even at injected powers well above the predicted L-H power threshold. A critical triangularly ($δ_\mathrm{crit}\simeq-0.15$), consistent with inherently ELM-free operation is identified, beyond which access to the second stability region for infinite-$n$ ballooning modes closes on DIII-D. It is also possible to close access to this region, and thereby prevent an H-mode transition, at weaker average triangularities ($δ\lesssimδ_\mathrm{crit}$) provided that at least one of the two x-points is still sufficiently negative. Enhanced low field side magnetic fluctuations during ELM-free operation are consistent with additional turbulence limiting the NT edge gradient. Despite the reduced upper limit on the pressure gradient imposed by ballooning stability, NT plasmas are able to support small pedestals and are typically characterized by an enhancement of edge pressure gradients beyond those found in traditional L-mode plasmas. Further, the pressure gradient inside of this small pedestal is unusually steep, allowing access to high core performance that is competitive with other ELM-free regimes previously achieved on DIII-D. Since ELM-free operation in NT is linked directly to the magnetic geometry, NT fusion pilot plants are predicted to maintain advantageous edge conditions even in burning plasma regimes, potentially eliminating reactor core-integration issues caused by ELMs. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Highest Fusion Performance without Harmful Edge Energy Bursts in Tokamak

Authors: SangKyeun Kim, Ricardo Shousha, SeongMoo Yang, Qiming Hu, SangHee Hahn, Azarakhsh Jalalvand, Jong-Kyu Park, Nikolas Christopher Logan, Andrew Oakleigh Nelson, Yong-Su Na, Raffi Nazikian, Robert Wilcox, Rongjie Hong, Terry Rhodes, Carlos Paz-Soldan, YoungMu Jeon, MinWoo Kim, WongHa Ko, JongHa Lee, Alexander Battey, Alessandro Bortolon, Joseph Snipes, Egemen Kolemen

Abstract: The path of tokamak fusion and 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 high-confinement plasmas. The application of 3D magnetic perturbations is the method in ITER and possibly in future fusion power plants to suppress this instability and avoid energy bus… ▽ More The path of tokamak fusion and 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 high-confinement plasmas. The application of 3D magnetic perturbations is the method in ITER and possibly in future fusion power plants to suppress this instability and avoid energy busts damaging the device. Unfortunately, the conventional use of the 3D field in tokamaks typically leads to degraded fusion performance and an increased risk of other plasma instabilities, two severe issues for reactor implementation. In this work, we present an innovative 3D field optimization, exploiting machine learning, real-time adaptability, and multi-device capabilities to overcome these limitations. This integrated scheme is successfully deployed on DIII-D and KSTAR tokamaks, consistently achieving reactor-relevant core confinement and the highest fusion performance without triggering damaging instabilities or bursts while demonstrating ITER-relevant automated 3D optimization for the first time. This is enabled both by advances in the physics understanding of self-organized transport in the plasma edge and by advances in machine-learning technology, which is used to optimize the 3D field spectrum for automated management of a volatile and complex system. These findings establish real-time adaptive 3D field optimization as a crucial tool for ITER and future reactors to maximize fusion performance while simultaneously minimizing damage to machine components. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Operational Space and Plasma Performance with an RMP-ELM Suppressed Edge

Authors: C. Paz-Soldan, S. Gu, N. Leuthold, P. Lunia, P. Xie, M. W. Kim, S. K. Kim, N. C. Logan, J. -K. Park, W. Suttrop, Y. Sun, D. B. Weisberg, M. Willensdorfer, the ASDEX-Upgrade, DIII-D, EAST, KSTAR Teams

Abstract: The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are foun… ▽ More The operational space and global performance of plasmas with edge-localized modes (ELMs) suppressed by resonant magnetic perturbations (RMPs) are surveyed by comparing AUG, DIII-D, EAST, and KSTAR stationary operating points. RMP-ELM suppression is achieved over a range of plasma currents, toroidal fields, and RMP toroidal mode numbers. Consistent operational windows in edge safety factor are found across devices, while windows in plasma shaping parameters are distinct. Accessed pedestal parameters reveal a quantitatively similar pedestal-top density limit for RMP-ELM suppression in all devices of just over 3x1019 m-3. This is surprising given the wide variance of many engineering parameters and edge collisionalities, and poses a challenge to extrapolation of the regime. Wide ranges in input power, confinement time, and stored energy are observed, with the achieved triple product found to scale like the product of current, field, and radius. Observed energy confinement scaling with engineering parameters for RMP-ELM suppressed plasmas are presented and compared with expectations from established H and L-mode scalings, including treatment of uncertainty analysis. Different scaling exponents for individual engineering parameters are found as compared to the established scalings. However, extrapolation to next-step tokamaks ITER and SPARC find overall consistency within uncertainties with the established scalings, finding no obvious performance penalty when extrapolating from the assembled multi-device RMP-ELM suppressed database. Overall this work identifies common physics for RMP-ELM suppression and highlights the need to pursue this no-ELM regime at higher magnetic field and different plasma physical size. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 22 pages, 11 figures

Assessment of vertical stability for negative triangularity pilot plants

Authors: S. Guizzo, A. O. Nelson, C. Hansen, F. Logak, C. Paz-Soldan

Abstract: Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When couple… ▽ More Negative triangularity (NT) tokamak configurations may be more susceptible to magneto-hydrodynamic instability, posing challenges for recent reactor designs centered around their favorable properties, such as improved confinement and operation free of edge-localized modes. In this work, we assess the vertical stability of plasmas with NT shaping and develop potential reactor solutions. When coupled with a conformal wall, NT equilibria are confirmed to be less vertically stable than equivalent positive triangularity (PT) configurations. Unlike PT, their vertical stability is degraded at higher poloidal beta. Furthermore, improvements in vertical stability at low aspect ratio do not translate to the NT geometry. NT equilibria are stabilized in PT vacuum vessels due to the increased proximity of the plasma and the wall on the outboard side, but this scenario is found to be undesirable due to reduced vertical gaps which give less spatial margin for control recovery. Instead, we demonstrate that informed positioning of passively conducting plates can lead to improved vertical stability in NT configurations on par with stability metrics expected in PT scenarios. An optimal setup for passive plates in highly elongated NT devices is presented, where plates on the outboard side of the device reduce vertical instability growth rates to 16% of their baseline value. For lower target elongations, integration of passive stabilizers with divertor concepts can lead to significant improvements in vertical stability. Plates on the inboard side of the device are also uniquely enabled in NT geometries, providing opportunity for spatial separation of vertical stability coils and passive stabilizers. △ Less

Submitted 26 January, 2024; originally announced January 2024.

Implications of Vertical Stability Control on the SPARC Tokamak

Authors: A. O. Nelson, D. T. Garnier, D. J. Battaglia, C. Paz-Soldan, I. Stewart, M. Reinke, A. J. Creely, J. Wai

Abstract: To achieve its performance goals, SPARC plans to operate in equilibrium configurations with a strong elongation of $κ_\mathrm{areal}\sim1.75$, destabilizing the $n=0$ vertical instability. However, SPARC also features a relatively thick conducting wall that is designed to withstand disruption forces, leading to lower vertical instability growth rates than usually encountered. In this work, we use… ▽ More To achieve its performance goals, SPARC plans to operate in equilibrium configurations with a strong elongation of $κ_\mathrm{areal}\sim1.75$, destabilizing the $n=0$ vertical instability. However, SPARC also features a relatively thick conducting wall that is designed to withstand disruption forces, leading to lower vertical instability growth rates than usually encountered. In this work, we use the TokSyS framework to survey families of accessible shapes near the SPARC baseline configuration, finding maximum growth rates in the range of $γ\lesssim100\,$s$^{-1}$. The addition of steel vertical stability plates has only a modest ($\sim25\%$) effect on reducing the vertical growth rate and almost no effect on the plasma controllability when the full vertical stability system is taken into account, providing flexibility in the plate conductivity in the SPARC design. Analysis of the maximum controllable displacement on SPARC is used to inform the power supply voltage and current limit requirements needed to control an initial vertical displacement of $5\%$ of the minor radius. From the expected spectra of plasma disturbances and diagnostic noise, requirements for filter latency and vertical stability coil heating tolerances are also obtained. Small modifications to the outboard limiter location are suggested to allow for an unmitigated vertical disturbance as large as $5\%$ of the minor radius without allowing the plasma to become limited. Further, investigations with the 3D COMSOL code reveal that strategic inclusion of insulating structures within the VSC supports are needed to maintain sufficient magnetic response. The workflows presented here help to establish a model for the integrated predictive design for future devices by coupling engineering decisions with physics needs. △ Less

Submitted 17 January, 2024; originally announced January 2024.

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

Authors: C. Hansen, I. G. Stewart, D. Burgess, M. Pharr, S. Guizzo, F. Logak, A. O. Nelson, C. Paz-Soldan

Abstract: 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 description… ▽ More 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 orders for analytic and cross-code test cases. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Journal ref: Computer Physics Communications 298, 109111 (2024)

Design of Passive and Structural Conductors for Tokamaks Using Thin-Wall Eddy Current Modeling

Authors: A. F. Battey, C. Hansen, D. Garnier, D. Weisberg, C. Paz-Soldan, R. Sweeney, R. A. Tinguely, A. J. Creely

Abstract: A new three-dimensional electromagnetic modeling tool ThinCurr has been developed using the existing PSI-Tet finite-element code in support of conducting structure design work for both the SPARC and DIII-D tokamaks. Within this framework a 3D conducting structure model was created for both the SPARC and DIII-D tokamaks in the thin-wall limit. This model includes accurate details of the vacuum vess… ▽ More A new three-dimensional electromagnetic modeling tool ThinCurr has been developed using the existing PSI-Tet finite-element code in support of conducting structure design work for both the SPARC and DIII-D tokamaks. Within this framework a 3D conducting structure model was created for both the SPARC and DIII-D tokamaks in the thin-wall limit. This model includes accurate details of the vacuum vessel and other conducting structural elements with realistic material resistivities. This model was leveraged to support the design of a passive runaway electron mitigation coil (REMC), studying the effect of various design parameters, including coil resistivity, current quench duration, and plasma vertical position, on the effectiveness of the coil. The REMC is a non-axisymmetric coil designed to passively drive large non-axisymmetric fields during the plasma disruption thereby destroying flux surfaces and deconfining RE seed populations. These studies indicate that current designs should apply substantial 3D fields at the plasma surface during future plasma current disruptions as well as highlight the importance of having the REMC conductors away from the machine midplane in order to ensure they are robust to off-normal disruption scenarios. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 33 pages and 14 figures. to be published in Nuclear Fusion

Simultaneous access to high normalized current, pressure, density, and confinement in strongly-shaped diverted negative triangularity plasmas

Authors: C. Paz-Soldan, C. Chrystal, P. Lunia, A. O. Nelson, K. E. Thome, M. E. Austin, T. B. Cote, A. W. Hyatt, A. Marinoni, T. H. Osborne, M. Pharr, O. Sauter, F. Scotti, T. M. Wilks, H. S. Wilson

Abstract: Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized current, pressure, density, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangu… ▽ More Strongly-shaped diverted negative triangularity (NT) plasmas in the DIII-D tokamak demonstrate simultaneous access to high normalized current, pressure, density, and confinement. NT plasmas are shown to exist across an expansive parameter space compatible with high fusion power production, revealing surprisingly good core stability properties that compare favorably to conventional positive triangularity plasmas in DIII-D. Non-dimensionalizing the operating space, edge safety factors below 3, normalized betas above 3, Greenwald density fractions above 1, and high-confinement mode (H-mode) confinement qualities above 1 are simultaneously observed, all with a robustly stable edge free from deleterious edge-localized mode instabilities. Scaling of the confinement time with engineering parameters reveals at least a linear dependence on plasma current although with significant power degradation, both in excess of expected H-mode scalings. These results increase confidence that NT plasmas are a viable approach to realize fusion power and open directions for future detailed study. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Robust avoidance of edge-localized modes alongside gradient formation in the negative triangularity tokamak edge

Authors: A. O. Nelson, L. Schmitz, C. Paz-Soldan, K. E. Thome, T. B. Cote, N. Leuthold, F. Scotti, M. E. Austin, A. Hyatt, T. Osborne

Abstract: In a series of high performance diverted discharges on DIII-D, we demonstrate that strong negative triangularity (NT) shaping robustly suppresses all edge-localized mode (ELM) activity over a wide range of plasma conditions: $\langle n\rangle=0.1-1.5\times10^{20}$m$^{-3}$, $P_\mathrm{aux}=0-15$MW and $|B_\mathrm{t}|=1-2.2$T, corresponding to $P_\mathrm{loss}/P_\mathrm{LH08}\sim8$. The full dataset… ▽ More In a series of high performance diverted discharges on DIII-D, we demonstrate that strong negative triangularity (NT) shaping robustly suppresses all edge-localized mode (ELM) activity over a wide range of plasma conditions: $\langle n\rangle=0.1-1.5\times10^{20}$m$^{-3}$, $P_\mathrm{aux}=0-15$MW and $|B_\mathrm{t}|=1-2.2$T, corresponding to $P_\mathrm{loss}/P_\mathrm{LH08}\sim8$. The full dataset is consistent with the theoretical prediction that magnetic shear in the NT edge inhibits access to ELMing H-mode regimes; all experimental pressure profiles are found to be at or below the infinite-$n$ ballooning stability limit. Importantly, we also report enhanced edge pressure gradients at strong NT that are significantly steeper than in traditional ELM-free L-mode plasmas and provide significant promise for NT reactor integration. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 5 pages, 5 figures

Self-consistent simulation of compressional Alfvén eigenmodes excited by runaway electrons

Authors: Chang Liu, Andrey Lvovskiy, Carlos Paz-Soldan, Stephen C. Jardin, Amitava Bhattacharjee

Abstract: Alfvénic modes in the current quench (CQ) stage of the tokamak disruption have been observed in experiments. In DIII-D the excitation of these modes is associated with the presence of high-energy runaway electrons, and a strong mode excitation is often associated with the failure of RE plateau formation. In this work we present results of self-consistent kinetic-MHD simulations of RE-driven compre… ▽ More Alfvénic modes in the current quench (CQ) stage of the tokamak disruption have been observed in experiments. In DIII-D the excitation of these modes is associated with the presence of high-energy runaway electrons, and a strong mode excitation is often associated with the failure of RE plateau formation. In this work we present results of self-consistent kinetic-MHD simulations of RE-driven compressional Alfvén eigenmodes (CAEs) in DIII-D disruption scenarios, providing an explanation of the CQ modes. Simulation results reveal that high energy trapped REs can have resonance with the Alfvén mode through their precession motion, and the resonance frequency is proportional to the energy of REs. The mode frequencies and their relationship with the RE energy are consistent with experimental observation. The perturbed magnetic fields from the modes can lead to spatial diffusion of runaway electrons including the nonresonant passing ones, thus providing the theoretical basis for a potential approach for runaway electron mitigation. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: 5 pages, 4 figures

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

Authors: R. A. Tinguely, I. Pusztai, V. A. Izzo, K. S'"arkimäki, T. Fülöp, D. T. Garnier, R. S. Granetz, M. Hoppe, C. Paz-Soldan, A. Sundström, R. Sweeney

Abstract: In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2… ▽ More In [V.A. Izzo et al 2022 Nucl. Fusion 62 096029], state-of-the-art modeling of thermal and current quench (CQ) MHD coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1 MA compared to ~8.7 MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.3 and with kinetic energies 0.2-15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 $\mathrm{m^2/s}$, is needed in the core to fully mitigate these REs. Properly accounting for (i) the CQ electric field's effect on RE transport in islands and (ii) the contribution of significant RE currents to disruption MHD may help achieve this. △ Less

Submitted 3 January, 2023; originally announced January 2023.

Runaway electron deconfinement in SPARC and DIII-D by a passive 3D coil

Authors: V. A. Izzo, I. Pusztai, K. Särkimäki, A. Sundström, D. Garnier, D. Weisberg, R. A. Tinguely, C. Paz-Soldan, R. S. Granetz, R. Sweeney

Abstract: The operation of a 3D coil--passively driven by the current quench loop voltage--for the deconfinement of runaway electrons is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear MHD modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code… ▽ More The operation of a 3D coil--passively driven by the current quench loop voltage--for the deconfinement of runaway electrons is modeled for disruption scenarios in the SPARC and DIII-D tokamaks. Nonlinear MHD modeling is carried out with the NIMROD code including time-dependent magnetic field boundary conditions to simulate the effect of the coil. Further modeling in some cases uses the ASCOT5 code to calculate advection and diffusion coefficients for runaway electrons based on the NIMROD-calculated fields, and the DREAM code to compute the runaway evolution in the presence of these transport coefficients. Compared with similar modeling in Tinguely, et al [2021 Nucl. Fusion 61 124003], considerably more conservative assumptions are made with the ASCOT5 results, zeroing low levels of transport, particularly in regions in which closed flux surfaces have reformed. Of three coil geometries considered in SPARC, only the $n=1$ coil is found to have sufficient resonant components to suppress the runaway current growth. Without the new conservative transport assumptions, full suppression of the RE current is maintained when the TQ MHD is included in the simulation or when the RE current is limited to 250kA, but when transport in closed flux regions is fully suppressed, these scenarios allow RE beams on the order of 1-2MA to appear. Additional modeling is performed to consider the effects of the close ideal wall. In DIII-D, the current quench is modeled for both limited and diverted equilibrium shapes. In the limited shape, the onset of stochasticity is found to be insensitive to the coil current amplitude and governed largely by the evolution of the safety-factor profile. In both devices, prediction of the q-profile evolution is seen to be critical to predicting the later time effects of the coil. △ Less

Submitted 25 July, 2022; originally announced July 2022.

H-mode inhibition in negative triangularity tokamak reactor plasmas

Authors: A. O. Nelson, C. Paz-Soldan, S. Saarelma

Abstract: Instability to high toroidal mode number ($n$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity ($δ$) shaping, preventing access to strong H-mode regimes when $δ\ll0$. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$n$ ballooning stability as a function of internal profiles a… ▽ More Instability to high toroidal mode number ($n$) ballooning modes has been proposed as the primary gradient-limiting mechanism for tokamak equilibria with negative triangularity ($δ$) shaping, preventing access to strong H-mode regimes when $δ\ll0$. To understand how this mechanism extrapolates to reactor conditions, we model the infinite-$n$ ballooning stability as a function of internal profiles and equilibrium shape using a combination of the CHEASE and BALOO codes. While the critical $δ$ required for avoiding $2^\mathrm{nd}$ stability to high-$n$ modes is observed to depend in a complicated way on various shaping parameters, including the equilibrium aspect ratio, elongation and squareness, equilibria with negative triangularity are robustly prohibited from accessing the $2^\mathrm{nd}$ stability region, offering the prediction that that negative triangularity reactors should maintain L-mode-like operation. In order to access high-$n$ $2^\mathrm{nd}$ stability, the local shear over the entire bad curvature region must be sufficiently negative to overcome curvature destabilization on the low field side. Scalings of the ballooning-limited pedestal height are provided as a function of plasma parameters to aid future scenario design. △ Less

Submitted 26 April, 2022; originally announced April 2022.

Latest results on quiescent and post-disruption runaway electron mitigation experiments at Frascati Tokamak Upgrade

Authors: D. Carnevale, P. Buratti, M. Baruzzo, W. Bin, F. Bombarda, L. Boncagni, C. Paz-Soldan, L. Calacci, M. Cappelli, C. Castaldo, S. Ceccuzzi, C. Centioli, C. Cianfarani, S. Coda, F. Cordella, O. D Arcangelo, J. Decker, B. Duval, B. Esposito, L. Gabellieri, S. Galeani, S. Garavaglia, C. Galperti, G. Ghillardi, G. Granucci , et al. (16 additional authors not shown)

Abstract: Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissi… ▽ More Results from the last FTU campaigns on the deuterium large (wrt FTU volume) pellet REs suppression capability, mainly due to the induced burst MHD activity expelling REs seed are presented for discharges with 0.5 MA and 5.3T. Clear indications of avalanche multiplication of REs following single pellet injection on 0.36 MA flat-top discharges is shown together with quantitative indications of dissipative effects in terms of critical electrical field increase due to fan-like instabilities. Analysis of large fan-like instabilities on post-disruption RE beams, that seem to be correlated with low electrical field and background density drops, reveal their strong RE energy suppression capability suggesting a new strategy for RE energy suppression controlling large fan instabilities. We demonstrate how such density drops can be induced using modulated ECRH power on post-disruption beams. △ Less

Submitted 25 May, 2021; v1 submitted 10 May, 2021; originally announced May 2021.

Comments: DRAFT that will be submitted to NF (the conference version is on FEC/IAEA conference 2020)

Self-consistent simulation of resistive kink instabilities with runaway electrons

Authors: Chang Liu, Chen Zhao, Stephen C. Jardin, Nathaniel Ferraro, Carlos Paz-Soldan, Yueqiang Liu, Brendan C. Lyons

Abstract: A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris… ▽ More A new fluid model for runaway electron simulation based on fluid description is introduced and implemented in the magnetohydrodynamics code M3D-C1, which includes self-consistent interactions between plasma and runaway electrons. The model utilizes the method of characteristics to solve the continuity equation for the runaway electron density with large convection speed, and uses a modified Boris algorithm for pseudo particle pushing. The model was employed to simulate magnetohydrodynamics instabilities happening in a runaway electron final loss event in the DIII-D tokamak. Nonlinear simulation reveals that a large fraction of runaway electrons get lost to the wall when kink instabilities are excited and form stochastic field lines in the outer region of the plasma. Plasma current converts from runaway electron current to Ohmic current, and get pinched at the magnetic axis. Given the good agreement with experiment, the simulation model provides a reliable tool to study macroscopic plasma instabilities in existence of runaway electron current, and can be used to support future studies of runaway electron mitigation strategies in ITER. △ Less

Submitted 21 April, 2021; originally announced April 2021.

Comments: 26 pages, 16 figures

Plasma Performance and Operational Space without ELMs in DIII-D

Authors: C. Paz-Soldan

Abstract: A database of DIII-D plasmas without ELMs compares the operating space and plasma performance of stationary no-ELM regimes: RMP-ELM suppression, QH-mode, I-mode, EDA H-mode, L-mode, and negative triangularity L-mode (Neg-D). Operational space is documented in terms of engineering and physics parameters, revealing divergent constraints. Some operational space discriminants (such as pedestal collisi… ▽ More A database of DIII-D plasmas without ELMs compares the operating space and plasma performance of stationary no-ELM regimes: RMP-ELM suppression, QH-mode, I-mode, EDA H-mode, L-mode, and negative triangularity L-mode (Neg-D). Operational space is documented in terms of engineering and physics parameters, revealing divergent constraints. Some operational space discriminants (such as pedestal collisionality) are well known, while others, such as low torque & safety factor, or high power & density, are less commonly emphasized. Normalized performance (confinement quality and normalized pressure) also discriminate the no-ELM regimes and favor the regimes tolerant to power in DIII-D: RMP, QH, and Neg-D. Absolute performance (triple product <p>$τ$) also discriminates no-ELM regimes and is found to rise linearly with IaB and also benefits from tolerance to power. The highest normalized performance using the metric <p>$τ$/IaB is found in QH and RMP regimes. Though comparable QH and RMP performance is found, the pedestal pressure (p_ped) is very different. p_ped in RMP plasmas is relatively low, and the best performance is found with a high core <p> fraction alongside high core rotation, consistent with an ExB shear confinement enhancement. p_ped of QH plasmas is significantly higher than RMP, and QH performance does not correlate with core rotation. However, the highest QH p_ped are found with high carbon fraction. While normalized performance of Neg-D plasmas is comparable to QH and RMP plasmas, the absolute confinement of Neg-D is lower, owing to low elongation and low p_ped achieved thus far. Only the EDA-H, Neg-D, and L-mode scenarios have approached divertor-compatible high separatrix density conditions, with Neg-D preserving the highest performance owing to its compatibility with both high power and density. △ Less

Submitted 6 December, 2020; originally announced December 2020.

Comments: 43 figures

Compressional Alfvén eigenmodes excited by runaway electrons

Authors: Chang Liu, Dylan P. Brennan, Andrey Lvovskiy, Carlos Paz-Soldan, Eric D. Fredrickson, Amitava Bhattacharjee

Abstract: Compressional Alfvén eigenmodes (CAE) driven by energetic ions have been observed in magnetic fusion experiments. In this paper, we show that the modes can also be driven by runaway electrons formed in post-disruption plasma, which may explain kinetic instabilities observed in DIII-D disruption experiments with massive gas injection. The mode-structure is calculated, as are the frequencies which a… ▽ More Compressional Alfvén eigenmodes (CAE) driven by energetic ions have been observed in magnetic fusion experiments. In this paper, we show that the modes can also be driven by runaway electrons formed in post-disruption plasma, which may explain kinetic instabilities observed in DIII-D disruption experiments with massive gas injection. The mode-structure is calculated, as are the frequencies which are in agreement with experimental observations. Using a runaway electron distribution function obtained from a kinetic simulation, the mode growth rates are calculated and found to exceed the collisional damping rate when the runaway electron density exceeds a threshold value. The excitation of CAE poses a new possible approach to mitigate seed runaway electrons during the current quench and surpassing the avalanche. △ Less

Submitted 1 October, 2020; v1 submitted 24 September, 2020; originally announced September 2020.

Comments: 20 pages, 8 figures

Non-Planar Coil Winding Angle Optimization for Compatibility with Non-Insulated High-Temperature Superconducting Magnets

Authors: Carlos Paz-Soldan

Abstract: The rapidly emerging technology of high-temperature superconductors (HTS) opens new opportunities for the development of non-planar non-insulated HTS magnets. This type of HTS magnet offers attractive features via its simplicity, robustness, and is well-suited for modest size steady-state applications such as a mid-scale stellarator. In non-planar coil applications the HTS tape may be subject to s… ▽ More The rapidly emerging technology of high-temperature superconductors (HTS) opens new opportunities for the development of non-planar non-insulated HTS magnets. This type of HTS magnet offers attractive features via its simplicity, robustness, and is well-suited for modest size steady-state applications such as a mid-scale stellarator. In non-planar coil applications the HTS tape may be subject to severe minor-axis bending strain ($ε_{bend}$), torsional strains ($ε_{tor}$) and transverse magnetic field components ($B_\perp$), all of which can limit the magnet operating space. A novel method of winding angle optimization is here presented to overcome these limitations. Essentially, this method: 1) calculates the peak εbend and B\perp for arbitary winding angle along an input coil filamentary trajectory, 2) defines a cost function including both, and then 3) uses tensioned splines to define a winding angle that reduces $ε$tor and optimizes the εbend and $B_\perp$ cost function. As strain limits are present even without $B_\perp$, this optimization is able to provide an assessment of the minimimum buildable size of an arbitrary non-planar non-insulating HTS coil. This optimization finds that for standard 4 mm wide HTS tapes the minimum size coils of the existing HSX, NCSX, and W7-X stellarator geometries are around 0.3 - 0.5 m in radius. For coils larger than this size, permitting a finite (yet tolerable) strain allows reduction of $B_\perp$. This enables a reduction of the HTS tape length required to achieve a given design magnetic field or equivalently an increase in the achievable magnetic field for fixed HTS tape length. The distinct considerations for optimizing a stellarator coilset to further ease compatibility with non-insulated HTS magnets are also discussed. △ Less

Submitted 3 March, 2020; originally announced March 2020.

Runaway Electron Seed Formation at Reactor-Relevant Temperature

Authors: C. Paz-Soldan, P. Aleynikov, E. M. Hollmann, A. Lvovskiy, I. Bykov, X. Du, N. W. Eidietis, D. Shiraki

Abstract: Systematic variation of the pre-disruption core electron temperature (Te) from 1 to 12 keV using an internal transport barrier scenario reveals a dramatic increase in the production of seed runaway electrons (REs), ultimately accessing near-complete conversion of the pre-disruption current into sub-MeV RE current. Injected Ar pellets are observed to ablate more intensely and promptly as Te rises.… ▽ More Systematic variation of the pre-disruption core electron temperature (Te) from 1 to 12 keV using an internal transport barrier scenario reveals a dramatic increase in the production of seed runaway electrons (REs), ultimately accessing near-complete conversion of the pre-disruption current into sub-MeV RE current. Injected Ar pellets are observed to ablate more intensely and promptly as Te rises. At high Te, the observed ablation exceeds predictions from published thermal ablation models. Simultaneously, the thermal quench (TQ) is observed to significantly shorten with increasing Te -- a surprising result. While the reason for the shorter TQ is not yet understood, candidate mechanisms include: insufficiently accurate thermal ablation models, enhanced ablation driven by the seed RE population, or significant parallel heat transport along stochastic fields. Kinetic modeling that self-consistently treats the plasma cooling via radiation, the induced electric field, and the formation of the seed RE is performed. Including the combined effect of the inherent dependence of hot-tail RE seeding on Te together with the shortened TQ, modeling recovers the progression towards near-complete conversion of the pre-disruption current to RE current as Te rises. Measurement of the HXR spectrum during the early current quench (CQ) reveals a trend of decreasing energy with pre-disruption Te. At the very highest Te ($\approx$ 12 keV), $\approx$ 100% conversion of the thermal current to runaway current is found. The energy of this peculiar RE beam is inferred to be sub-MeV as it emits vanishingly few MeV hard X-rays (HXRs). These measurements demonstrate novel TQ dynamics as Te is varied and illustrate the limitations of treating the RE seed formation problem without considering the inter-related dependencies of the pellet ablation, radiative energy loss, and resultant variations of the TQ duration. △ Less

Submitted 12 February, 2020; originally announced February 2020.

Wide operational windows of edge-localized mode suppression by resonant magnetic perturbations in the DIII-D tokamak

Authors: Q. M. Hu, R. Nazikian, B. A. Grierson, N. C. Logan, D. M. Orlov, C. Paz-Soldan, Q. Yu

Abstract: Edge-Localized-Mode (ELM) suppression by resonant magnetic perturbations (RMPs) generally occurs over very narrow ranges of the plasma current (or magnetic safety factor q95) in the DIII-D tokamak. However, wide q95 ranges of ELM suppression are needed for the safety and operational flexibility of ITER and future reactors. In DIII-D ITER Similar Shape (ISS) plasmas with n=3 RMPs, the range of q95… ▽ More Edge-Localized-Mode (ELM) suppression by resonant magnetic perturbations (RMPs) generally occurs over very narrow ranges of the plasma current (or magnetic safety factor q95) in the DIII-D tokamak. However, wide q95 ranges of ELM suppression are needed for the safety and operational flexibility of ITER and future reactors. In DIII-D ITER Similar Shape (ISS) plasmas with n=3 RMPs, the range of q95 for ELM suppression is found to increase with decreasing electron density. Nonlinear two-fluid MHD simulations reproduce the observed q95 windows of ELM suppression and the dependence on plasma density, based on the conditions for resonant field penetration at the top of the pedestal. When the RMP amplitude is close to the threshold for resonant field penetration, only narrow isolated magnetic islands form near the top of the pedestal, leading to narrow q95 windows of ELM suppression. However, as the threshold for field penetration decreases with decreasing density, resonant field penetration can take place over a wider range of q95. For sufficiently low density (penetration threshold) multiple magnetic islands form near the top of the pedestal giving rise to continuous q95 windows of ELM suppression. The model predicts that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure in DIII-D by shifting to higher toroidal mode number (n=4) RMPs. △ Less

Submitted 30 June, 2020; v1 submitted 13 December, 2019; originally announced December 2019.

Journal ref: Phys. Rev. Lett. 125, 045001 (2020)

The Density Dependence of Edge-Localized-Mode Suppression and Pump-out by Resonant Magnetic Perturbations in the DIII-D Tokamak

Authors: Q. M. Hu, R. Nazikian, B. Grierson, N. C. Logan, J-K. Park, C. Paz-Soldan, Q. Yu

Abstract: The density dependence of edge-localized-mode (ELM) suppression and density pump-out (density reduction) by n = 2 resonant magnetic perturbations (RMPs) is consistent with the effects of narrow well-separated magnetic islands at the top and bottom of the H-mode pedestal in DIII-D low-collisionality plasmas. Nonlinear two-fluid MHD simulations for DIII-D ITER Similar Shape (ISS) discharges show tha… ▽ More The density dependence of edge-localized-mode (ELM) suppression and density pump-out (density reduction) by n = 2 resonant magnetic perturbations (RMPs) is consistent with the effects of narrow well-separated magnetic islands at the top and bottom of the H-mode pedestal in DIII-D low-collisionality plasmas. Nonlinear two-fluid MHD simulations for DIII-D ITER Similar Shape (ISS) discharges show that, at low collisionality, low pedestal density is required for resonant field penetration at the pedestal top, consistent with the ubiquitous low density requirement for ELM suppression in these DIII-D plasmas. The simulations predict a drop in the pedestal pressure due to parallel transport across these narrow width (0.02) magnetic islands at the top of the pedestal that is stabilizing to Peeling-Ballooning-Modes (PBMs), and comparable to the pedestal pressure reduction observed in experiment at the onset of ELM suppression. The simulations predict density pump-out at experimentally relevant levels (-20%) at low pedestal collisionality due to very narrow (~0.01-0.02) RMP driven magnetic islands at the pedestal foot at ~0.99. The simulations show decreasing pump-out with increasing density, consistent with experiment, resulting from the inverse dependence of parallel particle transport on resistivity at the foot of the pedestal. The robust screening of resonant fields is predicted between the top and bottom of the pedestal during density pump-out and ELM suppression, consistent with the preservation of strong temperature gradients in the edge transport barrier as seen in experiment. △ Less

Submitted 13 December, 2019; v1 submitted 30 October, 2019; originally announced October 2019.

Journal ref: Physics of Plasmas, 2019

Experimental conditions to suppress edge localised modes by magnetic perturbations in the ASDEX Upgrade tokamak

Authors: W. Suttrop, A. Kirk, V. Bobkov, M. Cavedon, M. Dunne, R. M. McDermott, H. Meyer, R. Nazikian, C. Paz-Soldan, D. A. Ryan, E. Viezzer, M. Willensdorfer

Abstract: Access conditions for full suppression of Edge Localised Modes (ELMs) by Magnetic Perturbations (MP) in low density high confinement mode (H-mode) plasmas are studied in the ASDEX Upgrade tokamak. The main empirical requirements for full ELM suppression in our experiments are: 1. The poloidal spectrum of the MP must be aligned for best plasma response from weakly stable kink-modes, which amplify t… ▽ More Access conditions for full suppression of Edge Localised Modes (ELMs) by Magnetic Perturbations (MP) in low density high confinement mode (H-mode) plasmas are studied in the ASDEX Upgrade tokamak. The main empirical requirements for full ELM suppression in our experiments are: 1. The poloidal spectrum of the MP must be aligned for best plasma response from weakly stable kink-modes, which amplify the perturbation, 2. The plasma edge density must be below a critical value, $3.3 \times 10^{19}$~m$^{-3}$. The edge collisionality is in the range $ν^*_i = 0.15-0.42$ (ions) and $ν^*_e = 0.15-0.25$ (electrons). However, our data does not show that the edge collisionality is the critical parameter that governs access to ELM suppression. 3. The pedestal pressure must be kept sufficiently low to avoid destabilisation of small ELMs. This requirement implies a systematic reduction of pedestal pressure of typically 30\% compared to unmitigated ELMy H-mode in otherwise similar plasmas. 4. The edge safety factor $q_{95}$ lies within a certain window. Within the range probed so far, $q_{95}=3.5-4.2$, one such window, $q_{95}=3.57-3.95$ has been identified. Within the range of plasma rotation encountered so far, no apparent threshold of plasma rotation for ELM suppression is found. This includes cases with large cross field electron flow in the entire pedestal region, for which two-fluid MHD models predict that the resistive plasma response to the applied MP is shielded. △ Less

Submitted 29 June, 2018; v1 submitted 3 April, 2018; originally announced April 2018.

The effects of kinetic instabilities on the electron cyclotron emission from runaway electrons

Authors: Chang Liu, Lei Shi, Eero Hirvijoki, Dylan P. Brennan, Amitava Bhattacharjee, Carlos Paz-Soldan, Max E. Austin

Abstract: In this paper we show that the kinetic instabilities associated with runaway electron beams play an essential role for the production of high-level non-thermal electron-cyclotron-emission (ECE) radiation. Most of the non-thermal ECE comes from runaway electrons in the low-energy regime with large pitch angle, which are strongly scattered by the excited whistler waves. The power of ECE from runaway… ▽ More In this paper we show that the kinetic instabilities associated with runaway electron beams play an essential role for the production of high-level non-thermal electron-cyclotron-emission (ECE) radiation. Most of the non-thermal ECE comes from runaway electrons in the low-energy regime with large pitch angle, which are strongly scattered by the excited whistler waves. The power of ECE from runaway electrons is obtained using a synthetic diagnostic model based on the reciprocity method. The electron distribution function is calculated using a kinetic simulation model including the whistler wave instabilities and the quasilinear diffusion effects. Simulations based on DIII-D low-density discharge reproduces the rapid growth of the ECE signals observed in DIII-D experiments. Unlike the thermal ECE where radiation for a certain frequency is strongly localized inside the resonance region, the non-thermal ECE radiation from runaway electrons is nonlocal, and the emission-absorption ratio is higher than that of thermal electrons. The runaway electron tail is more significant for ECE with higher frequencies, and the ECE spectrum becomes flatter as RE population grows. The nonlinear behavior of the kinetic instabilities is illustrated in the osculations of the ECE waves. The good agreement with the DIII-D experimental observations after including the kinetic instabilities clearly illustrate the significance of the scattering effects from wave-particle interactions, which can also be important for runaway electrons produced in disruptions. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Role of kinetic instability in runaway electron avalanche and elevated critical electric fields

Authors: Chang Liu, Eero Hirvijoki, Guo-yong Fu, Dylan P. Brennan, Amitava Bhattacharjee, Carlos Paz-Soldan

Abstract: The effects of kinetic whistler-wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from DIII-D experiments, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results seen in a variety of tokamaks. We find an increase of the avalanche growth rate and threshold electric field, bringing the present model m… ▽ More The effects of kinetic whistler-wave instabilities on the runaway-electron (RE) avalanche is investigated. With parameters from DIII-D experiments, we show that RE scattering from excited whistler waves can explain several poorly understood experimental results seen in a variety of tokamaks. We find an increase of the avalanche growth rate and threshold electric field, bringing the present model much closer to observations than previous results. The excitation of kinetic instabilities and the scattering of resonant electrons are calculated self-consistently using a quasilinear model. We also explain the observed fast growth of electron cyclotron emission (ECE) signals and excitation of very low-frequency whistler modes observed in the quiescent RE experiments at DIII-D [D. A. Spong et al., submitted to Phys. Rev. Lett.]. These results indicate that by controlling the background thermal plasma temperature, the plasma wave can be excited spontaneously in tokamak disruptions and the avalanche generation of runaway electrons may be suppressed. △ Less

Submitted 5 January, 2018; originally announced January 2018.

Journal ref: Phys. Rev. Lett. 120, 265001 (2018)

Interpretation of runaway electron synchrotron and bremsstrahlung images

Authors: M. Hoppe, O. Embréus, C. Paz-Soldan, R. A. Moyer, T. Fülöp

Abstract: The crescent spot shape observed in DIII-D runaway electron synchrotron radiation images is shown to result from the high degree of anisotropy in the emitted radiation, the finite spectral range of the camera and the distribution of runaways. The finite spectral camera range is found to be particularly important, as the radiation from the high-field side can be stronger by a factor $10^6$ than the… ▽ More The crescent spot shape observed in DIII-D runaway electron synchrotron radiation images is shown to result from the high degree of anisotropy in the emitted radiation, the finite spectral range of the camera and the distribution of runaways. The finite spectral camera range is found to be particularly important, as the radiation from the high-field side can be stronger by a factor $10^6$ than the radiation from the low-field side in DIII-D. By combining a kinetic model of the runaway dynamics with a synthetic synchrotron diagnostic we see that physical processes not described by the kinetic model (such as radial transport) are likely to be limiting the energy of the runaways. We show that a population of runaways with lower dominant energies and larger pitch-angles than those predicted by the kinetic model provide a better match to the synchrotron measurements. Using a new synthetic bremsstrahlung diagnostic we also simulate the view of the Gamma Ray Imager (GRI) diagnostic used at DIII-D to resolve the spatial distribution of runaway-generated bremsstrahlung. △ Less

Submitted 12 December, 2017; originally announced December 2017.

Comments: 21 pages, 11 figures

Resistive and ferritic-wall plasma dynamos in a sphere

Authors: I. V. Khalzov, B. P. Brown, E. J. Kaplan, N. Katz, C. Paz-Soldan, K. Rahbarnia, E. J. Spence, C. B. Forest

Abstract: We numerically study the effects of varying electric conductivity and magnetic permeability of the bounding wall on a kinematic dynamo in a sphere for parameters relevant to Madison plasma dynamo experiment (MPDX). The dynamo is excited by a laminar, axisymmetric flow of von Karman type. The flow is obtained as a solution to the Navier-Stokes equation for an isothermal fluid with a velocity profil… ▽ More We numerically study the effects of varying electric conductivity and magnetic permeability of the bounding wall on a kinematic dynamo in a sphere for parameters relevant to Madison plasma dynamo experiment (MPDX). The dynamo is excited by a laminar, axisymmetric flow of von Karman type. The flow is obtained as a solution to the Navier-Stokes equation for an isothermal fluid with a velocity profile specified at the sphere's boundary. The properties of the wall are taken into account as thin-wall boundary conditions imposed on the magnetic field. It is found that an increase in the permeability of the wall reduces the critical magnetic Reynolds number Rm_cr. An increase in the conductivity of the wall leaves Rm_cr unaffected, but reduces the dynamo growth rate. △ Less

Submitted 21 May, 2012; originally announced May 2012.

Journal ref: Physics of Plasmas 19, 104501 (2012)