Brian K. Tran,Ben S. Southworth,Tommaso Buvoli
Abstract:Recently a new class of nonlinearly partitioned Runge-Kutta (NPRK) methods was proposed for nonlinearly partitioned systems of ordinary differential equations, $y' = F(y,y)$. The target class of problems are ones in which different scales, stiffnesses, or physics are coupled in a nonlinear way, wherein the desired partition cannot be written in a classical additive or component-wise fashion. Here we use rooted-tree analysis to derive full order conditions for NPRK$_M$ methods, where $M$ denotes the number of nonlinear partitions. Due to the nonlinear coupling and thereby mixed product differentials, it turns out the standard node-colored rooted-tree analysis used in analyzing ODE integrators does not naturally apply. Instead we develop a new edge-colored rooted-tree framework to address the nonlinear coupling. The resulting order conditions are enumerated, provided directly for up to 4th order with $M=2$ and 3rd-order with $M=3$, and related to existing order conditions of additive and partitioned RK methods.
What problem does this paper attempt to address?
This paper attempts to solve the order condition problems faced by the Non - linear Partitioned Runge - Kutta (NPRK) method when solving coupled systems with different scales, stiffness, and physical properties. Specifically, the paper mainly addresses the following issues:
1. **Limitations of Existing Methods**:
- Traditional methods such as the Additive Runge - Kutta (ARK) method are suitable for linear or component - decomposed systems, but for non - linearly coupled systems (i.e., systems that cannot be separated in the classical additive or component way), these methods are no longer applicable.
- The Non - linear Partitioned Runge - Kutta (NPRK) method can handle each parameter with different degrees of implicitness and can handle complex multi - physics problems more flexibly.
2. **New Theoretical Framework**:
- The paper introduces a new framework based on edge - colored rooted trees to derive the order conditions of the NPRK method. Due to the existence of non - linear coupling, the traditional node - colored rooted tree analysis is not applicable, so new tools need to be developed.
- Through this method, the paper is able to derive order conditions up to the fourth order (for two non - linear partitions) and third - order conditions (for three non - linear partitions).
3. **Calculation of Order Conditions**:
- Specific algorithms are provided to calculate the order conditions of the NPRK method, and the specific results are listed in the appendix.
- The order conditions of the NPRK method and the ARK method are compared, revealing the additional order conditions in the NPRK method caused by non - linear coupling.
4. **Application Prospects**:
- The NPRK method has shown good results in dealing with problems such as thermal radiation transfer and radiation hydrodynamics, and has the potential to be applied to other multi - physics systems.
- By studying these order conditions, a better understanding of the stability and accuracy of the NPRK method can be achieved, thus providing theoretical support for practical applications.
In summary, this paper aims to fill the gap in the theory of order conditions for the Non - linear Partitioned Runge - Kutta method and provides new theoretical tools and methods for solving the numerical integration problems of complex coupled systems.
### Formula Display
To ensure the correctness and readability of the formulas, the following are some key formulas mentioned in the paper:
- **General Form of the NPRK Method**:
\[
Y_i = y_n+h\sum_{j = 1}^s\sum_{k = 1}^s a_{ijk}F(Y_j, Y_k),\quad i = 1,\ldots,s,
\]
\[
y_{n + 1}=y_n+h\sum_{i = 1}^s\sum_{j = 1}^s b_{ij}F(Y_i, Y_j).
\]
- **Derivatives in Taylor Expansion**:
\[
y'=F(y, y),
\]
\[
y'' = D_1F(y, y)y'+D_2F(y, y)y',
\]
\[
y'''=D_{11}F(y, y)[y', y']+2D_{12}F(y, y)[y', y']+D_{22}F(y, y)[y', y']+D_1F(y, y)D_1F(y, y)y'+D_1F(y, y)D_2F(y, y)y'+D_2F(y, y)D_1F(y, y)y'+D_2F(y, y)D_2F(y, y)y'.
\]
- **Order Conditions**:
According to Theorem 3.1, the order conditions of the NPRK method are:
\[
\Phi(\tau)=\frac{1}{\gamma(\tau)}
\]
where $\tau$ is an edge - colored rooted tree, $\gamma(\tau)$ is the density of the tree, and $\Phi(\tau)$ is the elementary weight of the tree.
These formulas and theoretical frameworks together form the core content of this paper, helping to understand and derive the order conditions of the NPRK method.