Rujeko Chinomona

Rujeko Chinomona, Dibyendu Adak, William J. Barham et al.

We present a tensor-train discontinuous Galerkin (TT-DG) formulation for the Vlasov--Maxwell system that combines a modal DG discretization with low-rank tensor representations of the phase-space solution and discrete operators. The formulation exploits the tensor-product structure of the DG discretization to perform quadrature, differentiation, nonlinear upwind flux evaluation, and time integration directly in compressed form. The method is evaluated on several standard 1D2V Vlasov--Maxwell benchmark problems, including the streaming Weibel instability, weak Landau damping, and two-stream instability problems. Across these problems, the TT formulation reproduces the accuracy and conservation behavior of the underlying full-grid DG discretization while substantially reducing memory usage and runtime. For weakly nonlinear problems, compression ratios exceeding $10^4$ are obtained together with significant speedups relative to the full-grid solver. For the strongly nonlinear two-stream instability problem, the TT formulation remains effective despite reduced compressibility caused by fine-scale phase-space filamentation. These results demonstrate that tensor-train representations provide an effective approach for reducing the computational cost of deterministic DG-based kinetic plasma simulations while retaining the favorable numerical properties of the underlying discretization.

Vu Thai Luan, Rujeko Chinomona, Daniel R. Reynolds

This work focuses on the development of a new class of high-order accurate methods for multirate time integration of systems of ordinary differential equations. The proposed methods are based on a specific subset of explicit one-step exponential integrators. More precisely, starting from an explicit exponential Runge--Kutta method of the appropriate form, we derive a multirate algorithm to approximate the action of the matrix exponential through the definition of modified "fast" initial-value problems. These fast problems may be solved using any viable solver, enabling multirate simulations through use of a subcycled method. Due to this structure, we name these Multirate Exponential Runge--Kutta (MERK) methods. In addition to showing how MERK methods may be derived, we provide rigorous convergence analysis, showing that for an overall method of order $p$, the fast problems corresponding to internal stages may be solved using a method of order $p-1$, while the final fast problem corresponding to the time-evolved solution must use a method of order $p$. Numerical simulations are then provided to demonstrate the convergence and efficiency of MERK methods with orders three through five on a series of multirate test problems.