Engineering Papers⌕ Search

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

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↗
NASA NTRS1974

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.↗
DOE OSTI2025

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↗
NASA NTRS1984

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.↗
NASA NTRS1991

Collisionless reconnection in two-dimensional magnetotail equilibria

A two-dimensional particle simulation model based on the Darwin approximation to Maxwell's equations for studying collisionless reconnection in the magnetotail has been developed. Simulations of the pure ion tearing mode in a thin current sheet with normal B(z) field component demonstrate that in this limit this mode grows more slowly than expected based on previous analytic estimates. The saturation level of the tearing instability greatly surpasses estimates based on a simple trapping argument. The effect of the normal field component on the evolution of the tearing instability is considered. It is found that a normal field of even a few percent on axis strongly inhibits the growth of the instability.

Pritchett, P. L.↗
DOE OSTI2024

Macroscopic trends of linear tearing stability in cylindrical current profiles

Abstract The likelihood of realising tokamak power-plants will be greatly improved by the discovery of high-gain equilibria that resist the formation of small islands and hence avoid the disruptive neoclassical tearing mode. We propose a series of studies to understand how simple tokamak design can leverage aspects of tearing onset physics to maximise passive resistance to island formation. Here we investigate the variation that current profiles can bring about in preventing tearing onset through the cylindrical linear tearing stability parameter Δ ′ . A database of 159148 realistic pilot-plant current profiles was generated with Monte Carlo sampling, and the distribution of Δ ′ values was linked with interpretable profile characteristics. In agreement with prior theoretical and experimental studies, Δ ′ was found to be strongly correlated with the existence and steepness of a local toroidal current well or hill, with the former destabilising and the latter stabilising. In the absence of these two cases, the remaining Δ ′ values were linearly bounded by the toroidal current gradient at the rational surface.

Physics↗
NASA NTRS1984

Association between energetic particle bursts and Birkeland currents in the geomagnetic tail

IMP 8 observations in the geomagnetic tail of energetic particles at hundreds of keV and of the magnetic field are examined for a substorm period from 1700 to 1800 UT on Mar. 3, 1976. For this interval, it is found that bursts of energetic protons detected at a downstream distance of about 37 earth radii in the midnight tail region are accompanied by east-west magnetic field perturbations indicative of Birkeland current sheets. The magnitude of the magnetic field perturbations range from about 1 to 12 nT, corresponding to current densities of 1 to 9 mA/m when integrated over the current sheet thickness. These values are comparable to previous observations of Birkeland currents in the near-earth tail (X greater than -20 earth radii) and with the auroral-arc-associated Birkeland current densities at low altitudes scaled out to the tail region. The favorable comparison suggests that the closure of small-scale Birkeland current systems, at least in these instances, lies tailward of about 40 earth radii. The observed Birkeland current system and the magnitude of current densities are found to be consistent with ion tearing instability in the tail. The inferred wavelength of the tearing mode along the tail axis is about 2 to 17 earth radii. In the context of the tearing model, the presence of Birkeland currents and the occasional isotropic distribution of energetic electrons (of more than 0.22 MeV) suggest that the energetic particle events in this study occur on closed magnetic field lines, even though strong tailward streaming dominates the particle anisotropy in these events.

Lui, A. T. Y.↗
DOE OSTI2025

The dual-electron cyclotron emission based measurement of 3D structures on DIII-D tokamak

The measurements of the locked mode and other 3D magnetohydrodynamic modes are essential for real-time control of plasma stability in future fusion reactors. Current magnetic measurements present non-negligible challenges for the application on fusion pilot plants (FPPs) due to neutron radiation and pose challenging design issues. In this paper, we discuss an alternative solution with microwave-based measurements. The dual-electron cyclotron emission (ECE) method has been developed and applied on DIII-D. The dual-ECE measurement method demonstrates clear identification of magnetic islands compared to single ECE, which is highly beneficial for locked mode feedback control on the FPPs. Toroidally separated ECE measurements provide accurate three-dimensional structure of magnetic islands and tearing modes. This technique also offers deeper physical insights into the penetration of resonant magnetic perturbations. Numerical calculations indicate that small locked modes (magnetic island width > 2 cm) can be detected using the dual-ECE measurement method in ITER plasma. This will enhance measurement capabilities for the safe operation of FPP plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Can prominences form in current sheets

Two-dimensional numerical simulations of the formation of cold condensations in a vertical current sheet have been performed using the radiative, resistive MHD equations with line-tied boundary conditions at one end of the sheet. Prominence-like condensations are observed to appear above and below an X-line produced by the onset of the tearing-mode instability. Cooling in the sheet is initiated by Ohmic decay, with the densest condensations occurring in the region downstream of a fast-mode shock. This shock, which is due to the line-tied boundary conditions, terminates one of the two supermagnetosonic reconnection jets that develop when the tearing is fully developed. The condensation properties of shock waves, which may trigger or considerably enhance the conditions for thermal condensation are emphasized.

Malherbe, J. M.↗
DOE OSTI2025

Integrated modeling of alpha particle losses in JET DT plasmas

Alpha particle confinement is crucial for sustaining burning plasmas and designing future reactor concepts. Along with classical/prompt losses, various magnetohydrodynamic instabilities can lead to wave-particle interactions which can transport alpha particles outward from the plasma. This can result in a reduction in plasma heating/performance, and, at worst, damage in-vessel components. Joint European Torus’s recent deuterium–tritium campaigns in 2021–2023 have produced numerous alpha particle loss measurements with its scintillator probe and Faraday cup array fast ion loss detectors as discussed in Bonofiglo et al (2024 Nucl. Fusion 64 096038). This paper will report on integrated energetic particle transport modeling in support of those measurements. The modeling is accomplished with the TRANSP and ORBIT-kick codes with the use of recently developed reduced models which calculate mode structure, amplitude, and the evolving dynamics. When possible, constraints and comparisons to experiment are conducted. Case studies are performed on a variety of magnetohydrodynamic activity, including: fishbones, tearing modes (TMs), and sawtooth crashes. Additionally, a special case of an alpha-driven toroidal Alfvén eigenmode is briefly discussed, where modeling showed marginally weak alpha losses and was unable to support experimental observations. Coupled effects between a TM and toroidal field ripple are presented and were unable to replicate the observations in lost particle pitch but did duplicate the localized flattening of the measured neutron profile. Additional modeling results compare the magnitude of losses and energy/velocity-space sensitivities against experimental observations/measurements for each scenario. This work corroborates numeric alpha transport modeling while also identifying model deficiencies. While this report details alpha transport, it also presents open issues for discussion in assessing the validity of our numerical models towards burning plasmas.

Bonofiglo, Phillip J. [Princeton Plasma Physics La↗
DOE OSTI2024

Inter-ELM pedestal turbulence dynamics dependence on q 95 and temperature gradient

A series of dedicated experiments from DIII-D tokamak provide spatially and temporally resolved measurements of electron density and temperature, and multiscale and multichannel fluctuations over a wide range of conditions. Measurements of long wavelength density fluctuations in the type-I ELMing H-mode pedestal routinely revealed a coexistence of multiple instabilities that exhibit dramatic different dynamic behaviors as q 95 and temperature gradient are varied, apparently responsible for limiting pedestal temperature profiles. Two distinct frequency bands of density fluctuations are modulated with ELM cycle with frequency above 200kHz propagating in the electron diamagnetic direction in the lab frame (electron mode) and below 200kHz propagating in the ion diamagnetic direction (ion mode). The electron mode amplitude peaks near the electron temperature gradient region and increases with q 95 which seems to be correlated with the increased χ e at higher q 95 , similar to the characteristics expected for Micro-tearing Mode (MTM). At higher q 95 , during the inter-ELM period, the ion mode decays at later phase of the ELM cycle. Consistently, the poloidal correlation length of the ion mode is also found to reduce which suggests the possible E×B flow shear suppression of ion mode at later phase of the ELM cycle as the Er well recovers. In contrast, the electron mode grows during the ELM cycle and reaches saturation at around 50-60% of ELM period. Linear gyrokinetic simulations find the MTMs as the most unstable mode in the pedestal electron temperature gradient region. The higher q 95 and lower magnetic shear destabilize MTMs. These observations provide key insights of the underlying physics of multifield properties and rich dataset of experimental 'fingerprints' that enable new tests of the theoretical pedestal models and lead to developing a predictive model for pedestal formation on ITER and future burning plasma experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1989

Turbulent excitation of spontaneous reconnection

The long-term nonlinear evolution of a tearing-mode-unstable sheared-field plasma in a turbulent environment is explored. Two different physical configurations are modeled, and a different computational system is used for each. Results of both sets of calculations show that magnetic tearing arises spontaneously provided that the initial turbulence energy level is below the natural saturation level of the tearing instability. The relationship between these results and those of previous calculations are briefly discussed, concluding that there are no significant unexplainable disagreements.

Deeds, D.↗
NASA NTRS1994

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↗
NASA NTRS1992

Secondary instability in three-dimensional magnetic reconnection

We consider the transition to turbulence in three-dimensional reconnection of a magnetic neutral sheet. We find that the transition can occur via a three-step process. First, the sheet undergoes the usual tearing instability. Second, the tearing mode saturates to form a two-dimensional quasi-steady state. Third, this secondary equilibrium is itself unstable when it is perturbed by three-dimensional disturbances. Most of this paper is devoted to the analysis and simulation of the three-dimensional linear stability properties of the two-dimensional saturated tearing layer. The numerical simulations are performed with a semi-implicit, pseudospectral-Fourier collocation algorithm. We identify a three-dimensional secondary linear stability which grows on the ideal timescale. An examination of the modal energetics reveals that the largest energy transfer is from the mean field to the three-dimensional field, with the two-dimensional field acting as a catalyst.

Dahlburg, R. B.↗
NASA NTRS1985

The magnetohydrodynamics of current sheets

Examples of current sheets are summarized and their formation is described. A universal phenomenon in cosmic plasmas is the creation of sheets off intense current near X-type neutral points (where the magnetic field vanishes). These sheets are important as sites where the magnetic-field energy is converted efficiently into heat and bulk kinetic energy and where particles can be accelerated to high energies. Examples include disruptions in laboratory tokamaks, substorms in the earth's magnetosphere, and flares on the sun. The basic behavior of a one-dimensional sheet is presented, together with an account of the linear tearing-mode instability that can cause the field lines in such a sheet to reconnect. Such reconnection may develop in different ways: it may arise from a spontaneous instability or it may be driven, either from outside by motions or locally by a resistivity enhancement. Various processes are described that may occur during the nonlinear development of tearing, along with the many numerical and laboratory experiments that are aiding our understanding of this intriguing cosmical process.

Priest, E. R.↗
DOE OSTI2024

Verification of electromagnetic simulation capabilities in global gyrokinetic particle-in-cell code GTS

Recently, the numerical scheme presented by Mishchenko et al. enabled explicit gyrokinetic simulations of low-frequency electromagnetic instabilities in tokamaks at experimentally relevant values of plasma β⁠. This scheme resolved the long-standing cancellation problem that previously hindered gyrokinetic particle-in-cell code simulations of magnetohydrodynamic phenomena with inherently small parallel electric fields. Moreover, the scheme did not employ approximations that eliminate critical tearing-type instabilities. Here, we report on the implementation of this numerical scheme in the global gyrokinetic particle-in-cell code GTS. This implementation allows for a more complete and accurate picture of interaction between small scale turbulence and MHD modes in tokamaks. Additionally, we present a comprehensive set of verification simulations of numerous electromagnetic instabilities relevant to present-day tokamaks. These simulations encompass the kinetic ballooning mode, the internal kink mode, the tearing mode, the micro-tearing mode, and the toroidal Alfven eigenmode destabilized by energetic ions, which are all instrumental in understanding tokamak physics. We will also showcase the preliminary nonlinear simulations of kinetic ballooning instabilities and (2,1) island formation due to tearing mode instability. These simulations validate the accuracy of the scheme implementation and pave the way for studying how these instabilities affect plasma confinement and performance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Electron acceleration in solar flares and the transition from nonthermal to thermal hard X-ray phases

Observations are reviewed which indicate that hard X-rays during the impulsive phase of a flare typically start with a primarily nonthermal phase which undergoes a transition to a primarily thermal phase as the flare progresses. Recent theoretical work on the modified two-stream instability as an efficient electron accelerator and modeling of thermal hard X-ray sources is considered. A scenario which is termed the dissipative thermal model is proposed to explain the observations. Fast tearing modes occurring in a loop give rise to cross-field ion motion. This in turn excites the modified two-stream instability which converts about 50 percent of the ion energy into accelerated electrons along the loop as long as the plasma beta is less than 0.3. These electrons impact the chromosphere and boil off a part of it which rises up the loop. This density increase coupled with the temperature increase due to tearing causes the beta to increase beyond 0.3 and efficient electron acceleration ceases. This leads to the primarily thermal phase.

Smith, D. F.↗
DOE OSTI2024

Performance prediction applying different reduced turbulence models to the SMART tokamak

The SMall Aspect Ratio Tokamak (SMART) is currently being commissioned at the University of Seville and will be able to compare the performance of positive and negative triangularity plasmas at low aspect ratio. Predictive simulations have been performed for different machine scenarios and heating schemes using the TRANSP code. The objectives of these simulations are to predict the parameters expected in positive triangularity plasmas, to guide diagnostic development, and to validate transport models. Several reduced turbulence models have been used to predict electron and ion temperatures for the operational phase 2. All models provide similar results from approximately mid-radius to the separatrix but important discrepancies are found in the core region. These positive triangularity results are compared with experiments from a similar size machine like GLOBUS-M2. The multi-mode model (MMM) shows the best agreement. Simulations with different boundary conditions have been performed and no strong differences have been observed between them. The impact of neutral beam injection (NBI) on the predicted profiles has also been addressed. Rotation reduces turbulence levels so higher temperatures are achieved when included in the simulations. Studying the different contributions to the thermal diffusivities, it is observed that electron temperature gradient (ETG) turbulence dominates at the plasma core while micro-tearing modes (MTM) dominate at the edge in the electron channel. In the ion channel, the neoclassical contribution is dominant at the core and at the very edge while the Weiland component, which includes ion temperature gradient mode (ITG), trapped electron mode (TEM), kinetic ballooning mode (KBM), peeling mode (PM) and collisionless and collision dominated magnetohydrodynamic (MHD) modes governs the mid-radius region. For phase 3, two plasmas with different electron densities have been studied. The case with lower density matches well a specific discharge of GLOBUS-M2. The higher density plasma shows high performance with β N ≈ 3.8.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗