Prerpints.org logo
Preprint
Article

Exploring Experimental Isotope Scaling and Density Limit in Tokamak Transport

Altmetrics

Downloads

80

Views

39

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

08 August 2024

Posted:

12 August 2024

You are already at the latest version

Alerts
Abstract
As it turns out, both isotope scaling and density limits are phenomena closely linked to fluid closure. The necessity to include ion viscosity arises for both phenomena. Thus, we have added ion viscosity to our model. The experimental isotope scaling has been successfully recovered in our fluid model through parameter scans. Although ion viscosity typically exerts a small effect, the density limit is manifested by increasing the density by approximately tenfold from the typical experimental density. In our case, this increase originates from the density in the Cyclone base case. Notably, these phenomena would not manifest with a gyro-Landau fluid closure. The isotope scaling is nullified by the addition of a gyro-Landau term, while the density limit results from permitting ion viscosity to become comparable to the gyro-Landau term. The mechanism of zonal flows, demonstrated analytically for the Dimits upshift, yields insights into the isotope scaling observed in experiments. In our approach, ion viscosity is introduced in place of the Landau fluid resonances found in some fluid models. This implies that the mechanism of isotope scaling operates at the level of fluid closure in connection with the generation of zonal flows. The strength of zonal flows in our model has been verified, particularly in connection with the successful simulation of the nonlinear Dimits shift. Consequently, a role is played by our approach in the temperature perturbation part of the Reynolds stress.
Keywords: 
Subject: Physical Sciences  -   Fluids and Plasmas Physics

1. Introduction

The general problem of tokamak transport has been one of the main issues in fusion research for a long time [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,51]. However, the specific mechanism behind isotope scaling remains undetermined [1,2,3,4]. Similarly, the underlying reason for the density limit has not been clearly understood [11,36,51]. Several other papers [38,39,40,41,42,43] address the density limit. It is commonly associated with increased turbulent transport and collisions, which is also true in our model. In some cases (Gates et al. Ref. [40]), the coupling to magnetic perturbations is emphasized. Of course, magnetic field perturbations increase with the turbulence level in our model as well, but we do not see a causal relationship here. However, our model for the L-H transition [28], without ion viscosity, is electromagnetic and agrees very well with the model in Ref. [31]. The density limit in Ref. [31] is associated with collisions (as in our model) and was also discussed by Giacomin et al. Since we obtain the density limit by including ion viscosity, we are consistent with Ref. [31]. In the paper by Giacomin et al., it is also mentioned that the pressure gradient can become a considerable fraction of the minor radius, as in the case of MARFEs.
In this study, we link both phenomena to fluid closure, a mechanism that must also be present in kinetic formulations [10], though it is not always prominently highlighted. One aspect that appears to be important is the effects of classical dissipation [1,3]. As shown numerically in Ref. [6] and analytically in Ref. [7], using the reductive perturbation method [8], the Landau fluid resonance [14] changes the nonlinear Dimits upshift [10] strongly. Our fluid model results, in agreement with the kinetic Dimits shift [6,7], show that we have the correct strength of zonal flows. This means that the inverse turbulent cascade [16] is damped out so that there is no pileup of waves at the system size. This also verifies the choice of absorbing boundary for long wavelengths, as used in Ref. [17]. It is interesting to compare it with the resistive drift wave fluid model [9], which in Ref. [18] gave a kernel transport coefficient of the same form as ours. With kernel, we mean that we keep only diagonal elements of the transport matrix. This requires an outgoing turbulence boundary without pileup at the system size. In this work, both analytical and numerical methods are used. However, ion viscosity is typically at least an order of magnitude smaller than the Landau fluid resonance [14] and can thus be incorporated into our general approach to fluid closure. As it turns out, usually edge data give a stronger isotope effect (about 20% difference in transport between hydrogen and tritium). Thus, we also have an interest in the L-H transition [28,31,32,33] which is also due to zonal flows. Actually, our fluid model did very well in comparison with local transport simulations as made in Ref. [31]. Although our simulations were originally set up for EAST [28], we were able to investigate the same parameter regimes used in Refs. [31,32] for Alcator C-Mod. This allowed us to construct a stability diagram in the magnetohydrodynamic (MHD) parameter alpha and a corresponding diamagnetic resistive parameter alpha, showing good agreement. As found in Refs. [18,29], a quasilinear approach is quite adequate for studying drift wave turbulence, and in Ref. [29], a fluid approach was found to be accurate. Our scaling of temperature at the separatrix with the magnetic field and the power threshold were also consistent with those reported in Ref. [33]. Thus, we have secured the validity of our simulation of the L-H transition [28,44]. All our new findings are related to the excitation of zonal flows and the associated fluid closure aspects. This connects to our previous work on fast particles [12], where resonance broadening was shown to be important for both fast particle instabilities and drift waves, eliminating wave-particle resonances in the latter case. The foundation of this research lies in our drift wave fluid model [17] (usually referred to as the Weiland model), which still serves as a valid limit of the current model. The motivation behind implementing fluid closure was to incorporate all moments with sources observed in the experiment. This transport model still has the same normalization as in Ref. [17]. It has never been fitted to any other theoretical model or experiment, yet it consistently yields good agreement with conventional (JET, EAST, KSTAR, and DIII-D) and low aspect ratio tokamak (NSTX) plasma profiles [45,46,47].
As is clear from this introduction, zonal flows play an important role in tokamak transport. We here mention three more papers in this field [20,21,22]. However, none of these include resonance broadening [23,24], which is now the main argument for our fluid closure [12]. Contrary to flattening, which means that particles move out of resonance with waves, resonance broadening means that waves move out of resonance with particles.
The manuscript is structured as follows: Section 2 discusses the inclusion of all moments with sources in the experimental closure approach and the impact of resonance broadening on the fluid model. It concludes that resonance broadening stabilizes fast particle instabilities and validates the reactive closure for drift waves, allowing the use of Braginskii’s highest moment for successful modeling. Section 3 outlines the original motivation behind a fluid closure, driven by the aim to encompass all moments with sources and justified by a nonlinear Fokker-Planck equation. It further discusses the implications of the Fokker-Planck equation in understanding turbulence and resonance broadening effects in wave-particle interactions. In Section 4, we explore the incorporation of ion viscosity, emphasize the pivotal role of zonal flows generated by Reynolds stress—specifically, in the form of off-diagonal poloidal momentum flux—and stress the critical consideration of averaged resonances for accurate predictions. Section 5 highlights the relevance of edge data for understanding the H-mode barrier, presenting ion thermal diffusivity against normalized temperature gradients. The observed isotope mass effect leads to a confinement time scaling in agreement with experimental results, along with a global scaling involving heating power. Section 6 explores the strong dependence on system size and the density limit, highlighting the impact of gyro-Landau resonances on rotation dynamics. The interplay with E × B convection is elucidated, emphasizing the need for strong zonal flows to avoid wave pileup and reflections, particularly in finite radial systems. The discussion encompasses stabilizing nonlinearity, showcasing its critical role in preventing an inverse cascade towards longer wavelengths. The importance of such considerations is underscored in reactive systems, as demonstrated in previous works. Section 7 delves into successful applications of the transport model, focusing on isotope scaling and the density limit, while acknowledging challenges in the latter. The discussion underlines the importance of ion temperature perturbation resonance in turbulence-rotation balance and tokamak transport dynamics.

2. Fluid Model

The original approach to our fluid closure was to include all moments with sources in the experiment [17], the reason being that moments without sources would die out on the transport time scale. An important mechanism turned out to be resonance broadening [19]. This was applied to our model, with the conclusion that we would get a reactive closure with a diagonal part of the same type as in our fluid model [23]. We found that turbulent collisions have the same effect on a fluid description as classical collisions. A generalization to include fast particles verified both the stabilizing trend of resonance broadening on fast particle instabilities [12] and our previous result that our closure is valid for drift waves. Since resonance broadening only shows that we have a reactive closure, we again employed our original rule for the exact point of closure [52]. This closure means that we can take the highest moment at the Braginskii q * , and this has been extremely successful in previous studies in cases where the ion viscosity could be ignored. Resonance broadening means that nonlinear frequency shifts change the phase velocity of waves in such a way that waves move out of resonance with particles when we average over turbulent fluctuations. Thus, dissipative wave-particle resonances like Landau damping and magnetic drift resonances are averaged out. Of course, this will be true for all processes, so we can also use the Braginskii derivation for ion viscosity.

3. Formulation

The original motivation for our fluid closure was to include all moments with sources in the experiment. The elimination of higher moments was then motivated by a nonlinear Fokker-Planck equation. The Fokker-Planck equation for turbulent collisions [12,23] can be written as
∂ ∂ t + v · ∂ ∂ r W ( X , X ′ , t , t ′ ) = ∂ ∂ v β v + C i + D ν ∂ ∂ v W ( X , X ′ , t , t ′ ) + S v .
β = ∑ j β j | ϕ j | 2
D ν = ∑ j d j | ϕ j | 2 .
We are here interested in drift waves, so we take S v = 〈 S v 〉 = 0 . In the given context, W represents the transition probability between states X and X ′ over the time interval t − t ′ . The terms β , C i , D ν , and ϕ denote the nonlinear turbulent friction, classical friction term due to ion-ion collisions, nonlinear diffusion coefficient in velocity space, and electrostatic potential, respectively. The Fokker-Planck equation leads to the mean square of velocity deviation (velocity dispersion) developing as [23]
〈 ▵ v 2 〉 = D ν β ( 1 − e − β t ) .
This quantity measures the average deviation of particle velocities from their mean due to turbulent collisions over time. It helps in understanding the intensity and impact of turbulence on plasma behavior and illustrates how particles spread over time. A simplified derivation of Equation 3 was given in Ref. [34] Chapter 9, where the Chandrasekhar solution [53] of Equation 1 without viscosity was used as a weight function. Now, the classical collision term just gives the usual classical viscosity, so we will recover Equation 4 since Equation 3 means that we have a reactive closure as discussed in Section 4.
The Fokker-Planck equation can also be applied to fast particles if we add a fast particle source, and then we can recover equations for both fast particles and drift waves within appropriate limits [12]. In both cases, resonance broadening [19,23] reduces the wave-particle interaction. As it turns out, the resonance broadening is due to nonlinear frequency shifts, which would also remain in the coherent limit [26,27].
A solution for coherent explosive instability with a nonlinear frequency shift (Figure 1) was shown already in our book [24]. Similar to our recent work [48], we can explore the coherent three-wave interaction regime as also studied in Ref. [13] and Ref. [54] for drift waves [26,27].
The coherent, unstable three-wave system is stabilized by nonlinear frequency shifts. These work so as if the sign of the wave energy is shifted, thus we have alternatively Landau growth and damping. This is the way the effect of waves-particles [23] interaction is averaged out, leading to the absence of energy transfer between waves and particles, as shown by Equation 3 in the turbulent case.
In summary, this section motivates the use of a fluid closure approach, which utilizes a nonlinear Fokker-Planck equation to include all experimental moments and eliminate higher moments. The equation accounts for nonlinear turbulent friction, classical collision friction, and nonlinear diffusion, and can be extended to fast particles and drift waves, leading to resonance broadening and balanced wave-particle interactions. This approach provides a comprehensive framework for understanding turbulent collision dynamics and wave-particle interactions. By incorporating nonlinear effects, it offers a more accurate and robust description of plasma behavior, enabling better predictions and insights in various plasma physics applications.

4. Exploring Ion Viscosity Effects in Drift Wave Turbulence

In order to study isotope scaling, we add ion viscosity to our usual derivation.
3 2 n i ∂ ∂ t + v i · ∇ T i + P i ∇ · v i = − ∇ · q * i + i ν ii ,
As discussed after Equation 3, the model can use ion viscosity from Braginskii’s equations directly. Turbulent viscosity behaves like nonlinear friction (see Ref. [23]) and can be included without changing the model’s final results, where ion viscosity
ν ii = ν ee T e T i 3 / 2 ( m e m i ) 0.5 1 A 0.5 ,
where A is the isotope mass number and ν ee is the electron electron collision frequency. Our viscosity, however, has been taken from Ref. [13], Equation (1.8). This gives the scaling T i − 3 / 2 n / A 0 . 5 for the viscosity. Thus, we see that viscosity is smaller for heavier isotopes (favorable scaling) and larger for higher density (unfavorable scaling). Equation 5 leads to the ion temperature perturbation
δ T i T i = ω ω − 5 ω Di / 3 + i ν ii 2 3 δ n i n + ω * e ω ( η i − 2 3 ) e ϕ T e .
Note that the viscosity term enters at the fluid resonance in the ion energy equation, the most sensitive point in our fluid modeling. This sensitivity is the main improvement in our model. Here n i is the ion density, v i is the ion flow velocity, q * i is the diamagnetic ion heat flow, T e ( T i ) is the electron (ion) temperature, e is the electron charge, m i is the ion mass, m e is the electron mass, ϵ 0 is the permittivity of free space, ω = ω r + i γ , where ω r is the real frequency and γ is the mode’s growth rate, ω Di is the magnetic drift frequency, ω * e is the electron diamagnetic drift frequency, and η i is the ratio of ion temperature gradient to ion density gradient. It is important to recall that in our usual reactive limit, the linear threshold of Ion Temperature Gradient (ITG) modes appears exactly at the resonance in Equation 6. As it turns out, zonal flows play important roles in both Dimits shifts, the L-H transition, isotope scaling, and finally the density limit. Zonal flows are generated by the Reynolds stress.
We note that Equations (4), (5),(6) are a generalization of our fluid closure, where the principle is to omit dissipative kinetic resonances, keeping only moments with sources in the experiment. This can be motivated by resonance broadening, as explained in Ref. [12]. Resonance broadening and profile flattening are parallel phenomena that both remove wave particle resonances, although resonance broadening is strongly nonlinear. If you have a particle source that injects particles into the resonant region in phase space, it can balance the flattening so that the resonance survives. The same thing can happen with resonance broadening. Here we have the source of fast particles, which can balance the resonance broadening so that the resonance remains. In this case, it is the question of waves moving out of resonance with particles. However, as we clarify in our paper on fast particles [12], drift waves lack a source because their frequency is approximately two orders of magnitude lower than the source’s, leading to the averaging of the source. Then only resonance broadening remains, and the wave particle resonance is cancelled. We have shown this in Ref. [23], and our paper on fast particles [12] just uses this.
Γ p = 〈 V E r V θ 〉 = − 1 2 D B 2 k r k θ ϕ ^ * ϕ ^ + 1 τ P ^ i + c . c ,
In the expression, Γ p represents the off-diagonal poloidal momentum flux, V E r denotes the radial component of the E × B drift, V θ signifies the poloidal flow velocity, D B is defined as ρ s c s where c s represents the sound speed, k θ stands for the poloidal wavenumber, ϕ ^ is defined as e ϕ / T e , and normalized ion pressure, P ^ i is calculated as δ P i / P i . In Equation 7, it is noteworthy to observe the dependence on ion temperature through ion pressure ( P i ). This dependence is highly sensitive to the fluid closure, as indicated by Equation 6. Here, the last part, the E × B convection of the diamagnetic flow, is often ignored but is, in fact, the most important part.
In the absence of viscosity, marginal stability for the reactive system enters at the resonance in Equation 6. Because of the resonance, we get a particularly strong drive of rotation here (see Figure 2). Although ion viscosity is usually small, it can still be of importance here, and it turns out to generate the isotope scaling. We recall from Ref. [7] that the resonance in Equation 6 is completely smeared out by gyro-Landau fluid resonances, giving a significantly smaller rotation. We note that gyro-Landau resonances are averaged out by resonance broadening and are, without this averaging, typically an order of magnitude larger than ion viscosity but added in the same place. We may here add that nonlinear flattening means that particles move out of resonance with waves, while resonance broadening means that waves are moving out of resonance with particles. Thus, there is no way that we could obtain the isotope scaling if we kept the unaveraged gyro-Landau resonances.
In summary, this section examines the impact of ion viscosity on drift wave turbulence and isotope scaling, utilizing Braginskii’s formulation for ion viscosity and noting that turbulent viscosity behaves similarly to nonlinear friction. The derived equations illustrate the effect of ion viscosity on ion temperature perturbations and validate the fluid closure approach. The study highlights the importance of zonal flows and resonance broadening in achieving isotope scaling and maintaining plasma stability.

5. Isotope Effects on Transport and Confinement: Edge Considerations and Predictive Comparison

As pointed out above, edge data are relevant since the transport flux has to pass the edge. Since the H-mode barrier depends on the same type of zonal flows as the Dimits shift, ion thermal diffusivity as a function of the normalized temperature gradient is shown for the parameters relevant to the plasma edge in Figure 3. Figure 3 directly continues from Figure 2. We are extending the results of Ref. [7] by including ion viscosity.
It is known that the pedestal can be close to force balance, but the L-H transition process involves the creation of zonal flows via poloidal rotation. This allows the use of analytic results from Ref. [7] to determine transport calculations. When switching from the use of hydrogen to tritium, a 6% drop in transport is noted in typical edge data.
Using data from our local code and from Figure 3, the scaling leads to the confinement time ( τ E ) scaling with isotope mass (A),
τ E = A 0.2 .
This result agrees with the experimental result in Ref. [36]. We also recall our result for global scaling with heating power (P) [37]. Additionally, it is important to note the strong FLR stabilization for large gradients (see, e.g., Ref. [34], Equation 6.161).
In summary, these results may be stated as follows:
τ E = P − 0 . 67 A 0 . 2 .
Now, using our predictive code for ITER simulations of hydrogen and tritium, we instead obtain:
τ E = A 0 . 5 .
This result was obtained by comparing runs with different isotopes at different temperatures assuming that τ E scales as T − 2 / 3 . The ITER simulations, focusing on ELMY H-mode discharges and employing the Weiland [17,49] and NCLASS [50] neoclassical models for predictive simulation, are shown in Figure 4. The simulations initiate with prescribed sources and an assumed L-mode profile, progressing to the temperature of the L–H transition and pedestal. We do not make any assumptions about the occurrence of an L–H transition or the specific locations of barriers to temperature [44]. We note that these calculations are rather sensitive since we are close to a pole in the fluid equations. Thus, we expect that the difference in these energy confinement scalings (c.f. Ref. [4]) can be due to impurities in the ITER simulations, where Z eff = 1 . 65 . We did not include impurities in the local code.
In summary, this section examines isotope effects on transport and confinement at the plasma edge. It highlights the importance of edge data, as transport flux must pass through the edge, and shows that the H-mode barrier depends on zonal flows. Ion thermal diffusivity decreases with higher isotope mass. The L-H transition, driven by poloidal rotation, enables the use of analytic transport calculations. A switch from hydrogen to tritium results in a 6% reduction in transport, with confinement time ( τ E ) scaling as A 0 . 2 , consistent with experimental data. However, ITER simulations predict a different scaling of τ E = A 0 . 5 , possibly due to impurities not accounted for in local code simulations. The sensitivity of these calculations is noted, particularly near poles in the fluid equations.

6. Density Limit and System Size Dependence: Insights from Gyro-Landau Resonances and E × B Convection

We now recall from Ref. [7] that there is hardly any rotation if we include the unaveraged gyro-Landau resonances. Then we get much weaker zonal flow, which could hardly absorb the inverse turbulent cascade. Then we get a pileup of waves at the system size, leading to very strong transport, which is in accordance with the density limit Equation (11). Now we know that the gyro-Landau resonances are typically about an order of magnitude larger than the ion viscosity. However, ion viscosity will increase with density, so if we increase density by about a factor 10, we reach a density limit [38].
n G = I p π a 2 ,
where I p the plasma current, and a the minor radius. Our first observation is the strong dependence on the system size, meaning a strong potential dependence on perturbations with radial scale length approaching the system size. Now, the stabilization of an instability due to E × B convection is written:
γ δ T = v E · ∇ δ T .
leading to
e δ ϕ T e = γ ω * e 1 k x L n ,
where k x is the radial wavenumber and L n is the density gradient scale length. Note that in Equation 13, the diamagnetic drift accounts only for the density gradient. Thus, the density gradient is canceled (see, e.g., Ref. [34]). The saturation level in Equation (13) is in agreement with typical experiments and used both by us and in other works. A very critical point is here whether the nonlinearity in Equation (12) is entirely stabilizing. This is only the case if we are looking at the correlation length (corresponding to the mode number with the largest growth rate as normalized by the drift frequency) and if there are no reflections in k-space. This may be critical for the inverse cascade, leading to cascade towards longer wavelengths, i.e., towards the system size. Clearly, it is very important to take into account that we are here looking at systems with a finite radial size. To avoid getting a pileup of waves at the system size with reflections, we need strong zonal flows that absorb the inverse cascade. This is usually only available in reactive systems, as seen in Ref. [18] for only density transport and in Ref. [17] using our reactive fluid closure. We note that in both cases, we obtain a transport kernel of the type
D = γ 3 / k x 2 ω r 2 + γ 2 .
In Ref. [18], we have an exact fluid closure by taking a zero-ion temperature, while in our case, we have no kinetic dissipation because of resonance broadening [7,13]. In our case Equation (12) is applied to ion temperature. In our derivation, we also included pinch terms from a full quasilinear procedure [17]. The quasilinear approach was found in Ref. [29] to be valid within a few percent in the fluid case. We now recall that if the ion density is increased by about a factor 100, the ion viscosity in Equation 6 becomes comparable to an unnormalized gyro-Landau resonance that was, in fact, displayed in Ref. [7]. The result for the rotation was similar to that in Figure 5 below, where the density was only increased by a factor 10 above the Cyclone case. The result of the Cyclone work was an increase of up to a factor 3 in thermal conductivity for an unnormalized gyro-Landau fluid model. However, recent gyro-Landau fluid models have been normalized to the nonlinear gyrokinetic code [11].
Consequently, if the density is increased by approximately ten times compared to the density investigated in the Cyclone work [10], it is expected that there will be an increase in thermal transport, a decrease in flow shear, and the disappearance of the Dimits shift. The increase in transport will be accompanied by the breakdown of the H-mode barrier and the appearance of significant turbulence structures, known as MARFES (for multifaceted asymmetric radiation from the edge), similar to those observed in experimental studies at the density limit. Due to the normally larger sizes of MHD modes compared to drift waves, it is anticipated that MHD ballooning, kinetic ballooning modes, resistive ballooning modes, and peeling modes will also be present near the edge.
This section explores the density limit and system size dependence in plasma transport, focusing on gyro-Landau resonances and E × B convection. Including unaveraged gyro-Landau resonances results in weaker zonal flows, leading to increased transport and aligning with the density limit equation. Gyro-Landau resonances are significantly larger than ion viscosity, but ion viscosity increases with density, reaching the density limit when density is increased by about a factor of ten. Stabilization of instabilities due to E × B convection is crucial, with saturation levels matching experimental observations. Strong zonal flows are necessary to prevent wave pileup at system size, which is typically available in reactive systems. Increasing ion density significantly raises thermal transport, potentially breaking the H-mode barrier and causing substantial turbulence structures, similar to those seen in experiments at the density limit.

7. Summary

In this work, we have achieved two additional successful applications of our transport model: isotope scaling and the density limit. Regarding isotope scaling, we believe it fits perfectly. However, for the density limit [38,39], the evidence is less accurate. We found that increasing the density in our model to the gyro-Landau limit, giving viscosity comparable to the gyro-Landau resonance, would require a density far exceeding the density limit. We also discuss the implications of encountering MARFEs and ballooning instability in this scenario. MARFEs are very large-scale structures and can therefore be torn apart by zonal flows.
We recall the derivation and development of our fluid model, first established in 1988 using theory and simulations with an absorbing boundary for long wavelengths. The original model remains valid within the appropriate limit. The motivation for the absorbing boundary for long wavelengths is now supported by the damping of zonal flows. This first-principles model is not tailored to any other theoretical model or experiment. A notable extension of the model includes the incorporation of fast particles, applying resonance broadening for fast particle modes—an aspect utilized by other research groups.
From a broader perspective, transport in tokamaks is observed to be governed by a delicate balance between turbulence and rotation. In our model, similar to any fluid model, both turbulence and rotation (resulting from turbulence) originate at the resonance in the ion temperature perturbation. This resonance defines the linear threshold of the ITG and acts as the primary source of rotation. We assert that this interplay fundamentally controls transport. Given that both our isotope scaling and density limit rely predominantly on this resonance, we consider them fundamental. While acknowledging the existence of several other mechanisms contributing to both isotope scaling and density limits, we maintain confidence in the fundamentality of our proposed mechanisms. This confidence arises from their close relation to our fluid closure, previously shown to give excellent agreement with experiment and from their integral role in overall transport control, leading us to expect their dominance.

Author Contributions

Jan Weiland: Writing – original draft (equal); Writing - review and editing (equal). Tariq Rafiq: Funding acquisition (lead); Writing - original draft (equal); Writing - review and editing (equal). Eugenio Schuster: Writing - review and editing (support).

Funding

This material is based upon work supported by the U.S. Department of Energy, Office of Science, Office of Fusion Energy Sciences, under Award Numbers DE-SC0013977 and DE-SC0010661.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

One of us (JW) wants to express his gratitude to the Bogoliubov Institute of Theoretical Physics, Kiev, Ukraine, for developing the nonlinear kinetic theory on which our main derivations are based. The authors are grateful to Dr. HG Gustavsson for assistance with computation.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Conner, J.W.; Wilson, H.R. Survey of theories of anomalous transport. Plasma Physics and Controlled Fusion 1994, 36, 719. [CrossRef]
  2. Liewer, P.C. Measurements of microturbulence in tokamaks and comparisons with theories of turbulence and anomalous transport. Nuclear Fusion 1985, 25, 543. [CrossRef]
  3. Itoh, S.I.; Itoh, K. Hydrogen isotope effect on the Dimits shift. Nuclear Fusion 2016, 56, 106028. [CrossRef]
  4. Urano, H.; Takizuka, T.; Kikuchi, M.; Nakano, T.; Hayashi, N.; Oyama, N.; Kamada, Y. Publisher’s Note: Small ion-temperature-gradient scale length and reduced heat diffusivity at large hydrogen isotope mass in conventional H-mode plasmas [Phys. Rev. Lett. 109, 125001 (2012)]. Phys. Rev. Lett. 2012, 109, 149901. [CrossRef]
  5. Rosenbluth, M.N.; Hinton, F.L. Poloidal Flow Driven by Ion-Temperature-Gradient Turbulence in Tokamaks. Phys. Rev. Lett. 1998, 80, 724–727. [CrossRef]
  6. Wang, G.Q.; Ma, J.; Weiland, J. Zonal flows near marginal stability in drift wave transport. Physica Scripta 2015, 90, 065604. [CrossRef]
  7. Weiland, J.; Zagorodny, A. On the normalization of transport from ITG Modes. Physics of Plasmas 2016, 23, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.4964772/14027143/102307_1_online.pdf]. 102307, . [CrossRef]
  8. Nozaki, K.; Taniuti, T. Plasma Flow and a Soliton in a Theta Pinch. Journal of the Physical Society of Japan 1979, 46, 970–974. [CrossRef]
  9. Wakatani, M.; Hasegawa, A. A collisional drift wave description of plasma edge turbulence. The Physics of Fluids 1984, 27, 611–618, [https://pubs.aip.org/aip/pfl/article-pdf/27/3/611/12476138/611_1_online.pdf]. [CrossRef]
  10. Dimits, A.M.; Bateman, G.; Beer, M.A.; Cohen, B.I.; Dorland, W.; Hammett, G.W.; Kim, C.; Kinsey, J.E.; Kotschenreuther, M.; Kritz, A.H.; Lao, L.L.; Mandrekas, J.; Nevins, W.M.; Parker, S.E.; Redd, A.J.; Shumaker, D.E.; Sydora, R.; Weiland, J. Comparisons and physics basis of tokamak transport models and turbulence simulations. Physics of Plasmas 2000, 7, 969–983, [https://pubs.aip.org/aip/pop/article-pdf/7/3/969/12423001/969_1_online.pdf]. [CrossRef]
  11. Doyle, E.; Houlberg, W.; Kamada, Y.; Mukhovatov, V.; Osborne, T.; Polevoi, A.; Bateman, G.; Connor, J.; Cordey, J.; Fujita, T.; Garbet, X.; Hahm, T.; Horton, L.; Hubbard, A.; Imbeaux, F.; Jenko, F.; Kinsey, J.; Kishimoto, Y.; Li, J.; Luce, T.; Martin, Y.; Ossipenko, M.; Parail, V.; Peeters, A.; Rhodes, T.; Rice, J.; Roach, C.; Rozhansky, V.; Ryter, F.; Saibene, G.; Sartori, R.; Sips, A.; Snipes, J.; Sugihara, M.; Synakowski, E.; Takenaga, H.; Takizuka, T.; Thomsen, K.; Wade, M.; Wilson, H.; ITPA Transport Physics Topical Group and ITPA Confinement Database and Modelling Topical Group and ITPA Pedestal and Edge Topical Group. Chapter 2: Plasma confinement and transport. Nuclear Fusion 2007, 47, S18. [CrossRef]
  12. Weiland, J.; Rafiq, T.; Schuster, E. Fast particles in drift wave turbulence. Physics of Plasmas 2023, 30, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/5.0147320/17009722/042517_1_5.0147320.pdf]. 042517, . [CrossRef]
  13. Hasegawa, A. Plasma instabilities and nonlinear effects. Springer Verlag Springer Series on Physics Chemistry Space 1975, 8.
  14. Waltz, R.; Dominguez, R.; Hammett, G. Gyro-Landau fluid models for toroidal geometry. Physics of Fluids b: Plasma physics 1992, 4, 3138–3151.
  15. Waltz, R.E.; Staebler, G.M.; Dorland, W.; Hammett, G.W.; Kotschenreuther, M.; Konings, J.A. A gyro-Landau-fluid transport model. Physics of Plasmas 1997, 4, 2482–2496, [https://pubs.aip.org/aip/pop/article-pdf/4/7/2482/12447394/2482_1_online.pdf]. [CrossRef]
  16. Hasegawa, A.; Mima, K. Pseudo-three-dimensional turbulence in magnetized nonuniform plasma. The Physics of Fluids 1978, 21, 87–92, [https://pubs.aip.org/aip/pfl/article-pdf/21/1/87/12384514/87_1_online.pdf]. [CrossRef]
  17. Weiland, J.; Nordman, H. Proc. Varenna-Lausanne Workshop, Chexbres 1988 "Theory of Fusion Plasmas. (Bologna Editrice Compostori) 1988, p. 451.
  18. Connor, J.W.; Pogutse, O.P. On the relationship between mixing length and strong turbulence estimates for transport due to drift turbulence. Plasma Physics and Controlled Fusion 2001, 43, 155. [CrossRef]
  19. Dupree, T.H. A Perturbation Theory for Strong Plasma Turbulence. The Physics of Fluids 1966, 9, 1773–1782, [https://aip.scitation.org/doi/pdf/10.1063/1.1761932]. [CrossRef]
  20. Hahm, T.; Wang, L.; Wang, W.; Yoon, E.; Duthoit, F. Isotopic dependence of residual zonal flows. Nuclear Fusion 2013, 53, 072002. [CrossRef]
  21. Sugama, H.; Watanabe, T.H.; Horton, W. Collisionless kinetic-fluid model of zonal flows in toroidal plasmas. Physics of Plasmas 2007, 14, 022502, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.2435329/16121570/022502_1_online.pdf]. [CrossRef]
  22. Kim, S.S.; Jhang, H. A conservative gyrofluid model: Effect of closure on energetics. Physics of Plasmas 2020, 27, 092305, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/5.0020332/15811813/092305_1_online.pdf]. [CrossRef]
  23. Zagorodny, A.; Weiland, J. Statistical theory of turbulent transport (non-Markovian effects). Physics of Plasmas 1999, 6, 2359–2372.
  24. Weiland, J.; Wilhelmsson, H. Coherent non-linear interaction of waves in plasmas. Oxford Pergamon Press International Series on Natural Philosophy 1977, 88.
  25. Weiland, J.; Wilhelmsson, H. Repetitive Explosive Instabilities. Physica Scripta 1973, 7, 222. [CrossRef]
  26. Mattor, N.; Parker, S.E. Nonlinear Kinetic Fluid Equations. Phys. Rev. Lett. 1997, 79, 3419–3422. [CrossRef]
  27. Holod, I.; Weiland, J.; Zagorodny, A. Nonlinear fluid closure:Three mode slab ion temperature gradient problem with diffusion. Physics of Plasmas 2002, 9, 1217–1220, [. [CrossRef]
  28. Weiland, J. Simulations of the L-H transition on experimental advanced superconducting Tokamak. Physics of Plasmas 2014, 21, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.4901597/15780223/122501_1_online.pdf]. 122501, . [CrossRef]
  29. Weiland, J. The role of zonal flows in reactive fluid closures. Plasma Science and Technology 2018, 20, 074007. [CrossRef]
  30. Weiland, J.; Zagorodny, J. Drift wave theory for transport in tokamaks. Rev. Mod. Plasma Phys 2019, 3, 8. [CrossRef]
  31. Rogers, B.N.; Drake, J.F.; Zeiler, A. Phase Space of Tokamak Edge Turbulence, the L-H Transition, and the Formation of the Edge Pedestal. Phys. Rev. Lett. 1998, 81, 4396–4399. [CrossRef]
  32. Hubbard, A.E.; Boivin, R.L.; Drake, J.F.; Greenwald, M.; In, Y.; Irby, J.H.; Rogers, B.N.; Snipes, J.A. Local variables affecting H-mode threshold on Alcator C-Mod. Plasma Physics and Controlled Fusion 1998, 40, 689. [CrossRef]
  33. Hubbard, A.E.; Hughes, J.W.; Bespamyatnov, I.O.; Biewer, T.; Cziegler, I.; LaBombard, B.; Lin, Y.; McDermott, R.; Rice, J.E.; Rowan, W.L.; Snipes, J.A.; Terry, J.L.; Wolfe, S.M.; Wukitch, S.; The Alcator C-Mod Group. H-mode pedestal and threshold studies over an expanded operating space on Alcator C-Mod. Physics of Plasmas 2007, 14, 056109, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.2714297/13888117/056109_1_online.pdf]. [CrossRef]
  34. Weiland, J. Stability and Transport in Magnetic Confinement Systems; Springer, New York, Heidelberg, 2012.
  35. Snyder, P.; Groebner, R.; Hughes, J.; Osborne, T.; Beurskens, M.; Leonard, A.; Wilson, H.; Xu, X. A first-principles predictive model of the pedestal height and width: development, testing and ITER optimization with the EPED model. Nuclear Fusion 2011, 51, 103016. [CrossRef]
  36. ITER Physics Expert Group on Confinement and Transport and ITER Physics Expert Group on Confinement Modelling and Database and ITER Physics Basis Editors. Chapter 2: Plasma confinement and transport. Nuclear Fusion 1999, 39, 2175. [CrossRef]
  37. Weiland, J.; Nordman, H. Enhanced confinement regimes in transport code simulations of toroidal drift wave transport. Nuclear Fusion 1991, 31, 390. [CrossRef]
  38. Greenwald, M.; Terry, J.; Wolfe, S.; Ejima, S.; Bell, M.; Kaye, S.; Neilson, G. A new look at density limits in tokamaks. Nuclear Fusion 1988, 28, 2199. [CrossRef]
  39. Borrass, K.; Loarte, A.; Maggi, C.; Mertens, V.; Monier, P.; Monk, R.; Ongena, J.; Rapp, J.; Saibene, G.; Sartori, R.; Schweinzer, J.; Stober, J.; Suttrop, W.; Workprogramme collaborators, E. Recent H-mode density limit studies at JET. Nuclear Fusion 2004, 44, 752. [CrossRef]
  40. Gates, D.A.; Delgado-Aparicio, L. Origin of Tokamak Density Limit Scalings. Phys. Rev. Lett. 2012, 108, 165004. [CrossRef]
  41. Zanca, P.; Sattin, F.; Escande, D.; Pucella, G.; Tudisco, O. A unified model of density limit in fusion plasmas. Nuclear Fusion 2017, 57, 056010. [CrossRef]
  42. Giacomin, M.; Pau, A.; Ricci, P.; Sauter, O.; Eich, T.; the ASDEX Upgrade team.; Contributors, J.; the TCV team. First-Principles Density Limit Scaling in Tokamaks Based on Edge Turbulent Transport and Implications for ITER. Phys. Rev. Lett. 2022, 128, 185003. [CrossRef]
  43. Singh, R.; Diamond, P.H. Zonal shear layer collapse and the power scaling of the density limit: old L-H wine in new bottles. Plasma Physics and Controlled Fusion 2022, 64, 084004. [CrossRef]
  44. Rafiq, T.; Weiland, J. Self-consistent core-pedestal ITER scenario modeling. Nuclear Fusion 2021, 61, 116005. [CrossRef]
  45. Rafiq, T.; Kritz, A.H.; Tangri, V.; Pankin, A.Y.; Voitsekhovitch, I.; Budny, R.V.; JET EFDA Contributors. Integrated modeling of temperature profiles in L-mode tokamak discharges. Physics of Plasmas 2014, 21, 122505, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.4903464/15780800/122505_1_online.pdf]. [CrossRef]
  46. Rafiq, T.; Wang, Z.; Morosohk, S.; Schuster, E.; Weiland, J.; Choi, W.; Kim, H.T. Validating the Multi-Mode Model’s Ability to Reproduce Diverse Tokamak Scenarios. Plasma 2023, 6, 435–458. [CrossRef]
  47. Rafiq, T.; Wilson, C.; Clauser, C.; Schuster, E.; Weiland, J.; Anderson, J.; Kaye, S.; Pankin, A.; LeBlanc, B.; Bell, R. Predictive modeling of NSTX discharges with the updated multi-mode anomalous transport module. Nuclear Fusion 2024, 64, 076024. [CrossRef]
  48. Weiland, J.; Rafiq, T.; Schuster, E. Nonlinearities in magnetic confinement, ionospheric physics, and population explosion leading to profile resilience. submitted to Physics of Plasmas 2024.
  49. Rafiq, T.; Kritz, A.H.; Weiland, J.; Pankin, A.Y.; Luo, L. Physics basis of Multi-Mode anomalous transport module. Physics of Plasmas 2013, 20, 032506, [https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.4794288/14795596/032506_1_online.pdf]. [CrossRef]
  50. Houlberg, W.A.; Shaing, K.C.; Hirshman, S.P.; Zarnstorff, M.C. Bootstrap current and neoclassical transport in tokamaks of arbitrary collisionality and aspect ratio. Physics of Plasmas 1997, 4, 3230–3242, [https://pubs.aip.org/aip/pop/article-pdf/4/9/3230/12664245/3230_1_online.pdf]. [CrossRef]
  51. Ding, S.; Garofalo, A.M.; Wang, H.Q.; Weisberg, D.B.; Li, Z.Y.; Jian, X.; Eldon, D.; Victor, B.S.; Marinoni, A.; Hu, Q.M.; Carvalho, I.S.; Odstrčil, T.; Wang, L.; Hyatt, A.W.; Osborne, T.H.; Gong, X.Z.; Qian, J.P.; Huang, J.; McClenaghan, J.; Holcomb, C.T.; Hanson, J.M. A high-density and high-confinement tokamak plasma regime for fusion energy. Nature 2024, 629, 555–560.
  52. Zagorodny, A.; Weiland, J. Closure at the Irreducible Part of the Fourth Moment for the Case of Constant Coefficients in the Fokker-Planck Equation. AIP Conference Proceedings 2011, 1392, 24–32, [https://pubs.aip.org/aip/acp/article-pdf/1392/1/24/11579434/24_1_online.pdf]. [CrossRef]
  53. Chandrasekhar, S. Stochastic Problems in Physics and Astronomy. Rev. Mod. Phys. 1943, 15, 1–89. [CrossRef]
  54. Coppi, B.; Rosenbluth, M.; Sudan, R. Nonlinear interactions of positive and negative energy modes in rarefied plasmas (I). Annals of Physics 1969, 55, 207–247. [CrossRef]
Preprints 114717 g001
Figure 1. This figure exhibits similarity to the suppression of explosive instability by a nonlinear frequency shift [27]. Reproduced from [I. Holod, J. Weiland, and A. Zagorodny Physics of Plasmas 9, 1217 (2002)], with the permission of AIP Publishing.
Figure 1. This figure exhibits similarity to the suppression of explosive instability by a nonlinear frequency shift [27]. Reproduced from [I. Holod, J. Weiland, and A. Zagorodny Physics of Plasmas 9, 1217 (2002)], with the permission of AIP Publishing.
Preprints 114717 g001
Preprints 114717 g002
Figure 2. (a) Normalized growth rate and normalized flow shear (b) transport (ion thermal diffusivity) are derived for Cyclone parameters based on the findings in Ref. [7], incorporating ion viscosity corresponding to hydrogen in Equation 6. The strong flow shear at marginal stability arises from the fluid resonance described in Equation 6. Waves with wavelengths approaching the system size will inevitably reach marginal stability at some point, leading to strong damping. Consequently, transport is heavily influenced by the fluid closure. The reactive closure results from the detuning of wave-particle resonances due to resonance broadening or nonlinear frequency shifts.
Figure 2. (a) Normalized growth rate and normalized flow shear (b) transport (ion thermal diffusivity) are derived for Cyclone parameters based on the findings in Ref. [7], incorporating ion viscosity corresponding to hydrogen in Equation 6. The strong flow shear at marginal stability arises from the fluid resonance described in Equation 6. Waves with wavelengths approaching the system size will inevitably reach marginal stability at some point, leading to strong damping. Consequently, transport is heavily influenced by the fluid closure. The reactive closure results from the detuning of wave-particle resonances due to resonance broadening or nonlinear frequency shifts.
Preprints 114717 g002
Preprints 114717 g003
Figure 3. Transport as a function of normalized ion temperature gradients for both hydrogen and tritium shows a decrease in transport with a higher isotope mass. These results are obtained by adding ion viscosity to the calculations conducted in Ref. [7].
Figure 3. Transport as a function of normalized ion temperature gradients for both hydrogen and tritium shows a decrease in transport with a higher isotope mass. These results are obtained by adding ion viscosity to the calculations conducted in Ref. [7].
Preprints 114717 g003
Preprints 114717 g004
Figure 4. The ion temperature is depicted by the full line before and after the L-H transition in a simulation of ITER, utilizing the Weiland fluid model along with the neoclassical module.
Figure 4. The ion temperature is depicted by the full line before and after the L-H transition in a simulation of ITER, utilizing the Weiland fluid model along with the neoclassical module.
Preprints 114717 g004
Preprints 114717 g005
Figure 5. (a) Normalized growth rates and normalized flow shear (b) transport (ion thermal diffusivity) illustrate similar variations as seen in Figure 2. Both the rotation and Dimits shift exhibit a decreasing trend, consistent with the observations in Ref. [7] for the gyro-Landau fluid model. This particular case corresponds to a density of 4 . 5 × 10 20 m − 3 .
Figure 5. (a) Normalized growth rates and normalized flow shear (b) transport (ion thermal diffusivity) illustrate similar variations as seen in Figure 2. Both the rotation and Dimits shift exhibit a decreasing trend, consistent with the observations in Ref. [7] for the gyro-Landau fluid model. This particular case corresponds to a density of 4 . 5 × 10 20 m − 3 .
Preprints 114717 g005
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.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated