Abstract:This paper introduces a wavelet Galerkin method for solving two-dimensional elliptic interface problems of the form in $-\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.

Wavelet-galerkin Method for the Singular Perturbation Problem with Boundary Layers

High Accuracy Non-Equidistant Method for Singular Perturbation Reaction-Diffusion Problem

Positive Solutions to a Singular Second Order Three-Point Boundary Value Problem

A higher order numerical method for singularly perturbed elliptic problems with characteristic boundary layers

Numerical Solution of Singularly Perturbed Problems via both Galerkin and Subdomain Galerkin methods

A wavelet method for solving a class of nonlinear boundary value problems

A Finite-Element Method for Laplace- and Helmholtz-Type Boundary Value Problems with Singularities

Singularly Perturbed Boundary Value Problems Via the Fixed Point Theorem

A meshless method using wavelets

A sixth-order wavelet integral collocation method for solving nonlinear boundary value problems in three dimensions

Non-singular boundary integral methods for fluid mechanics applications

Boundary elements with mesh refinements for the wave equation

A high-precision method for general nonlinear initial-boundary value problems

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

On the numerical solution of a hyperbolic initial boundary value problem by hypersingular boundary integral equations

Singular perturbation boundary and interior layers problems with multiple turning points

Boundary value problems initial condition identification by a wavelet-based Galerkin method

A Tailored Finite Point Method for a Singular Perturbation Problem on an Unbounded Domain

A Weakly Singular Boundary Integral Formulation Of The External Helmholtz Problem Valid For All Wavenumbers

Wavelet-based Edge Multiscale Finite Element Methods for Singularly Perturbed Convection-Diffusion Equations

Numerical and mathematical analysis of nonlocal singular Emden–Fowler type BVPs by improved Taylor-wavelet method