Abstract:In this work we consider a relativistic drift-kinetic model for runaway electrons along with a Fokker-Planck operator for small-angle Coulomb collisions, a radiation damping operator, and a secondary knock-on (Boltzmann) collision source. We develop a new scalable fully implicit solver utilizing finite volume and conservative finite difference schemes and dynamic mesh adaptivity. A new data management framework in the PETSc library based on the p4est library is developed to enable simulations with dynamic adaptive mesh refinement (AMR), distributed memory parallelization, and dynamic load balancing of computational work. This framework and the runaway electron solver building on the framework are able to dynamically capture both bulk Maxwellian at the low-energy region and a runaway tail at the high-energy region. To effectively capture features via the AMR algorithm, a new AMR indicator prediction strategy is proposed that is performed alongside the implicit time evolution of the solution. This strategy is complemented by the introduction of computationally cheap feature-based AMR indicators that are analyzed theoretically. Numerical results quantify the advantages of the prediction strategy in better capturing features compared with nonpredictive strategies; and we demonstrate trade-offs regarding computational costs. The robustness with respect to model parameters, algorithmic scalability, and parallel scalability are demonstrated through several benchmark problems including manufactured solutions and solutions of different physics models. We focus on demonstrating the advantages of using implicit time stepping and AMR for runaway electron simulations.

PETSc/TS: A Modern Scalable ODE/DAE Solver Library

PETSc/TAO Developments for GPU-Based Early Exascale Systems

Direct Approaches for PDE-constrained Dynamic Optimization Based on Space–time Orthogonal Collocation on Finite Elements

A Quasi-Sequential Algorithm for PDE-constrained Optimization Based on Space–time Orthogonal Collocation on Finite Elements

Efficiently and easily integrating differential equations with JiTCODE, JiTCDDE, and JiTCSDE

Application of PETSc in Soil-Water Coupled Geotechnical Problems

Flexible, Scalable Mesh and Data Management using PETSc DMPlex

AFEPack: A General-Purpose C++ Library for Numerical Solutions of Partial Differential Equations

Open source Differentiable ODE Solving Infrastructure

A C++ implementation of the discrete adjoint sensitivity analysis method for explicit adaptive Runge-Kutta methods enabled by automatic adjoint differentiation and SIMD vectorization

SUNDIALS time integrators for exascale applications with many independent systems of ordinary differential equations

An Efficient Local Optimizer-Tracking Solver for Differential-Algebriac Equations with Optimization Criteria

Extending a serial 3D two-phase CFD code to parallel execution over MPI by using the PETSc library for domain decomposition

Scalability of high-performance PDE solvers

SUNDIALS Time Integrators for Exascale Applications with Many Independent ODE Systems

A multi-physics compiler for generating numerical solvers from differential equations

A robust solver for large-scale heat transfer topology optimization

PETScML: Second-order solvers for training regression problems in Scientific Machine Learning

Scalable Implicit Solvers with Dynamic Mesh Adaptation for a Relativistic Drift-Kinetic Fokker-Planck-Boltzmann Model

Extreme-scale Multigrid Components within PETSc

Automated MPI-X code generation for scalable finite-difference solvers