Abhijit Biswas,David I. Ketcheson,Steven Roberts,Benjamin Seibold,David Shirokoff
Abstract:Explicit Runge--Kutta (RK) methods are susceptible to a reduction in the observed order of convergence when applied to initial-boundary value problem with time-dependent boundary conditions. We study conditions on explicit RK methods that guarantee high-order convergence for linear problems; we refer to these conditions as weak stage order conditions. We prove a general relationship between the method's order, weak stage order, and number of stages. We derive explicit RK methods with high weak stage order and demonstrate, through numerical tests, that they avoid the order reduction phenomenon up to any order for linear problems and up to order three for nonlinear problems.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of order reduction that occurs when explicit Runge - Kutta (RK) methods are used to handle initial - boundary value problems (IBVP) with time - dependent boundary conditions. Specifically:
1. **Order reduction phenomenon**: When explicit RK methods are applied to stiff problems or initial - boundary value problems with time - dependent boundary conditions, the actually observed convergence order is usually lower than the order predicted by classical analysis. This is mainly due to the low "stage order" of the method.
2. **Linear and nonlinear problems**: For many problems, such as hyperbolic partial differential equations (PDE), explicit RK methods are the first choice, but since the stage order of explicit RK methods is at most 1, order reduction becomes an important challenge.
3. **Limitations of existing solutions**: Existing methods for avoiding order reduction usually involve modifying the boundary conditions, but these methods have many disadvantages, such as requiring derivatives of the boundary conditions, reducing the stable time step, introducing additional numerical parameters, and making complex modifications to the solver code.
4. **Proposed new method**: This paper proposes a method to avoid order reduction by selecting RK coefficients that satisfy specific weak stage order conditions. This method does not require complex modifications to the boundary conditions, but is achieved by optimizing the coefficients of the RK method.
5. **Main contributions**:
- Provide theoretical results that enable the construction of explicit RK methods with high weak stage order.
- Characterize the linear and nonlinear stability properties of these methods.
- Construct explicit RK methods with arbitrary high order and weak stage order, and demonstrate their effectiveness and superiority in a series of problems.
### Markdown representation of formulas
To ensure the correctness and readability of the formulas, the following is the Markdown - format representation of the key formulas in the paper:
- **Definition of weak stage order**:
\[
Y := \text{span}\left\{ b, A^Tb, \ldots, (A^T)^{s - 1}b \right\}
\]
\[
K_m := \left\{ \tau(1), A\tau(1), \ldots, A^{s - 1}\tau(1), \tau(2), A\tau(2), \ldots, A^{s - 1}\tau(m) \right\}
\]
\[
\tau(k) := Ac_k-\frac{1}{k}c_k
\]
- **Weak stage order conditions**:
\[
0 = b^T(I - zA)^{-1}\tau(k),\quad \forall z\in\mathbb{C},\quad 1\leq k\leq q
\]
- **Stability function**:
\[
R(z) := 1+\sum_{j = 1}^s z^j b^T A^{j - 1}e
\]
Through these formulas, the paper explores in detail how to construct efficient explicit RK methods by satisfying the weak stage order conditions, thereby avoiding the order reduction problem.