Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MARS-F”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Chasing the multi-modal plasma response in MAST-U

Abstract Achieving edge localized modes (ELMs) suppression in spherical tokamaks by applying resonant magnetic perturbations (RMPs) has proven challenging. The poloidal spectrum of the applied RMP is a key parameter that has an impact on the capability to mitigate and eventually suppress ELMs. In this work the resistive magnetohydrodynamic code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is used to evaluate the possibility of directly measuring the plasma response in MAST-U, and particularly its variation as function of the applied poloidal spectrum, in order to guide the experimental validation of the predicted best RMP configuration for ELM suppression. Toroidal mode number n = 2 RMP is considered to minimize the presence of sidebands, and to avoid the deleterious core coupling of n = 1. Singular Value Decomposition is used to highlight linearly independent structures in the simulated magnetic 3D fields and how those structures can be measured at the wall where the magnetic sensors are located. Alternative ways to measure the multimodal plasma response and how they can be used to infer the best RMP configuration to achieve ELM suppression are also presented, including the plasma displacement and the 3D footprints at the divertor plates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

First observation of RMP ELM mitigation on MAST Upgrade

Abstract The first experimental attempts at controlling edge localised modes (ELMs) via the application of resonant magnetic perturbations on the MAST Upgrade tokamak are reported. Using the linear MHD model MARS-F, the phase shift between the upper and lower coil rows ΔΦ was optimised for toroidal mode number n = 1 and n = 2 fields, to provide forward guidance to experiments. In low β N discharges, the application of n = 1 3D fields caused the ELM frequency f E L M to increase by over a factor 20 relative to the reference, and also induced a locked mode, which did not cause a plasma termination nor an H-L back transition. However when β N was raised, this induced locked mode caused plasma termination which precluded mitigation access. Initially, applying a numerically optimised n = 2 field had no effect. However applying a rigid toroidal shift to this field caused a locked mode disruption, demonstrating the presence of a substantial n = 2 error field. Coil current ramps were conducted with ΔΦ set at 6 different values, resulting in either locked mode disruptions or no effect, but mitigation with n = 2 fields was not established.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D

Post-disruption runaway electron (RE) kinetic energy K and pitch angle sin$\vartheta$ are critical parameters for determining resulting first wall material damage during wall strikes, but are very challenging to measure experimentally. During the final loss instability, confined RE K and sin$\vartheta$ are reconstructed during center-post wall strikes for both high impurity (high-Z) and low impurity (low-Z) plasmas by combining soft x-ray, hard x-ray, synchrotron emission, and total radiated power measurements. Deconfined (wall impacting) RE sin$\vartheta$ is then reconstructed for these shots by using time-decay analysis of infra-red imaging. Additionally, deconfined RE K and sin$\vartheta$ are reconstructed for a low-Z downward loss shot by analyzing resulting damage to a sacrificial graphite dome limiter. The damage analysis uses multi-step modeling simulating plasma instability, RE loss orbits, energy deposition, and finally material expansion (MARS-F, KORC, GEANT-4, and finally COMSOL). Overall, mean kinetic energies are found to be in the range $\langle$$K$$\rangle$ $≈ 3-4$ MeV for confined REs. KORC simulations indicate that the final loss instability process does not change individual RE kinetic energy K. In conclusion, confined RE pitch angles are found to be fairly low initially pre-instability, $\langle$sin$\vartheta$$\rangle$ $≈ 0.1 – 0.2$, but appear to increase roughly 2x, to $\langle$sin$\vartheta$$\rangle$ $≈ 0.3 – 0.4$ for both confined and deconfined REs during instability onset in the low-Z case; this increase is not observed in the high-Z case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Modeling fast ion losses due to tearing and internal kink perturbations in MAST-U

Fast ion (FI) loss properties in the presence of tearing mode and internal kink perturbations are numerically investigated for discharges in the MAST-U spherical tokamak, utilizing the MARS-F magnetohydrodynamic stability code and the REORBIT test particle guiding-center orbit-following module. Here, tracing about 100 000 particle markers sampled from the equilibrium distribution of the neutral-beam injection induced FIs, it is found that about 10% out of the total strike the limiting surface (including the divertor surface) in MAST-U discharge 46943, assuming a maximum perturbation of 100 G inside the plasma (corresponding to ~6 G at the Mirnov probe location at the outboard mid-plane). Detailed particle tracing, assuming a uniform initial distribution in the 2D phase space (at given radial locations), reveals that initially counter-current FIs launched near the plasma edge are subject to significant prompt losses, while almost all initially co-current ions remain well confined at the assumed perturbation level. Most lost FIs strike the lower-half of the limiting surface. Finite gyro-radius effects prevent lost ions from striking the top-outer corner of the super-X divertor chamber. A scan of the perturbation level (based on discharge 45163) reveals, not surprisingly, an approximately linear scaling of the particle loss fraction (for counter-current FIs) with respect to the perturbation amplitude.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of equilibrium pressure on plasma response to RMPs in a spherical tokamak

This study presents a comprehensive analysis of the equilibrium pressure on the plasma response to resonant magnetic perturbations (RMPs) in the spherical tokamak (ST) MAST-U, employing both single-fluid and MHD-kinetic hybrid models (implemented via the MARS-F/K codes). As a key finding, the study identifies two different pressure-driven eigenmodes, exhibiting Sturmian property, that affect the Troyon no-wall limits for the onset of the n = 1 and n = 2 ( n is the toroidal mode number) ideal external kink instabilities as well as the corresponding plasma response to the applied RMP. With increasing equilibrium pressure, the plasma response to RMPs is significantly enhanced in the ST plasma, particularly in the high-pressure regime where kinetic effects strongly stabilize the external kink instability. The Troyon no-wall limit divides the plasma response into two regions: well below the limit, the response amplitudes and trends (versus pressure) are similar between the fluid and kinetic models; as the equilibrium pressure approaches the Troyon limit, the kinetic model predicts significant amplification of the RMP field, up to 30 times for cases considered. A relatively weak dependence of the optimal coil phasing on the equilibrium pressure is computed in this ST plasma, similar to the trend obtained for the conventional aspect ratio devices. These findings underscore the importance of incorporating kinetic effects in accurate prediction of the plasma response to RMPs in high-pressure ST tokamak plasmas and provide a theoretical basis for optimizing RMP-based control of the edge-localized modes in future ST devices.

RMP↗
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↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental investigation of beam-ion losses induced by magnetic perturbations using the light ion beam probe technique in the ASDEX Upgrade tokamak

The impact of externally applied magnetic perturbations (MPs) on fast-ion losses has been investigated by means of the light ion beam probe (LIBP) technique in the ASDEX Upgrade tokamak. The LIBP technique allows to experimentally infer the fast-ion orbit displacement induced by MPs via first-orbit losses using scintillator based fast-ion loss detector (FILD) measurements. The fast-ion orbit displacement against different applied MP spectra has been studied. These shots were conducted in ELM mitigated H-mode plasmas. A rigid rotation of the MP coils was applied with a frequency of 1 Hz, with an n = 2 configuration and changing the differential phase between the upper and lower set of coils (ΔΦ ul ) on a shot-to-shot basis. Beam sources Q7 (tangential) and Q8 (radial) were used to probe different fast-ion orbits with FILD1. The measured fast-ion orbit displacement ranges from 3 to 20 mm approximately, and no qualitative difference is observed between ions from beam sources Q7 and Q8. The minimum is found for a ΔΦ ul ~ 50°, which is shifted with respect to the minimum of the plasma boundary displacement, found at ΔΦ ul ~ 0°. A first attempt to validate the orbit following code ASCOT—including the plasma response calculated with the MARS-F code—against these experimental measurements is performed. Finally, while the dependence of the first-orbit fast-ion displacement with ΔΦ ul does not match the experimental measurements, these simulations do capture other features such as the order of magnitude of the orbit displacement and the importance of the toroidal spectrum of the applied perturbation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Surrogate models for plasma displacement and current in 3D perturbed magnetohydrodynamic equilibria in tokamaks

Abstract A numerical database of over one thousand perturbed three-dimensional (3D) equilibria has been generated, constructed based on the MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) computed plasma response to the externally applied 3D field sources in multiple tokamak devices. Perturbed 3D equilibria with the n = 1–4 ( n is the toroidal mode number) toroidal periodicity are computed. Surrogate models are created for the computed perturbed 3D equilibrium utilizing model order reduction (MOR) techniques. In particular, retaining the first few eigenstates from the singular value decomposition (SVD) of the data is found to produce reasonably accurate MOR-representations for the key perturbed quantities, such as the perturbed parallel plasma current density and the plasma radial displacement. SVD also helps to reveal the core versus edge plasma response to the applied 3D field. For the database covering the conventional aspect ratio devices, about 95% of data can be represented by the truncated SVD-series with inclusion of only the first five eigenstates, achieving a relative error (RE) below 20%. The MOR-data is further utilized to train neural networks (NNs) to enable fast reconstruction of perturbed 3D equilibria, based on the two-dimensional equilibrium input and the 3D source field. The best NN-training is achieved for the MOR-data obtained with a global SVD approach, where the full set of samples used for NN training and testing are stretched and form a large matrix which is then subject to SVD. The fully connected multi-layer perceptron, with one or two hidden layers, can be trained to predict the MOR-data with less than 10% RE. As a key insight, a better strategy is to train separate NNs for the plasma response fields with different toroidal mode numbers. It is also better to apply MOR and to subsequently train NNs separately for conventional and low aspect ratio devices, due to enhanced toroidal coupling of Fourier spectra in the plasma response in the latter case.

3D equilibrium↗
DOE OSTI2022

Toroidal modeling of 3D perturbations generated by current filaments in scrape-off layer in tokamak with biased divertor targets

The divertor biasing technique offers a promising alternative to control the edge localized mode (ELM) as well as the divertor heat load in tokamaks, as compared to the resonant magnetic perturbation (RMP) generated by magnetic coils. The linear resistive MHD code MARS-F (Liu et al 2000 Phys. Plasmas 7 3681) is employed to study the plasma response to the n = 2 (n is the toroidal mode number) field perturbation in HL-2A, produced by the divertor biasing current filaments flowing in the scrape-off layer region. The response field due to biasing currents is compared with the RMP field utilized for ELM control in HL-2A. The strength of the radial resonant field perturbation, produced by the biasing currents of 100 A level, is found to be comparable with RMP generated by several kA of ELM control coil currents for the reference plasma. The plasma normal displacement near the X-point and the associated neoclassical toroidal viscosity torque are also computed to be similar between these two techniques. The modeling results thus strongly suggest that the biasing technique can be applied to control ELMs. Moreover, the biasing currents produce field perturbations, including the plasma response, that are localized more near the plasma edge (compared to the RMP counterpart), thus reducing the chance of mode locking associated with core perturbations. Particle orbit tracing also reveals that the biasing current produced magnetic perturbation tends to widen the heat deposition region and induce the strike point splitting of the ion saturation flow on the outer divertor surface, consistent with experimental observations in HL-2A. These toroidal modeling results confirm the possibility of ELM control and plasma exhaust solution by the divertor biasing technique.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effects of external kink and fishbone-like modes on energetic particle transport in tokamak plasmas

Transport and loss of beam injected energetic particles (EPs) due to three-dimensional perturbations, associated with the external kink (XK) instability and fishbone-like mode (FLM), are numerically investigated utilizing the guiding center following code ORBIT for static toroidal plasmas in HL-2A. The perturbation structure for the XK is computed by the MARS-F code and then mapped to the Boozer coordinates as defined in ORBIT. The simulation shows that the EP profile experiences a significant change in the middle of the plasma column, when the XK-induced radial magnetic field perturbation amplitude, normalized by the equilibrium field, exceeds a threshold value of about 10 -2 . The EP transport is found to be dominated by a diffusion process instead of convection. Furthermore, by scanning the perturbation frequency as a free parameter while maintaining the XK mode structure (thus mimicking the FLM as observed in DIII-D and JT-60U tokamaks), redistribution and loss of EPs are found to be substantially enhanced due to strong resonances between the FLM and EPs, when the mode frequency exceeds a threshold value of ~2 kHz for the case considered. For either XK or FLM, the response of passing EPs to the perturbation is dominant due to the assumed tangential neutral beam injection. Most lost EPs due to these instabilities are initially passing particles but are eventually lost through trapped orbits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Resonant mode effects on rotation braking induced by n = 1 resonant magnetic perturbations in the EAST tokamak

Abstract The spectrum effects on toroidal rotation braking, induced by n = 1 resonant magnetic perturbations (RMPs) in the discharges with q 95 = 4.1 and q 95 = 5.1 , are studied in the EAST tokamak. Here n is the toroidal mode number, RMP spectrum is varied by scanning δ ϕ U L , the phase difference between the upper and lower rows of RMP coils. The toroidal rotation changes periodically with the periodic δ ϕ U L scanning and such an effect is stronger in the discharge with lower q 95 = 4.1 . The spectrum dependence of the neoclassical toroidal viscosity (NTV) torque, modeled by NTVTOK based on the magnetic perturbation obtained from MARS-F calculation, agrees well with that of the experimentally observed braking torques in both discharges. The modeled NTV torque is stronger in the discharge with lower q 95 , which also agrees with the observations. The comparisons between the spectrum dependence of the NTV and magnetic perturbations show that the resonant mode of magnetic perturbations near the plasma edge mainly contribute the NTV torque. These agreements between modeling and experiments highlight the capability of NTV theory in explaining the experimental observation in the EAST tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Effect of resonant magnetic perturbations including toroidal sidebands on magnetic footprints and fast ion losses in HL-2M

Externally applied resonant magnetic perturbations (RMPs), generated by magnetic coils located outside the plasma (referred to as RMP coils), provide an effective way to control the edge localized mode (ELM) in tokamak devices. Due to the discrete nature of the toroidal distribution of these window-frame coils, toroidal sidebands always exist together with the fundamental harmonics designed for ELM control. In this work, the MARS-F code (Liu et al 2000 Phys. Plasmas7 3681) is applied to investigate the detailed features of the RMP spectra considering both the dominant harmonic (n = 2) and the associated sideband (n = 6), and the impact of the combined fields on magnetic footprints as well as on the fast ion losses for a reference double-null scenario in the HL-2M device. It is found that the sum of the n = 2 and n = 6 RMP fields splits the footprint and widens the footprint area, as compared to the single-n (n = 2) harmonic case. The resistive plasma response breaks the up–down symmetry of the footprint pattern on the outer divertor plates, which is otherwise symmetric assuming vacuum RMP fields. Considering fast ion losses, a threshold value exists for the initially launched radial position of test particles, as well as for the RMP coil current, before the loss occurs. When the threshold criterion is satisfied, the combined n = 2 and n = 6 RMP fields enhance the fast ion loss rate by , as compared to that of the n = 2 component alone. These results illustrate the important role of the sideband of RMP fields on the magnetic footprints and fast ion losses in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Real time detection of multiple stable MHD eigenmode growth rates towards kink/tearing modes avoidance in DIII-D tokamak plasmas

Abstract Real time detection of time evolving growth rates of multiple stable magnetohydrodynamic (MHD) eigenmodes has been achieved in DIII-D tokamak experiments via multi-mode three-dimensional (3D) active MHD spectroscopy. The measured evolution of the multi-modes’ growth rates is in good accordance with the variation of the plasma β N . Using experimental equilibria, resistive MARS-F simulations found the two least stable modes to have comparable growth rates to those experimentally measured. Real time and offline calculations of the modes’ growth rates show comparable results and indicate that cleaner system input and output signals will improve the accuracy of the real time stability detection. Moreover, the shortest real time updating time window of multi-mode eigenvalues can be about 2 ms in DIII-D experiments. This real time monitoring of stable, macroscopic kink and tearing modes thus provides an effective tool for avoidance of the most common causes of tokamak disruption.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

MHD stability analysis against pressure and current-driven modes in the SMall Aspect Ratio Tokamak

Linear magnetohydrodynamic (MHD) simulations for the SMall Aspect Ratio Tokamak (SMART) have been carried out for the first time, for both positive (PT) and negative triangularity (NT) shaped plasmas using the MARS-F code. The MHD stability of projected SMART plasmas against internal kinks, infernal modes and edge peeling-ballooning modes have been analyzed for a wide range of realistic equilibria. A stabilization of internal kinks and infernal modes is observed when increasing the safety factor profile and reducing plasma beta. PT shaped plasmas are more stable against both internal kinks and infernal modes than their counterpart NT shaped plasmas. Toroidal flows have little impact on the MHD stability of the internal kinks, but they have a strong stabilizing effect on infernal modes, which can be further mitigated in NT shaped plasmas. The MHD stability of peeling-ballooning modes is reduced in NT shaped plasmas, as observed in conventional tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Comparison of plasma response models for RMP effects on the divertor and scrape-off layer in KSTAR

Resonant magnetic perturbations (RMPs) are beneficial for control of edge localized modes (ELMs) in tokamaks. Nevertheless, a side effect of RMPs is the appearance of helical striations in the particle and heat loads onto divertor targets. The extent and field line connection of these striations is significantly altered by the plasma response to external perturbations. For an ELM suppressed high confinement (H-mode) plasma at KSTAR, magnetic footprints are computed by FLARE based on plasma response from GPEC, MARS-F, M3D-C1 and JOREK with substantial differences in the resulting footprints (from 2 cm to 14 cm). This is reflected in EMC3-EIRENE simulations of the resulting heat loads: it is found that either the peak value or the extent of the striations appears to be overestimated compared to infrared camera measurements. Reasonable agreement can only be achieved for the smallest footprint for lower input power and lower cross-field transport, or for higher upstream density and radiative power losses.

Frerichs, H. [Univ. of Wisconsin, Madison, WI (Uni↗
DOE OSTI2026

Toroidal phase of post-disruption runaway electron loss to wall in presence of applied 3-D fields in DIII-D

The toroidal distribution of runaway electrons (REs) striking the centre post in DIII-D post-disruption RE final loss events is simulated utilizing the MARS-F code combined with the REORBIT module, by tracing the guiding-centre drift orbits of test REs in the presence of external (applied) plus internal (intrinsic instability) 3-D fields. To better recover the experimental results, three different equilibria with different safety factor profiles are adopted. Dominant resistive magnetohydrodynamic instabilities are found to be m/n=2/1, 3/2 and 1/1 modes (where m and n are the poloidal and toroidal mode numbers). Externally applied n=1 resonant magnetic perturbations (RMPs), with both even and odd parity configuration, are considered in each equilibrium. Locking of MHD instability to the applied RMPs is simulated by minimizing the total perturbed magnetic energy in the plasma. Here, the simulated toroidal locations, where most REs strike in the presence of an edge m/n= 2/1 instability and the n=1 RMPs in even parity, are found to agree well with experiments. Modeling also captures the measured n = 2 RE loss pattern on the wall for experiments where an edge m/n=3/2 instability is dominant instead of the m/n=2/1 mode. The key overall findings include (a) the measured global RE toroidal impact pattern in experiments can be explained by a locked resistive kink instability, (b) the toroidal phase of the locked mode is well predicted by minimizing the total perturbed magnetic energy, and (c) this edge locked mode determines the toroidal phase of the peak RE wall impact.

Bai, Xue [University of California San Diego, La J↗
DOE OSTI2025

A 3D helical filament surrogate model for 3D tokamak equilibria

A novel approach for efficient representation of three-dimensional (3D) tokamak equilibria is investigated, where a set of helical current filaments occupying the plasma region are employed to resolve deviations from the two-dimensional (2D) axi-symmetric state. A discrete set of 3D filaments, located at rational surfaces for a given toroidal mode number n and following the 2D equilibrium field lines (thus forming closed current loops), are found to provide a surrogate model of 3D equilibria with reasonable accuracy. Specifically, application of the filament model to 3D perturbed equilibria, due to the resonant magnetic perturbation (RMP) in DIII-D and MAST-U discharges, reveals that (1) a single helical filament per rational surface is sufficient; (2) 21 such helical filaments are capable of representing the n = 2 3D response field in MAST-U with less than 10% relative error as compared to that computed by a full magnetohydrodynamic code; (3) optimizing currents (both amplitude and phase) flowing in 3D filaments with fixed geometry, the highest accuracy fitting is found to depend on the characteristics of the 3D equilibria such as the coil current phasing of the RMP coils in our case studies. Here, whis filament approach is also applicable for generating surrogate models of other type of 3D tokamak equilibria, including those during the initial phase of the plasma disruption.

MARS-F↗