Abstract:The high resolution numerical perturbation (NP) algorithm is analyzed and tested using various convective-diffusion equations. The NP algorithm is constructed by splitting the second order central difference schemes of both convective and diffusion terms of the convective-diffusion equation into upstream and downstream parts, then the perturbation reconstruction functions of the convective coefficient are determined using the power-series of grid interval and eliminating the truncated errors of the modified differential equation. The important nature, i.e. the upwind dominance nature, which is the basis to ensuring that the NP schemes are stable and essentially oscillation free, is firstly presented and verified. Various numerical cases show that the NP schemes are efficient, robust, and more accurate than the original second order central scheme.
[1] Fletcher C A J 1997 Computational Techniques for Fluid Dynamics (Berlin: Springer) [2] Gao Z 2000 Adv. Mech.30 200 (in Chinese) [3] Gao Z 2010 Adv. Mech.40 607 (in Chinese) [4] Shen Y Q, Gao Z and Yang D H 2003 ACTA Aerodynamica Sin.21 342 (in Chinese) [5] Yang M Y, Shu S and Li M J 2010 Chin. J. Hydrodyn.25 307 (in Chinese) [6] Zhu K and Li M J 2011 Chin. J. Theor. Appl. Mech.43 55 (in Chinese) [7] Gao Z and Yang G W 2004 Acta Mech. Sin.20 580 [8] Gao Z, Shen Y Q and Yang G W 2004 CFD J.13 204 [9] Gao Z 2010 Chin. J. Theor. Appl. Mech.42 811 (in Chinese)