Wavelet optimized finite difference method using interpolating wavelets for self-adjoint singularly perturbed problems

We design a wavelet optimized finite difference (WOFD) scheme for solving self-adjoint singularly perturbed boundary value problems. The method is based on an interpolating wavelet transform using polynomial interpolation on dyadic grids. Small dissipation of the solution is captured significantly using an adaptive grid. The adaptive feature is performed automatically by thresholding the wavelet coefficients. Numerical examples have been solved and compared with non-standard finite difference schemes in [J.M.S. Lubuma, K.C. Patidar, Uniformly convergent non-standard finite difference methods for self-adjoint singular perturbation problems, J. Comput. Appl. Math. 191 (2006) 228-238]. The proposed method outperforms the non-standard finite difference for studying singular perturbation problems for small dissipations (very small @e) and effective grid generation. Therefore, the proposed method is better for studying the more challenging cases of singularly perturbed problems.

[1]  David L. Donoho,et al.  Interpolating Wavelet Transforms , 1992 .

[2]  R. Bellman,et al.  Quasilinearization and nonlinear boundary-value problems , 1966 .

[3]  Gilles Deslauriers,et al.  Symmetric Iterative Interpolation Processes , 1989 .

[4]  Mani Mehra,et al.  Cubic spline Adaptive Wavelet Scheme to Solve singularly perturbed Reaction Diffusion Problems , 2007, Int. J. Wavelets Multiresolution Inf. Process..

[5]  Lutz Tobiska,et al.  Numerical Methods for Singularly Perturbed Differential Equations , 1996 .

[6]  A. Harten Adaptive Multiresolution Schemes for Shock Computations , 1994 .

[7]  Naoki Saito,et al.  Multiresolution representations using the autocorrelation functions of compactly supported wavelets , 1993, IEEE Trans. Signal Process..

[8]  Wil H. A. Schilders,et al.  Uniform Numerical Methods for Problems with Initial and Boundary Layers , 1980 .

[9]  Gilles Deslauriers,et al.  Symmetric iterative interpolation processes , 1989 .

[10]  JohnM . Miller,et al.  Robust Computational Techniques for Boundary Layers , 2000 .

[11]  J. J. Miller,et al.  Fitted Numerical Methods for Singular Perturbation Problems , 1996 .

[12]  J. Lubuma,et al.  Uniformly convergent non-standard finite difference methods for self-adjoint singular perturbation problems , 2006 .

[13]  Leland Jameson,et al.  A Wavelet-Optimized, Very High Order Adaptive Grid and Order Numerical Method , 1998, SIAM J. Sci. Comput..

[14]  H. Weinberger,et al.  Maximum principles in differential equations , 1967 .

[15]  Mohan K. Kadalbajoo,et al.  Fitted mesh B-spline collocation method for solving self-adjoint singularly perturbed boundary value problems , 2005, Appl. Math. Comput..