Abstract:We introduce a class of unconditionally energy stable, high order accurate schemes for gradient flows in a very general setting. The new schemes are a high order analogue of the minimizing movements approach for generating a time discrete approximation to a gradient flow by solving a sequence of optimization problems. In particular, each step entails minimizing the associated energy of the gradient flow plus a movement limiter term that is, in the classical context of steepest descent with respect to an inner product, simply quadratic. A variety of existing unconditionally stable numerical methods can be recognized as (typically just first order accurate in time) minimizing movement schemes for their associated evolution equations, already requiring the optimization of the energy plus a quadratic term at every time step. Therefore, our approach gives a painless way to extend these to high order accurate in time schemes while maintaining their unconditional stability. In this sense, it can be viewed as a variational analogue of Richardson extrapolation.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a class of unconditionally energy - stable, high - order accurate numerical schemes for handling gradient flow problems. Specifically, the author aims to generalize the existing minimizing movements methods so that they can be applied to a wider range of gradient flow problems and improve the accuracy of time discretization while maintaining unconditional stability.
### Problem Background
Gradient flow is an important class of partial differential equations (PDEs), which describes the dissipation process of system energy over time. A typical gradient flow equation can be expressed as:
\[
u' = -\nabla_H E(u)
\]
where \( u \) is the state variable of the system, \( E:H \to \mathbb{R} \) is the energy functional of the system, \( H \) is a Hilbert space, and \( \langle \cdot, \cdot \rangle \) is the inner product on it. An important property of the gradient flow is that it dissipates energy over time:
\[
\frac{d}{dt} E(u) = \langle \nabla_H E(u), u' \rangle = - \| \nabla_H E(u) \|^2 \leq 0
\]
### Existing Methods and Their Limitations
Existing numerical methods, such as the Backward Euler method, can ensure unconditional stability, but usually have only first - order time accuracy. These methods achieve time discretization at each step by solving an optimization problem:
\[
u^{n + 1} = \arg \min_u \left( E(u) + \frac{1}{2k} \| u - u^n \|^2 \right)
\]
This method is called minimizing movements. However, for higher - order time accuracy, existing methods often find it difficult to simultaneously ensure unconditional stability and high accuracy.
### Core Contributions of the Paper
This paper proposes a new variational extrapolation method, which can generalize the existing minimizing movements methods to high - order accuracy while maintaining unconditional stability. Specifically, the new method has the following features:
1. **Fully General**: There are no additional assumptions (such as convexity) on the energy functional \( E \), only sufficient differentiability is required.
2. **Unconditional Energy Stability**: Ensure energy dissipation at each step.
3. **High - Order Accuracy**: Reach at least third - order time accuracy.
4. **Only a Few Standard Minimizing Movements Solves are Required at Each Time Step**: Equivalent to Backward Euler sub - steps or optimizing the associated energy plus quadratic terms.
### Method Overview
The new method is a linear M - order scheme in the following form:
1. Initialize \( U_0 = u^n \)
2. For \( m = 1,\ldots,M \), solve:
\[
U_m = \arg \min_u \left( E(u) + \sum_{i = 0}^{m - 1} \frac{\gamma_{m,i}}{2k} \| u - U_i \|^2 \right)
\]
3. Set \( u^{n + 1} = U_M \)
By reasonably choosing the coefficients \( \gamma_{m,i} \), the unconditional energy stability and high - order accuracy of the new method can be ensured.
### Numerical Experiments
The paper shows the application of the new method to several classical gradient flow problems, including ordinary differential equations (ODEs) and partial differential equations (PDEs), and verifies the expected convergence rate and energy dissipation properties.
In conclusion, this paper solves the key problem of how to improve the time discretization accuracy of gradient flow problems while maintaining unconditional stability.