A New Numerical Scheme for Convection-Dominated Transport Equations.

Abstract

In order to obtain numerical solutions stable and accurate for convection-dominated transport equations, we propose a criterion in constructing numerical schemes for the convection term that roots of the characteristic equation for resulting difference equation have poles. By imposing this criterion on the difference coefficients for the convection term, we construct a new numerical scheme robust for convection-dominated equations. The present new scheme coincides with the QUICK scheme when the mesh Reynolds number (Rm) is 8/3, which is the critical value for its stability, while it approaches the second-order upwind scheme as Rm goes to infinity. Hence the present scheme interpolates a stable scheme between the QUICK scheme at Rm=8/3 and the second-order upwind scheme at Rm=infinity. This new scheme shows good numerical solutions for one-dimensional, convection-diffusion equations.