Erik Burman,Rebecca Durst,Miguel A. Fernández,Johnny Guzmán,Sijing Liu
Abstract:We consider a loosely coupled, non-iterative Robin-Robin coupling method proposed and analyzed in [J. Numer. Math., 31(1):59--77, 2023] for a parabolic-parabolic interface problem and prove estimates for the discrete time derivatives of the scalar field in different norms. When the interface is flat and perpendicular to two of the edges of the domain we prove error estimates in the $H^2$-norm. Such estimates are key ingredients to analyze a defect correction method for the parabolic-parabolic interface problem. Numerical results are shown to support our findings.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is: For the parabolic - parabolic interface problem, when using the loosely - coupled, non - iterative Robin - Robin coupling method, how to estimate the error of the discrete - time derivative in different norms and prove the validity of these estimates.
Specifically, the paper focuses on the following points:
1. **Estimation of discrete - time derivative**:
- The paper aims to prove the error estimates in different norms of the time - discrete derivative of the scalar field when using the Robin - Robin coupling method to solve the parabolic - parabolic interface problem.
- In particular, when the interface is flat and perpendicular to two sides of the domain, the author proves the second - order convergence in the \(H^2\) norm.
2. **Application of error estimates**:
- These error estimates are the key to analyzing the defect correction method, which aims to improve the first - order time convergence to second - order.
- The defect correction method depends on the prediction step and the correction step. The prediction step is exactly the Robin - Robin splitting method analyzed in this paper, and the second - order accuracy of the correction step depends on the second - order accuracy of the time difference in the prediction step.
3. **Numerical verification**:
- The validity of the theoretical results is verified through numerical experiments, showing the convergence rates of error estimates in different configurations.
### Summary of mathematical formulas
- **Discrete - time derivative**:
\[
\partial_{\Delta t} v^{n + 1}=\frac{v^{n + 1}-v^{n}}{\Delta t}
\]
\[
\partial^2_{\Delta t} v^{n + 1}=\frac{v^{n + 1}-2v^{n}+v^{n - 1}}{(\Delta t)^2}
\]
- **Error estimate**:
- For the general case, the error \(E\) satisfies:
\[
Z_N(W, U, \Lambda)+\sum_{n = 0}^{N - 1}S_{n + 1}(W, U, \Lambda)\leq C(\Delta t)^2Y
\]
- In the special case (such as when the interface is perpendicular to the boundary), the error estimate in the \(H^2\) norm is:
\[
\Delta t\sum_{n = 0}^{N - 1}\nu_f\|D^2(\partial_{\Delta t}U^{n + 1})\|^2_{L^2(\Omega_f)}\leq C(\Delta t)^2(Y + Y'+\nu_f\|\partial^3_tw\|^2_{L^2((0, T), L^2(\Omega_s))})
\]
Here, \(Y\) and \(Y'\) are norms defined in specific Sobolev spaces, representing the boundedness of the higher - order derivatives of the solution.
Through these estimates, the author not only proves the high - order convergence in time of the Robin - Robin method, but also provides a theoretical basis for the subsequent defect correction method.