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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2010

Current-interchange tearing modes: Conversion of interchange-type modes to tearing modes

It is shown that, in addition to usual neoclassical tearing modes, another type of nonclassical tearing mode exists in tokamaks: viz., current-interchange tearing modes (CITMs). CITMs are directly driven by unstable pressure-driven electromagnetic or electrostatic modes of the interchange type (e.g., interchange/ballooning modes, drift waves, etc.) due to the current gradient in tokamaks. Interchange-type modes exchange not only thermal and magnetic energies between flux tubes but also current. In a plasma with a current (or resistivity) gradient, such an interchange can create a current sheet at a mode resonance surface and result in the excitation of CITMs. Note that the interchange mode (i.e., Rayleigh-Taylor instability) is fundamental to tokamak physics. This new theory has an effect on both resistive magnetohydrodynamic stability and transport theories. Instabilities of the interchange type could be directly converted into CITMs, alternative to forming turbulent eddies through nonlinear coupling as in conventional transport theories. In particular, our CITM theory fills in the component in the transport theory of Rechester and Rosenbluth [Phys. Rev. Lett. 40, 38 (1978)] for the origin of magnetic island structure in axisymmetric tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Zonal flow dynamics in the double tearing mode with antisymmetric shear flows

The generation dynamics and the structural characteristics of zonal flows are investigated in the double tearing mode (DTM) with antisymmetric shear flows. Two kinds of zonal flow oscillations are revealed based on reduced resistive magnetohydrodynamics simulations, which depend on the shear flow amplitudes corresponding to different DTM eigen mode states, elaborated by Mao et al. [Phys. Plasmas 20, 022114 (2013)]. For the weak shear flows below an amplitude threshold, v{sub c}, at which two DTM eigen states with antisymmetric or symmetric magnetic island structure are degenerated, the zonal flows grow oscillatorily in the Rutherford regime during the nonlinear evolution of the DTMs. It is identified that the oscillation mechanism results from the nonlinear interaction between the distorted islands and the zonal flows through the modification of shear flows. However, for the medium shear flows above v{sub c} but below the critical threshold of the Kelvin-Helmholtz instability, an oscillatory growing zonal flow occurs in the linear phase of the DTM evolution. It is demonstrated that the zonal flow oscillation originates from the three-wave mode coupling or a modulation instability pumped by two DTM eigen modes with the same frequency but opposite propagating direction. With the shear flows increasing, the amplitude of zonal flow oscillation increases first and then decreases, whilst the oscillation frequency as twice of the Doppler frequency shift increases. Furthermore, impacts of the oscillatory zonal flows on the nonlinear evolution of DTM islands and the global reconnection are also discussed briefly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Transport and Dynamics in Toroidal Fusion Systems

The study entitled, "Transport and Dynamics in Toroidal Fusion Systems," (TDTFS) applied analytical theory and numerical computation to investigate topics of importance to confining plasma, the fourth state of matter, with magnetic fields. A central focus of the work is how non-thermal components of the ion particle distribution affect the "sawtooth" collective oscillation in the core of the tokamak magnetic configuration. Previous experimental and analytical research had shown and described how the oscillation frequency decreases and amplitude increases, leading to "monster" or "giant" sawteeth, when the non-thermal component is increased by injecting particle beams or by exciting ions with imposed electromagnetic waves. The TDTFS study applied numerical computation to self-consistently simulate the interaction between macroscopic collective plasma dynamics and the non-thermal particles. The modeling used the NIMROD code [Sovinec, Glasser, Gianakon, et al., J. Comput. Phys. 195, 355 (2004)] with the energetic component represented by simulation particles [Kim, Parker, Sovinec, and the NIMROD Team, Comput. Phys. Commun. 164, 448 (2004)]. The computations found decreasing growth rates for the instability that drives the oscillations, but they were ultimately limited from achieving experimentally relevant parameters due to computational practicalities. Nonetheless, this effort provided valuable lessons for integrated simulation of macroscopic plasma dynamics. It also motivated an investigation of the applicability of fluid-based modeling to the ion temperature gradient instability, leading to the journal publication [Schnack, Cheng, Barnes, and Parker, Phys. Plasmas 20, 062106 (2013)]. Apart from the tokamak-specific topics, the TDTFS study also addressed topics in the basic physics of magnetized plasma and in the dynamics of the reversed-field pinch (RFP) configuration. The basic physics work contributed to a study of two-fluid effects on interchange dynamics, where "two-fluid" refers to modeling independent dynamics of electron and ion species without full kinetic effects. In collaboration with scientist Ping Zhu, who received separate support, it was found that the rule-of-thumb criteria on stabilizing interchange has caveats that depend on the plasma density and temperature profiles. This work was published in [Zhu, Schnack, Ebrahimi, et al., Phys. Rev. Lett. 101, 085005 (2008)]. An investigation of general nonlinear relaxation with fluid models was partially supported by the TDTFS study and led to the publication [Khalzov, Ebrahimi, Schnack, and Mirnov, Phys. Plasmas 19, 012111 (2012)]. Work specific to the RFP included an investigation of interchange at large plasma pressure and support for applications [for example, Scheffel, Schnack, and Mirza, Nucl. Fusion 53, 113007 (2013)] of the DEBS code [Schnack, Barnes, Mikic, Harned, and Caramana, J. Comput. Phys. 70, 330 (1987)]. Finally, the principal investigator over most of the award period, Dalton Schnack, supervised a numerical study of modeling magnetic island suppression [Jenkins, Kruger, Hegna, Schnack, and Sovinec, Phys. Plasmas 17, 12502 (2010)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

The q-profile effect on high-order harmonic q = 1 tearing mode generation during sawtooth crashes

The effect of q-profiles on the excitation of high-order harmonic q=1 tearing modes during sawtooth crashes is investigated by a collisionless fluid model with the electron inertia term in Ohm's law. It is found that for a flat q-profile in the core region, the high-order harmonics, such as m/n=2/2 and/or m/n=3/3 modes, comparable to or stronger than the m/n=1/1 component, can be excited during tokamak sawteeth. The stronger the magnetic shear on the q=1 surface is, the more unstable the higher-m modes are. For smoothly monotonously increased q-profiles, a lower q value on the plasma edge tends to easily excite higher-m harmonics at the same level as the m = 1 mode simultaneously. The spatial characteristics of the eigenmodes in the cases with the typical q-profiles are also discussed. In addition, the basic feature of the magnetic island structures in the nonlinear evolution is numerically obtained, which is consistent qualitatively with the experimentally reconstructed phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

TWO-DIMENSIONAL MAGNETOHYDRODYNAMIC SIMULATIONS OF BARRED GALAXIES

Barred galaxies are known to possess magnetic fields that may affect the properties of bar substructures such as dust lanes and nuclear rings. We use two-dimensional high-resolution magnetohydrodynamic (MHD) simulations to investigate the effects of magnetic fields on the formation and evolution of such substructures, as well as on the mass inflow rates to the galaxy center. The gaseous medium is assumed to be infinitesimally thin, isothermal, non-self-gravitating, and threaded by initially uniform, azimuthal magnetic fields. We find that there exists an outermost x{sub 1}-orbit relative to which gaseous responses to an imposed stellar bar potential are completely different between inside and outside. Inside this orbit, gas is shocked into dust lanes and infalls to form a nuclear ring. Magnetic fields are compressed in dust lanes, reducing their peak density. Magnetic stress removes further angular momentum of the gas at the shocks, temporarily causing the dust lanes to bend into an 'L' shape and eventually leading to a smaller and more centrally distributed ring than in unmagnetized models. The mass inflow rates in magnetized models correspondingly become larger, by more than two orders of magnitude when the initial fields have an equipartition value with thermal energy, than in the unmagnetized counterparts. Outside the outermost x{sub 1}-orbit, on the other hand, an MHD dynamo due to the combined action of the bar potential and background shear operates near the corotation and bar-end regions, efficiently amplifying magnetic fields. The amplified fields shape into trailing magnetic arms with strong fields and low density. The base of the magnetic arms has a thin layer in which magnetic fields with opposite polarity reconnect via a tearing-mode instability. This produces numerous magnetic islands with large density that propagate along the arms to turn the outer disk into a highly chaotic state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Targeted Optimization of Quasi-Symmetric Stellarators

The proposed research focuses on targeted areas of plasma physics dedicated to improving the stellarator concept. Research was pursued in the technical areas of edge/divertor physics in 3D configurations, magnetic island physics in stellarators, the role of 3D shaping on microinstabilities and turbulent transport and energetic ion confinement in stellarators.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Closure of the single fluid magnetohydrodynamic equations in presence of electron cyclotron current drive

In the presence of electron cyclotron current drive (ECCD), the Ohm's law of single fluid magnetohydrodynamics is modified as E + v × B = η(J – J{sub EC}). This paper presents a new closure relation for the EC driven current density appearing in this modified Ohm's law. The new relation faithfully represents the nonlocal character of the EC driven current and its main origin in the Fisch-Boozer effect. The closure relation is validated on both an analytical solution of an approximated Fokker-Planck equation as well as on full bounce-averaged, quasi-linear Fokker-Planck code simulations of ECCD inside rotating magnetic islands. The new model contains the model put forward by Giruzzi et al. [Nucl. Fusion 39, 107 (1999)] in one of its limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Linear gyrokinetic theory and computation of the gyrocenter motion based on the exact canonical variables for axisymmetric tokamaks

Linear gyrokinetic theory based on the Lie-transform perturbation method is presented in terms of the exact canonical variables. In the linear drift approximation, it is shown that the gyrocenter equations of motion based on the canonical variables are equivalent to the usual guiding-center equations of motion. A numerical code is developed to advance the gyrocenter motion in terms of the exact canonical variables with arbitrary perturbations. It is found that a static magnetic island in a tokamak has little effect on the trapped particle orbits due to the conservation of the longitudinal invariant; and it induces the island structure of passing particle orbits due to the fact that the longitudinal invariant for the passing particles is broken by the asymmetric perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Vanishing Magnetic Shear And Electron Transport Barriers In The RFX-Mod Reversed Field Pinch

We define the safety factor q for the helical plasmas of the experiment RFX-mod by accounting for the actual three-dimensional nature of the magnetic flux surfaces. Such a profile is not monotonic but goes through a maximum located in the vicinity of the electron transport barriers measured by a high resolution Thomson scattering diagnostic. Helical states with a single axis obtained in visco resistive magnetohydrodynamic numerical simulations exhibit similar nonmonotonic q profiles provided that the final states are preceded by a magnetic island phase, like in the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

On the behavior of plasma with a negative current in the L-2M stellarator

The behavior of plasma with a negative current in the range of 5 kA ≤ vertical bar J vertical bar ≤ 16 kA in the L-2M stellarator is studied. Plasma is created and heated by means of electron cyclotron resonance heating (ECRH). Here, the definition “negative” stands for a toroidal current directed in such a way that it reduces the initial rotational transform µ in absolute value. It is shown that, only at relatively small values of the negative current ( vertical bar J vertical bar ∼ 6 kA), the increase in the shear leads to a reduction in electromagnetic turbulence. As vertical bar J vertical bar rises, the average level of magnetic fluctuations in the stationary phase of the discharge gradually increases. After ECRH is switched off (i.e., when plasma cools down), there is a time interval in which the amplitude of magnetic field fluctuations increases. No coherent three-dimensional modes (e.g., double tearing modes) that could accompany the penetration of the current into the plasma were observed. As for axisymmetric perturbations with m ≠ 0, analysis of the experimental data shows that they are absent at vertical bar J vertical bar ≤ 8 kA, appear in some discharges at vertical bar J vertical bar = 9 kA, and are observed in each discharge at vertical bar J vertical bar ≥ 10 kA. The frequency of these oscillations is on the order of 10 kHz, the poloidal mode number is m = 2, and the duration is 4–6 ms. The absence of oscillations with m = 1 and the small value of the external magnetic field for the m = 2 perturbation indicate that the radial size of the magnetic island is less than 1 cm. This, in turn, suggests that the radial profile of the plasma current changes with increasing vertical bar J vertical bar.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical experiments on the detailed energy conversion and spectrum studies in a corona current sheet

In this paper, we study the energy conversion and spectra in a corona current sheet (CS) by 2.5 dimensional MHD numerical simulations. Numerical results show that many Petschek-like fine structures with slow-mode shocks mediated by plasmoid instabilities develop during the magnetic reconnection process. The termination shocks can also be formed above the primary magnetic island and at the head of secondary islands. These shocks play important roles in generating thermal energy in a corona CS. For a numerical simulation with initial conditions close to the solar corona environment, the ratio of the generated thermal energy to the total dissipated magnetic energy is around 1/5 before secondary islands appear. After secondary islands appear, the generated thermal energy starts to increase sharply and this ratio can reach a value of about 3/5. In an environment with a relatively lower plasma density and plasma β, the plasma can be heated to a much higher temperature. After secondary islands appear, the one-dimensional energy spectra along the CS do not behave as a simple power law and the spectrum index increases with the wave number. The average spectrum index for the magnetic energy spectrum along the CS is about 1.8. The two-dimensional spectra intuitively show that part of the high energy is cascaded to large kx and ky space after secondary islands appear. The plasmoid distribution function calculated from numerical simulations behaves as a power law closer to f(ψ)∼ψ{sup −1} in the intermediate ψ regime. By using η{sub eff}=v{sub inflow}⋅L, the effective magnetic diffusivity is estimated to be about 10{sup 11} ∼ 10{sup 12} m{sup 2} s{sup −1}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

The tokamak density limit: A thermo-resistive disruption mechanism

The behavior of magnetic islands with 3D electron temperature and the corresponding 3D resistivity effects on growth are examined for islands with near-zero net heating in the island interior. We refer to the resulting class of non-linearities as thermo-resistive effects. In particular, the effects of varying impurity mix on the previously proposed local island onset threshold [Gates and Delgado-Aparicio, Phys. Rev. Lett. 108, 165004 (2012)] are examined and shown to be consistent with the well established experimental scalings for tokamaks at the density limit. A surprisingly simple semi-analytic theory is developed which imposes the effects of heating/cooling in the island interior as well as the effects of island geometry. For the class of current profiles considered, it is found that a new term that accounts for the thermal effects of island asymmetry is required in the modified Rutherford equation. The resultant model is shown to exhibit a robust onset of a rapidly growing tearing mode—consistent with the disruption mechanism observed at the density limit in tokamaks. A fully non-linear 3D cylindrical calculation is performed that simulates the effect of net island heating/cooling by raising/suppressing the temperature in the core of the island. In both the analytic theory and the numerical simulation, the sudden threshold for rapid growth is found to be due to an interaction between three distinct thermal non-linearities which affect the island resistivity, thereby modifying the growth dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerical study on the influence of electron cyclotron current drive on tearing mode

Controlling tearing modes by localized current drive is explored by using numerical simulation with a set of compressible magnetohydrodynamics equations. By examining the effects of different characteristics of driven current, such as density distribution, duration time, and deposition location, it is found that a driven current with larger magnitude and more focused deposition region shows a better suppression effect on the tearing modes. Meanwhile destabilizing effects are also observed when a driven current over a certain magnitude is applied continuously. In comparison with those on the X-point of the magnetic island, the results are better when the current deposition is targeted on the O-point. In addition, the timing control of the current deposition will be also addressed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Linear stability and nonlinear dynamics of the fishbone mode in spherical tokamaks

Extensive linear and nonlinear simulations have been carried out to investigate the energetic particle-driven fishbone instability in spherical tokamak plasmas with weakly reversed q profile and the q{sub min} slightly above unity. The global kinetic-MHD hybrid code M3D-K is used. Numerical results show that a fishbone instability is excited by energetic beam ions preferentially at higher q{sub min} values, consistent with the observed appearance of the fishbone before the “long-lived mode” in MAST and NSTX experiments. In contrast, at lower q{sub min} values, the fishbone tends to be stable. In this case, the beam ion effects are strongly stabilizing for the non-resonant kink mode. Nonlinear simulations show that the fishbone saturates with strong downward frequency chirping as well as radial flattening of the beam ion distribution. An (m, n) = (2, 1) magnetic island is found to be driven nonlinearly by the fishbone instability, which could provide a trigger for the (2, 1) neoclassical tearing mode sometimes observed after the fishbone instability in NSTX.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The magnetic furnace: intense core dynamos in B stars

The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through three-dimensional (3D) global simulations of convective core dynamos operating within a young 10 M{sub ⊙} B-type star, using the anelastic spherical harmonic code. These simulations capture the inner 65% of this star by radius, encompassing the convective nuclear-burning core (about 23% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05% to 16% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE ROLE OF KELVIN–HELMHOLTZ INSTABILITY FOR PRODUCING LOOP-TOP HARD X-RAY SOURCES IN SOLAR FLARES

We propose a model for the formation of loop-top hard X-ray (HXR) sources in solar flares through the inverse Compton mechanism, scattering the surrounding soft X-ray (SXR) photons to higher energy HXR photons. We simulate the consequences of a flare-driven energy deposit in the upper chromosphere in the impulsive phase of single loop flares. The consequent chromosphere evaporation flows from both footpoints reach speeds up to hundreds of kilometers per second, and we demonstrate how this triggers Kelvin–Helmholtz instability (KHI) in the loop top, under mildly asymmetric conditions, or more toward the loop flank for strongly asymmetric cases. The KHI vortices further fragment the magnetic topology into multiple magnetic islands and current sheets, and the hot plasma within leads to a bright loop-top SXR source region. We argue that the magnetohydrodynamic turbulence that appears at the loop apex could be an efficient accelerator of non-thermal particles, which the island structures can trap at the loop-top. These accelerated non-thermal particles can upscatter the surrounding thermal SXR photons emitted by the extremely hot evaporated plasma to HXR photons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Interacting tilt and kink instabilities in repelling current channels

We present a numerical study in resistive magnetohydrodynamics (MHD) where the initial equilibrium configuration contains adjacent, oppositely directed, parallel current channels. Since oppositely directed current channels repel, the equilibrium is liable to an ideal magnetohydrodynamic tilt instability. This tilt evolution, previously studied in planar settings, involves two magnetic islands or flux ropes, which on Alfvénic timescales undergo a combined rotation and separation. This in turn leads to the creation of (near) singular current layers, posing severe challenges to numerical approaches. Using our open-source grid-adaptive MPI-AMRVAC software, we revisit the planar evolution case in compressible MHD, as well as its extension to two-and-a-half-dimensional (2.5D) and full three-dimensional (3D) scenarios. As long as the third dimension can be ignored, pure tilt evolutions result that are hardly affected by out of plane magnetic field components. In all 2.5D runs, our simulations do show secondary tearing type disruptions throughout the near singular current sheets in the far nonlinear saturation regime. In full 3D runs, both current channels can be liable to additional ideal kink deformations. We discuss the effects of having both tilt and kink instabilities acting simultaneously in the violent, reconnection-dominated evolution. In 3D, both the tilt and the kink instabilities can be stabilized by tension forces. As a concrete space plasma application, we argue that interacting tilt-kink instabilities in repelling current channels provide a novel route to initiate solar coronal mass ejections, distinctly different from the currently favored pure kink or torus instability routes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2010

Initial performance of fully-coupled AMG and approximate block factorization preconditioners for solution of implicit FE resistive MHD.

This brief paper explores the development of scalable, nonlinear, fully-implicit solution methods for a stabilized unstructured finite element (FE) discretization of the 2D incompressible (reduced) resistive MHD system. The discussion considers the stabilized FE formulation in context of a fully-implicit time integration and direct-to-steady-state solution capability. The nonlinear solver strategy employs Newton-Krylov methods, which are preconditioned using fully-coupled algebraic multilevel (AMG) techniques and a new approximate block factorization (ABF) preconditioner. The intent of these preconditioners is to enable robust, scalable and efficient solution approaches for the large-scale sparse linear systems generated by the Newton linearization. We present results for the fully-coupled AMG preconditioner for two prototype problems, a low Lundquist number MHD Faraday conduction pump and moderately-high Lundquist number incompressible magnetic island coalescence problem. For the MHD pump results we explore the scaling of the fully-coupled AMG preconditioner for up to 4096 processors for problems with up to 64M unknowns on a CrayXT3/4. Using the island coalescence problem we explore the weak scaling of the AMG preconditioner and the influence of the Lundquist number on the iteration count. Finally we present some very recent results for the algorithmic scaling of the ABF preconditioner.

99 GENERAL AND MISCELLANEOUS↗