Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear tearing mode theory”

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 19 records

Modeling q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

A toroidal asymptotic matching model of the response of a tokamak plasma to a static resonant magnetic perturbation (RMP) is used to simulate the n = 3 RMP-induced edge-localized-mode-suppression windows in q 95 that are evident when the plasma current is slowly ramped in DIII-D discharge #145380. All quantities employed in the simulation are derived from experimental measurements, apart from the neutral particle data. Three cases are considered. In the first case, the natural frequencies of tearing modes resonant in the plasma are determined by the ion flows at the corresponding resonant surfaces, which is the prediction of nonlinear tearing mode theory. In the second case, the natural frequencies are determined by the local E x B velocities at the resonant surfaces. In the third case, the natural frequencies are determined by the electron flows at the resonant surfaces, which is the prediction of linear tearing mode theory. The second case gives the best agreement between the simulations and the experimental observations. The first and third cases only lead to partial agreement between the simulations and the observations. In the first case, the lack of complete agreement may be a consequence of using an inaccurate assumption for the neutral particle distribution in the pedestal. In the third case, the lack of complete agreement is probably due to the fact that the response of a tokamak plasma to an RMP is not accurately described by linear tearing mode theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Global linear and nonlinear gyrokinetic simulations of tearing modes

Abstract To better understand multi-scale interactions between global tearing modes and microturbulence in the Madison Symmetric Torus (MST) reversed-field pinch (RFP), the global gyrokinetic code Gene is modified to describe global tearing mode instability via a shifted Maxwellian (SM) distribution consistent with experimental equilibria. The implementation of the SM is tested and benchmarked by comparisons with different codes and models. Good agreement is obtained in code-code and code-theory comparisons. Linear stability of tearing modes of a non-reversed MST discharge is studied. A collisionality scan is performed to the lowest order unstable modes ( n = 5, n = 6) and shown to behave consistently with theoretical scaling. The nonlinear evolution is simulated, and saturation is found to arise from mode coupling and transfer of energy from the most unstable tearing mode to small-scale stable modes mediated by the m = 2 tearing mode. The work described herein lays the foundation for nonlinear simulation and analysis of the interaction of tearing modes and gyroradius-scale instabilities in RFP plasmas.

Physics↗

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Does the resistive tearing instability nonlinearly evolve to a fast reconnection mode

A fundamental problem in applying linear tearing instability theory to the rapid processes (particle acceleration, heating) in flares was the characteristically slow rate of reconnection. This problem can be at least partially overcome if the tearin mode nonlinearly evolves to a regime in which the reconnection rate is substantially enhanced, such as that for the Petschek configuration. This possibility was often suggested, and some numerical simulations appear to provide support for such a view. Numerical simulation are used to study the nonlinear evolution of the tearing stability and show that a fast Petschek-like regime may not be achieved. This conclusion follows when there are sufficient grid points within the diffusion region to completely resolve the nonlinear dynamic interactions in the diffusion layer. When the numerical resolution is not adequate, the solution does appear to approach a Petschek configuration. The resolved solution contains reverse flow vortices and current sheets, terminated with a current reversal, similar to those obtained by Syrovatsky (JEPT, 33, 933, 1971).

Steinolfson, R. S.↗

Ballooning theory for micro-tearing mode in tokamak

This paper aims to investigate the impact of magnetic drift on the linear micro-tearing mode by using a kinetic approach to derive a reduced two-field eigen system in real space. Here, since the magnetic drift in real space has derivatives, it is more efficient to solve the mode equations in a Fourier-ballooning representation using the two-dimensional (2D) ballooning transform. The lowest-order eigen system in the Fourier-ballooning representation consists of two integral equations, which are numerically solved using the finite difference method for both eigenvalues and wave functions. The main results will be presented through graphical eigenvalue scans for each parameter. Furthermore, we present a graphical comparison between the predictions of the ballooning theory and GENE gyrokinetic code simulation in the pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Further modeling of q 95 windows for the suppression of edge localized modes by resonant magnetic perturbations in the DIII-D tokamak

An improved resonant plasma response model that more accurately captures the physics of the interaction between a tokamak plasma and a resonant magnetic perturbation (RMP) is developed. The model interpolates between the linear and the nonlinear response regimes and takes into account the fact that the slip-frequency is non-zero in the nonlinear regime. The improved model is incorporated into the extended perturbed equilibrium code (EPEC) toroidal asymptotic matching code. The modified EPEC code is used to investigate RMP-induced edge-localized-mode (ELM) suppression in DIII-D H-mode discharge #145380. Somewhat surprisingly, allowing for a finite slip-frequency (i.e., relaxing the so-called no-slip constraint) is found to only slightly facilitate the locking of driven magnetic island chains to the RMP, and, hence, to only slightly facilitate RMP-induced ELM suppression. This is true despite the fact that the nature of non-locked island solutions is radically different when the no-slip constraint is imposed compared to when it is relaxed (in the first case, the widths of the island chains driven at the rational surfaces pulsate, and in the second case, they remain steady). The previously obtained conclusion that the response of a typical H-mode tokamak plasma to an RMP cannot be accurately modeled by linear theory is confirmed. Here, the previously obtained conclusion that the best agreement between theory and observations is achieved by assuming that the natural frequencies of tearing modes, in the absence of the RMP, are determined by the local equilibrium E × B velocity is also confirmed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The electrostatic effect for the collisionless tearing mode

Electron dynamics has not been self-consistently considered in collisionless tearing mode theories to date because of the mathematical complexity of the Vlasov-Maxwell equations. It has been found, using computer simulations, that electrostatic fields play an important role in the tearing mode. Vlasov theory, including the electrostatic field, is investigated for topologies with both antiparallel and nonantiparallel magnetic field lines. The electrostatic field influences the resonant current in the neutral sheet which is a non-MHD effect, and modifies the linear growth rate. At the magnetopause, where the field lines are not antiparallel, the electrostatic effect acts to raise the linear growth rate of the tearing mode. On the other hand, in the magnetotail, where magnetic field lines are antiparallel, the electrostatic effect reduces the tearing mode growth rate.

Hoshino, Masahiro↗

A new approach to the linear theory of single-species tearing in two-dimensional quasi-neutral sheets

We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Brittnacher, M.↗

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisionless tearing instability in magnetotail plasmas

The problem of the linear stability of collisionless tearing modes in the earth's magnetotail is revisited. It is found that the collisionless tearing mode is linearly unstable with wavelengths of the order of 10 R(E). It is shown that an important feature neglected in earlier theories is a nonzero equilibrium B(y)-field. The physics of the instability is elucidated in the context of a simple slab model and a sheared parabolic model which is representative of the magnetotail in which all three components of the magnetic field are nonzero. The dispersion equation for the instability is obtained by a boundary-layer analysis. The implications of the theory for recent observations on current disruption and diversion during substorms is discussed.

Wang, Xiaogang↗

New linear stability parameter to describe low-β electromagnetic microinstabilities driven by passing electrons in axisymmetric toroidal geometry

In magnetic confinement fusion devices, the ratio of the plasma pressure to the magnetic field energy, β, can become sufficiently large that electromagnetic microinstabilities become unstable, driving turbulence that distorts or reconnects the equilibrium magnetic field. In this paper, a theory is proposed for electromagnetic, electron-driven linear instabilities that have current layers localised to mode-rational surfaces and binormal wavelengths comparable to the ion gyroradius. The model retains axisymmetric toroidal geometry with arbitrary shaping, and consists of orbit-averaged equations for the mode-rational surface layer, with a ballooning space kinetic matching condition for passing electrons. The matching condition connects the current layer to the large scale electromagnetic fluctuations, and is derived in the limit that β is comparable to the square root of the electron-to-ion-mass ratio. Electromagnetic fluctuations only enter through the matching condition, allowing for the identification of an effective β that includes the effects of equilibrium flux surface shaping. The scaling predictions made by the asymptotic theory are tested with comparisons to results from linear simulations of micro-tearing and electrostatic microinstabilities in MAST discharge #6252, showing excellent agreement. In particular, it is demonstrated that the effective β can explain the dependence of the local micro-tearing mode (MTM) growth rate on the ballooning parameter θ 0 – possibly providing a route to optimise local flux surfaces for reduced MTM-driven transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Linear tearing modes of a forced current-sheet equilibrium

Linear tearing modes are studied in a nearly singular forced current-sheet equilibrium, such as could result from global magnetic forces in the solar corona. Growth rates for the tearing modes, determined by solving the linearized reduced MHD (Strauss) equations numerically, were found to scale as (gamma)tau(d) = about S exp 4/5 k(y) exp 4/5, where S is the Lundquist number, k(y) is the wavenumber, and tau(d) is the classical resistive diffusion time. This scaling is in agreement with predictions from analytical theory. Because of the faster S scaling of these modes compared to the tearing modes of a diffuse current-sheet equilibrium, the modes have much higher growth rates (by a factor of about 10,000) for coronal values of S (about 10 to the 12th). For coronal parameters, the growth times of these new modes are estimated to be on the order of several hours to days, as compared to growth times of months to years for the tearing modes in a diffuse current-shear equilibrium. The growth times are comparable to reconnection times scales required in models of coronal heating by magnetic field dissipation.

Liewer, Paulett C.↗

Predicted thresholds for RMP ELM suppression access in double-null configurations

Modeling of DIII-D plasmas spanning shapes from single to double-null (DN) reveals new insights into the nature of resonant magnetic perturbation (RMP) conditions necessary for edge-localized mode (ELM) suppression. The suppression of ELMs with RMPs has proven difficult in DN configurations, where no device has thus far reported any hints of suppression. Modeling using the GPEC code finds a reduced high-field side response closer to DN shaping. The resulting synthetic diagnostic measurements are consistent with what has been observed in experiments on DIII-D, validating the plasma response model in this regime. While common metrics for suppression do not illustrate a clear distinction between single-null (SN) and DN cases, the pedestal top resonant field does show a ∼20% decrease at DN shaping in modeling. Field penetration is assessed using linear tearing theory, which demonstrates a lack of sufficient pedestal top resonant flux in the DN shape, requiring at least 1.5 × greater RMP coil currents than what was used in experiment. Analysis from drift kinetic simulations further indicate up to 2 × larger critical island widths are required at the pedestal top for tearing mode growth compared to SN cases. Effective island widths inferred from 3D ideal MHD are also analyzed, where maximum widths in lab coordinates indicate a threshold of ∼18–24× the ion gyroradius for sufficient profile flattening for ELM suppression. These results suggest that ELM suppression may be possible in DN with sufficiently large RMP coil amplitude. Future prospects of achieving RMP ELM suppression in the DN configuration may involve going to lower triangularity, which is also highlighted in this work.

RMP ELM suppression↗

Collisionless drift-tearing modes in the magnetopause

The linear stability properties of collisionless drift-tearing modes are analyzed in a modified Harris equilibrium model of the magnetopause. Particular attention is paid to the relevance of the parametric behavior of growth rates to the 'magnetic percolation' theory of flux transfer event formation (Galeev et al., 1986). Numerical methods are used to solve the drift-tearing eigenmode equations and the results are compared with those previously obtained by analytical methods. The analytical results are found to correctly model important parametric dependencies but to typically overestimate the rate of growth. The eigenmode equations are numerically difficult, and an integration scheme utilizing Ricatti transforms is developed to affect their solution.

Gladd, N. T.↗

Linear theory of tearing in a high-beta plasma

The linear dispersion relation for a collisionless plasma in a sheared one-dimensional current sheet is calculated with reference to conditions in the daytime magnetopause. Calculations are extended to include plasmas with beta approximately equal to 1. It is found that the tearing mode eigenstructure and temporal growth rate are a sensitive function of the ratios l sub s/l sub G, l sub s/l sub T, and l sub s/l sub n, where l sub s is the shearing length of the magnetic field, l sub G is the gradient scale length, and l sub T is the temperature gradient scale length. In particular, if beta is approximately equal to 1, and l sub s is less than l sub G, l sub n, and l sub T, then the thickness of the layer over which particles are resonantly accelerated by the induction magnetic field is approximately rho, a thermal gyroradius. If the above conditions are not satisfied, plasma gradients may electrostatically stabilize the mode.

Quest, K. B.↗

Lower-hybrid-drift wave turbulence in the distant magnetotail

Recent satellite observations of electrostatic and magnetic noise in the distant magnetotail (Gurnett et al., 1976) can be explained by the excitation of the lower-hybrid-drift instability. In particular, it is shown that (1) existence conditions for the lower-hybrid-drift instability can be met, (2) the observed frequency spectra and polarization are in good agreement with the predictions of linear theory, and (3) the observed amplitudes of fluctuations are consistent with the nonlinear theory of this mode. Moreover, the observation of this instability suggests that the anomalous transport properties associated with these waves, which are important in many laboratory devices, may play a crucial role in the macroscopic evolution of magnetotail processes such as field line merging, tearing instabilities, or 'fireballs'.

Huba, J. D.↗

Gyrokinetic simulation of pedestal degradation correlated with enhanced magnetic turbulence in a DIII-D ELMy H-mode discharge

Gyrokinetic simulation of a dedicated pedestal density ramping-up discharge on DIII-D can reproduce the enhancement of magnetic turbulence in the pedestal, which is identified to be caused by micro-tearing modes (MTMs). Here, an increase of MTM amplitude results in higher electron thermal diffusivity, consistent with experimentally observed lower electron temperature gradient and degraded pedestal height. Gyrokinetic simulation identifies the major cause of MTM enhancement to be the increase of collisionality, which has a significant impact on the MTM intensity and is beyond the description of any (quasi-)linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetotail dynamics under isobaric constraints

Using linear theory and nonlinear MHD simulations, we investigate the resistive and ideal MHD stability of two-dimensional plasma configurations under the isobaric constraint dP/dt = 0, which in ideal MHD is equivalent to conserving the pressure function P = P(A), where A denotes the magnetic flux. This constraint is satisfied for incompressible modes, such as Alfven waves, and for systems undergoing energy losses. The linear stability analysis leads to a Schroedinger equation, which can be investigated by standard quantum mechanics procedures. We present an application to a typical stretched magnetotail configuration. For a one-dimensional sheet equilibrium characteristic properties of tearing instability are rediscovered. However, the maximum growth rate scales with the 1/7 power of the resistivity, which implies much faster growth than for the standard tearing mode (assuming that the resistivity is small). The same basic eigen-mode is found also for weakly two-dimensional equilibria, even in the ideal MHD limit. In this case the growth rate scales with the 1/4 power of the normal magnetic field. The results of the linear stability analysis are confirmed qualitatively by nonlinear dynamic MHD simulations. These results suggest the interesting possibility that substorm onset, or the thinning in the late growth phase, is caused by the release of a thermodynamic constraint without the (immediate) necessity of releasing the ideal MHD constraint. In the nonlinear regime the resistive and ideal developments differ in that the ideal mode does not lead to neutral line formation without the further release of the ideal MHD constraint; instead a thin current sheet forms. The isobaric constraint is critically discussed. Under perhaps more realistic adiabatic conditions the ideal mode appears to be stable but could be driven by external perturbations and thus generate the thin current sheet in the late growth phase, before a nonideal instability sets in.

Birn, Joachim↗