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.

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

Efficient Chebyshev Spectral Method for Solving Linear Elliptic Pdes Using Quasi-Inverse Technique

Application of thin plate splines for solving a class of boundary integral equations arisen from Laplace's equations with nonlinear boundary conditions

A Periodic Wavelet Method for the Second Kind of the Logarithmic Integral Equation

On the Approximate Solution of Integral Equations with Logarithmic Kernels Using the Third Kind of Chebyshev Polynomials

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

A new Laguerre wavelets‐based method for solving Fredholm integral equations with weakly singular logarithmic kernel

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

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

Galerkin-Chebyshev Spectral Method and Block Boundary Value Methods for Two-Dimensional Semilinear Parabolic Equations

Application of linear Legendre multi-wavelets collocation method for solution of fourth order boundary value problems

Multiwavelet Compression for the Boundary Integral Equation on an Open Wedge

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

Numerical solution of a weakly singular integral equation by the method of orthogonal polynomials

Trigonometric Hermite Wavelet Approximation for the Integral Equations of Second Kind with Weakly Singular Kernel

Spectral Monic Chebyshev Approximation for Higher Order Differential Equations

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

CALCULATION OF LOGARITHMIC SINGULARITIES INTEGRAL IN DUAL BOUNDARY ELEMENT METHODS

Fast Fourier-Galerkin Methods for First-Kind Logarithmic-Kernel Integral Equations on Open Arcs

Implicit-Explicit Finite Difference Approximations of a Semilinear Heat Equation with Logarithmic Nonlinearity

On the boundary element method for some nonlinear boundary value problems