Accelerated convergence of the numerical simulation of incompressible flow in general curvilinear co-ordinates by discretizations on overset grids

The convergence rate of a methodology for solving incompressible flow in general curvilinear co-ordinates is analyzed. Overset grids (double-staggered grids type), each defined by the same boundaries as the physical domain are used for discretization. Both grids are Marker and Cell (MAC) quadrilateral meshes with scalar variables (pressure, temperature, etc.) arranged at the center and the Cartesian velocity components at the middle of the sides of the mesh. The problem was checked against benchmark solutions of natural convection in a squeezed cavity, heat transfer in concentric and eccentric horizontal cylindrical annuli and hot cylinder in a duct. Convergence properties of Poisson's pressure equations which arise from application of the SIMPLE-like procedure are analyzed by several methods: successive overrelaxation, symmetric successive overrelaxation, modified incomplete factorization, and conjugate gradient. A genetic algorithm was developed to solve problems of numerical optimization of calculation time, in a space of iteration numbers and relaxation factors. The application provides a means of making an unbiased comparison between the double-staggered grids method and the standard interpolation method. Furthermore, the convergence rate was demonstrated with the well-known calculation of natural convection heat transfer in concentric horizontal cylindrical annuli. Calculation times when double staggered grids were used were 6-10 times shorter than those achieved by interpolation.

[1]  S. Patankar.  Numerical Heat Transfer and Fluid Flow , 2018, Lecture Notes in Mechanical Engineering.

[2]  I. Demirdzic,et al.  Fluid flow and heat transfer test problems for non‐orthogonal grids: Bench‐mark solutions , 1992 .

[3]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[4]  Richard J Goldstein,et al.  An Experimental Study of Natural Convection Heat Transfer in Concentric and Eccentric Horizontal Cylindrical Annuli , 1978 .

[5]  G. D. Raithby,et al.  A method for computing three dimensional flows using non‐orthogonal boundary‐fitted co‐ordinates , 1984 .

[6]  A. Shklyar,et al.  Numerical method for calculation of the incompressible flow in general curvilinear co‐ordinates with double staggered grid , 2003 .

[7]  R. Temam,et al.  Calcul des Ecoulements d'un Fluide Visqueux Incompressibile , 1971 .

[8]  C. W. Hirt,et al.  An Arbitrary Lagrangian-Eulerian Computing Method for All Flow Speeds , 1997 .

[9]  Gene H. Golub,et al.  Matrix computations , 1983 .

[10]  R. P. Kendall,et al.  An Approximate Factorization Procedure for Solving Self-Adjoint Elliptic Difference Equations , 1968 .

[11]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[12]  Wei Shyy,et al.  A numerical study of annular dump diffuser flows , 1985 .

[13]  W. T. Sha,et al.  A SEMI-IMPLICIT CALCULATION PROCEDURE FOR FLOWS DESCRIBED IN BOUNDARY-FITTED COORDINATE SYSTEMS , 1980 .

[14]  A. T. Prata,et al.  Finite-difference solutions of convection-diffusion problems in irregular domains, using a nonorthogonal coordinate transformation , 1984 .

[15]  David M. Young,et al.  Applied Iterative Methods , 2004 .

[16]  G. P. Boerstoel,et al.  The Krylov accelerated SIMPLE(R) method for flow problems in industrial furnaces , 2000 .

[17]  Yankui Liu,et al.  Modeling data envelopment analysis by chance method in hybrid uncertain environments , 2010, Math. Comput. Simul..

[18]  H. L. Stone.  ITERATIVE SOLUTION OF IMPLICIT APPROXIMATIONS OF MULTIDIMENSIONAL PARTIAL DIFFERENTIAL EQUATIONS , 1968 .

[19]  Michael Schäfer,et al.  Efficient discretization of pressure-correction equations on non-orthogonal grids , 2003 .