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.

Euler wavelets operational matrix of integration and its application in the calculus of variations

Solving nonlinear Volterra integro-differential equations of fractional order by using Euler wavelet method

An Efficient Numerical Method Based on Bell Wavelets for Solving the Fractional Integro-Differential Equations with Weakly Singular Kernels

Solving the Euler equations by wavelet interpolation

A Wavelet-based Method for Numerical Calculation of Singular Integrals in Eddy Current Analysis

Weak Formulation of Finite Element Method Using Wavelet Basis Functions

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

Solution of linear and nonlinear singular value problems using operational matrix of integration of Fibonacci wavelets

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

Water Wave Solutions Obtained by Variational Method

Waveform Relaxation Methods of Nonlinear Integral-Differential-Algebraic Equations

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

A new generalized Bell wavelet and its applications for solving linear and nonlinear integral equations

Shannon wavelet approximations of linear differential operators

Variational Integrals of a Class of Nonhomogeneous 𝒜-Harmonic Equations

An adaptive wavelet method for nonlinear differential-algebraic equations

A Daubechies wavelet-based method for elastic problems

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

On the generalized wavelet-Galerkin method

An Application of the Gegenbauer Wavelet Method for the Numerical Solution of the Fractional Bagley-Torvik Equation

Numerical Inversion of Laplace Transform via Wavelet Operational Matrix and Its Applications to Fractional Differential Equations