Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETIC ISLANDS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

First demonstration of radiative power exhaust with impurity seeding in the island divertor at Wendelstein 7-X

Radiative power exhaust by impurity seeding was demonstrated for the first time in island divertor configurations at the stellarator Wendelstein 7-X. Feasibility of stable plasma operation was shown during seeding with both neon (Ne) and nitrogen (N 2 ). High radiative power losses (≥ 80%) were found to reduce the divertor heat loads globally by 2/3 with both seeding gases injected at a single toroidal location into one of five magnetic islands. Heat flux detachment was achieved for the price of a loss of (≥ -15%) in the stored energy. Ne seeding allows for sustained enhancement of edge radiation with a very slow decay of line emission of several tens of seconds after the end of the injection indicating a high recycling of this noble gas at the carbon main plasma facing components. In N 2 seeded discharges it is shown that a response of line emission and plasma parameters is in correlation to the puff duration which indicates a higher level of absorption of this seeding gas in the wall. Continuous N 2 seeding results in global cooling of the scrape-off layer (SOL) and decay of radiation over several seconds after the injection. Damping of counter-streaming SOL flows, and divertor particle fluxes induced by Ne and N 2 seeding have been measured and provide evidence for a reduction of the convective part of the divertor heat fluxes. Losses in density in response to seeding can be compensated by feedback controlled divertor fueling. The controlled reduction of heat fluxes within this complex 3D edge island geometry is a very promising finding concerning detachment control in a future all-metal divertor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2009

The Time-Dependent Structure of the Electron Reconnection Layer

Collisionless magnetic reconnection is often associated with time-dependent behavior. Specifically, current layers in the diffusion region can become unstable to tearing-type instabilities on one hand, or to instabilities with current-aligned wave vectors on the other. In the former case, the growth of tearing instabilities typically leads to the production of magnetic islands, which potentially provide feedback on the reconnection process itself, as well as on the rate of reconnection. The second class of instabilities tend to modulate the current layer along the direction of the current flow, for instance generating kink-type perturbations, or smaller-scale turbulence with the potential to broaden the current layer. All of these processes contribute to rendering magnetic reconnection time-dependent. In this presentation, we will provide a summary of these effects, and a discussion of how much they contribute to the overall magnetic reconnection rate.

Hesse, Michael↗
DOE OSTI2021

Skyrmions in an oblique field

Here we study skyrmion states induced by the oblique field in a square magnetic island. To describe stable states and dynamics of the skyrmion, we employ the lattice model that accounts for all relevant interactions and uses parameters in the range of practical interest. Depending on the orientation of the field, two or four spatially separated energy minima emerge in the oblique field. The depth of the energy minima can be controlled by the strength and orientation of the magnetic field. We study the dynamics of the skyrmion and show that it can be moved between any two states by the field gradient. Possible applications for developing binary and quaternary skyrmion-based memory are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The General Atomics Center for Magnetohydrodynamic Modeling (CEMM) Final Report

The Center for Extended Magnetohydrodynamic Modeling (CEMM) is built around the further development, verification, validation, and other application of two major extended magnetohydrodynamics (MHD) codes NIMROD and M3D-C1. The General Atomics (GA) CEMM Program focuses on three main research areas: extension of the M3D-C1 MHD code plasma response calculations to new regimes, coupling of the DK4D drift kinetic equation solver to a temperature evolution equation, and simulations of disruption mitigation for unstable plasmas with pre-existing magnetic islands using massive gas injection (MGI) with the NIMROD MHD code. A systematic scan of the zero-crossing of the rotation profile in a DIII-D ITER-similar shape equilibrium using linear, time-independent modeling with M3D-C1 was performed. A new Python utility, autoC1, has been developed to automate the running of M3D-C1 for linear simulations. A number of ELM suppression experiments were analyzed including a validation and cross-code verification with the MARS-F MHD code of magnetic signal phase shifts across ELM-suppression bifurcations. A linearized, time-dependent temperature equation was tightly coupled to the electron and ion drift-kinetic equations (DKEs) in the DK4D drift-kinetic equation solver. Various explicit and implicit implementations were attempted, but all methods proved numerically unstable, even with the addition of numerical dissipation. Several simulations were completed investigating the effects of pre-existing islands on disruption mitigation by MGI in DIII-D using NIMROD.

3D MHD Simulation↗
DOE OSTI2020

Block–spiral magnetism: An exotic type of frustrated order

Significance Magnetic frustration in spin or electronic models emerges from the failure of the system to fulfill simultaneously conflicting local requirements. The latter typically arise from the lattice geometry or are induced by special spin–spin interactions in the system at various distances. Here we show that the competing energy scales of the seemingly nonfrustrated orbital-selective Mott phase of the low-dimensional multiorbital Hubbard model can originate a “block–chiral magnetism,” i.e., a state with rigidly rotating spin–magnetic islands. Furthermore, we show how such an exotic spin state influences the electronic properties of the system, revealing parity-breaking quasiparticles.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Pedestal collapse by resonant magnetic perturbations

Pedestal collapse (i.e., the complete loss of the edge transport barrier (ETB)) in DIII-D H-mode plasmas occurs when resonant magnetic perturbations (RMPs) penetrate the steep gradient region at the plasma edge. Normally, RMP driven magnetic islands can occur at the top and bottom of the H-mode pedestal and these islands generate conditions consistent with edge-localized-mode (ELM) suppression and density pump-out, respectively, based on nonlinear two-fluid MHD simulations. In contrast, MHD simulations show that the steep pressure gradient region between the top and bottom of the DIII-D pedestal is generally immune to resonant field penetration due to large local E × B and diamagnetic flows. By this fortuitous circumstance, the edge-transport-barrier and H-mode confinement can be maintained while achieving ELM suppression. However, pedestal collapse can occur in DIII-D when the screening flows are inadequate to prevent field penetration in the steep gradient region of the pedestal. Non-linear two-fluid MHD simulations support the role of resonant field penetration in pedestal collapse for DIII-D H-mode plasmas with weak edge E × B and diamagnetic screening flows. ITER will likely have weaker edge screening flows than present experiments due to its much larger size, making it more susceptible to resonant field penetration in the steep gradient region of the pedestal. Analysis of model ITER equilibria demonstrates that resonant field penetration in the steep pressure gradient region is possible for RMP levels of the order required for ELM suppression. The effect of such penetration on the ITER pedestal will depend sensitively on the resulting degree of island overlap.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Onset of Plasmoid Reconnection during Magnetorotational Instability

The evolution of current sheets in accretion flows undergoing magnetorotational instability (MRI) is examined through two- and three-dimensional numerical modeling of the resistive MHD equations in global cylindrical geometry. With an initial uniform magnetic field aligned in the vertical (z) direction, MRI produces radially extended toroidal (azimuthal) current sheets. In both 2D and 3D when axisymmetric modes dominate, these current sheets attract each other and merge in the poloidal (rz) plane, driving magnetic reconnection when the Lundquist number S > 3 × 10 2 , making it a possible source of plasmoids (closed magnetic loops) in accretion disks. At high Lundquist numbers in the 2D regime, starting at S = 5 × 10 3 , self-consistent MRI-generated current sheets become thin and subject to plasmoid instability, and therefore spontaneous magnetic reconnection. When nonaxisymmetric 3D modes dominate, turbulence makes the azimuthal current sheets more unstable and stretch vertically. Toroidally extended vertical current sheets in the inner region, as well as larger 3D magnetic islands in the outer regions of the disks are also formed. These findings have strong ramifications for astrophysical disks as potential sources of plasmoids that could cause local heating, particle acceleration, and high energy EM radiation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Mode penetration induced by transient magnetic perturbations

Three dimensional magnetic fields in tokamaks can induce forced magnetic reconnection (FMR) and produce magnetic islands on resonant surfaces. Conventional analytic solutions to FMR focus on describing the time asymptotic state given a steady-state field error. The purpose of this report is to understand the nonlinear dynamics of mode penetration, an evolution from a high-slip, flow-screened metastable equilibrium into a low-slip, field-penetrated metastable equilibrium. Here, we extend previous work by incorporating a temporally varying external magnetic field as a simple model for a magnetohydrodynamic (MHD) event that produces resonant magnetic perturbations. Proof-of-principle, extended-MHD, NIMROD computations vary parameterizations of the transient external perturbation to probe the threshold for mode penetration. We test these computational results against analytical theory that captures the temporal evolution properties of the electromagnetic and viscous forces during and after a transient. We find qualitative agreement between computational and analytical results. However, computational tools are necessary to accurately capture the threshold conditions for mode penetration induced by an MHD transient.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Numerical study of $\delta$-function current sheets arising from resonant magnetic perturbations

We report general three-dimensional toroidal ideal magnetohydrodynamic equilibria with a continuum of nested flux surfaces are susceptible to forming singular current sheets when resonant perturbations are applied. The presence of singular current sheets indicates that, in the presence of non-zero resistivity, magnetic reconnection will ensue, leading to the formation of magnetic islands and potentially regions of stochastic field lines when islands overlap. Numerically resolving singular current sheets in the ideal magnetohydrodynamics (MHD) limit has been a significant challenge. This work presents numerical solutions of the Hahm–Kulsrud–Taylor (HKT) problem, which is a prototype for resonant singular current sheet formation. The HKT problem is solved by two codes: a Grad–Shafranov (GS) solver and the Stepped Pressure Equilibrium Code (SPEC) code. The GS solver has built-in nested flux surfaces with prescribed magnetic fluxes. The SPEC code implements multi-region relaxed magnetohydrodynamics (MRxMHD), whereby the solution relaxes to a Taylor state in each region while maintaining force balance across the interfaces between regions. As the number of regions increases, the MRxMHD solution appears to approach the ideal MHD solution assuming a continuum of nested flux surfaces. We demonstrate agreement between the numerical solutions obtained from the two codes through a convergence study.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Access to stable, high pressure tokamak pedestals using local electron cyclotron current drive

Abstract Experiments on the DIII-D tokamak demonstrate that edge localized mode (ELM) stability can be manipulated using localized electron cyclotron current drive (ECCD) in conjunction with resonant magnetic perturbations (RMPs). The injection of counter-plasma-current edge ECCD reduces the RMP amplitude required to suppress ELMs and bifurcates the pedestal into a high-confinement regime with 7 kPA pedestal pressure. This is the first time such a high confinement regime has been accessed through the bifurcation from the ballooning stability branch predicted by existing models. These observations are consistent with modeled ECCD manipulation of magnetic islands induced by the RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Effect of toroidal plasma currents on the Wendelstein 7-X Scrape-Off Layer

The role of toroidal plasma currents for the island divertor scrape-off layer in the stellarator Wendelstein 7-X is investigated using reciprocating electric probes. Experiments show that small amounts (of a few kA) of plasma current are sufficient to significantly affect the scrape-off layer plasma conditions, whereas higher plasma currents above 10kA result in more drastic changes. This behavior is linked to the effect of the plasma current on the rotational transform profile, which can result in significant shifts of the edge magnetic islands. These shifts affect the interaction of the islands with the divertor and can eventually result in a transition from a diverted to a limited plasma configuration. The probe observations are complemented by further edge diagnostics including plasma flow measurements, divertor Langmuir probes, divertor thermography and impurity spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A divergence-free constrained magnetic field interpolation method for scattered data

An interpolation method to evaluate magnetic fields, given its unstructured and scattered magnetic data, is presented. The method is based on the reconstruction of the global magnetic field using a superposition of orthogonal functions. The coefficients of the expansion are obtained by minimizing a cost function defined as the L2 norm of the difference between the ground truth and the reconstructed magnetic field evaluated on the training data. The divergence-free condition is incorporated as a constraint in the cost function, allowing the method to achieve arbitrarily small errors in the magnetic field divergence. An exponential decay of the approximation error is observed and compared with the less favorable algebraic decay of local splines. Compared to local methods involving computationally expensive search algorithms, the proposed method exhibits a significant reduction of the computational complexity of the field evaluation, while maintaining a small error in the divergence even in the presence of magnetic islands and stochasticity. Finally, applications to the computation of Poincaré sections using data obtained from numerical solutions of the magnetohydrodynamic equations in toroidal geometry are presented and compared with local methods currently in use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of hydrogen fueling, recycling, and confinement at Wendelstein 7-X via a single-reservoir particle balance

Abstract A single-reservoir particle balance for the main plasma species hydrogen has been established for W7-X. This has enabled the quantitative characterization of the particle sources in the standard island divertor configuration for the first time. Findings from attached scenarios with two different island sizes with a boronized wall and turbo molecular pumping are presented. Fueling efficiencies, particle flows and source locations were measured and used to infer the total particle confinement time τ p . Perturbative gas injection experiments served to measure the effective particle confinement time τ p * . Combining both confinement times provides access to the global recycling coefficient R ¯ . Hydrogen particle inventories have been addressed and the knowledge of particle sources and sinks reveals the core fueling distribution and provides insight into the capability of the magnetic islands to control exhaust features. Measurements of hydrogen fueling efficiencies were sensitive to the precise fueling location and measured between 12% and 31% with the recycling fueling at the strike line modeled at only 6%, due to much higher densities. 15% of the total 5.2 × 10 22 a/s recycling flow ionizes far away from the recycling surfaces in the main chamber. It was shown that 60% of recycled particles ionize above the horizontal and 18% above the vertical divertor target, while the remainder of the recycling flow ionizes above the baffle (7%). Combining these source terms with their individual fueling efficiencies resolves the core fueling distribution. Due to the higher fueling efficiency in the main chamber, up to 51% of the total 5.1 × 10 21 s −1 core fueling particles are entering the confined plasma from the main chamber. τ p values in the range of 260 ms were extracted for these discharges. Together with τ p , the global recycling coefficient R ¯ was resolved for every τ p * measurement and a typical value close to unity was obtained. An increase of the island size, resulted in no change of τ p , but doubled τ p * , indicating the feasibility of the control coils as an actuator to control exhaust features without affecting core confinement properties.

recycling↗
DOE OSTI2021

Plasma–surface interaction in the stellarator W7-X: conclusions drawn from operation with graphite plasma-facing components

W7-X completed its plasma operation in hydrogen with island divertor and inertially cooled test divertor unit (TDU) made of graphite. A substantial set of plasma-facing components (PFCs), including in particular marker target elements, were extracted from the W7-X vessel and analysed post-mortem. The analysis provided key information about underlying plasma–surface interactions (PSI) processes, namely erosion, transport, and deposition as well as fuel retention in the graphite components. The net carbon (C) erosion and deposition distribution on the horizontal target (HT) and vertical target (VT) plates were quantified and related to the plasma time in standard divertor configuration with edge transform ι = 5/5, the dominant magnetic configuration of the two operational phases (OP) with TDU. The operation resulted in integrated high net C erosion rate of 2.8 mg s -1 in OP1.2B over 4809 plasma seconds. Boronisations reduced the net erosion on the HT by about a factor 5.4 with respect to OP1.2A owing to the suppression of oxygen (O). In the case of the VT, high peak net C erosion of 11μm at the strike line was measured during OP1.2B which converts to 2.5 nm s -1 or 1.4 mg s -1 when related to the exposed area of the target plate and the operational time in standard divertor configuration. PSI modelling with ERO2.0 and WallDYN-3D is applied in an interpretative manner and reproduces the net C erosion and deposition pattern at the target plates determined by different post-mortem analysis techniques. This includes also the 13 C tracer deposition from the last experiment of OP1.2B with local 13 CH 4 injection through a magnetic island in one half module. The experimental findings are used to predict the C erosion, transport, and deposition in the next campaigns aiming in long-pulse operation up to 1800 s and utilising the actively cooled carbon-fibre composite (CFC) divertor currently being installed. The CFC divertor has the same geometrical design as the TDU and extrapolation depends mainly on the applied plasma boundary. Extrapolation from campaign averaged information obtained in OP1.2B reveals a net erosion of 7.6 g per 1800 s for a typical W7-X attached divertor plasma in hydrogen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence

Here using 2.5 dimensional kinetic particle-in-cell simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasisteady reconnection rates. Note that in these simulations, the ion Larmor radius is comparable to the ion inertial length. Moreover, the transition to a more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasisteady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with the exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfvén speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments, we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Electron cyclotron emission detection of neoclassical tearing modes for control for ITER

Successful operation of ITER requires control of magnetic instabilities including neoclassical tearing modes (NTMs) that can degrade confinement and lead to disruption. Low latency detection by electron cyclotron emission (ECE) diagnostics has been demonstrated in a few current experiments. Using a synthetic diagnostic, we demonstrate low latency NTM detection for ITER with plasmas described by ITER IMAS database scenarios and with realistic limitations imposed on the instrumentation by these high temperature scenarios. 2/1 NTMs are detected 430 ms after magnetic island seeding and before island locking. The radiometer configuration was optimized using simulation, and the smallest detectable island size was explored. Island sizes of ∼3 cm are detectable at the 2/1 surface. The simulated signals incorporate recent physics models for island growth and rotation, which show early locking and continued island growth after locking and before disruption. This work determines limits for ITER ECE spatial resolution imposed by relativistic broadening of channels, which informs hardware design. Real-time detection is demonstrated in hardware that is required by ITER, including on an NI PXI-7853R FPGA system. Development of a synthetic diagnostic and details of the hardware will be discussed.

Cyclotron radiation↗
DOE OSTI2021

Passive deconfinement of runaway electrons using an in-vessel helical coil

A helical coil designed to passively generate non-axisymmetric fields during a plasma disruption is shown (via electromagnetic analysis, linear MHD modeling, and relativistic drift orbit tracing) to be effective at deconfining runaway electrons (REs) on a time scale significantly faster than the plasma current quench. Magnetic equilibria from DIII-D RE-producing scenarios are used to calculate the toroidal electric field generated during the current quench phase of a disruption, which in turn drives current in the proposed n = 1 in-vessel helical coil, without the need for any external power supplies or disruption detection or prediction techniques. Simulations of the plasma evolution using the TokSys GS Evolve code predict the inductive coupling of coil currents up to 12% of the pre-disruption plasma current into the helical coil. The coil geometry is parametrically varied to maximize both the non-resonant and resonant components of the 3D magnetic perturbation, resulting in δB/B ≈ 10 –2 and a vacuumisland overlapwidth of up to 0.7ψ N . The REORBIT module of the MARS-F code is used to model the full non-axisymmetric magnetic field and trace RE drift orbits to determine the effect on RE deconfinement, with up to 70% of the RE orbits lost after 0.2 ms. A two-stage evolution of the RE orbit loss fraction is observed to be caused by resonant trapping between multiple magnetic island chains. Finally, electromagnetic and thermal stresses on the coil are calculated to be within operational limits for installation in DIII-D, and scale favorably to a reactor-size device. Furthermore, these findings motivate future experimental study of the helical coil concept in DIII-D or other tokamaks.

3D coil↗