Subsolar magnetopause observation and kinetic simulation of a tripolar guide magnetic field perturbation consistent with a magnetic island
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
Abstract not provided
Magnetic reconnection is a primary mechanism for particle energization in space and astrophysical plasmas. By carrying out two-dimensional (2D) fully kinetic simulations, we study particle acceleration during magnetic reconnection in plasmas with different plasma β (the ratio between the thermal pressure and the magnetic pressure). For the high-β cases, we do not observe significant particle acceleration. In the low-β regime (β<0.1), we find that reconnection is efficient at energizing both electrons and ions. While the distribution of accelerated particles integrated over the whole simulation box appears highly non-thermal, it is actually the superposition of a series of distributions in different sectors of a 2D magnetic island. Each of those distributions has only a small non-thermal component compared with its thermal core. By tracking a large number of particles, we show that particles get energized in X-line regions, contracting magnetic islands, and magnetic island coalescence regions. We obtain the particle energization rate j⋅E by averaging over particle drift motions and find that it agrees well with the particle kinetic energy change. We quantify the contribution of curvature drift, gradient drift, polarization drift, magnetization, non-gyrotropic effect, and parallel electric field in different acceleration regions. We find that the major energization is due to particle curvature drift along the motional electric field. The other particle motions contribute less but may become important in different acceleration regions. The highly efficient particle energization in low-β plasmas may help us understand the strong particle energization in solar flares and accretion disk coronae.
Abstract The safe operation of most tokamaks, especially the large ones, relies on the feedback control of vertical displacement events (VDEs). However, most of these feedback control systems are based on axisymmetric VDE models. In this study, we use NIMROD simulations to study the role of non-axisymmetric perturbations in free drift vertical displacement events on EAST. The high- n modes in the non-axisymmetric VDE grow first, which drives the formation of high- n magnetic island chains. Subsequently, the magnetic island chains grow and overlap with each other, leading to the destruction of the magnetic flux surface, which induces a minor disruption and accelerates the start of the following major disruption. The magnetic island and the stochastic magnetic field allow the toroidally non-axisymmetric poloidal plasma current to jet towards the hoop force direction, forming finger-like and filamentary structures. Such a plasma current non-axisymmetry strongly depends on the anisotropy in the thermal transport coefficients.
Shock waves, as shown by simulations and observations, can generate high levels of downstream vortical turbulence, including magnetic islands. We consider a combination of diffusive shock acceleration (DSA) and downstream magnetic-island-reconnection-related processes as an energization mechanism for charged particles. Observations of electron and ion distributions downstream of interplanetary shocks and the heliospheric termination shock (HTS) are frequently inconsistent with the predictions of classical DSA. We utilize a recently developed transport theory for charged particles propagating diffusively in a turbulent region filled with contracting and reconnecting plasmoids and small-scale current sheets. Particle energization associated with the anti-reconnection electric field, a consequence of magnetic island merging, and magnetic island contraction, are considered. For the former only, we find that (i) the spectrum is a hard power law in particle speed, and (ii) the downstream solution is constant. For downstream plasmoid contraction only, (i) the accelerated spectrum is a hard power law in particle speed; (ii) the particle intensity for a given energy peaks downstream of the shock, and the distance to the peak location increases with increasing particle energy, and (iii) the particle intensity amplification for a particular particle energy, f(x,c/c{sub 0})/f(0,c/c{sub 0}), is not 1, as predicted by DSA, but increases with increasing particle energy. The general solution combines both the reconnection-induced electric field and plasmoid contraction. The observed energetic particle intensity profile observed by Voyager 2 downstream of the HTS appears to support a particle acceleration mechanism that combines both DSA and magnetic-island-reconnection-related processes.
An emerging paradigm for the dissipation of magnetic turbulence in the supersonic solar wind is via localized quasi-2D small-scale magnetic island reconnection processes. An advection-diffusion transport equation for a nearly isotropic particle distribution describes particle transport and energization in a region of interacting magnetic islands [1; 2]. The dominant charged particle energization processes are 1) the electric field induced by quasi-2D magnetic island merging, and 2) magnetic island contraction. The acceleration of charged particles in a “sea of magnetic islands” in a super-Alfvénic flow, and the energization of particles by combined diffusive shock acceleration (DSA) and downstream magnetic island reconnection processes are discussed.
The international ITER tokamak has the objective of demonstrating the scientific feasibility of magnetic confinement fusion as a source of energy. A concern towards the achievement of this goal is represented by major disruptions: complete losses of confinement often initiated by a non-rotating ('locked') magnetic island created by magnetic reconnection. During disruptions, energy and particles accumulated in the plasma volume over many seconds are lost in a few milliseconds and released on the plasma-facing materials. In addition, multi-MA level currents flowing in the tokamak plasma for its sustainment and confinement are lost, also in milliseconds, thus terminating the plasma discharge and causing electromagnetic stresses that, if unmitigated, could lead to excessive device wear. Moreover it is shown that magnetic perturbations can be used to avoid disruptions by "guiding" the magnetic island to lock in a position where it is accessible to millimetre wave beams that fully stabilize it.
The aim of this project with international collaborations (PI: J.-K. Park of Princeton Plasma Physics Laboratory) is to develop a unified physics basis and predictive capability for the control of edge-localized modes (ELMs) with optimized non-axisymmetric (3D) fields, by leveraging the unique research capabilities of international KSTAR tokamak in Korea. At the University of California, Irvine, we use a fusion simulation code GTC to study the effects of the 3D magnetic fields (including both magnetic islands and stochastic magnetic field lines) on microturbulence. We focus on the effects of 3D equilibrium on turbulent and neoclassical transport that leads to ELM suppression. During the project period, the UCI team of this collaboration has made significant progress in understanding turbulent transport in 3D geometry including both tokamaks with RMP and stellarators. Four UCI postdoctoral researchers, Drs. Javier Nicolau, Gyungjin Choi, Pengfei Liu, and Xishuo Wei were partially supported by this project and collaborating with the entire project team. We have published 4 journal papers and gave several invited talks. This project advance physics understand of 3D effects on microturbulence, specially effects of magnetic islands. It will contribute to the DOE SciDAC ISEP project in the area of GTC development for general 3D geometry capability, rigorous verification and validation, and readiness for the emerging exascale computing. Through collaborations with computational scientists, the project helps to advance high performance computing in the area of GPU acceleration, linear solver, and data management.
The interactions between magnetic islands are considered to play an important role in electron acceleration during magnetic reconnection. In this paper, two-dimensional particle-in-cell simulations are performed to study electron acceleration during multiple X line reconnection with a guide field. Because the electrons remain almost magnetized, we can analyze the contributions of the parallel electric field, Fermi, and betatron mechanisms to electron acceleration during the evolution of magnetic reconnection through comparison with a guide-center theory. The results show that with the magnetic reconnection proceeding, two magnetic islands are formed in the simulation domain. Next, the electrons are accelerated by both the parallel electric field in the vicinity of the X lines and the Fermi mechanism due to the contraction of the two magnetic islands. Then, the two magnetic islands begin to merge into one, and, in such a process, the electrons can be accelerated by both the parallel electric field and betatron mechanisms. During the betatron acceleration, the electrons are locally accelerated in the regions where the magnetic field is piled up by the high-speed flow from the X line. At last, when the coalescence of the two islands into one big island finishes, the electrons can be further accelerated by the Fermi mechanism because of the contraction of the big island. With the increase of the guide field, the contributions of the Fermi and betatron mechanisms to electron acceleration become less and less important. When the guide field is sufficiently large, the contributions of the Fermi and betatron mechanisms are almost negligible.
In this work, we compare gyrokinetic (GK) with fully kinetic Particle-in-Cell (PIC) simulations of magnetic reconnection in the limit of strong guide field. In particular, we analyze the limits of applicability of the GK plasma model compared to a fully kinetic description of force free current sheets for finite guide fields (b{sub g}). Here, we report the first part of an extended comparison, focusing on the macroscopic effects of the electron flows. For a low beta plasma (β{sub i} = 0.01), it is shown that both plasma models develop magnetic reconnection with similar features in the secondary magnetic islands if a sufficiently high guide field (b{sub g} ≳ 30) is imposed in the kinetic PIC simulations. Outside of these regions, in the separatrices close to the X points, the convergence between both plasma descriptions is less restrictive (b{sub g} ≳ 5). Kinetic PIC simulations using guide fields b{sub g} ≲ 30 reveal secondary magnetic islands with a core magnetic field and less energetic flows inside of them in comparison to the GK or kinetic PIC runs with stronger guide fields. We find that these processes are mostly due to an initial shear flow absent in the GK initialization and negligible in the kinetic PIC high guide field regime, in addition to fast outflows on the order of the ion thermal speed that violate the GK ordering. Since secondary magnetic islands appear after the reconnection peak time, a kinetic PIC/GK comparison is more accurate in the linear phase of magnetic reconnection. For a high beta plasma (β{sub i} = 1.0) where reconnection rates and fluctuations levels are reduced, similar processes happen in the secondary magnetic islands in the fully kinetic description, but requiring much lower guide fields (b{sub g} ≲ 3)
We report a two-stage crash process in edge localized mode (ELM) driven by resistive drift-ballooning modes (RDBMs) numerically simulated in a full annular torus domain with a scale-separated four-field reduced MHD (RMHD) model using the BOUT++ framework. In the early nonlinear phase, the small first crash is triggered by linearly unstable RDBMs, and m/n=2/1 magnetic islands are nonlinearly excited by nonlinear coupling of RDBMs as well as their higher harmonics. Here, m is the poloidal mode number, n is the toroidal mode number, the q = 2 rational surface exists near the pressure gradient peak, and q is the safety factor. Simultaneously, middle-n RDBM turbulence develops but is poloidally localized around X-points of the magnetic islands, leading to the small energy loss. The second large crash occurs in the late nonlinear phase. Higher harmonics of m/n=2/1 magnetic islands well develop around the q = 2 surface via nonlinear coupling and make the magnetic field stochastic by magnetic island overlapping. Turbulence heat transport develops at X-points of higher harmonics of m/n=2/1 magnetic islands, resulting in the turbulence spreading in the poloidal direction. The large second crash is triggered when the turbulence covers the whole poloidal region so that the magnetic island generation and magnetic field stochastization before the large crash can be interpreted as ELM precursors. It is concluded that the ELM trigger is attributed to the turbulent spreading in the poloidal direction in synchronization with the magnetic field stochastization and the crash is driven by E × B convection rather than the conventional Rechester–Rosenbluth anomalous electron heat transport.
A new method to control the tearing mode (TM) in tokamaks has been proposed [Q. Hu and Q. Yu, Nucl. Fusion 56, 034001 (5pp.) (2016)], according to which, the external resonant magnetic perturbation needs to be applied in certain magnetic island phase regions. Therefore, it is very important to identify the helical phase of magnetic islands in real time. The TM in tokamak plasmas is normally rotating and carries magnetic oscillations, which are known as Mirnov oscillations and can be detected by Mirnov probes. When the O-point or X-point of the magnetic island passes through the probe, the signal will experience a zero-crossing. A poloidal Mirnov probe array and a corresponding island phase identification method are presented. A field-programmable gate array is used to provide the magnetic island helical phase in real time by using multichannel zero crossing detection. This system has been developed on the J-TEXT tokamak and works well. This paper introduces the establishment of the fast magnetic island phase identifying system.
We perform 2D resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible to the formation of blobs that are ejected in the jet. Our simulation with low plasma β (Case I) shows that magnetic islands form easily and propagate upward in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 Å and 211 Å) and the cool channel (304 Å) in observations by the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory . The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s{sup −1}. These parameters are similar to the properties of blobs observed in extreme-ultraviolet (EUV) jets by AIA. The Kelvin–Helmholtz instability develops in our simulation with moderately high plasma β (Case II) and leads to the formation of bright vortex-like blobs above the multiple high magnetosonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magnetosonic Mach number regions disappear. In the lower plasma β case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current-sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin–Helmholtz instability along the jet are both possible causes of the formation of blobs observed at EUV wavelengths.
Abstract According to recent DIII-D experiments (Logan et al 2024 Nucl. Fusion 64 014003), injecting edge localized electron cyclotron current drive (ECCD) in the counter-plasma-current (counter- I p ) direction reduces the n = 3 resonant magnetic perturbation (RMP) current threshold for edge-localized mode (ELM) suppression, while co- I p ECCD during the suppressed ELM phase causes a back transition to ELMing. This paper presents nonlinear two-fluid simulations on the ECCD manipulation of edge magnetic islands induced by RMP using the TM1 code. In the presence of a magnetic island chain at the pedestal-top, co- I p ECCD is found to decrease the island width and restore the initially degraded pedestal pressure when its radial deposition location is close to the rational surface of the island. With a sufficiently strong co- I p ECCD current, the RMP-driven magnetic island can be healed, and the pedestal pressure fully recovers to its initial ELMing state. On the contrary, counter- I p ECCD is found to increase the island width and further reduce the pedestal pressure to levels significantly below the peeling-ballooning-mode limited height, leading to even stationary ELM suppression. These simulations align with the results from DIII-D experiments. However, when multiple magnetic island chains are present at the pedestal-top, the ECCD current experiences substantial broadening, and its effects on the island width and pedestal pressure become negligible. Further simulations reveal that counter- I p ECCD enhances RMP penetration by lowering the penetration threshold, with the degree of reduction proportional to the amplitude of ECCD current. For the ∼1 MW ECCD in DIII-D, the predicted decrease in the RMP penetration threshold for ELM suppression is approximately 20%, consistent with experimental observations. These simulations indicate that edge-localized ECCD can be used to either facilitate RMP-driven ELM suppression or optimize the confinement degradation.
Disruption affected by locked modes (LMs) is one of the most critical issues to be solved for tokamak fusion reactors. This paper aimed to understand the plasma response to external rotating resonant magnetic perturbations (RMPs), which are applied at a few kilohertz for controlling the pre-existing LMs. In the J-TEXT tokamak, the plasma response to a rotating RMP showed a feature of travelling wave, if no tearing mode (TM) existed. With a LM, the application of a rotating RMP led to unique features of plasma responses, which could be roughly described by the formula δb θ (θ, t)=2Asin(mθ–mθ O )cos(2πf RRMP t), where δb θ , A, m, θ O and f RRMP are the perturbed poloidal magnetic field, the amplitude, the poloidal mode number, the poloidal location of the TM's O-point and the frequency of the rotating RMP. These features are analogous to that of a standing wave in the TM rest frame, with the nodes locating at the O-/X- points. The perturbations of electron temperature, T e , due to plasma responses were zero or minimal inside the magnetic island or outside the magnetic island but at the same poloidal positions as the O-/X- points. Both δb θ and the T e perturbations increased linearly with the rotating RMP. A phenomenological model is presented that the rotating RMP drives the phase oscillation of the TM by applying a periodical electromagnetic torque. Then the phase oscillation of TM, combined with the poloidal gradient of the plasma parameters, induces the plasma response measured experimentally. The nonlinear numerical modelling confirmed the forced phase oscillation of the TM, and qualitatively verified the amplitude and phase dependences among the plasma response, the phase oscillation of the TM and the rotating RMP, which have been predicted by the model.
In the Wendelstein 7-X (W7-X) stellarator, the vacuum rotational transform, ι, has a flat radial profile and does not cross any major rational resonance. Nevertheless, during plasma operation the ι–profile can be strongly modified by electron cyclotron current drive in such a way that the resulting ι-profile passes through low-order rational values, and this can trigger magnetohydrodynamic (MHD) events. Indeed, W7-X plasmas are sometimes subject to repetitive collapses of core confinement, which can be observed regardless of the direction in which the EC current is driven. Even though the origin of these MHD instabilities is under investigation, the crashes may be connected to the formation of magnetic islands and magnetic reconnection. In the present work, we try to shed light on the dynamics of different events happening during the course of sawtooth cycles in W7-X by proposing a model that combines a slow current diffusion with a recipe for fast relaxation that conserves the corresponding helical flux (Kadomtsev 1975 Fiz. Plazmy 1 710–15). We also propose a simple model based on Taylor relaxation (Taylor 1974 Phys. Rev. Lett. 33 1139), (Taylor 1986 Rev. Mod. Phys. 58 741) to predict the nonlinear redistribution of plasma current caused by the largest of the observed events.
Numerical simulations of runaway electrons (REs) with a particular emphasis on orbit dependent effects in 3-D magnetic fields are presented. The simulations were performed using the recently developed Kinetic Orbit Runaway electron Code (KORC) that computes the full-orbit relativistic dynamics in prescribed electric and magnetic fields including radiation damping and collisions. The two main problems of interest are synchrotron radiation and impurity-based RE dissipation. Synchrotron radiation is studied in axisymmetric fields and in 3-D magnetic configurations exhibiting magnetic islands and stochasticity. For passing particles in axisymmetric fields, neglecting orbit effects might underestimate or overestimate the total radiation power depending on the direction of the radial shift of the drift orbits. For trapped particles, the spatial distribution of synchrotron radiation exhibits localized “hot” spots at the tips of the banana orbits. In general, the radiation power per particle for trapped particles is higher than the power emitted by passing particles. The spatial distribution of synchrotron radiation in stochastic magnetic fields, obtained using the MHD code NIMROD, is strongly influenced by the presence of magnetic islands. 3-D magnetic fields also introduce a toroidal dependence on the synchrotron spectra, and neglecting orbit effects underestimates the total radiation power. In the presence of magnetic islands, the radiation damping of trapped particles is larger than the radiation damping of passing particles. Results modeling synchrotron emission by RE in DIII-D quiescent plasmas are also presented. The computation uses EFIT reconstructed magnetic fields and RE energy distributions fitted to the experimental measurements. Qualitative agreement is observed between the numerical simulations and the experiments for simplified RE pitch angle distributions. However, it is noted that to achieve quantitative agreement, it is necessary to use pitch angle distributions that depart from simplified 2-D Fokker-Planck equilibria. Finally, using the guiding center orbit model (KORC-GC), a preliminary study of pellet mitigated discharges in DIII-D is presented. The dependence of RE energy decay and current dissipation on initial energy and ionization levels of neon impurities is studied. Furthermore, the computed decay rates are within the range of experimental observations.