Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electrostatic”

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Gyrokinetic simulations of core turbulence and thermal transport in the high-β P discharge on EAST

The properties of core turbulence and thermal transport are investigated for EAST high-β P (β P ~ 3.1) plasmas with dominant electron heating (Te/Ti > 1) via gyrokinetic simulation with the NLT code. Linear simulations identify that the electrostatic η e -driven trapped electron mode (η e -TEM) dominates in the core region (ρ< 0.7) and the ion temperature gradient (ITG) mode dominates in the out region (ρ≥ 0.7), consistent with the linear threshold analysis of micro-instabilities. Sensitivity analysis shows that the normalized electron density gradient (R/L ne ) and ITG (R/L Ti ) are two effective parameters to stabilize TEM instability. Nonlinear simulations are also carried out and compared with the experimental results, which show that the electron thermal internal transport barrier (ITB) in EAST high-β P plasma is determined by the TEM-induced turbulence. A higher zonal flow shearing rate is observed in the ITB region (0.2 < ρ < 0.34), which can regulate energy transport induced by TEM turbulence and facilitate the formation of e-ITB. A plausible positive feedback mechanism to mitigate turbulence transport and improve energy confinement via enhanced ion heating is proposed for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Linear instabilities in the hot-ion regime in a high-field spherical tokamak

Abstract The present operating high-field compact spherical tokamak ST40 is an important step toward an ST-based fusion reactor (Gryaznevich and Asunta 2017 Fusion Eng. Des. 123 177; McNamara 2023 Nucl. Fusion 63 054002). Temperature and density profiles and their uncertainties in recent hot ion ST40 plasmas with central ion temperatures exceeding 8.6 keV, have been inferred using an integrated analysis of several diagnostics including line-of-sight integration and volume average measurements, as well as limited profile information from a charge-exchange-recombination spectrometer. A linear gyrokinetic stability analysis has been carried out to identify the most unstable micro-instabilities in these hot ion plasmas. In one of these plasmas, it is found that linear growth rates of both ion- and electron-scale ( k θ ρ s ⩾ 0.2 where k θ is the poloidal wavenumber and ρ s is the ion gyro-radius with sound speed) modes decrease from the edge toward the core of the plasma (i.e. from ρ = 0.8 to 0.3 where ρ is the square root of the normalized magnetic flux), while the change of linear growth rates at k θ ρ s < 0.2 is non-monotonic. In particular, at ρ = 0.3 no unstable mode was found at k θ ρ s > 5 , and about an order of magnitude reduction in the maximum linear growth rate in the ion-scale wavenumber range of 0.2 ⩽ k θ ρ s < 1 is seen at ρ = 0.3 compared with the ρ = 0.8 location. It is found through parametric scans that while the ion temperature gradient (ITG) mode/trapped electron mode (TEM) and kinetic ballooning mode (KBM)/kinetic shear Alfvén (KSA) mode are more important in the plasma core, ubiquitous mode (UM) is found to be the important ion-scale instability at ρ ⩾ 0.6 . At ρ = 0.8, two instabilities are found at the electron scale, with one being the typical electrostatic electron temperature gradient (ETG) mode at ρ e scale and beyond and the other being an unidentified low-real-frequency instability (at the intermediate scale between electron and ion gyroradii) which can propagate in the ion direction and is driven by ETG. The existence of UM at the ion scale and the low-real-frequency instability and ETG mode at the electron scale at ρ = 0.8 is supported by an extensive linear gyrokinetic stability analysis of another ST40 hot ion plasma in a wider parametric range, taking into account the uncertainties in the experimental profiles.

ubiquitous mode↗
DOE OSTI2024

Effects of hydrogen isotope species on ITG microturbulence in LHD

Abstract The linear and nonlinear effects of hydrogen isotope species on ion temperature gradient (ITG) instability in the Large Helical Device (LHD) stellarator are studied using radially global gyrokinetic simulation. We found that the coupling range of linear toroidal harmonics depends on the ion mass of the hydrogen isotope. The growth rates of ITG mode are almost the same for H, D, and T plasmas, indicating a gyro-Bohm scaling of ion-mass dependence. The nonlinear electrostatic simulations show that the zonal flow breaks the radially elongated eigenmode structures and reduces the size of the turbulence eddies, which suppresses the turbulence and the ion heat transport in the LHD. The turbulence amplitude without the zonal flow is almost the same for H, D, and T plasmas, while it decreases with increasing ion mass of the hydrogen isotope when the zonal flow is present. The reduction of the turbulent transport with larger ion mass is mostly due to the enhancement of zonal flows by larger ion mass. The ion heat conductivity deviates from the gyro-Bohm scaling for both cases with and without the zonal flow.

Physics↗
DOE OSTI2024

Sheath constraints on turbulent magnetised plasmas

Abstract A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale λ D . In magnetised plasmas where the magnetic field is inclined at an oblique angle α with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale ρ S cos ⁡ α , caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness l ms ∼ max ( λ D , ρ S cos ⁡ α ) is smaller than the characteristic length scale L of spatial fluctuations in the bulk plasma, such that the limit l ms / L → 0 is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle α ∼ δ ≡ ρ S / L ≪ 1 , electric fields tangential to the target transport ions towards the target via E × B drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm–Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.

kinetic Bohm–Chodura condition↗
DOE OSTI2025

Bifurcated merging operations of two spherical tokamak plasmas for reconnection heating and magnetic helicity injection

In tokamak plasma merging experiments, high-power ion heating is achieved only when the current sheet thickness is compressed to ion gyroradius ρ i . The magnetic reconnection converts 40%–50% of the reconnecting magnetic field energy (B rec 2 /2μ 0 ) into the ion heating energy (W i ) with the scaling law of W i proportional to B rec 2 . If the compressed current sheet thickness is larger than ρ i the magnetic reconnection converts only 5%–10% of B rec 2 /2μ 0 into W i . Finally, the high-power ion heating can be understood from the experimental finding that the poloidal electric field due to electrostatic quadrupolar potential in the downstream region increases linearly with both B rec and the guide field B t .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Role of perturbed parallel magnetic field effects in predicting turbulent transport in NSTX

This study presents analysis of gyrokinetic simulations on the National Spherical Torus Experiment (NSTX) to investigate the effects of electromagnetic fields on plasma turbulence and transport. The simulations, performed with varying levels of fidelity using the gyrokinetic CGYRO code, include electrostatic (ES), single-field electromagnetic (EM1), and two-field electromagnetic (EM2) models. A detailed comparison across the simulation database reveals that electromagnetic effects increase both predicted growth rates and quasilinear fluxes, with EM2 simulations producing stronger turbulence than ES and EM1 cases. Quasilinear modeling using QLGYRO demonstrates that while the perturbed parallel magnetic field (δB ∥ ) does not drastically affect the total flux at experimental gradients, it leads to a shift in the dominant instability, altering mode structures from microtearing to kinetic ballooning modes (KBMs). The proximity of the plasma profiles to the KBM threshold is explored, with the experimental conditions being near the onset of KBM-driven transport. The KBM, with its large growth rates, is identified as a potential driver of electron temperature flattening, as it can rapidly transport heat across flux surfaces. Performing stability analysis shows core-localized unstable a low- mode that could contribute to the flattening at the early times of the discharge. TGYRO predictive modeling, incorporating both TGLF and QLGYRO, indicates that the inclusion of δB ∥ significantly improves the accuracy of temperature profile predictions in NSTX high-beta plasmas, although challenges remain in modeling the sharp flux discontinuities caused by KBM-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗
DOE OSTI2025

Investigation of radiated-power for low aspect ratio fusion plasmas

Investigating radiated power in fusion plasmas is of utmost importance to understand the effect of undesired impurities, such as metals, in present devices, or also desired impurities, such as noble gasses, to purposefully radiate a large fraction of power in future devices. These studies are especially important for high Z impurities which will play a crucial role in future generation fusion pilot plants. In this work, we have developed a power radiation analysis module, which is used to investigate 2D distributions of impurity densities and radiated power asymmetries caused due to both plasma rotation and the cooling rate dependence on temperature profile for the cases of experimental NSTX plasmas and designed scenarios for the spherical tokamak advanced reactor (STAR), both being low aspect ratio tokamaks. Two different atomic databases have been tested during this work to study their impact on the radiated power distribution, especially for high Z impurities. Also in PRAM, self-consistent calculations of two-dimensional electron, main ion, and impurity ion densities are carried out using one-dimensional input density, temperature, and rotation profiles. In the case of NSTX, discharges with high rotation of ∼ 170 km s, measured with charge exchange recombination spectroscopy, have been investigated. Rotation-induced charge separation, leading to an electrostatic potential, is calculated iteratively to a self-consistent solution while testing high Z impurities to observe any 2D asymmetry in the core radiated power due to centrifugal forces. The STAR design, being much larger (R = 4 m), is projected to have a much lower rotation, and is shown to have low rotation-induced asymmetries, on the order of ten percent or less, between the low field and high field sides. However, another effect not due to rotation but to the dependence of impurity cooling rates on temperature can lead to radiation peaking off-axis, near the plasma edge. This effect is noticeable for argon in NSTX, for example, but can also be enhanced for certain impurities at much higher temperatures projected for STAR (T e0 ~ 32 keV), for example for undesired tungsten or possibly desired xenon.

NSTX↗
DOE OSTI2020

Experimental studies of the difference between plasma potentials measured by Langmuir probes and emissive probes in presheaths

The plasma potential measured by cylindrical and planar Langmuir probes has been shown to differ from the plasma potential measured by emissive probes in the neighborhood of the presheath near a negatively biased electrode immersed in a weakly collisional low temperature argon plasma. There are two principal results demonstrated in this paper. First, while it is well known that Langmuir probes cannot reliably measure plasma potentials inside of sheaths, results presented here demonstrate that the problem persists in presheaths, the quasineutral plasma bordering sheaths. It is known that emissive probes analyzed in the limit of zero emission accurately measure the plasma potential in the sheath. It is now clear that they are the only known electrostatic probe technique able to measure the plasma potential accurately throughout the presheath. Second, it is shown that the difference between potential measurements made by Langmuir probes and emissive probes in the body of the plasma, farther than a presheath distance from the boundary, is not proportional to T e , as has been previously claimed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electron dynamics in radio frequency magnetron sputtering argon discharges with a dielectric target

Abstract We demonstrate a self-consistent and complete description of electron dynamics in a typical electropositive radio frequency magnetron sputtering (RFMS) argon discharge with a dielectric target. The electron dynamics, including the electron power absorption dynamics in one radio frequency (RF) period, is studied via a fully kinetic 2d3v particle-in-cell/Monte Carlo collision (PIC/MCC) electrostatic simulation. The interplay between the fundamental plasma parameters is analyzed through their spatiotemporal dynamics. Due to the influence of magnetic trap on the electron transport, a spatially dependent charging that perturbs the electric potential is observed on the dielectric target surface, resulting in a spatially dependent ion energy distribution along the target surface. The E × B drift-to-discharge current ratio is in approximate agreement with Bohm diffusion. The electron power absorption can be primarily decoupled into the positive Ohmic power absorption in the bulk plasma region and the negative pressure-induced power absorption near the target surface. Ohmic power absorption is the dominant electron power absorption mechanism, mostly contributed by the azimuthal electron current. The power absorption due to electron inertial effects is negligible on time-average. Both the maximum power absorption and dissipation of electrons appear in the bulk plasma region during the second half of the RF period, implying a strong electron trapping in magnetron discharges. The contribution of secondary electrons is negligible under typical RFMS discharge conditions.

Physics↗
DOE OSTI2023

The transport dynamics of tens of micrometer-sized water droplets in RF atmospheric pressure glow discharges

Abstract Charging of particles having diameters of tens of microns has been extensively studied at atmospheric pressure in the context of, for example, electrostatic precipitators where the focus was on unipolar charging. The ambipolar charging of particles in atmospheric pressure plasmas, and of droplets in particular, has received less attention. The plasma activation of droplets is of interest for water purification, fertilizer production and materials synthesis, all of which depend on the transport of the droplets through the plasma, which in turn depends on their charging. In this paper, we report on the transport dynamics of water droplets, tens of microns in diameter, carried by the gas flow through an atmospheric pressure radiofrequency glow discharge sustained in helium. The droplets pass through the plasma with minimal evaporation and without reaching the Rayleigh limit. The droplet trajectory in the presence and absence of the plasma provides insights on the forces acting on the droplet. The measurements were analyzed using results from a three-dimensional fluid model and a two-dimensional plasma hydrodynamics model. We found that the transport dynamics as the droplet enters and leaves the plasma are due to differential charging of the droplet in the plasma gradients of the bounding sheaths to the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Study of synchronous RF pulsing in dual frequency capacitively coupled plasma

Low-pressure multi-frequency capacitively coupled plasmas (CCPs) are used for numerous etch and deposition applications in the semiconductor industry. Pulsing of the radio-frequency (RF) sources enables control of neutral and charged species in the plasma on a millisecond timescale. The synchronous (i.e. simultaneous, in-phase) pulsing of both power sources in a dual frequency CCP is examined in this article. Due to the low gas pressure, modeling has been done using the electrostatic particle-in-cell/Monte Carlo collision method. Here, the objective of this work is to investigate the sensitivity of the plasma properties to small changes in timing during synchronous pulsing of the two RF sources. It is demonstrated that small deviations in the on and off times of the two RF sources can lead to major changes in the plasma characteristics. This high sensitivity is of concern for process repeatability but can be utilized to enable better control of the dynamics of plasma-surface interaction. In the simulations, the pulsing parameters (on and off times and ramp rates) are varied and the temporal evolution of plasma characteristics such as electron density (n e ), species current at the electrode, and electron temperature are examined. It is demonstrated that if the low-frequency (LF) source is turned off a few μs before (or after) the high-frequency source, n e during the off-state is significantly higher (or lower) due to the frequency coupling effect. Similarly, turning on the LF source with a small delay results in a sharp increase in the plasma density when the RF sources are turned on.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Understanding atom probe’s analytical performance for iron oxides using correlation histograms and ab initio calculations

Abstract Field evaporation from ionic or covalently bonded materials often leads to the emission of molecular ions. The metastability of these molecular ions, particularly under the influence of the intense electrostatic field (10 10 Vm −1 ), makes them prone to dissociation with or without an exchange of energy amongst them. These processes can affect the analytical performance of atom probe tomography (APT). For instance, neutral molecules formed through dissociation may not be detected at all or with a time of flight no longer related to their mass, causing their loss from the analysis. Here, we evaluated the changes in the measured composition of FeO, Fe 2 O 3 and Fe 3 O 4 across a wide range of analysis conditions. Possible dissociation reactions are predicted by density-functional theory calculations considering the spin states of the molecules. The energetically favoured reactions are traced on to the multi-hit ion correlation histograms, to confirm their existence within experiments, using an automated Python-based routine. The detected reactions are carefully analyzed to reflect upon the influence of these neutrals from dissociation reactions on the performance of APT for analysing iron oxides.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Impact of divertor material on neutral recycling and discharge fueling in DIII-D

We report that experiments with the lower divertor of DIII-D during the Metal Rings Campaign (MRC) show that the fraction F of atomic D in the total recycling flux is material-dependent and varies through the ELM cycle, which may affect divertor fueling. Between ELMs, F C ~ 10% and FW ~ 40%, consistent with expectations if all atomic recycling is due to reflections. During ELMs, F C increases to 50% and F W to 60%. In contrast, the total D recycling coefficient including atoms and molecules R stays close to unity near the strike point where the surface is saturated with D. During ELMs, R can deviate from unity, increasing during high energy ELM-ion deposition (net D release) and decreasing at the end of the ELM which leads to ability of the target to trap the ELM-deposited D. The increase of R > 1 in response to an increase in ion impact energy Ei has been studied with small divertor target samples using Divertor Materials Evaluation System (DiMES). An electrostatic bias was applied to DiMES to change E i by 90 eV. On all studied materials including C, Mo, uncoated and W-coated TZM (>99% Mo, Ti, and Zr alloy), W, and W fuzz, an increase of E i transiently increased the D yield (and R) by ~10%. On C there was also an increase in the molecular D 2 yield, probably due to ion-induced D 2 desorption. Despite the measured increase in F on W compared to C, attached H-mode shots with OSP on W during MRC did not demonstrate a higher pedestal density. About 8% increase in the edge density could be seen only in attached L-mode scenarios. The difference can be explained by higher D trapping in the divertor and lower divertor fueling efficiency in H- versus L-mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Self-consistent interaction of linear gravitational and electromagnetic waves in non-magnetized plasma

Here, this paper explores the hybridization of linear metric perturbations with linear electromagnetic (EM) perturbations in non-magnetized plasma for a general background metric. The local wave properties are derived from first principles for inhomogeneous plasma, without assuming any symmetries of the background metric. First, we derive the effective (“oscillation-center”) Hamiltonian that governs the average dynamics of plasma particles in a prescribed quasimonochromatic wave that involves metric perturbations and EM fields simultaneously. Then, using this Hamiltonian, we derive the backreaction of plasma particles on the wave itself and obtain gauge-invariant equations that describe the resulting self-consistent gravito-electromagnetic (GEM) waves in a plasma. The transverse tensor modes of gravitational waves are found to have no interaction with the plasma and the EM modes in the geometrical-optics limit. However, for $longitudinal$ GEM modes with large values of the refraction index, the interplay between gravitational and EM interactions in plasma can have a strong effect. In particular, the dispersion relation of the Jeans mode is significantly affected by electrostatic interactions. As a spin-off, our calculation also provides an alternative resolution of the so-called Jeans swindle.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Testing predictions of electron scale turbulent pedestal transport in two DIII-D ELMy H-modes

In this study, we present linear and nonlinear gyrokinetic analyses in the pedestal region of two DIII-D ELMy H-mode discharges using the CGYRO code. The otherwise matched discharges employ different divertor configurations to investigate the impact of varying recycling and particle source on pedestal profiles. Linear gyrokinetic simulations find electrostatic ion-scale instabilities (ion temperature gradient and trapped electron modes, ITG-TEM) are present just inside the top of the pedestal with growth rates that are enhanced significantly by parallel velocity shear. In the sharp gradient region, E x B shearing rates are comparable or larger than ion scale growth rates, suggesting the suppression of ITG-TEM modes in this region. Instead, the electron temperature profiles are found to be correlated with and just above the electron temperature gradient (ETG) instability thresholds. Using gradients varied within experimental uncertainties, nonlinear electron-scale gyrokinetic simulations predict electron heat fluxes from ETG turbulence, that when added to neoclassical ion thermal transport simulated by NEO, account for 30-60% of the total experimental heat flux. In addition, the neoclassical electron particle flux is found to contribute significantly to the experimental fluxes inferred from SOLPS-ITER analysis. Additional nonlinear gyrokinetic simulations are run varying input gradients to develop a threshold-based reduced model for ETG transport, finding a relatively simple dependence on η e = L ne /L Te . Predictive transport simulations are used to validate this pedestal-specific ETG model, in conjunction with a model for neoclassical particle transport. In both discharges, the predicted electron temperatures are always overpredicted, indicative of the insufficient stiffness in the ETG pedestal model to account for all of the experimental electron thermal transport. In the case of the closed divertor discharge with lower particle source, the predicted electron density is close to the experiment, consistent with the magnitude of neoclassical particle transport in that discharge. However, the density profiles are overpredicted in the open divertor discharge (larger particle source), due to insufficient model transport. In conclusion, the implications for other mechanisms accounting for the remainder of transport in the sharp gradient region in the two discharges are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

W7-X and the sawtooth instability: towards realistic simulations of current-driven magnetic reconnection

Magnetic reconnection in W7-X is studied by means of global numerical simulations in a series of models of increasing complexity. The magnetic geometry ranges from that of a cylinder to the full three-dimensional field of W7-X, and the equations solved range from ideal magnetohydrodynamics (MHD) to gyrokinetics. We simulate plasmas from the first operation phase with electron cyclotron current drive (ECCD). These are characterized by an equilibrium magnetic field featuring an ECCD-distorted 'humped' profile of the rotational transform ι, with ι = 1 in two radial locations. Such plasmas generally show sawtooth activity, hence motivating the present study. We pay particular attention to the role of equilibrium current density gradients in the destabilization of reconnecting modes. When the equilibrium temperature and density gradients are artificially suppressed (to eliminate the pressure gradient drive), the perturbed electrostatic potential is radially localized between the locations at which ι = 1. This is shown with a purely collisionless gyrokinetic model, in cylindrical geometry. In the real toroidal geometry of W7-X, for a non-ideal MHD model including a uniform resistivity, electron inertia and (numerical) viscosity, the same qualitative behaviour is observed. In particular, even if a resonant (m, n) = (1, -1) perturbation is initialized, the most unstable mode is the (m, n) = (-4, 4), where m and n are the poloidal and toroidal mode numbers, respectively. Other modes are destabilized due to geometric coupling. The growth rate of this instability scales as η 1/3 , where η is the plasma resistivity, thus suggesting that ECCD drives ideal MHD stable W7-X plasmas towards non-ideal marginality. An ideal magnetohydrodynamic analysis confirms the result. A fluid-kinetic hybrid version of the EUTERPE code shows that gyrokinetic ions have a stabilizing effect on these modes. For W7-X relevant collisionalities, the growth rate scales linearly with the electron skin depth, d e . Implications of our results for sawtoothing W7-X operation are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗