Wh. Calhoon et Rl. Roach, A NATURALLY UPWINDED CONSERVATIVE PROCEDURE FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON NON-STAGGERED GRIDS, Computers & fluids, 26(5), 1997, pp. 525-545
Upwind methods have been previously adopted in incompressible flow cal
culations in an attempt to eliminate the need for artificial viscosity
while maintaining high accuracy. Standard upwind methods found in the
literature are primarily concerned with accurate modeling of the conv
ection terms and generally neglect the viscous terms in the analysis.
A new procedure is developed here which seeks to include the viscous t
erms so that regions of significant gradients will not be over-dissipa
ted. The approach taken here is an extension of the work of Refs [1-5]
in which interpolating functions are derived from direct integration
of linearized forms of the governing equations. The resulting function
s are then used as interpolants for differencing the full equations. D
uring the process, no attempt is made to deliberately upwind bias the
differencing. However, because of the manner in which coefficients of
the differencing are functions of the cell Reynolds number, the scheme
introduces its own upwind bias. This property was first observed by R
oscoe [1] in which the observation that this class of schemes preserve
s the hyperbolic/elliptic nature of the original system of equations w
as made. The present work yields, for the first time, a fully conserva
tive Roscoe-type method. Although the method approaches the developmen
t of the finite difference equations in a different manner than tradit
ional upwind procedures, it achieves higher accuracy in a manner analo
gous to the MUSCL approach [6] used in standard upwind schemes. The pr
ocedure also uses a unique discretization of the continuity equation w
hich involves pressure terms as in Ref. [5]. This equation is fully co
nservative and satisfies mass conservation to machine zero. The new sc
heme is first demonstrated on the 1-D, viscous, Burger's equation and
results are compared with a first and second-order Roe scheme. It is t
hen applied to the 2-D, incompressible Navier-Stokes equations for non
-staggered grids. Results for the driven cavity problem are compared w
ith proven methods and are found to be in excellent agreement. (C) 199
7 Elsevier Science Ltd.