Abstract:This work presents a numerical investigation of different approximation techniques for the temporal weights used in the Dual Weighted Residual (DWR) method applied to a time-dependent convection-diffusion equation which is assumed to be convection-dominated. It is a continuation of a previous work by the authors where spatial weights were compared for a steady-state case. A higher-order finite elements approach is compared to a more cost-efficient higher-order reconstruction approach. Numerical examples point out the results regarding accuracy, efficiency and stability reasons.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to study the performance of different time - weight approximation techniques in the Dual Weighted Residual (DWR) method for time - dependent convection - diffusion equations dominated by convection. Specifically, the author focuses on the differences in accuracy, efficiency and stability of these approximation techniques and evaluates the performance of these techniques through numerical experiments.
### Core problems of the paper
1. **Convection - dominated problem**: The paper studies the situation where the convection term is much larger than the diffusion term, which will lead to non - physical oscillations in the numerical solution, etc. Therefore, effective error estimation and adaptive mesh refinement methods are required.
2. **Approximation of time weights**: The time weights in the DWR method are crucial for error estimation, and different approximation techniques will affect the final error estimation results. The paper compares two main time - weight approximation techniques:
- **High - order reconstruction method (hoRe)**: The order of the time - discrete solution is increased through the time reconstruction operator.
- **High - order finite element method (hoFE)**: A higher - order finite element is used to approximate the dual solution.
3. **Error estimation and adaptive algorithms**: The a posteriori error estimator obtained by the DWR method can guide the selection of space - time adaptive algorithms to improve computational efficiency and accuracy.
### Mathematical model
The paper considers the time - dependent convection - diffusion - reaction equation in the following form:
\[
\frac{\partial u}{\partial t}-\nabla\cdot(\nu\nabla u)+\mathbf{b}\cdot\nabla u + \sigma u=f\quad\text{in }Q = \Omega\times(0, T),
\]
where:
- \(u\) is the unknown function,
- \(\nu\) is the diffusion coefficient,
- \(\mathbf{b}\) is the convection velocity vector,
- \(\sigma\) is the reaction coefficient,
- \(f\) is the source term,
- \(\Omega\subset\mathbb{R}^d\) is the spatial region,
- \((0, T)\) is the time interval.
### Main contributions
1. **Numerical study**: Through numerical experiments on different time - weight approximation techniques, their performance in convection - dominated problems is evaluated.
2. **Adaptive algorithm**: A space - time adaptive algorithm based on the DWR method is proposed, which can automatically adjust the space - time grid according to the error estimator.
3. **Comparative analysis**: The differences in accuracy, efficiency and stability between the high - order reconstruction method and the high - order finite element method are compared in detail.
### Conclusion
The paper verifies the effectiveness of the two time - weight approximation techniques through two numerical examples. The results show that although the high - order finite element method shows better accuracy and stability in some cases, the difference between the two methods is not as significant as in the steady - state problem. In addition, the high - order reconstruction method has a lower computational cost and may be more advantageous in practical applications.
In general, this paper provides an effective error estimation and adaptive mesh refinement method for convection - dominated problems, which helps to improve the accuracy and efficiency of numerical simulations.