Abstract:In this paper, we investigate 1D elliptic equations $-\nabla\cdot (a\nabla u)=f$ with rough diffusion coefficients $a$ that satisfy $0<a_{\min}\le a\le a_{\max}<\infty$ and $f\in L_2(\Omega)$. To achieve an accurate and robust numerical solution on a coarse mesh of size $H$, we introduce a derivative-orthogonal wavelet-based framework. This approach incorporates both regular and specialized basis functions constructed through a novel technique, defining a basis function space that enables effective approximation. We develop a derivative-orthogonal wavelet multiscale method tailored for this framework, proving that the condition number $\kappa$ of the stiffness matrix satisfies $\kappa\le a_{\max}/a_{\min}$, independent of $H$. For the error analysis, we establish that the energy and $L_2$-norm errors of our method converge at first-order and second-order rates, respectively, for any coarse mesh $H$. Specifically, the energy and $L_2$-norm errors are bounded by $2 a_{\min}^{-1/2} \|f\|_{L_2(\Omega)} H$ and $4 a_{\min}^{-1}\|f\|_{L_2(\Omega)} H^2$. Moreover, the numerical approximated solution also possesses the interpolation property at all grid points. We present a range of challenging test cases with continuous, discontinuous, high-frequency, and high-contrast coefficients $a$ to evaluate errors in $u, u'$ and $a u'$ in both $l_2$ and $l_\infty$ norms. We also provide a numerical example that both coefficient $a$ and source term $f$ contain discontinuous, high-frequency and high-contrast oscillations. Additionally, we compare our method with the standard second-order finite element method to assess error behaviors and condition numbers when the mesh is not fine enough to resolve coefficient oscillations. Numerical results confirm the bounded condition numbers and convergence rates, affirming the effectiveness of our approach.

Hierarchical Method for Elliptic Problems Using Wavelet

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

A Derivative-Orthogonal Wavelet Multiscale Method for 1D Elliptic Equations with Rough Diffusion Coefficients

Adaptive Multiwavelet-Hierarchical Method for Multiscale Computation

Novel Multilevel Preconditioners for the Systems Arising from Plane Wave Discretization of Helmholtz Equations with Large Wave Numbers

Wavelet Galerkin Method for an Electromagnetic Scattering Problem

Weak Formulation of Finite Element Method Using Wavelet Basis Functions

Two-Level Hierarchical Basis Preconditioner for Elliptic Equations with Jump Coefficients

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

High-order adaptive multiresolution wavelet upwind schemes for hyperbolic conservation laws

A Daubechies wavelet-based method for elastic problems

Solving the Euler equations by wavelet interpolation

Multi-Scale Daubechies Wavelet-Based Method for 2-D Elastic Problems

A Coiflet Wavelet Homotopy Technique for Nonlinear PDEs: Application to the Extreme Bending of Orthotropic Plate with Forced Boundary Constraints

Adaptive BDDC Algorithms for the System Arising from Plane Wave Discretization of Helmholtz Equations

A Dynamic Multiscale Lifting Computation Method Using Daubechies Wavelet

On Wavelet-Based Numerical Homogenization

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

A Novel Least Squares Method for Helmholtz Equations with Large Wave Numbers

A HAM-based wavelet approach for nonlinear partial differential equations: Two dimensional Bratu problem as an application

Numerical Solution of Poisson Equation with Wavelet Bases of Hermite Cubic Splines on the Interval