We propose a new high order finite difference discretization strategy, which is based on the Richardson extrapolation technique and an operator interpolation scheme, to solve convection diffusion equations. For a particular implementation, we solve a fine grid equation and a coarse grid equation by using a fourth order compact difference scheme. Then we combine the two approximate solutions and use the Richardson extrapolation to compute a sixth order accuracy coarse grid solution. A sixth order accuracy fine grid solution is obtained by interpolating the sixth order coarse grid solution using an operator interpolation scheme. Numerical results are presented to demonstrate the accuracy and efficacy of the proposed finite difference discretization strategy, compared to the sixth order combined compact difference (CCD) scheme, and the standard fourth order compact difference (FOC) scheme.
Mathematics Subject Classification: 65N06, 65N55, 65F10.
Technical Report 347-02, Department of Computer Science, University of Kentucky, Lexington, KY, 2002.
This research was supported in part by the U.S. National Science Foundation under the grant CCR-9902022, CCR-9988165, CCR-0092532, and ACI-0202934, in part by the U.S. Department of Energy under grant DE-FG02-02ER45961, in part by the Japan Research Organization for Information Science & Technology, and in part by the University of Kentucky Research Committee.