Abstract:We present reduced basis approximations and rigorous a posteriori error bounds for the instationary Stokes equations. We shall discuss both a method based on the standard formulation as well as a method based on a penalty approach, which combine techniques developed in our previous work on parametrized saddle point problems with current reduced basis techniques for parabolic problems. The analysis then shows how time integration affects the development of reduced basis a posteriori error bounds as well as the construction of computationally efficient reduced basis approximation spaces. To demonstrate their performance in practice, the methods are applied to a Stokes flow in a two-dimensional microchannel with a parametrized rectangular obstacle; evolution in time is induced by a time-dependent velocity profile on the inflow boundary. Numerical results illustrate (i) the rapid convergence of reduced basis approximations, (ii) the performance of a posteriori error bounds with respect to sharpness, and (iii) computational efficiency.
What problem does this paper attempt to address?
This paper attempts to solve the problem of efficient numerical solution of the unsteady Stokes equations (transient Stokes equations), especially to provide reliable reduced - basis (RB) approximations and strict a posteriori error bounds for parametric problems. Specifically, the authors propose two methods:
1. **Method based on the standard form**: Directly deal with the standard mixed variational form of the unsteady Stokes equations.
2. **Penalty - function - based method**: Introduce a penalty term in the continuity equation, thus avoiding the problem of estimating the inf - sup stability constant.
These two methods combine the existing reduced - basis techniques and the reduced - basis techniques for parabolic problems, aiming to improve computational efficiency and ensure the strictness and accuracy of the error bounds.
### Main research contents
- **Time discretization**: Use the backward Euler method to discretize time, and regard each time step as a saddle - point problem.
- **Reduced - basis approximation**: Construct an efficient approximate solution by Galerkin projection into a low - dimensional reduced - basis space.
- **A posteriori error bounds**: Develop strict a posteriori error bounds, which take into account not only the influence of the diffusion term and the incompressibility constraint, but also the influence of the mass term.
- **Computational efficiency**: Through the off - line - on - line computing strategy, ensure that the computational cost in the on - line stage does not depend on the dimension of the high - dimensional finite - element space, but only on the smaller dimension of the reduced - basis space.
- **Adaptive sampling**: Use the POD greedy algorithm to construct the reduced - basis space to capture the parameter dependence and time - evolution characteristics and ensure the accuracy of the approximation.
### Application examples
To verify the effectiveness of the proposed methods, the authors apply them to the Stokes flow problem in a two - dimensional microchannel, where the time evolution is induced by the time - varying velocity profile at the inlet boundary. The numerical results show the following points:
1. **Fast convergence**: The fast convergence of the reduced - basis approximate solution.
2. **Accuracy of error bounds**: The effectiveness and accuracy of the a posteriori error bounds.
3. **Computational efficiency**: The high efficiency of the method in practical applications.
### Formula representation
The formulas involved in the paper are represented in Markdown format as follows:
- Unsteady Stokes equations with time discretization:
\[
\frac{1}{\Delta t} m(u^k - u^{k - 1}, v; \mu) + a(u^k, v; \mu) + b(v, p^k; \mu) = f^k(v; \mu) \quad \forall v \in X
\]
\[
b(u^k, q; \mu) - \varepsilon c(p^k, q; \mu) = g^k(q; \mu) \quad \forall q \in Y
\]
- Definition of a posteriori error bounds (for the case of \(\varepsilon = 0\)):
\[
\Delta_k^N(\mu) = \left[ \Delta t \sum_{j = 1}^k \left( \frac{\|r_1^j\|^2_{X'}}{\alpha_a^{LB}(\mu)} + \frac{2}{\beta^{LB}(\mu)} \left( 1 + \frac{\gamma_a^{UB}(\mu)}{\alpha_a^{LB}(\mu)} \right) \|r_1^j\|_{X'} \|r_2^j\|_{Y'} + \left( \frac{\gamma_m^{UB}(\mu)}{\Delta t} + \frac{(\gamma_a^{UB}(\mu))^2}{\alpha_a^{LB}(\mu)} \right) \|r_2^j\|^2_{Y'} \right) \right]^{1/2}
\]
These formulas show how to construct a posteriori error bounds through the dual norms of residuals and the continuity, coercivity and inf - sup conditions of the true solution.