No Arabic abstract
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.
The role of turbulence in current generation and self-excitation of magnetic fields has been studied in the geometry of a mechanically driven, spherical dynamo experiment, using a three dimensional numerical computation. A simple impeller model drives a flow which can generate a growing magnetic field, depending upon the magnetic Reynolds number, Rm, and the fluid Reynolds number. When the flow is laminar, the dynamo transition is governed by a simple threshold in Rm, above which a growing magnetic eigenmode is observed. The eigenmode is primarily a dipole field tranverse to axis of symmetry of the flow. In saturation the Lorentz force slows the flow such that the magnetic eigenmode becomes marginally stable. For turbulent flow, the dynamo eigenmode is suppressed. The mechanism of suppression is due to a combination of a time varying large-scale field and the presence of fluctuation driven currents which effectively enhance the magnetic diffusivity. For higher Rm a dynamo reappears, however the structure of the magnetic field is often different from the laminar dynamo; it is dominated by a dipolar magnetic field which is aligned with the axis of symmetry of the mean-flow, apparently generated by fluctuation-driven currents. The fluctuation-driven currents have been studied by applying a weak magnetic field to laminar and turbulent flows. The magnetic fields generated by the fluctuations are significant: a dipole moment aligned with the symmetry axis of the mean-flow is generated similar to those observed in the experiment, and both toroidal and poloidal flux expulsion are observed.
Using fully kinetic simulations, we study the suppression of asymmetric reconnection in the limit where the diamagnetic drift speed >> Alfven speed and the magnetic shear angle is moderate. We demonstrate that the slippage between electrons and the magnetic flux facilitates reconnection, and can even result in fast reconnection that lacks one of the outflow jets. Through comparing a case where the diamagnetic drift is supported by the temperature gradient with a companion case that has a density gradient instead, we identify a robust suppression mechanism. The drift of the x-line is slowed down locally by the asymmetric nature of the current sheet and the resulting tearing modes, then the x-line is run over and swallowed by the faster-moving following flux.
Sustained ELM mitigation has been achieved using RMPs with a toroidal mode number of n=4 and n=6 in lower single null and with n=3 in connected double null plasmas on MAST. The ELM frequency increases by up to a factor of eight with a similar reduction in ELM energy loss. A threshold current for ELM mitigation is observed above which the ELM frequency increases approximately linearly with current in the coils. A comparison of the filament structures observed during the ELMs in the natural and mitigated stages shows that the mitigated ELMs have the characteristics of type I ELMs even though their frequency is higher, their energy loss is reduced and the pedestal pressure gradient is decreased. During the ELM mitigated stage clear lobe structures are observed in visible-light imaging of the X-point region. The size of these lobes is correlated with the increase in ELM frequency observed. The RMPs produce a clear 3D distortion to the plasma and it is likely that these distortions explain why ELMs are destabilised and hence why ELM mitigation occurs.
The application of resonant magnetic perturbations (RMPs) produces splitting of the divertor strike point due to the interaction of the RMP field and the plasma field. The application of a rotating RMP field causes the strike point splitting to rotate, distributing the particle and heat flux evenly over the divertor. The RMP coils in MAST have been used to generate a rotating perturbation with a toroidal mode number n=3. The ELM frequency is doubled with the application of the RMP rotating field, whilst maintaining the H mode. During mitigation, the ELM peak heat flux is seen to be reduced by 50% for a halving in the ELM energy and motion of the strike point, consistent with the rotation of the applied RMP field, is seen using high spatial resolution (1.5mm at the target) heat flux profiles measured using infrared (IR) thermography.
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.