Spatially partitioned embedded Runge-Kutta methods

From MaRDI portal



Abstract: We study spatially partitioned embedded Runge--Kutta (SPERK) schemes for partial differential equations (PDEs), in which each of the component schemes is applied over a different part of the spatial domain. Such methods may be convenient for problems in which the smoothness of the solution or the magnitudes of the PDE coefficients vary strongly in space. We focus on embedded partitioned methods as they offer greater efficiency and avoid the order reduction that may occur in non-embedded schemes. We demonstrate that the lack of conservation in partitioned schemes can lead to non-physical effects and propose conservative additive schemes based on partitioning the fluxes rather than the ordinary differential equations. A variety of SPERK schemes are presented, including an embedded pair suitable for the time evolution of fifth-order weighted non-oscillatory (WENO) spatial discretizations. Numerical experiments are provided to support the theory.


The paper deals with spatially partitioned embedded Runge-Kutta (SPERK) methods as an appropriate alternative to more traditional schemes in the time integration of the ordinary differential equation (ODE) resulting from a previous space discretization of partial differential equations (PDEs) when the smoothness of their solution or the magnitude of their coefficients vary strongly in space. The attention is focused on embedded pairs, since they offer a better efficiency and avoid the order reduction which is present in other non-embedded schemes.NEWLINENEWLINEWith SPERK integrators, each component method of the overall scheme is applied over different parts of the spatial domain, and this allows, in particular, to use larger time steps. SPERK methods are applied in two different ways. In the first (the equation-based partitioning), the ODE is partitioned, and at each point in space a particular time-stepping scheme is selected. In this case, the order of accuracy is shown to be equal to the minimum of the order of accuracy of the schemes composing the embedded pair, but on the other hand it is not conservative when applied to a scalar conservation law.NEWLINENEWLINEThe second partitioning method operates on fluxes, and is called flux-based partitioning. This has the advantage of being conservative when applied to a conservation law. It is shown that it may suffer from order reduction, although this loss of accuracy is not observed in practice.NEWLINENEWLINEBoth procedures are positivity preserving under a suitable time step restriction when the component schemes are strong stability preserving methods.NEWLINENEWLINESeveral explicit SPERK schemes are designed and numerically illustrated on numerical examples. One of them, in particular, is suitable to be used in the time integration of PDEs, previously discretized in space with fifth-order weighted non-oscillatory techniques.




Cited in
(19)








This page was built for publication: Spatially partitioned embedded Runge-Kutta methods

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2870632)