Jump-preserving polynomial interpolation in non-manifold polyhedra

Martin Averseng
2024-04-30
Abstract:We construct a piecewise-polynomial interpolant $u \mapsto \Pi u$ for functions $u:\Omega \setminus \Gamma \to \mathbb{R}$, where $\Omega \subset \mathbb{R}^d$ is a Lipschitz polyhedron and $\Gamma \subset \Omega$ is a possibly non-manifold $(d-1)$-dimensional hypersurface. This interpolant enjoys approximation properties in relevant Sobolev norms, as well as a set of additional algebraic properties, namely, $\Pi^2 = \Pi$, and $\Pi$ preserves homogeneous boundary values and jumps of its argument on $\Gamma$. As an application, we obtain a bounded discrete right-inverse of the "jump" operator across $\Gamma$, and an error estimate for a Galerkin scheme to solve a second-order elliptic PDE in $\Omega$ with a prescribed jump across $\Gamma$.
Numerical Analysis,Analysis of PDEs
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to construct a piecewise polynomial interpolation operator $\Pi_h$ for the interpolation of the function $u: \Omega \setminus \Gamma \to \mathbb{R}$, where $\Omega \subset \mathbb{R}^d$ is a Lipschitz polyhedron and $\Gamma \subset \Omega$ is a possibly non - manifold $(d - 1)$ - dimensional hypersurface. This interpolation operator has approximation properties in the relevant Sobolev norms, as well as a set of additional algebraic properties, specifically including: 1. **Idempotency**: $\Pi_h^2=\Pi_h$. 2. **Preserving homogeneous boundary conditions**: For all $u \in H_0^1(\Omega \setminus \Gamma)$, we have $\Pi_h u \in H_0^1(\Gamma, \Omega)$. 3. **Preserving jumps**: $\Pi_h$ preserves the jumps of its arguments on $\Gamma$. In addition, the paper also provides a bounded discrete right - inverse "jump" operator and gives a Galerkin scheme to solve the error estimate of the second - order elliptic partial differential equations (PDEs) with predetermined jumps. ### Specific problem description 1. **Construction of the interpolation operator**: - The paper first defines the interpolation operator $\Pi_h: H^1(\Omega \setminus \Gamma)\to V_p(\Omega_h; \Gamma)$, where $V_p(\Omega_h; \Gamma)$ is a finite - dimensional subspace consisting of $p$-th degree polynomials on each mesh element of $\Omega_h$. - The interpolation operator $\Pi_h$ needs to satisfy the properties mentioned above, especially preserving jumps and homogeneous boundary conditions. 2. **Approximation properties**: - The paper proves the continuity and approximation properties of the interpolation operator $\Pi_h$ in Sobolev norms. Specifically, for a given integer $p\geq1$ and a real number $\frac{1}{2}<t\leq p + 1$, there exists a constant $C(\gamma_0, p, t)>0$ such that if $\Omega_h$ is a $\gamma_0$-shape - regular mesh, then the interpolation operator $\Pi_h$ satisfies the following continuity properties: \[ \|\Pi_h u\|_{H^t(K)}^2\leq C(\gamma_0, p, t)\sum_{K'\in\omega_K}\|u\|_{H^t(K')}^2, \] \[ \|u - \Pi_h u\|_{H^s(K)}^2\leq C(\gamma_0, p, t)h_K^{2(t - s)}\sum_{K'\in\omega_K}|u|_{H^t(K')}^2, \] where $s\in[0, t]$, $h_K$ is the diameter of the mesh element $K$, and $\omega_K$ is the set of all mesh elements that share at least one vertex with $K$. 3. **Applications**: - **Bounded discrete right - inverse jump operator**: The paper constructs a bounded discrete right - inverse jump operator to handle jump problems. - **Error estimate of the Galerkin scheme**: The paper gives a Galerkin scheme to solve the second - order elliptic PDEs with predetermined jumps and provides error estimates. ### Summary The paper solves the problem of numerical solution of partial differential equations in complex geometric structures (especially polyhedra containing non - manifold hypersurfaces) by constructing an interpolation operator $\Pi_h$ with specific properties. This method not only provides theoretical guarantees but also provides an effective tool for numerical simulation in practical applications.