Mixed convection flow and heat transfer across a square cylinder under the influence of aiding buoyancy at low Reynolds numbers

In this paper, mixed convection flow and heat transfer around a long cylinder of square cross-section under the influence of aiding buoyancy are investigated in the vertical unconfined configuration (Reynolds number, Re = 1–40 and Richardson number, Ri = 0–1). The semi-explicit finite volume method implemented on the collocated grid arrangement is used to solve the governing equations along with the appropriate boundary conditions. The onset of flow separation occurs between Re = 1–2, between Re = 2–3 and between Re = 3–4 for Ri = 0, 0.5 and 1, respectively. The flow is found to be steady for the range of conditions studied here. The friction, pressure and total drag coefficients are found to increase with Richardson number, i.e., as the influence of aiding buoyancy increases drag coefficients increase at the constant value of the Reynolds number. The temperature field around the obstacle is presented by isotherm contours at the Prandtl number of 0.7 (air). The local and average Nusselt numbers are calculated to give a detailed study of heat transfer over each surface of the square cylinder and an overall heat transfer rate and it is found that heat transfer increases with increase in Reynolds number and/or Richardson number. The simple expressions for the wake length and average cylinder Nusselt number are obtained for the range of conditions covered in this work.     

Numerical prediction of flow and heat transfer of power-law fluids in a plane channel with a built-in heated square cylinder

Structure of unsteady laminar flow and heat transfer of power-law fluids in two-dimensional horizontal plane channel with a built-in heated square cylinder is studied numerically. The governing equations are solved using a control volume finite element method (CVFEM) adapted to the staggered grid. Computations are performed over a range of Reynolds and Richardson numbers from Re = 20 to 200 and from Ri = 0 to 8, respectively at fixed Prandtl number Pr = 50 and blockage ratio value β′ = 1/8. Three different values of the power-law index (n = 0.5, 1 and 1.4) are considered in this study to show its effect on the value of the critical Reynolds number defining the transition between two different flow regimes (symmetrical and periodic flows), the variations of Strouhal number, drag and lift coefficients and the heat transfer from the square cylinder as function of Reynolds number. Heat transfer correlations are obtained through forced convection. A discussion about the buoyancy effect on the flow pattern and the heat transfer for different power-law index is also presented.     

Numerical study of swirling flows in a cylindrical container with co-/counter-rotating end disks under stable temperature difference

A numerical study was conducted to investigate swirling flows of a Boussinesq fluid confined in a cylindrical container with co-/counter-rotating end disks. A vertically stable temperature gradient is imposed, with the stationary sidewall assumed as adiabatic. Flows are studied for a range of parameters: the Reynolds number, Re, 100 ⩽ Re ⩽ 2000; the Richardson number, Ri, 0 ⩽ Ri ⩽ 1.0; and the Prandtl number, Pr, Pr = 1.0. The ratio of the angular velocity of the top disk to the bottom disk, s, −1.0 ⩽ s ⩽ 1.0. The cylinder aspect ratio: h = 2.0. For the case of negligibly small temperature difference (Ri ∼ 0) and high Re, interior fluid rotates with an intermediate angular velocity of both end disks when they are co-rotating (s > 0). When end disks are counter-rotating (s < 0), shearing flow with meridional recirculation is created. For the case of large temperature difference (Pr · Ri ∼ O(1)), the Ekman suction is suppressed and the sidewall boundary layer disappears at mid-height of the cylinder. For all the values of s considered in the present study, the bulk of the fluid is brought close to rest with the fluid in the vicinity of both end disks rotating in each direction. The secondary flow in the meridional section of the cylindrical container exhibits various types of vortices as the governing parameters are varied. These flow patterns are presented in the form of diagrams on the (s, Re) plane and (s, Ri) plane. The average Nusselt number is computed and presented as functions of Ri, Re and s.     

Numerical simulation of double diffusive mixed convection in an open enclosure with different cylinder locations

This article reports numerical simulation of the double diffusive mixed convection around a cylinder in an open enclosure with an inlet and exit ports. The temperature and mass concentration of the cylinder are higher than those of the inlet flow and the cylinder can be at three different locations (lower, middle and upper) in the enclosure. The inlet flow with low temperature and mass concentration is located at the lower-left wall of the enclosure and the exit is at the upper-right wall. Other walls are assumed to be adiabatic. Effects of Lewis number Le, buoyancy ratio Br, and cylinder locations on the double diffusive mixed convection are investigated at Richardson number Ri = 1.0 and 0.01 while Prandtl number Pr is kept at 0.7. Streamlines, isotherms, isoconcentrations, and the average and local Sherwood number at different parameters are reported to characterize the double diffusive mixed convection phenomena in the open enclosure.     

Assessment of dissipative particle dynamics to simulate combined convection heat transfer: Effect of compressibility

The current study explored the capability of a discrete particle method known as dissipative particle dynamics with energy conservation (eDPD) to simulate combined convection heat transfer in a vertical lid driven cavity. The study investigated two cases of aiding and opposing buoyancy mechanisms in the lid driven cavity. The eDPD results were compared against the finite volume solutions for the range of Richardson number, 10^− 2 ≤ Ri ≤ 10². The method showed good comparison for the range of Richardson number 10^− 2 ≤ Ri ≤ 10¹. However, the eDPD method showed deviation from the FV solutions for a high value of Richardson number, Ri = 10², and this deviation is attributed to the compressibility of eDPD system experienced at such high value of Richardson number. Parametric study on the influence of the Richardson number (Ri) on the eDPD compressibility was conducted and presented via temperature isotherms, streamlines, velocity contours, velocity vectors, temperature and velocity profiles.     

Unsteady wake dynamics and heat transfer in forced and mixed convection past a circular cylinder in cross flow for high Prandtl numbers

This paper investigates the combined effect of Prandtl number and Richardson number on the wake dynamics and heat transfer past a circular cylinder in crossflow using a SUPG based finite element method. The computations are carried out for 80 < Re < 180, 0.7 < Pr < 100 and $0?\mathrm{Ri}?2$. The results have been presented for both forced and mixed convection flows. In the case of forced convection, crowding of temperature contours with reduced spatial spread is observed for increasing Prandtl numbers. The local and average Nusselt numbers are found to increase with increasing Reynolds number and Prandtl number. The average Nusselt number and Colburn factor are found to vary as Re^0.548 Pr ^0.373 and Re^?0.452, respectively. The extrapolated results of the average Nusselt number for low and high Reynolds numbers are found to match quite well with the available results in literature. Effect of Prandtl number shows various interesting phenomena for the mixed convective flows. Increasing the Prandtl numbers resulted in decreasing deflection and strength in the wake structures. The effect of baroclinic vorticity production during vortex shedding has been demonstrated at the vicinity of the cylinder and near wake. The Strouhal number is found to decrease with increasing Prandtl number, in the case of buoyancy induced flow. The effect of increasing Prandtl number is manifested as the stabilizing effect in the flow. This is, perhaps, the first time that such behavior for the Prandtl number is being reported. Additionally it is observed that the average Nusselt number decreases with increasing Richardson number.     

A numerical study of fluid flow passes two heated/cooled square cylinders in a tandem arrangement via lattice Boltzmann method

In this paper, the fluid flow pass two heated/cooled square cylinders in a tandem arrangement is simulated via the Multiple-Relaxation-Time lattice Boltzmann method. The distance between the upstream and downstream cylinder varies from the rear of the upstream one to 5 times of the cylinder width. The numerical experiments are done under different Richardson numbers (Ri, represents the effect of the buoyancy force) for two typical Re = 100, 60. The buoyancy effect on the flow and heat transfer around the two cylinders is mainly investigated. As is shown, if the force is in the same direction of incoming flow, the vortex street is always suppressed and no critical spacing seems to exist. However, if the force is in the opposite direction of the incoming flow, the vortex street can always be generated and the critical spacing always seems to exist. Correspondingly, the heat transfer around the cylinders measured by the Nusselt number on the surfaces of the two cylinders also shows different characteristics for various Ri s.     

Mixed convective transport in a lid-driven cavity containing a nanofluid and a rotating circular cylinder at the center

The mixed convective transport of Cu-H₂O nanofluid in a differentially heated and lid-driven square enclosure in the presence of a rotating circular cylinder is investigated numerically. The top wall of the enclosure is sliding from left to right at a uniform speed while all other walls are stationary. A thermally insulated circular cylinder is placed centrally within the enclosure. The cylinder can rotate about its centroidal axis. The top and bottom walls are kept isothermal at different temperatures while the side walls are assumed adiabatic. Simulations are performed for, Richardson number 1 ≤ Ri ≤ 10, dimensionless rotational speed 0 ≤ Ω ≤ 5 and nanoparticle concentration 0 ≤ ϕ ≤ 0.20 keeping the Grashof number fixed as Gr = 10⁴. The flow and thermal fields are analyzed through streamline and isotherm plots for various Ω and Ri. Furthermore, the drag coefficient of the moving lid and Nusselt number of the hot wall are also computed to understand the effects of Ω and Ri on them. It is observed that the heat transfer greatly depends on the rotational speed of the cylinder, mixed convective strength and the nanoparticle concentration.     

Momentum and heat transfer characteristics of a semi-circular cylinder immersed in power-law fluids in the steady flow regime

The continuity, momentum and energy equations describing the flow and heat transfer of power-law fluids over a semi-circular cylinder have been solved numerically in the two-dimensional steady flow regime. The influence of the Reynolds number (Re), Prandtl number (Pr) and power-law index (n) on the local and global flow and heat characteristics have been studied over wide ranges of conditions as follows: 0.01 ? Re ? 30, 1 ? Pr ? 100 and 0.2 ? n ? 1.8. The variation of drag coefficient and Nusselt number with the Reynolds number, Prandtl number and power-law index is shown over the aforementioned ranges of conditions. In addition, streamline and isotherm profiles along with the recirculation length and distribution of pressure coefficient and Nusselt number over the surface of the semi-circular cylinder are also presented to gain further insights into the nature of the underlying kinematics. The wake size (recirculation length) shows almost linear dependence on the Reynolds number (Re ? 1) for all values of power-law index studied herein. The drag values show the classical inverse variation with the Reynolds number, especially for shear-thinning fluids at low Reynolds numbers. The point of maximum pressure coefficient is found slightly displaced from the front stagnation point for highly shear-thinning fluids, whereas for shear-thickening and Newtonian fluids, it coincides with the front stagnation point. For fixed values of the Prandtl number and Reynolds number, the rate of heat transfer decreases with the gradual increase in power-law index; this effect is particularly striking at high Prandtl numbers due to the thinning of the thermal boundary layer. Conversely, as expected, shear-thinning behavior facilitates heat transfer and shear-thickening impedes it. The effect of power-law index on both momentum and heat-transfer characteristics is seen to be appreciable at low Reynolds numbers and it gradually diminishes with the increasing Reynolds number.     

Laminar mixed convection in a channel with a built-in semi-circular cylinder under the effect of cross-buoyancy

Effects of cross-buoyancy mixed convection on flow and heat transfer characteristics of a long semi-circular cylinder (long in neutral direction) in a confined channel have been investigated in the laminar regime. The numerical results have been presented and discussed for the range of conditions as Reynolds number (Re) = 1–40, Richardson number (Ri) = 0–4, Prandtl number (Pr) = 0.71–50 and blockage ratio (β) = 16.67%–50%. The drag coefficient increases with increasing Richardson number and/or blockage ratio. The average Nusselt number is showing a maximum relative enhancement of approximately 45% for Ri = 4 with respect to corresponding forced convection value (Ri = 0). The average Nusselt number increases with increase in Prandtl number and shows a maximum relative enhancement of approximately 1136% for Pr = 50 with respect to corresponding value at Pr = 0.71. On the other hand, the maximum relative variation of the total drag coefficient is found to be approximately 55% for Ri = 4 with respect to corresponding value at Ri = 0. Finally, the simple heat transfer correlation is obtained for the proceeding range of control parameters.     

Flow over and forced convection heat transfer in Newtonian fluids from a semi-circular cylinder

Momentum and heat transfer characteristics of a semi-circular cylinder immersed in unconfined flowing Newtonian fluids have been investigated numerically. The governing equations, namely, continuity, Navier–Stokes and energy, have been solved in the steady flow regime over wide ranges of the Reynolds number (0.01 ? Re ? 39.5) and Prandtl number (Pr ? 100). Prior to the investigation of drag and heat transfer phenomena, the critical values of the Reynolds number for wake formation (0.55 < Re^c < 0.6) and for the onset of vortex shedding (39.5 < Re_c < 40) have been identified. The corresponding values of the lift coefficient, drag coefficient, and Strouhal number are also presented. After establishing the limit of the steady flow regime, the influence of the Reynolds number (0.01 ? Re ? 39.5) and Prandtl number (Pr = 0.72, 1, 10, 50 and 100) on the global flow and heat transfer characteristics have been elucidated. Detailed kinematics of the flow is investigated in terms of the streamline and vorticity profiles and the variation of pressure coefficient in the vicinity of the cylinder. The functional dependence of the individual and total drag coefficients on the Reynolds number is explored. The Nusselt number shows an additional dependence on the Prandtl number. In addition, the isotherm profiles, local Nusselt number (Nu_L) and average Nusselt number (Nu) are also presented to analyze the heat transfer characteristic of a semi-circular cylinder in Newtonian media.     

Characteristics for flow and heat transfer around a circular cylinder near a moving wall in wide range of low Reynolds number

The present study numerically investigates two-dimensional laminar fluid flow and heat transfer past a circular cylinder near a moving wall. Numerical simulations to calculate the fluid flow and heat transfer past a circular cylinder are performed for different Reynolds numbers varying in the range of 60–200 and a fixed Prandtl numbers of 0.7 (air) in the range of 0.1 ? G/D ? 4, where G/D is the ratio of the gap between the cylinder and a moving wall, G and the cylinder diameter, D. The flow and thermal fields become the steady state below the critical gap ratios of 0.8, 0.4 and 0.2 for the Reynolds numbers of 60, 80 and 100, respectively. As the gap ratio decreases, the magnitude of lift coefficient for all Reynolds numbers increased significantly with diminishing G/D due to the ground effect. The cases of Reynolds numbers of 60, 80 and 100 revealed the sharp slope of drag coefficient in the range of the gap ratio where the flow transfers from the unsteady state to the steady state. As the Reynolds number decreases, the variation of Nusselt is much significant and increases considerably with decreasing G/D.     

Two-dimensional unsteady forced convection heat transfer in power-law fluids from a cylinder

Forced convection heat transfer characteristics of a cylinder (maintained at a constant temperature) immersed in a streaming power-law fluids have been studied numerically in the two-dimensional (2-D), unsteady flow regime. The governing equations, namely, continuity, momentum and thermal energy, have been solved using a finite volume method based solver (FLUENT 6.3) over wide ranges of conditions (power law index, 0.4 ? n ? 1.8; Reynolds number, 40 ? Re ? 140; Prandtl number, 1 ? Pr ? 100). In particular, extensive numerical results elucidating the influence of Reynolds number, Prandtl number and power-law index on the isotherm patterns, local and average Nusselt numbers and their evolution with time are discussed in detail. Over the ranges of conditions considered herein, the nature of flow is fully periodic in time. The heat transfer characteristics are seen to be influenced in an intricate manner by the value of the Reynolds number (Re), Prandtl number (Pr) and the power-law index (n). Depending upon the value of the power-law index (n), though the flow transits from being steady to unsteady somewhere in the range ～33 < Re < 50, the fully periodic behavior is seen only beyond the critical value of the Reynolds number (Re). As expected, the average Nusselt number increases with an increase in the values of Reynolds and/or Prandtl numbers, irrespective of the value of the flow behavior index. A strong influence of the power-law index on both local and time-averaged Nusselt numbers was observed. Broadly, all else being equal, shear-thinning behavior (n < 1) promotes heat transfer whereas shear-thickening behavior (n > 1) impedes it. Furthermore, this effect is much more pronounced in shear-thinning fluids than that in shear-thickening fluids.     

Steady forced convection heat transfer from a heated circular cylinder to power-law fluids

Forced convection heat transfer to incompressible power-law fluids from a heated circular cylinder in the steady cross-flow regime has been investigated numerically by solving the momentum and thermal energy equations using a finite volume method and the QUICK scheme on a non-uniform Cartesian grid. The dependence of the average Nusselt number on the Reynolds number (5 ⩽ Re ⩽ 40), power-law index (0.6 ⩽ n ⩽ 2) and Prandtl number (1 ⩽ Pr ⩽ 1000) has been studied in detail. The numerical results are used to develop simple correlations as functions of the pertinent dimensionless variables. In addition to the average Nusselt number, the effects of Re, Pr and n on the local Nusselt number distribution have also been studied to provide further physical insights. The role of the two types of thermal boundary conditions, namely, constant temperature and uniform heat flux on the surface of the cylinder has also been presented.     

Combined convection heat transfer for thermally developing aiding flow in an inclined circular cylinder with constant heat flux

Experiments have been conducted to study the local and average heat transfer by mixed convection for hydrodynamically fully developed, thermally developing and thermally fully developed laminar air flow in an inclined circular cylinder. The experimental setup consists of aluminum cylinder as test section with 30 mm inside diameter and 900 mm heated length (L/D = 30), is subjected to a constant wall heat flux boundary condition. The investigation covers Reynolds number range from 400 to 1600, heat flux is varied from 70 W/m² to 400 W/m² and cylinder angles of inclination including 30°, 45° and 60°. The hydrodynamically fully developed condition has been achieved by using aluminum entrance section pipes (calming sections) having the same inside diameter as test section pipe but with variable lengths. The entrance sections included two long calming sections, one with length of 180 cm (L/D = 60), another one with length of 240 cm (L/D = 80) and two short calming sections with lengths of 60 cm (L/D = 20), 120 cm (L/D = 40). The results present the surface temperature distribution along the cylinder length, the local and average Nusselt number distribution with the dimensionless axial distance Z⁺. For all entrance sections, the results showed an increase in the Nusselt number values as the heat flux increases and as the angle of cylinder inclination moves from θ = 60° inclined cylinder to θ = 0° horizontal cylinder. The mixed convection regime has been bounded by the convenient selection of Re number range and the heat flux range, so that the obtained Richardson numbers (Ri) is varied approximately from 0.13 to 7.125. The average Nusselt numbers have been correlated with the (Rayleigh numbers/Reynolds numbers) in empirical correlations.     

Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids

The behaviour of nanofluids is investigated numerically inside a two-sided lid-driven differentially heated square cavity to gain insight into convective recirculation and flow processes induced by a nanofluid. A model is developed to analyze the behaviour of nanofluids taking into account the solid volume fraction χ. The transport equations are solved numerically with finite volume approach using SIMPLE algorithm. Comparisons with previously published work on the basis of special cases are performed and found to be in excellent agreement. The left and the right moving walls are maintained at different constant temperatures while the upper and the bottom walls are thermally insulated. Three case were considered depending on the direction of the moving walls. Governing parameters were 0.01 < Ri < 100 but due to space constraints only the results for 0.1 < Ri < 10 are presented. It is found that both the Richardson number and the direction of the moving walls affect the fluid flow and heat transfer in the cavity. Copper–Water nanofluid is used with Pr = 6.2 and solid volume fraction χ is varied as 0.0%, 8%, 16% and 20%. Detailed results are presented for flow pattern and heat transfer curves.     

Experimental investigation of mixed convection heat transfer for thermally developing flow in a horizontal circular cylinder

An experimental investigation is presented on mixed (free and forced) convection to study the local and average heat transfer for hydrodynamically fully developed, thermally developing and thermally fully developed laminar air flow in a horizontal circular cylinder. The experimental setup consists of aluminum cylinder as test section with 30 mm inside diameter and 900 mm heated length (L/D = 30), is subjected to a constant wall heat flux boundary condition. The investigation covers Reynolds number range from 400 to 1600, the heat flux varied from 60 W/m² to 400 W/m² and with cylinder inclination angle of θ = 0° (horizontal). The hydrodynamically fully developed condition is achieved by using an aluminum entrance section pipes (calming sections) having the same inside diameter as test section pipe but with variable lengths. The entrance sections included two long calming sections, one with length of 180 cm (L/D = 60), another one with length of 240 cm (L/D = 80) and two short calming sections with lengths 60 cm (L/D = 20), 120 cm (L/D = 40). The surface temperature variation along the cylinder surface, the local and average Nusselt number variation with the dimensionless axial distance Z⁺ were presented. For all entrance sections, it was found an increase in the Nusselt number values as the heat flux increases. It was concluded that the free convection effects tended to decrease the heat transfer results at low Re while to increase the heat transfer results for high Re. The combined convection regime could be bounded by a suitable selection of Re number ranges and the heat flux ranges. The obtained Richardson numbers (Ri) range varied approximately from 0.13 to 7.125. The average Nusselt numbers were correlated with the (Rayleigh numbers/Reynolds numbers). The proposed correlation has been compared with available literature and showed satisfactory agreement.     

Heat transfer in a small gap between co-axial rotating cylinders

This paper investigates the local heat transfer of a co-axial rotating cylinder. In the inner flow field of the rotating cylinder, the dimensionless parameters include the rotational Reynolds number (Re_Ω) and buoyancy parameter (Gr). The test rig is designed to make the rotating in the inner cylinder and stationary in the outer cylinder. The local temperature distributions of the inner and outer cylinder on axial direction were measured. Under the experimental condition, whereas the ranges of the rotational Reynolds number are 2400 ≤ Re_Ω ≤ 45,000. Experimental results reveal that the rotational Reynolds number's increase is with the heat transfer coefficient distributions increase types. Finally, the local heat transfer rate on the wall are correlated and compared with that in the existing literature.     

Effects of Reynolds and Prandtl numbers on heat transfer from a square cylinder in the unsteady flow regime

The effects of the Reynolds and Prandtl numbers on the rate of heat transfer from a square cylinder are investigated numerically in the unsteady two-dimensional periodic flow regime, for the range of conditions 60 ? Re ? 160 and 0.7 ? Pr ? 50 (the maximum value of Peclet number being 4000). A semi-explicit finite volume method has been used on a non-uniform collocated grid arrangement to solve the governing equations. Using the present numerical results, simple heat transfer correlations are obtained for the constant temperature and constant heat flux conditions on the solid square cylinder. In addition, the variation of the time averaged local Nusselt number on the each face of the obstacle and representative isotherm plots are presented to elucidate the role of Prandtl number on heat transfer in the unsteady flow regime.     

Forced convection heat transfer from an elliptical cylinder to power-law fluids

Forced convection heat transfer to incompressible power-law fluids from a heated elliptical cylinder in the steady, laminar cross-flow regime has been studied numerically. In particular, the effects of the power-law index (0.2 ? n ? 1.8), Reynolds number (0.01 ? Re ? 40), Prandtl number (1 ? Pr ? 100) and the aspect ratio of the elliptic cylinder (0.2 ? E ? 5) on the average Nusselt number (Nu) have been studied. The average Nusselt number for an elliptic cylinder shows a dependence on the Reynolds and Prandtl numbers and power-law index, which is qualitatively similar to that for a circular cylinder. Thus, heat transfer is facilitated by the shear-thinning tendency of the fluid, while it is generally impeded in shear-thickening fluids. The average Nusselt number values have also been interpreted in terms of the usual Colburn heat transfer factor (j). The functional dependence of the average Nusselt number on the dimensionless parameters (Re, n, Pr, E) has been presented by empirically fitting the numerical results for their easy use in process design calculations.