Superconvergence of solution derivatives for the Shortley-Weller difference approximation of Poisson's equation. Part I: smoothness problems

摘要

The finite difference method (FDM) using the Shortley-Weller approximation can be viewed as a special kind of the finite element methods (FEMs) using the piecewise bilinear and linear functions, and involving some integration approximation. When u ∈ C~3(S-bar) (i.e., u ∈C~(3,0)(S-bar)) and f ∈ C~2(S-bar), the superconvergence rate O(h~2) of solution derivatives in discrete H~1 norms by the FDM is derived for rectangular difference grids, where h is the maximal mesh length of difference grids used, and the difference grids are not confined to be quasiunform. Comparisons are made on the analysis by the maximum principle and the FEM analysis, conversions between the FDM and the linear and bilinear FEMs are discussed, and numerical experiments are provided to support superconvergence analysis made.