Search NASA⌕ Search

SEARCH · Search NASA

Results for “Boundary plasma”

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 OSTI2021

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radio frequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto (J. Plasma Phys., vol. 86, 815860302, 2020). The behavior is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston (Phys. Fluids, vol. 14, 1971, pp. 2719-2726) and Auerbach (Phys. Fluids, vol. 20, 1977, pp. 1836-1844). For a wave of parallel wave number k|| interacting with weakly collisional plasma species of collision frequency ν and thermal speed vth , the effective collision frequency becomes of order ν(k_||v_th /ν)^2/3>> ν . The narrow boundary layers that arise because of the diffusive nature of the collisions allows a physically meaningful wave-particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behavior of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave-particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Identification of a window for quiescent H-mode operation in MHD stability diagram of DIII-D plasmas

A window for quiescent H-mode (QH-mode) operation in DIII-D plasmas in the diagram of pedestal magnetohydrodynamic (MHD) stability was identified through linear MHD stability analysis, accounting for the effects of plasma rotation and ion diamagnetic drift. The operation window lies between the stability boundaries of the kink/peeling mode (K/PM) identified with and without the inclusion of plasma rotation effects. In this region, the mode remains unstable unless rotation effects are considered alongside the ion diamagnetic drift, which is consistently included in the analysis. The stabilization of the MHD mode, facilitated by the coupled effects of plasma rotation and ion diamagnetic drift, plays a crucial role in widening the window, enabling the attainment of the QH-mode state. Furthermore, the suppression of edge localized modes (ELM) can be achieved by controlling the pedestal structure to ensure the plasma state remains within the operation window. The location of the operation window in the stability diagram depends on the K/PM stability properties. Therefore, optimizing conditions for QH-mode requires adjustments based on changes in stability characteristics. A pressure pedestal and its associated bootstrap current density near the last closed flux surface are advantageous for situating the plasma state within the window. However, excessive current density can trigger ELMs. This trend was confirmed through comparisons of MHD stability diagrams between QH-mode and ELMy H-mode plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Collisionless plasma transport mechanisms in stochastic open magnetic field lines in tokamaks

The collisionless plasma transport in stochastic magnetic fields connecting to wall boundary has been studied to understand the underlying mechanisms of thermal quench physics in tokamak disruption. In this study, we present a comprehensive understanding of the plasma dynamics in three-dimensional stochastic open magnetic field lines, taking into account the consistent coupling of electron and ion dynamics through the ambipolar electric field. The open stochastic field lines act as a 3D magnetic mirror consisting of the magnetic well and hill regions. It was found that the magnetic hill regions play a critical role in enhancing electron thermal transport. When the plasma collapses to the wall, the 3D ambipolar potential arises in the stochastic layer to maintain the quasi-neutrality of the plasma. The E × B vortices induced by the 3D ambipolar potential mix the plasma across stochastic field lines and enhance the radial transport. Particularly, the E × B mixing between the magnetic well and hill regions provides a collisionless detrapping mechanism that plays a major role in the loss of high-energy trapped electrons. As a result, the electron temperature steadily decreases at a rate of about -0.5 keV ms -1 , comparable to the rate observed in thermal quench experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Neural net modeling of equilibria in NSTX-U

Neural networks (NNs) offer a path towards synthesizing and interpreting data on faster timescales than traditional physics-informed computational models. In this work we develop two NNs relevant to equilibrium and shape control modeling, which are part of a suite of tools being developed for the National Spherical Torus Experiment-Upgrade for fast prediction, optimization, and visualization of plasma scenarios. The networks include Eqnet, a free-boundary equilibrium solver trained on the EFIT01 (Equilibrium FITtting 01) reconstruction algorithm, and Pertnet, which is trained on the Gspert code and predicts the non-rigid plasma response, a nonlinear term that arises in shape control modeling. The NNs are trained with different combinations of inputs and outputs in order to offer flexibility in use cases. In particular, Eqnet can use magnetic diagnostics as inputs and act as an EFIT-like reconstruction algorithm, or, by using pressure and current profile information the NN can act as a forward Grad–Shafranov equilibrium solver. This forward-mode version is envisioned to be implemented in the suite of tools for simulation of plasma scenarios. The reconstruction-mode version gives some performance improvements compared to the online reconstruction code real-time EFIT, especially when vessel eddy currents are significant. Here, we report strong performance for all NNs indicating that the models could reliably be used within closed-loop simulations or other applications. Some limitations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulations of divertor plasmas with inverse sheaths

The effect of strong electron emission from material surfaces has been proposed to form an “inverse sheath”: a region with a positive potential relative to the near-wall plasma which prevents the flow of ions to the wall [M. D. Campanell, “Negative plasma potential relative to electronemitting surfaces,” Phys. Rev. E. 88, 033103 (2013); M. D. Campanell and M. V. Umansky, “Strongly emitting surfaces unable to float below plasma potential,” Phys. Rev. Lett. 116, 1–5 (2016); M. D. Campanell and G. R. Johnson, “Thermionic cooling of the target plasma to a sub-ev temperature,” Phys. Rev. Lett. 122, 1–5 (2019)]. We assess the viability of this regime in a tokamak device using the 2D edge plasma transport code UEDGE [T. Rognlien et al ., “A fully implicit, time dependent 2-D fluid code for modeling tokamak edge plasmas,” J. Nucl. Mater. 196–198, 347–351 (1992)]. Since the UEDGE code does not consider the sheath region directly, we apply boundary conditions at the divertor targets which emulate the physics of both “standard” and “inverse” sheath regimes [R. Masline et al ., “Influence of the inverse sheath on divertor plasma performance in tokamak edge plasma simulations,” Contrib. Plasma Phys. 60, e201900097 (2020)]. Using these boundary conditions, we perform scoping studies to assess plasma parameters near the target by varying the density at the core-edge interface. We observe a smooth transition in the resultant profiles of plasma parameters for the standard sheath, and a bifurcation across the simulation set for plasmas with an inverse sheath. The cause of this bifurcation is assessed by performing the parameter scan both with and without impurity radiation; we observe that the bifurcation persists in both cases, indicating that this bifurcation is caused by plasma recombination.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Reimagining full wave rf quasilinear theory in a tokamak

The velocity dependent resonant interaction of particles with applied radiofrequency (rf) waves during heating and current drive in the presence of pitch angle scattering collisions gives rise to narrow collisional velocity space boundary layers that dramatically enhance the role of collisions as recently shown by Catto ( J. Plasma Phys. , vol. 86, 2020, 815860302). The behaviour is a generalization of the narrow collisional boundary layer that forms during Landau damping as found by Johnston ( Phys. Fluids , vol. 14, 1971, pp. 2719–2726) and Auerbach ( Phys. Fluids , vol. 20, 1977, pp. 1836–1844). For a wave of parallel wave number ${k_{||}}$ interacting with weakly collisional plasma species of collision frequency $\nu$ and thermal speed ${v_{\textrm{th}}}$ , the effective collision frequency becomes of order $\nu {({k_{||}}{v_{th}}/\nu )^{2/3}} \gg \nu $ . The narrow boundary layers that arise because of the diffusive nature of the collisions allow a physically meaningful wave–particle interaction time to be defined that is the inverse of this effective collision frequency. The collisionality implied by the narrow boundary layer results in changes in the standard quasilinear treatment of applied rf fields in tokamaks while remaining consistent with causality. These changes occur because successive poloidal interactions with the rf are correlated in tokamak geometry and because the resonant velocity space dependent interactions are controlled by the spatial and temporal behaviour of the perturbed full wave fields rather than just the spatially local Landau and Doppler shifted cyclotron wave–particle resonance condition associated with unperturbed motion of the particles. The correlation of successive poloidal circuits of the tokamak leads to the appearance in the quasilinear operator of transit averaged resonance conditions localized in velocity space boundary layers that maintain negative definite entropy production.

Physics↗
DOE OSTI2021

Experiment and model for a Stokes layer in a strongly coupled dusty plasma

A Stokes layer, which is a flow pattern that arises in a viscous fluid adjacent to an oscillatory boundary, was observed in an experiment using a two-dimensional strongly coupled dusty plasma. Liquid conditions were maintained using laser heating, while a separate laser manipulation applied an oscillatory shear that was localized and sinusoidal. The evolution of the resulting flow was analyzed using space-time diagrams. These figures provide an intuitive visualization of a Stokes layer, including features such as the depth of penetration and wavelength. Another feature, the characteristic speed for the penetration of the oscillatory flow, also appears prominently in space-time diagrams. To model the experiment, the Maxwell-fluid model of a Stokes layer was generalized to describe a two-phase liquid. In our experiment, the phases were gas and dust, where the dust cloud was viscoelastic due to strong Coulomb coupling. Furthermore, the model is found to agree with the experiment, in the appearance of the space-time diagrams, and in the values of the characteristic speed, depth of penetration, and wavelength.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of CMAS viscosity on the infiltration depth in thermal barrier coatings of different microstructures

Calcium magnesium aluminum silicate (CMAS) is one of the leading concerns for the gas turbine industry. The effects of CMAS viscosity and the coating microstructure on CMAS infiltration depth were explored by conducting a time dependent interaction study. Three CMAS compositions were used from literature, and their viscosities predicted through FactSage viscosity module were drastically different. The interaction was carried out on three different TBCs synthesized using the solution precursor plasma spray process (SPPS): two of the TBCs were made of yttrium aluminum garnet (YAG) having different microstructures that promote different modes of CMAS infiltration, and one TBC was made of gadolinium zirconate (GZO). All samples had stress relieving vertical cracks and different intensities of horizontally banded porosity, (inter pass boundaries IPBs). A concentration of 100 mg/cm 2 of CMAS was applied on the TBCs which were then subjected to a 5-minute interaction at 1300 °C. Samples were analyzed using scanning electron microscopy (SEM), electron dispersive X-ray spectroscopy (EDXS), and transmission electron microscopy (TEM). Low viscosity CMAS readily penetrated the TBCs while more viscous CMAS showed less penetration. The depth of CMAS infiltration depended on the coating microstructure. In the YAG with IPBs, the CMAS spread horizontally in the IPBs before infiltrating deeper, resulting in reduced infiltration depth compared to other samples in spite of having wider vertical cracks. TEM and EDXS analysis were performed to investigate the phases present in the CMAS-TBC interaction region in YAG. Two regions were chosen, the top TBC surface in direct contact with the sea of CMAS, and the region at the CMAS penetration ended within the coating. Finally, the results showed that no secondary phases like apatite were observed in YAG, thus it can be concluded that the arrest of CMAS happened solely because of CMAS viscosity and the short infiltration time.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Parametric dependencies of resonant layer responses across linear, two-fluid, drift-MHD regimes

Non-axisymmetric magnetic fields arising in a tokamak either by external or internal perturbations can induce complex non-ideal MHD responses in their resonant surfaces while remaining ideally evolved elsewhere. This layer response can be characterized in a linear regime by a single parameter called the inner-layer Δ, which enables outer-layer matching and the prediction of torque balance to non-linear island regimes. In this work, we follow strictly one of the most comprehensive analytic treatments including two-fluid and drift MHD effects and keep the fidelity of the formulation by incorporating the numerical method based on the Riccati transformation when quantifying the inner-layer Δ. The proposed scheme reproduces not only the predicted responses in essentially all asymptotic regimes but also with continuous transitions as well as improved accuracies. In particular, the Δ variations across the inertial regimes with viscous or semi-collisional effects have been further resolved, in comparison with additional analytic solutions. The results imply greater shielding of the electromagnetic torque at the layer than what would be expected by earlier work when the viscous or semi-collisional effects can compete against the inertial effects, and also due to the intermediate regulation by kinetic Alfvén wave resonances as rotation slows down. These are important features that can alter the non-axisymmetric plasma responses including the field penetration by external fields or island seeding process in rotating tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Effects of wave damping and finite perpendicular scale on three-dimensional Alfvén wave parametric decay in low-beta plasmas

Shear Alfvén wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how the finite perpendicular wave scale—as found in both laboratory and space plasmas—affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfvén wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron–ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfvén wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role: waves of different perpendicular scales exhibit similar instability excitation as long as the magnitude of the parallel ponderomotive force remains unchanged. These findings are corroborated by theoretical analysis and estimates. This new understanding of 3D kinetic Alfvén wave PDI physics is essential for laboratory study of the basic plasma process and may also aid future evaluation of the relevance/role of PDI in low-beta space plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Parallel Cut-Cell Algorithm for the Free-Boundary Grad--Shafranov Problem

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad--Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. Here the algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad--Shafranov solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

First snowflake divertor experiments in MAST-U tokamak

First snowflake (SF) divertor experiments in the MAST-U tokamak demonstrated steady-state snowflake-plus divertor configurations in 450 kA ohmic L-mode plasmas. The SF divertor configuration features a second poloidal field (PF) null in the divertor region close by or overlapping with the main X-point. The resulting low PF region and two additional divertor legs (strike points) may lead to additional power and particle flux sharing via a hypothesized convective cell, and increased plasma-wetted area and radiation. The free-boundary Grad–Shafranov equilibrium code FIESTA was used to design SF configurations with several inter-null distances and orientations. In the experiment, the SF configurations with inter-null distances 0.13–0.20 m and lasting 0.2–0.3 s were obtained. Parallel connection lengths between the midplane and the outer strike point in the SF configurations (evaluated at field lines 1–2 mm from the separatrix in the midplane) were 25–30 m, higher than in the standard divertor (20–25 m) or the Super-X divertor (25 m). Diagnostic measurements highlighted salient SF features. The infra-red video bolometer diagnostic showed that the radiated power peaking in the PF null region was not as pronounces as in the standard divertor. Divertor ion fluxes measured by target Langmuir probes showed increased ion flux in the plate region where a secondary SF strike point landed, concomitantly with the SF formation. These measurements may suggest that some particle and heat redistribution was taking place in the convective SF zone. The first SF experiments provide a basis for future SF studies in MAST-U tokamak with higher input power, improved plasma control and diagnostic measurements, to be compared with the modeling predictions of plasma convective SF mixing and lower density strike point detachment threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

The boundary density profile of a Coulomb droplet. Freezing at the edge

In this paper, we revisit the problem of computing the boundary density profile of a droplet of two-dimensional one-component plasma (2D OCP) with logarithmic inter- action between particles in a confining harmonic potential. At a sufficiently low temperature, but still in the liquid phase, the density exhibits oscillations as a function of the distance to the boundary of the droplet. We obtain the density profile numerically using Monte-Carlo simulations of the 2D OCP. We argue that the decay and period of those oscillations can be explained within a picture of the Wigner crystallization near the boundary, where the crystal is gradually melted with the increasing distance to the boundary.

2D OCP↗
DOE OSTI2023

Using the Stix finite element RF code to investigate operation optimization of the ICRF antenna on Alcator C-Mod

Abstract As the Ion Cyclotron Radio Frequency range (ICRF) heating becomes more favorable in fusion devices, the urgency of predicting and mitigating impurity generation that arises from it becomes more pressing. In the ICRF regime, rectified Radio Frequency (RF) sheaths are known to form at antenna and material edges that influence negative effects like sputtering and a decrease in heating efficiency. Methods to mitigate the formation of these RF sheaths through RF image currents cancellation have been experimentally studied. A power-phasing scan done on Alcator C-Mod in which the amount of power on the two inner straps ( P in ) versus the total 4 straps ( P tot ) was varied showed a minimization of enhanced potentials between P in / P tot ∼ 0.7–0.9 while impurities were minimized for P in / P tot ∼ 0.5–0.8. New capabilities in the realm of representing the RF sheath numerically now allow for these experiments to be simulated. Given the size of the sheath relative to the scale of the device, it can be approximated as a Boundary Condition (BC). A new parallelized cold-plasma wave equation solver called Stix implements a non-linear sheath impedance model BC formulated by Myra et al (2015 Phys. Plasmas 22 062507) through the method of finite elements using the MFEM library [ http://mfem.org ]. It is seen that Stix shows qualitative agreement with the measured C-Mod enhanced potentials.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The role of electric dominance for particle injection in relativistic reconnection

ABSTRACT Magnetic reconnection in relativistic plasmas – where the magnetization $\sigma \gg 1$ – is regarded as an efficient particle accelerator, capable of explaining the most dramatic astrophysical flares. We employ two-dimensional (2D) particle-in-cell simulations of relativistic pair-plasma reconnection with vanishing guide field and outflow boundaries to quantify the impact of the energy gain occurring in regions of electric dominance ($E\gt B$) for the early stages of particle acceleration (i.e. the ‘injection’ stage). Given an injection threshold energy $\epsilon ^\ast =\sigma /4$ for the particles that eventually attain energy $\epsilon _{\rm T}$ by time T, we calculate the mean fractional contribution $\zeta (\epsilon ^\ast ,\epsilon _{\rm T})$ by $E\gt B$ fields to particle energization at the time when the threshold $\epsilon ^\ast$ is reached. We find that $\zeta$ monotonically increases with $\sigma$ and $\epsilon _{\rm T}$; for $\sigma \gtrsim 50$ and $\epsilon _{\rm T}/\sigma \gtrsim 8$, we find that $\gtrsim 80~{{\ \rm per\ cent}}$ of the energy gain obtained before reaching $\epsilon ^\ast =\sigma /4$ occurs in $E\gt B$ regions. We find that $\zeta$ is independent of simulation box size $L_x$, as long as $\epsilon _{\rm T}$ is normalized to the maximum particle energy, which scales as $\epsilon _{\rm max}\propto L_{\rm x}^{1/2}$ in 2D. The distribution of energy gains $\epsilon _{\chi }$ acquired in $E\gt B$ regions can be modelled as $\mathrm{ d}N/\mathrm{ d}\epsilon _{\chi }\propto \epsilon _{\chi }^{-0.35}\exp [-(\epsilon _{\chi }/0.06\, \sigma)^{0.5}]$. Our results help assess the role of electric dominance in relativistic reconnection with vanishing guide fields, which is realized in the magnetospheres of black holes and neutron stars.

Gupta, Sanya (ORCID:0000000151944384)↗
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗