Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
10-15 May 2021
Virtual Event
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Theory-based models for the control of W7-X divertor plasmas

13 May 2021, 14:00
4h 45m
Virtual Event

Virtual Event

Regular Poster Magnetic Fusion Experiments P6 Posters 6


Andreas Dinklage (Max-Planck-Institut für Plasmaphysik)


Both the exploration of previously unattained plasma conditions in stellarators and high-reliability requirements (for safe, long-pulse operation at higher plasma performance) pose new challenges for plasma control in Wendelstein 7-X (W7-X). With upgrades being installed at present, the development of predictive tools for safe long-pulse operation of the optimized stellarator W7-X becomes more and more an urgent task. But the pace for the achievement of reliable plasma operation enjoying full benefit from inherent advantages of stellarators, i.e.~quiescent, stable steady-state operation will critically depend on control capabilities. This paper describes a novel approach for stellarators to develop theory-based, simplified models for control - tailored to stellarator specific requirements. A result of the proposed approach to control localized heat loads on the divertor shows the feasibility of the concept. Since the models originate from theory or scalings, it is expected that they can be applied to future experiments. Therefore, the validity range of the models are explored and result in a comprehensive characterization of physics regimes attained in W7-X so far. Both aspects of this study – validated physics models for operation and the characterization of physics regimes in optimized stellarators - contribute to the physics basis for potential future fusion power plant based on optimized stellarators.

Wendelstein 7-X (W7-X) is a superconducting magnetic confinement fusion device built as a proof-of-concept to assess the reactor potential of optimized stellarators. The plasma volume ($30m^{-3}$) is large enough to expect temperatures leading to low collision frequencies (also collisionalities $\nu^*$). In 3D-fields of classical stellarators, $\nu^*$ as low as for fusion leads to unacceptable losses. The optimized stellarator W7-X is built to demonstrate that detrimental low-$\nu^*$ effects can be mitigated by optimized magnetic fields. In its forthcoming campaigns (“OP2”), W7-X aims to demonstrate safe and reliable steady-state operation. Equally important, the plasma beta and collisionalities should meet requirements that allow one to examine the physics of fusion-grade plasmas. At present, the machine is being equipped with water-cooled plasma-facing components. Heating and fueling upgrades are being installed along with extended wall surveillance and measures to detect high, localized power loads.

Operation space of W7-X test-divertor discharges (OP1). $\nu^*$ are collisionalities, $\langle \beta \rangle$ is the volume averaged plasma beta.

Results shown in this paper are based on systematic transport modelling of W7-X discharges from recent campaigns with an uncooled divertor (“OP1”: limited in pulse lengths ($\tau_{pulse} < 100$s) and heating power $P<7\mathrm{MW}$). The modelled discharges encompass a density range from $2\times 10^{19}\mathrm{m}^{-3}$ to $1.2\times 10^{20}m^{-3}$. The plasmas were heated with electron cyclotron heating (ECH) in X2- and O2-polarization. Fig. 1 summarizes the analyzed data representatively covering discharges in the first campaigns. The representation is chosen to address electron-ion equilibration ($T_e/T_i$), long-mean-free-path conditions (collisionality $\nu^*$) and plasma pressure ($\beta$) since those quantities are key for scenarios in the forthcoming campaigns of W7-X. The arrows indicate gaps from conducted experiments to target values for high performance in W7-X. Fig. 1 clearly reveals different stellarator-specific transport regimes. Core-electron-root-confinement (with positive radial electric fields and $T_e < T_i$ require low electron collisionalities and depend on the field ripples as predicted. Ion-root conditions were attained by density increase but not yet entering low-$\nu^*$ regimes due to a lack of ion heating. At high densities, divertor detachment shows up. Highest $\beta$ plasmas were obtained transiently with pellet injection.

Signatures for stellarator specific, long-mean-free-path transport is consistently found below a threshold of thermal collisionality in the plasma center of about $\nu^*(0) < 0.01$. Electron root conditions (size of the circles in Fig. 1 shows the extent of the positive radial electric fields in the plasma center) are obtained at high electron temperatures but lacking thermal equilibration. The transport modelling reveals that neoclassical energy losses are smaller than 40% of the total energy losses even in the core (electron-root: $\rho < 1/3 \ldots 1/2$). Anomalous transport, radiation and charge-exchange losses (showing up at the edge and SOL) dominate the energy transport otherwise. In ion-root conditions, neoclassical contributions to fluxes are small (as expected at higher $\nu^*$).

Profiles of discharge 20180927.033 at $t_{ref} = 3.5\mathrm{s}$. (TS: Thomson scattering, XICS: x-ray imaging spectroscopy). (d) and (e) are results from a neoclassical transport modelling code.

In contrast to the radial energy flows, neoclassical mechanisms apparently dominate the parallel transport. The time response of the plasma current is modelled by an L/R-response model ($I(t) = I^\infty[1-\exp(-t/\tau_{L/R})]$). While fits are frequently used to determine the stationary current $I^\infty$ and the L/R-response time $\tau_{L/R}$ from experimental data, the approach in this paper is different: theory and measurements are used to forward model (potentially in real time) $I^\infty$ and $\tau_{L/R}$ in order to predict the plasma current. To this end, stationary bootstrap current profiles were calculated (Fig. 2(e)) from measurements of densities and temperatures (Fig. 2(a) and (b)). The required neoclassical modelling included precalculated, theoretical transport coefficients and the (not intrinsically ambipolar) radial electric field (Fig. 2(c)). The resistivity was determined from the theoretical parallel transport coefficients (Fig. 2(d)), the plasma inductance was estimated from the mean major and minor radius, respectively. Errors (broken lines in Fig. 2) include uncertainties of the profile fit and $Z_{eff}$.

Measured plasma current and forward modelled bootstrap current from profile modelling at 3.5s in the W7-X discharge 20180927.033. Dashed lines indicate expected model uncertainties.

Fig. 3 shows results of the L/R-response model with pre-calculated parameters from neoclassical theory. To demonstrate predictive capabilities, measurements (Fig. 2 (a) and (b)) at a reference time ($t_{ref} = 3.5\mathrm{s}$ early in the discharge) are used for predictive calculations. Good agreement of the model prediction (using information at $t_{ref}$) with the plasma current is obtained for periods $t \gg t_{ref}$. It is concluded that the L/R-response model using measured profiles provides a real-time prediction for plasma control. For discharge simulations (prior experiments), profiles of relevant quantities are reproducibly set by the heating power, the densities and the magnetic configuration. A parametrization of $I^\infty$ and $\tau_{L/R}$ is therefore an option for pre-calculating the plasma current response or for simplified (but life) measurements of profile proxies.

Limitations of the model are found in transient phases of plasma current redistribution, e.g. during plasma start-up or changes of the heating power and current drive with minor effects on the long-term current evolution. Moreover, if magnetic islands are intentionally shifted inside the plasma, differences are found both in $I^\infty$ and $\tau_{L/R}$ (in the order of a few tens of percent). The difference due to magnetic islands (even beyond the flattening of the pressure profiles) is investigated. For O2-heated, long-pulse discharges, systematic deviations in the predicted current occur. Differences may be due to electron cyclotron current drive (ECCD) from multi-pass reflections of strong microwave heating.

The findings demonstrate the potential of theory-based, but reasonably simple time response models to predict the plasma current. The presented results can be applied for in-situ control of changes in the edge rotational transform (affecting strike-line positions) being computationally much faster than full time-dependent modelling. Current seeding with ECCD could be used as an actuator. Further applications for the general approach lie in the prediction of stored energy and density control. Since the control scheme is proposed to rely on physics-based modelling or scalings, the approach can be extrapolated to upgraded heating capabilities expected in the next campaigns of W7-X. Ultimately, the models therefore contribute to the development of reliable plasma scenarios for next-step stellarators.

Affiliation Max-Planck-Institut für Plasmaphysik
Country or International Organization Germany

Primary authors

Andreas Dinklage (Max-Planck-Institut für Plasmaphysik) Dr Golo Fuchert (Max-Planck Institut für Plasmaphysik) Robert Wolf (Max-Planck-Institute for Plasma Physics) Arturo Alonso (Laboratorio Nacional de Fusión - CIEMAT) Dr Tamara Andreeva (Max-Planck Institut für Plasmaphysik) Dr Craig Beidler (Max-Planck-Institut für Plasmaphysik) Dr Marco de Baar (DIFFER) Dr Yu Gao (Max-Planck Institut für Plasmaphysik) Joachim Geiger (Max-Planck-Institute for Plasma Physics, Greifswald, Germany) Dr Marcin Jakubowski (Max-Planck Institut für Plasmaphysik) Mrs Heike Laqua (Max Planck Institut für Plasmaphysik) Dr Nikolai Marushchenko (Max-Planck Institut für Plasmaphysik) Dr Ulrich Neuner (Max-Planck Institut für Plasmaphysik) Dr Novimir Pablant (Princeton Plasma Physics Laboratory) Dr Andrea Pavone (Max-Planck Institut für Plasmaphysik) Kian Rahbarnia (Max-Planck Institut für Plasmaphysik) John Schmitt (Princeton Plasma Physics Laboratory) Dr Håkan Smith (Max-Planck-Institut für Plasmaphysik) Dr Torsten Stange (Max-Planck Institut für Plasmaphysik) Dr Yuriy Turkin (Max-Planck Institut für Plasmaphysik) W7-X Team (Max Planck Institut für Plasmaphysik)

Presentation Materials