Numerical investigation of natural convection of Al<Subscript>2</Subscript>O<Subscript>3</Subscript> nanofluid in vertical annuli

This paper presents the numerical study of internal free convection of Al₂O₃ water nanofluid in vertical annuli. Vertical walls are maintained at constant temperatures and horizontal walls are adiabatic. Results are validated by experimental data. Effect of nanofluids on natural convection is investigated as a function of geometrical and physical parameters and particle fractions for aspect ratio of 1 ≤ H/L ≤ 5, Grashof number of 10³ ≤ Gr ≤ 10⁵ and concentration of 0 ≤ ϕ ≤ 0.06. More than 330 different numerical cases are investigated to develop a new correlation for the Nusselt number. This correlation is presented as a function of Nusselt number of base fluid and particle fraction which is a linear decreasing function of particle fraction. The developed correlation for annuli is also valid for the natural convection of Al₂O₃ water nanofluid in a square cavity. Furthermore, the effect of the viscosity and conductivity models on the Nusselt number of nanofluids in cylindrical cavities are discussed.     

Convection in a Lid-Driven Heat-Generating Porous Cavity with Alternative Thermal Boundary Conditions

Mixed convection flow in a two-sided lid-driven cavity filled with heat-generating porous medium is numerically investigated. The top and bottom walls are moving in opposite directions at different temperatures, while the side vertical walls are considered adiabatic. The governing equations are solved using the finite-volume method with the SIMPLE algorithm. The numerical procedure adopted in this study yields a consistent performance over a wide range of parameters that were 10⁻⁴ ≤ Da ≤ 10⁻¹ and 0 ≤ Ra _I ≤ 10⁴. The effects of the parameters involved on the heat transfer characteristics are studied in detail. It is found that the variation of the average Nusselt number is non-linear for increasing values of the Darcy number with uniform or non-uniform heating condition.     

Natural convection in inclined two dimensional rectangular cavities

Steady two-dimensional natural convection in fluid filled cavities is numerically investigated. The channel is heated from below and cooled from the top with insulated side walls and the inclination angle is varied. The field equations for a Newtonian Boussinesq fluid are solved numerically for three cavity height based Rayleigh numbers, Ra = 10⁴, 10⁵ and 10⁶, and several aspect ratios. The calculations are in excellent agreement with previously published benchmark results. The effect of the inclination of the cavity to the horizontal with the angle varying from 0° to 180° and the effect of the startup conditions on the flow pattern, temperature distribution and the heat transfer rates have been investigated. Flow admits different configurations at different angles as the angle of inclination is increased depending on the initial conditions. Regardless of the initial conditions Nusselt number Nu exhibits discontinuities triggered by gradual transition from multiple cell to a single cell configuration. The critical angle of inclination at which the discontinuity occurs is strongly influenced by the assumed startup field. The hysteresis effect previously reported is not always present when the calculations are reversed from 90° to 0°. A comprehensive study of the flow structure, the Nu variation with varying angle of inclination, the effect of the initial conditions and the hysteresis effect are presented.     

Heat transfer by laminar flow in a vertical and horizontal pipe with twisted-tape inserts

 This article presents the results of laboratory research on heat exchange while heating water in horizontal and vertical tubes with twisted-tape inserts. The scope of the research: 70 ≤ Re ≤ 4000 3.6 ≤ Pr ≤ 5.9 8.6 ≤ Gz ≤ 540 The research was held for three cases: – horizontal experimental tube – vertical experimental tube, the direction of flow according to the free convection vector – vertical experimental tube, the direction of flow not in accordance with the free convection vector For such cases the correlation equation was defined Nu_T=f(Gz; y), Nu = f(Gz) and the proportion Nu_T/Nu was analysed. Received on 30 March 2000     

Average and extremal properties of heat transfer and shear stress on a wall surface in Rayleigh–Bénard convection

In this study surface-averaged and extremal properties of heat transfer and shear stress on the upper wall surface of Rayleigh–Bénard convection are numerically examined. The Prandtl number was raised up to 10³, and the Rayleigh number was changed between 10⁴ and 10⁷. As a result, average Nusselt number Nu and shear rate τ/Pr depends on Pr, Ra, and the entire numerical results are distributed between two correlation equations corresponding to small and large Pr. The small and large Pr equations are closely related to steady and unsteady flow regimes, respectively. Nevertheless, a single relation τ/Pr ~ Nu ^3.0 exists to explain the entire results. Similarly the change of local maximal properties Nu _max and τ _max/Pr depends on Pr, Ra, and these values are also distributed between two correlation equations corresponding to small and large Pr cases. Despite such complicated dependence we can obtain a correlation equation as a form of τ _max/Pr ~ Nu _max^2.6, which has not been obtained theoretically.     

Non-Darcian and Anisotropic Effects on Free Convection in a Porous Enclosure

The free convective flow and heat transfer, within the framework of Boussinesq approximation, in an anisotropic fluid filled porous rectangular enclosure subjected to end-to-end temperature difference have been investigated using Brinkman extended non-Darcy flow model. The studies involve simultaneous consideration of hydrodynamic and thermal anisotropy. The flow and temperature fields in general are governed by, Ra, the Rayleigh number, AR, the aspect ratio of the slab, K*, the permeability ratio and k*, the thermal conductivity ratio, and Da, Darcy number. Numerical solutions employing the successive accelerated replacement (SAR) scheme have been obtained for 100 ≤ Ra ≤ 1000, 0.5 ≤ AR ≤ 5, 0.5 ≤ K* ≤ 5, 0.5 ≤ k* ≤ 5, and 0 ≤ Da ≤ 0.1. It has been found that [`(Nu)]{\overline {Nu}}, average Nusselt number increases with increase in K* and decreases as k* increases. However, the magnitude of the change in [`(Nu)]{\overline {Nu}} depends on the parameter Da, characterizing the Brinkman extended non-Darcy flow.     

Transient natural convection cooling of a high Prandtl number fluid in a cubical cavity

This work presents a numerical analysis of the effects of thermal boundary conditions, fluid variable viscosity and wall conduction on transient laminar natural convection of a high Prandtl number (Pr=4×10⁴) fluid (Golden Syrup) in a cubical cavity. The simulations consider physical situations realizable at laboratory scale using a cavity with Plexiglas walls of 1 cm of thickness, and inside dimension of L=20 cm. The initial Rayleigh (Ra) number is 10⁶. The cavity is initially full of fluid at rest and at constant temperature (T _i =45°C) higher than the temperature of the walls (T _w =25°C). The time evolution of the flow patterns, the temperature contours, the mean temperature of the fluid and the Nusselt number (Nu) of eight different cases of cooling are presented and analyzed.     

Characterization of fluid flow patterns and heat transfer in horizontal channel mixed convection

Two mechanisms of roll initiation are highlighted in a horizontal channel flow, uniformly heated from below, at constant heat flux (Γ = 10, Pr = 7, 50 ≤ Re ≤ 100, 0 ≤ Ra ≤ 10⁶). The first mechanism is the classical one, it occurs for low Rayleigh numbers and is initiated by the lateral wall effect. The second occurs for higher Rayleigh numbers and combines the previous effect with a supercritical vertical temperature gradient in the lower boundary layer, which simultaneously triggers pairs of rolls in the whole zone in between the two lateral rolls. We have found that in the present configuration, the transition between the two roll initiation mechanisms occurs for Ra/Re ² ≈ 18. Consequently, the heat transfer is significantly enhanced compared to the pure forced convection case owing to the flow pattern responsible of the continuous flooding the heated wall with cold fluid.     

Experimental and numerical study on laminar natural convection in a cavity heated from bottom due to an inclined fin

Natural convection heat transfer in an inclined fin attached square enclosure is studied both experimentally and numerically. Bottom wall of enclosure has higher temperature than that of top wall while vertical walls are adiabatic. Inclined fin has also adiabatic boundary conditions. Numerical solutions have been done by writing a computer code in Fortran platform and results are compared with Fluent commercial code and experimental method. Governing parameters are Rayleigh numbers (8.10⁵ ≤ Ra ≤ 4 × 10⁶) and inclination angle (30° ≤ and ≤ 120°). The temperature measurements are done by using thermocouples distributed uniformly at the wall of the enclosure. Remarkably good agreement is obtained between the predicted results and experimental data. A correlation is also developed including all effective parameters on heat transfer and fluid flow. It was observed that heat transfer can be controlled by attaching an inclined fin onto wall.     

Effect of the surface thermal radiation on turbulent natural convection in tall cavities of façade elements

The effect of the surface thermal radiation in tall cavities with turbulent natural convection regime was analyzed and quantified numerically. The parameters considered were: the Rayleigh number 10⁹–10¹², the aspect ratio 20, 40 and 80 and the emmisivity 0.0–1.0. The percentage contribution of the radiative surface to the total heat transfer has a maximum value of  15.19% (Ra = 10⁹, A = 20) with emissivity equal to 1.0 and a minimum of 0.5% (Ra = 10¹², A = 80) with ε* = 0.2. The average radiative Nusselt number for a fixed emissivity is independent of the Rayleigh number, but for a fixed Rayleigh number diminishes with the increase of the aspect ratio. The results indicate that the surface thermal radiation does not modify significantly the flow pattern in the cavity, just negligible effects in the bottom and top of the cavity were observed. Two different temperature patterns were observed a conductive regime Ra = 10⁹ and a boundary layer regime Ra = 10¹².     

Buoyancy and cross flow effects on heat transfer of multiple impinging slot air jets cooling a flat plate at different orientations

The present article reports on heat transfer characteristics associated with multiple laminar impinging air jet cooling a hot flat plat at different orientations. The work aims to study the interactions of the effects of cross flow, buoyancy induced flow, orientation of the hot surface with respect to gravity, Reynolds numbers and Rayleigh numbers on heat transfer characteristics. Experiments have been carried out for different values of jet Reynolds number, Rayleigh number and cross flow strength and at different orientations of the air jet with respect to the target hot plate. In general, the effective cooling of the plate has been observed to be increased with increasing Reynolds number and Rayleigh number. The results concluded that the hot surface orientation is important for optimum performance in practical applications. It was found that for Re ≥ 400 and Ra ≥ 10,000 (these ranges give 0.0142 ≤ Ri ≤ 1.59 the Nusselt number is independent on the hot surface orientation. However, for Re ≤ 300 and Ra ≥ 100,000 (these ranges give 1.59 ≤ Ri ≤ 42.85): (i) the Nusselt number for horizontal orientation with hot surface facing down is less that that of vertical orientation and that of horizontal orientation with hot surface facing up, and (ii) the Nusselt number of vertical orientation is approximately the same as that of horizontal orientation with hot surface facing up. For all surfaces orientations and for the entire ranges of Re and Ra, it was found that increasing the cross flow strength decreases the effective cooling of the surface.     

Experimental investigation of natural convection heat transfer in horizontal and inclined annular fluid layers

In the present study, an experimental investigation of heat transfer and fluid flow characteristics of buoyancy-driven flow in horizontal and inclined annuli bounded by concentric tubes has been carried out. The annulus inner surface is maintained at high temperature by applying heat flux to the inner tube while the annulus outer surface is maintained at low temperature by circulating cooling water at high mass flow rate around the outer tube. The experiments were carried out at a wide range of Rayleigh number (5 × 10⁴ < Ra < 5 × 10⁵) for different annulus gap widths (L/D _o = 0.23, 0.3, and 0.37) and different inclination of the annulus (α = 0°, 30° and 60°). The results showed that: (1) increasing the annulus gap width strongly increases the heat transfer rate, (2) the heat transfer rate slightly decreases with increasing the inclination of the annulus from the horizontal, and (3) increasing Ra increases the heat transfer rate for any L/D _o and at any inclination. Correlations of the heat transfer enhancement due to buoyancy driven flow in an annulus has been developed in terms of Ra, L/D _o and α. The prediction of the correlation has been compared with the present and previous data and fair agreement was found.     

Secondary flow patterns and mixing in laminar pulsating flow through a curved pipe

Mixing by secondary flow is studied by particle image velocimetry (PIV) in a developing laminar pulsating flow through a circular curved pipe. The pipe curvature ratio is η = r ₀/r _c  = 0.09, and the curvature angle is 90°. Different secondary flow patterns are formed during an oscillation period due to competition among the centrifugal, inertial, and viscous forces. These different secondary-flow structures lead to different transverse-mixing schemes in the flow. Here, transverse mixing enhancement is investigated by imposing different pulsating conditions (Dean number, velocity ratio, and frequency parameter); favorable pulsating conditions for mixing are introduced. To obviate light-refraction effects during PIV measurements, a T-shaped structure is installed downstream of the curved pipe. Experiments are carried out for the Reynolds numbers range 420 ≤ Re_st ≤ 1,000 (Dean numbers 126.6 ≤ Dn ≤ 301.5) corresponding to non-oscillating flow, velocity component ratios 1 ≤ (β = U _max,osc/U _m,st) ≤ 4 (the ratio of velocity amplitude of oscillations to the mean velocity without oscillations), and frequency parameters 8.37 < (α = r ₀(ω/ν)^0.5) < 24.5, where α² is the ratio of viscous diffusion time over the pipe radius to the characteristic oscillation time. The variations in cross-sectional average values of absolute axial vorticity (|ζ|) and transverse strain rate (|ε|) are analyzed in order to quantify mixing. The effects of each parameter (Re_st, β, and α) on transverse mixing are discussed by comparing the dimensionless vorticities (|ζ _P |/|ζ _S |) and dimensionless transverse strain rates (|ε _P |/|ε _S |) during a complete oscillation period.     

Combined Laminar Mixed Convection and Surface Radiation using Asymptotic Computational Fluid Dynamics (ACFD)

This paper reports the use of the technique of combining asymptotics with computational fluid dynamics (CFD), known as asymptotic computational fluid dynamics (ACFD), to handle the problem of combined laminar mixed convection and surface radiation from a two dimensional, differentially heated lid driven cavity. The fluid under consideration is air, which is radiatively transparent, and all the walls are assumed to be gray and diffuse and having the same hemispherical, total emissivity (ɛ). The computations have been performed on FLUENT 6.2. The full radiation problem (i.e. all the walls are radiatively black corresponding to ɛ = 1) is first taken up and the method of “perturbing and blending” is used wherein, first, limiting solutions of natural and forced convection are perturbed, to obtain correlations for the weighted average convective Nusselt numbers for the full radiation case. These correlations are then blended suitably in order to obtain a composite correlation for the weighted average convective Nusselt number that is valid for the entire mixed convection range, i.e., 0 ≤ Ri ≤ ∞. This correlation is then expanded in terms of ɛ to obtain an expression for the average convective Nusselt number that is valid for any ɛ in the range 0 ≤ ɛ ≤ 1. In so far as radiation heat transfer is concerned, using asymptotic arguments, a new weighted average radiation Nusselt number is defined such that this quantity can be expanded just in terms of ɛ. Hence, by the use of ACFD, the number of solutions required to obtain reasonably accurate correlations for both the convective and radiative heat transfer rates and hence the total heat transfer rate (Nu _total = Nu _C + Nu _R), is substantially reduced. More importantly, the correlations for convection and radiation are asymptotically correct at their ends. The effect of secondary variables like aspect ratio and the case of unequal wall emissivities can also be included without significant additional effort.     

Experimental laminar Rayleigh-Bénard convection in a cubical cavity at moderate Rayleigh and Prandtl numbers

Rayleigh-Bénard convection in a cubical cavity with adiabatic or conductive sidewalls is experimentally analyzed at moderate Rayleigh numbers (Ra ≤ 8 × 10⁴) using silicone oil (Pr=130) as the convecting fluid. Under these conditions the flow is steady and laminar. Three single-roll-type structures and an unstable toroidal roll have been observed inside the cavity with nearly adiabatic sidewalls. The sequence from the conductive state consists of a toroidal roll that evolves to a diagonally oriented single roll with increasing Rayleigh number. This diagonal roll, which is stabilized by the effect of the small but finite conductivity of the walls, shifts its axis of rotation towards to two opposite walls, and back to the diagonal orientation to allow for the increase in circulation that occurs as the Rayleigh number is further increased. Conduction at the sidewalls modifies this sequence in the sense that the two initial single rolls finally evolve into a four-roll structure. Once formed, this four-roll structure remains stable when decreasing the Rayleigh number until the initial single diagonally oriented roll is again recovered. The topology and the velocity fields of all structures, characterized with visualization and particle image velocimetry, respectively, are in good agreement with numerical results reported previously for the cavity with adiabatic walls, as well as with the numerical predictions obtained in the present study for perfectly conducting lateral walls. Received: 10 August 1998/Accepted: 1 August 2000     

Streamwise evolution of an inclined cylinder wake

The streamwise evolution of an inclined circular cylinder wake was investigated by measuring all three velocity and vorticity components using an eight-hotwire vorticity probe in a wind tunnel at a Reynolds number Re_d of 7,200 based on free stream velocity (U _∞) and cylinder diameter (d). The measurements were conducted at four different inclination angles (α), namely 0°, 15°, 30°, and 45° and at three downstream locations, i.e., x/d = 10, 20, and 40 from the cylinder. At x/d = 10, the effects of α on the three coherent vorticity components are negligibly small for α ≤ 15°. When α increases further to 45°, the maximum of coherent spanwise vorticity reduces by about 50%, while that of the streamwise vorticity increases by about 70%. Similar results are found at x/d = 20, indicating the impaired spanwise vortices and the enhancement of the three-dimensionality of the wake with increasing α. The streamwise decay rate of the coherent spanwise vorticity is smaller for a larger α. This is because the streamwise spacing between the spanwise vortices is bigger for a larger α, resulting in a weak interaction between the vortices and hence slower decaying rate in the streamwise direction. For all tested α, the coherent contribution to [`(v<sup>2</sup>)] \overline{{v^{2}}} is remarkable at x/d = 10 and 20 and significantly larger than that to [`(u²)] \overline{{u^{2}}} and [`(w<sup>2</sup>)]. \overline{{w^{2}}}. This contribution to all three Reynolds normal stresses becomes negligibly small at x/d = 40. The coherent contribution to [`(u²)] \overline{{u^{2}}} and [`(v²)] \overline{{v^{2}}} decays slower as moving downstream for a larger α, consistent with the slow decay of the coherent spanwise vorticity for a larger α.     

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.     

Effect of the use of the balance and gradient methods as a result of experimental investigations of natural convection action with regard to the conception and construction of measuring apparatus

Measurement apparatus designed and constructed according to conceptions of the authors, enabled a more precise calculation of the heat transfer coefficient with the balance and gradient methods. Construction and use of the apparatus and devices are described below, results of experimental investigations for horizontal and vertical, isothermal, flat plates obtained independently with the balance and gradient methods, are also presented. The following equations were found:Nu=0.612 · (Ra)^1/4 10⁴ ≦Ra ≦ 10⁸ for vertical platesNu=0.766 · (Ra)^1/5 10⁴ ≦Ra ≦ 10⁷ Nu=0.173 · (Ra)^1/3 10⁵ ≦Ra≦ 10⁸ for horizontal plates. On the basis of the results obtained from both these methods, differences of natural convection acting from vertical and horizontal plates are discussed. The usefulness of the balance and gradient methods have been considered for qualitative and quantitative investigations of heat transfer by natural convection.     

Prediction of natural convection flow using network model and numerical simulations inside enclosure with distributed solid blocks

Steady state natural convection of a fluid with Pr ≈ 1 within a square enclosure containing uniformly distributed, conducting square solid blocks is investigated. The side walls are subjected to differential heating, while the top and bottom ones are kept adiabatic. The natural convection flow is predicted employing the nondimensional volumetric flow rate (Q_max^* Q_{\max }^{*} ) by using a network model and also using numerical simulations. For identical solid and fluid thermal conductivities (i.e. k _s  = k _f ), a parametric study of the effect of number of blocks (N ²), gap size (δ) and enclosure Rayleigh number (Ra) on Q_max^* Q_{\max }^{*} is performed using the two approaches. Network model predictions are observed to agree well with that from the simulations until Raδ³ ~ 12. Considering the enclosure with blocks as a porous medium, for a fixed enclosure Ra number, increasing the number of blocks for a fixed volumetric porosity leads to a decrease in enclosure permeability, which in turn reduces the flow rate. When the number of blocks is fixed, and for a given Ra number, the flow rate increases as the porosity increases by widening the gap between the blocks.     

Optimal design of angled rib turbulators in a cooling channel

In the present study, an optimal design for enhancement of heat transfer and thermal performance for a stationary channel with angled rib turbulators was investigated to find the most suitable rib geometry. Among various design parameters, two design variables, rib angle of attack (α) and pitch-to-rib height (p/e), were chosen. The ranges of two design variables were set as 30° ≤ α ≤ 80° and 3.0 ≤ p/e ≤ 15.0. Approximations for design of the best rib turbulators were obtained using the advanced response surface method with functional variables. The second-order response surfaces (or correlations) within the ranges of two design variables were completed by this method. As for the optimized results, maximum averaged heat transfer value was obtained at α = 53.31° and p/e = 6.50, while the highest thermal performance value was presented at α = 54.67° and p/e = 6.80.