A comparative study of central and upwind difference schemes using the primitive variables

The use of the velocity-pressure formulation of the Navier-Stokes equations for the numerical solution of fluid flow problems is favoured for free-surface problems, more involved flow configurations, and three-dimensional flows. Many engineering problems involve such features in addition to strong inertial effects. The computational instabilities arising from central-difference schemes for the convective terms of the governing equations impose serious limitations on the range of Reynolds numbers that can be investigated by the numerical method. Solutions for higher Reynolds numbers Re > 1000 could be reached using upwind-difference schemes. A comparative study of both schemes using a method based on the primitive variables is presented. The comparison is made for the model problem of the driven flow in a square cavity. Using a central scheme stable solutions of the pressure and velocity fields were obtained for Reynolds numbers up to 5000. The streamfunction and vorticity fields were calculated from the resulting velocity field and compared with previous solutions. It is concluded that total upwind differencing results in a considerable change in the flow pattern due to the false diffusion. For practical calculations, by a proper choice of a small amount of partial upwind differencing the vorticity diffusion near the walls and the global features of the solutions are not sigificantly altered.