h‐adaptive finite element solution of high Rayleigh number thermally driven cavity problem

An h‐adaptive finite element code for solving coupled Navier‐Stokes and energy equations is used to solve the thermally driven cavity problem. The buoyancy forces are represented using the Boussinesq approximation. The problem is characterised by very thin boundary layers at high values of Rayleigh number (>106). However, steady state solutions are achievable with adequate discretisation. This is where the auto‐adaptive finite element method provides a powerful means of achieving optimal solutions without having to pre‐define a mesh, which may be either inadequate or too expensive. Steady state and transient results are given for six different Rayleigh numbers in the range 103 to 108 for a Prandtl number of 0.71. The use of h‐adaptivity, based on a posteriori error estimation, is found to ensure a very accurate problem solution at a reasonable computational cost.

[1]  John S. Campbell,et al.  Local and global smoothing of discontinuous finite element functions using a least squares method , 1974 .

[2]  T. J.R. Hughes,et al.  ANALYSIS OF TRANSIENT ALGORITHMS WITH PARTICULAR REFERENCE TO STABILITY BEHAVIOR. , 1983 .

[3]  Asif Usmani,et al.  Finite element modelling of natural convection controlled change of phase , 1992 .

[4]  T. Saitoh,et al.  High-accuracy bench mark solutions to natural convection in a square cavity , 1989 .

[5]  P. LeQuéré,et al.  Accurate solutions to the square thermally driven cavity at high Rayleigh number , 1991 .

[6]  O. C. Zienkiewicz,et al.  Adaptivity and mesh generation , 1991 .

[7]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[8]  S. Paolucci,et al.  Natural convection in an enclosed vertical air layer with large horizontal temperature differences , 1986, Journal of Fluid Mechanics.

[9]  J. Donea.  A Taylor–Galerkin method for convective transport problems , 1983 .

[10]  Jörg Imberger,et al.  Unsteady natural convection in a rectangular cavity , 1980, Journal of Fluid Mechanics.

[11]  O. C. Zienkiewicz,et al.  Error estimation and adaptivity in flow formulation for forming problems , 1988 .

[12]  A. Usmani.  Solution of steady and transient advection problems using an h-adaptive finite element method , 1999 .

[13]  O. Zienkiewicz,et al.  Natural Convection in a Square Enclosure by a Finite-Element, Penalty Function Method Using Primitive Fluid Variables , 1978 .

[14]  P. Le Quéré,et al.  Computation of natural convection in two-dimensional cavities with Chebyshev polynomials , 1985 .

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

[16]  R. Sani,et al.  On the time-dependent solution of the incompressible Navier-Stokes equations in two and three dimensions , 1980 .

[17]  O. C. Zienkiewicz,et al.  A simple error estimator and adaptive procedure for practical engineerng analysis , 1987 .

[18]  S. Paolucci,et al.  Transition to chaos in a differentially heated vertical cavity , 1989, Journal of Fluid Mechanics.

[19]  J. T. Cross,et al.  Finite element analysis of heat transfer and flow problems using adaptive remeshing including application to solidification problems , 1991 .

[20]  G. de Vahl Davis,et al.  Natural convection in a square cavity: A comparison exercise , 1983 .

[21]  R. L. Sani,et al.  Solution of the time-dependent incompressible Navier-Stokes and Boussinesq equations using the Galerkin finite element method , 1980 .

[22]  Robert L. Lee,et al.  Don''t suppress the wiggles|they''re telling you something! Computers and Fluids , 1981 .