Engineering PapersSearch

DOE OSTI · 2537932

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

Hollmann, E. M. [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000262676589)·Marini, Claudio [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000348197765)·Rudakov, D. L. [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000252664269)·Martinez-Loran, E. [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000338414981)·Beidler, Matthew [Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)] (ORCID:0000000273853886)·Herfindal, Jeffrey L. [Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)] (ORCID:000000032846597X)·Shiraki, Daisuke [Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)]·Liu, Yueqiang Q. [General Atomics, San Diego, CA (United States)] (ORCID:0000000281928411)·Bykov, I. [General Atomics, San Diego, CA (United States)]·Lvovskiy, Andreas [General Atomics, San Diego, CA (United States)] (ORCID:0000000236491169)·Lasnier, C. [Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)] (ORCID:0000000271092278)·Ren, Jun [University of Tennessee, Knoxville, TN (United States)] (ORCID:0000000222561007)·Ratynskaia, Svetlana [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0000000267123625)·Rizzi, Tommaso [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0009000521957260)·Paschalidis, Konstantinos [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0009000173335544)·Tolias, Panagiotis [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0000000196328104)·Kulachenko, Artem [KTH Royal Institute of Technology, Stockholm (Sweden)]·Pitts, R. A. [ITER Organization, Saint-Paul-lez-Durance (France)] (ORCID:0000000194552698)

Abstract

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

Explore related subjects

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

Hollmann, E. M. [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000262676589), Marini, Claudio [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000348197765), Rudakov, D. L. [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000252664269), Martinez-Loran, E. [University of California, San Diego, La Jolla, CA (United States)] (ORCID:0000000338414981), Beidler, Matthew [Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)] (ORCID:0000000273853886), Herfindal, Jeffrey L. [Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)] (ORCID:000000032846597X), Shiraki, Daisuke [Oak Ridge National Laboratory (ORNL), Oak Ridge, TN (United States)], Liu, Yueqiang Q. [General Atomics, San Diego, CA (United States)] (ORCID:0000000281928411), Bykov, I. [General Atomics, San Diego, CA (United States)], Lvovskiy, Andreas [General Atomics, San Diego, CA (United States)] (ORCID:0000000236491169), Lasnier, C. [Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)] (ORCID:0000000271092278), Ren, Jun [University of Tennessee, Knoxville, TN (United States)] (ORCID:0000000222561007), Ratynskaia, Svetlana [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0000000267123625), Rizzi, Tommaso [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0009000521957260), Paschalidis, Konstantinos [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0009000173335544), Tolias, Panagiotis [KTH Royal Institute of Technology, Stockholm (Sweden)] (ORCID:0000000196328104), Kulachenko, Artem [KTH Royal Institute of Technology, Stockholm (Sweden)], Pitts, R. A. [ITER Organization, Saint-Paul-lez-Durance (France)] (ORCID:0000000194552698). 2025-02-21. Measurement of post-disruption runaway electron kinetic energy and pitch angle during final loss instability in DIII-D. https://doi.org/10.1088/1361-6587%2Fadb5b6

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related reports

Fast solvers for tokamak fluid models with PETSc

Multigrid (MG) is widely recognized as a highly effective solver for the model problem, the Laplacian, but textbook MG fails on most problems of interest. MG methods have been applied to complex, real-world applications with careful consideration of the physical model and discretization. In this work we develop the first step in applying MG methods to science and engineering relevant magnetohydrodynamics (MHD) tokamak models in the M3D-C1 (https://m3dc1.pppl.gov) fusion energy science code. The semi-implicit time integrator in M3D-C1 is composed of many linear solves. The implicit advance of the momentum equation is the most challenging and is the focus of this work. The current production solver in M3D-C1 is a block Jacobi (BJ) preconditioner within a Krylov solver, where blocks group degrees of freedom on planes of constant toroidal coordinate. BJ convergence degrades as the number of planes increases due to the spectral properties of the matrix preconditioned with BJ. The partially magnetic field-aligned, regular toroidal grid structure in M3D-C1 is amenable to semi-coarsening geometric MG in the toroidal direction. This paper develops such a solver and demonstrates competitive performance on a runaway electron model of a SPARC (https://cfs.energy/technology/sparc) disruption, and superior robustness on a stellarator model on which the BJ solver fails to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Final Report: A Multi-Channel Fusion Product

The goal of this project was to measure charged fusion products from the d(d,p)t reaction in MAST-U plasmas as a function of time and position with good energy resolution using a system of up to six charged particle detectors. The data from this new diagnostic will make it possible to determine the neutral beam ion density profile as a function of R, z, and t with reduced model dependency and contribute new information to a global analysis of fast ion diagnostic data needed for the determination of the fast ion distribution function (velocity space tomography).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Unitary Qubit Lattice Algorithms for Plasma Physics

This final technical report summarizes research conducted under DOE Award DE-SC0021653 to develop unitary Quantum Lattice Algorithms for modeling electromagnetic wave propagation and scattering in complex media, including plasmas. The project developed and validated quantum-inspired formulations of Maxwell's equations that preserve unitary evolution and can be evaluated on classical high-performance computing systems while providing a foundation for future quantum-computing implementations. Major accomplishments include the development of two- and three-dimensional algorithms for electromagnetic scattering; scalable, distributed-memory implementations demonstrated on the Perlmutter supercomputer; formulations for nonlinear lossless fluid dynamics and cold, lossless, inhomogeneous magnetized plasmas; and an explicit quantum algorithm for a time-discretized Lorenz model. Simulations reproduced a range of characteristic wave phenomena, including transient effects that are not readily apparent in conventional frequency-domain studies, demonstrating the effectiveness of the proposed approach for modeling complex electromagnetic and plasma systems. The work establishes a unified theoretical and computational framework for quantum and quantum-inspired simulation and provides a foundation for future implementation on fault-tolerant quantum systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY