The numerical solution of diffusion problems in strongly heterogeneous non-isotropic materials

摘要

A new second-order finite-difference algorithm for the numerical solution of diffusion problems in strongly heterogeneous and non- isotropic media is constructed. On problems with rough coefficients or highly nonuniform grids, the new algorithm is superior to all other algorithms we have compared it with. For problems with smooth coefficients on smooth grids, the method is comparable with other second-order methods. The new algorithm is formulated for logically rectangular grids and is derived using the support- operators method.