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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analytic Method to Estimate Particle Acceleration in Flux Ropes

The mechanism that accelerates particles to the energies required to produce the observed high-energy emission in solar flares is not well understood. Drake et al. (2006) proposed a kinetic mechanism for accelerating electrons in contracting magnetic islands formed by reconnection. In this model, particles that gyrate around magnetic field lines transit from island to island, increasing their energy by Fermi acceleration in those islands that are contracting. Based on these ideas, we present an analytic model to estimate the energy gain of particles orbiting around field lines inside a flux rope (2.5D magnetic island). We calculate the change in the velocity of the particles as the flux rope evolves in time. The method assumes a simple profile for the magnetic field of the evolving island; it can be applied to any case where flux ropes are formed. In our case, the flux-rope evolution is obtained from our recent high-resolution, compressible 2.5D MHD simulations of breakout eruptive flares. The simulations allow us to resolve in detail the generation and evolution of large-scale flux ropes as a result of sporadic and patchy reconnection in the flare current sheet. Our results show that the initial energy of particles can be increased by 2-5 times in a typical contracting island, before the island reconnects with the underlying arcade. Therefore, particles need to transit only from 3-7 islands to increase their energies by two orders of magnitude. These macroscopic regions, filled with a large number of particles, may explain the large observed rates of energetic electron production in flares. We conclude that this mechanism is a promising candidate for electron acceleration in flares, but further research is needed to extend our results to 3D flare conditions.

Energetic particle acceleration↗

CICART Center For Integrated Computation And Analysis Of Reconnection And Turbulence

CICART is a partnership between the University of New Hampshire (UNH) and Dartmouth College. CICART addresses two important science needs of the DoE: the basic understanding of magnetic reconnection and turbulence that strongly impacts the performance of fusion plasmas, and the development of new mathematical and computational tools that enable the modeling and control of these phenomena. The principal participants of CICART constitute an interdisciplinary group, drawn from the communities of applied mathematics, astrophysics, computational physics, fluid dynamics, and fusion physics. It is a main premise of CICART that fundamental aspects of magnetic reconnection and turbulence in fusion devices, smaller-scale laboratory experiments, and space and astrophysical plasmas can be viewed from a common perspective, and that progress in understanding in any of these interconnected fields is likely to lead to progress in others. The establishment of CICART has strongly impacted the education and research mission of a new Program in Integrated Applied Mathematics in the College of Engineering and Applied Sciences at UNH by enabling the recruitment of a tenure-track faculty member, supported equally by UNH and CICART, and the establishment of an IBM-UNH Computing Alliance. The proposed areas of research in magnetic reconnection and turbulence in astrophysical, space, and laboratory plasmas include the following topics: (A) Reconnection and secondary instabilities in large high-Lundquist-number plasmas, (B) Particle acceleration in the presence of multiple magnetic islands, (C) Gyrokinetic reconnection: comparison with fluid and particle-in-cell models, (D) Imbalanced turbulence, (E) Ion heating, and (F) Turbulence in laboratory (including fusion-relevant) experiments. These theoretical studies make active use of three high-performance computer simulation codes: (1) The Magnetic Reconnection Code, based on extended two-fluid (or Hall MHD) equations, in an Adaptive Mesh Refinement (AMR) framework, (2) the Particle Simulation Code, a fully electromagnetic 3D Particle-In-Cell (PIC) code that includes a collision operator, and (3) GS2, an Eulerian, electromagnetic, kinetic code that is widely used in the fusion program, and simulates the nonlinear gyrokinetic equations, together with a self-consistent set of Maxwell’s equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Wide Operational Windows of Edge-Localized Mode Suppression by Resonant Magnetic Perturbations in the DIII-D Tokamak

Edge-Localized-Mode (ELM) suppression by resonant magnetic perturbations (RMPs) generally occurs over very narrow ranges of the plasma current (or magnetic safety factor $q_{95}$) in the DIII-D tokamak. However, wide $q_{95}$ ranges of ELM suppression are needed for the safety and operational flexibility of ITER and future reactors. In DIII-D ITER Similar Shape (ISS) plasmas with $n$=3 RMPs, the range of $q_{95}$ for ELM suppression is found to increase with decreasing electron density. Nonlinear two-fluid MHD simulations reproduce the observed $q_{95}$ windows of ELM suppression and the dependence on plasma density, based on the conditions for resonant field penetration at the top of the pedestal. When the RMP amplitude is close to the threshold for resonant field penetration, only narrow isolated magnetic islands form near the top of the pedestal, leading to narrow $q_{95}$ windows of ELM suppression. However, as the threshold for field penetration decreases with decreasing density, resonant field penetration can take place over a wider range of $q_{95}$. For sufficiently low density (penetration threshold) multiple magnetic islands form near the top of the pedestal giving rise to continuous $q_{95}$ windows of ELM suppression. The model predicts that wide $q_{95}$ windows of ELM suppression can be achieved at substantially higher pedestal pressure in DIII-D by shifting to higher toroidal mode number ($n$=4) RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron and Proton Acceleration in Trans-relativistic Magnetic Reconnection: Dependence on Plasma Beta and Magnetization

Nonthermal electron acceleration via magnetic reconnection is thought to play an important role in powering the variable X-ray emission from radiatively inefficient accretion flows around black holes. The trans-relativistic regime of magnetic reconnection—where the magnetization σ, defined as the ratio of magnetic energy density to enthalpy density, is ~1—is frequently encountered in such flows. By means of a large suite of two-dimensional particle-in-cell simulations, we investigate electron and proton acceleration in the trans-relativistic regime. Additionally, we focus on the dependence of the electron energy spectrum on σ and the proton β (the ratio of proton thermal pressure to magnetic pressure). We find that the electron spectrum in the reconnection region is nonthermal and can be modeled as a power law. At low β, the slope, p, is independent of β and hardens with increasing σ as $p\simeq 1.8+0.7/\sqrt{\sigma }$. Electrons are primarily accelerated by the nonideal electric field at X-points, either in the initial current layer or in current sheets generated between merging magnetic islands. At higher values of β, the electron power law steepens, and the electron spectrum eventually approaches a Maxwellian distribution for all values of σ. At values of β near β max ≈ 1/4σ, when both electrons and protons are relativistically hot prior to reconnection, the spectra of both species display an additional component at high energies, containing a few percent of particles. Finally, these particles are accelerated via a Fermi-like process by bouncing between the reconnection outflow and a stationary magnetic island

43 PARTICLE ACCELERATORS↗

First observation of plasma healing via helical equilibrium in tokamak disruptions

A prompt transition of equilibrium states in a disruptive plasma is observed in the presence of three-dimensional (3D) resonant magnetic fields. A dominantly stochastic field line state heals into a 3D equilibrium with intact flux surfaces, resembling a single helicity magnetic island in the plasma core. As a consequence, the current quench process is bifurcated, i.e., from a completed quench accompanied by a vertically unstable disruption to a rapid termination of current decay followed by a full recovery to its predisruption value. Finally, this observation implies prompt stochastic fields healing through the rapid non-linear growth of magnetic island in plasma core may provide a plausible path for a soft-landing scenario of tokamak disruption in high plasma current regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energetic Particles of keV–MeV Energies Observed near Reconnecting Current Sheets at 1 au

We provide evidence for particle acceleration up to ∼5 MeV at reconnecting current sheets in the solar wind based on both case studies and a statistical analysis of the energetic ion and electron flux data from the five Advanced Composition Explorer Electron, Proton, and Alpha Monitor (EPAM) detectors. The case study of a typical reconnection exhaust event reveals (i) a small-scale peak of the energetic ion flux observed in the vicinity of the reconnection exhaust and (ii) a long-timescale atypical energetic particle event (AEPE) encompassing the reconnection exhaust. AEPEs associated with reconnecting strong current sheets last for many hours, even days, as confirmed by statistical studies. The case study shows that time-intensity profiles of the ion flux may vary significantly from one EPAM detector to another partially because of the local topology of magnetic fields, but mainly because of the impact of upstream magnetospheric events; therefore, the occurrence of particle acceleration can be hidden. The finding of significant particle energization within a time interval of ±30 hr around reconnection exhausts is supported by a superposed epoch analysis of 126 reconnection exhaust events. We suggest that energetic particles initially accelerated via prolonged magnetic reconnection are trapped and reaccelerated in small- or medium-scale magnetic islands surrounding the reconnecting current sheet, as predicted by the transport theory of Zank et al. Other mechanisms of initial particle acceleration can contribute also.

79 ASTRONOMY AND ASTROPHYSICS↗

Performance of Wendelstein 7-X stellarator plasmas during the first divertor operation phase

Wendelstein 7-X is the first comprehensively optimized stellarator aiming at good confinement with plasma parameters relevant to a future stellarator power plant. Plasma operation started in 2015 using a limiter configuration. After installing an uncooled magnetic island divertor, extending the energy limit from 4 to 80 MJ, operation continued in 2017. For this phase, the electron cyclotron resonance heating (ECRH) capability was extended to 7 MW, and hydrogen pellet injection was implemented. The enhancements resulted in the highest triple product (6.5 × 10 19 keV m ₋3 s) achieved in a stellarator until now. Plasma conditions [T e (0) ≈ T i (0) ≈ 3.8 keV, τ E > 200 ms] already were in the stellarator reactor-relevant ion-root plasma transport regime. Stable operation above the 2nd harmonic ECRH X-mode cutoff was demonstrated, which is instrumental for achieving high plasma densities in Wendelstein 7-X. Further important developments include the confirmation of low intrinsic error fields, the observation of current-drive induced instabilities, and first fast ion heating and confinement experiments. The efficacy of the magnetic island divertor was instrumental in achieving high performance in Wendelstein 7-X. Symmetrization of the heat loads between the ten divertor modules could be achieved by external resonant magnetic fields. Full divertor power detachment facilitated the extension of high power plasmas significantly beyond the energy limit of 80 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Endogenous and asymmetric magnetic reconnection with associated processes of relevance to fusion burning plasmas

We report an endogenous magnetic reconnection process is characterized by a driving factor that lays within the layer where a drastic change of magnetic field topology occurs. This kind of process is shown to take place in the presence of an electron temperature gradient in a well-confined plasma where, referring to quasi-collisionless regimes, the resulting electron temperature fluctuations can be anisotropic. Then a class of (radially) localized reconnecting modes is identified. These involve a transverse generated field $\tilde{B}_x$ of odd parity (as a function of the radial variable), and have finite (phase) velocities of propagation contrary to commonly considered reconnecting modes. The widths of the relevant reconnection layers remain significant even when large macroscopic distances are considered. Given that there are plasmas in the Universe with considerable electron thermal energy contents, these features can be relied upon in order to produce magnetic field generation, or conversion of magnetic energy into particle energy when the coupling of the localized odd modes to extended even modes can be significant. In any case, the resulting magnetic islands are not symmetric. With their excitation these modes can extract momentum from the main body of the plasma column which should recoil in the opposite direction. The excitation of antisymmetric endogenous modes is shown to be relevant to the electron temperature heating due to the reaction products in a fusion burning plasma as, in this case, the longitudinal thermal conductivity on selected rational magnetic surfaces can be decreased, relative to its collisional value, by the effects of reconnection. This kind of steepening is proposed to have a role in enhancing the growth rate of the instability involved in disruption events of the plasma column. The best agreement between theory and experiments concerning the onset of magnetic reconnection is (probably) represented by the theory of the internal kink mode. The observed accelerated reconnection rate following the onset is suggested as being explained by the formation of a relatively large magnetic island with a local steepening of the electron temperature gradient. A new kind of odd 'thermonuclear heating' mode associated with symmetric reconnection is identified.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Equilibrium β-limits in classical stellarators

Here, a numerical investigation is carried out to understand the equilibrium β-limit in a classical stellarator. The stepped-pressure equilibrium code is used in order to assess whether or not magnetic islands and stochastic field-lines can emerge at high β. Two modes of operation are considered: a zero-net-current stellarator and a fixed-iota stellarator. Despite the fact that relaxation is allowed, the former is shown to maintain good flux surfaces up to the equilibrium β-limit predicted by ideal-magnetohydrodynamics (MHD), above which a separatrix forms. The latter, which has no ideal equilibrium β-limit, is shown to develop regions of magnetic islands and chaos at sufficiently high β, thereby providing a ‘non-ideal β-limit’. Perhaps surprisingly, however, the value of β at which the Shafranov shift of the axis reaches a fraction of the minor radius follows in all cases the scaling laws predicted by ideal-MHD. We compare our results to the High-Beta-Stellarator theory of Freidberg and derive a new prediction for the non-ideal equilibrium β-limit above which chaos emerges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Role of fast ions in spontaneous neoclassical tearing mode instabilities in NSTX

Spontaneous neoclassical tearing modes (spontaneous NTMs) have been observed in the National Spherical Torus Experiment (NSTX). Recently, a TRANSP based model for the modified Rutherford equation analysis was developed to accurately predict the effect of fast ions on the growth of magnetic islands. The goal of this paper is to utilize the TRANSP NTM model to understand the role of fast ions in NSTX spontaneous NTMs. Our analysis shows that when the passing fast ion driven kinetic neoclassical polarization current term is larger than 1% of the bootstrap current term in the modified Rutherford equation, fast ions play a decisive role in the early growth phase of spontaneous NTM. This result is applicable to any tokamak plasmas with a fast ion drift orbit width, or the fast ion poloidal Larmor radius, comparable to the magnetic island width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗