Explicit Stabilized Multirate Method For Stiff Differential Equations

Stabilized Runge???Kutta methods are especially efficient for the numerical solution of large systems of stiff nonlinear differential equations because they are fully explicit. For semi-discrete parabolic problems, for instance, stabilized Runge???Kutta methods overcome the stringent stability condition of standard methods without sacrificing explicitness. However, when stiffness is only induced by a few components, as in the presence of spatially local mesh refinement, their efficiency deteriorates. To remove the crippling effect of a few severely stiff components on the entire system of differential equations, we derive a modified equation, whose stiffness solely depends on the remaining mildly stiff components. By applying stabilized Runge???Kutta methods to this modified equation, we then devise an explicit multirate Runge???Kutta??? Chebyshev (mRKC) method whose stability conditions are independent of a few severely stiff components. Stability of the mRKC method is proved for a model problem, whereas its efficiency and usefulness are demonstrated through a series of numerical experiments.

MATHICSE technical Report : Stabilized explicit multirate methods for stiff differential equations

Assyr Abdulle, Giacomo Rosilho De Souza

Stabilized Runge–Kutta (aka Chebyshev) methods are especially efficient for the numerical solution of large systems of stiff differential equations because they are fully explicit; hence, they are inherently parallel and easily accommodate nonlinearity. For semi-discrete parabolic (or diffusion dominated) problems, for instance, stabilized Runge–Kutta methods overcome the stringent stability condition of standard methods without sacrificing explicitness. However, when much of the stiffness is only induced by a few components, as in the presence of spatially local mesh refinement, their efficiency deteriorates. To remove the crippling effect of a few severely stiff components on the entire system of differential equations, we derive a modified equation, whose stiffness solely depend on the remaining mildly stiff components. By applying stabilized Runge–Kutta methods to this modified equation, we then devise an explicit multirate Runge–Kutta–Chebyshev (mRKC) method whose stability conditions are independent of a few severely stiff components. Stability of the mRKC method is proved for a model problem, whereas its efficiency and usefulness are demonstrated through a series of numerical experiments.

MATHICSE2020

Explicit Stabilized Multirate Method For Stiff Stochastic Differential Equations

Assyr Abdulle, Giacomo Rosilho De Souza

Stabilized explicit methods are particularly efficient, for large systems of stiff stochastic differential equations (SDEs) due to their extended stability domain. However, they lose their efficiency when a severe stiffness is induced by very few "fast" degrees of freedom, as the stiff and nonstiff terms are evaluated concurrently. Therefore, inspired by [A. Abdulle, M. J. Grote, and G. Rosilho de Souza, Explicit stabilized multirate method for stiff differential equations, Math. Comp., in press, 2022], we introduce a stochastic modified equation whose stiffness depends solely on vi the "slow" terms. By integrating this modified equation with a stabilized explicit scheme, we devise a multirate method which overcomes the bottleneck caused by a few severely still terms and recovers the efficiency of stabilized schemes for large systems of nonlinear SDEs. The scheme is not based on any scale separation assumption of the SDE. Therefore, it is employable for problems stemming from the spatial discretization of stochastic parabolic partial differential equations on locally refined grids. The multirate scheme has strong Order 1/2, weak order 1, and its stability is proved on a model problem. Numerical experiments confirm the efficiency and accuracy of the scheme.

SIAM PUBLICATIONS2022

A posteriori error estimation for partial differential equations with random input data

Diane Sylvie Guignard

This thesis is devoted to the derivation of error estimates for partial differential equations with random input data, with a focus on a posteriori error estimates which are the basis for adaptive strategies. Such procedures aim at obtaining an approximation of the solution with a given precision while minimizing the computational costs. If several sources of error come into play, it is then necessary to balance them to avoid unnecessary work. We are first interested in problems that contain small uncertainties approximated by finite elements. The use of perturbation techniques is appropriate in this setting since only few terms in the power series expansion of the exact random solution with respect to a parameter characterizing the amount of randomness in the problem are required to obtain an accurate approximation. The goal is then to perform an error analysis for the finite element approximation of the expansion up to a certain order. First, an elliptic model problem with random diffusion coefficient with affine dependence on a vector of independent random variables is studied. We give both a priori and a posteriori error estimates for the first term in the expansion for various norms of the error. The results are then extended to higher order approximations and to other sources of uncertainty, such as boundary conditions or forcing term. Next, the analysis of nonlinear problems in random domains is proposed, considering the one-dimensional viscous Burgers' equation and the more involved incompressible steady-state Navier-Stokes equations. The domain mapping method is used to transform the equations in random domains into equations in a fixed reference domain with random coefficients. We give conditions on the mapping and the input data under which we can prove the well-posedness of the problems and give a posteriori error estimates for the finite element approximation of the first term in the expansion. Finally, we consider the heat equation with random Robin boundary conditions. For this parabolic problem, the time discretization brings an additional source of error that is accounted for in the error analysis. The second part of this work consists in the analysis of a random elliptic diffusion problem that is approximated in the physical space by the finite element method and in the stochastic space by the stochastic collocation method on a sparse grid. Considering a random diffusion coefficient with affine dependence on a vector of independent random variables, we derive a residual-based a posteriori error estimate that controls the two sources of error. The stochastic error estimator is then used to drive an adaptive sparse grid algorithm which aims at alleviating the so-called curse of dimensionality inherent to tensor grids. Several numerical examples are given to illustrate the performance of the adaptive procedure.

EPFL2016

A posteriori error estimation for partial differential equations with random input data

Diane Sylvie Guignard

EPFL2016