Wavelet Galerkin method for fourth-order multi-dimensional elliptic partial differential equations

Abstract

We describe a wavelet Galerkin method for numerical solutions of fourth-order linear and nonlinear partial differential equations (PDEs) in 2D and 3D based on the use of Daubechies compactly supported wavelets. Two-term connection coefficients have been used to compute higher-order derivatives accurately and economically. Localization and orthogonality properties of wavelets make the global matrix sparse. In particular, these properties reduce the computational cost significantly. Linear system of equations obtained from discretized equations have been solved using GMRES iterative solver. Quasi-linearization technique has been effectively used to handle nonlinear terms arising in nonlinear biharmonic equation. To reduce the computational cost of our method, we have proposed an efficient compression algorithm. Error and stability estimates have been derived. Accuracy of the proposed method is demonstrated through various examples.