GDQ METHOD FOR NATURAL CONVECTION IN A SQUARE CAVITY USING VELOCITY–VORTICITY FORMULATION

ABSTRACT This article describes a compact numerical algorithm based on the generalized differential quadrature (GDQ) method for the numerical analysis of natural convection in a differentially heated square cavity. The velocity–vorticity form of the Navier–Stokes equations and energy equation are used to represent the mass, momentum, and energy conservations of the fluid medium in the cavity. The GDQ form of the governing equations and the vorticity definition at the boundaries are solved by a coupled solution algorithm using a global matrix scheme for all the field variables. The vorticity values at the boundary are correctly imposed using the GDQ method, which approximates a given space derivative with higher-order accuracy compared to the existing schemes based on Taylor's series expansion. This has assured a divergence-free solution for the flow field by satisfying the continuity constraint, though the pressure term is not used directly in the present formulation. The proposed algorithm is validated for a lid-driven cavity flow for Reynolds number Re = 100, 400, and 1,000, and the predicted velocity profiles are in excellent agreement with the benchmark solutions. The algorithm is then used to compute the average Nusselt number and flow parameters for natural convection in a square cavity for Rayleigh number Ra = 103, 104, 105, and 106. These results are in better agreement with the benchmark solutions than the results obtained by other numerical schemes, which used much finer grids compared to the present scheme.

[1]  William H. Press,et al.  Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing , 1989 .

[2]  G. D. Davis.  Natural convection of air in a square cavity: A bench mark numerical solution , 1983 .

[3]  Takeo S. Saitoh,et al.  Benchmark solutions for natural convection in a cubic cavity using the high-order time–space method , 2004 .

[4]  D. Young,et al.  Two-Dimensional Incompressible Flows by Velocity-Vorticity Formulation and Finite Element Method , 2001 .

[5]  Marc S. Ingber,et al.  A Galerkin implementation of the generalized Helmholtz decomposition of vorticity formulations , 2001 .

[6]  Wagdi G. Habashi,et al.  Finite Element Solution of the 3D Compressible Navier-Stokes Equations by a Velocity-Vorticity Method , 1993 .

[7]  Wen-Zhong Shen,et al.  Numerical method for unsteady 3D Navier-Stokes equations in velocity-vorticity form , 1997 .

[8]  Numerical simulation of natural convection in a square cavity by SIMPLE-generalized differential quadrature method , 2002 .

[9]  Chang Shu,et al.  AN EFFICIENT APPROACH TO SIMULATE NATURAL CONVECTION IN ARBITRARILY ECCENTRIC ANNULI BY VORTICITY-STREAM FUNCTION FORMULATION , 2000 .

[10]  M. Napolitano,et al.  A multigrid solver for the vorticity‐velocity Navier‐Stokes equations , 1991 .

[11]  Hermann F. Fasel,et al.  Investigation of the stability of boundary layers by a finite-difference model of the Navier—Stokes equations , 1976, Journal of Fluid Mechanics.

[12]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[13]  O. Daube.  Resolution of the 2D Navier-Stokes Equations in Velocity-Vorticity Form by Means of an Influence Matrix Technique , 1992 .

[14]  Chang Shu,et al.  A 3D incompressible thermal lattice Boltzmann model and its application to simulate natural convection in a cubic cavity , 2004 .

[15]  C. Davies,et al.  A novel velocity-vorticity formulation of the Navier-Stokes equations with applications to boundary layer disturbance evolution , 2001 .

[16]  L. Quartapelle,et al.  A review of vorticity conditions in the numerical solution of the ζ–ψ equations , 1999 .

[17]  A. J. Baker,et al.  A 3D incompressible Navier–Stokes velocity–vorticity weak form finite element algorithm , 2002 .

[18]  Bakhtier Farouk,et al.  A numerical study of three-dimensional natural convection in a differentially heated cubical enclosure , 1991 .

[19]  D. L. Young,et al.  Arbitrary Lagrangian-Eulerian finite element analysis of free surface flow using a velocity-vorticity formulation , 2004 .

[20]  R. Sani,et al.  On pressure boundary conditions for the incompressible Navier‐Stokes equations , 1987 .

[21]  G. Mallinson,et al.  Three-dimensional natural convection in a box: a numerical study , 1977, Journal of Fluid Mechanics.

[22]  C. Shu.  Differential Quadrature and Its Application in Engineering , 2000 .

[23]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[24]  Leopold Škerget,et al.  Mixed boundary elements for laminar flows , 1999 .

[25]  Chang Shu,et al.  Numerical solutions of incompressible Navier—Stokes equations by generalized differential quadrature , 1994 .

[26]  F. Stella,et al.  A vorticity-velocity method for the numerical solution of 3D incompressible flows , 1993 .

[27]  Chang Shu,et al.  Numerical computation of three‐dimensional incompressible Navier–Stokes equations in primitive variable form by DQ method , 2003 .

[28]  Robert B. Kinney,et al.  Time‐dependent natural convection in a square cavity: Application of a new finite volume method , 1990 .

[29]  J. C. Kalita,et al.  Fully compact higher-order computation of steady-state natural convection in a square cavity. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  G. Labrosse,et al.  A FIRST INCURSION INTO THE 3D STRUCTURE OF NATURAL CONVECTION OF AIR IN A DIFFERENTIALLY HEATED CUBIC CAVITY, FROM ACCURATE NUMERICAL SOLUTIONS , 2000 .

[31]  W. Habashi,et al.  Finite element solution of the navier-stokes equations by a velocity-vorticity method , 1990 .

[32]  R. Bellman,et al.  DIFFERENTIAL QUADRATURE: A TECHNIQUE FOR THE RAPID SOLUTION OF NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS , 1972 .

[33]  C. Shu,et al.  APPLICATION OF GENERALIZED DIFFERENTIAL QUADRATURE TO SOLVE TWO-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS , 1992 .