Since 18 of December 2019 uses Nucleus credentials. Visit our help pages for information on how to Register and Sign-in using Nucleus.
10-15 May 2021
Nice, France
Europe/Vienna timezone
The Conference will be held virtually from 10-15 May 2021

Verification of an Adaptive Sparse Grid Approach to Simulation of Runaway Electron Dynamics

11 May 2021, 08:30
Nice, France

Nice, France

Regular Poster Magnetic Fusion Theory and Modelling P1 Posters 1


David Green (Oak Ridge National Laboratory)


We demonstrate the reduced computational cost of applying adaptive sparse grids to the problem of linear kinetic simulation of runaway electron dynamics. This work is motivated by the need to routinely perform high-dimensional kinetic simulation in a noise-free manner, with the implication being that a reduced computational cost for kinetic simulation allows for more routine investigation. Here we verify our solver with a series of tests (e.g., Fig. 1), including demonstrating that it yields the expected runaway electron production rate. We show that a non-adapted sparse grid approach yields significant savings in the computational resources required to resolve the runaway electron dynamics to sufficient accuracy. We also address a common criticism of sparse grids, in that we show that adding adaptivity to sparse grids recovers the advantage of sparse grids when the features of the solution are not coordinate aligned.
(left) 2D (0x-2v) Maxwellian electron distribution function and sparse grid locations (white crosses) for a lev=5 hierarchical, deg=5 basis. **(right):** Time dynamics of a cos(pitch angle)=0 slice of fe relaxing to a Maxwellian (dotted) showing accuracy to at least 10-12.

Our linear model includes terms for electric field acceleration, Coulomb collisions, and radiation damping. Figure 1 shows a test where we simulate on a sparse grid for a sufficiently small accelerating electric field that the distribution relaxes to a known Maxwellian (dotted line) in steady-state. There are several other efforts focused on the phase space simulation of runaway electron dynamics (e.g., 1, 2, [3], [4]). However, of these, it is only the Monte-Carlo approaches that have begun looking at including the full spatial variation required for quantitative analysis and prediction. This is because the Monte-Carlo methods require less computational resources than the continuum (mesh/grid-based) methods. As these approaches are advanced to include a self-consistent coupling of the distribution function evolution and field solutions to capture how the runaway population feeds back on the accelerating electric field, a continuum approach that can incorporate the spatial dimensions at reasonable compute cost is motivated. Here we investigate if adaptive sparse grids [5,6] allow for calculating the dynamics of the phase space distribution of electrons in a robust and accurate manner. Our approach is constructing a high-order, discontinuous Galerkin solver using a "sparse" polynomial space where the sparsity is introduced via a truncated tensor product of 1D bases. We employ upwinding for advection terms and LDG [7] for diffusion operators.
We further demonstrate (cf., Fig. 2) that adding adaptivity to the sparse grid discretization recovers the advantages of sparse grids in such where it is traditionally lost. The truncated tensor product drops, for example, those elements from the space whose total degree is less than the maximum degree in any one dimension. This will lead to cases where the advantages of the sparse grid are lost when a solution relies on the parts of the tensor product which were thrown away. Figure 2 demonstrates how adaptivity makes sparse grids robust to this issue.

Demonstration of adaptive sparse grids retaining a reduction in degrees of freedom for solutions with non-coordinate aligned features. Panels (a)-(c) represent manufactured solutions with a parameter C whose value represents the degree of coordinate misalignment. Panel (a) also shows the full-grid (FG) locations, panel (b) the sparse grid (SG) locations, and panel (c) the adaptive-sparse-grid (ASG) locations. Panels (d)-(f) show the error analysis for the 3 grid types for a deg=3 (d/e) and deg=5 (f) basis.

In summary, the adaptive sparse grid approach appears to be a viable approach to enabling more routine kinetic simulation of fusion plasmas. We will continue to develop and apply this approach to incorporate more of the spatial dimensionality typically lacking from continuum simulation of runaway electron dynamics.
This research used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-00OR22725.

1 D. del-Castillo-Negrete et al., Physics of Plasmas 25.5 (2018): 056104.
2 M. Landreman, et al., Comp. Physics Comm. 185.3 (2014): 847-855.
[3] R. W. Harvey et al., Physics of Plasmas 7.11 (2000): 4590-4599.
[4] D. Daniel et al., arXiv preprint arXiv:1902.10241 (2019).
[5] D. Pflüger et al., Journal of Complexity 26.5 (2010): 508-522.
[6] W. Guo et al., SIAM Journal on Scientific Computing 38.6 (2016): A3381-A3409.
[7] B. Cockburn et al., SIAM Journal on Numerical Analysis 35.6 (1998): 2440-2463.

Affiliation Oak Ridge National Laboratory
Country or International Organization United States

Primary authors

David Green (Oak Ridge National Laboratory) Dr Lin Mu (University of Georgia) Dr Eduardo D'Azevedo (Oak Ridge National Laboratory) Diego del-Castillo-Negrete (Fusion Energy Division. Oak Ridge National Laboratory) Dr Matthew Graham Lopez (Oak Ridge National Laboratory) Mr Benjamin Tyler McDaniel (Oak Ridge National Laboratory) Wael Elwasif (Oak Ridge National Laboratory) Timothy R. Younkin (University of Tennessee ) Dr Louis Wonnell (Oak Ridge National Laboratory) Mr Harry Hughes (Oak Ridge National Laboratory)

Presentation Materials

There are no materials yet.