基于区间双正交小波的二维椭圆界面问题Galerkin格式
Galerkin Scheme Using Biorthogonal Wavelets on Intervals for 2D Elliptic Interface Problems
摘要 Abstract
本文介绍了一种求解二维椭圆界面问题的Wavelet-Galerkin方法,其形式为$-\nabla\cdot(a\nabla u)=f$在$\Omega\backslash \Gamma$中,其中$\Gamma$是在$\Omega$内的光滑界面。变量标量系数$a>0$和源项$f$可能在$\Gamma$处出现不连续性。通过利用由双线性有限元导出的双正交小波基,该基作为$H^1_0(\Omega)$的Riesz基,我们设计了一种策略,使其在$L_2(\Omega)$-范数下的收敛率接近最优:$O(h^2|\log(h)|^2)$,而在$H^1(\Omega)$-范数下为$O(h|\log(h)|)$(相对于逼近阶次)。为了处理$\Gamma$的几何形状以及解$u$的奇异性(其梯度在$\Gamma$处不连续),在界面上引入了额外的小波元素。双正交小波基的对偶部分在证明这些收敛率方面起着关键作用。我们发展了小波的加权Bessel性质,推导了分数Sobolev空间中的各种不等式,并采用有限元论证建立了理论收敛结果。为了实现更高的精度并有效处理高对比度系数$a$,我们的方法类似于无网格方法,依赖于在整个域内及靠近界面处增加小波元素的数量,无需像有限元方法那样重新网格化。与其他求解椭圆界面问题的方法不同,使用$H^1_0(\Omega)$的Wavelet Riesz基确保了系数矩阵的条件数保持较小且一致有界,无论矩阵大小如何。
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.