Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userObservations of Heteroclinic Bifurcations in Resistive MHD Simulations of the Plasma Response to Resonant Magnetic Perturbations
123
0
0.0
(
0
)
Authors
T. E. Evans
Ask ChatGPT about the research
No Arabic abstract
A new class of static magnetohydrodynamic (MHD) magnetic island bifurcations is identified in rotating spherical tokamak plasmas during single- and two-fluid resistive MHD simulations. As the magnitude of an externally applied non-axisymmetric magnetic field perturbation is increased in these simulations, the internal flux surfaces that make up a sub-set of the resonant helical magnetic islands in the plasma gradually elongate and undergo heteroclinic bifurcations. The bifurcation results in the creation of a new set of hyperbolic-elliptic fixed points as predicted by the Poincare-Birkoff fixed point theorem. Field line calculations without including the resistive MHD plasma response to the applied perturbation field do not undergo this class of bifurcations indicating the importance of plasma self-organization in the bifurcation process.
rate research
Read More
Plasma response to resonant magnetic perturbations (RMPs) near the pedestal top is crucial for accessing edge localized modes (ELMs) suppression in tokamaks. Since radial location of rotation zero-crossing plays a key role in determining the threshold for field penetration of RMP, plasma response may be different in low input torque plasmas. In this work, the linear MHD code MARS-F is applied to reveal the dependence of plasma response to RMP on rotation zero-crossing by a scan of rotation profiles based on an EAST equilibrium. It is shown that the plasma response is enhanced when zero-crossing occurs near rational surfaces. The dependence of plasma response on the location of rotation zero-crossing is well fitted by a double Gaussian, indicating two effects in this enhancement. One is induced by rotation screening effect shown as a wide base (with a width around 10-20 krad/s), and the other is related to resistive singular layer effect characterized by a localized peak (with a width around 3-4 krad/s). The width of the peak scales with the resistive singular layer width. The plasma displacement suggests the response is tearing like when zero-crossing is within the singular layer, while it is kink like when zero-crossing is far from the layer. The enhancement of magnetic islands width at the peak is only around a factor of two, when the absolute value of local rotation is not larger than 10-20 krad/s. It is further confirmed in a modeling of plasma response in an EAST ELM suppression discharge. Though there is a zero-crossing in $Etimes B$ rotation but not in electron perpendicular rotation, no significant difference in plasma response is obtained using these two rotation profiles. This suggests that the rotation near pedestal top should not be far away from zero but may not be necessary to have zero-crossing for accessing ELM suppression.
The plasma response from an external n = 2 magnetic perturbation field in ASDEX Upgrade has been measured using mainly electron cyclotron emission (ECE) diagnostics and a rigid rotating field. To interpret ECE and ECE-imaging (ECE-I) measurements accurately, forward modeling of the radiation transport has been combined with ray tracing. The measured data is compared to synthetic ECE data generated from a 3D ideal magnetohydrodynamics (MHD) equilibrium calculated by VMEC. The measured amplitudes of the helical displacement around the low field side midplane are in reasonable agreement with the one from the synthetic VMEC diagnostics. Both exceed the prediction from the vacuum field calculations and indicate the presence of a kink response at the edge, which amplifies the perturbation. VMEC and MARS-F have been used to calculate the properties of this kink mode. The poloidal mode structure of the magnetic perturbation of this kink mode at the edge peaks at poloidal mode numbers larger than the resonant components |m| > |nq|, whereas the poloidal mode structure of its displacement is almost resonant |m| ~ |nq|. This is expected from ideal MHD in the proximity of rational surfaces. The displacement measured by ECE-I confirms this resonant response.
The effect of magnetic perturbations (MPs) on the helical self-organization of shaped tokamak plasmas is discussed in the framework of the nonlinear 3D MHD model. Numerical simulations performed in toroidal geometry with the textsc{pixie3d} code [L. Chacon, Phys. Plasmas {bf 15}, 056103 (2008)] show that $n=1$ MPs significantly affect the spontaneous quasi-periodic sawtoothing activity of such plasmas. In particular, the mitigation of sawtooth oscillations is induced by $m/n=1/1$ and $2/1$ MPs. These numerical findings provide a confirmation of previous circular tokamak simulations, and are in agreement with tokamak experiments in the RFX-mod and DIII-D devices. Sawtooth mitigation via MPs has also been observed in reversed-field pinch simulations and experiments. The effect of MPs on the stochastization of the edge magnetic field is also discussed.
The impact of resonant magnetic perturbations (RMPs) on the power required to access H-mode is examined experimentally on MAST. Applying RMP in n=2,3,4 and 6 configurations causes significant delays to the timing of the L-H transition at low applied fields and prevents the transition at high fields. The experiment was primarily performed at RMP fields sufficient to cause moderate increases in ELM frequency, f mitigated/f natural~3. To obtain H-mode with RMPs at this field, an increase of injected beam power is required of at least 50% for n=3 and n=4 RMP and 100% for n=6 RMP. In terms of power threshold, this corresponds to increases of at least 20% for n=3 and n=4 RMPs and 60% for n=6 RMPs. This RMP affected power threshold is found to increase with RMP magnitude above a certain minimum perturbed field, below which there is no impact on the power threshold. Extrapolations from these results indicate large increases in the L-H power threshold will be required for discharges requiring large mitigated ELM frequency.
General three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal MHD limit has been a significant challenge. This work presents numerical solutions of the Hahm-Kulsrud-Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad-Shafranov (GS) solver and the SPEC code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), where the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution approaches the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate excellent agreement between the numerical solutions obtained from the two codes through a thorough convergence study.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
