Abstract:Higher-order numerical methods are used to find accurate numerical solutions to hyperbolic partial differential equations and equations of transport type. Limiting is required to either converge to the correct type of solution or to adhere to physically motivated local maximum principles. Less restrictive limiting procedures are required so as to not severely decrease the accuracy.
In this paper, we develop an existing slope limiter framework, to achieve different local boundedness principles for higher-order schemes on unstructured meshes. Quadrature points contributing to numerical fluxes can be limited based on face defined maximum principles, and the resulting cell mean at the next timestep can satisfy a cell mean maximum principle but with less limiting. We demonstrate the practical application of the introduced framework to a second-order finite volume scheme as well as a fourth-order finite volume scheme, in the context of the advection equation.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to maintain the Local Maximum Principles (LMPs) in high - order numerical methods while reducing the impact of the limiting process on accuracy. Specifically, the author has developed an existing slope limiter framework to achieve different types of local boundedness principles on unstructured grids, thereby ensuring that the numerical solution satisfies the physically reasonable local maximum principle while guaranteeing high - order accuracy. This includes:
1. **Limiting problems in high - order numerical methods**: Limiting is required to ensure that the numerical solution converges to the correct type or follows the physically motivated local maximum principle. However, traditional limiting methods may significantly reduce the accuracy of the solution.
2. **Maintenance of the local maximum principle**: The paper proposes a new slope limiter framework that can maintain the local maximum principle without severely reducing accuracy. In particular, by limiting the maximum principle defined for faces, the cell - average value in the next time step satisfies the maximum principle for cell - average values, but with less limitation.
3. **Application examples**: The paper shows the practical application of this framework in the second - order and fourth - order finite - volume methods for the two - dimensional advection equation.
### Specific problem description
- **Background**: High - order numerical methods are widely used to solve hyperbolic partial differential equations and transport equations. These methods can provide more accurate numerical solutions, but in order to ensure the physical rationality of the solutions, limiting is usually required.
- **Problem**: Traditional limiting methods are often too strict, resulting in a decline in the accuracy of the solution. Therefore, a new limiting framework needs to be developed that can maintain the local maximum principle while reducing the impact on the accuracy of the solution.
### Main contributions of the paper
1. **New slope limiter framework**: The author has extended the existing slope limiter framework and proposed new limiting conditions, enabling high - order finite - volume methods to maintain the local maximum principle on unstructured grids.
2. **Theoretical analysis**: Through Theorem 2.1, the author gives sufficient conditions for methods of any order to maintain the local boundedness principle on unstructured grids, and introduces some symbols and diagrams to illustrate these conditions.
3. **Practical application**: The paper shows the application of this framework in the second - order and fourth - order finite - volume methods, especially in the two - dimensional transport equation. Through these examples, the effectiveness and practicality of the new framework are verified.
### Mathematical formulas
- **Local maximum principle**:
\[
\bar{u}_K^{n + 1} \in [m_K^n, M_K^n], \quad \forall K \in \mathcal{M}
\]
where \(m_K^n\) and \(M_K^n\) are the time - dependent maximum and minimum values defined locally at time step \(n\).
- **Evolution equation of high - order finite - volume methods**:
\[
\bar{u}_K^{n + 1} = \bar{u}_K^n - \frac{\Delta t}{|K|} \sum_{L \in N(K)} |\sigma_{KL}| \sum_{q \in \sigma_{KL}} w_{\sigma_{KL}}^q f_{KL}(p_K(x_q), p_L(x_q), \mathbf{v} \cdot \mathbf{n}_{KL})
\]
- **Conditions for the local maximum principle**:
\[
\Delta t \frac{|\sigma_{KL}| w_{\sigma_{KL}}^q}{w_K^q |K|} \frac{\partial f_{KL}}{\partial p_K} \leq 1, \quad \forall q \in \sigma_{KL}, \forall L \in N(K), \forall K \in \mathcal{M}
\]
Through these formulas and theoretical analysis, the paper successfully solves the problem of maintaining the local maximum principle in high - order numerical methods while reducing the impact of limiting on accuracy.