Numerical Homogenization of Fractal Interface Problems

Ralf Kornhuber,Joscha Podlesny,Harry Yserentant

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

2020-07-22

Abstract:We consider the numerical homogenization of a class of fractal elliptic interface problems inspired by related mechanical contact problems from the geosciences. A particular feature is that the solution space depends on the actual fractal geometry. Our main results concern the construction of projection operators with suitable stability and approximation properties. The existence of such projections then allows for the application of existing concepts from localized orthogonal decomposition (LOD) and successive subspace correction to construct first multiscale discretizations and iterative algebraic solvers with scale-independent convergence behavior for this class of problems.

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to perform numerical homogenization for a class of fractal interface problems inspired by mechanical contact problems in geoscience. The characteristic of this type of problems is that the solution space depends on the actual fractal geometric structure. Specifically, the paper focuses on how to construct projection operators with appropriate stability and approximation properties, and use these projection operators to apply the existing Localized Orthogonal Decomposition (LOD) and successive subspace correction methods, thereby constructing multi - scale discretizations and iterative algebraic solvers for this type of problems, ensuring that their convergence behavior is independent of the scale.

### Main Research Contents

1. **Problem Background**:
- The paper considers a class of elliptic fractal interface problems inspired by mechanical contact problems in geoscience.
- The solution spaces of these problems depend on the actual fractal geometric structure, so traditional homogenization methods (such as the classical periodicity assumption) are no longer applicable.

2. **Objectives**:
- Construct projection operators $\Pi_k: H \to S_k$ with appropriate stability and approximation properties.
- Use these projection operators to apply the existing LOD and successive subspace correction methods, construct multi - scale discretizations and iterative algebraic solvers, and ensure that their convergence behavior is independent of the scale.

3. **Methods**:
- Through local Poincaré inequalities, trace lemmas and quasi - interpolation methods, construct projection operators that satisfy the required properties.
- Prove the existence of these projection operators and analyze their stability and approximation properties.

4. **Applications**:
- Apply the LOD method and subspace correction method to construct multi - scale discretizations and iterative solvers.
- Verify the effectiveness of the proposed method through numerical experiments and show its potential applications in geology.

### Mathematical Formulas

- Stability and approximation properties of projection operators:
\[
\|\Pi v\|_H \leq c \|v\|_H, \quad \|v - \Pi v\|_{L^2} \leq Ch \|v\|_H \quad \forall v \in H
\]

- Local Poincaré inequality:
\[
\|v - \overline{v}\|_{0,B}^2 \leq C r_G \left( r_G |\nabla v|^2_{B \setminus \Gamma} + \sum_{j = k+1}^\infty (1 + c)^{j - k} C_{k,j} \|\llbracket v \rrbracket\|^2_{0, \Gamma_j \cap B(G, r_G)} \right)
\]
where $\overline{v}=\frac{1}{|B|}\int_B v \, dx$ represents the average value of $v$ on the ball $B$.

Through these methods and tools, the paper aims to provide effective numerical solutions for complex fractal interface problems, especially for multi - scale problems in geoscience.