A derivative-orthogonal wavelet multiscale method for elliptic equations with rough diffusion coefficients

In this paper, we investigate 1D elliptic equations −∇·(a∇u)=f with rough diffusion coefficients a satisfying 0 < a _min ≤ a ≤ a _max < ∞ and rough source terms f ∈ L ₂(Ω). 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 κ of the stiffness matrix satisfies κ ≤ a _max/ a _min, independent of H . For the error analysis, we establish that the energy and L ₂-norm errors of our method converge at first-order and second-order rates, respectively, for any coarse mesh H . Specifically, the energy and L ₂-norm errors are bounded by 2a_min^−1/2∥f∥_L2(Ω)H and 4a_min⁻¹∥f∥_L2(Ω)H². 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 au ′ in both l ₂ and l _∞ 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 (FEM) and the special FEM in [6] 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. Thus, our method is capable of handling both the rough diffusion coefficient a and the rough source term f . In the special case that a is constant but f is rough, our method achieves the optimal condition number κ=1 and second-order L ₂ convergence without solving any equations.