Vortex formation and heat transfer in turbulent flow past a transverse cavity with inclined frontal and rear walls

The process of vortex formation, distributions of pressure coefficients, and convective heat transfer in a turbulent flow past a cavity with a low aspect ratio and inclined frontal and rear walls were experimentally studied. The angle of wall inclination φ was varied in the interval from 30° to 90°. Visualization techniques were applied to trace the evolution of the flow with the angle φ as the transverse cavity became more open. Pressure fields in the longitudinal and transverse sections on the bottom wall of the cavity, and on its frontal and rear walls, were measured. The measured distributions of temperature in the longitudinal and transverse sections on the three heated walls, and the obtained thermographic fields over the whole heated surface, were used to calculate local and average heat-transfer coefficients. It is found that in the interval of wall inclination angles φ = 60–70° the flow in the cavity becomes unstable, with the primary vortex changing its structure from single-cellular to double-cellular. As a result, the distributions of static pressure and surface temperature across and along the cavity suffer dramatic changes. At smallest angles φ the flow re-attachment point gets displaced into the cavity to cause an abrupt growth of pressure and heat-transfer coefficients on the rear wall, which leads to a slight increase of the surface-mean pressure and heat transfer inside the cavity. At the angle of instability, φ = 60°, the local heat-transfer coefficient decreases markedly over the cavity span from the end faces of the cavity toward its center, and a most pronounced intensification of heat transfer is observed.     

The effect of the top and bottom wall temperatures on the laminar natural convection in an air-filled square cavity

The effect of the top and bottom wall temperatures on the natural convection heat transfer characteristics in an air-filled square cavity driven by a difference in the vertical wall temperatures was investigated by measuring the temperature distributions along the heated vertical wall and visualizing the flow patterns in the cavity. The experiments were performed at a horizontal Grashof number of 1.9 × 10⁸. Increasing the top wall temperature resulted in a separated flow region on the top wall, which caused a secondary flow between the separated flow and the boundary layer on the heated vertical wall. This secondary flow had a significant effect on the heat transfer in this region. Changes in the top and bottom wall temperatures changed the temperature gradient and the average temperature of the air outside the thermal boundary layers in the cavity. The local heat transfer along much of the heated vertical wall could be correlated by Nu = C · Ra^0.32, but the constant C increased when the average of the top and bottom wall temperatures increased.     

Experimental investigation of opposing mixed convection hear transfer in the vertical flat channel in a laminar–turbulent transition region

In this paper we present the results on experimental investigation of the local opposing mixed convection heat transfer in the vertical flat channel with symmetrical heating in a laminar–turbulent transition region. The experiments were performed in airflow (p = 0.1–1.0 MPa) in the range of Re from 1.5 × 10³ to 6.6 × 10⁴ and Gr_q up to 1 × 10¹¹ at the limiting condition q_w1 = q_w2 = const. The analysis of the results revealed significant increase in the heat transfer with increasing of air pressure (Gr number). Also sharp increase in heat transfer was noticed in the region with vortex flow in comparison with the turbulent flow region.     

Analysis of combined thermal and magnetic convection ferrofluid flow in a cavity

A numerical analysis of the magnetic gradient and thermal buoyancy induced cavity ferrofluid flow is conducted by a semi-implicit finite element method. The physical model for a square cavity containing two different temperature side walls and a magnet near bottom wall is described by mass, momentum and energy equations. Conditions for the fixed Prandtl number, Rayleigh number and different ferro-hydrodynamic interaction parameter are studied for 5 × 10⁸ ≤ β ≤ 1.6 × 10¹⁰. Results show the flow strength increases with the strengthening magnetic field. However, the side-wall heat transfer rate presents a decrease trend with the increase in magnetic field strength, for the magnet located near the bottom central area evokes the circulation to move toward the central portion. In summary, a proper choice of magnet strength and location can adjust the flow field and local heat transfer rate to fit the practical application.     

Experimental study of thermal behaviors in a rectangular channel with baffle of pores

This experimental study attempts to explore the local heat transfer in rectangular channel with baffles, and analyzes the experimental results of baffles with different heights and pores in the event of five Reynolds numbers and three heating quantities. Apart from increasing the perturbation of flow field, the channel's flow field with baffles, which is similar to a backward-facing step flow field, is very helpful to heat transfer. To obtain an optimized baffle and increase the perturbation of flow field, this experiment employed baffles with five heights (H = 10–50 mm) and different numbers of pores (N = 1–3), as well as heat flux: Q = 40–100 l/min, Reynolds number: 702–1752, and heating quantity: q_in = 90–750 W/m². In addition to measurement of overall temperature distribution, emphasis is also placed on analysis of local heat transfer coefficient. Furthermore, heat transfer distribution of channel can be applied to explain how the baffles of pores have an influence upon backward-facing step flow field, shear layer, recirculation region, reattachment region and redeveloped boundary layer. Finally, some empirical formulas derived form experimental results may provide a reference for future design.     

Heatlines based natural convection analysis in tilted isosceles triangular enclosures with linearly heated inclined walls: effect of various orientations

Natural convection in isosceles triangular enclosures with various configurations (case 1 — inverted, case 2 — straight and case 3 — tilted) is studied via heatline analysis for linear heating of inclined walls. Detailed analysis and comparison for various base angles (φ = 45°, 60°) of triangular enclosures have been carried out for a range of fluids (Pr = 0.015 − 1000) within Ra = 10³ − 10⁵ using Galerkin finite element method. The heat flow distributions indicate conduction dominant heat transfer at low Ra (Ra = 10³) for case 1 and case 2 whereas in case 3, convective heat flow is observed due to high buoyancy force. As Ra increases, enhanced thermal mixing is observed at the core of the cavity. Wall to wall heat transfer occurs at walls AB and AC due to linear heating boundary condition in all the cases. Although the distributions of fluid flow and heat flow are qualitatively similar for φ = 45° and 60°, the intensity of fluid flow and heat flow decreases as φ increases. Strength of fluid flow and heat flow circulation cells is found to be higher in case 3 for identical parameters. Results show that upper side wall (AC) for case 3 exhibits higher heat transfer rates whereas heat transfer rates for walls AB and AC are the same for case 1 and case 2. Also Nu_AB is higher for case 2 followed by case 1 and case 3 at the middle portion of wall AB. Thus to achieve high heat transfer from fluid to wall at the central region, case 2 and case 3 configurations may be recommended at high Ra (Ra = 10⁵) and Pr, irrespective of φ.     

Finite element based heatline approach to study mixed convection in a porous square cavity with various wall thermal boundary conditions

A penalty finite element method based simulation is performed to analyze the influence of various walls thermal boundary conditions on mixed convection lid driven flows in a square cavity filled with porous medium. The relevant parameters in the present study are Darcy number (Da = 10^?5 ? 10^?3), Grashof number (Gr = 10³ ? 10⁵), Prandtl number (Pr = 0.7–7.2), and Reynolds number (Re = 1–10²). Heatline approach of visualizing heat flow is implemented to gain a complete understanding of complex heat flow patterns. Patterns of heatlines and streamlines are qualitatively similar near the core for convection dominant flow for Da = 10^?3. Symmetric distribution in heatlines, similar to streamlines is observed irrespective of Da at higher Gr in natural convection dominant regime corresponding to smaller values of Re. A single circulation cell in heatlines, similar to streamlines is observed at Da = 10^?3 for forced convection dominance and heatlines are found to emanate from a large portion on the bottom wall illustrating enhanced heat flow for Re = 100. Multiple circulation cells in heatlines are observed at higher Da and Gr for Pr = 0.7 and 7.2. The heat transfer rates along the walls are illustrated by the local Nusselt number distribution based on gradients of heatfunctions. Wavy distribution in heat transfer rates is observed with Da ? 10^?4 for non-uniformly heated walls primarily in natural convection dominant regime. In general, exponential variation of average Nusselt numbers with Grashof number is found except the cases where the side walls are linearly heated. Overall, heatlines are found to be a powerful tool to analyze heat transport within the cavity and also a suitable guideline on explaining the Nusselt number variations.     

Heat flow visualization for natural convection in rhombic enclosures due to isothermal and non-isothermal heating at the bottom wall

Analysis has been carried out for the energy distribution and thermal mixing in steady laminar natural convective flow through the rhombic enclosures with various inclination angles, φ for various industrial applications. Simulations are carried out for various regimes of Prandtl (Pr) and Rayleigh (Ra) numbers. Dimensionless streamfunctions and heatfunctions are used to visualize the flow and energy distribution, respectively. Multiple flow circulations are observed at Pr = 0.015 and 0.7 for all φs at Ra = 10⁵. On the other hand, two asymmetric flow circulation cells are found to occupy the entire cavity for φ = 75° at higher Pr (Pr = 7.2 and 1000) and Ra (Ra = 10⁵). Heatlines are found to be parallel circular arcs connecting the cold and hot walls for the conduction dominant heat transfer at Ra = 10³. The enhanced convective heat transfer is explained with dense heatlines and convective loop of heatlines at Ra = 10⁵. Heatlines clearly demonstrate that the left wall receives heat from the bottom wall as heatlines directly connect both the walls whereas the convective heat circulation cells play lead role to distribute the heat along the right wall, especially for smaller φs. On the other hand, the heat flow is evenly distributed to both side walls at higher φs via convection as well as direct conductive transport. Significant convective heat transfer from the bottom hot wall to the left cold wall occurs for φ = 30° cavity whereas the heat transfer to the right cold wall is maximum for φ = 75° irrespective of Pr. Average Nusselt number studies also show that φ = 30° cavity gives maximum heat transfer rate from the bottom to left wall irrespective of Pr in isothermal heating case. On the other hand, enhanced thermal mixing occurs at φ = 75° for both isothermal and non-isothermal heating strategies except at Pr = 0.015 in isothermal heating case.     

Analysis of Bejan’s heatlines on visualization of heat flow and thermal mixing in tilted square cavities

This article analyzes the detailed heat transfer phenomena during natural convection within tilted square cavities with isothermally cooled walls (BC and DA) and hot wall AB is parallel to the insulated wall CD. A penalty finite element analysis with bi-quadratic elements has been used to investigate the results in terms of streamlines, isotherms and heatlines. The present numerical procedure is performed over a wide range of parameters (10³ ? Ra ? 10⁵,0.015 ? Pr ? 1000,0° ? φ ? 90°). Secondary circulations cells are observed near corner regions of cavity for all φ’s at Pr = 0.015 with Ra = 10⁵. Two asymmetric flow circulation cells are found to occupy the entire cavity for φ = 15° at Pr = 0.7 and Pr = 1000 with Ra = 10⁵. Heatlines indicate that the cavity with inclination angle φ = 15° corresponds to large convective heat transfer from the wall AB to wall DA whereas the heat transfer to wall BC is maximum for φ = 75°. Heat transfer rates along the walls are obtained in terms of local and average Nusselt numbers and they are explained based on gradients of heatfunctions. Average Nusselt number distributions show that heat transfer rate along wall DA is larger for lower inclination angle (φ = 15°) whereas maximum heat transfer rate along wall BC occur for higher inclination angle (φ = 75°).     

Investigation of coated tubes in cross-flow boiling

The boiling in cross-flow is investigated for coated tubes (low-porosity, flame-sprayed) in this paper. The effect of surface roughness on flow boiling heat transfer for a horizontal tube surface in cross-flow is studied for saturated boiling of water at atmospheric pressure. The parameters varied were for flow velocity up to 3.24 kg/s (G = 258.49 kg/m² s), heat flux from 12 to 45 kW/m², surface roughness (Ra) from 0.3296 to 4.731 μm. Nominal enhancement in heat transfer coefficient at higher mass flux may be attributed to the continued nucleation at the uppermost surfaces (in the wake region of the flow) of the rougher tubes thereby increasing the overall heat transfer rate. The flow boiling data was found to best fit the Kutateladze asymptotic equation h = h_l[1 + (h_npb/h_l)ⁿ]^1/n with the value of n = 2.258 (which is close to the value of n = 2 suggested by Kutateladze).     

Experimental study of R-134a evaporation heat transfer in a narrow annular duct

An experiment is carried out here to investigate the evaporation heat transfer and associated evaporating flow pattern for refrigerant R-134a flowing in a horizontal narrow annular duct. The gap of the duct is fixed at 1.0 and 2.0 mm. In the experiment, the effects of the duct gap, refrigerant vapor quality, mass flux and saturation temperature and imposed heat flux on the measured evaporation heat transfer coefficient h_r are examined in detail. For the duct gap of 2.0 mm, the refrigerant mass flux G is varied from 300 to 500 kg/m² s, imposed heat flux q from 5 to 15 kW/m², vapor quality x_m from 0.05 to 0.95, and refrigerant saturation temperature T_sat from 5 to 15 °C. While for the gap of 1.0 mm, G is varied from 500 to 700 kg/m² s with the other parameters varied in the same ranges as that for δ = 2.0 mm. The experimental data clearly show that the evaporation heat transfer coefficient increases almost linearly with the vapor quality of the refrigerant and the increase is more significant at a higher G. Besides, the evaporation heat transfer coefficient also rises substantially at increasing q. Moreover, a significant increase in the evaporation heat transfer coefficient results for a rise in T_sat, but the effects are less pronounced in the narrower duct at a low imposed heat flux and a high refrigerant mass flux. Furthermore, the evaporation heat transfer coefficient increases substantially with the refrigerant mass flux except at low vapor quality. We also note that reducing the duct gap causes a significant increase in h_r. In addition to the heat transfer data, photos of R-134a evaporating flow taken from the duct side show the change of the dominant two-phase flow pattern in the duct with the experimental parameters. Finally, an empirical correlation for the present measured heat transfer coefficient for the R-134a evaporation in the narrow annular ducts is proposed.     

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.     

The performance evaluation of Al2O3/water nanofluid flow and heat transfer in microchannel heat sink

The numerical modeling of the conjugate heat transfer and fluid flow of Al₂O₃/water nanofluid through the microchannel heat sink is presented in the paper. The laminar flow regime was considered along with viscous dissipation effect. The microchannel heat sink with square microchannels and D_h = 50 μm is considered. The heat flux was fixed to q = 35 W/m² with heating and cooling cases. The water based Al₂O₃ nanofluid was encountered with various volume concentrations of Al₂O₃ particles $?=1''–''9\%$ and three diameters of the particle d_p = 13, 28 and 47 nm. The analysis is performed on the results obtained for the local heat transfer coefficients based on a fixed pumping power. The results reveal a different local heat transfer behavior compared to the analysis made on a basis of the constant Re.     

Numerical simulation of turbulent natural convection in a rectangular enclosure having finite thickness walls

Turbulent natural convection in a rectangular enclosure having finite thickness heat-conducting walls at local heating at the bottom of the cavity provided that convective-radiative heat exchange with an environment on one of the external borders has been numerically studied. Mathematical simulation has been carried out in terms of the dimensionless Reynolds averaged Navier–Stokes (RANS) equations in stream function–vorticity formulations. The formulation comprises the standard two equation k–ε turbulence model with wall functions, along with the Boussinesq approximation, for the flow and heat transfer. The special attention was paid to the effects of the Grashof number 10⁸ ? Gr < 10¹⁰, the transient factor 0 < τ < 1000 and the thermal conductivity ratio k_2,1 = 5.7 × 10^?4, 6.8 × 10^?5 both on local and on integral problem parameters. Detailed results including stream lines, temperature profiles and correlation for the average Nusselt number in terms of Grashof number have been obtained.     

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.     

Flow and heat transfer across a confined square cylinder in the steady flow regime: Effect of Peclet number

The flow and heat transfer characteristics of an isolated square cylinder in crossflow placed symmetrically in a planar slit have been investigated for the range of conditions as 1 ⩽ Re ⩽ 45, 0.7 ⩽ Pr ⩽ 4000 (Pe ⩽ 4000) and β = 1/8, 1/6 and 1/4. Heat transfer correlations have been obtained in the steady flow regime for the constant temperature and constant heat flux boundary conditions on the solid square cylinder in crossflow. In addition, variation of the local Nusselt number on each face of the obstacle and representative isotherm plots are presented to elucidate the role of Prandtl number and blockage ratio on drag coefficient and heat transfer.     

CFD analysis of convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer

Steady Reynolds-Averaged Navier–Stokes (RANS) CFD is used to evaluate the forced convective heat transfer at the surfaces of a cube immersed in a turbulent boundary layer, for applications in atmospheric boundary layer (ABL) wind flow around surface-mounted obstacles such as buildings. Two specific configurations are analysed. First, a cube placed in turbulent channel flow at a Reynolds number of 4.6 × 10³ is considered to validate the numerical predictions by comparison with wind-tunnel measurements. The results obtained with low-Reynolds number modelling (LRNM) show a satisfactory agreement with the experimental data for the windward surface. Secondly, a cube exposed to high-Reynolds number ABL flow is considered. The heat transfer in the boundary layer is analysed in detail. The dimensionless parameter y¹, which takes into account turbulence, is found to be more appropriate for evaluating heat transfer than the commonly used y⁺ value. Standard wall functions, which are frequently used for high-Reynolds number flows, overestimate the convective heat transfer coefficient (CHTC) significantly (±50%) compared to LRNM. The distribution of the CHTC–U₁₀ correlation over the windward surface is reported for Reynolds numbers of 3.5 × 10⁴ to 3.5 × 10⁶ based on the cube height and U₁₀, where U₁₀ is the wind speed in the undisturbed flow at a height of 10 m.     

Mixed convection from an open cavity in a horizontal channel

Heat transfer results for mixed convection from a bottom heated open cavity subjected to an external flow are reported in this study for a wide range of the governing parameters (i.e., 1 ≤ Re ≤ 2000, 0 ≤ Gr ≤ 10⁶) over cavities with various aspect ratios (A = 0.5, 1, 2 and 4). It has been found that the Reynolds number and Garshof number control the flow pattern and the occurrence of recirculating cells while the aspect ratio has a significant influence on the orientation of these cells. Heat transfer from the cavity base approaches that of natural convection at a low Reynolds number (i.e., the asymptotic natural convection regime) and approaches that of forced convection at a high Reynolds number (i.e., the asymptotic forced convection regime). In the mixed convection regime, the heat transfer rate is reduced and the flow may become unstable. A unique heat transfer correlation which covers all three convection regimes is also presented.     

Numerical study of hydrodynamic and heat transfer of nanofluid flow in microchannels containing micromixer

In this study heat transfer and fluid flow of Al₂O₃/water nanofluid in two dimensional parallel plate microchannel without and with micromixers have been investigated for nanoparticle volume fractions of ϕ = 0, ϕ = 4%  and base fluid Reynolds numbers of Re_f = 5, 20, 50. One baffle on the bottom wall and another on the top wall work as a micromixer and heat transfer enhancement device. A single-phase finite difference FORTRAN code using Projection method has been written to solve governing equations with constant wall temperature boundary condition. The effect of various parameters such as nanoparticle volume fraction, base fluid Reynolds number, baffle distance, height and order of arrangement have been studied. Results showed that the presence of baffles and also increasing the Re number and nanoparticle volume fraction increase the local and averaged heat transfer and friction coefficients. Also, the effect of nanoparticle volume fraction on heat transfer coefficient is more than the friction coefficient in most of the cases. It was found that the main mechanism of enhancing heat transfer or mixing is the recirculation zones that are created behind the baffles. The size of these zones increases with Reynolds number and baffle height. The fluid pushing toward the wall by the opposed wall baffle and reattaching of separated flow are the locations of local maximum heat transfer and friction coefficients.     

Heat transfer study in a discoidal system: The influence of rotation and space between disks

This article presents an experimental study of the local heat transfer on the rotor surface in a discoidal rotor–stator system air-gap in which an air jet comes through the stator and impinges the rotor. To determine the surface temperatures, measurements were taken on the rotor, using an experimental technique based on infrared thermography. A thermal balance was used to identify the local convective heat transfer coefficient. The influence of the dimensionless spacing interval G between the disks and of the rotational Reynolds number Re was measured and compared with the data available in bibliography. Local convective heat transfer coefficients were obtained for an axial Reynolds number Re_j = 41.6 × 10³, a rotational Reynolds number Re between 0.2 × 10⁵ and 5.16 × 10⁵, and a dimensionless spacing interval G ranging from 0.01 to 0.16.