A NATURALLY UPWINDED CONSERVATIVE PROCEDURE FOR THE INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON NON-STAGGERED GRIDS

Citation
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
Citations number
33
Categorie Soggetti
Computer Application, Chemistry & Engineering",Mechanics,"Computer Science Interdisciplinary Applications
Journal title
ISSN journal
00457930
Volume
26
Issue
5
Year of publication
1997
Pages
525 - 545
Database
ISI
SICI code
0045-7930(1997)26:5<525:ANUCPF>2.0.ZU;2-Y
Abstract
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.