Sarah K. Locke,John R. Singler
Abstract:In a recent work [B. Koc et al., <a class="link-https" data-arxiv-id="2010.03750" href="https://arxiv.org/abs/2010.03750">arXiv:2010.03750</a>, SIAM J. Numer. Anal., to appear], the authors showed that including difference quotients (DQs) is necessary in order to prove optimal pointwise in time error bounds for proper orthogonal decomposition (POD) reduced order models of the heat equation. In this work, we introduce a new approach to including DQs in the POD procedure. Instead of computing the POD modes using all of the snapshot data and DQs, we only use the first snapshot along with all of the DQs and special POD weights. We show that this approach retains all of the numerical analysis benefits of the standard POD DQ approach, while using a POD data set that has half the number of snapshots as the standard POD DQ approach, i.e., the new approach is more computationally efficient. We illustrate our theoretical results with numerical experiments.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: when using the proper orthogonal decomposition (POD) method to construct reduced order models (ROMs), how to optimize the pointwise - in - time error bounds by introducing difference quotients (DQs) while improving computational efficiency. Specifically, the author proposes a new method that uses only the first data snapshot and all difference quotients, instead of using all data snapshots and difference quotients, to avoid data redundancy and reduce computational cost.
### Background and Problem Description of the Paper
1. **Limitations of Existing Methods**
- In the traditional POD method, in order to obtain the optimal pointwise - in - time error bounds, all data snapshots and difference quotients are usually used together to calculate the POD modes. This will lead to linear correlation of the data set and increase the computational complexity.
- The size of the data set in the traditional method is almost twice that of the new method, resulting in lower computational efficiency.
2. **Research Motivation**
- Researchers hope to find a method that can reduce data redundancy and improve computational efficiency without losing the advantages of numerical analysis.
- The specific goal is to verify whether it is possible to use difference quotients to retain all the numerical analysis advantages of the POD method while reducing the computational burden by selecting specific data sets and weights.
### Core Idea of the New Method
The author proposes a new POD method that uses only the first data snapshot and all difference quotients and introduces special POD weights. The key features of this method are as follows:
- **Data Set Selection**: Only use the first data snapshot \( u_1 \) and all difference quotients \( \partial u_j \), instead of all data snapshots and difference quotients.
- **Linear Independence**: If the original data set is linearly independent, then the data set of the new method is also linearly independent, thus avoiding redundant information.
- **Computational Efficiency**: The size of the data set used in the new method is \( N \), while the size of the data set in the traditional method is \( 2N - 1 \), so the new method is more efficient.
### Mathematical Formula Representation
The error function of the new method is defined as:
\[ E_{\text{DQ1}}^r=\|u_1 - \Pi_X^r u_1\|^2_X+\sum_{j = 1}^{N - 1}\Delta t\|\partial u_j-\Pi_X^r\partial u_j\|^2_X \]
where:
- \( u_1 \) is the first data snapshot,
- \( \partial u_j=\frac{u_{j + 1}-u_j}{\Delta t} \) is the difference quotient,
- \( \Pi_X^r \) is the orthogonal projection onto the subspace \( X_r \).
### Main Conclusions
Through theoretical analysis and numerical experiments, the author proves that the new method not only retains all the numerical analysis advantages of the traditional POD method but also significantly improves computational efficiency. Specifically:
- **Point Error Bounds**: The new method can obtain similar pointwise - in - time error bounds as the traditional method.
- **Computational Efficiency**: The size of the data set used in the new method is halved, reducing the computational complexity.
In conclusion, the main contribution of this paper is to propose a more efficient POD method that can reduce computational cost while maintaining the advantages of numerical analysis.