Engineering PapersSearch

SEARCH · Engineering Papers

Results for “tearing modes”

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 73 records · Page 4

Neural network-based classification and regression of magnetohydrodynamic modes in tokamaks

We present a machine learning-based magnetohydrodynamic (MHD) classifier and regressor that utilizes real or complex-valued 3D magnetic sensor array data to determine neoclassical tearing mode (NTM) onset times in tokamaks with millisecond accuracy. The input dataset consists of poloidal profiles of complex Fourier amplitudes with an n = 1 toroidal mode number from 144 human-labeled ITER Baseline Scenario discharges in the DIII-D tokamak, spanning both tearing-dominated and sawtooth-dominated regimes. Since m, n = 2,1 NTMs frequently emerge alongside sawteeth at the same frequency in this scenario, the focus is on isolating the m = 1 and m = 2 components of the n = 1 MHD mode near the tearing onset. To improve model regularization and prediction stability, singular value decomposition was applied to balance the sawtooth and tearing datasets. The enriched datasets facilitated training neural networks that learn the key distinguishing features of sawtooth and tearing modes in the poloidal profiles of their magnetic amplitude and phase. When the modes occur independently, the networks achieve perfect classification due to the modes’ distinct characteristics and low measurement noise. In the more experimentally relevant case where both modes coexist, the networks maintain exceptional performance across key metrics. Tests on synthetic data with known ground truth demonstrate the superior accuracy of the neural network trained on complex-valued input compared to models using real amplitude, phase, or pseudo-complex data, achieving both a mean time delay and standard deviation below 1 ms. Notably, standard linear regression methods fitting the dominant singular modes to the data closely match the neural network’s performance. Applying these methods across a broad range of H-mode scenarios will enable future studies to systematically identify dominant NTM triggers as scenario-specific variables, paving the way for more effective tearing mode avoidance strategies in future fusion reactor designs.

machine learning

Magnetotail energy storage and the variability of the magnetotail current sheet

The antiparallel field configuration in the earth magnetotail is a prime location to search for the important astrophysical process known as magnetic reconnection. The magnetic evidence that reconnection occurs in the earth's magnetotail is that the energy of the tail invariably decreases at the time of global substorm onsets, and that an increased number of closed field lines are seen in the near-earth magnetotail after substorms. This increased northward flux is often convecting earthward from the expected site of the X line after the substorm onset, while southward flux is often convecting tailward. Prior to a substorm, energy accumulates in the tail as the field assumes a more taillike configuration with a smaller field component across the equatorial plane. The high-beta plasma sheet becomes thin and has been observed to be at least as thin as a few thousand km near substorm onset times. The small Bn component and thin plasma sheet are the conditions which favor the onset of the tearing-mode instability, which is thought to lead to neutral-line formation and the reconnection associated with substorm onset. Intervals of weak and highly variable fields in the plasma sheet provide evidence for the tearing-mode instability.

Fairfield, D. H.

Kinetic aspects of tail dynamics - Theory and simulation

Kinetic theories relevant to the geomagnetic tail are reviewed. The topics discussed include kinetic instabilities, simulations, and current-sheet particle acceleration. Tearing mode and reconnection theories are emphasized. Kinetic treatment is appropriate for these topics since the tail plasma is collisionless. Fluid calculations are appropriate when stochastic processes dominate and for studies where long wavelengths are important. However, fluid treatments of tearing modes and reconnection require a finite resistivity in the diffusion region. Thus, although 'anomalous resistivity' can be guessed or in some cases calculated, ideally the kinetic treatment is often to be preferred. Particle motion and acceleration in the current sheet can give rise to beam-like distributions in the plasma-sheet boundary layer. Studies of current-sheet particle motion have also been used as the basis for 'kinetic' tail equilibrium models. Furthermore, quite recently current-sheet particle motion is used directly in Coroniti's explosive tail reconnection model. The 'inertial conductivity' from the equilibrium models provides the 'dissipation' necessary for reconnection.

Speiser, T. W.

Transient brightenings of interconnecting loops. III Interpretation

Two alternative interpretations of the sudden X-ray brightenings observed in loops that interconnect active regions are presented. A fast tearing mode may be excited in those newly formed interconnecting loops within which sufficient magnetic free energy is stored to drive the mode. Alternatively, anomalous Joule heating driven by an inductive electric field parallel to the magnetic field, varying on a time scale of order of a minute, may cause the brightenings. It is suggested that it is plausible that the fast tearing mode may be the cause of brightenings in the young newly formed interconnecting loops, whereas the anomalous Joule heating might occur in old loop connections when an exteral disturbance propagates through them.

Spicer, D. S.

The root cause of disruptive NTMs and paths to stable operation in DIII-D ITER baseline scenario plasmas

Analyses of the DIII-D ITER Baseline Scenario database support that the disruptive m,n=2,1 magnetic islands are pressure gradient driven, non-linear instabilities seeded in a sequence of stochastic transient magnetic perturbations, and that the current profile relaxation does not affect the m,n=2,1 island onset rate. At low torque, these Neoclassical Tearing Modes are most commonly seeded by non-linear 3-wave coupling when the differential rotation between the q=1 & q=2 rational surfaces approaches zero. Lack of statistically significant difference between the current profiles of stable and unstable states, as well as lack of correlation between the tearing mode onset rate and the current profile relaxation both reject causality between the current profile evolution and the 2,1 magnetic island onsets in these plasmas. These support that preserving the differential rotation between the q=1 and q=2 rational surfaces is key to long pulse stable operation in the plasma scenario planned for ITER, while optimization of the current profile within the explored parameter space may lead to much weaker improvements than sustaining the differential rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Quantification of locked mode instability triggered by a change in confinement

Abstract This work presents the first analysis of the disruptive locked mode (LM) triggered by the dynamics of a confinement change. It shows that, under certain conditions, the LM threshold during the transient is significantly lower than expected from steady states. We investigate the sensitivity to a controlled n = 1 error field (EF) activated prior to the L-H transition in the COMPASS tokamak, at q 95 ∼ 3, β N ∼ 1, and using EF coils on the high-field side of the vessel. A threshold for EF penetration subsequent to the L-H transition is identified, which shows no significant trend with density or applied torque, and is an apparent consequence of the reduced intrinsic rotation of the 2/1 mode during this transient phase. This finding challenges the assumption made in theoretical and empirical works that natural mode rotation can be predicted by global plasma parameters and urges against using any parametric EF penetration scaling derived from steady-state experiments to define the EF correction strategy in the entire discharge. Furthermore, even at EFs below the identified penetration threshold, disruptive locking of sawtooth-seeded 2/1 tearing modes is observed after about 30% of L-H transitions without external torque.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Tearing at the dayside magnetopause

Initial observations of the dayside magnetopause thickness using magnetometers on board ISEE 1 (International Sun-Earth Explorer) and ISEE 2 indicate that this boundary may only be a few ion gyro-radii thick during periods of observed southward interplanetary magnetic field (IMF). Such a thin current layer can destabilize the collisionless tearing mode, a probable first step in the reconnection sequence. The linear growth rate is calculated assuming typical magnetopause parameters and it is found that it is a sensitive function of IMF orientation, magnetopause thickness, and electron number density. It is also shown that the finite size of the dayside magnetopause necessitates the existence of a two-dimensional tearing mode wave vector spectrum, a consequence of the inability of a thermal electron to maintain Landau resonance with the wave for a growth period. Implications regarding reconnection are discussed.

Quest, K. B.

Tearing Stable Stationary ITER Baseline Operation in DIII-D

We present a tearing stable, stationary and reproducible operating solution for DIII-D plasmas characterized by the ITER Baseline Scenario normalized parameter set and shape. In these plasmas low differential rotation (Δf1,2) between the core and edge is identified as the single root cause of the onset of the primary limiting instability, the 2,1 Neoclassical Tearing Mode (NTM), while a group of other examined parameters have a much weaker impact on the stability. Explanation is offered by drift kinetic simulations, which show a six-fold reduction in the NTM onset threshold due to diminished stabilizing ion polarization currents when the seed island drift frequency ceases in the local plasma frame. Based on the experimental observation and the explanation provided by the theory, a new method is developed to sustain Δf1,2 by modifying the edge neoclassical potential through neutral gas fueling. This method enables stationary ITER baseline operation free of disruptive 2,1 tearing modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Impulsive phase heating by uni-directional current systems in solar flares

The possibility of producing rapid heating of a substantial coronal volume by current dissipation, driven by a tearing instability associated with a single unidirectional current system, is investigated. For the fast tearing mode, it is found that, for initial growth rates gamma(f) equal or less than 0.3/s, the electron heating is offset by convective losses, resulting in a very slow temperature rise, and hard X-ray emitting temperatures are never realized. For the larger growth rates, heating from 10 to the 7th to 10 to the 8th K can be achieved in a few seconds; however, in this regime, the maximum volume heated is some three to five orders of magnitude less than the volume of heated material that is inferred from hard X-ray emission measures. The results suggest that to achieve the required heating for the fast tearing mode, a more complicated geometry involving multiple small-scale, oppositely directed current channels may be necessary.

La Rosa, T. N.

Analysis of the impact of parallel magnetic fluctuations on linear gyrokinetic stability in NSTX-U and verification of gyro-fluid models

In this work, we use the CGYRO gyrokinetic code to analyze two L- and one H-mode discharges from the National Spherical Torus Experiment (NSTX) and NSTX-Upgrade (NSTX-U) selected due to their different mix of ion-scale driftwaves, ion temperature gradient (ITG) mode and trapped electron mode (TEM), and electromagnetic instabilities, kinetic ballooning mode (KBM), and micro-tearing mode (MTM) in the plasma core. It is found that the effect of parallel magnetic fluctuations is strongly destabilizing to the unstable KBMs compared to calculations with only perpendicular magnetic fluctuations. Two discharges have a mix of ITG/TEM and MTMs that are predicted to be dominant instability across the plasma radius. The parallel magnetic fluctuations are found to have little effect on the MTM stability but are destabilizing to ITG/TEM modes. To test the validity of the gyro-fluid linear stability codes TGLF and GFS at low aspect ratio, a database of linear growth rates has been created using the CGYRO gyrokinetic code. The database is comprised of various parameter scans around a standardized set of NSTX-U core parameters. It contains a group of electrostatic cases and an electromagnetic group that includes the effects of perpendicular and parallel magnetic fluctuations. Comparing the results from the GFS and TGLF models, we find that GFS exhibits the best agreement with the database of CGYRO linear growth rates. Comparing the model results for the electromagnetic scans shows that GFS captures the effects of parallel magnetic fluctuations accurately, while the TGLF model does not, as it lacks sufficient perpendicular energy resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL

The effects of Ohmic heating and stable radiation on magnetic tearing

A study is made of the effect of a temperature-dependent Coulomb-like resistivity on the planar tearing mode. The local evolution of the temperature is described by an energy equation which includes Joule heating and optically thin radiation. The resulting system of coupled linear magnetohydrodynamic equations is solved numerically, and eigenfunctions and growth rates are obtained. In the absence of radiation, there are two distinct solutions above a critical value of the magnetic Reynolds number S, a tearing-like mode and a Joule-heating mode. Below this point, the growth rates coalesce into a conjugate-complex pair. When stable radiation (dR/dT greater than 0) is added, the heating mode disappears and a modified tearing excitation exists to much lower values of S before its growth is cut off by Ohmic heating. Examples are given for solar coronal parameters, and for those characteristic of fusion-research devices. The introduction of an effective value for the resistivity, in the presence of energy transport, allows a simple qualitative discussion of the different modes.

Tachi, T.

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

NSTX-U research advancing the physics of spherical tokamaks

Abstract The objectives of NSTX-U research are to reinforce the advantages of STs while addressing the challenges. To extend confinement physics of low- A , high beta plasmas to lower collisionality levels, understanding of the transport mechanisms that set confinement performance and pedestal profiles is being advanced through gyrokinetic simulations, reduced model development, and comparison to NSTX experiment, as well as improved simulation of RF heating. To develop stable non-inductive scenarios needed for steady-state operation, various performance-limiting modes of instability were studied, including MHD, tearing modes, and energetic particle instabilities. Predictive tools were developed, covering disruptions, runaway electrons, equilibrium reconstruction, and control tools. To develop power and particle handling techniques to optimize plasma exhaust in high performance scenarios, innovative lithium-based solutions are being developed to handle the very high heat flux levels that the increased heating power and compact geometry of NSTX-U will produce, and will be seen in future STs. Predictive capabilities accounting for plasma phenomena, like edge harmonic oscillations, ELMs, and blobs, are being tested and improved. In these ways, NSTX-U researchers are advancing the physics understanding of ST plasmas to maximize the benefit that will be gained from further NSTX-U experiments and to increase confidence in projections to future devices.

74 ATOMIC AND MOLECULAR PHYSICS

The Onset of Magnetic Reconnection in Tail-Like Equilibria

Magnetic reconnection is a fundamental mode of dynamics in the magnetotail, and is recognized as the basic mechanisms converting stored magnetic energy into kinetic energy of plasma particles. The effects of the reconnection process are well documented by spacecraft observations of plasmoids in the distant magnetotail, or bursty bulk flows, and magnetic field dipolarizations in the near Earth region. Theoretical and numerical analyses have, in recent years, shed new light on the way reconnection operates, and, in particular, which microscopic mechanism supports the dissipative electric field in the associated diffusion region. Despite this progress, however. the question of how magnetic reconnection initiates in a tail-like magnetic field with finite flux threading the current i.sheet remains unanswered. Instead, theoretical studies supported by numerical simulations support the point-of-view that such plasma and current sheets are stable with respect to collisionless tearing mode. In this paper, we will further investigate this conclusion, with emphasis on the question whether it remains valid in plasma sheets with embedded thin current sheets. For this purpose, we perform particle-in-cell simulations of the driven formation of thin current sheets, and their subsequent evolution either to equilibrium or to instability of a tearing-type mode. In the latter case we will pay particular attention to the nature of the electric field contribution which unmagnetizes the electrons.

Hesse, Michael

Fracture under combined modes in 4340 steel

An experimental investigation was conducted to study the interaction of combined modes of loading on crack instability in the presence of the opening and sliding modes of stress intensity factors, the opening and tearing modes of stress intensity factors, and all three modes of stress intensity factors. Through-cracked and surface-cracked flat and round specimens, and round notched bar specimens fabricated from high strength 4340 steel were used for the investigation. The results are evaluated to determine fracture criteria under the combined modes of stress intensity factors for the 4340 steel. These results are compared with the results of other investigators obtained for different materials.

Shah, R. C.

A study of resistive peeling–ballooning modes across spherical tokamaks

We investigate how non-ideal-magnetohydrodynamics (MHD) effects, in particular plasma resistivity, impact the peeling–ballooning stability thresholds in spherical tokamaks. This analysis follows the discovery of resistive kink-peeling modes in ELMing National Spherical Torus Experiment (NSTX) discharges. In the present study we extend this modeling to ELMing pulses in the Mega Ampere Spherical Tokamak (MAST) and MAST—Upgrade (MAST-U), where we find a clear resistive scaling for peeling–ballooning modes. While in NSTX ideal-MHD predicts stability for ELMing discharges, in MAST-U we find that the plasma is slightly unstable to peeling–ballooning modes, but is fully stabilized once diamagnetic effects are considered in terms of a growth rate normalization. A resistive power law scaling is calculated for these modes on MAST-U, which lies in between that of tearing modes and resistive interchange modes. A comparison between M3D-C1 and NIMROD shows reasonable agreement for this scaling. Resistivity destabilizes the modes and the peeling–ballooning unstable domain is considerably expanded in both, MAST and MAST-U. In addition to the MAST/-U pulses we also analyze resistive PB stability in a NSTX-similarity discharge on DIII-D. While having a different aspect ratio from NSTX, this discharge uses NSTX-like shaping parameters, toroidal field and plasma current. By considering these discharges alongside NSTX cases, we identify conditions influencing the onset of resistive peeling–ballooning modes. Furthermore, our findings indicate that magnetic shear in the pedestal region is closely linked to the emergence of resistive edge modes.

MAST-U

Radiative tearing - Magnetic reconnection on a fast thermal-instability time scale

Two energy modification mechanisms which are known to occur in sheared magnetic fields are the tearing and thermal instabilities. These processes can be studied separately with formalisms incorporating just the effective driving mechanism of interest (finite resistivity for the tearing mode and unstable radiation for the thermal mode). A model which includes both effects, and a temperature-dependent resistivity, indicates that modified forms of these two instabilities may coexist for identical physical conditions. When they are isolated computationally, one can show that their limiting growth rates are approximately those of the uncoupled instabilities. The spatial structure and energy content of these two new hybrid processes are then individually examined and are found to differ considerably from those obtained from separate treatments of the driving mechanisms. The faster radiative instability, which has a hydromagnetically scaled growth rate like the condensation mode of the thermal instability, is shown to involve a substantial amount of magnetic field reconnection. This can be partially explained by a large temperature drop (or resistivity rise) at the X-point. The island width of the Coulomb-coupled radiative mode is 30 percent of that produced by a comparable level of the slower tearing instability. In addition, the perturbed magnetic energy in the radiative instability is 5 times that of the perturbed thermal energy, indicating an appreciable modification of the initial magnetic structure.

Steinolfson, R. S.