Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETIC ISLANDS”

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.

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Plasma Energization in Colliding Magnetic Flux Ropes

Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized with two magnetic islands embedded in a reconnecting current sheet. The two islands collide and eventually merge into a single island. Particles are accelerated during this process as the magnetic energy is released and converted to the plasma energy, including bulk kinetic energy increase by the ideal electric field, and thermal energy increase by the fluid compression and the nonideal electric field. We find that contributions from these different energization mechanisms are all important and comparable with each other. Fluid shear and a nongyrotropic pressure tensor also contribute to the energy conversion process. For simulations with different box sizes ranging from $L_x$ ~ 25–100$d_i$ and ion-to-electron mass ratios $m_i /m_e$ = 25,100, and 400, we find that the general evolution is qualitatively the same for all runs, and the energization depends only weakly on either the system size or the mass ratio. The results may help us understand plasma energization in solar and heliospheric environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Electron thermal confinement in a partially stochastic magnetic structure

Using a high-repetition-rate Thomson scattering diagnostic, we observe a peak in electron temperature $T_e$ coinciding with the location of a large magnetic island in the Madison Symmetric Torus. Magnetohydrodynamic modeling of this quasi-single helicity plasma indicates that smaller adjacent islands overlap with and destroy the large island flux surfaces. The estimated stochastic electron thermal conductivity (≈30 m 2 /s) is consistent with the conductivity inferred from the observed $T_e$ gradient and ohmic heating power. Island-shaped $T_e$ peaks can result from partially stochastic magnetic islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Bifurcations of axisymmetric plasma equilibrium in a tokamak

Bifurcation of solutions to the Grad–Shafranov-type equation for helically symmetric plasma near the threshold for tearing instability are analyzed. Quadratic and cubic nonlinearities were added to the linear dependence of the current density on the helical flux. Depending on the character of nonlinearity, two types of bifurcation can be observed, the “small” and the “large” ones. The small bifurcation is typical of cubic nonlinearity and reveals itself in the growth of the magnetic island from zero as the profile parameter increases above the instability threshold. The large bifurcation is typical of quadratic nonlinearity and causes jumplike formation of a large-scale magnetic island upon exceeding the instability threshold. As the profile parameter decreases below the instability threshold, the large-scale island continues to persist for some time (the hysteresis effect) and then suddenly disappears.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On the resonance amplification of magnetic perturbations near the threshold of tearing instability in a tokamak

Using a cylindrical model, a relatively simple description is presented of how a magnetic field perturbation stimulated by a low external helical current or a small helical distortion of the boundary and generating magnetic islands penetrates into a plasma column with a magnetic surface q=m/n to which tearing instability is attached. Linear analysis of the classical instability with an aperiodic growth of the perturbation in time shows that the perturbation amplitude in plasma increases in a resonant manner as the discharge parameters approach the threshold of tearing instability. In a stationary case, under the assumption on the helical character of equilibrium, which can be found from the two-dimensional nonlinear equation for the helical flux, there is no requirement for the small size of the island. Examples of calculations in which magnetic islands are large near the threshold of tearing instability are presented. The bifurcation of equilibrium near the threshold of tearing instability in plasma with a cylindrical boundary, i.e., the existence of helical equilibrium (along with cylindrical equilibrium) with large islands, is described. Moreover, helical equilibrium can also exist in the absence of instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

MAGNETIC NULL POINTS IN KINETIC SIMULATIONS OF SPACE PLASMAS

We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic particle-in-cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind, and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly (LMA) and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3–9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and LMA simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Full suppression of runaway electron generation by the mode penetration of resonant magnetic perturbations during disruptions on J-TEXT

The avoidance and suppression of runaway electron (RE) generation during disruptions is of great importance for the safe operation of tokamaks. Massive gas injection is enacted to suppress the generation of REs, but the poor gas mixing efficiency and extremely high density required to suppress RE generation make the full RE suppression unreliable. The magnetic perturbations provide an alternative RE suppression during disruptions. The use of mode penetration induced by resonant magnetic perturbations (RMPs) to suppress RE generation has been investigated on the J-TEXT tokamak. For a sufficiently long mode penetration duration, robust runaway suppression has been reached during the disruptions. The m/n = 2/1 mode RMP with high amplitude excites large magnetic islands inside the plasma and leads to the large-scale destruction of magnetic surfaces during disruptions, which results in RE loss and runaway-free disruptions. The essential island width required for runaway suppression is estimated to be larger than 0.16 as the minor radius. This value might be slightly underestimated because of the misalignment between the electron cyclotron emission diagnostic and the island O-point. NIMROD simulations are used to investigate the effect of magnetic islands on RE generation during disruption, showing that the large magnetic islands have the ability to enhance RE seed loss during disruptions. RMP can excite large magnetic islands in the target plasma without tearing mode and might be a way to prevent RE generation during disruptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Beaming of particles and synchrotron radiation in relativistic magnetic reconnection

Relativistic reconnection has been invoked as a mechanism for particle acceleration in numerous astrophysical systems. According to idealized analytical models, reconnection produces a bulk relativistic outflow emerging from the reconnection sites (X-points). The resulting radiation is therefore highly beamed. Using two-dimensional particle-in-cell simulations, we investigate particle and radiation beaming, finding a very different picture. Instead of having a relativistic average bulk motion with an isotropic electron velocity distribution in its rest frame, we find that the bulk motion of the particles in X-points is similar to their Lorentz factor γ, and the particles are beamed within ∼5/γ. On the way from the X-point to the magnetic islands, particles turn in the magnetic field, forming a fan confined to the current sheet. Once they reach the islands they isotropize after completing a full Larmor gyration and their radiation is no longer strongly beamed. The radiation pattern at a given frequency depends on where the corresponding emitting electrons radiate their energy. Lower-energy particles that cool slowly spend most of their time in the islands and their radiation is not highly beamed. Only particles that quickly cool at the edge of the X-points generate a highly beamed fan-like radiation pattern. The radiation emerging from these fast cooling particles is above the burn-off limit (∼100 MeV in the overall rest frame of the reconnecting plasma). This has significant implications for models of gamma-ray bursts and active galactic nuclei that invoke beaming in that frame at much lower energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Error field identification through torque balance on a saturated island in DIII-D

Measurement of the electromagnetic torque on a magnetic island could be an attractive method for error field identification in the early phase of ITER operation. Previous DIII-D experiments (Strait 2014 Nucl. Fusion 54 073004; Shiraki et al 2015 Plasma Phys. Control. Fusion 57 025016) have demonstrated the principle of this approach using a stationary or slowly rotating island, while recent developments in magnetic data analysis (Sweeney and Strait 2019 Phys. Plasmas 26 012509) allow the field of a rapidly rotating island to be readily distinguished from that of the wall currents induced by its rotation. In a recent experiment, a rotating n = 1 magnetic perturbation forced a saturated magnetic island to rotate, thus sampling all toroidal phases periodically in a single discharge. The phase and amplitude of the error field are inferred from analysis of the time-dependent torque balance on the island, including torques from the error field, the applied magnetic perturbation, and the wall currents induced by rotation of the applied perturbation and the island. Furthermore, results agree well with those from more conventional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analytical model of plasma response to external magnetic perturbation in absence of no-slip condition

Recent simulation and experimental results suggest that the magnetic island and flow on a resonant surface often do not satisfy the “no-slip” condition in the steady state. A new theoretical model on a nonlinear plasma response to external magnetic perturbation in the absence of a no-slip condition is proposed. The model is composed of the equations for the evolution of both the width and phase of the magnetic island due to forced reconnection driven by the external magnetic perturbation, and the force-balance equation for the plasma flow. When the island width is much less than the resistive layer width, the island growth is governed by the linear Hahm–Kulsrud–Taylor solution in the presence of time-dependent plasma flow. In the other regime when the island width is much larger than the resistive layer width, the evolution of both the island width and phase can be described using the Rutherford theory. The island solution is used to construct the quasilinear electromagnetic force, which, together with the viscous one, contributes to the nonlinear variation in plasma flow. The no-slip condition assumed in the conventional error field theory is not imposed here, where the island oscillation frequency depends on but does not necessarily equal to the plasma flow frequency at the rational surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On the magnetic reconnection of resistive tearing mode with the dynamic flow effects

Magnetic reconnection usually occurs in turbulent environments, which may not only provide anomalous resistivity to enhance reconnection rates but also significantly modify the reconnection process through direct nonlinear interaction with magnetic islands. This study presents numerical simulations investigating the effects of an imposed dynamic flow on magnetic reconnection, based on a two-dimensional reduced resistive MHD model. Results show that while the linear stability properties of the resistive tearing mode are moderately affected by the dynamic flow, nonlinear evolution is significantly modified by radial parity, amplitude, and frequency of the dynamic flow. After the slowly evolving nonlinear Rutherford stage, the reconnection process is found to progress in two phases by including the dynamic flow. A Sweet-Parker like current sheet is formed in the first phase. Afterwards, plasmoid instability is triggered in the second phase, where multiple plasmoids are continuously generated and ejected along the current sheet, leading to an impulsive bursty reconnection. The reconnection rate is considerably enhanced in the range of low resistivity as compared to without flow. We found that plasmoid instability onset and evolution are strongly influenced by the frequency and radial parity of the dynamic flows. The scaling of effective reconnection rates with the flow is found to be independent of resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Effects of plasmoid formation on sawtooth process in a tokamak

For realistic values of the Lundquist number in tokamak plasmas, 1/1 magnetic islands lead to the formation of a secondary thin current sheet, which breaks up into a chain of small magnetic islands, called plasmoids. The role of plasmoid dynamics during the sawtooth reconnection process in fusion plasmas remains an unresolved issue. In this study, systematic simulations are performed to investigate the resistive internal kink mode using the full resistive magnetohydrodynamics equations implemented in the NIMROD code in a simplified tokamak geometry. For Lundquist number S≥1.6×10 7 , the secondary current sheet is found to be unstable to plasmoids during the nonlinear resistive kink mode evolution with a critical aspect ratio of the current sheet of ~70. The merging of small plasmoids leads to the formation of a monster plasmoid that can significantly affect the primary island evolution. This may provide an explanation for the partial reconnection observed in sawtooth experiments.

Computer simulation↗
DOE OSTI2021

Stellarator optimization for good magnetic surfaces at the same time as quasisymmetry

A method is demonstrated to optimize a stellarator's geometry to eliminate magnetic islands and achieve other desired physics properties at the same time. For many physics quantities that have been used in stellarator optimization, including quasisymmetry, neoclassical transport, and magnetohydrodynamic stability, it is convenient to use a magnetic equilibrium representation that assures the existence of magnetic surfaces. However, this representation hides the possible presence of magnetic islands, which are typically undesirable. To include both surface-based objectives and island widths in a single optimization, two fixed-boundary equilibrium calculations are run at each iteration of the optimization: one that enforces the existence of magnetic surfaces (the Variational Moments Equilibrium Code) [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)] and one that does not (the Stepped Pressure Equilibrium Code) [Hudson et al., Phys. Plasmas 19, 112502 (2012)]. By penalizing the island residues in the objective function, the two magnetic field representations are brought into agreement during the optimization. Here, an example is presented in which, particularly on the surface where quasisymmetry was targeted, quasisymmetry is achieved more accurately than in previously published examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

An asymptotic-preserving Lagrangian algorithm for the time-dependent anisotropic heat transport equation

We propose a Lagrangian numerical algorithm for a time-dependent, anisotropic temperature transport equation in magnetized plasmas in the large guide field regime. The approach is based on an analytical integral formal solution of the parallel (i.e., along the magnetic field) transport equation with sources, and it is able to accommodate both local and non-local parallel heat flux closures. The numerical implementation is based on an operator-split formulation, with two straightforward steps: a perpendicular transport step (including sources), and a Lagrangian (field-line integral) parallel transport step. Algorithmically, the first step is amenable to the use of modern iterative methods, while the second step has a fixed cost per degree of freedom (and is therefore scalable). Accuracy-wise, the approach is free from the numerical pollution introduced by the discrete parallel transport term when the perpendicular to parallel transport coefficient ratio X ⊥ /X ∥ becomes arbitrarily small, and is shown to capture the correct limiting solution when ε = X⊥L 2 ∥ /X1L 2 ⊥ → 0 (with L∥∙ L⊥ , the parallel and perpendicular diffusion length scales, respectively). Therefore, the approach is asymptotic-preserving. We demonstrate the capabilities of the scheme with several numerical experiments with varying magnetic field complexity in two dimensions, including the case of transport across a magnetic island.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Electron acceleration in three-dimensional magnetic reconnection with a guide field

Kinetic simulations of 3D collisionless magnetic reconnection with a guide field show a dramatic enhancement of energetic electron production when compared with 2D systems. In the 2D systems, electrons are trapped in magnetic islands that limit their energy gain, whereas in the 3D systems the filamentation of the current layer leads to a stochastic magnetic field that enables the electrons to access volume-filling acceleration regions. The dominant accelerator of the most energetic electrons is a Fermi-like mechanism associated with reflection of charged particles from contracting field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

ELECTRON ACCELERATION BY CASCADING RECONNECTION IN THE SOLAR CORONA. I. MAGNETIC GRADIENT AND CURVATURE DRIFT EFFECTS

We investigate the electron acceleration by magnetic gradient and curvature drift effects in cascading magnetic reconnection of a coronal current sheet via a test particle method in the framework of the guiding center approximation. After several Alfvén transit times, most of the electrons injected at the current sheet are still trapped in the magnetic islands. A small fraction of the injected electrons precipitate into the chromosphere. The acceleration of trapped electrons is dominated by the magnetic curvature drifts, which change the parallel momentum of the electron, and appears to be more efficient than the acceleration of precipitating electrons, which is dominated by the perpendicular momentum change caused by the magnetic gradient drifts. With the resulting trapped energetic electron distribution, the corresponding hard X-ray (HXR) radiation spectra are calculated using an optically thin Bremsstrahlung model. Trapped electrons may explain flare loop top HXR emission as well as the observed bright spots along current sheets trailing coronal mass ejections. The asymmetry of precipitating electrons with respect to the polarity inversion line may contribute to the observed asymmetry of footpoint emission.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Energetic particle transport in 3D equilibrium. Final report

In this project, we have further developed the state-of-the-art fusion simulation code, gyrokinetic toroidal code (GTC), for studying the energetic particle (EP) transport in the 3D equilibrium. We have applied the new capabilities to simulate the EP transport by the Alfven eigenmodes in stellarators and by 3D magnetic islands in tokamaks. As part of V&V efforts for 3D capability, the 3D capabilities of GTC have been applied to study thermal transport in stellarators and tokamak with 3D magnetic fields by microturbulence, which could cause additional EP loss. Effects magnetic islands on neoclassical bootstrap current have also been studied. The computational capability developed and physics learned in this project are in support of fusion SciDAC Center for Integrated Simulation of Energetic Particles in Burning Plasmas (ISEP). During this project, the PI has supervised 6 postdoctoral researchers and 4 PhD students at the University of California, Irvine. A complete list of 49 publications and 23 invited talks partially supported by this grant since 7/2015 appears at the end of this report. Five key publications closely related to this project are highlighted in this report.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Reconnection dynamics with secondary tearing instability in compressible Hall plasmas

The dynamics of a secondary tearing instability is systematically investigated based on compressible Hall magnetohydrodynamic. It is found that in the early nonlinear phase of magnetic reconnection before onset of the secondary tearing instability, the geometry of the magnetic field in the reconnection region tends to form a Y-type structure in a weak Hall regime, instead of an X-type structure in a strong Hall regime. A new scaling law is found that the maximum reconnection rate in the early nonlinear stage is proportional to the square of the ion inertial length (γ∝d{sub i}{sup 2}) in the weak Hall regime. In the late nonlinear phase, the thin elongated current sheet associated with the Y-type geometry of the magnetic field breaks up to form a magnetic island due to a secondary tearing instability. After the onset of the secondary tearing mode, the reconnection rate is substantially boosted by the formation of the X-type geometries of magnetic field in the reconnection regions. With a strong Hall effect, the maximum reconnection rate linearly increases with the increase of the ion inertial length (γ∝d{sub i})

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗