Search NASA⌕ Search

SEARCH · Search NASA

Results for “islanded mode”

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 127 records · Page 7

Global gyrokinetic simulations of the impact of magnetic island on ion temperature gradient driven turbulence

The effect of island width on the multi-scale interactions between magnetic island (MI) and ion temperature gradient (ITG) turbulence has been investigated based on the global gyrokinetic approach. It is found that the coupling between the island and turbulence is enhanced when the MI width (w) becomes larger. A vortex flow that is highly sensitive to the width of the MI can be triggered, ultimately resulting in a potent E × B shear flow and a consequent reduction in turbulent transport. The shearing rate induced by the vortex flow is minimum at the O-point while it is maximum at both of the two reconnection points of the island, i.e. the X-points, regardless of the island width. There exists a nonmonotonic relationship between zonal flow (ZF) amplitude and island width, showing that the ZF is partially suppressed by medium-sized MIs whereas enhanced in the case of large island. A larger MI can tremendously damage the ITG mode structure, resulting in higher turbulent transport at the X-point whereas a lower one at the O-point, respectively. Such phenomenon will be less distinct at very small island widths below w/a ~ 8% ≈ 12ρ i (a is the minor radius and ρ i the ion gyroradius), where it shows that turbulence near the X-point is hardly affected although it is still suppressed inside the island. Furthermore, the influence of different island sizes on turbulence transport level is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling Solar Flare Hard X-ray Images and Spectra Observed with RHESSI

Observations obtained with the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of a flare on February 20, 2002 indicate a hard X-ray (HXR) coronal source at or near the top of a flare loop (called a HXR looptop source). The existence of the HXR looptop source suggests that magnetic reconnection, which is believed to power flares, occurs above the loop. In order to explain this HXR looptop source, I created a steady-state particle transport model, in which high-energy electrons are continuously injected at the top of a semicircular flare loop. Based on the simulation results, I find that the model predictions are consistent with the RHESSI observations in many respects, but the spectrum of the looptop source obtained from the model is steeper than that from the RHESSI data. This suggests that, instead of being accelerated above the loop as generally believed, the particles might be accelerated in the looptop itself. RHESSI-observations of three other homologous flares that occurred between April 14 and 16, 2002, provide strong evidence for the presence of a large scale current sheet above a flare loop, which is the basis of standard flare models. The most convincing finding is the presence of the temperature distribution of a separate coronal source above the flare loops: the hotter part of the coronal source was located lower in altitude the cooler part. Together with the fact that the hotter flare loops are higher than the cooler loops, the observations support the existence of a large-scale current sheet between the top of the flare loops and the coronal source above. Blob-Like sources along a line above the loop in the decay phase of the April 15, 2002, flare, which are suggestive of magnetic islands initiated by the tearing-mode instability, and the observation of a cusp structure in microwaves, further support the presence of the current sheet. The observations of the three homologous flares reveal two other features which are beyond the predictions of the standard flare models: the downward motion of flare loops in the early impulsive phase of each flare, and an initially stationary coronal source above the loops. These features me believed to be related to the formation and development of a current sheet. In particular, the downward loop motion seem to be a common phenomenon in flares, suggesting the necessity for modifications to the existing standard flare. models. Finally, thanks to the broad energy coverage of the RHESSI spectra, a low-energy cutoff of 28 (plus or minus 2) keV in the nonthermal electron distribution was determined for the April 15, 2002, flare. As a result, the energy carried by the nonthermal electrons is found to be comparable to the thermal energy of the flare, but one order of magnitude larger than the kinetic energy of the associated coronal mass ejection. The method used to deduce the electron low-energy cutoff will be useful in the analyses of similar events.

Sui, Linhui↗

Macroscopic trends of neoclassical tearing stability in high-field H-mode tokamak pilot plants

The neoclassical tearing mode (NTM) stability metric—minimum marginally stable island width $w$$^{*}_{m}$—was compared across 14651 inductive high-field tokamak pilot plant equilibria. Larger devices with reduced elongation and/or increased minor radius demonstrated an order-of-magnitude increase in $w$$^{*}_{m}$, primarily due to a reduction in bootstrap drive. This work is part of an ongoing effort to ensure passive NTM-stability in the ARC tokamak, in which the technology to achieve active tearing-suppression with localised electron cyclotron current drive does not yet exist. The equilibrium scenarios in the database were Monte Carlo generated and normalised to the same >400MW fusion power, minimum pressure scenario at a range of plasma currents, before tearing analysis using the modified Rutherford equation was applied for all resonant poloidal and toroidal m, n modes up to n = 4. Single-helicity toroidal Δ' calculations in resistive DCON set the minimum marginally stable island width, and a simple modal scaling proportional to –m 2 n –1 was identified for high-m Δ' values. The dominant correlates of $w$$^{*}_{m}$ and Δ' across the database were analysed using interpretable machine learning techniques.

NTM seeding↗

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↗

The role of edge resonant magnetic perturbations in edge-localized-mode suppression and density pump-out in low-collisionality DIII-D plasmas

Two-fluid nonlinear MHD simulations using the TM1 code demonstrate that the formation of magnetic islands at the top and bottom of the H-mode pedestal, together with the strong screening of resonant fields in the gradient region of the pedestal, can account for ELM suppression and density pump-out by n = 2 resonant magnetic perturbations (RMPs) in low-collisionality DIII-D ITER similar shape (ISS) plasmas. Using experimentally relevant transport coefficients, neoclassical resistivity, electron collisionality, and RMP amplitudes, nonlinear MHD simulations reproduce the observed level of density reduction (density pump-out) in DIII-D due to the formation of narrow magnetic islands in the resistive foot of pedestal. For large amplitude RMPs ( B r / B t > 1 × 10 -4 ) simulations predict resonant field penetration and resulting in pressure reduction at the top of the pedestal, consistent with experimental observations at the onset of ELM suppression. The predicted reduction in the height and width of the pedestal by magnetic island enhanced collisional transport provides a quantitative mechanism for the stabilization of the peeling-ballooning mode (PBM). Importantly, these simulations predict strong screening of resonant fields in the steep gradient region of the pedestal due to strong E × B rotation and diamagnetic flows. Yet, if the plasma resistivity is made artificially larger (~ 10X) than neoclassical, the simulations predict magnetic stochasticity throughout the plasma edge and the collapse of the pedestal due to the reduction in the penetration threshold with increasing resistivity. A scaling relation for resonant field penetration at the pedestal top, using several hundred nonlinear simulations, reproduces the density and E × B rotation dependence of the ELM suppression threshold observed in DIII-D. Simulations using ITER model equilibria indicate that the penetration threshold at the top of the ITER pedestal will be lower than in DIII-D due to the anticipated lower perpendicular flow velocities in ITER, resulting in the weaker screening of resonant fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neoclassical Tearing Mode Seeding by Nonlinear Three-Wave Interactions in Tokamaks

We report the experimental observation of seed magnetic island formation by nonlinear three-wave coupling of magnetic island triplets. In this experiment, disruptive 2,1 islands are seeded by the coupling of 4,3 and 3,2 tearing modes to a central 1,1 sawtooth precursor. Three-wave interactions between these modes are conclusively identified by bispectral analysis, indicating fixed phase relationships in agreement with theory. Furthermore, this new observation of this seeding mechanism has important implications for future reactors that must operate in stable plasma equilibria, free of disruptive 2,1 islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the stabilisation of locked tearing modes in ITER and other large tokamaks

Tearing modes in tokamaks typically rotate while small and then lock at a fixed location when larger. Research on present-day devices has focused almost exclusively on stabilisation of rotating modes, as it has been considered imperative to avoid locked modes. However, in larger devices, stabilisation during the rotating phase is made difficult by fast locking at small island widths, and large broadening of the stabilising wave-driven current profile. In contrast, the smaller island width at locking not only mitigates the deleterious consequences of locked modes, but also permits their efficient stabilisation. On large devices, it thus becomes surprisingly advantageous to allow the mode to grow and lock naturally before stabilising it, challenging the mainstream strategy of neoclassical tearing mode stabilisation during the rotating phase. Furthermore, calculations indicate that a locked island stabilisation strategy should be adopted in the ITER tokamak, with a large potential impact on the fusion gain and disruptivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

Experimental evidence of electron-cyclotron current drive-based neoclassical tearing mode suppression threshold reduction during mode locking on DIII-D

The repetitive locking-unlocking process of a large $m/n = 2/1$ neoclassical tearing mode (NTM) is studied on the DIII-D tokamak. Cyclical periods of partial suppression are achieved when β p and the expected bootstrap current $J_{\mathrm{BS}}$ drop upon back-transitions to L-mode. Diminished L-mode bootstrap current is required to explain decreased island size as modeled by the generalized Rutherford equation. At EC power levels too low to fully suppress the NTM, partial suppression is enough to reduce the torque on the mode to the point where viscous drag from the background plasma overcomes the electromagnetic torque and restores fast plasma rotation and H-mode. Additional current drive is shown to further reduce the size of the mode until full suppression is achieved. These results suggest that the suppression of locked modes may be aided by the lower bootstrap current in L-mode plasmas, possibly improving chances of disruption avoidance in ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Islands far outside the horizon

Information located in an entanglement island in semiclassical gravity can be nonperturbatively reconstructed from distant radiation, implying a radical breakdown of effective field theory. We show that this occurs well outside of the black hole stretched horizon. We compute the island associated to large-angular momentum Hawking modes of a four-dimensional Schwarzschild black hole. These modes typically fall back into the black hole but can be extracted to infinity by relativistic strings or, more abstractly, by asymptotic boundary operators constructed using the timelike tube theorem. Remarkably, we find that their island can protrude a distance of order $\sqrt{\ell_p{r}_{\textrm{hor}}}$ outside the horizon. This is parametrically larger than the Planck scale ℓ p and is comparable to the Bohr radius for supermassive black holes. Therefore, in principle, a distant observer can determine experimentally whether the black hole information paradox is resolved by complementarity, or by a firewall.

AdS-CFT correspondence↗

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

Error field measurements with rotating RMP fields for DIII-D H-mode

3D magnetic sensors are employed to identify the amplitude and toroidal phase of error fields (EF) by analyzing the torque balance for magnetic islands entrained by rotating resonant magnetic perturbations (RMPs) in DIII-D H-mode plasmas. This technique of torque balance allows for efficient error field identification, offering a valuable tool for scenario-specific and optimized error field compensation (EFC) and requires only magnetic diagnostics. The torque balance used in this work includes the contribution from electromagnetic torque due to error fields, wall response, RMP fields, and viscous. Results show that viscous torque plays a crucial role, particularly during locked modes and H-mode plasmas, ensuring accurate data fits with lower residuals. The torque balance technique reveals that the L- and H-mode plasmas have distinct EF configurations, and consistent with a model-based EF analysis including MHD response in IPEC and the SURFMN EF simulation. This technique shows great robustness in measuring the intrinsic EF amplitude regardless of its amplitude or toroidal phase. Repeated discharges with EFC disparities exhibit consistent results of intrinsic error field within a reasonable range near the “standard” error field compensation. Additionally, the use of a rotating n = 1 resonant magnetic perturbation offers the advantage of reducing disruption risks by entraining saturated magnetic islands. These findings are instrumental for optimizing EF correction in fusion devices, thereby enhancing tearing mode suppression and overall plasma stability.

3D magnetic sensors↗

Predicting operational windows of ELMs suppression by resonant magnetic perturbations in the DIII-D and KSTAR tokamaks

A newly developed plasma response model, combining the nonlinear two-fluid MHD code TM1 and toroidal MHD code GPEC run in ideal mode, quantitatively predicts the narrow isolated q 95 windows (Δq 95 ~ 0.1) of ELM suppression by n = 1, 2 and 3 resonant magnetic perturbations (RMPs) in both DIII-D and KSTAR tokamaks across a wide range of plasma parameters. The key physics that unites both experimental observations and our simulations is the close alignment of essential resonant q-surfaces and the location of the top of the pedestal prior to an ELM. This alignment permits an applied RMP to produce field penetration due to the lower E × B rotation at the pedestal top rather than being screened. The model successfully predicts that narrow magnetic islands form when resonant field penetration occurs at the top of pedestal, and these islands are easily screened when q 95 moves off resonance, leading to very narrow windows of ELM suppression (typically Δq 95 ~ 0.1). Furthermore, the observed reduction in the pedestal height is also well captured by the calculated classical collisional transport across the island. We recover observed q 95 , β N and plasma shape dependence of ELM suppression due to the effect of magnetic islands on pedestal transport and Peeling-Ballooning- Mode (PBM) stability. Importantly, experiments do occasionally observe wide windows of ELM suppression (Δq 95 > 0.5). Our model reveals that at low pedestal-top density multiple islands open, leading to wide operational windows of ELM suppression consistent with experiment. The model indicates that wide q95 windows of ELM suppression can be achieved at substantially higher pedestal pressure with less confinement degradation in DIII-D by operating at higher toroidal mode number (n = 4) RMPs. As a result, this can have significant implications for the operation of the ITER ELM control coils for maintaining high confinement together with ELM suppression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tri-Level Linear Programming Model for Automatic Load Shedding Using Spectral Clustering

Traditional load shedding schemes can be inadequate in grids with high renewable penetration, leading to unstable events and unnecessary grid islanding. Although for both manual and automatic operating modes load shedding areas have been predefined by grid operators, they have remained fixed, and may be sub-optimal due to dynamic operating conditions. In this work, a distributed tri-level linear programming model for automatic load shedding to avoid system islanding is presented. Preventing islanding is preferred because it reduces the need for additional load shedding besides the disconnection of transmission lines between islands. This is crucial as maintaining the local generation-demand balance is necessary to preserve frequency stability. Furthermore, uneven distribution of generation resources among islands can lead to increased load shedding, causing economic and reliability challenges. This issue is further compounded in modern power systems heavily dependent on non-dispatchable resources like wind and solar. The upper-level model uses complex power flow measurements to determine the system areas to shed load depending on actual operating conditions using a spectral clustering approach. The mid-level model estimates the area system state, while the lower-level model determines the locations and load values to be shed. The solution is practical and promising for real-world applications.

Baquedano-Aguilar, Mario D.↗

Drag exerted by drift-acoustic radiation on a magnetic island

The Island Equilibrium and Transport (IslET) model is extended to allow for the effects of ion compressibility albeit for a constant electron temperature. The model describes fully relaxed, saturated, steady-state magnetic islands arising from the evolution of a tearing mode caused either by an instability or an externally imposed resonant perturbation. It assumes that the collision frequency lies in the semi-collisional regime. A salient new feature of the extended model is an equation describing the transport of parallel momentum. The model is applied to evaluate the drag exerted on a thin island by the emission of drift-acoustic waves.

Physics↗

Epitaxial columnar growth of strain-free antiferromagnetic Weyl semimetal Mn 3 Sn on wurtzite c -plane GaN/Al 2 O 3 (0001)

Weyl semimetal thin films with excellent crystalline quality are of great interest for antiferromagnetic spintronics. Mn 3 Sn is one Weyl semimetal with great properties and promise for exciting science and applications. It has proven very challenging, however, to grow Mn 3 Sn thin films with smooth surfaces, negligible strain, and excellent crystallinity. In this work, we discuss the successful preparation of epitaxial Mn 3 Sn (0001)-oriented thin films via molecular beam epitaxial growth on c -plane wurtzite GaN which was grown by MBE on Al 2 O 3 (0001). We present the reflection high energy electron diffraction analysis along with x-ray diffraction in order to demonstrate the crystalline quality of the film, and we give atomic models to explain the epitaxial orientation relationships between the crystal lattices of the substrate, GaN layer, and Mn 3 Sn layer. Importantly, we discuss the film lattice parameters as compared to expected values, demonstrating negligible strain both in-plane and out-of-plane . Atomic force microscopy reveals an epitaxial columnar growth mode characterized by flat-top-mesa islands, while scanning tunneling microscopy shows the atomically smooth surfaces of the mesa-top structures. Finally, Rutherford backscattering informs the stoichiometry of the film as well as the layer thicknesses.

Chiral antiferromagnetic material↗

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↗