LES-based filter-matrix lattice Boltzmann model for simulating turbulent natural convection in a square cavity

In this paper, a novel thermal filter-matrix lattice Boltzmann model based on large eddy simulation (LES) is proposed for simulating turbulent natural convection. In this study, the Vreman subgrid-scale eddy-viscosity model is introduced into the present framework of LES to accurately predict the flow in near-wall region. Two dimensional numerical simulations of natural convection in a square cavity were performed at high Rayleigh number varying from 10⁷ to 10¹⁰ with a fixed Prandtl number of Pr = 0.71. The influences of the higher-order terms upon the present results at high Rayleigh numbers are examined, taking Ra = 10⁷ and 10⁸ as the example, revealing that the proper minimization of the higher-order terms can improve numerical accuracy of present model for high Rayleigh convective flow. For the turbulent convective flow, the time-averaged quantities in the median lines are presented and compared against those available results from previous studies. The general structure of turbulent boundary layers is well predicted. All numerical results exhibit good agreement with the benchmark solutions available in the previous literatures.     

Turbulent natural convection in vertical parallel-plate channels

The problem of buoyancy driven turbulent flow in parallel-plate channels is investigated. The investigation is limited to vertical channels of uniform cross-section with different modes of heating. The details of the flow and thermal fields are obtained from the solution of the conservation equations of mass, momentum, and energy in addition to equations of the low Reynolds number turbulence model. The study covers Rayleigh number ranging from 10⁵ to 10⁷ and focuses on the effect of channel geometry on the characteristic of the flow and thermal fields as well as the local and average Nusselt number variation. A Nusselt number correlation has been developed in terms of a modified Rayleigh number and channel aspect ratio for the cases of symmetrically heated isothermal and isoflux conditions.     

Laminar natural convection in inclined enclosures bounded by a solid wall

Laminar natural convection has been studied in enclosures bounded by a solid wall with its outer boundary at constant temperature while the opposing side has a constant heat flux. Two-dimensional equations of conservation of mass, momentum and energy, with the Boussinesq approximation are solved using a finite difference method. The numerical procedure adopted is based on the SIMPLER algorithm. Various parameters were: Rayleigh number (from 10³ to 10⁶), dimensionless conductivity of bounding wall (from 1 to 10) and dimensionless wall width (from 0.15 to 0.5), aspect ratio (from 0.5 to 1) and the inclination angle (from 30^° to 180^°). The results are reduced in terms of the normalized Nusselt number as a function of the Rayleigh number, and other dimensionless parameters. The isotherms and streamlines are produced for various Rayleigh numbers and geometrical conditions. It is found that the heat transfer is an increasing function of the Rayleigh number, wall to fluid conductivity ratio, enclosure aspect ratio and a decreasing function of the wall thickness. It passes from a maximum for the inclination angle of about 80^°.     

Time-dependent natural convection in a square cavity: Application of a new finite volume method

A new finite volume (FV) approach with adaptive upwind convection is used to predict the two-dimensional unsteady flow in a square cavity. The fluid is air and natural convection is induced by differentially heated vertical walls. The formulation is made in terms of the vorticity and the integral velocity (induction) law. Biquadratic interpolation formulae are used to approximate the temperature and vorticity fields over the finite volumes, to which the conservation laws are applied in integral form. Image vorticity is used to enforce the zero-penetration condition at the cavity walls. Unsteady predictions are carried sufficiently forward in time to reach a steady state. Results are presented for a Prandtl number (Pr) of 0-71 and Rayleigh numbers equal to 10³, 10⁴ and 10⁵. Both 11 × 11 and 21 × 21 meshes are used. The steady state predictions are compared with published results obtained using a finite difference (FD) scheme for the same values of Pr and Ra and the same meshes, as well as a numerical bench-mark solution. For the most part the FV predictions are closer to the bench-mark solution than are the FD predictions.     

Numerical analysis of free convective flows in partially open enclosure

 Laminar steady state buoyancy induced flows in a two-dimensional, air filled partial open enclosure with a discrete flush mounted iso-flux heat source on one of its walls is investigated numerically. The transport equations for energy and vorticity are solved with the aid of the ADI finite difference scheme on uniform mesh. Because of the specific application of the present study in the air cooling of electronic equipments, results are obtained only for a Prandtl number of 0.71 with an aspect ratio of 1.0 for a range of Rayleigh numbers, Ra (≤10⁵), heat flux parameter, Q and opening parameter, A ₀ using constant properties and Boussinesq approximation by imposing approximate conditions at the opening. Results of flow and temperature patterns, velocity and temperature profiles shows that the outgoing flow is governed by strong characteristics of the cavity condition whereas the incoming flow influenced by outside conditions. It is observed that Rayleigh number considerably affects the flow and thermal fields within the open enclosure when compared with intensity of heat flux and size of the opening. Received on 22 January 2001 / Published online: 29 November 2001     

Numerical study of steady magneto-convection around an adiabatic body inside a square enclosure in low Prandtl numbers

In this paper, the effects of Prandtl number on the steady magneto-convection around a centrally located adiabatic body inside a square enclosure are numerically investigated. Two-dimensional nonlinear governing equations are discretized using the control volume method and hybrid scheme. The equations are solved using SIMPLER algorithm. The results are displayed in the form of streamlines and isotherms when the Rayleigh number varies between 10³ and 10⁶, the Hartmann number changes between 0 and 100 and the Prandtl number ranges between 0.005 and 0.1. The ratio of the buoyancy force to the Lorentz force (Ra/Ha ²) is introduced as an index to compare the contribution of natural convection and magnetic field strength on heat transfer. The results obtained from numerical modeling show that the Prandtl number has not considerable effect on heat transfer at low Rayleigh numbers. The effect of magnetic field strength on convection is increased by increasing Prandtl number. The effect of Prandtl number on the average Nusselt number in the presence of a magnetic field is less than the case without a magnetic field.     

Transient natural convection between concentric and vertically eccentric spheres with mixed boundary conditions

Transient analysis has been investigated numerically to determine heat transfer by natural convection between concentric and vertically eccentric spheres with constant heat flux on the inner wall and a specified isothermal temperature on the outer wall. The governing equations, in terms of vorticity, stream function and temperature are expressed in a spherical polar coordinate system. The alternating direction implicit method and the successive over-relaxation techniques are applied to solve the finite difference form of governing equations. A physical model is introduced which accounts for the effects of fluid buoyancy as well as eccentricity of the outer sphere. Transient solutions of the entire flow field are obtained for a range of modified Rayleigh number (10³<Ra?<5×10⁵), for a Prandtl number of 0.7 and a radius ratio of 2.0, with the outer sphere near the top and bottom of the inner sphere (ε=±0.625). Results of the parametric study conducted further reveal that the heat and flow fields are primarily dependent on the modified Rayleigh number and the eccentricity of the spherical annulus. The results of average Nusselt numbers are also compared with the results obtained for flow between two isothermal spheres.     

Numerical study of simultaneous natural convection heat transfer from both surfaces of a uniformly heated thin plate with arbitrary inclination

 Numerical studies were conducted to investigate the natural convection heat transfer around a uniformly heated thin plate with arbitrary inclination in an infinite space. The numerical approach was based on the finite volume technique with a nonstaggered grid arrangement. For handling the pressure–velocity coupling the SIMPLE-algorithm was used. QUICK scheme and first order upwind scheme were employed for discretization of the momentum and energy convective terms respectively. Plate width and heating rate were used to vary the modified Rayleigh number over the range of 4.8×10⁶ to 1.87×10⁸. Local and average heat transfer characteristics were compared with regarding to the inclination angle. The empirical expressions for local and average Nusselt number were correlated. It has been found that for inclination angle less than 10^∘, the flow and heat transfer characteristics are complicated and the average Nusselt number can not be correlated by one equation while for inclination angle larger than 10^∘, the average Nusselt number can be correlated into an elegant correlation. Received on 18 April 2001 / Published online: 29 November 2001     

Analytical and numerical study of natural convection heat transfer in an inclined composite enclosure

Laminar natural convection heat transfer in inclined fluid layers divided by a partition with finite thickness and conductivity is studied analytically and numerically. The governing equations for the fluid layers are solved analytically in the limit of a thin layered system with constant flux boundary conditions. The study covers of the range of Ra from 10³ to 10⁷, from 0° to 180° and the thermal conductivity ratio of partition to fluid ratioK from 10^–2 to 10⁶. The Prandtl number was 0.72 (for air). Results are obtained in terms of an overall Nusselt number as a function of Rayleigh number, angle of inclination of the system, mid layer thickness, and mid layer thermal conductivity. The critical Rayleigh number for the onset of convection in a bottom-heated horizontal system is predicted. The results are compared with the numerical results obtained by solving the complete system of governing equations, using SIMPLER method, as well as with the limiting cases in the literature.     

The numerical solution of the unsteady natural convection flow in a square cavity at high Rayleigh number using SADI method

The unsteady natural convection flow in a square cavity at high Rayleigh number Ra=10 ⁷ and 2×10 ⁷ has been computed using cubic spline integration. The required solutions to the two dimensional Navier-Stokes and energy equations have been obtained using two alternate numerical formulations on non-uniform grids. The main features of the transient flow have been briefly discussed. The results obtained by using the present method are in good agreement with the theoretical predictions ^[1,2].The steady state results have been compared with accurate solutions presented recently for Ra=10 ⁷.     

Natural convection with surface radiation from a planar heat generating element mounted freely in a vertical channel

Experiments and numerical simulations have been conducted to study the conjugate heat transfer by natural convection and surface radiation from a planar heat generating element placed centrally between two adiabatic vertical plates. The relevant problem dependent parameters considered in this study are modified Rayleigh number, channel aspect ratio, stream-wise location of the heat generating element, and surface emissivities of the heat generating element and the adiabatic side plates. Experiments are conducted for different values of modified Rayleigh number ranging from 3.2 × 10⁵ to 1.6 × 10⁷ and surface emissivities 0.05, 0.55, 0.75 and 0.85. The interdependence between the heat transfer mechanism and the flow field under the influence of surface radiation on natural convection is explored and discussed. Experimental correlations for total and convective Nusselt number, and dimensionless temperature in terms of relevant parameters have been developed. The mathematical model governing the problem has been numerically solved using a commercial computational fluid dynamics package FLUENT 6.3 and the numerical predictions substantiate the experimental observations.     

A least-squares finite element method for time-dependent incompressible flows with thermal convection

The time-dependent Navier–Stokes equations and the energy balance equation for an incompressible, constant property fluid in the Boussinesq approximation are solved by a least-squares finite element method based on a velocity–pressure–vorticity–temperature–heat-flux ( u –P–ω–T– q ) formulation discretized by backward finite differencing in time. The discretization scheme leads to the minimization of the residual in the l²-norm for each time step. Isoparametric bilinear quadrilateral elements and reduced integration are employed. Three examples, thermally driven cavity flow at Rayleigh numbers up to 10⁶, lid-driven cavity flow at Reynolds numbers up to 10⁴ and flow over a square obstacle at Reynolds number 200, are presented to validate the method.     

Non-Darcy buoyancy driven flows in a fluid saturated porous medium: the use of asymptotic computational fluid dynamics (ACFD) approach

Natural convection in a fluid saturated porous medium has been numerically investigated using a generalized non-Darcy approach. The governing equations are solved by using Finite Volume approach. First order upwind scheme is employed for convective formulation and SIMPLE algorithm for pressure velocity coupling. Numerical results are presented to study the influence of parameters such as Rayleigh number (10⁶ ≤Ra ≤10⁸), Darcy number (10⁻⁵ ≤ Da ≤ 10⁻²), porosity (0.4 ≤ ɛ ≤ 0.9) and Prandtl number (0.01 ≤ Pr ≤ 10) on the flow behavior and heat transfer. By combining the method of matched asymptotic expansions with computational fluid dynamics (CFD), so called asymptotic computational fluid dynamics (ACFD) technique has been employed to generate correlation for average Nusselt number. The technique is found to be an attractive option for generating correlation and also in the analysis of natural convection in porous medium over a fairly wide range of parameters with fewer simulations for numerical solutions.     

Laminar natural convection in a partially divided rectangular cavity at high Rayleigh number

Finite element predictions of two-dimensional laminar natural convection in a partially divided rectangular cavity at high Rayleigh number are presented. The walls are differentially heated, the horizontal surfaces are insulated and the cavity contains a partial vertical divider which is centrally located and whose height is varied. Detailed results are presented for an aluminium half-divider in water for Rayleigh number up to 10¹¹ and compared directly with recent experiments in a cavity of aspect ratio 1/2. The predicted flow and heat transfer are in good agreement with the measurements and confirm the existence of a high Rayleigh number regime with characteristic behaviour that differs significantly from that found at lower Rayleigh number. In addition, the effects of the divider height, the divider conductivity, the fluid Prandtl number and the cavity aspect ratio are studied. The results show that a direct simulation of the complex flow and heat transfer that occurs in partially divided cavities is possible for realistic physical conditions.     

Thermal effects of the micropolar fluid parameters on the laminar free convection in spherical annuli with mixed boundary conditions

The effects of micro-rotation and vortex viscosity in micropolar fluids have been investigated numerically to determine heat transfer by natural convection between concentric and vertically eccentric spheres with specified mixed boundary conditions. Calculations were carried out systematically for several different eccentricities and a range of modified Rayleigh numbers to determine the average Nusslet numbers which are affected by the micropolar parameters (F) of the flow and temperature fields. The skin friction stress on the walls has also been studied and discussed. The governing equations, in terms of vorticity, stream function, temperature and angular momentum are expressed in a spherical polar coordinate system. Results were obtained for steady heat-transfer in spherical annuli at a Prandtl number of 0.7, with the modified Rayleigh number ranging from 10³ to 5 × 10⁵, for a radius ratio of 2.0 and eccentricities varying from −0.625 to +0.625. Comparisons are attempted between the Newtonian fluid and micropolar fluid.     

Free convection from one thermal and solute source in a confined porous medium

This paper reports a numerical study of double diffusive natural convection in a vertical porous enclosure with localized heating and salting from one side. The physical model for the momentum conservation equation makes use of the Darcy equation, and the set of coupled equations is solved using the finite-volume methodology together with the deferred central difference scheme. An extensive series of numerical simulations is conducted in the range of −10 ⩽ N ⩽ + 10, 0 ⩽ R _t ⩽ 200, 10⁻² ⩽ Le ⩽ 200, and 0.125 ⩽ L ⩽ 0.875, where N, R _t , Le, and L are the buoyancy ratio, Darcy-modified thermal Rayleigh number, Lewis number, and the segment location. Streamlines, heatlines, masslines, isotherms, and iso-concentrations are produced for several segment locations to illustrate the flow structure transition from solutal-dominated opposing to thermal dominated and solutal-dominated aiding flows, respectively. The segment location combining with thermal Rayleigh number and Lewis number is found to influence the buoyancy ratio at which flow transition and flow reversal occurs. The computed average Nusselt and Sherwood numbers provide guidance for locating the heating and salting segment.     

Thermo-Bioconvection in a Square Porous Cavity Filled by Oxytactic Microorganisms

This paper studies the thermo-bioconvection in a square porous cavity filled by oxytactic microorganisms. The Darcy model with Boussinesq approximation has been used to solve the flow and heat and mass transfer in the porous region. The governing equations formulated in terms of the dimensionless stream function, temperature and concentration have been solved using the finite difference method. Comparison with results from the open literature of the mean Nusselt number for a square cavity filled with a regular porous medium is made. It is shown that the results are in very good agreement. The main objective was to investigate the influence of the traditional Rayleigh number Ra = 10, 100, bioconvection Rayleigh number Rb = 10, 100, Lewis number Le = 1, 10, and Péclet number Pe = 0.1, 1 on the fluid flow and heat and mass transfer. Comprehensive analysis of an effect of these key parameters on the Nusselt and Sherwood numbers at the vertical walls has been conducted.     

Natural convection in a rectangular cavity with wall temperature decreasing at a uniform rate

The transient natural convection in a fluid contained in a rectangular enclosure, the wall of which is maintained at a uniform temperature which changes at a steady rate, is approached by a numerical method. Numerical solutions are obtained forPr=0.73, 7.3 and 73 and a range of Rayleigh numbersRa=10² ～ 10⁸. At relatively low Rayleigh numbers the flow is characterized by the development of double cells with flow up the center and down the sidewalk However it was found that an increase of the Rayleigh number leads to the development of strong secondary circulation on the axis of symmetry of the cavity near the top wall. Thus, as the Rayleigh number is increased the secondary cells grow in size. The effects of the secondary cells on the temperature field and heat transfer coefficients are discussed. Most results are obtained for the case of a square cavity (E=2) but the influence of the aspect ratio of the cavity is also studied forE=1 and 4.     

A numerical study of transient natural convection under time-dependent gravitational acceleration field

A transient natural convection caused by abruptly heating of a vertical wall in a square enclosure which is under time-dependent gravitational acceleration field is studied numerically. During the computing process, a penalty finite element method is adopted to solve the governing equations. The factor of the time-dependent gravitational acceleration is included in Rayleigh number, therefore two variation of Rayleigh number are considered, one is from 10⁵ via 10⁶ to ? 10⁵ and the other is from 10⁴ via 10⁵ to ? 10⁴. Since the variations of flow field and temperature distribution of the fluid in the square enclosure are complicated and unpredictable, the isotherms, streamlines and time dependent heat transfer rate are examined for investigating the heat transfer mechanism during the transient process. The results show that the agreement of the heat transfer rate of this study with that of quasi-steady state is mainly dependent on the response of the variation rate of boundary layer along the vertical walls to that of time-dependent gravitational acceleration.     

Role of the inclination of an inverted‐V upper plate on the heat and flow behavior of trapped gases in a modified Rayleigh–Bénard cavity

Thermal buoyant air inside a modified Rayleigh–Bénard (RB) cavity bounded by a lower flat plate and an inverted‐V upper plate has been investigated numerically using the finite‐volume method. The second‐order‐accurate QUICK and SIMPLE schemes were used for the discretization of the convective terms and the pressure–velocity coupling in the set of conservation equations, respectively. The problem under study is controlled by two parameters: (1) the Rayleigh number ranging from 10³ to 10⁶ and (2) the relative height of the vertical sidewalls d. In reference to the latter, it varies from one limiting case corresponding to the standard RB cavity (a rectangle with d = 1) to another limiting case represented by an isosceles triangular cavity where d = 0. The numerical results for the velocity and temperature fields are presented in terms of streamlines, isotherms, local and mean heat fluxes. An additional effort was devoted to determine the critical Ra values characterizing the transition from symmetrical to asymmetrical buoyant airflow responsive to incremental changes in Ra. For purposes of engineering design, a general correlation equation for the Nusselt number in terms of the pertinent Ra and d was constructed using nonlinear multiple regression theory. Copyright © 2007 John Wiley & Sons, Ltd.