No Arabic abstract
The role of anisotropic thermal diffusivity on tearing mode stability is analysed in general toroidal geometry. A dispersion relation linking the growth rate to the tearing mode stability parameter, Delta, is derived. By using a resistive MHD code, modified to include such thermal transport, to calculate tearing mode growth rates, the dispersion relation is employed to determine Delta in situations with finite plasma pressure that are stabilised by favourable average curvature in a simple resistive MHD model. We also demonstrate that the same code can also be used to calculate the basis-functions [C J Ham, et al, Plasma Phys. Control. Fusion 54 (2012) 105014] needed to construct Delta.
The stabilization of tearing modes with rf driven current benefits from the cooperative feedback loop between rf power deposition and electron temperature within the island. This effect, termed rf current condensation, can greatly enhance and localize current driven within magnetic islands. It has previously been shown that the condensation effect opens the possibility of passive stabilization with broad rf profiles, as would be typical of LHCD for steady state operation. Here we show that this self-healing effect can be dramatically amplified by operation in a hot ion mode, due to the additional electron heat source provided by the hotter ions.
In large hot tokamaks like JET, the width of the reconnecting layer for resistive modes is determined by semi-collisional electron dynamics and is much less than the ion Larmor radius. Firstly a dispersion relation valid in this regime is derived which provides a unified description of drift-tearing modes, kinetic Alfven waves and the internal kink mode at low beta. Tearing mode stability is investigated analytically recovering the stabilising ion orbit effect, obtained previously by Cowley et al. [Phys. Fluids (29) 3230 1986], which implies large values of the tearing mode stability parameter Delta prime are required for instability. Secondly, at high beta it is shown that the tearing mode interacts with the kinetic Alfven wave and that there is an absolute stabilisation for all Delta prime due to the shielding effects of the electron temperature gradients, extending the result of Drake et. al [Phys. Fluids (26) 2509 1983] to large ion orbits. The nature of the transition between these two limits at finite values of beta is then elucidated. The low beta formalism is also relevant to the m=n=1 tearing mode and the dissipative internal kink mode, thus extending the work of Pegoraro et al. [Phys. Fluids B (1) 364 1989] to a more realistic electron model incorporating temperature perturbations, but then the smallness of the dissipative internal kink mode frequency is exploited to obtain a new dispersion relation valid at arbitrary beta. A diagram describing the stability of both the tearing mode and dissipative internal kink mode, in the space of Delta prime and beta, is obtained. The trajectory of the evolution of the current profile during a sawtooth period can be plotted in this diagram, providing a model for the triggering of a sawtooth crash.
Calculations of tearing mode stability in tokamaks split conveniently into an external region, where marginally stable ideal MHD is applicable, and a resonant layer around the rational surface where sophisticated kinetic physics is needed. These two regions are coupled by the stability parameter. Pressure and current perturbations localized around the rational surface alter the stability of tearing modes. Equations governing the changes in the external solution and - are derived for arbitrary perturbations in axisymmetric toroidal geometry. The relationship of - with and without pressure flattening is obtained analytically for four pressure flattening functions. Resistive MHD codes do not contain the appropriate layer physics and therefore cannot predict stability directly. They can, however, be used to calculate -. Existing methods (Ham et al. 2012 Plasma Phys. Control. Fusion 54 025009) for extracting - from resistive codes are unsatisfactory when there is a finite pressure gradient at the rational surface and favourable average curvature because of the Glasser stabilizing effect. To overcome this difficulty we introduce a specific pressure flattening function that allows the earlier approach to be used. The technique is first tested numerically in cylindrical geometry with an artificial favourable curvature. Its application to toroidal geometry is then demonstrated using the toroidal tokamak tearing mode stability code T7 (Fitzpatrick et al. 1993 Nucl. Fusion 33 1533) which uses an approximate analytic equilibrium. The prospects for applying this approach to resistive MHD codes such as MARS-F (Liu et al. 2000 Phys. Plasmas 7 3681) which utilize a fully toroidal equilibrium are discussed.
The RF stabilization of tearing modes with current condensation has the potential to increase stabilization efficiency and loosen power localization requirements. Such benefits stem from the cooperative feedback between the RF deposition and resulting island temperature perturbation governed by diffusion. A self consistent treatment of the damping of an rf ray as it traverses the island shows that low damping scenarios can require unfavorably high powers to overcome initial power leakage and effectively capitalize on the nonlinear effect. In this work it is demonstrated that for such regimes,modulated stabilization schemes can achieve significant improvements in heating and current drive contributions to stabilization for the same average power as a continuous wave scheme. The impact of modulation frequency and duty cycle on the performance is explored, the results of which suggest modulation strategies in which the pulsing periods are kept on the order of a diffusive time.
Direct evidence of micro-turbulence effect on the onset of neoclassical tearing mode (NTM) is reported for the first time in this letter. A puzzling positive correlation between critical width of seed island of NTM and normalized plasma pressure beta_p is first observed employing a novel method for clearly separating the processes of seed island and the onset of NTM in the EAST tokamak. Different from the methods developed before, the width of the seed island is well controlled by slowly ramping up the current in resonant magnetic perturbation coils. It is revealed that the positive correlation is mainly attributed to the enhancement of perpendicular transport by micro-turbulence, which overcomes the destabilizing effect of beta_p on the onset of NTM. Reduced magnetohydrodynamics (MHD) modeling well reproduced the two states of nonlinear bifurcations observed in this experiment by including the finite transport effect. This result provides a new route for understanding multi-scale interaction in plasma physics.