Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Disruptions”

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,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Environmental exposure to industrial air pollution is associated with decreased male fertility

Objective: To understand how chronic exposure to industrial air pollution is associated with male fertility through semen parameters. Design: Retrospective cohort study. Subjects: Men in the Subfertility, Health and Assisted Reproduction cohort who underwent a semen analysis 2005-2017 with ≥1 measured semen parameter (N=21,563). Intervention(s): Residential histories for each man were constructed using locations from administrative records linked through the Utah Population Database. Industrial facilities with air emissions of nine endocrine disrupting compound chemical classes were identified from the Environmental Protection Agency Risk-Screening Environmental Indicators microdata. Chemical levels were linked with residential histories for the 5 years prior to each semen analysis. Main Outcome Measures: Semen analyses were classified as azoospermic or oligozoospermic (< 15 M/mL) using World Health Organization cutoffs for concentration. Bulk semen parameters such as concentration, total count, ejaculate volume, total motility, total motile count, and total progressive motile count were also measured. Multivariable regression models with robust standard errors were used to associate exposure quartiles for each of the nine chemical classes with each semen parameter, adjusting for age, race, and ethnicity, as well as neighborhood socioeconomic disadvantage. Results: After adjustment for demographic covariates, several chemical classes were associated with azoospermia and decreased total motility and volume. For exposure in the 4th relative to 1st quartile, significant associations were observed for acrylonitrile (β total motility = -0.87 pp), aromatic hydrocarbons (odds ratio [OR]azoospermia = 1.53; β volume = -0.14 mL), dioxins (OR azoospermia = 1.31; β volume = -0.09 mL; β total motility = -2.65 pp), heavy metals (β total motility = -2.78pp), organic solvents (OR azoospermia = 1.75; β volume = -0.10 mL), organochlorines (OR azoospermia = 2.09; β volume = -0.12 mL), phthalates (OR azoospermia = 1.44; β volume = -0.09 mL; β total motility = -1.21 pp), and silver particles (OR azoospermia = 1.64; β volume = -0.11 mL). All semen parameters significantly decreased with increasing socioeconomic disadvantage. Men who lived in the most disadvantaged areas had concentration, volume, and total motility of 6.70 M/mL, 0.13 mL, and 1.79 pp lower, respectively. Count, motile count, and total progressive motile count all decreased by 30–34 M.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

Risk-averse optimization for resilience enhancement of complex engineering systems under uncertainties

With the growth of complexity and extent, large scale interconnected network systems, e.g., transportation networks or infrastructure networks, become more vulnerable to external disturbances. Hence, managing potential disruptive events during the design, operating, and recovery phase of an engineered system and therefore improving the system’s resilience is an important yet challenging task. Here, to ensure system resilience after the occurrence of failure events, this study proposes a mixed-integer linear programming (MILP) based restoration framework using heterogeneous dispatchable agents. The scenario-based stochastic optimization (SO) technique is adopted to deal with the inherent uncertainties imposed on the recovery process from nature. Moreover, different from conventional SO using deterministic equivalent formulations, the CVaR risk measure is implemented for this study because of the temporal sparsity of the decision making in applications such as the recovery from extreme events. The resulting restoration framework involves a large-scale MILP problem and thus an adequate decomposition technique i.e. modified Lagrangian dual decomposition, is also employed to achieve tractable computational complexity. Case study results based on the IEEE 37-bus test feeder demonstrate the benefits of using the proposed framework for resilience improvement as well as the advantages of adopting SO formulations.

42 ENGINEERING↗
DOE OSTI2023

A physics-informed deep learning model of the hot tail runaway electron seed

A challenging aspect of the description of a tokamak disruption is evaluating the hot tail runaway electron seed that emerges during the thermal quench. This problem is made challenging due to the requirement of describing a strongly non-thermal electron distribution, together with the need to incorporate a diverse range of multiphysics processes, including magnetohydrodynamic instabilities, impurity transport, and radiative losses. Here this work develops a physics-informed neural network (PINN) tailored to the solution of the hot tail seed during an axisymmetric thermal quench. Here, a PINN is developed to identify solutions to the adjoint relativistic Fokker–Planck equation in the presence of a rapid quench of the plasma's thermal energy. It is shown that the PINN is able to accurately predict the hot tail seed across a range of parameters, including the thermal quench timescale, initial plasma temperature, and local current density, in the absence of experimental or simulation data. The hot tail PINN is verified by comparison with direct Monte Carlo simulations, with excellent agreement found across a broad range of thermal quench conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Resurfacing promotes antibacterial activity of a lipid A–binding nanobody

Nanobodies have been pursued as candidates for antimicrobial design due to their small size and versatile binding capacities, but direct antibacterial activity of a nanobody has yet to be described. Here, we employed a bacterial surface display platform to screen a synthetic library of nanobody variants for antimicrobial potential. We identified a candidate that binds the essential lipid A component of gram-negative lipopolysaccharide. Nonetheless, this nanobody required a weakened outer membrane to access its target and elicit its toxic activity. Borrowing from observations of innate immune proteins, we found that resurfacing nanobodies with positively charged residues enabled them to bind and perturb the gram-negative outer membrane, but this alone was not sufficient for toxic activity. However, when we resurface our lipid A-targeting nanobody, it gained the ability to disrupt the outer membrane and enact its antibacterial function against wild-type bacteria. This development of a dual-function nanobody that can reach and bind previously inaccessible gram-negative targets introduces a route for antimicrobial biologic advancement.

antibacterial↗
DOE OSTI2020

Analysis of the Shattered Pellet Injection Fragment Plumes Generated by Machine Specific Shatter Tube Designs

Shattered pellet injector systems have been installed on DIII-D, JET, and KSTAR and used to experimentally determine the effectiveness of the shattered pellet injection (SPI) process in mitigating the deleterious effects of a tokamak plasma disruption. Pellets are fired, and before entering the plasma, strike a bent tube known as a shatter tube causing the pellet to shatter. The process of pellet fragmentation is a chaotic process that can be described in terms of fragment size distribution through a statistical model that incorporates the effects of the pellet material and impact characteristics. In addition to the fragment size distribution, the shatter plume has other characteristics of interest, such as a fragment velocity distribution and temporal mass evolution. The fragment velocity distribution is important because it is needed to accurately model the spread and location of the ablation and the deposition of impurities in the plasma over time. The temporal mass evolution is necessary to determine the time-resolved delivery of mass to the plasma.Due to installation constraints, the shatter tube currently installed on JET has a unique geometry with a modest S-bend followed by a 20-deg bend at the end of the tube. The DIII-D and KSTAR shatter tube design is a simple tube bent through an angle of 20 deg followed by a straight section. The resulting shatter sprays from the JET shatter tube and a 20-deg miter bend shatter tube were experimentally characterized for various pellet materials and speeds. Laboratory testing of these shatter tubes allows for the use of fast cameras to capture the fragment spray traveling through a large vacuum chamber. These high-speed videos of the shatter plumes allow the fragment size distribution, temporal mass evolution, and velocity distribution of the fragments within the plume to be determined. Overall, this paper presents a comparison of the unique geometry of the JET shatter tube to the miter bend geometries used for shattering and some insight into the variables that may be adjusted to produce the optimal shatter spray. The impact of entrained propellant gas on the resulting shatter spray was examined during testing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Testing of Fusion Machine Components Against Cryogenic Pellet Impacts

Shattered pellet injection (SPI), currently the most effective method of disruption mitigation, is currently implemented on tokamaks worldwide for experimental purposes. Cryogenic pellets are formed and fired into an angled surface before entering the plasma. The impact with the angled surface causes the pellets to fragment into a cloud of particles with the purpose of increasing the surface area for ablation. As pellets traverse guide tubes, depending on design, there is a chance of an off-normal pellet impact. Pellet impacts are also, depending on design, possible in the plasma chamber if not fully ablated and assimilated, or if the fragment plume is not directed in the proper direction. This paper outlines a series of pellet impact tests on various tiles and components relevant to the ITER, JET, and ASDEX Upgrade SPI systems. In conclusion, testing was done to assess the potential damage from pellet and fragment plume impacts through high-speed imaging and the visual inspection of components.

Cryogenic pellets↗
DOE OSTI2020

Projected global stability of high beta MAST-U spherical tokamak plasmas

Assessment of the limits of stability of tokamak plasmas is essential to operation in high fusion performance ranges without disruption of the plasma current. Projected equilibria have been generated for the MAST-U spherical tokamak experiment, an upgrade of the previous MAST device, in order to prepare for operation. These equilibria are scanned in pressure and current profiles, and assessed with the DCON and MARS-F stability codes to find the so-called 'no-wall' beta limit, above which resistive wall mode instabilities can be expected in the absence of other stabilising effects. The no-wall limit was generally found to increase as plasma internal inductance increased. The equilibria are also assessed for the 'with-wall' limit, theoretically the highest achievable performance point, again with the DCON and MARS-F codes, including different approximate axisymmetric walls, and with the VALEN code which includes a 3D model of the surrounding conducting structure. Similar limits were found, despite the difference between the 2D and 3D codes in the treatment of the wall. Conducting passive stabilisation plates, which were newly installed in MAST-U, are in a region of significant mode perturbation when the plasma β N is sufficiently high and eddy currents are driven in these structures. Due to the increased stabilising effect of the wall in MAST-U vs. MAST, a significant gap exists between the approximate no-wall limits of β N / l i = 6.71 and 7.13, found from DCON and MARS-F respectively, and the with-wall limits of β N / l i = 8.23 and 8.53 for the equilibrium profiles analysed in this study. This opens a region of high beta operating space in MAST-U for potentially stable operation if non-ideal effects or active control can stabilise the resistive wall mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Simulation of shell pellet injection strategies for ITER-scale tokamaks

Dispersive shell pellet (DSP) injection is considered as an alternative to shattered pellet injection as a disruption mitigation system for ITER, and strategies for penetration of a shell pellet into ITER are modeled with the 3D MHD code NIMROD. Because the high plasma temperatures lead to rapid ablation of the shell, delivery of the dispersive payload to the core of ITER will be very challenging. Two strategies to increase payload delivery depth are modeled: first, multiple staggered pellets are simulated in DIII-D, to assess the ability for one DSP to "piggy-back" on another to reach deeper into the core; second, DSP injection after pre-dilution-cooling with deuterium is simulated in ITER, in order to reduce the plasma temperature before shell pellet arrival. The DIII-D simulations show that a second, slower pellet can penetrate much deeper once the release of the first payload strongly cools the mid-radius region. When the pellets are staggered, deeper penetration of the second pellet leads to higher radiation fraction and larger runaway electron loss fraction, consistent with single pellet results. However, simultaneously released pellets at mid-radius that do not trigger a large n=1 mode produce an even higher radiation fraction. The ITER simulations show that an inside-out TQ can be produced with a payload release just inside of the q=2 surface, which is achieved at a speed of 800 m/s after pre-dilution cooling. Although stochastization of the core leads to a complete thermal quench, the edge flux surfaces are surprisingly robust in the ITER simulations, regardless of payload release location. As a result, runaway electron losses would not be expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Interpretive MHD modeling of dispersive shell pellet injection for rapid shutdown in tokamaks

Dispersive shell pellet (DSP) injection is modeled with the extended-MHD code NIMROD for interpretive insight into the results of recent DIII-D DSP experiments and to explore the dynamics of an inside-out thermal quench for disruption mitigation in tokamaks. Simulations of the pre-thermal quench (TQ) phase indicate that the upper bound for the quantity of ablated carbon shell material that will not perturb the flux surfaces is in the ballpark of, but somewhat below the experimental quantity. Even below this quantity, sufficient electrons are added to the plasma by the shell material to produce significant dilution cooling before the TQ is triggered. Simulations carried through the end of the TQ have very large amplitude MHD fluctuations (δB/B>10 -2 ) at the time of the plasma current spike associated with current profile redistribution. Finally, after the plasma current spike, which is of comparable amplitude to that measured in DIII-D experiments, none of the runaway electron test-particles whose orbits are tracked throughout the simulation remain confined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A semi-supervised machine learning detector for physics events in tokamak discharges

Databases of physics events have been used in various fusion research applications, including the development of scaling laws and disruption avoidance algorithms, yet they can be time-consuming and tedious to construct. This paper presents a novel application of the label spreading semi-supervised learning algorithm to accelerate this process by detecting distinct events in a large dataset of discharges, given few manually labeled examples. A high detection accuracy (> 85%) for H-L back transitions and initially rotating locked modes is demonstrated on a dataset of hundreds of discharges from DIII-D with manually identified events for which only 3 discharges are initially labeled by the user. Lower yet reasonable performance (~75%) is also demonstrated for the core radiative collapse, an event with a much lower prevalence in the dataset. Additionally, analysis of the performance sensitivity indicates that the same set of algorithmic parameters is optimal for each event. Furthermore, this suggests that the method can be applied to detect a variety of other events not included in this paper, given that the event is well described by a set of 0D signals robustly available on many discharges. Procedures for analysis of new events are demonstrated, showing automatic event detection with increasing fidelity as the user strategically adds manually labeled examples. Detections on Alcator C-Mod and EAST are also shown, demonstrating the potential for this to be used on a multi-tokamak dataset.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Lagrangian particle model for 3D simulation of pellets and SPI fragments in tokamaks

A 3D numerical model for the ablation of pellets and shattered pellet injection (SPI) fragments in tokamaks in the plasma disruption mitigation and fueling parameter space has been developed based on the Lagrangian particle code [R. Samulyak, X. Wang, H.-S. Chen, Lagrangian Particle Method for Compressible Fluid Dynamics, J. Comput. Phys., 362 (2018), 1-19]. The pellet code implements the low magnetic Reynolds number MHD equations, kinetic models for the electronic heating, a pellet surface ablation model, an equation of state that supports multiple ionization states, radiation, and a model for grad-B drift of the ablated material across the magnetic field. The Lagrangian particle algorithm is highly adaptive, capable of simulating a large number of fragments in 3D while eliminating numerical difficulties of dealing with the tokamak background plasma. The code has achieved good agreement with theory for spherically symmetric ablation flows. Axisymmetric simulations of neon and deuterium pellets in magnetic fields ranging from 1 to 6 Tesla have been compared with previous simulations using the FronTier code, and very good agreement has also been obtained. Furthermore, the main physics contribution of the paper is a detailed study of the influence of 3D effects, in particular grad-B drift, on pellet ablation rates and properties of ablation clouds. Smaller reductions of ablation rates in magnetic fields compared to axially symmetric simulations have been demonstrated because the ablated material is not confined to narrowing channels in the presence of grad-B drift. Contribution of various factors in the grad-B drift model has also been quantified.

97 MATHEMATICS AND COMPUTING↗
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

Advances in physics and applications of 3D magnetic perturbations on the J-TEXT tokamak

In the last two years, three major technical improvements have been made on J-TEXT in supporting of the expanded operation regions and diagnostic capabilities. (1) The successful commission of the 105 GHz/500 kW/1 s electron cyclotron resonance heating (ECRH) system increasing the core electron temperature from 0.9 keV up to around 1.5 keV. (2) The poloidal divertor configuration with an X-point in the high-field side has been achieved. In particular, the 400 kW electron cyclotron wave has also been successfully injected into the diverted plasma. (3) A 256-channel electron cyclotron emission imaging diagnostic system and two sets of four-channel Doppler backscattering diagnostics have been successfully developed on J-TEXT, allowing detailed measurement of the electron temperature and density fluctuations for turbulence and MHD research. The locked mode (LM), especially the 2/1 LM, is one of the biggest threats to the plasma operation. Both the thresholds of 2/1 and 3/1 LM are observed to vary non-monotonically on electron density. The electrode biasing was applied successfully to unlock the LM from either a rotating or static resonant magnetic perturbation (RMP) field. In the presence of 2/1 LM, three kinds of standing wave (SW) structures have been observed to share a similar connection to the island structure, i.e. the nodes of the SWs locate around the O- or X-points of the 2/1 island. The control and mitigation of disruption is essential to the safe operation of ITER, and it has been systematically studied by applying a RMP field, massive gas injection (MGI) and shattered pellet injection on J-TEXT. When the RMP-induced 2/1 LM is larger than a critical width, the MGI shutdown process can be significantly influenced. If the phase difference between the O-point of LM and the MGI valve is +90° (or –90°), the penetration depth and the assimilation of impurities can be enhanced (or suppressed) during the pre-thermal quench (TQ) phase and result in a faster (or slower) TQ. A secondary MGI can also suppress the runaway electron (RE) generation, if the additional high-Z impurity gas arrives at the plasma edge before TQ. When the secondary MGI has been applied after the formation of the RE current plateau, the RE current can be dissipated, and the dissipation rate increases with the injected impurity quantity but saturates with a maximum of 28 MA s –1 . The non-local transport is experimentally observed in the ion transport channel. The electron thermal diffusivity significantly increases with the ECRH power. Theoretical work shows that significant intrinsic current can be driven by electromagnetic turbulence, and the robust formation mechanism of the E × B staircase is identified from the Hasegawa–Wakatani system.

divertor↗
DOE OSTI2021

MHD simulations of cold bubble formation from 2/1 tearing mode during massive gas injection in a tokamak

Massive gas injection (MGI) experiments have been carried out in many tokamaks to study disruption dynamics and mitigation schemes. Two events often observed in those experiments are the excitation of the m = 2, n = 1 magnetohydrodynamic mode, and the formation of cold bubble structure in the temperature distribution before the thermal quench (TQ). Here m is the poloidal mode number, n the toroidal mode number. The physics mechanisms underlying those phenomena, however, have not been entirely clear. In this work, our recent NIMROD simulations of the MGI process in a tokamak have reproduced the main features of both events, which has allowed us to examine and establish the causal relation between them. In these simulations, the 3/1 and 2/1 islands are found to form successively after the arrival of impurity ion cold front at the corresponding q = 3 and q = 2 rational surfaces. At the interface between impurity and plasma, a local thin current sheet forms due to an enhanced local pressure gradient and moves inward following the gas cold front, this may contribute to the formation of a dominant 2/1 mode. Following the growth of the 2/1 tearing mode, the impurity penetration into the core region inside the q = 2 surface gives rise to the formation of the cold bubble temperature structure and initiates the final TQ. Here, a subdominant 1/1 mode developed earlier near the q = 1 surface alone does not cause such a cold bubble formation, however, the exact manner of the preceding impurity penetration depends on the nature of the 1/1 mode: kink-tearing or quasi-interchange.

2/1 tearing mode↗
DOE OSTI2022

Hot-tail electrons’ impact on assimilation and injection penetration of D2 shattered pellet injections

The assimilation and the penetration of injected materials play significant roles in the mitigation efficiency of shattered pellet injection (SPI) as a disruption mitigation system (DMS), hence an accurate description of the fragment ablation rate is desirable. Current mainstream 3D MHD codes modelling SPIs mostly assume instantaneous thermalization between the previously hot ambient electrons and the newly released cold electrons. This could result in underestimation of the ablation rate if the hot electron thermalization time is comparable or even longer than the fragment flying time across flux tubes, putting the validity of these 3D SPI simulations into question. To resolve this doubt, we, in this study, investigate the thermalization dynamics and the overall hot-electron impact. The finite-time collisional thermalization of hot-tail electrons in a rapidly cooling plasma, as well as the so-called ‘self-limiting’ effect are considered. The former effect tends to deplete the colder population within a hot-tail species, while the latter is found to preferentially deplete the higher energy population. The combined result is found to cause an almost self-similar decay of the hot electron distribution function, while its shape does not deviate much from that of Maxwellian distribution and the mean energy does not change much during the thermalization process. Based on this observation, axisymmetric JOREK D2 SPI simulations were carried out with additional hot-tail contribution to evaluate their overall impact onto the injection assimilation and penetration. It is found that the hot-tail effect indeed causes enhanced assimilation and shallower penetration, although the overall effect depends on the exact injection configuration, with the slow injection showing negligible hot-tail effect while the fast single non-shattered pellet case shows drastic hot-tail ablation enhancement. For ITER-like SPI parameters, there is no significant deviation in the total assimilation, but some deviation in the injection penetration is observed for the fast injection velocity cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

MHD modeling of shattered pellet injection in JET

Abstract Nonlinear 3D MHD simulations of shattered-pellet injection (SPI) in JET show prototypical SPI-driven disruptions using the M3D-C1 and NIMROD extended-MHD codes. Initially, radiation-driven thermal quenches are accelerated by MHD activity as the pellet crosses rational surfaces, leading to a radiation spike, global stochasticization of the magnetic field, and a complete thermal quench. Eventually, current quenches, preceded by a current spike are seen as the Ohmic heating becomes equal to the radiative cooling. The results are qualitatively similar for both a single monolithic pellet, pencil-beam model, and a realistic shatter to represent the SPI plume. A scan in viscosity from 500 to 2000 m 2 s −1 for MHD simulations finds that reducing viscosity increases MHD activity and decreases thermal quench time slightly. A realistic cloud of fragments modeling shows that mixed-D–Ne pellet travels deeper into the plasma core before the thermal quench. At the slow pellet speeds, the pellet is found to be moving slowly enough inward that even the 5% neon in the mixed pellet is enough to effectively radiate the thermal energy available. Radiation toroidal peaking is predicted to be at levels consistent with experimental observations and reduced as the pellet travels deeper into the plasma. These simulations lay the ground work for more-sophisticated validative and predictive modeling of SPI in JET using both M3D-C1 and NIMROD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Staged cooling of a fusion-grade plasma in a tokamak thermal quench

Abstract In tokamak disruptions where the magnetic connection length becomes comparable to or even shorter than the plasma mean-free-path, parallel transport can dominate the energy loss and the thermal quench of the core plasma goes through four phases (stages) that have distinct temperature ranges and durations. The main temperature drop occurs while the core plasma remains nearly collisionless, with the parallel electron temperature T e ∥ dropping in time t as T e ∥ ∝ t − 2 and a cooling time that scales with the ion sound wave transit time over the length of the open magnetic field line. These surprising physics scalings are the result of effective suppression of parallel electron thermal conduction in an otherwise bounded, quasineutral, and collisionless plasma, which is different from what are known to date on electron thermal conduction along the magnetic field in a nearly collisionless and quasineural plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Wall heating by subcritical energetic electrons generated by the runaway electron avalanche source *

Abstract Subcritical energetic electrons (SEEs) produced by the runaway electron (RE) avalanche source at energies below the runaway threshold are found to be the primary contributor to surface heating of plasma-facing components (PFCs) during final loss events. This finding is supported by theoretical analysis, computational modeling with the Kinetic Orbit Runaway electrons Code (KORC), and qualitative agreement with DIII-D experimental observations. The avalanche source generates significantly more secondary electrons below the runaway threshold, which thermalize rapidly when well-confined. However, during a final loss event, the RE beam impacts the first wall, and SEEs are deconfined before they can thermalize. Additionally, because the energy deposition length decreases faster than energy, the deposited energy density, and thus the maximum PFC surface temperature change, is larger for SEEs than REs. KORC simulations employ an analytic first wall to model particle deconfinement onto a non-axisymmetric wall composed of individual tiles. PFC surface heating is calculated using a 1D model extended to include an energy-dependent deposition length scale. Simulations of DIII-D qualitatively agree with infrared (IR) imaging only when SEEs from the avalanche source are included. These results demonstrate that SEEs are the dominant contributor to PFC surface heating and indicate that the avalanche source plays a critical role in the PFC damage caused during final loss events. The prominence of SEEs also has important implications for interpreting IR imaging, one of the primary diagnostics for RE-wall interaction diagnosis, despite REs dominating the energy and current density. This result improves predictions of wall damage due to post-disruption REs to estimate material lifetime and design RE mitigation systems for ITER and future reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗