Diffusion coefficients at control-volume surfaces are required for numerical simulation of 10 convection-diffusion equations. Various interpolation methods for these diffusion coefficients are briefly discussed and compared with the analytical solutions of both pure diffusion and convection-diffusion problems. It is found that the harmonic mean method is not as accurate and reliable as is supposed. A new method is thus developed. Extensive numerical comparisons are given for both the harmonic mean and the present method, and the results 15 show that the method proposed is accurate, reliable, and easy to use.