Application of factor difference scheme to solving discrete flow equations based on unstructured grid

Abstract

A second-order mixing difference scheme with a limiting factor is deduced with the reconstruction gradient method and applied to discretizing the Navier-Stokes equation in an unstructured grid. The transform of nonorthogonal diffusion items generated by the scheme in discrete equations is provided. The Delaunay triangulation method is improved to generate the unstructured grid. The computing program based on the SIMPLE algorithm in an unstructured grid is compiled and used to solve the discrete equations of two types of incompressible viscous flow. The numerical simulation results of the laminar flow driven by lid in cavity and flow behind a cylinder are compared with the theoretical solution and experimental data respectively. In the former case, a good agreement is achieved in the main velocity and drag coefficient curve. In the latter case, the numerical structure and development of vortex under several Reynolds numbers match well with that of the experiment. It is indicated that the factor difference scheme is of higher accuracy, and feasible to be applied to Navier-Stokes equation.

Keywords

unstructured grid / mixing difference scheme / limiting factor / numerical simulation / unsteady viscous flow

References

[1] American Committee of AcademicApplied Mechanics. Research of fluid mechanics: Face towards the demand of nation[J]. Advances in Mechanics, 2006, 36(4): 619-625.

[2] Jameson A., Mavriplis D. Finite volume solution of two dimensional Euler equations on a regular triangular mesh[J]. AIAA J, 1986, 24(4): 611-618.

[3] Mavriplis D J, Jameson A. Multigrid Solution of the Euler Equations on Unstructured and Adaptive Meshes [R]. ICASE Rept. No 87-53, NASA CR-178346, 1988.

[4] Patanker S. V., Spalding D. B. A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows[J]. Int Journal Heat Mass Transfer, 1972, 15 1787-1806.

[5] Taniguchi N., Kobayashi T. Finite volume method on the unstructured grid system[J]. Computer and Fluids, 1991, 19(3): 287-295.

[6] Tao W. The Modern Progress of Computational Heat Transfer[M]. 2000, Beijing: Science Press.

[7] Wang Y. The Study on Generation Techniques of Unstructured Grid and Application in SIMPLE Algorithm [D]. 2005, Tianjin: Tianjin University.

[8] Davidson L. A pressure correction method for unstructured meshes with arbitrary control volumes[J]. International Journal for Numerical Methods in Fluids, 1996, 22(4): 265-281.

[9] Zhang C. Numerical Solution of 3-D Turbulent Flow in Unstructured Grid and the Application Study in Centrifuge Blower[D]. 1999, Xi’an: School of Energy and Power Engineering, Xi’an Jiaotong University.

[10] Ghia U., Chia K. N., Shin C. T. High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method[J]. J of Comp Phys, 1982, 48(3): 387-411.

[11] Bouard R., Coutanceau M. The early stage of development of the wake behind an impulsively started cylinder for 40<Re<104[J]. J of Fluid Mechanics, 1980, 101 975-987.

[12] Zhang Z., Dong Z. Viscous Fluid Mechanics[M]. 1998, Beijing: Tsinghua University Press.

[13] Song G. The PIV Measurement of Unsteady Flow Around Cylinder and Steady Flow in Vaneless Diffuser[D]. 2002, Xi’an: School of Energy and Power Engineering, Xi’an Jiaotong University.