Abstract
It is a well-known phenomenon called superconvergence in the mathematical literature that the error level of an integral quantity can be much smaller than the magnitude of the local errors involved in the computation of this quantity. When discretizing an integrated form of Fick's second law of diffusion the local errors reflect the accuracy of individual concentration points while the integral quantity has the physical meaning of the flux. This article demonstrates how an extraordinary fast exponential convergence towards zero can be achieved for the simulated flux error on the basis of finite–difference approximations that are only second-order (Box 2 method) or fourth-order (Box 4 method) accurate as far as the level of local errors is concerned. © 2005 Wiley Periodicals, Inc. J Comput Chem 26: 619–632, 2005

This publication has 24 references indexed in Scilit: