Skip to content

Author

Stefan R. Schnake

2 papers indexed here

We haven’t gathered this author’s papers yet. Follow them and we’ll fetch their work.

Not the right person? Other researchers publish under this name.

Preprint Sep 2026

Adaptive Sparse-grid Discontinuous Galerkin Approximations the Bhatnagar--Gross--Krook Model

This work studies adaptive sparse-grid discontinuous Galerkin (DG) discretizations for the Bhatnagar--Gross--Krook (BGK) model, a kinetic equation posed in four- and six-dimensional phase-space. Standard DG methods are rendered impractical for the BGK model by the curse of dimensionality, motivating compressed representations that adapt to the solution in time. Using the adaptive sparse-grid DG method, we quantify accuracy and compression by comparing the adaptive degrees of freedom to full-grid DG methods and by assessing the resulting kinetic and fluid quantities in both fluid and rarefied regimes. Test cases include a relaxation problem, a multidimensional Sod shock tube, and shear/expansion flows used in prior low-rank BGK studies. To build an efficient Maxwellian evaluation without violating conservation, a central obstacle for structure-perserving BGK simulations, we introduce a hybrid interpolation strategy that exploits velocity separability to recover the correct discrete collision invariants and prove conservation of the resulting discrete collision operator on adaptive sparse grids. Our results show that the adaptive sparse-grid strategy can recover accurate and physically relevant solutions with sharp gradients, and the method reduces the active degrees of freedom by factors ranging from several-fold to several orders of magnitude, with the largest reductions occurring in the six-dimensional examples. All computations are performed with the open-source ASGarD adaptive sparse-grid DG library.

Stefan R. Schnake, M. Stoyanov, E. Endeve et al. · 0 citations
Preprint Aug 2026

Sweep-based, implicit solutions of the multidimensional BGK equation on unstructured grids

We present a nodal discontinuous Galerkin method for solving the Bhatnagar-Gross-Krook (BGK) kinetic equation on multi-dimensional, unstructured grids. The method uses implicit, sweep-based solvers and a moment-preserving projection of the Maxwellian source to enable high-order accuracy in time while avoiding restrictive time steps imposed by boundary layers and other geometry-induced features. We verify that the method is correct in the continuum limit by comparing to closed-form and high-order solutions of the Sod shock problem on 2 and 3D unstructured grids. Linear L2 stability is demonstrated for a B-stable diagonally implicit Runge-Kutta method of third order. The solver uses a hybrid parallel scheme based on spatial domain decomposition with local sweeps performed on CPU and GPU hardware. Platform-portability is demonstrated through the development of new GPU-friendly, graph-based sweep algorithms that are implemented using the Kokkos performance portability library and achieve greater than 20 times speedup on NVIDIA H100 GPUs compared to 64-core AMD EPYC 9654 CPUs. Finally, we show results on the Frontier supercomputer at the Oak Ridge Leadership Computing Facility for a boundary value problem with 2.77 trillion phase space degrees of freedom that executed on 1536 nodes utilizing 6144 AMD MI250X GPUs.

T. Evans, Ryan S. Glasby, Cory D. Hauck et al. · 0 citations

We use cookies to run the site and, with your consent, for analytics and to show ads. See our Cookie Policy.