Dipartimento di Matematica Guido Castelnuovo, Università Sapienza Roma
Yulong Xing (Ohio State University)
In the common practice of the method-of-lines (MOL) approach for discretizing a time-dependent partial differential equation (PDE), one first applies spatial discretization to convert the PDE into an ordinary differential equation system. Subsequently, a time integrator is used to discretize the time variable. When a multi-stage Runge-Kutta (RK) method is used for time integration, by default, the same spatial operator is used at all RK stages. However, recent studies on perturbed RK methods indicate that not all RK stages are born equal – breaking the MOL structure and applying rough approximations at specific RK stages may not affect the overall accuracy of the numerical scheme. In this talk, we present two of our recent explorations on blending rough stage operators in RK discontinuous Galerkin (DG) methods for solving hyperbolic conservation laws. In our first work, we mix the DG operator with the local derivative operator, yielding an RKDG method featuring compact stencils and simple boundary treatment. In our second work, we mix the DG operators with polynomials of degrees k and k-1, and the resulting method may allow larger time step sizes and fewer floating-point operations per time step.
giuseppe.visconti@uniroma1.it