Next Article in Journal
A Review of Experimental Investigations into the Time Evolution of Low-Pressure Capacitively Coupled Plasmas in Their Early Stages of Development
Next Article in Special Issue
Exploring Experimental Isotope Scaling and Density Limit in Tokamak Transport
Previous Article in Journal
Positive- and Negative-Polarity Nanosecond-Pulsed Cryogenic Plasma in Liquid Argon
Previous Article in Special Issue
Spatial Distribution Analyses of Axially Long Plasmas under a Multi-Cusp Magnetic Field Using a Kinetic Particle Simulation Code KEIO-MARC
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of ICRF Heating Schemes in ITER Non-Active Plasmas Using PION+ETS Integrated Modeling

1
Barcelona Supercomputing Center, 08034 Barcelona, Spain
2
Catalan Institution for Research and Advanced Studies (ICREA), 08010 Barcelona, Spain
3
Royal Institute of Technology (KTH), School of Electrical Engineering and Computer Science, Department of Electrical Engineering, Division of Electromagnetic Engineering and Fusion Science, 10044 Stockholm, Sweden
*
Author to whom correspondence should be addressed.
Plasma 2024, 7(3), 517-530; https://doi.org/10.3390/plasma7030028
Submission received: 30 April 2024 / Revised: 8 July 2024 / Accepted: 11 July 2024 / Published: 19 July 2024
(This article belongs to the Special Issue New Insights into Plasma Theory, Modeling and Predictive Simulations)

Abstract

:
The PION code has been integrated into the European Transport Solver (ETS) transport workflow, and we present the first application to model Ion Cyclotron Resonance Frequency (ICRF) heating scenarios in the next-step fusion reactor ITER. We present results of predictive, self-consistent and time-dependent simulations where the resonant ion concentration is varied to study its effects on the performance, with a special emphasis on the resulting bulk ion heating and thermal ion temperature. We focus on two ICRF heating schemes, i.e., fundamental H minority heating in a 4He plasma at 2.65 T/7.5 MA and a three-ion ICRF scheme consisting of fundamental 3He heating in a H-4He plasma at 3.3 T/ 8.8 MA. The H minority heating scenario is found to result in strong absorption by resonant H ions as compared to competing absorption mechanisms and dominant background electron heating for H concentrations up to 10%. The highest H absorption of ∼80% of the applied ICRF power and highest ion temperature of ∼15 keV are obtained with an H concentration of 10%. For the three-ion scheme in 85%:15% H:4He plasma, PION+ETS predicts 3He absorption in the range of 21–65% for 3He concentrations in the range of 0.01–0.20%, with the highest 3He absorption at a 3He concentration of 0.20%.

1. Introduction

Auxiliary heating systems are essential for the plasma in ITER to reach the necessary temperature for fusion to take place. According to the ITER Research Plan in 2018 [1], these systems will be comprised of Ion Cyclotron Resonance Frequency (ICRF) heating, Electron Cyclotron Resonance Heating (ECRH), and Neutral Beam Injection (NBI). The operations at ITER will follow a staged approach, the details of which are currently under review [2]. The 2018 Plan proposed starting with two Pre-Fusion Power Operation (PFPO I and II) phases, where the auxiliary heating systems would be tested on non-active hydrogen (H) and helium (He) plasmas, and continuing with the Fusion Power Operation (FPO) phase, where a transition would be made to deuterium (D), tritium (T), and deuterium–tritium (DT) plasmas at full plasma current and in a toroidal field.
In this paper, we study ICRF heating in ITER non-active plasmas with q 95 = 3 at half field and current (2.65 T/7.5 MA), where q 95 is the value of the safety factor at the radial location where 95% of the magnetic flux is contained within that surface. This heating method uses fast magnetosonic waves that are launched by external antennas, propagated in the plasma and are absorbed by different plasma species. The wave–particle resonance with plasma ions, ω = k v + l ω c , takes place when the wave frequency ω matches the Doppler-shifted ion cyclotron frequency ω c or its harmonic. Here, k is the wavenumber parallel to the background magnetic field, v is the ion parallel velocity, and l = 1 , 2 , 3 represents the harmonics of the wave. The l = 1 resonance corresponds to the fundamental resonance, the l = 2 resonance to the second harmonic resonance, and so on. Typically, only a small number of plasma ions resonate with the launched waves. In most cases, these are either minority ions for l = 1 and majority ions with a finite Larmor radius with respect to the wavelength of the launched wave for l > 1 . In addition, there is direct electron damping of the launched waves by electron Landau damping and transit time magnetic pumping, which competes with ion cyclotron damping [3].
The different options for ICRF heating in ITER plasmas have been investigated in [4]. It was concluded that strong central absorption can be achieved in half-field placing the fundamental resonance of H in the plasma center using an ICRF frequency of 40 MHz in 4He plasmas. Furthermore, it was found that schemes with good absorption in H plasmas are lacking at half-field. To alleviate this situation, it was proposed to use a three-ion ICRF scheme at a higher magnetic field in the range of 3.3 T. In a three-ion scheme, a small concentration of a third ion species Z is added to the plasma that consists of two other ion species Y and X. The fundamental cyclotron resonance layer of the species Z is tuned so that it coincides with the mode conversion layer of the other two species [5]. It was proposed in [5] to add a small amount of 3He to an H plasma with a 4He minority at a magnetic field of 3.3 T, which was shown to give rise to good absorption at the fundamental 3He resonance at 40 MHz.
In this paper, we study the above two heating scenarios, i.e., fundamental H minority heating in H plasma at an ITER half-field of 2.65 T and the three-ion scheme 3He-H-4He in H plasma at a magnetic field of 3.3 T using the ICRF modeling code PION [6] integrated into the transport modeling workflow European Transport Solver (ETS) [7]. The aim is to predict the time evolution of the ICRF heating characteristics and plasma behavior when ICRF heating is applied in these conditions. While PION and ETS have been used separately for plasma heating modeling in the past, this paper reports the first results using the PION+ETS integration. We present results of predictive, self-consistent, time-dependent simulations and carry out H and 3He density scans to determine the effects of the concentration on the time evolution of the absorbed power density profiles, collisional power transfer, and the ion temperature. Particular attention is given to the resulting bulk ion heating and thermal ion temperature given their importance to overall fusion plasma performance and the fact that ICRF heating will be the only auxiliary heating scheme in ITER that can provide dominant bulk ion heating and access to the H mode [8]. The scans carried out in the concentration of resonant ion species are relevant because the concentration of resonant ion species affects ICRF power partitioning between the competing absorption mechanisms, the average energy of the resonant ions and, thereby, the fraction of ICRF power transferred from resonant ions to bulk ions and electrons in collisions with the background plasma. When the resonant ion energy is above the critical energy E crit = 14.8 A T e [ i n i Z i 2 / ( n e A i ) ] 2 / 3 (given in eV), where A is the atomic mass of the resonant species, T e is the electronic temperature (given in eV), n i is the ionic density (given in cm−3), Z i is the atomic number (with the sum over i referring to all thermal ion species) and n e is the electronic density (given in cm−3) [3], the ions are more likely to collide with electrons, resulting in dominant collisional bulk electron heating. If the energy is below the critical energy, bulk ion heating will dominate.
This paper is organized as follows. A brief introduction of the tools used in the modeling is given in Section 2. In Section 3, the results of the simulations are presented and compared to earlier studies using a number of ICRF modeling codes with fixed plasma parameters. Finally, Section 4 contains some final remarks and conclusions.

2. Methodology

The simulations in this paper have been conducted with the PION+ETS integration. In the following, these tools are briefly described.

2.1. PION

PION is a relatively fast ICRF modeling code based on simplified models [6]. It calculates the time evolution of the ICRF power absorption together with the velocity distribution functions of the resonant ions in a self-consistent way. The main features of PION can be divided into two procedures, i.e., the power deposition model and the Fokker–Planck calculations [9]. The power deposition model computes how the ions and electrons absorb the ICRF waves. The power deposition model of PION approximates full-wave codes like LION [10]. It decomposes the wave field into two parts: one part in the strong absorption limit and another part in the weak absorption limit. The time evolution of the velocity distribution function is calculated using a one-dimensional Fokker–Planck equation for the pitch angle-averaged distribution function. Both procedures are coupled by the computation of the dielectric tensor. The Fokker–Planck calculation provides the contributions to the dielectric tensor by the ICRF-accelerated resonant ions. These contributions are used for the calculation of the power deposition and modify the power partitioning and power absorption profiles.
A typical PION simulation involves several time steps. In each time step, the background plasma parameters are read, and the absorbed power deposition is calculated. Then, the result of this calculation is used as the input of the one-dimensional Fokker–Planck solver, which outputs a distribution function. This distribution function is then used to calculate the absorbed power in the next time step [9].
PION has been extensively validated against experimental data on JET [6,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32], AUG [28,33,34,35,36,37,38], DIII-D [39], and WEST (formerly Tore Supra) [40,41] for many minority and majority heating schemes. At JET, it is part of the second data processing chain. Recently, PION was integrated [42] into the ITER Modeling and Analysis Suite (IMAS) [43]. IMAS is the computational platform that supports plasma operations and research activities at ITER. It uses a data model that can describe both experimental and simulated data with the same representation. As part of the IMAS integration, the PION results in IMAS were compared with the results of the standalone PION code (outside IMAS), and overall, good agreement was found. Furthermore, verification of PION’s data handling within the ETS framework was carried out in the initial stages of the integration, including checks on retrieved ion species and flux surface mapping, which influence many other variables in PION such as temperature and density. Regarding the handling of input and output of variables calculated by PION, there are a number of flags, such as in the output files that were used to verify the input and output data.
In this work, we integrate PION into the transport modeling workflow of the European Transport Solver (ETS) within IMAS. This integration is relevant because it provides the capabilities to simulate the evolution of a plasma discharge. Both interpretative and predictive simulations are possible. Furthermore, the integration of PION into IMAS has allowed collective development of integrated modeling tools and workflows.

2.2. ETS

The ETS is one of the transport modeling workflows within the IMAS framework. It is a modular package of physics modules combined into a workflow, where PION has been included as one of the Heating and Current Drive (H&CD) codes. Other heating codes integrated into ETS as H&CD modules include CYRANO [44] or StixRedist [45]. The ETS has been developed with the objective of building the capabilities to compute the full discharge evolution of a power plant at the tokamak scale. To meet this objective, the ETS counts with a high degree of modularity, a separation of physics and numeric parts, a flexible workflow, and the ability to both treat several ion components (including impurities) and use stiff transport models [7]. The ETS functions as a suite of validated codes containing several transport actors integrated into its workflow which use different models, such as neoclassical or gyrokinetic transport models.
The ETS is formed out of coupled codes where the information exchange takes place through well-defined generalized Interface Data Structures (IDSs) acting as standardized interfaces [46]. This is the standard layered structure used in IMAS, where the reactor data model is mapped to the solver data model through the IDSs. The IDSs contain physics definitions of the equilibrium, neoclassical transport and core sources, transport, impurities, and profile. The IDS equilibrium represents the equilibrium at the previous time step, and the IDS core source represents the sources of current, electron and ion energy, momentum, and particles, whereas the IDS core transport represents the transport coefficients of these quantities. In each iteration, the ETS core code reads the input from the equilibrium, core source and transport IDSs, as well as the core profile representation of the plasma state at the previous time step. The ETS then passes the input to the physics modules, which solve the transport equations for the density of current, the ion and electron density, the ion and electron temperature, and the toroidal velocity. The core code will then output a core profile containing a new plasma state consisting of the new plasma composition, ionic and electronic temperatures and densities, and various global, time-dependent quantities calculated from the 1D profiles.

3. Modeling and Results

In this section, we use PION+ETS integration to study H minority heating in 4He plasma and three-ion-scheme H:4He (85%:15%) with a 3He minority at 2.65 T and 3.3 T, respectively, in ITER. As discussed in the introduction, these schemes have been identified [4,5] as promising ICRF schemes for heating ITER non-active plasmas.

3.1. Fundamental H Minority Heating in 4He Plasma at 2.65 T/7.5 MA

For the study of fundamental H ( ω = ω H ) minority heating in ITER 4He plasma, synthetic ITER shot 110005 with a pulse duration of 647 s was used. The plasma was pre-heated using 20 MW ECRF heating. The PION+ETS simulation was started during the flat-top phase, at t i = 300 s, with the plasma and ICRF parameters as shown in Table 1 and Figure 1. An ICRF frequency of 40 MHz was chosen to place the fundamental H resonance in the plasma center. The ICRF power of 20 MW was chosen based on the ITER Research Plan [1]. In the simulation set-up, the temperature evolution of the ions and electrons was set to be predictive, whilst the density evolution was set to be interpretative. The chosen transport model uses a combination of a neoclassical model and a multi-ion model; a simplified analytical expression that uses transport coefficients. The H concentration was constant both in radius and in time and was varied between 1.0% and 10% from simulation to simulation as in [42]. These percentages are calculated in terms of the electron density shown in Table 1. The simulation was performed with a full antenna spectrum discretized by 100 toroidal mode numbers N. Finally, the duration of the simulation was chosen to be 8 s to allow ample time for the evolution of the plasma toward a new steady state.
Figure 2 shows the power density absorption profile as given by PION+ETS at t = 308 s. There are three competing absorption mechanisms, i.e., fundamental H damping, 2nd harmonic 4He damping and direct electron damping. Among them, the fundamental H absorption is the dominating mechanism. As shown in Figure 2, most of the ICRF power is absorbed by resonant H ions independently of the H concentration. The H resonance is located at s = 0.03 on the low-field side (LFS) of the plasma. Here, s stands for the square root of the normalized poloidal flux surface. Figure 2 also shows the single-pass H absorption (SPA) coefficients, which measure how much of the ICRF wave is absorbed the first time it travels through the plasma [3]. The SPA coefficients given here are considered at the ICRF resonance with respect to the dominant toroidal mode number N = 54. As we can see from Figure 2, the SPA on H increases with the H concentration in the range of H concentration considered (1.0–10%). For the lowest H concentration of 1.0%, the SPA on H is low enough (80%) for significant 4He absorption (20% of the ICRF power, cf. Table 2), to take place. The 4He absorption arises because the chosen ICRF frequency also coincides with the 2nd harmonic 4He resonance. Consequently, when H damping decreases as the H concentration is decreased, 4He absorption by majority 4He damping starts to compete with the dominant H damping. As shown by Figure 2, direct electron damping is rather small and also takes place mainly off-axis for all H concentrations considered.
Figure 3 shows the collisional power density transfer profiles at t = 308 s. In minority heating, the ICRH wave accelerates the minority species, in this case H, to non-thermal velocities, larger than the velocities of the thermal H population, thereby building up a fast ion tail in the H distribution function. The heat transfer from this fast ion population to the thermal ions and electrons has a strong dependence on this fast ion tail, as discussed above in the introduction. The average energy of the H ions in the tail of the fast ion distribution function is higher (in the range of 459–1719 keV, cf. Table 2) than their critical energy (158–173 keV) for all H concentrations considered. Therefore, the H ions will collide mainly with the electron population and transfer their energy to them. As the H concentration is increased to 10%, the average energy of the fast H decreases until E fast , H 2.5 E crit , resulting in a balanced power equipartition and P ci / P c 0.5 , which agrees with the equation for total energy given up by the resonant particles transferred into the thermal ions of the plasma in [3].
Figure 4 shows the time evolution of the thermal ion temperature at the location of the ICRF resonance at the normalized flux surface s = 0.03 . The results suggest that at lower H concentrations, the time evolution of the plasma is more nonlinear due to the evolution of the ICRF power deposition in time when the fast resonant ion populations build up. We note, however, that the transient evolution in the first seconds of the simulations shown in Figure 4 depends on the parameters used in the transport model and may not be observed experimentally. Also, it takes a surprisingly long time for the plasma to reach a new steady state. The temperatures at the end of the simulation at t = 308 s are in the range of 11.9–15.0 keV, with the lowest T i = 11.9 keV with an H concentration of 2.5% and the highest T i = 15.0 keV with an H concentration of 10%. As shown in Figure 4, the ion temperature increases as the H concentration increases in the range of 2.5 to 10%. The exception to this trend is the case with an H concentration of 1.0% for which the ion temperature is close to that for an H concentration of 5.0% and higher than that for an H concentration of 2.5%. At the low H concentration of 1.0%, 2nd harmonic 4He damping is significant, as discussed above, and contributes mainly to collisional bulk ion heating. Moreover, it has a more peaked profile around the resonance located at s = 0.03 . This explains its relatively good performance in terms of bulk ion heating and ion temperature.
Prior to this work, the H minority scenario in an ITER 4He plasma was modeled in [4,42,47,48] using fixed background plasma parameters. In [4], good absorption was found with combined SPA coefficients summed over electrons and ion species close to 1 for H concentrations in the range of 2–8%, which is consistent with the SPA values shown in Figure 1. In [42], PION was used as a stand-alone code within IMAS (PION+IMAS) to study H minority heating in 4He plasmas with H minority concentrations ranging from 1% to 9.2%. An increase in H power absorption with an increasing H concentration from 1% to 5% was found, with a slight decrease in H absorption when the H concentration was increased to 10%. In [47], a wider range of H concentrations up to 19% at a lower magnetic field of 2.5 T and a higher electron density of 5.1 × 10 19 m−3 was explored to benchmark full-wave ICRF codes. Significant differences in the ICRF power partitioning between different full-wave codes were reported, with some of the codes suggesting a maximum H absorption around an H concentration of 5%. The simulations carried out by [48] using the 1D TOMCAT code [49] provided power absorption percentages for fundamental H, 2nd harmonic 4He, and direct electron damping at 72%, 6%, and 22%, respectively, for an H concentration of 5%. These were, again, very consistent with the results from PION+ETS, with percentages for those same absorption mechanisms at 78%, 2%, and 20%, respectively, although TOMCAT calculated a slightly higher absorption for the 2nd harmonic 4He at the expense of fundamental H absorption. The ion absorption was notably central, whereas electron absorption tended to be flatter and more off-axis. In this aspect, PION+ETS, PION+IMAS, and TOMCAT coincided. However, these earlier results for fixed plasma parameters in [42,47] are in contrast to the findings of our PION+ETS simulations regarding an optimum H concentration to maximize H absorption. Previous studies reported maximum H absorption below 10%, whereas in PION+ETS simulations, H absorption increased across the range of 1.0–10%, albeit minimally from 5.0% to 10%, reflecting the plasma temperature’s self-consistent evolution with ICRF heating.

3.2. Three-Ion Scheme 3He-H-4He in H Plasma at 3.3 T/8.8 MA

The three-ion scheme 3He-H-4He was studied with the PION+ETS integration using the synthetic ITER discharge 104010 at 3.3 T with the main ion composition H:4He of 85%:15% and a small 3He concentration in the range of 0.01–0.20%. These concentrations were chosen to allow a direct comparison with the modeling performed in [4]. This synthetic discharge has only one time slice in the discharge flat-top phase at t i = 500 s when the plasma is pre-heated using 30 MW of ECRF heating. The plasma density and the ion and electron temperatures are as shown in Table 3 and Figure 5. The PION+ETS simulation was started at t = 500 s, and the duration of the simulation was 5 s to let the plasma evolve towards its new steady state. The fundamental resonance of the 3He minority ( ω = ω 3 H e ) was targeted using an ICRF frequency of 40 MHz and 20 MW of ICRF power. The simulation setup of the ETS transport model is the same as that employed in the case of the H minority scheme at 2.65 T studied above in Section 3.1.
The competing absorption mechanisms as given by PION+ETS are fundamental 3He damping and direct electron damping. Figure 6 shows that the 3He absorption takes place much more off-axis than the power deposition in the H minority scenario at 2.65 T (see Figure 1). The 3He resonance is located in the high-field side (HFS) of the plasma at s = 0.62. The electron damping takes place over the whole radial profile, but it is more prominent on-axis, where the electron temperature and density are highest [50]. In Table 4, it can be seen that most of the power is absorbed by the 3He resonance for the 3He concentrations above 0.05%. Below this concentration, the dominating mechanism is direct electron damping. This can also be understood in terms of the SPA, which is very low (7%) for the 3He concentration of 0.01%, resulting in most of the power in the ICRF wave not being absorbed at the resonance by the resonant 3He ions but at the center of the plasma by the electrons. From the results in Table 4, we can see the trend of increasing power absorbed with increasing 3He concentration, with a maximum of 65% at a 3He concentration of 0.20%. The SPA is also highest (89%) at this concentration.
In terms of the collisional power transfer shown in Figure 7, there is dominant electron heating for the lowest and highest 3He concentrations of 0.01% and 0.20%, and dominant bulk ion heating for 0.05% and 0.10%. The average energy of the fast 3He ions in the tail of the fast ion distribution function is higher in the case of the lowest concentration, but lower in every other case (in the range of 775–71.9 keV c.f. Table 4) than their critical energy (601–486 keV). Therefore, the 3He ions will collide mainly with the bulk ions and transfer their energy to them. It should be noted that the collisional transfer-related values in Table 4 do not seem to follow the trend seen in the H minority scheme at 2.65 T studied above in Section 3.1, where a high E fast value and a low E c r i t value result in dominant electron heating, and as these values converge, bulk ion heating becomes more relevant, according to [3]. The reason for this might be that the collisional power deposition profiles are broader in this case with a HFS resonance and a lower SPA. Therefore, the volume-integrated collisional powers P ci and P ce do not follow the scaling expected from the local E fast and E c r i t at the resonance, as given in Table 4. If we consider the local values of p ci and p ce at the resonance s = 0.62 in Figure 7, it can be seen that, apart from the lowest minority concentration case, the p ce and p ci profiles seem to meet close to the resonance, resulting in a relatively equivalent local power partitioning between bulk ion and electrons, which is consistent with E fast and E c r i t values in Table 4.
Figure 8 shows the time evolution of the thermal ion temperature at different 3He concentrations. In this scenario, the temperatures range from 5.78 keV to 5.92 keV at t f = 505 s, with the lowest T i = 5.78 keV at a 3He concentration of 0.01% and the highest T i = 5.92 keV at a 3He concentration of 0.20%. As shown in Figure 8, the ion temperature increases as the 3He concentration increases in the range of 0.01 to 0.20%. It should be noted that the transient evolution in the first seconds of the simulation is quite similar to the evolution in the H minority scheme at 2.65 T and also depends on the parameters used in the transport model. The highest 3He concentration yields the largest ICRF absorption, SPA and, consequently, the highest final ion temperature. It is also worth noting that the profiles for all concentrations have a lower temperature at the end of the simulation than at the beginning. This can be expected since the plasma was initially heated by a higher ECRH power than the ICRF power that was then applied during the simulation. Even though a percentage of the ICRF power has been absorbed by the resonant ions, it has not been enough to increase the ionic temperature.
Previously, the three-ion ICRF scheme 3He-H-4He at 3.3 T in H plasma was modeled [4] using TORIC [51], with a lower electron density of 3.4 × 10 19 m 3 and fixed plasma parameters. Effective damping was noted for 4He (10–15%) and 3He (0.05%). Competing absorption mechanisms included 3He (absorbing 70% ICRF power), direct electron damping (30%), and minor absorption by H ions at the plasma edge (not simulated in PION+ETS here). A combined SPA of 0.8–0.9 for direct electron heating and off-axis 3He ion heating was quoted [4]. Our time-evolving ETS+PION results at n e = 4.9 × 10 19 m 3 show comparable, albeit somewhat smaller, 3He absorption of ∼45% of ICRF power at a 3He concentration of 0.05% (c.f. Table 4), and stronger 3He absorption of 65% of ICRF power (SPA = 89%) at a higher 3He concentration of 0.20% (c.f. Figure 4).

4. Conclusions

In this paper, we presented the results of predictive simulations of two ICRF heating schemes in the ITER non-activated phase using for the first time integrated modeling with the ICRF modeling code PION integrated into the European Transport Solver (ETS) workflow within the ITER Analysis and Modeling Suite (IMAS). For the H minority scenario in the 4He plasma at 2.65 T, H damping was found to be a strong ICRF power absorption mechanism compared to the competing absorption mechanisms. The H power absorption and SPA increased with the H concentration and reached a maximum of ∼80% and a single pass absorption coefficient of SPA = 99 % at an H concentration of 10%, which was the highest value of the H concentration considered in this study. Electron heating dominated over bulk ion heating for every concentration except for an H concentration of 10%, where the collisional power equipartition was balanced. For the three-ion scheme 3He-H-4He at 3.3 T, the 3He resonance was the main absorption mechanism for concentrations above 0.05%, absorbing up to 65% at a 3He concentration of 0.20%. However, SPA values were generally lower compared to H minority heating, with a maximum SPA = 89 % at the 3He concentration of 0.20%, and the increase in ionic temperature for the concentration scan was almost negligible (∼0.2 keV). Overall, our results suggest that the plasma response to changes in ICRF heating, e.g., via changes in minority ion concentration as studied here, can be highly non-linear. Self-consistent simulations of ICRF and plasma transport physics are required in order to optimize the resulting plasma performance. There are plans for more comparisons between codes to keep validating the results, such as between PION as an H&CD actor within ETS and the ICRF heating code StixReDist [45] within ETS. Further work is in progress and will be reported elsewhere to further exploit the capabilities of the workflow integration and study the effect of ICRF heating on the plasma performance in ITER.

Author Contributions

Conceptualization, T.B., M.J.M. and D.G.; methodology, T.B. and M.J.M.; software, X.S. and J.M.; validation, T.B., J.M. and T.J.; formal analysis, T.B.; investigation, T.B. and T.J.; resources, T.B.; data curation, T.B.; writing—original draft preparation, T.B.; writing—review and editing, M.J.M. and D.G.; visualization, T.B. and T.J.; supervision, M.J.M.; project administration, M.J.M.; funding acquisition, M.J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union via the Euratom Research and Training Programme (Grant Agreement No 101052200 — EUROfusion).

Data Availability Statement

The data generated and/or analysed during the current study are not publicly available for legal/ethical reasons but are available from the corresponding author on reasonable request.

Acknowledgments

This work was carried out within the framework of the EUROfusion Consortium. Views and opinions expressed are, however, those of the author(s) only and do not necessarily reflect those of the European Union or the European Commission. Neither the European Union nor the European Commission can be held responsible for them.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. ITER Organization. ITER Technical Report ITR-18-003; ITER Organization: Saint-Paul-lez-Durance, France, 2018. [Google Scholar]
  2. Coblentz, L. ITER Council Meeting: ITER Project Making Progress, Preparing Updated Baseline; ITER Organization: Saint-Paul-lez-Durance, France, 2023. [Google Scholar]
  3. Stix, T. Waves in Plasmas; Springer: New York, NY, USA, 1992. [Google Scholar]
  4. Schneider, M.; Artaud, J.-F.; Bonoli, P.; Kazakov, Y.; Lamalle, P.; Lerche, E.; Van Eester, D.; Wright, J. ICRF heating schemes for the ITER non-active phase. EPJ Web Conf. 2017, 157, 03046. [Google Scholar] [CrossRef]
  5. Kazakov, Y.O.; Van Eester, D.; Dumont, R.; Ongena, J. On resonant ICRF absorption in three-ion component plasmas: A new promising tool for fast ion generation. Nucl. Fusion 2015, 55, 032001. [Google Scholar] [CrossRef]
  6. Eriksson, L.-G.; Hellsten, T.; Willen, U. Comparison of time dependent simulations with experiments in ion cyclotron heated plasmas. Nucl. Fusion 1993, 33, 1037. [Google Scholar] [CrossRef]
  7. Coster, D.P.; Basiuk, V.; Pereverzev, G.; Kalupin, D.; Zagórski, R.; Stankiewicz, R.; Huynh, P.; Imbeaux, F. The European Transport Solver. IEEE Trans. Plasma Sci. 2010, 38, 2085–2092. [Google Scholar] [CrossRef]
  8. Ryter, F.; Orte, L.B.; Kurzan, B.; McDermott, R.M.; Tardini, G.; Viezzer, E.; Bernert, M.; Fischer, R.; The ASDEX Upgrade Team. Experimental evidence for the key role of the ion heat channel in the physics of the L–H transition. Nucl. Fusion 2014, 54, 083003. [Google Scholar] [CrossRef]
  9. Eriksson, L.G.; Hellsten, T. A model for calculating ICRH power deposition and velocity distribution. Phys. Scr. 1995, 52, 70. [Google Scholar] [CrossRef]
  10. Hellsten, T.; Eriksson, L.G. A modelling scheme for the direct electron heating profiles during ion cyclotron resonance heating. Nucl. Fusion 1989, 29, 2165. [Google Scholar] [CrossRef]
  11. Start, D.F.H.; Jacquinot, J.; Bergeaud, V.; Bhatnagar, V.P.; Cottrell, G.A.; Clement, S.; Eriksson, L.-G.; Fasoli, A.; Gondhalekar, A.; Gormezano, C.; et al. D-T Fusion with Ion Cyclotron Resonance Heating in the JET Tokamak. Phys. Rev. Lett. 1998, 80, 4681. [Google Scholar] [CrossRef]
  12. Start, D.F.H.; Jacquinot, J.; Bergeaud, V.; Bhatnagar, V.P.; Conroy, S.W.; Cottrell, G.A.; Clement, S.; Ericsson, G.; Eriksson, L.-G.; Fasoli, A.; et al. Bulk ion heating with ICRH in JET DT plasmas. Nucl. Fusion 1999, 39, 321. [Google Scholar] [CrossRef]
  13. Mantsinen, M.J.; Eriksson, L.-G.; Bhatnagar, V.P.; Cottrell, G.A.; Gondhalekar, A.; Gormezano, C.; König, R.; Lomas, P.; Righi, E.; Rimini, F.G.; et al. Analysis of bulk ion heating with ICRH in JET high-performance plasmas. Plasma Phys. Control. Fusion 1999, 41, 843. [Google Scholar] [CrossRef]
  14. Mantsinen, M.J.; Eriksson, L.-G.; Gondhalekar, A.; Hellsten, T. Evidence for regions of nearly suppressed velocity space diffusion caused by finite Larmor radius effects during ICRF heating. Nucl. Fusion 1999, 39, 459. [Google Scholar] [CrossRef]
  15. Mantsinen, M.J.; Jarvis, O.N.; Kiptily, V.G.; Sharapov, S.E.; Alper, B.; Eriksson, L.-G.; Gondhalekar, A.; Heeter, R.F.; McDonald, D.C. First observation of pT fusion in JET tritium plasmas with ICRF heating of protons. Nucl. Fusion 2001, 41, 1815. [Google Scholar] [CrossRef]
  16. Eriksson, L.-G.; Mantsinen, M.J.; Rimini, F.G.; Nguyen, F.; Gormezano, C.; Start, D.F.H.; Gondhalekar, A. ICRF heating of JET plasmas with the third harmonic deuterium resonance. Nucl. Fusion 1998, 38, 265. [Google Scholar] [CrossRef]
  17. Eriksson, L.-G.; Mantsinen, M.J.; Bhatnagar, V.P.; Gondhalekar, A.; Gormezano, C.; Harbour, P.J.; Hellsten, T.; Jacquinot, J.; Jäckel, H.J.; Lawson, K.; et al. Theoretical analysis of ICRF heating in JET DT plasmas. Nucl. Fusion 1999, 39, 337–352. [Google Scholar] [CrossRef]
  18. Sharapov, S.E.; Borba, D.; Fasoli, A.; Kerner, W.; Eriksson, L.-G.; Heeter, R.F.; Huysmans, G.T.A.; Mantsinen, M.J. Stability of alpha particle driven Alfvén eigenmodes in high performance JET DT plasmas. Nucl. Fusion 1999, 39, 373. [Google Scholar] [CrossRef]
  19. Mantsinen, M.J.; Mayoral, M.-L.; Kiptily, V.G.; Sharapov, S.E.; Alper, B.; Bickley, A.; de Baar, M.; Eriksson, L.-G.; Gondhalekar, A.; Hellsten, T.; et al. Alpha-Tail Production with Ion-Cyclotron-Resonance Heating of 4He-Beam Ions in JET Plasmas. Phys. Rev. Lett. 2002, 88, 105002. [Google Scholar] [CrossRef]
  20. Salmi, A.; Mantsinen, M.J.; Beaumont, P.; de Vries, P.; Eriksson, L.-G.; Gowers, C.; er Hel, P.; Laxåback, M.; Noterdaeme, J.-M.; Testa, D.; et al. JET experiments to assess the clamping of the fast ion energy distribution during ICRF heating due to finite Larmor radius effects. Plasma Phys. Control. Fusion 2006, 48, 717. [Google Scholar] [CrossRef]
  21. Mantsinen, M.J.; Eriksson, J.; Challis, C.; Frigione, D.; Garcia, J.; Giroud, C.; Hjalmarsson, A.; King, D.B.; JET Contributors. ICRH analysis of high-performance JET hybrid discharges using PION modelling and neutron spectrometry measurements. In Proceedings of the 42nd EPS Conference on Plasma Physics, Lisbon, Portugal, 22–26 June 2015; p. P2.171. [Google Scholar]
  22. Mantsinen, M.; Challis, C.; Frigione, D.; Graves, J.; Hobirk, J.; Belonohy, E.; Czarnecka, A.; Eriksson, J.; Gallart, D.; Goniche, M.; et al. The Role of Combined ICRF and NBI Heating in JET Hybrid Plasmas in Quest for High D-T Fusion Yield. EPJ Web Conf. 2017, 157, 03032. [Google Scholar] [CrossRef]
  23. Gallart, D.; Mantsinen, M.; Challis, C.; Frigione, D.; Graves, J.; Hobirk, J.; Belonohy, E.; Czarnecka, A.; Eriksson, J.; Goniche, M.; et al. Modelling of combined ICRF and NBI heating in JET hybrid plasmas. EPJ Web Conf. 2017, 157, 03015. [Google Scholar] [CrossRef]
  24. Garcia, J.; Challis, C.; Gallart, D.; Garzotti, L.; Görler, T.; King, D.; Mantsinen, M.; JET Contributors. Challenges in the extrapolation from DD to DT plasmas: Experimental analysis and theory based predictions for JET-DT. Plasma Phys. Control. Fusion 2017, 59, 014023. [Google Scholar] [CrossRef]
  25. Mantsinen, M.J.; Gallart, D.; Belonohy, E.; Challis, C.; Czarnecka, A.; Eriksson, J.; Frigione, D.; Graves, J.; Goniche, M.; Hellesen, C.; et al. Optimising the use of ICRF waves in JET hybrid plasmas for high fusion yield. In Proceedings of the 44th EPS Conference on Plasma Physics, Belfast, UK, 26–30 June 2017; p. O3.110. [Google Scholar]
  26. Gallart, D.; Mantsinen, M.J.; Challis, C.; Frigione, D.; Graves, J.; Belonohy, E.; Casson, F.; Czarnecka, A.; Eriksson, J.; Garcia, J.; et al. Modelling of JET hybrid plasmas with emphasis on performance of combined ICRF and NBI heating. Nucl. Fusion 2018, 58, 106037. [Google Scholar] [CrossRef]
  27. Hellesen, C.; Mantsinen, M.; Conroy, S.; Ericsson, G.; Eriksson, J.; Kiptily, V.G.; Nabais, F.; JET Contributors. Analysis of resonant fast ion distributions during combined ICRF and NBI heating with transients using neutron emission spectroscopy. Nucl. Fusion 2018, 58, 056021. [Google Scholar] [CrossRef]
  28. Mantsinen, M.J.; Kazakov, Y.O.; Bobkov, V.; Craciunescu, T.; Gallart, D.; Giacomelli, L.; Kappatou, A.; Nocente, M.; Weil, M. Modelling of three-ion ICRF schemes with PION. In Proceedings of the 46th EPS Conference on Plasma Physics, Milan, Italy, 8–12 July 2019; p. O5.102. [Google Scholar]
  29. Kazakov, Y.O.; Nocente, M.; Mantsinen, M.J.; Ongena, J.; Baranov, Y.; Craciunescu, T.; Dreval, M.; Dumont, R.; Eriksson, J. Plasma heating and generation of energetic D ions with the 3-ion ICRF + NBI scenario in mixed H-D plasmas at JET-ILW. Nucl. Fusion 2020, 60, 112013. [Google Scholar] [CrossRef]
  30. Taylor, D.M.A.; Mantsinen, M.J.; Gallart, D.; Manyer, J.; Sirén, P.; JET Contributors. Effect of inclusion of pitch-angle dependence on a simplified model of RF deposition in tokamak plasma. Plasma Phys. Control. Fusion 2022, 64, 055015. [Google Scholar] [CrossRef]
  31. Gallart, D.; Mantsinen, M.J.; Manyer, J.; Planas, E.; Taylor, D.M.A.; Garcia, J.; Frigione, D.; Garzotti, L.; Kim, H.-T.; Nocente, M.; et al. Prediction of ICRF minority heating schemes for JET D–T experiments. Plasma Phys. Control. Fusion 2022, 64, 125006. [Google Scholar] [CrossRef]
  32. García-Muñoz, M.; Fahrbach, H.-U.; Günter, S.; Igochine, V.; Mantsinen, M.J.; Maraschek, M.; Martin, P.; Piovesan, P.; Sassenberg, K.; Zohm, H. Fast-Ion Losses due to High-Frequency MHD Perturbations in the ASDEX Upgrade Tokamak. Phys. Rev. Lett. 2008, 100, 055005. [Google Scholar] [CrossRef]
  33. Mantsinen, M.J.; Bilato, R.; Bobkov, V.V.; Kappatou, A.; McDermott, R.M.; Nocente, M.; Odstrčil, T.; Tardini, G.; Bernert, M.; Dux, R.; et al. Bulk ion heating with ICRF waves in tokamaks. AIP Conf. Proc. 2015, 1689, 030005. [Google Scholar] [CrossRef]
  34. Mantsinen, M.J.; Bobkov, V.; Gallart, D.; Geiger, B.; Johnson, T.; Meyer, H.; Nocente, M.; Ochoukov, R.; Odstrcil, T.; Perelli, E.; et al. Third harmonic ICRF heating of deuterium beam ions on ASDEX Upgrade. In Proceedings of the 43rd EPS Conference on Plasma Physics, Leuven, Belgium, 4–8 July 2016; p. P1.035. [Google Scholar]
  35. Sharapov, S.E.; Garcia-Munoz, M.; Van Zeel, M.A.; Bobkov, B.; Classen, I.G.J.; Ferreira, J.; Figueiredo, A.; Fitzgerald, M.; Galdon-Quiroga, J.; Gallart, D.; et al. The effects of electron cyclotron heating and current drive on toroidal Alfvén eigenmodes in tokamak plasmas. Plasma Phys. Control. Fusion 2017, 60, 014026. [Google Scholar] [CrossRef]
  36. Galdon-Quiroga, J.; Garcia-Munoz, M.; Sanchis-Sanchez, L.; Mantsinen, M.; Fietz, S.; Igochine, V.; Maraschek, M.; Rodriguez-Ramos, M.; Sieglin, B.; Snicker, A.; et al. Velocity space resolved absolute measurement of fast ion losses induced by a tearing mode in the ASDEX Upgrade tokamak. Nucl. Fusion 2018, 58, 036005. [Google Scholar] [CrossRef]
  37. Mantsinen, M.J.; Bobkov, V.; Gallart, D.; Pütterich, T.; Sauter, O.; the EUROfusion MST1 Team; the ASDEX Upgrade Team. Modelling of ICRF heating in ASDEX Upgrade discharges with pure wave heating relevant to the ITER baseline scenario. In Proceedings of the 45th EPS Conference on Plasma Physics, Prague, Czech Republic, 2–6 July 2018; p. P1.1072. [Google Scholar]
  38. Manyer, J.; Mantsinen, M.; Gallart, D.; Florido-Llinas, M.; Taylor, D.; Bobkov, V.; Pütterich, T.; Sauter, O.; Weil, M.; the EUROfusion MST1 Team; et al. Modelling of dual-frequency ICRF heating in ASDEX Upgrade discharges relevant to the ITER baseline scenario. In Proceedings of the 47th EPS Conference on Plasma Physics, Virtual, 21–25 June 2021; p. P.5.1055. [Google Scholar]
  39. Mantsinen, M.J.; Petty, C.C.; Eriksson, L.-G.; Mau, T.K.; Pinsker, R.I.; Porkolab, M. Analysis of combined fast wave current drive and neutral beam injection in the DIII-D tokamak. Phys. Plasmas 2002, 9, 1318. [Google Scholar] [CrossRef]
  40. Eriksson, L.-G.; Hoang, G.T.; Bergeaud, V. On the role of ion heating in ICRF heated discharges in Tore Supra. Nucl. Fusion 2001, 41, 91. [Google Scholar] [CrossRef]
  41. Becoulet, A.; the Equipe Tore Supra. Heating and Current Drive System for Tore Supra Steady-State Operation. J. Plasma Fusion Res. 2000, 3, 51. [Google Scholar]
  42. Arbina, I.L.; Mantsinen, M.J.; Sáez, X.; Gallart, D.; Gutiérrez, A.; Taylor, D.; Johnson, T.; Pinches, S.D.; Schneider, M.; the EUROfusion-IM Team. First applications of the ICRF modelling code PION in the ITER Integrated Modelling and Analysis Suite. In Proceedings of the 46th EPS Conference on Plasma Physics, Milan, Italy, 8–12 July 2019; p. P4.1079. [Google Scholar]
  43. Imbeaux, F.; Pinches, S.D.; Lister, J.B.; Burav, Y.; Casper, T.; Duval, B.; Guillerminet, B.; Hosokawa, M.; Houlberg, W.; Huynh, P.; et al. Design and first applications of the ITER integrated modelling & analysis suite. Nucl. Fusion 2015, 55, 123006. [Google Scholar] [CrossRef]
  44. Lamalle, P.U. LPP-ERM/KMS Report 101. Ph.D. Thesis, Université de Mons, Mons, Belgium, 1994. [Google Scholar]
  45. Van Eester, D.; Lerche, E. Simple 1D Fokker-Planck modelling of ion cyclotron resonance frequency heating at arbitrary cyclotron harmonics accounting for Coulomb relaxation on non-Maxwellian populations. Plasma Phys. Control. Fusion 2011, 53, 9. [Google Scholar] [CrossRef]
  46. Kalupin, D. Construction of the European transport solver under the European Integrated Tokamak Modelling Task Force. In Proceedings of the 35th EPS Conference on Plasma Physics, Hersonissos, Greece, 9–13 June 2008; Volume 32D, p. P-.5.027. [Google Scholar]
  47. Bilato, R.; Coster, D.; Dumont, R.; Johnson, T.; Klingshirn, H.-J.; Lerche, E.; Sauter, O.; Brambilla, M.; Figini, L.; Van Eester, D.; et al. ICRF-code benchmark activity in the framework of the European task-force on integrated Tokamak Modelling. AIP Conf. Proc. 2014, 1580, 291–294. [Google Scholar] [CrossRef]
  48. Koch, R.; Lerche, E.; Van Eester, D.; Nightingale, M. High frequency fast wave current drive for DEMO. AIP Conf. Proc. 2011, 1406, 349–352. [Google Scholar] [CrossRef]
  49. Van Eester, D.; Koch, R. A variational principle for studying fast-wave mode conversion. Plasma Phys. Control. Fusion 1998, 40, 11. [Google Scholar] [CrossRef]
  50. Porkolab, M.; Bécoulet, A.; Bonoli, P.T.; Gormezano, C.; Koch, R.; Majeski, R.J.; Messiaen, A.; Noterdaeme, J.M.; Petty, C.; Pinsker, R.; et al. Recent progress in ICRF physics. Plasma Phys. Control. Fusion 1998, 40, A35. [Google Scholar] [CrossRef]
  51. Brambilla, M. Numerical simulation of ion cyclotron waves in tokamak plasmas. Plasma Phys. Control. Fusion 1999, 41, 1. [Google Scholar] [CrossRef]
Figure 1. Overview of ITER discharge 110005 for H concentrations of [H] = 1.0, 2.5, 5.0 and 10%. (a) Temperature of the ions at the start of the simulation ( T i 0 ), (b) temperature of the ions at the end of the simulation ( T i ), (c) density of the H ions at the start of the simulation ( n H 0 ), (d) temperature of the electrons at the start of the simulation ( T e 0 ), (e) temperature of the electrons at the end of the simulation ( T e ), (f) density of the electrons at the start of the simulation n e 0 .
Figure 1. Overview of ITER discharge 110005 for H concentrations of [H] = 1.0, 2.5, 5.0 and 10%. (a) Temperature of the ions at the start of the simulation ( T i 0 ), (b) temperature of the ions at the end of the simulation ( T i ), (c) density of the H ions at the start of the simulation ( n H 0 ), (d) temperature of the electrons at the start of the simulation ( T e 0 ), (e) temperature of the electrons at the end of the simulation ( T e ), (f) density of the electrons at the start of the simulation n e 0 .
Plasma 07 00028 g001
Figure 2. ICRF power density absorbed by fundamental H (red), 2nd 4He harmonic (green), and direct electron damping (blue) as a function of the normalized flux surface, s, at t f = 308 s for H concentrations of 1.0, 2.5, 5.0 and 10%. SPA coefficients on H at the resonance and toroidal mode number N = 54 are given in the legends.
Figure 2. ICRF power density absorbed by fundamental H (red), 2nd 4He harmonic (green), and direct electron damping (blue) as a function of the normalized flux surface, s, at t f = 308 s for H concentrations of 1.0, 2.5, 5.0 and 10%. SPA coefficients on H at the resonance and toroidal mode number N = 54 are given in the legends.
Plasma 07 00028 g002
Figure 3. Orbit redistributed collisional power transferred by the resonant ion species to bulk ions (red) and electrons (blue) at t f = 308 s.
Figure 3. Orbit redistributed collisional power transferred by the resonant ion species to bulk ions (red) and electrons (blue) at t f = 308 s.
Plasma 07 00028 g003
Figure 4. Time evolution of the thermal ions temperature at the normalized flux surface s = 0.03 where the ICRF resonance is located, from t i = 300 s to t f = 308 s, for H concentrations of 1.0, 2.5, 5.0 and 10%.
Figure 4. Time evolution of the thermal ions temperature at the normalized flux surface s = 0.03 where the ICRF resonance is located, from t i = 300 s to t f = 308 s, for H concentrations of 1.0, 2.5, 5.0 and 10%.
Plasma 07 00028 g004
Figure 5. Overview of ITER discharge 104010 for 3He concentrations of [3He] = 0.01, 0.05, 0.10 and 0.20%. (a) Temperature of the ions at the start of the simulation ( T i 0 ), (b) temperature of the ions at the end of the simulation ( T i ), (c) density of the 3He ions at the start of the simulation ( n 3 H e 0 ), (d) temperature of the electrons at the start of the simulation ( T e 0 ), (e) temperature of the electrons at the end of the simulation ( T e ), (f) density of the electrons at the start of the simulation n e 0 .
Figure 5. Overview of ITER discharge 104010 for 3He concentrations of [3He] = 0.01, 0.05, 0.10 and 0.20%. (a) Temperature of the ions at the start of the simulation ( T i 0 ), (b) temperature of the ions at the end of the simulation ( T i ), (c) density of the 3He ions at the start of the simulation ( n 3 H e 0 ), (d) temperature of the electrons at the start of the simulation ( T e 0 ), (e) temperature of the electrons at the end of the simulation ( T e ), (f) density of the electrons at the start of the simulation n e 0 .
Plasma 07 00028 g005
Figure 6. RF-power density absorbed by the fundamental 3He (green) and direct electron damping (blue) as a function of the normalized flux surface, s, at the end of the simulation at t f = 505 s, for [3He] = 0.01, 0.05, 0.10 and 0.20%. SPA coefficients on 3He at the resonance and toroidal mode number N = 54 are given in the legends.
Figure 6. RF-power density absorbed by the fundamental 3He (green) and direct electron damping (blue) as a function of the normalized flux surface, s, at the end of the simulation at t f = 505 s, for [3He] = 0.01, 0.05, 0.10 and 0.20%. SPA coefficients on 3He at the resonance and toroidal mode number N = 54 are given in the legends.
Plasma 07 00028 g006
Figure 7. Orbit redistributed collisional power transfer from the resonant ion species to bulk ions (red) and electrons (blue) at t f = 505 s, for [3He] = 0.01, 0.05, 0.10 and 0.20%.
Figure 7. Orbit redistributed collisional power transfer from the resonant ion species to bulk ions (red) and electrons (blue) at t f = 505 s, for [3He] = 0.01, 0.05, 0.10 and 0.20%.
Plasma 07 00028 g007
Figure 8. Time evolution of the thermal ion temperature at the location of the 3He resonance at s = 0.62 from t i = 500 s to t f = 505 s for [3He] = 0.01, 0.05, 0.10 and 0.20%.
Figure 8. Time evolution of the thermal ion temperature at the location of the 3He resonance at s = 0.62 from t i = 500 s to t f = 505 s for [3He] = 0.01, 0.05, 0.10 and 0.20%.
Plasma 07 00028 g008
Table 1. Plasma composition (with the minority ion species in brackets), heating scheme, magnetic field ( B 0 ), ICRF frequency (f), ICRF power ( P I C R F ), central electron density ( n e 0 ) and temperature ( T e 0 ).
Table 1. Plasma composition (with the minority ion species in brackets), heating scheme, magnetic field ( B 0 ), ICRF frequency (f), ICRF power ( P I C R F ), central electron density ( n e 0 ) and temperature ( T e 0 ).
PlasmaHeating B 0 f P ICRF n e 0 × 10 19 T e 0
Scheme (T) (MHz) (MW) (m−3) (keV)
4He-(H) ω = ω H 2.6540203.310
Table 2. H concentrations, power absorbed by resonant H and 4He ions ( P a b s , H , P a b s , 4 H e ), fractions of collisional power transferred from resonant ions to bulk ions ( P c i / P c ) and background electrons ( P c e / P c ), critical energy ( E c r i t ) and average energy of the fast H ions ( E f a s t , H ) as given by PION+ETS.
Table 2. H concentrations, power absorbed by resonant H and 4He ions ( P a b s , H , P a b s , 4 H e ), fractions of collisional power transferred from resonant ions to bulk ions ( P c i / P c ) and background electrons ( P c e / P c ), critical energy ( E c r i t ) and average energy of the fast H ions ( E f a s t , H ) as given by PION+ETS.
H P abs , H P abs , 4 He P ci / P P ce / P E crit E fast , H
(%)(MW)(MW)(%)(%)(keV)(keV)
1.011.53.9234.265.81581719
2.515.20.7326.873.2153568
5.015.60.4039.960.1154366
1015.70.2548.852.2173459
Table 3. Plasma composition (with the minority ion species in brackets), heating scheme, magnetic field ( B 0 ), ICRF frequency (f), ICRF power ( P I C R F ), central electron density ( n e 0 ) and temperature ( T e 0 ).
Table 3. Plasma composition (with the minority ion species in brackets), heating scheme, magnetic field ( B 0 ), ICRF frequency (f), ICRF power ( P I C R F ), central electron density ( n e 0 ) and temperature ( T e 0 ).
PlasmaHeating B 0 f P ICRF n e 0 × 10 19 . T e 0
Scheme (T) (MHz) (MW) (m−3) (keV)
H-4He-(3He) ω = ω 3 H e 3.340204.912
Table 4. 3He concentrations, power absorbed by resonant 3He ions and electrons ( P a b s , 3 H e , P a b s , e ), fractions of collisional power transferred from resonant ions to bulk ions ( P c i / P c ) and background electrons ( P c e / P c ), critical energy ( E c r i t ) and average energy of the fast 3He ions ( E f a s t , 3 H e ) as given by PION+ETS.
Table 4. 3He concentrations, power absorbed by resonant 3He ions and electrons ( P a b s , 3 H e , P a b s , e ), fractions of collisional power transferred from resonant ions to bulk ions ( P c i / P c ) and background electrons ( P c e / P c ), critical energy ( E c r i t ) and average energy of the fast 3He ions ( E f a s t , 3 H e ) as given by PION+ETS.
3He P abs , 3 He P abs , e P ci / P c P ce / P c E c r i t E fast , 3 He
(%)(MW)(MW)(%)(%)(keV)(keV)
0.014.2614.017.182.9601775
0.058.569.3951.948.1510169
0.1010.57.1855.744.348671.9
0.2011.36.0843.856.2496102
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bensadon, T.; Mantsinen, M.J.; Jonsson, T.; Gallart, D.; Sáez, X.; Manyer, J. Analysis of ICRF Heating Schemes in ITER Non-Active Plasmas Using PION+ETS Integrated Modeling. Plasma 2024, 7, 517-530. https://doi.org/10.3390/plasma7030028

AMA Style

Bensadon T, Mantsinen MJ, Jonsson T, Gallart D, Sáez X, Manyer J. Analysis of ICRF Heating Schemes in ITER Non-Active Plasmas Using PION+ETS Integrated Modeling. Plasma. 2024; 7(3):517-530. https://doi.org/10.3390/plasma7030028

Chicago/Turabian Style

Bensadon, Tomas, Mervi J. Mantsinen, Thomas Jonsson, Dani Gallart, Xavier Sáez, and Jordi Manyer. 2024. "Analysis of ICRF Heating Schemes in ITER Non-Active Plasmas Using PION+ETS Integrated Modeling" Plasma 7, no. 3: 517-530. https://doi.org/10.3390/plasma7030028

APA Style

Bensadon, T., Mantsinen, M. J., Jonsson, T., Gallart, D., Sáez, X., & Manyer, J. (2024). Analysis of ICRF Heating Schemes in ITER Non-Active Plasmas Using PION+ETS Integrated Modeling. Plasma, 7(3), 517-530. https://doi.org/10.3390/plasma7030028

Article Metrics

Back to TopTop