Search NASA⌕ Search

SEARCH · Search NASA

Results for “edge 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

MARS-F/K modeling of plasma response and fast ion losses due to RMP in KSTAR

Abstract The toroidal single-fluid magnetohydrodynamic (MHD) code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) and the MHD-kinetic hybrid code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503) are utilized to study the plasma response to the n = 1 ( n is the toroidal mode number) resonant magnetic perturbation (RMP), applied to suppress the type-I edge localized mode (ELM) in a KSTAR discharge. Both the resistive-rotating and ideal-static plasma models identify strong screening of the resonant radial field harmonics of the applied RMP due to the plasma response, and predict a strong edge-peeling response of the plasma which is consistent with the optimal ELM control coil current configuration adopted in experiment. The RMP-induced radial displacement of the plasma, computed by the resistive-rotating plasma model, agrees reasonably well with that reconstructed from the measured data in the plasma core. Taking into account the drift kinetic response of fast ions, MARS-K hybrid modeling also finds quantitative agreement of the plasma core fluid pressure perturbation with experiment. Based on the MARS-F computed plasma response, a guiding-center orbit-tracing simulation finds about 0.3% of fast ion losses due the n = 1 RMP in the KSTAR ELM control experiment considered. Most losses are associated with counter-current fast ions located near the plasma edge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Toroidal modeling of runaway electron loss due to 3D fields in ITER

Mitigation of runaway electrons (REs) by three-dimensional (3D) magnetic field perturbations is numerically investigated for the ITER 15 MA baseline D–T scenario, utilizing the MARS-F code (Liu et al Phys. Plasmas 7 3681) with a drift orbit test particle tracing module. Considered are two types of 3D fields: the n = 3 (n is the toroidal mode number) resonant magnetic perturbation (RMP) utilized for the purpose of controlling the edge localized modes in ITER, and perturbations generated by the n = 1 magneto-hydrodynamic (MHD) instabilities in a post-disruption plasma. The RMP field, applied to a pre-disruption plasma, is found to be moderately effective in mitigating the RE seeds in ITER when vacuum field model is assumed. Up to ~40% loss fraction is possible at 90 kA-turn coil current. The mitigation efficiency is however substantially reduced, down to less than 5%, when the plasma response is taken into account. This is due to strong screening of the resonant magnetic field components by the plasma response resulting in much less field line stochasticity. On the other hand, the MARS-F modeling, based on the DINA-simulated post-disruption equilibria, shows that the n = 1 resistive kink instabilities develop in these plasmas, as the edge safety factor qa evolves and drops below integer numbers. RE mitigation by these MHD instabilities is sensitive to the eigenmode structure. The best mitigation is achieved as q a drops below 3, when a global kink instability occurs that encompasses both internal and external components. This global instability is found to be capable of mitigating over 80% MeV-level passing RE orbits at a field perturbation |δB|/B 0 that is comparable to that observed in DIII-D experiments, and full mitigation if the perturbation amplitude is doubled. The 'wetted' area on the ITER limiting surface, due to MHD instability induced RE loss, generally increases with the perturbation amplitude (together with increasing loss fraction). At the highest perturbation level assumed in this study, the wetted area reaches ~60% of the total limiting surface area. Lastly, the lost RE orbits mainly strike the outer divertor region of the limiting surface, with some fraction also hitting a wide area along the inboard side of the surface.

3D fields↗
DOE OSTI2024

Non-linear MHD modelling of transients in tokamaks: a review of recent advances with the JOREK code

Transient magneto-hydrodynamic (MHD) events like edge localized modes (ELMs) or disruptions are a concern for magnetic confinement fusion power plants. Research with the MHD code JOREK towards understanding control of such instabilities is reviewed here in a concise way to provide a complete overview, while we refer to the original publications for details. Experimental validation for unmitigated vertical displacement events progressed. The mechanism of vertical force mitigation by impurity injection was identified. Two-way eddy current coupling to CARIDDI was completed. Shattered pellet injection was simulated in JET, KSTAR, ASDEX Upgrade (AUG) and ITER. Benign runaway electron beam termination in JET and ITER was studied. Coupling of kinetic REs to the MHD is ongoing and a virtual RE synchrotron radiation diagnostic was developed. Regarding pedestal physics, regimes devoid of large ELMs in AUG were simulated and predictive JT60-SA simulations are ongoing. For ELM suppression by resonant magnetic perturbations (RMPs), AUG, ITER and EAST simulations were performed. A free boundary RMP model was validated against experiments. Evidence for penetrated magnetic islands at the pedestal top based on AUG experiments and simulations was found. Simulations of the naturally ELM-free quiescent H-mode in AUG and HL-3 show external kink mode formation prevents pedestal build-up towards an ELM within windows of the edge safety factor. With kinetic neutral particles, high field side high density formation in ITER was simulated and with kinetic impurities, tungsten transport in AUG RMP plasmas was studied. To capture turbulent transport, electro-static full-f particle in cell models for ion temperature gradient and trapped electron modes were established and benchmarked. Application to RMP plasmas shows enhanced turbulence in comparison to unperturbed states. Energetic particle interactions with MHD were studied. Flux pumping that prevents the safety factor on axis from dropping below unity was simulated. First non-linear stellarator applications include current relaxation in $l$ = 2 stellarators, while verification for advanced stellarators progresses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Understanding the formation of a low-pressure pedestal in the presence of a strong internal transport barrier in DIII-D high β p plasmas

As a promising scenario for fusion reactors, the high poloidal-beta (β P ) scenario is characterized by a strong large radius internal transport barrier (ITB), which significantly enhances the overall confinement quality and the bootstrap current fraction for fully non-inductive operation. It is frequently observed that in the presence of a strong ITB, the pedestal height is lower and is accompanied by small edge localized modes (ELMs), which further improves the compatibility of a high performance core with an edge solution. A mechanism for the formation of the low pedestal is proposed in this paper. It is found that the strong ITB creates an off-axis bootstrap current to clamp the local safety factor q, and thus the magnetic shear in the outer core/pedestal region is increased. Gyrokinetic simulations with the CGYRO code show that the higher magnetic shear brings the experimental profiles into the range where the growth rate of drift-wave instabilities and thus transport is higher, and therefore a lower pedestal gradient is expected. Here, the combination of low pedestal and high magnetic shear further enhances the turbulent transport across the whole pedestal, consistent with power balance analysis. Such a positive feedback mechanism ultimately results in a lower pressure pedestal as observed in experiments. Under such a low pedestal, linear simulations with BOUT++ predict the growth rates of peeling–ballooning modes to be lower across the whole toroidal mode number spectra, and the nonlinear BOUT++ simulation exhibits lower saturated fluctuation intensity as well, consistent with the experimentally observed lower ELM size.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Tailoring tokamak error fields to control plasma instabilities and transport

A tokamak relies on the axisymmetric magnetic fields to confine fusion plasmas and aims to deliver sustainable and clean energy. However, misalignments arise inevitably in the tokamak construction, leading to small asymmetries in the magnetic field known as error fields (EFs). The EFs have been a major concern in the tokamak approaches because small EFs, even less than 0.1%, can drive a plasma disruption. Meanwhile, the EFs in the tokamak can be favorably used for controlling plasma instabilities, such as edge-localized modes (ELMs). Here we show an optimization that tailors the EFs to maintain an edge 3D response for ELM control with a minimized core 3D response to avoid plasma disruption and unnecessary confinement degradation. We design and demonstrate such an edge-localized 3D response in the KSTAR facility, benefiting from its unique flexibility to change many degrees of freedom in the 3D coil space for the various fusion plasma regimes. This favorable control of the tokamak EF represents a notable advance for designing intrinsically 3D tokamaks to optimize stability and confinement for next-step fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Impact of toroidal rotation on the resistive ballooning modes in ASDEX Upgrade tokamak

Abstract In this work, we investigate the behavior of instabilities appearing between type-I edge localized modes (ELMs), with increasing neutral beam injection (NBI) power concomitant increase in toroidal rotation, and compare it to the modeling result of the linear magneto-hydrodynamic (MHD) code CASTOR3D. An injection of one NBI beam, increasing toroidal rotation, results in the mode slowing down from 12 kHz to 7 kHz, and its associated radial displacement decreases from 5 mm to 3.5 mm. In addition, modes shift radially outwards towards higher q , decreasing their poloidal mode numbers. The mode velocity is measured to be close to the E × B velocity with significant uncertainties. Through a set of CASTOR3D simulations with varying profiles, resistivity has been identified as the primary contributor to the growth rates. Only a small stabilizing effect due to toroidal rotation has been observed. While experimental results show a decrease of mode frequency with rotation, the opposite trend is observed in modeling. Reasons for discrepancies between modeling and experiment are discussed. Nevertheless, a main contributor to the mode frequency has been identified to be rotation velocity. CASTOR3D classifies modes as resistive ballooning modes as they do not appear unstable in ideal MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Drift-kinetic effects of tungsten on plasma response to RMP in ITER

Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.

ITER↗
DOE OSTI2020

Modeling plasma toroidal flow profile control via NTV torque with n = 2 3D fields in MAST-U

Toroidal modeling utilizing the MARS-Q code shows that large (compared to electromagnetic torque) neoclassical toroidal viscous (NTV) torque can be achieved in the plasma core region in reference MAST(-U) L-mode plasmas, by applying the n=2 (n is the toroidal mode number) 3-D fields generated by the magnetic coils used for controlling the edge localized mode and/or by the error field correction coils. Large NTV torque, occurring at relatively slow plasma flow, in combination with strong variation of the torque amplitude versus the coil phasing of the 3-D coils, offers a tool to control the plasma toroidal rotation profile in spherical tokamaks. To effectively employ the toroidal NTV torque, the potentially fast initial flow needs to first be damped by 3-D fields. In conclusion, MARS-Q quasi-linear initial value simulations demonstrate that this is achievable in MAST-U, within the designed 3-D coil current capability for the device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Disruptive neoclassical tearing mode seeding in DIII-D with implications for ITER

New studies identify the critical parameters and physics governing disruptive neoclassical tearing mode (NTM) onset. An m/n = 2/1 mode in DIII-D that begins to grow robustly after a seeding event (edge localized mode ELM or sawtooth precursor and crash) causes the mode rotation to drop close to the plasma's Er = 0 rest frame; this condition opens the stabilizing ion-polarization current 'gate' and destabilizes an otherwise marginally stable NTM. Our new experimental and theoretical insights and novel toroidal theory-based modeling are benchmarked and scalable to ITER and other future experiments. Here, the nominal ITER rotation at q = 2 is found to be stabilizing ('gate closed') except for MHD-induced transients that could 'open the gate'. Extrapolating from the DIII-D ITER baseline scenario (IBS) discharges, MHD transients are much more likely to destabilize problematic robustly growing 2/1 NTMs in ITER; this makes predictions of seeding and control of both ELMs and sawteeth imperative for more than just minimizing divertor pulsed-heat loading.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2008

System Estimates Radius of Curvature of a Segmented Mirror

A system that estimates the global radius of curvature (GRoC) of a segmented telescope mirror has been developed for use as one of the subsystems of a larger system that exerts precise control over the displacements of the mirror segments. This GRoC-estimating system, when integrated into the overall control system along with a mirror-segment- actuation subsystem and edge sensors (sensors that measure displacements at selected points on the edges of the segments), makes it possible to control the GROC mirror-deformation mode, to which mode contemporary edge sensors are insufficiently sensitive. This system thus makes it possible to control the GRoC of the mirror with sufficient precision to obtain the best possible image quality and/or to impose a required wavefront correction on incoming or outgoing light. In its mathematical aspect, the system utilizes all the information available from the edge-sensor subsystem in a unique manner that yields estimates of all the states of the segmented mirror. The system does this by exploiting a special set of mirror boundary conditions and mirror influence functions in such a way as to sense displacements in degrees of freedom that would otherwise be unobservable by means of an edge-sensor subsystem, all without need to augment the edge-sensor system with additional metrological hardware. Moreover, the accuracy of the estimates increases with the number of mirror segments.

Rakoczy, John↗
DOE OSTI2023

Impact of Soret effect on hydrogen and helium retention in PFC tungsten under ELM-like conditions

Abstract In our previous work, we have demonstrated using nonequilibrium molecular-dynamics simulations that the fluxes of helium and self-interstitial atoms in the presence of a thermal gradient in tungsten are directed opposite to the heat flux, indicating that species transport is governed by a Soret effect, namely, thermal-gradient-driven diffusion, characterized by a negative heat of transport that drives species transport uphill, i.e. from the cooler to the hot regions of the tungsten sample. In this work, the findings of our thermal and species transport analysis have been implemented in our cluster-dynamics code, Xolotl, which has been used to compute temperature and species profiles over spatiotemporal scales representative of plasma-facing component (PFC) tungsten under typical reactor operating conditions, including extreme heat loads at the plasma-facing surface characteristic of plasma instabilities that induce edge localized modes (ELMs). We demonstrate that the steady-state species profiles, when properly accounting for the Soret effect, vary significantly from those where temperature-gradient-driven transport is not accounted for and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten. Although our cluster-dynamics simulations do not yet include self-clustering of helium or hydrogen blister formation, our simulation results show that the Soret effect substantially reduces helium and hydrogenic species retention inside PFC tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effects beyond ideal MHD on stability of wide and enhanced pedestal regimes in NSTX

Stability of edge-localized modes (ELMs) in spherical tokamaks is explored using the extended MHD model. Linear NIMROD simulations have been performed for three NSTX discharges 132543, 132588, and 141133, to investigate the role of resistivity, diffusivity, and shear flows on the onset of ELMs. The first discharge represents the wide pedestal regime and the later two discharges represent ELM-free enhanced pedestal H-mode. We first present the effect of toroidal rotation shear and find a flow shear destabilizing effect in these NSTX discharges. Simulations are also extended to include the two-fluid and ion gyroviscosity effects. Simulations show that the flow shear can shift the mode spectrum and alter the critical condition of ELM onset. We also uncover that ELM onset prediction in spherical tokamaks requires effects beyond MHD, in particular gyroviscosity and diamagnetic terms could stabilize Peeling-Ballooning modes consistent with the experimental observation of ELM-free regimes in NSTX. The findings give new insight into the nature of the interplay between resistivity, flows, and diamagnetic stabilization in ELM suppression and have potential applications to ELM control schemes in NSTX-U and next-generation spherical tokamaks. This study identifies the essential physical effects that must be included in future predictive and validation simulations.

NSTX↗
DOE OSTI2023

Pedestal particle balance studies in JET-ILW H-mode plasmas

JET-ILW type I ELMy H-modes at 2.5 MA/2.8 T with constant NBI heating (23 MW) and gas fuelling rate were performed, utilising edge localised mode (ELM) pacing by vertical kicks and plasma shaping (triangularity, δ) as tools to disentangle the effects of ELMs, inter-ELM transport and edge stability on the pedestal particle balance. In agreement with previous studies, the pedestal confinement improves with increasing δ, mostly due to a significant increase in pedestal density while the ELM frequency ($f_{\mathrm{ELM}}$) is decreased. Improved pedestal confinement with increasing δ was observed even when the pedestal MHD stability was degraded artificially by vertical kicks, implying that increased triangularity may favourably affect the inter-ELM pedestal recovery. The workflow developed to quantify the pedestal particle balance uses high time-resolution profile reflectometry to characterise the inter-ELM evolution of the plasma particle content ($\mathrm{d} N/\mathrm{d}t$), the NEO drift-kinetic solver to evaluate the neoclassical fluxes and interpretative EDGE2D-EIRENE simulations to estimate the edge particle source. The edge particle source is then constrained by deuterium Balmer-α line intensity measurements in the main chamber, which are, however, strongly affected by reflections from the metal walls. The reflections are accounted for by the CHERAB code taking the divertor emission (the brightest light source in the torus) distribution from imaging spectroscopy measurements as input. Our analysis shows that in the second half of the ELM cycle, the volume-integrated particle source is larger than $\mathrm{d} N/\mathrm{d}t$, indicating that transport plays a key role in the inter-ELM pedestal recovery.

fuelling↗
DOE OSTI2024

HL-2A's ELM cycle simulations by integrating BOUT++'s drift MHD and transport code

A new integrating model has been developed to couple tokamak edge multiscale magnetohydrodynamic (MHD) events and transport simulations, such as edge-localized mode (ELM) cycles. As a proof of principle, we first start from a set of three-field two-fluid model equations, which includes the pressure, current, and vorticity. Here, the equations are separated into the slowly evolving part of the axisymmetric component by taking a time average of the axisymmetric component. The time-averaged fluxes, which are quadratic in fluctuating quantities, act as driven terms for the time-averaged axisymmetric quantities that determine the plasma transport, and therefore the large-scale evolution of the plasma profiles. Then the HL-2A's ELM cycles are simulated using the model. Good agreements of ELM size and pedestal recovery time have been achieved for the solutions obtained from the coupled simulation compared with experiment. For one ELM cycle simulation, the coupled code can achieve a speedup of a factor of up to 30 over standalone code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effect of detachment on Magnum-PSI ELM-like pulses: direct observations and qualitative results

Abstract Conditions similar to those at the end of the divertor leg in a tokamak were replicated in the linear plasma machine Magnum-PSI. The neutral pressure in the target chamber is then increased to cause the target to transition from an attached to a detached state. Superimposed to this steady state regime, edge localised mode (ELM)-like pulses are reproduced, resulting in a sudden increase in plasma temperature and density, such that the heat flux increases transiently by half an order of magnitude. Visible light emission, target thermography, and Thomson scattering are used to demonstrate that the higher the neutral pressure the more energy is removed from the ELM-like pulse in the volume. If the neutral pressure is sufficiently high, the ELM-like pulse can be prevented from affecting the target and the plasma energy is fully dissipated in the volume instead (ID 4 in table 1). The visible light images allow the division of the ELM-plasma interaction process of ELM energy dissipation into 3 ‘stages’ ranging from no dissipation to full dissipation (the target plasma is detached). In the second publication related to this study, spectroscopic data is analysed with a Bayesian approach, to acquire insights into the significance of molecular processes in dissipating the plasma energy and particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Design and experimental demonstration of feedback adaptive RMP ELM controller toward complete long pulse ELM suppression on KSTAR

Operation of a fusion power plant requires robust edge localized mode (ELM) suppression simultaneously with high plasma performance. In this study, we describe a novel feedback adaptive resonant magnetic perturbation (RMP) ELM controller designed to address this problem by achieving optimized ELM suppression through the advanced application of 3D RMPs. From real-time $D_α$ data, the controller is achievable robust ELM suppression while simultaneously minimizing the applied RMP in order to enhance plasma performance. In real time, the instantaneous ELM-frequency is analyzed with an adaptive feedback algorithm to determine amplitudes and phases of RMP coil currents that will maximize plasma performance while maintaining ELM suppression. When applied through the KSTAR plasma control system in several experiments using $textiti{n}$ = 1 RMPs, robust ELM suppression is achieved and sustained in feedback while reducing the RMP strength to ~65% of its initial value. Minimization of the RMP strength in this manner not only allows for operation of longer discharges due to a decrease in flux consumption but also allows for a strong recovery of up to ~60% of $β_N$ throughout the ELM-free period.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Quasi-linear toroidal simulations of resonant magnetic perturbations in eight ITER H-mode scenarios

Abstract Both linear and quasi-linear aspects of the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for various H-mode scenarios in ITER, covering the pre-fusion power operation and the fusion power operation phases. Linear response computations for eight ITER scenarios, with varying plasma current and toroidal magnetic field, reveal that the best coil current phasing for controlling the type-I edge localized modes (ELMs) scales roughly linearly with the edge safety factor. The coil phasing is defined as the relative toroidal phase of the coil currents between different rows, for a given toroidal harmonic. Quasi-linear initial value simulation, which is the focus of the present study, shows that application of the n = 3 ( n is the toroidal mode number) RMP field has a minimum side effect on the plasma core momentum confinement but potentially a large effect on the global particle transport. Generally, the RMP field with the best (worst) coil phasing for ELM control produces the strongest (weakest) effect on the plasma edge flow and the overall density. This robustly holds for all eight ITER scenarios. Consequently, in order to minimize the RMP induced side effects while achieving ELM control (suppression) in ITER, a compromise is necessary in choosing the coil current configuration.

Physics↗
DOE OSTI2022

Overview of recent progress in 3D field physics in KSTAR

Various 3D field physics challenges of magnetically confined plasmas arise when the driving source comes from either externally applied non-axisymmetric 3D magnetic perturbations or plasma instabilities inside the plasma. Recently, several key outstanding topics of 3D field physics have been extensively studied in the Korean Superconducting Tokamak Advanced Research (KSTAR), such as edge-localized-mode (ELM) control by resonant magnetic perturbation (RMP), error field (EF) control, 3D field effects on rotation and transport, and RMP-induced alteration of divertor heat flux and detachment. KSTAR has a few physically unique features (i.e., high rotation and long-pulse plasmas with a low intrinsic EF) and machine/diagnostic capabilities (i.e., 3-row in-vessel control coil and state-of-the-art 2D/3D imaging diagnostics), which have been taken advantage of until now to address critical 3D field physics issues relevant to ITER and K-DEMO. Among many remarkable achievements are the robust access to and control of n = 1 RMP ELM suppression, along with a development of its physics basis tools, parameter expansion, optimization, and long-pulse control techniques. Nonetheless, a series of unresolved 3D physics themes, as well as limited coverage of 3D field operating regimes, have also been identified as future works for the 3D field research in KSTAR. In this paper, we provide an overview about the recent progress of KSTAR 3D field physics and present future plans of KSTAR 3D research toward a future fusion reactor.

3D field physics↗