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.

Compression of boundary integral operators discretized by anisotropic wavelet bases

Boundary Knot Method for Some Inverse Problems Associated with the Helmholtz Equation

Wavelet compressed, modified Hilbert transform in the space-time discretization of the heat equation

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

Boundary elements with mesh refinements for the wave equation

A fast direct solver for nonlocal operators in wavelet coordinates

Hybrid matrix compression for high-frequency problems

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

Multi-level, Multi-frequency, compressed Full Waveform Inversion

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

Weak Formulation of Finite Element Method Using Wavelet Basis Functions

A Fast Fourier-Galerkin Method for Solving Boundary Integral Equations on Torus-Shaped Surfaces

Application of dual-Chebyshev wavelets for the numerical solution of boundary integral equations with logarithmic singular kernels

A Fast Fourier-Galerkin Method Solving Boundary Integral Equations for the Helmholtz Equation with Exponential Convergence

A Fully Discrete Fast Fourier–Galerkin Method Solving a Boundary Integral Equation for the Biharmonic Equation

On the boundary element method for some nonlinear boundary value problems

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

Uniform $\mathcal{H}$-matrix Compression with Applications to Boundary Integral Equations

On the multistep time discretization of linear\newline initial-boundary value problems and their boundary integral equations

Hierarchical Method for Elliptic Problems Using Wavelet

Convolution quadrature for the wave equation with a nonlinear impedance boundary condition