Engineering topics

Loarte, A.

Publications and source records attributed to Loarte, A..

15 records

DOE OSTI2024

Initial design concepts for solid boron injection in ITER

As part of ITER’s consideration to change its first wall material from beryllium to tungsten, the ITER Organization has proposed studying the feasibility of real-time solid boron injection (SBI) into the plasma to coat the walls and divertor to supplement glow discharge boronization (GDB). Boron deposits getter oxygen and reduce sputtering of tungsten from plasma facing components (PFCs). Particularly in areas with significant plasma wall interactions, boron coatings are expected to be short-lived under high performance plasma conditions. The proposed SBI system aims to maintain boron layers in these areas to avoid excessive radiation from tungsten in the plasma as a risk mitigation to ensure ITER will be able to reach and sustain Q = 10 conditions. The system will be used sparingly, as redeposition of boron can lead to significant tritium retention, which must be minimized in ITER to comply with nuclear safety concerns. SBI is proposed to limit and precisely control the amount of boron injected in real-time during plasma operation. Here, some of the design requirements and initial concepts for an SBI system in ITER are presented based on previous results carried out with SBI systems on a number of tokamaks and stellarators around the world. Previous results using SBI systems installed by PPPL have injected boron particles from 5 µm–2 mm diameter at calibrated rates of 2–200 mg/s in real-time during plasma operation on AUG, DIII-D, EAST, KSTAR, LHD, TFTR, WEST, and W7-X, leading to improved wall conditions with reduced plasma impurity concentrations and radiated power and improved plasma performance. The boron is ionized in the plasma edge and then deposited on plasma-wetted surfaces. On AUG, EAST and WEST reduced tungsten sputtering sources were observed following several discharges with SBI. Extrapolation of these SBI results are presented to estimate the amount of boron needed for wall conditioning in ITER. Real-time SBI control requirements and plasma operation scenarios for ITER are also described.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Overview of the KSTAR experiments toward fusion reactor

The Korean Superconducting Tokamak Advanced Research has been focused on exploring the key physics and engineering issues for future fusion reactors by demonstrating the long pulse operation of high beta steady-state discharge. Advanced scenarios are being developed with the goal for steady-state operation, and significant progress has been made in high ℓ i , hybrid and high beta scenarios with β N of 3. In the new operation scenario called fast ion regulated enhanced (FIRE), fast ions play an essential role in confinement enhancement. GK simulations show a significant reduction of the thermal energy flux when the thermal ion fraction decreases and the main ion density gradient is reversed by the fast ions in FIRE mode. Optimization of 3D magnetic field techniques, including adaptive control and real-time machine learning control algorithm, enabled long-pulse operation and high-performance ELM-suppressed discharge. Symmetric multiple shattered pellet injections (SPIs) and real-time disruption event characterization and forecasting are being performed to mitigate and avoid the disruptions associated with high-performance, long-pulse ITER-like scenarios. Finally, the near-term research plan will be addressed with the actively cooled tungsten divertor, a major upgrade of the NBI and helicon current drive heating, and transition to a full metallic wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Helium plasma operations on ASDEX Upgrade and JET in support of the non-nuclear phases of ITER

For its initial operational phase, ITER has until recently considered using non-nuclear hydrogen (H) or helium (He) plasmas to keep nuclear activation at low levels. To this end, the Tokamak Exploitation Task Force of the EUROfusion Consortium carried out dedicated experimental campaigns in He on the ASDEX Upgrade (AUG) and JET tokamaks in 2022, with particular emphasis put on the ELMy H-mode operation and plasma-wall interaction processes as well as comparison to H or deuterium (D) plasmas. Both in pure He and mixed He + H plasmas, H-mode operation could be reached but more effort was needed to obtain a stable plasma scenario than in H or D. Even if the power threshold for the LH transition was lower in He, entering the type-I ELMy regime appeared to require equally much or even more heating power than in H. Suppression of ELMs by resonant magnetic perturbations was studied on AUG but was only possible in plasmas with a He content below 19%; the reason for this unexpected behaviour remains still unclear and various theoretical approaches are being pursued to properly understand the physics behind ELM suppression. The erosion rates of tungsten (W) plasma-facing components were an order of magnitude larger than what has been reported in hydrogenic plasmas, which can be attributed to the prominent role of He 2+ ions in the plasma. For the first time, the formation of nanoscale structures (W fuzz) was unambiguously demonstrated in H-mode He plasmas on AUG. However, no direct evidence of fuzz creation on JET was obtained despite the main conditions for its occurrence being met. The reason could be a delicate balance between W erosion by ELMs, competition between the growth and annealing of the fuzz, and coverage of the surface with co-deposits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

PFPO plasma scenarios for exploration of long pulse operation in ITER

Long Pulse Scenarios (LPS) in ITER foreseen during the Pre-Fusion Power Operation (PFPO) phase of the ITER Research Plan (IRP) are assessed using 1.5D transport simulations within the ASTRA framework. Such assessment is required to predict the operational space for LPS operation in PFPO, as well as to evaluate which physics processes for LPS operation during Fusion Power Operation (FPO) could be studied during PFPO. An important aspect in the development of LPSs in PFPO is to minimize lifetime consumption of the Central Solenoid (CS) for these scenarios. The maximum pulse length achievable for LPSs in PFPO with no consumption of CS lifetime (currents in CS coils $&LessSlantEqual;$30 kA per turn) has been assessed for a range of heating schemes and heating mixes, confinement regimes (L-mode and H-mode) and for helium and hydrogen plasmas. The operational space of LPS and pulse length has been explored through density scans with the Heating and Current Drive mix required for the FPO Q $&GreaterSlantEqual;$ 5 steady-state plasma scenario (namely Neutral Beam Injection and Electron Cyclotron Heating) including acceptable shine through losses on the first wall for both helium and hydrogen plasmas. Fast particle physics aspects that are common between FPO plasmas and LPS PFPO H-mode plasmas at low densities are studied including MHD stability analysis with the KINX code and non-perturbative critical gradient model based on high-n Toroidal Alfven Eigenmodes (TAE) stability kinetic ballooning code HINST calculations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Toward holistic understanding of the ITER-like resonant magnetic perturbation (RMP) ELM control on KSTAR

Here, KSTAR has clarified a set of unresolved 3D physics issues utilizing the ITER-like in-vessel, three-row, resonant magnetic perturbation (RMP) configurations. Since RMP-driven, edge-localized-modes (ELMs)-crash control elevates the divertor heat flux peak through its impact on edge plasma parameters and transport, a series of intentionally misaligned RMP configurations (IMCs) have been explored to investigate the relationship between RMP ELM control and divertor heat fluxes, while searching for an ideal IMC that could be favorable in both aspects. First of all, the contrasting influence of kink vs anti-kink phasing on the ELM-crash suppression has been articulated, demonstrating the synergistic benefit of 'kink' phasing on ELM-crash-suppression. On the other hand, the three-row IMC in the anti-kink phasing becomes more insensitive to the ELM-crashes at the sub-marginal level of RMP, consistent with theory. Meanwhile, the divertor 'wetted' area of ELM-crash-suppression gets narrower than that of ELM-crash-mitigation, suggesting that ELM-crash-mitigation remains advantageous over ELM-crash-suppression in terms of time-averaged divertor thermal loading. In comparison, based on a set of two-row IMCs, no evidence of divertor heat flux broadening was found during ELM-crash-suppression, supporting a hypothesis that the dispersal of the divertor heat flux in three-row IMCs cannot be driven by helically structured two-row RMPs alone. Among ITER-like three-rows, lower two-row RMPs have been found to be much more effective in suppressing the ELM-crashes than upper two-row RMPs. Although it is quite preliminary, the up/down asymmetric dependence of RMP coupling may be generically attributed to lower-single-null plasmas. Such a holistic understanding of RMP-driven, ELM-crash-control in KSTAR is expected not only to elucidate various subtle points in the vicinity of ELM-crash-suppression, but also to clarify the relevant divertor thermal loading issues for ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Loss of energetic particles due to resistive wall mode instability in ITER

Effects of an unstable n = 1 (n is the toroidal mode number) resistive wall mode (RWM) on the energetic particle (EP) confinement and loss are numerically investigated, for an ITER steady state scenario with 10 MA plasma current and 5.3 T toroidal field. The eigenfunction of the RWM is computed, with the associated three-dimensional magnetic field perturbation superposed with the 2D equilibrium field for tracing the EP drift orbits. Considered are mono-energetic EPs at 0.5 MeV and 1 MeV for deuterium ions, and 3.5 MeV for fusion-born alphas, with a range of distribution in the particle pitch angle. Modeling finds that less than 20% of EPs can be lost to the limiting surface in ITER assuming a source distribution uniform in minor radius, due to an unstable RWM that produces 100 Gauss poloidal field perturbation at the outboard mid-plane just inside the (effective) resistive wall surface. On top of the initial prompt drift orbit loss for counter-current EPs, the RWM induced particle loss occurs on a one second time scale, which is comparable to the RWM growth time in ITER. The 'wetted' area, due to the lost EPs striking the limiting surface, is generally found to be large due to the RWM. This is a favorable prediction for ITER. Here, the loss distribution in the poloidal angle is more uniform for co-current EPs. Counter-current EPs experience outward orbit drift when launched from the low-field side, and tend to more often hit the bottom region of the limiting surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Divertor detachment in the pre-fusion power operation phase in ITER during application of resonant magnetic perturbations

Detachment of the divertor plasma during application of resonant magnetic perturbation (RMP) fields is evaluated for hydrogen H-mode plasma during the first pre-fusion power operation (PFPO-1) phase in ITER by 3D plasma boundary modelling with EMC3–EIRENE. Plasma response effects from a linearized, resistive, single fluid MHD model are discussed, which includes partial screening of the externally applied field—but also field amplification near the separatrix. Furthermore, this field amplification is found to play a pivotal role for the magnetic footprint on the divertor targets, but is sensitive to model parameters. Extensions of the footprint beyond the straight portions of the ITER vertical divertor targets, optimized for high stationary heat flux handling, may be possible depending on the level of toroidal rotation in the plasma. Exhaust from the bulk plasma is guided by the helical corrugations (lobes) of the perturbed separatrix, and this results in an upstream heat flux that is distributed over these lobes with lower peak values than in the typical radial heat flux profiles seen in the absence of magnetic perturbations. As a consequence, an earlier onset (with respect to the upstream density) of detachment is found in the traditional strike zone when RMPs are applied, but secondary, non-axisymmetric strike locations appear—and those remain attached at temperatures above 10 eV. Neon seeding can mitigate these non-axisymmetric heat loads, but this becomes less efficient for large magnetic footprints.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Drift orbit islands of energetic particles due to 3D fields in ITER

In this study, we report a comparative study of the orbit islands formed by the guiding center drift motion of test energetic particles (EPs), including both high-energy deuterium ions and fusion-born α-particles, with respect to the magnetic islands formed by tracing the field lines, for an ITER plasma representing the 15 MA baseline scenario. The particle drift orbit is modified by the presence of static resonant magnetic perturbations (RMPs) of different toroidal mode numbers n (n = 1 and 3 in this study), forming orbit islands in the Poincaré plane by passing EPs. The key findings are as follows. (i) With an n = 1 RMP field, the size of the orbit islands combined with the total RMP field, including the plasma response, is about three times smaller than that obtained by assuming the corresponding vacuum field. (ii) The orbit island size is not sensitive to the EP energy, and is comparable to those of the magnetic islands. (iii) Passing EPs with outward radial orbital drift form orbit islands that shift inward with respect to the corresponding magnetic islands, and vice versa. This radial shift, ΔΨ, measured in the normalized equilibrium poloidal flux, is quantified as ΔΨ ≃ 0.03X for the ITER baseline plasma and assuming n = 1 RMPs, with $X\equiv {(M/{M}_{\mathrm{p}})}^{1/2}{(E[\mathrm{M}\mathrm{e}\mathrm{V}])}^{1/2}{Z}^{-1}$, where M is the particle mass, Mp is the proton mass, E is the particle energy, and Z is the particle charge number. (iv) Trapped EPs do not form drift orbit islands, even in the presence of the 3D fields produced by the 90 kAt ELM control coil current in ITER, and thus possess good confinement properties. Nevertheless, the deformation of the trapped particle drift orbits in the Poincaré plane shows that the canonical toroidal angular momentum of EPs is no longer conserved in the presence of the RMP fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Property of neoclassical GAMs induced by pellet generated plasma perturbations in the gyrokinetic code XGC

A pellet source model has been implemented into the gyrokinetic code XGC and applied to neoclassical geodesic acoustic mode (GAM) study as the first step to investigate potentially rapid kinetic spread of pellet-born particles following their injection into plasma. Here in this study, GAM oscillations of the radial electric field, as well as of their effect on the radial particle and energy fluxes, are studied and the observed frequencies agree with theoretical expectations. It is verified that (i) GAM oscillations are driven both within the pellet-fuelled region and for locations at a larger radius than this region and not for locations radially inside the pellet-fuelled region, and that, (ii) even though the poloidal spread along the magnetic field lines is fast, the time-averaged radial plasma transport under GAM oscillations is kept to roughly the same level as the neoclassical transport in the absence of plasma turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Multi-device Studies of Pedestal Physics and Confinement in the I-mode Regime

This paper describes joint ITPA studies of the I-mode regime, which features an edge thermal barrier but with L-mode-like particle and impurity transport, and without ELMs. The regime has been demonstrated on Alcator C-Mod, ASDEX Upgrade and DIII-D tokamaks, over a wide range of device parameters and pedestal conditions, including low collisionality. Access to I-mode is independent of heating method (neutral beam injection, ICRH and/or ECRH). Normalized energy confinement H98y,2 ≥ 1 has been achieved for a range of 3≤q95≤4.9 and scales favourably with power. Changes in turbulence in the pedestal region accompany the transition from L-mode to I-mode. The L-I threshold increases with plasma density and current, and with device size, but weakly with BT. The upper limit of power for I-modes increases with toroidal field, and is largest on Alcator C-Mod at B > 5 T. Issues for extrapolation to ITER and other future fusion devices are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Alcator C-Mod: Research in support of ITER and steps beyond

This paper presents an overview of recent highlights from research on Alcator C-Mod. Significant progress has been made across all research areas over the last two years, with particular emphasis on divertor physics and power handling, plasma-material-interaction studies, ELM-suppressed pedestal dynamics, core transport and turbulence, and RF heating and current drive utilizing Ion Cyclotron and Lower Hybrid tools. Specific results of particular relevance to ITER include: inner wall SOL transport studies that have led, together with results from other experiments, to the change of the detailed shape of the inner wall in ITER; runaway electron studies showing that the critical electric field required for runaway generation is much higher than predicted from collisional theory; core tungsten impurity transport studies reveal that tungsten accumulation is naturally avoided in typical C-Mod conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Detachment in Fusion Plasmas with Symmetry Breaking Magnetic Perturbation Fields

Power exhaust from the bulk plasma is significantly altered by symmetry breaking magnetic perturbation fields, because these create direct connections (perturbed field lines) from the confined high temperature plasma to solid surfaces. Additionally, the same amount of power is distributed among those new exhaust channels as for a symmetric magnetic configuration, which reduces the local upstream heat flux flowing down the perturbed field lines, thereby making access to detachment easier (i.e., at lower upstream density) for the divertor plasma near the location corresponding to the symmetric magnetic separatrix. However, the divertor plasma regions with connection to the bulk plasma are extended nonaxisymmetrically further outside, where significant heat loads occur, unlike in the symmetric configuration. The temperature remains high at those locations, which reduces the divertor plasma dissipation capacity, making the mitigation of heat loads more difficult to achieve.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

ELM control optimization for various ITER scenarios based on linear and quasi-linear figures of merit

For the purpose of better understanding type-I edge localized mode (ELM) control in ITER with resonant magnetic perturbation (RMP) fields, the plasma response to RMP is computed by a resistive full magneto-hydrodynamic model in toroidal geometry. Five scenarios designed for ITER are considered, ranging from the pre-nuclear to nuclear phases. The plasma response to RMP is quantified by the plasma surface displacement near the X-point of the divertor plasma and at the outboard mid-plane. The optimal coil configurations between two high- Q deuterium-tritium (DT) scenarios (at the same plasma current of 15 MA and the same magnetic field of 5.3 T but different fusion gains, Q = 5 and 10) are predicted to be similar. For the other ITER scenarios with similar edge safety factor q 95 ~ 3 to that of the baseline scenario, the optimal coil phasing is also similar. The optimization results are different for a half-current full-field (7.5 MA/5.3 T) scenario, largely due to the difference in q 95 . The RMP coil currents are also optimized to tailor the core vs edge toroidal torques exerted by the 3D RMP fields on the plasma column. Torque optimization, with various objective functions proposed in the study, is useful for minimizing the side effects of RMP on the plasma core flow in ITER, while still maintaining the ELM control capability. Full utilization of three rows of ELM control coils in ITER is found to be essential to ensure both flexibility and robustness of ELM control, in terms of both linear and quasilinear plasma responses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗