Multiderivative Runge-Kutta Flux Reconstruction for Hyperbolic Conservation Laws

Abstract We extend the fourth-order, two-stage multiderivative Runge-Kutta (MDRK) scheme to the flux reconstruction (FR) framework by writing both stages in terms of a time-averaged flux and then using the approximate Lax-Wendroff (LW) procedure to compute the time-averaged flux. Numerical flux is carefully constructed to enhance Fourier CFL stability and accuracy. A subcell-based blending limiter is developed for the MDRK scheme which ensures that the limited scheme is provably admissibility preserving. Along with being admissibility preserving, the blending scheme is constructed to minimize dissipation errors using Gauss-Legendre (GL) solution points and performing the MUSCL-Hancock (MH) reconstruction on subcells. The accuracy enhancement of the blending scheme is numerically verified on compressible Euler equations, with test cases involving shocks and small-scale structures.