1. Introduction
The main goal of nuclear fusion is to generate electric energy, reproducing on earth the nuclear reactions that power the stars. Nuclear fusion is among the energy sources able to guarantee world sustainability without CO2 production, contributing to meeting the quick growth of global energy demand, which is expected to be doubled by 2050 due to the combined effect of the increase of world population and energy need per person in developing countries. However, different challenges must be addressed to make fusion a robust form of electric energy production.
Magnetically confined fusion uses a magnetic field to maintain the plasma (fully ionized gas) away from the wall. The most promising configuration of a magnetic fusion device is the tokamak, a toroidal device in which plasma plays the role of the secondary winding of a transformer that has to be kept in equilibrium by a set of poloidal field currents [
1].
In 2012, the EUROFUSION consortium published the “European Fusion Roadmap” [
1], in which eight missions have been identified and addressed to realize nuclear fusion by the end of the century. In particular, mission 2 refers to the power exhaust problem, which probably represents the main challenge towards the realization of magnetic confinement fusion. According to the roadmap, the “risk exists that the baseline strategy pursued in ITER cannot be extrapolated to a fusion power plant”, hence “a dedicated Divertor Tokamak Test (DTT) facility will be necessary”. The final aim of the DTT device, whose construction is starting at ENEA C.R. Frascati—Italy, is the test of different solutions for the divertor in view of a DEMOnstration fusion reactor (DEMO), such as alternative plasma configurations and innovative materials for the divertor able to cope with very large power fluxes (>10 MW/m
2). In 2015, the first proposal of the DTT project was published [
2], while in 2019, a revised version was released [
3] after a careful revision of the project according to the suggestions of the EUROFUSION consortium.
This paper is focused on the description of DTT plasma scenarios, showing a revised version with respect to [
4], while in [
5], a detailed description of the Single Null (SN) scenario fluid modeling is provided. Even if a DTT PF (Poloidal Field)/CS (Central Solenoid) coil system is able to accommodate a comprehensive set of divertor configurations, including SN, X-Divertor (XD), Single Null with Negative-Triangularity (SN-NT), Double Null (DN) and Snowflake (SF), for the design of the DTT first divertor, priority has been given to the SN (reference configuration), XD and SN-NT scenarios [
6,
7].
The paper is organized as follows.
Section 2 presents the DTT machine configuration, showing the position, dimension and number of turns of the poloidal field coils, highlighting the limits in terms of magnetic fields, forces, voltages and currents.
Section 3 reports the plasma scenario for the SN configuration, including the ramp-up, flat-top and ramp-down phases.
Section 4 shows the alternative configurations for a flat top, whereas
Section 5 draws the main conclusions of the work.
2. DTT Machine Configuration and Constraints
In this section, the present design of a DTT device is shown, and the main constraints considered for the development of the plasma scenarios are summarized.
2.1. Machine Description
The current DTT design, whose poloidal section is illustrated in
Figure 1, foresees a major radius
= 2.19 m, a minor radius
= 0.70 m with a maximum plasma current
= 5.5 MA and a toroidal magnetic field
= 6 T at
R = 2.14 m.
DTT is characterized by two main passive structures, namely the vacuum vessel and the stabilizing plates. The vacuum vessel consists of two 15 mm thick shells in stainless steel 316LN. The in-vessel stabilizing plates consist of two loops with a thickness of 30 mm made of copper on the upper and lower side of the vessel, which are connected by a bridge so that the currents in the lower and upper part are anti-parallel and it only reacts to the net flux change of the plasma. The stabilizing plates allow DTT to have elongated plasmas with a relatively low growth rate of the vertical instability, with tolerable effects during the breakdown phase and robustness in the case of plasma disruptions.
The poloidal field magnetic coil system of DTT is characterized by 12 superconducting coils plus 6 in-vessel copper coils. The Central Solenoid (CS) is characterized by six modules, each of them consisting of the series of three conductors at high, medium and low fields realized in
and capable of withstanding a peak magnetic field of
= 13.6 T in the inner conductor. The current limit in the CS modules is
= 31.34 kA. The six PF coils are realized in
(PF1 and PF6) and
(PF2, PF3, PF4 and PF5). In
Table 1, the peak fields and the maximum currents of each PF coil are reported.
DTT is equipped with six in-vessel copper coils. Two up-down symmetric coils (IVCU and IVCL) are dedicated to the vertical stabilization of the plasma column and to the fast radial control of plasma centroid. Four additional coils (IVC1, IVC2, IVC3 and IVC4) located in the divertor region are dedicated to the control of quantities related to power exhaust and to the implementation of control strategies, such as strike-points sweeping [
8]. It is important to remark that the nominal plasma scenarios are obtained without the contribution of the in-vessel coils.
Table 2 reports the coordinates, the dimension and the number of turns of the superconducting magnets.
2.2. Plasma Scenarios Constraints
For the design of the plasma scenarios, the following constraints have been considered.
2.2.1. Vertical Forces
The force limits adopted on the CS/PF coils are:
Maximum vertical force on the CS stack not exceeding 20 MN;
Maximum separation force in the CS stack not exceeding 30 MN;
Maximum vertical force not exceeding 30 MN for PF1/PF6;
Maximum vertical force not exceeding 40 MN for PF2/PF5 and PF3/PF4.
2.2.2. Magnetic Fields
The maximum magnetic field at the location of the CS coils shall not exceed 13.6 T in the high field section, 10.8 T in the medium field section and 8.9 T in the low field section. The maximum magnetic field at the location of the PF coils changes as a function of the radial coordinate of the coils. Peak magnetic field values are reported in
Table 1.
2.2.3. Plasma
Minimum distance of 40 mm between the plasma and the inner side of the first wall;
Minimum distance of 60 mm between the plasma and the outer side of the first wall;
The maximum plasma current is 5.5 MA, whereas the plasma shape parameters should be similar to the present EU design of DEMO, with a (plasma triangularity) and (plasma elongation) ; with these parameters, taking into account the presence of the passive conductor, the plasma control system is able to stabilize the plasma vertically by means of the in-vessel coils.
2.2.4. CS/PF Voltages
The plasma scenarios have been designed according to the following voltage limits for CS/PF coils:
Maximum voltage of 800 V for CS circuits;
Maximum voltage of 2 kV for PF1 and PF6;
Maximum voltage of 3 kV for PF2, PF3, PF4 and PF5.
3. SN Plasma Scenario
The reference SN scenario is characterized by a 5.5 MA plasma current at flat-top and a total duration of about 100 s. In the following, the details of the pre-magnetization, the ramp-up, flat top and ramp-down phases are presented. The plasma configurations have been produced by the CREATE-NL code [
9] and optimized using the CREATE-L code [
10].
3.1. Pre-Magnetization
A good magnetic field null during the plasma breakdown (i.e., a large central hexapolar region even at the low field) can also be guaranteed in the most demanding case (ohmic breakdown without EC assistance) with a toroidal electric field of 0.8 V/m. This configuration can be obtained with various combinations of currents satisfying the constraints of
Section 2.2, starting from a premagnetization characterized by a poloidal flux of 16.2 Vs (
Figure 2).
3.2. Ramp-Up Phase
The ramp-up phase has been designed according to premagnetization flux of
= 16.2 Vs (
Figure 2) and applying the Ejima formula with
= 0.35 [
11]. Indeed, the flux at the start of the flat-top (SOF) turns out to be:
After the breakdown phase, with a hexapolar null point in the center of the chamber, at
= 0.2 s a first circular plasma with
= 200 kA,
(poloidal beta) = 0.1 and
(internal inductance) = 0.9 is obtained. Progressively, the plasma current
is assumed to increase with a ramp rate of 200 kA/s, while the plasma shape evolves from a circular to an elliptical shape leaning on the inner side of the first wall. Between
= 9.2 s and
= 14.2 s, the plasma current rises up to 3 MA, achieving the X-point formation keeping
and
constant. In this ramp-up scenario, the plasma remains in a limited configuration for about 15 s. From
= 14.2 s to
= 27 s, the current reaches its flat-top value of 5.5 MA, and the plasma shape parameters achieve their target values while
remains at a low value. It is worth mentioning that, according to the voltage limits presented in
Section 2.2, the ramp-up phase can be realized with a maximum plasma current rate value of about 650 kA/s due to the 800 V limitation of CS power supplies.
Figure 3 shows the time evolution of plasma shapes during the ramp-up phase, while
Table 3 provides the evolution of CS/PF currents. In
Figure 4, the time traces of CS/PF voltages during the ramp-up phase are reported, and in
Table 4, the time evolution of main plasma parameters, such as radial and vertical position of plasma centroid (
and
), boundary flux (
), safety factor (
), plasma perimeter and volume (
and
), is reported.
3.3. Flat-Top Phase
Once the current ramp-up phase is terminated at = 27 s, the plasma current has reached the reference flat top value of 5.5 MA while the poloidal beta is still 0.1.
The transition from low (L-mode) to high (H-mode) energy confinement time is obtained by introducing thermal energy into the plasma above a certain threshold value. The back H–L transition is also observed in some cases. These transitions are characterized by a significant change of .
Waiting for a more detailed description of the L–H transition, we assume that at
= 36 s full additional heating is applied, causing an increase of the internal kinetic energy on a time scale longer than the plasma energy confinement time with an increase of
from 0.1 to the target value of 0.65. After
= 36 s, all plasma physical parameters are assumed to remain nearly constant up to the end of the plasma current plateau at
= 84 s. The duration of the flat-top is determined by the estimation of the loop voltage that in DTT is fixed at 0.16 Vs, assuming full additional power. It is worth mentioning that it would be possible to reach an end of flat-top configuration (EOF) at
= 89 s, fulfilling all the engineering constraints reported in
Section 2, but a margin needs to be stored to allow a safe H–L transition. In
Figure 5, the time evolution of plasma shapes during the flat-top phase is illustrated, while in
Figure 6 the behavior of voltages in the flat-top phase is presented.
Table 5 and
Table 6 report the behavior of currents and main plasma parameters.
3.4. Ramp-Down Phase
For the definition of the plasma snapshots during the ramp-down phase, the H–L transition has been modeled as a drop of from the flat-top value of 0.65 to 0.1 within a time interval of 1 s with fixed plasma current, assuming a mean value of among flat top case (0.16 Vs) and the ramp-down case (0.8 Vs). The main assumptions used for the design of the plasma ramp-down phase are hereafter reported:
We assume a decrease of
from
V at 5.5 MA to
V at 1 MA (75% drop during ramp-down). In DTT, we have assumed three possible ramp rates for the ramp-down phase of 250 (natural ramp rate for DTT), 500 and 750 kA/s. In
Figure 7, the behavior of the boundary flux
, according to Formula (2), for different plasma current ramp rates is illustrated. For the sake of simplicity, in the following, we will focus on ramp-down with a 250 kA/s plasma current ramp rate.
After the H–L transition, the plasma volume gradually reduces while
is kept constant to 0.1 and
linearly increases. Plasma shape evolves from an elliptical shape to a circular one, and in the time interval between
97 and
101 s, the divertor to limiter transition occurs. After
103 s, the plasma discharge is terminated. For future development, the possibility to delay the divertor–limiter transition is under analysis. In
Figure 8 and
Figure 9, the time evolution of plasma shapes and CS/PF voltages during the ramp-down phase are reported, respectively, while in
Table 7 and
Table 8, the behavior of currents and plasma parameters are presented.
4. Alternative Plasma Configurations
The flexibility of the DTT machine allows it to produce a wide range of alternative configurations, which increase divertor radiation without excessive core performance degradation, to study different solutions to the power exhaust problem. DTT tokamak is, in fact, capable of realizing XD at 4.5 MA, SN-NT at 4 MA, DN at 5 MA and SN with a long leg at 5.5 MA.
Impurity seeding is likely to be required for dissipation of the high values of power flux both for SN [
5] and alternative configurations [
12]. Impurity seeding increases the effective ion charge Z
eff, which, in turn, has an effect on the flat-top duration.
Power flux density reduction on the divertor target is obtained by various means, e.g., plasma detachment, flaring of the magnetic field lines and strike point sweeping. Plasma detachment is characterized by the dissipation of power, momentum and particle flux along the open field lines from the midplane to the divertor, where the neutral particles are not obliged to spiral around magnetic field lines. The XD concept [
13] seeks to flare flux surfaces near the divertor target by placing a second X-point behind the divertor plate, with a consequent increase of flux expansion. This flaring may introduce a mechanism that counteracts the upward movement of the detachment front and result in a more robust detachment [
14], similar to the stability of a stationary shockwave in a supersonic nozzle flow. The DTT device is able to develop an XD divertor scenario at 4.5 MA characterized by a flat-top duration of 40 s. Furthermore, since DTT is equipped with four coils located in the divertor region, the local magnetic topology can be easily modified (e.g., the distance between the null points) in order to study possible beneficial effects on power exhaust-related quantities.
The main feature of SN-NT [
15] configuration is to reach H-mode-level confinement, maintaining a low edge pressure so as to have ELM-free regimes. Encouraging results have been obtained on TCV [
16] and DIII-D [
17], with a consequent increasing interest towards the negative triangularity concept. In the DTT device an SN-NT scenario at 4 MA characterized by an upper triangularity
= −0.30 and a flat-top duration of 70 s can be realized. According to the particular topology of the SN-NT plasma shape, a dedicated divertor able to accommodate such a configuration must be realized.
The DN configuration is characterized by a second null point located in the upper part of the main chamber. This diverts a significant fraction of the heat load to the inner divertor of a SN to a second target at a larger radius, which increases the wall interaction area and decreases the peak heat load reaching the targets. However, it also decreases the connection length to the target. In DTT tokamak, it is possible to realize a DN configuration at 5 MA thanks to the up-down symmetry of the CS/PF coils system. It is worth mentioning that, in order to realize the DN configuration, an upper divertor is mandatory.
The concept of SN with a long-leg is very similar to the Super-X (SX) [
18], which aims to increase the total flux expansion by increasing the major radius of the divertor target. SN with long-leg at 5.5 MA has been realized in DTT by modifying the SN configuration. In particular, the outer leg length has been increased up to 55 cm while the upper part of the plasma has been moved toward the first wall in order to keep
= 1.66.
In
Figure 10, the SOF plasma shapes of alternative configurations of DTT tokamak are shown, while in
Table 9, the main plasma parameters are reported.
5. Conclusions
In this paper, the SN and alternative plasma configurations for DTT tokamak have been presented. The SN plasma scenario has been analyzed in detail, showing the ramp-up, flat-top and ramp-down phases. The ramp-up phase has been optimized with a plasma current rate variation of 200 kA/s reaching the SOF configuration at
= 27 s. However, a maximum plasma current ramp of 650 kA/s can be imposed, fulfilling all engineering constraints. Concerning the flat-top phase, to realize an H-L transition compatible with the coil voltage limits, the EOF equilibrium has been fixed at
= 84 s (even the flat top might be extended up to
= 89 s). Regarding the ramp-down phase, according to the constraints of
Section 2, the current can decrease at a rate from 250 to 750 kA/s.
Furthermore, alternative plasma configurations have been presented, proving the good flexibility of the DTT device, which is capable of obtaining DN at 5 MA, XD at 4.5 MA, SN-NT at 4 MA and SN with a long leg at 5 MA. Although these configurations are feasible from an electromagnetic point of view, some modifications must be applied to the machine geometry in order to accommodate such plasma shapes. The need to develop and study alternative divertor concepts is the main feature of DTT, whose aim is to test different solutions in terms of power exhaust in view of DEMO design and implementation.