An adaptive non-symmetric finite volume and boundary element coupling method for a fluid mechanics interface problem

Christoph Erath,Robert Schorr

DOI: https://doi.org/10.48550/arXiv.1605.07031

2016-05-23

Abstract:We consider an interface problem often arising in transport problems: a coupled system of partial differential equations with one (elliptic) transport equation on a bounded domain and one equation (in this case the Laplace problem) on the complement, an unbounded domain. Based on the non-symmetric coupling of the finite volume method and boundary element method of [Erath et al., <a class="link-https" data-arxiv-id="1509.00440" href="https://arxiv.org/abs/1509.00440">arXiv:1509.00440</a>, 2015] we introduce a robust residual error estimator. The upper bound of the error in an energy (semi)norm is robust against variation of the model data. The lower bound, however, additionally depends on the Peclet number. In several examples we use the local contributions of the a~posteriori error estimator to steer an adaptive mesh-refining algorithm. The adaptive FVM-BEM coupling turns out to be an efficient method especially to solve problems from fluid mechanics, mainly because of the local flux conservation and the stable approximation of convection dominated problems.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the numerical simulation and solution of interface problems in fluid mechanics. Specifically, it focuses on the coupling system between the flow and transport problems described by the elliptic transport equation in a bounded domain and the uniform diffusion process in the corresponding unbounded external region. In order to accurately approximate such problems, the author introduced an asymmetric finite volume method (FVM) and boundary element method (BEM) coupling method and developed a robust residual a posteriori error estimator.

### Specific description of the problem

1. **Model problem**:
- The transport equation in the internal region \(\Omega\) is an elliptic equation.
- The diffusion equation in the external region \(\Omega_e=\mathbb{R}^d\setminus\Omega\) is the Laplace equation.
- On the interface \(\Gamma = \partial\Omega=\partial\Omega_e\), there are inflow and outflow parts \(\Gamma_{\text{in}}\) and \(\Gamma_{\text{out}}\).

2. **Objectives**:
- Develop an effective numerical method to solve this coupling system.
- Provide a robust a posteriori error estimator to guide the adaptive mesh refinement algorithm, thereby improving the solution accuracy.

### Main contributions

1. **Asymmetric FVM - BEM coupling method**:
- Introduced an asymmetric FVM - BEM coupling method, which uses FVM in the internal region and BEM in the external region.
- By transforming the external problem into an integral equation, the problem of truncating the unbounded external region is avoided.

2. **Robust a posteriori error estimator**:
- Proposed a robust residual a posteriori error estimator to ensure the reliability (upper bound) and efficiency (lower bound) of error estimation.
- The error estimator is robust in the energy (semi -) norm, especially in convection - dominated problems.

3. **Adaptive mesh refinement**:
- Use the local contribution of the a posteriori error estimator to guide the adaptive mesh refinement algorithm, which improves the solution efficiency.

### Mathematical formula representation

- **Energy (semi -) norm**:
\[
|||v|||^2_\Omega:=\|A^{1/2}\nabla v\|^2_{L^2(\Omega)}+\left\|\left(\frac{1}{2}\text{div}\,b + c\right)^{1/2}v\right\|^2_{L^2(\Omega)}
\]

- **Bilinear form**:
\[
B((u,\varphi);(v,\psi)) := A(u,v)-\langle\varphi, v\rangle_\Gamma+\langle\psi,(\frac{1}{2}-K)u\rangle_\Gamma+\langle\psi, V\varphi\rangle_\Gamma
\]

- **Residual**:
\[
R := f-\text{div}(-A\nabla u_h + bu_h)-cu_h
\]

- **Jump term**:
\[
J|_E :=
\begin{cases}
[(-A\nabla u_h)|_{E,T}-(-A\nabla u_h)|_{E,T'}]\cdot n & \text{if } E\in E_I\\
(-A\nabla u_h + bu_h)\cdot n+\varphi_h + t_0 & \text{if } E\in E_{\Gamma_{\text{in}}} \\
-A\na