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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
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 OSTI2011

SIESTA: A scalable iterative equilibrium solver for toroidal applications

A new solver for rapidly obtaining magnetohydrodynamic (MHD) equilibria in toroidal systems in the presence of islands and stochastic regions is described. It is based on the Kulsrud-Kruskal MHD energy minimization principle. To carry out this minimization, small displacements are made around a convenient set of curvilinear coordinates obtained from a nearby three-dimensional equilibrium that assumes nested surfaces. Because the changes of the magnetic fields and pressure are small, corresponding to small changes in the initial magnetic and kinetic energies, solutions for the linearized perturbations can be used to rapidly and iteratively find lower energy states with magnetic islands. A physics-based preconditioner is developed to accelerate the convergence of the iterative procedure to obtain an ideal MHD equilibrium with broken magnetic surfaces (islands).

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 OSTI2013

Ideal magnetohydrodynamic stability of configurations without nested flux surfaces

Existing numerical tools for calculating the MHD stability of magnetically confined plasmas generally assume the existence of nested flux surfaces. These tools are therefore not immediately applicable to configurations with magnetic islands or regions with an ergodic magnetic field. However, in practice, these islands or ergodic regions are often small, and their effect on MHD stability can then be evaluated using a perturbation theory developed in the present paper. This procedure allows the effect of the broken magnetic topology on the stability of each eigenmode to be calculated without requiring any knowledge about the perturbed eigenfunctions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulation of non-resonant internal kink mode with toroidal rotation in the National Spherical Torus Experiment

Plasmas in spherical and conventional tokamaks, with weakly reversed shear q profile and minimum q above but close to unity, are susceptible to an non-resonant (m,n) = (1,1) internal kink mode. This mode can saturate and persist and can induce a (2,1) seed island for Neoclassical Tearing Mode. [Breslau et al. Nucl. Fusion 51, 063027 (2011)]. The mode can also lead to large energetic particle transport and significant broadening of beam-driven current. Motivated by these important effects, we have carried out extensive nonlinear simulations of the mode with finite toroidal rotation using parameters and profiles of an NTSX plasma with a weakly reversed shear profile. The numerical results show that, at the experimental level, plasma rotation has little effect on either equilibrium or linear stability. However, rotation can significantly influence the nonlinear dynamics of the (1,1) mode and the induced (2,1) magnetic island. The simulation results show that a rotating helical equilibrium is formed and maintained in the nonlinear phase at finite plasma rotation. In contrast, for non-rotating cases, the nonlinear evolution exhibits dynamic oscillations between a quasi-2D state and a helical state. Furthermore, the effects of rotation are found to greatly suppress the (2,1) magnetic island even at a low level.

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 OSTI2013

Eigenmode characteristics of the double tearing mode in the presence of shear flows

The double tearing mode (DTM) is characterized by two eigen states with antisymmetric or symmetric magnetic island structure, referred to as the even or odd DTM. In this work, we systematically revisit the DTM instabilities in the presence of an antisymmetric shear flow with a focus on eigenmode characteristics as well as the stabilization or destabilization mechanism in a wide parameter region. Both initial value simulation and eigenvalue analysis are performed based on reduced resistive MHD model in slab geometry. A degenerated eigen state is found at a critical flow amplitude v{sub c}. The even (or odd) DTM is stabilized (or destabilized) by weak shear flow below v{sub c} through the distortion of magnetic islands mainly due to the global effect of shear flow rather than the local flow shear. The distortion can be quantified by the phase angles of the perturbed flux, showing a perfect correspondence to the growth rates. As the shear flow increases above v{sub c}, the degenerated eigen state bifurcates into two eigen modes with the same growth rate but opposite propagating direction, resulting in an oscillatory growth of fluctuation energy. It is identified that two eigen modes show the single tearing mode structure due to the Alfven resonance (AR) occurring on one current sheet. Most importantly, the AR can destabilize the DTMs through enhancing the plasma flow exerting on the remaining island. Meanwhile, the local flow shear plays a remarkable stabilizing role in this region. In addition, the eigenmode characteristic of the electromagnetic Kelvin-Helmholtz instability is also discussed.

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↗