Abstract:The aim of this paper is to design the explicit radial basis function (RBF) Runge-Kutta methods for the initial value problem. We construct the two-, three- and four-stage RBF Runge-Kutta methods based on the Gaussian RBF Euler method with the shape parameter, where the analysis of the local truncation error shows that the s-stage RBF Runge-Kutta method could formally achieve order s+1. The proof for the convergence of those RBF Runge-Kutta methods follows. We then plot the stability region of each RBF Runge-Kutta method proposed and compare with the one of the correspondent Runge-Kutta method. Numerical experiments are provided to exhibit the improved behavior of the RBF Runge-Kutta methods over the standard ones.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to design explicit radial basis function (RBF) Runge - Kutta methods for solving initial - value problems. Specifically, the author aims to construct two - order, three - order and four - order RBF Runge - Kutta methods based on the Gaussian RBF Euler method, and shows through the analysis of local truncation errors that the s - order RBF Runge - Kutta method can theoretically achieve s + 1 - order accuracy. In addition, the paper also explores the convergence of these RBF Runge - Kutta methods, and plots the stability regions of each method, comparing them with the corresponding standard Runge - Kutta methods. The improved performance of the RBF Runge - Kutta methods in practical applications is demonstrated through numerical experiments.
### Background of the Paper
The Runge - Kutta method is a widely used method for solving initial - value problems, in the following form:
\[
u' = f(t, u), \quad a < t \leqslant b
\]
\[
u(a) = u_0
\]
Traditional s - order explicit Runge - Kutta methods can improve accuracy by calculating multiple intermediate values within a single step. However, for s ≥ 5, s - order accuracy cannot be achieved using only s stages. To overcome this limitation, researchers have proposed various improvement methods, such as Total Variation Diminishing (TVD) Runge - Kutta methods, Strong Stability Preserving (SSP) Runge - Kutta methods, exponential Runge - Kutta methods, etc.
### Main Contributions of the Paper
1. **Design of RBF Runge - Kutta Methods**: The paper proposes a new explicit RBF Runge - Kutta method, which improves the order of the method by introducing shape parameters to eliminate the main error terms.
2. **Theoretical Analysis**: Through local truncation error analysis, it is proved that the s - order RBF Runge - Kutta method can achieve s + 1 - order accuracy.
3. **Convergence Proof**: The convergence of the proposed RBF Runge - Kutta method is proved in detail.
4. **Stability Analysis**: The stability regions of each RBF Runge - Kutta method are plotted and compared with the standard Runge - Kutta methods.
5. **Numerical Experiments**: The superior performance of the RBF Runge - Kutta method in practical applications, especially in terms of accuracy improvement, is verified through numerical experiments.
### Specific Methods
- **Two - order RBF Runge - Kutta Method**:
\[
v_{n + 1}=v_n + h(b_1 k_1 + b_2 k_2)
\]
\[
k_1 = f(t_n, v_n)
\]
\[
k_2 = f(t_n + c_2 h, v_n e^{-\epsilon^2_n (a_{21} h)^2}+h a_{21} k_1)
\]
- **Three - order RBF Runge - Kutta Method**:
\[
v_{n + 1}=v_n + h(b_1 k_1 + b_2 k_2 + b_3 k_3)
\]
\[
k_1 = f(t_n, v_n)
\]
\[
k_2 = f(t_n + c_2 h, v_n e^{-\epsilon^2_{n2} (c_2 h)^2}+h a_{21} k_1)
\]
\[
k_3 = f(t_n + c_3 h, v_n e^{-\epsilon^2_{n3} (c_3 h)^2}+h (a_{31} k_1 + a_{32} k_2))
\]
- **Four - order RBF Runge - Kutta Method**:
\[
v_{n + 1}=v_n + h(b_1 k_1 + b_2 k_2 + b_3 k_3 + b_4 k_4)
\]
\[
k_1 = f(t_n, v_n)
\]
\[
k_2 = f(t_n + c_2 h, v_n e^{-\epsilon^2_{n2} (c_2 h)^2}+h a_{21} k_1