Magnetic Islands Research Articles

A semi-analytical model of RF condensation that can handle localized power depositions

A nonlinear effect, RF (radio frequency) condensation, can be used to facilitate RF stabilization of magnetic islands. Previously studied semi-analytical models for RF condensation are suited mainly for broad deposition profiles and are unable to handle power depositions that are localized in the interior of a magnetic island. The power deposition profile of electron cyclotron waves, in particular, can be highly localized near a cyclotron resonance. Here, a model is developed that can handle both localized profiles and broad profiles. This allows a comparison of RF condensation for narrow vs broad deposition profiles, and it allows a study of the dependence of RF condensation of localized deposition profiles on key parameters.

Upstream Plasma Waves and Downstream Magnetic Reconnection at a Reforming Quasi-parallel Shock

With the help of a two-dimensional particle-in-cell simulation model, we investigate the long-time evolution (near 100Ωi0−1 , where Ω i0 is the ion gyrofrequency in the upstream) of a quasi-parallel shock. Some of the upstream ions are reflected by the shock front, and their interactions with the incident ions excite low-frequency magnetosonic waves in the upstream. Detailed analyses have shown that the dominant wave mode is caused by the resonant ion–ion beam instability, and the wavelength can reach tens of ion inertial lengths. Although these plasma waves are directed toward the upstream in the upstream plasma frame, they are brought by the incident plasma flow toward the shock front, and their amplitude is enhanced during the approaching. The interaction of the upstream plasma waves with the shock leads to the cyclic reformation of the shock front, and the reformation period is slightly larger than 10 Ωi0−1 . When crossing the shock front, these large-amplitude plasma waves are compressed and evolve into current sheets in the transition region of the shock. At last, magnetic reconnection occurs in these current sheets, accompanying the generation of magnetic islands. Simultaneously, there still exist plasma waves of another kind, with the wavelength of several ion inertial lengths in the ramp of the shock, which are excited by the nonresonant ion–ion beam instability. The current sheets in the transition region are distorted and broken into several segments when the plasma waves of this kind are transmitted into the downstream, where magnetic reconnection and the generated islands have a much smaller size. No obvious ion flow can be observed around some X-lines produced in the magnetic reconnection, and this implies that electron-only reconnection may occur.

Coherent structures of beam-driven whistler mode in the presence of magnetic islands in the magnetopause

The Magnetospheric Multiscale Mission (MMS) has perceived whistler wave generation, coherent structures, and related turbulence close to the magnetopause reconnection zones. The current research examines coherent structure of whistler wave driven by an intense electron beam at the magnetopause’s magnetic reconnection sites as well as by the dynamic growth of magnetic islands. A nonlinear model of high-frequency whistler wave and low-frequency magnetosonic wave has been developed by using the two-fluid approximation. Nonlinear dynamics of 3D whistler wave and magnetosonic wave have been solved by the pseudo spectral method along with the predictor-corrector method and finite difference method. The simulation’s outcomes demonstrate the temporal and spatial development of the whistler localized structures and current sheets as a witness to the turbulence’s existence. Moreover, the turbulent power spectra have been investigated. The formation of the thermal tail of energetic electrons has been studied using the power-law scaling of turbulence development. We determined the scale sizes of current sheets and localized structures using a semi-analytic model and showed that these scale sizes rely on the power of whistler wave. We predict that the acceleration of the energetic electrons and heating in the Magnetopause may be caused by whistler wave.

Intermittent Electron-only Reconnection at Lunar Mini-magnetospheres

Lunar crustal magnetic anomalies (LCMA) are sub-ion-gyroradius structures that have been shown to stand off the solar wind (SW) plasma from the Moon’s surface, forming shock-like discontinuities and reflecting incident SW protons. In this Letter, the results of high-resolution, 2D fully kinetic simulations show a bursty electron-only magnetic reconnection in the SW-LCMA interaction region, characterized by the quasiperiodic formation and ejection of magnetic islands and strong parallel electron flows along the X-point separator lines. The islands are observed to modify the magnetic pressure pileup and Hall electric field above the LCMA, leading to sharp increases in reflected protons that drive electromagnetic fluctuations downstream and short distances upstream in the SW.

Toxoplasmosis epidemic in a population of urbanised allied rock-wallabies (Petrogale assimilis) on Magnetic Island (Yunbenun), North Queensland.

A mortality event involving 23 allied rock-wallabies (Petrogale assimilis) displaying neurological signs and sudden death occurred in late April to May 2021 in a suburban residential area directly adjacent to Magnetic Island National Park, on Magnetic Island (Yunbenun), North Queensland, Australia. Three allied rock-wallabies were submitted for necropsy, and in all three cases, the cause of death was disseminated toxoplasmosis. This mortality event was unusual because only a small, localised population of native wallabies inhabiting a periurban area on a tropical island in the Great Barrier Reef World Heritage Area were affected. A disease investigation determined the outbreak was likely linked to the presence of free-ranging feral and domesticated cats inhabiting the area. There were no significant deaths of other wallabies or wildlife in the same or other parts of Magnetic Island (Yunbenun) at the time of the outbreak. This is the first reported case of toxoplasmosis in allied rock-wallabies (Petrogale assimilis), and this investigation highlights the importance of protecting native wildlife species from an infectious and potentially fatal parasitic disease.

Perturbed Ion Temperature and Toroidal Flow Profile Measurements in Rotating Neoclassical Tearing Mode Magnetic Islands.

The perturbed ion temperature and toroidal flow were measured in rotating neoclassical tearing modes (NTM) in a tokamak for the first time. These toroidally and radially resolved profiles were obtained by impurity ion spectroscopy in a 2,1 NTM in DIII-D. In agreement with drift-kinetic simulations, the electron temperature profile is flat, while the ion temperature gradient is restored across the magnetic island O point in the presence of fast ions; the perturbed flow has minima in the O points and maxima at the X points. These measurements provide the first confirmation of the theoretically expected ion temperature and flow response to a magnetic island needed to predict the NTM onset threshold scaling for ITER and other future tokamaks.

Numerical MHD simulations of solar flares and their associated small-scale structures

ABSTRACT Using numerical simulations, we study the formation and dynamics of solar flares in a local region of the solar atmosphere. The magnetohydrodynamics (MHD) equations describe the dynamic evolution of flares, including space-dependent and anomalous magnetic resistivity and highly anisotropic thermal conduction on a 2.5 D slice. We adopt an initial solar atmospheric model in magnetohydrostatic equilibrium, with a magnetic configuration consisting of a vertical current sheet, which helps trigger the magnetic reconnection process. Specifically, we study three scenarios, two with only resistivity and the third with resistivity plus thermal conduction. The main results of the numerical simulations show differences in the global morphology of the flares, including the post-flare loops and the current sheet in three cases. In particular, localized resistivity produces more substructure around the post-flare loops that could be related to the Ritchmyer–Meshkov Instability (RMI). Furthermore, in the scenario of anomalous resistivity, we identify the formation of a plasmoid and a jet at coronal heights. On the other hand, in the scenario with resistivity plus thermal conduction, the post-flare loops are smooth, and no apparent substructures develop. Besides, in the z-component of the current density for the Res + TC case, we observe the development of multiple magnetic islands generated due to the Tearing instability in the non-linear regime.

Unveiling non-flat profiles within magnetic islands in tokamaks

The presence of non-flat profiles on magnetic island is studied for the first time through gyrokinetic simulations alongside a simplified Lagrangian model. We have identified that inside a magnetic island, the non-flatness of density and temperature profiles is controlled by a dimensionless parameter α≡w*ŝϵ/qρ*, which is a function of normalized island width w*=w/a0, magnetic shear ŝ, inverse aspect ratio ϵ=a0/R, safety factor q, and normalized gyroradius ρ*=ρ/a0. The gyroradius ρ* dependence of the control parameter α leads to a species-selective transition of profiles from flat to concave only for electrons having high α∼O(1). The finding elucidates that electron profiles tend to increasingly deviate from the flat state for a larger magnetic island, in contrast to the conventional wisdom.

Generalized kinetic equation for tokamak plasma equilibrium distribution function

A generalized kinetic equation for the equilibrium distribution function in a finite beta, arbitrary tokamak plasma is derived. The equation is correct to second order in ρ/L (ρ is the particle Larmor radius and L is the system size). Resolving the finite Larmor radius length scales with no restriction on the ratio of poloidal to total equilibrium magnetic field, Bϑ/B, it generalizes the drift kinetic theory of Hazeltine [Phys. Plasmas 15, 77 (1973)] to the limit of Bϑ/B∼1 (e.g., to ensure validity for spherical tokamaks). Two cases are considered. The first provides the equilibrium distribution function, consistent with the generalized gyrokinetic formalism of Dudkovskaia et al. [Plasma Phys. Controlled Fusion 65, 045010 (2023)], derived specifically to capture neoclassical equilibrium currents in the gyrokinetic stability analyses in strong gradient regions. The second assumes short length scales in the direction perpendicular to the magnetic field, which can occur as a result of small coherent magnetic structures in the plasma, such as neoclassical tearing mode magnetic islands close to threshold. This then extends the drift island equations of Dudkovskaia et al. [Nucl. Fusion 63, 016020 (2023)] for the plasma response to magnetic islands to a spherical tokamak plasma configuration. Resolving ρ∼ρϑ (or Bϑ∼B), where ρϑ is the particle poloidal Larmor radius, is also expected to influence calculations of the magnetic island propagation frequency and the associated contributions to the island onset conditions.

Effects of electron viscosity on resonant layer responses to non-axisymmetric magnetic perturbations

The resonant field penetration to magnetic islands is the central MHD mechanism of non-axisymmetric plasma responses in a tokamak such as disruptive locking or favorable ELM stabilization. The resonant field penetration can be induced by any non-ideal processes as manifested in the delicate balance under the generalized Ohm's law. Here, we show that the viscous effects by electrons are not ignorable in the field penetration unlike previous presumption, even if the electron viscosity is as small as the square root of its mass compared to the ions. It is clear that its effects become only bigger if the electron viscosity becomes anomalously large. The work strictly follows the three-field model in the linear regime targeting the prediction of the onset of the field penetration and successfully extending it with electron viscosity and identifying new regimes. The results also indicate that the error field thresholds become more strongly dependent on plasma density than ones predicted in the linear regimes without the electron viscosity, which is consistent with experimental observations and thus a significant implication.

Role of nonlinear structures and associated turbulence generation dayside magnetosphere reconnection sites

A numerical simulation is implemented to investigate the role of whistler wave in the presence of the ponderomotive force-driven density modification and the magnetic field perturbation due to preexisting magnetic islands. The dynamical equation governing the whistler waves is derived and solved using numerical methods to assess their role. The simulation produced results that revealed the nonlinear structures, confirming the system's turbulent characteristics. Moreover, we utilized a semi-analytical model, applying the paraxial approximation, to estimate the scale size of the current sheet, which was found to be between 2λe and λe. This measurement provides evidence of the presence of a superthin electron current sheet within reconnection sites. Qualitative consistency is established by comparing the simulation results with findings reported in the literature.

Effects of diamagnetic drift on nonlinear interaction between multi-helicity neoclassical tearing modes

A numerical study of the diamagnetic drift effect on the nonlinear interaction between multi-helicity neoclassical tearing modes (NTMs) is carried out using a set of four-field equations including two-fluid effects. The results show that, in contrast to the single-fluid case, 5/3 NTM cannot be completely suppressed by 3/2 NTM with diamagnetic drift flow. Both modes exhibit oscillation and coexist in the saturated phase. To better understand the effect of the diamagnetic drift flow on multiple-helicity NTMs, the influence of typical relevant parameters is investigated. It is found that the average saturated magnetic island width increases with increasing bootstrap current fraction f b but decreases with the ion skin depth δ. In addition, as the ratio of parallel to perpendicular transport coefficients χ ∥/χ ⊥ increases, the average saturated magnetic island widths of the 3/2 and 5/3 NTMs increase. The underlying mechanisms behind these observations are discussed in detail.

Localization of beam generated whistler wave and turbulence generation in reconnection region of magnetopause

Whistler waves have been studied for many years in relation to turbulence and particle heating, and observations show that they are crucial to magnetic reconnection. Recent research has revealed a close relationship between magnetic reconnection and turbulence. The current work investigates the whistler turbulence caused by the energetic electron beam in the magnetic reconnection sites of magnetopause and also due to dynamic evolution of magnetic islands. For this, we develop a model based upon the two-fluid approximation to study whistler dynamics, propagating in the medium with the pre-existing chain of magnetic islands and under the influence of background density perturbation originating from ponderomotive nonlinearity of wave. Dynamics of nonlinear whistler have been solved with pseudo-spectral approach and a finite difference method with a modified predictor–corrector method and a Runge Kutta method for the semianalytical model. In the current research, we study how the nonlinear whistler wave contributes to the significant space phenomenon, i.e., turbulence, localization, and magnetic reconnection. We have also investigated the formation of a current sheet in a magnetopause region of the order of few-electron inertial length. We analyzed the power spectrum at the magnetopause when the system reached a quasi-steady condition. Our new approach to study whistler turbulence by an energetic electron beam at the magnetic reconnection sites has extensive applications to space plasmas, shedding a new light on the study of magnetic reconnection in nature.

Effect of nonlinearity, magnetic islands on turbulence and observation of electron energization, temperature anisotropy at Earth’s magnetopause (magnetosphere)

Space missions’ observations have shown that waves such as lower hybrid waves (LHWs), whistler waves, and kinetic Alfven waves play a vital role in magnetic reconnection, turbulence, and particle acceleration. This paper studies the effect of nonlinearity and the magnetic islands on lower hybrid turbulence and current sheets in Earth’s magnetopause region. The evolution of electromagnetic LHW has been studied with numerical model using pseudo-spectral method for spatial integration and finite difference method with modified predictor-corrector approach for temporal integration. We have considered both ion and electron dynamics and included electromagnetic and warm plasma effects in our model. The study outcomes reveal that both the nonlinear effects and magnetic islands are responsible for the evolution of LHWs and current sheets to a chaotic and turbulent state. We have also used the semi-analytical model to elaborate on the physics behind the localization. Finally, the nonlinear model with field perturbations (magnetic islands) is further used to elaborate on the electron energization and temperature anisotropy near reconnection regions. We have also discussed the relevance of model predictions in the context of the MMS mission observations at Earth’s magnetopause.

Enhancement of ECCD by the current condensation effect for stabilizing large magnetic islands caused by neoclassical tearing modes in tokamak plasmas

The radio frequency current condensation effect reported in Reiman and Fisch (2018 Phys. Rev. Lett. 121 225001) is modeled in the nonlinear resistive magnetohydrodynamic code. A series of numerical investigations have been performed to investigate the enhancement of electron cyclotron current drive (ECCD) by the current condensation effect during the control of neoclassical tearing mode (NTM) in tokamak plasmas. In the numerical model, both the parallel transport and the perpendicular transport of electron temperature are considered. The EC driven current and driven perturbed electron temperature can nonlinearly evolve within the given magnetic configuration and eventually reach saturation states. The input power threshold of ECCD and the fold bifurcation phenomenon are numerically verified via nonlinear simulations. The numerical results show good agreements with the analytical results. Moreover, spatial distributions of EC current for the two solutions at different condensed level are displayed. The control effectiveness of ECCD for large NTM islands has been evaluated while considering the current condensation effect. While taking into account current condensation effect, for a sufficiently large input power, a larger island can be more effectively stabilized than a smaller one, which suggests a reassessment of the previous idea that the ECCD should always be turned on as early as possible. The potential physics mechanism behind the ECCD control have all been discussed in detail. Furthermore, the condensation effect is found to have favorable effects on the radial misalignment of ECCD. In the consideration of the situation for extremely localized control needs, a highly peaked heating profile is adopted to verify that the fold bifurcation phenomenon still exists and the current condensation effect can still take effect in this extreme condition.

Impact of microwave beam scattering by density fluctuations on the electron–cyclotron power deposition profile in tokamaks

Electron–cyclotron waves are a tool commonly used in tokamaks, in particular to drive current. Their ability to drive current in a very localized manner renders them an optimal tool for MHD mode mitigation. However, such applications require high accuracy and good control of the power deposition location to efficiently target the magnetic islands. It has been indirectly observed that the suprathermal electron distribution, resulting from the wave absorption, is broader than what is expected from experimentally-constrained forward drift-kinetic modeling. The present paper explores the possibility that beam scattering through the turbulent edge of the plasma may explain this observed discrepancy. In particular, full-wave studies exhibit three beam broadening regimes, from superdiffusive to diffusive, with an intermediate regime characterized by a Lorentzian beam profile with a slightly increased full-width at half maximum with respect to the quiet plasma case. In the tokamak à configuration variable, dedicated plasma scenarios have been developed to test this hypothesis. A realistic worst-case fluctuation scenario falls into this intermediate beam broadening regime. By comparing the experimental hard x-ray emission from suprathermal electron Bremmstrahlung with the emission calculated by coupling a full-wave model to a Fokker–Planck solver, it is shown that, in the tested cases, the beam broadening is not sufficient to explain the aforementioned discrepancy between simulation and experiment and that another mechanism must play the main role in broadening the suprathermal electron distribution.

Effect of 2/1 tearing mode on radiation asymmetry during disruptions on J-TEXT

Plasma disruption in the tokamak is a huge threat to the safe operation of ITER, which will cause serious damage to the tokamak device. At present, massive impurity injection is the main disruption mitigation method, include massive gas injection (MGI) and shattered pellet injection (SPI) which can dissipate most of the internal thermal energy via radiation. However, radiation asymmetry will still cause damage to the device when it results the concentration of radiation in a small area. It has been observed on many devices that n = 1 mode has a strong influence on radiation asymmetry. In this paper, the experimental results of radiation asymmetry on J-TEXT tokamak have been introduced, and there is a close to sine function relationship between toroidal peaking factor (TPF) and the phase of 2/1 tearing mode. It is found that when the injection position is aligned with the X-point of the pre-existing magnetic island, TPF has the maximum value, and while the injection position is between the X-point and the O-point, TPF has the minimum value. This will be helpful for the study of disruption mitigation of future large devices.

Anisotropic heat diffusion in stochastic magnetic fields

AbstractThe magnetic topology is a critical issue in fusion plasma research. An example is the Resonant Magnetic Perturbation (RMP), which controls the edge transport in tokamaks. However, the physics of how the RMP affects edge transport is not clear. One reason is the transport process on the stochastic magnetic field is poorly understood. This study examines anisotropic heat diffusion numerically to understand heat transport in stochastic magnetic fields. We developed a numerical model of an anisotropic temperature diffusion model, where the significant deviation of the parallel and perpendicular thermal conductivity exists. We applied this implementation to the realistic stellarator geometry with the stochastic magnetic field in the edge. The smooth temperature profile is obtained for the large perpendicular diffusion, although the magnetic field is stochastic. However, for another case of significant parallel diffusion, the small flattening of the temperature on the magnetic island in the stochastic region appears. That result suggests that the stochastic magnetic field can keep the finite temperature gradient if the connection length of the magnetic field line in the stochastic region is sufficiently long.

First observation of beta-induced Alfvén eigenmode inside the edge magnetic island on the J-TEXT tokamak

Experiments conducted on the J-TEXT tokamak have provided the first evidence that the Beta-induced Alfvén Eigenmode (BAE) is localized inside the isolated helical flux tube of its edge m/n = 3/1 magnetic island. The observations show that the BAE forms a standing wave inside the magnetic island, with its nodes located at the X- and O-points of the magnetic island. When the island is cut open by contact with the limiter plates, the BAE is found to remain inside the remnant closed island in the Scrape-Off Layer, but its amplitude decreases as the width of the remnant island becomes smaller.

Gyrokinetic simulation of magnetic-island-induced electric potential vortex mode

Ion temperature gradient (ITG)-driven turbulence with embedded static magnetic islands is simulated by utilizing a gyrokinetic theory-based global turbulence transport code (GKNET) in this work. Different from the traditional equilibrium circular magnetic-surface average (EMSA) method, an advanced algorithm that calculates the perturbed magnetic-surface average (PMSA) of the electric potential has been developed to precisely deal with the zonal flow component in a non-circular magnetic surface perturbed by magnetic islands. Simulations show that the electric potential vortex structure inside islands induced by the magnetic islands is usually of odd parity when using the EMSA method. It is found that the odd symmetry vortex can transfer into an even one after a steep zonal flow gradient, i.e. the flow shear has been built in the vicinity of the magnetic islands by adopting the PMSA algorithm. The phase of the potential vortex in the poloidal cross section is coupled with the zonal flow shear. Such an electric potential vortex mode may be of essential importance in wide topics, such as the turbulence spreading across magnetic islands, neoclassical tearing mode physics, and also the interaction dynamics between the micro-turbulence and MHD activities.