Galerkin Scheme Using Biorthogonal Wavelets on Intervals for 2D Elliptic Interface Problems

Abstract: This paper introduces a wavelet Galerkin method for solving two-dimensional elliptic interface problems of the form $-\nabla\cdot(a\nabla u)=f$ in $\Omega\backslash \Gamma$, where $\Gamma$ is a smooth interface within $\Omega$. The variable scalar coefficient $a>0$ and source term $f$ may exhibit discontinuities across $\Gamma$. By utilizing a biorthogonal wavelet basis derived from bilinear finite elements, which serves as a Riesz basis for $H^1_0(\Omega)$, we devise a strategy that achieves nearly optimal convergence rates: $O(h^2 |\log(h)|^2)$ in the $L_2(\Omega)$-norm and $O(h |\log(h)|)$ in the $H^1(\Omega)$-norm with respect to the approximation order. To handle the geometry of $\Gamma$ and the singularities of the solution $u$, which has a discontinuous gradient across $\Gamma$, additional wavelet elements are introduced along the interface. The dual part of the biorthogonal wavelet basis plays a crucial role in proving these convergence rates. We develop weighted Bessel properties for wavelets, derive various inequalities in fractional Sobolev spaces, and employ finite element arguments to establish the theoretical convergence results. To achieve higher accuracy and effectively handle high-contrast coefficients $a$, our method, much like meshfree approaches, relies on augmenting the number of wavelet elements throughout the domain and near the interface, eliminating the need for re-meshing as in finite element methods. Unlike all other methods for solving elliptic interface problems, the use of a wavelet Riesz basis for $H^1_0(\Omega)$ ensures that the condition numbers of the coefficient matrices remain small and uniformly bounded, regardless of the matrix size.