Numerical investigation of heat transfer enhancement on a porous vortex-generator applied to a block-heated channel

The numerical simulation is used to obtain the unsteady laminar flow and convective heat transfer in the block-heated channel with the porous vortex-generator. The general Darcy–Brinkman–Forchheimer model is adopted for the porous vortex-generator. The parameters studies including porosity, Darcy number, width-to-height ratio of porous vortex-generator and Reynolds number have been explored on heat transfer enhancement and vortex-induced vibration in detail. The results indicate that heat transfer enhancement and vortex-induced vibration increase with increasing Reynolds number and width-to-height ratio. However, the porosity has slight influence on heat transfer enhancement and vortex-induced vibration. When Darcy number is 10^?3 or 10^?4, installing a porous vortex-generator with B/h = 1.0 improves overall heat transfer the best along heated blocks, and has a strong reduction of vortex-induced vibration.     

Numerical study of transient conjugate heat transfer of a turbulent impinging jet

This study presents the numerical study of transient conjugate heat transfer in a high turbulence air jet impinging over a flat circular disk. The numerical simulation of transient, two-dimensional cylindrical coordinate, turbulent flow and heat transfer is adopted to test the accuracy of the theoretical model. The turbulent governing equations are resolved by the control-volume based finite-difference method with a power-low scheme, and the well-known low-Re κ–ω turbulence model to describe the turbulent structure. The SIMPLE algorithm is adopted to solve the pressure–velocity coupling. The parameters studied include turbulent flow Reynolds number (Re = 16,100–29,600), heated temperature of a circular disk (T_h = 373 K) or heat flux (q″ = 63–189 kW/m²), and orifice to heat-source spacing (H/D = 4–10). The numerical results of the transient impinging process indicate that the jet Reynolds number has a significant effect on the hydrodynamics and heat transfer, particularly in the stagnation region of an impinging jet. High turbulence values lead to greater heat transfer coefficients in the stagnation region and cause a bypass of the laminar-to-turbulent transition region in the wall jet region. Induced turbulence from the environment around the jet also influences the variation of the stagnation heat transfer. The modeling approach used here effectively captures both the stagnation region behavior and the transition to turbulence, thus forming the basis of a reliable turbulence model.     

Numerical investigation of pressure drop reduction without surrendering heat transfer enhancement in partially porous channel

The present study is to investigate the numerical simulation of steady laminar forced convection in a partially porous channel, with four dissimilar porous-blocks, attached to the strip heat sources at the bottom wall. The analysis is based on the Navier–Stokes equation in the fluid field, the Darcy–Brinkman–Forchheimer flow model in the porous field, and the energy equations for two thermal fields. The effects of variations of different parameters such as porous blocks Darcy numbers, arrangements of dissimilar blocks, Forchheimer coefficient, Reynolds number, thermal conductivity and Prandtl number are investigated and the velocity and temperature fields are presented and discussed. In the dissimilar partially porous channel, it is found that when the blocks sorted from the lowest to the highest Da in the flow direction, the total heat transfer enhancement is almost the same as in the similar porous channel (Nu/Nu_sim = 92%), while the total pressure drop is considerably lower (P/P_sim = 28%). In addition, reverse arrangement of porous blocks is suggested to prepare more uniform temperature gradient in all heat sources.     

The effect of confinement on the flow and turbulent heat transfer in a mist impinging jet

Numerical study of the effect of confinement on a flow structure and heat transfer in an impinging mist jets with low mass fraction of droplets (M_L1 ? 1%) were presented. The turbulent mist jet is issued from a pipe and strikes into the center of the flat heated plate. Mathematical model is based on the steady-state RANS equations for the two-phase flow in Euler/Euler approach. Predictions were performed for the distances between the nozzle and the target plate x/(2R) = 0.5–10 and the initial droplets size (d₁ = 5–100 μm) at the varied Reynolds number based on the nozzle diameter, Re = (1.3–8) × 10⁴. Addition of droplets causes significant increase of heat transfer intensity in the vicinity of the jet stagnation point compared with the one-phase air impinging jet. The presence of the confinement upper surface decreases the wall friction and heat transfer rate, but the change of friction and heat transfer coefficients in the stagnation point is insignificant. The effect of confinement on the heat transfer is observed only in very small nozzle-to-plate distances (H/(2R) < 0.5) both in single-phase and mist impinging jets.     

Performance assessment in a heat exchanger tube with alternate clockwise and counter-clockwise twisted-tape inserts

The article presents an experimental study of turbulent heat transfer and flow friction characteristics in a circular tube equipped with two types of twisted tapes: (1) typical twisted tapes and (2) alternate clockwise and counterclockwise twisted tapes (C–CC twisted tapes). Nine different C–CC twisted tapes are tested in the current work; they included the tapes with three twist ratios, y/w = 3.0, 4.0 and 5.0, each with three twist angles, θ = 30^o, 60^o and 90^o. The experiments have been performed over a Reynolds number range of 3000–27,000 under uniform heat flux conditions, using water as working fluid. The obtained results reveal that the C–CC twisted-tapes provide higher heat transfer rate, friction factor and heat transfer enhancement index than the typical twisted-tapes at similar operating conditions. The results also show that the heat transfer rate of the C–CC tapes increases with the decrease of twist ratio and the increase of twist angle values. Depending on Reynolds number, twist ratio and twist angle values, the mean Nusselt numbers in the tube fitted with the C–CC twisted tapes are higher than those with the typical ones and the plain tube around 12.8–41.9% and 27.3–90.5%, respectively. The maximum heat transfer enhancement indexes of the C–CC twisted tapes with θ = 90^o for y/w = 3.0, 4.0 and 5.0, are 1.4, 1.34 and 1.3, respectively. In addition, correlations of the Nusselt number and the friction factor for using the C–CC twisted tapes are also determined. Both predicted Nusselt number and friction factor are within ±15% and ±15% deviation compared to the experimental data.     

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.     

Measurement of local heat/mass transfer coefficients on a dimple using naphthalene sublimation

Local and average heat/mass transfer characteristics on a single dimple were investigated using a naphthalene sublimation technique. The dimple depth in this study ranged from 20% to 40% of the channel height. The experimental conditions covered the range from laminar to low-velocity turbulent flow regimes, 500 ? Re_H ? 5000. Secondary flows from the dimple were clearly observed in the transient flow regime of Re_H = 2000–3000. The velocity fluctuation in the mixing layer over the dimple increased with the dimple depth and the Reynolds number. The impingement of the mixing layer and the induced secondary flows augmented the Sherwood number around the rear rim of the dimple and in the rear plateau region, respectively. For a Reynolds number of 3000, the Sherwood number increased significantly due to the increased fluctuation in the mixing layer and the intensified secondary flows from the dimple. The heat/mass transfer augmentation factors increased as the Reynolds number increased, reaching 1.5 at a Reynolds number of 5000.     

Experimental studies on heat transfer and friction factor characteristics of turbulent flow through a circular tube with small pipe inserts

Heat transfer performance and pressure drop tests were performed on a circular tube with small pipe inserts. These inserts with different spacer lengths (S = 100, 142.9 and 200 mm) and arc radii (R = 5, 10 and 15 mm) were tested at Reynolds numbers between 4000 and 18,000. Tap water was used as working fluid. The use of pipe inserts allowed for a high heat transfer coefficient with relatively low flow resistance. The Nusselt number and friction factor increase with the decrease in spacer length. Optimal results were obtained for S = 100 mm (R = 10 mm). Heat transfer rates and friction factors were enhanced by 2.09–2.67 and 1.59–1.85 times, respectively, to those in the plain tube. Performance evaluation criterion (PEC) values were approximately 1.79–2.17. The Nusselt number and friction factor increase with the decrease in arc radius. Small pipe inserts with R = 5 mm and S = 100 mm show maximal heat transfer rates of 2.61–3.33 and friction factors of 1.6–1.8 times those of the empty tube. The PEC values were 2.23–2.7. Compared with other inserts, pipe inserts can transfer more heat for the same pumping power for their unique structure.     

A porous model for the square pin-fin heat sink situated in a rectangular channel with laminar side-bypass flow

This work illustrates the compact heat sink simulations in forced convection flow with side-bypass effect. Conventionally, the numerical study of the fluid flow and heat transfer in finned heat sinks employs the detailed model that spends a lot of computational time. Therefore, some investigators begin to numerically study such problem by using the compact model (i.e. the porous approach) since the regularly arranged fin array can be set as a porous medium. The computations of the porous approach model will be faster than those of the detailed mode due to the assumption of the volume-averaging technique. This work uses the Brinkman–Forchheimer model for fluid flow and two-equation model for heat transfer. A configuration of in-line square pin-fin heat sink situated in a rectangular channel with fixed height (H = 23.7 mm), various width and two equal-spacing bypass passages beside the heat sink is successfully studied. The pin-fin arrays with various porosities (ε = 0.358–0.750) and numbers of pin-fins (n = 25–81), confined within a square spreader whose side length (L) is 67 mm, are employed. The numerical results suggest that, within the range of present studied parameters (0.358 ? ε ? 0.750, 25 ? n ? 81 and 1 ? W/L ? 5), the pin-fin heat sink with ε = 0.750 and n = 25 is the optimal cooling configuration based on the maximum ratio of Nusselt number to dimensionless pumping power (Nu/(ΔP × Re³)). Besides, based on medium Nu/(ΔP × Re³) value and suitable channel size, W/L = 2–3 is suggested as the better size ratio of channel to heat sink.     

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.     

The investigation of groove geometry effect on heat transfer for internally grooved tubes

An experimental study of surface heat transfer and friction characteristics of a fully developed turbulent air flow in different grooved tubes is reported. Tests were performed for Reynolds number range 10,000–38,000 and for different geometric groove shapes (circular, trapezoidal and rectangular). The ratio of tube length-to-diameter is 33. Among the grooved tubes, heat transfer enhancement is obtained up to 63% for circular groove, 58% for trapezoidal groove and 47% for rectangular groove, in comparison with the smooth tube at the highest Reynolds number (Re = 38,000). Correlations of heat transfer and friction coefficient were obtained for different grooved tubes. In evaluation of thermal performance, it is seen that the grooved tubes are thermodynamically advantageous (Ns, a < 1) up to Re = 30,000 for circular and trapezoidal grooves and up to Re = 28,000 for rectangular grooves. It is observed that there is an optimum value of the entropy generation number at about Re = 17,000 for all investigated grooves.     

Unsteady conjugate natural convection in a square enclosure filled with a porous medium

Mathematical simulation of unsteady natural convection modes in a square cavity filled with a porous medium having finite thickness heat-conducting walls with local heat source in conditions of heterogeneous heat exchange with an environment at one of the external boundaries has been carried out. Numerical analysis was based on Darcy–Forchheimer model in dimensionless variables such as a stream function, a vorticity vector and a temperature. The special attention was given to analysis of Rayleigh number effect Ra = 10⁴, 10⁵, 10⁶, of Darcy number effect Da = 10^?5, 10^?4, 10^?3, ∞, of the transient factor effect 0 < τ < 1000 and of the heat conductivity ratio k_2,1 = 3.7 × 10^?2, 5.7 × 10^?4, 6.8 × 10^?5 on the velocity and temperature fields. The influence scales of the defining parameters on the average Nusselt number have been detected.     

Enhancement of an impingement heat transfer between turbulent mist jet and flat surface

The work presents the results of numerical investigation of the flow structure and heat transfer of impact mist jet with low concentration of droplets (M_L1 ? 1%). The downward gas-droplets jet issued from a pipe and strikes into at a center of the circular target wall. Mathematical model is based in the solution to RANS equations for the two-phase flow in Euler approximation. For the calculation of the fluctuation characteristics of the dispersed phase equations of Zaichik et al. (1997) [35] model were applied. Predictions were performed for the distances between the nozzle and target plate x/(2R) = 1–10 and the initial droplets size (d₁ = 5–100 μm) at the fixed Reynolds number based on the nozzle diameter, Re = 26,600. Addition of droplets causes significant increase of heat transfer intensity in the vicinity of the jet stagnation point compared with the one-phase air impact jet.     

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 study of mixed convection within porous square cavities using Bejan’s heatlines: Effects of thermal aspect ratio and thermal boundary conditions

The present numerical study deals with mixed convection flows within square enclosures filled with porous media. The influence of various thermal boundary conditions on bottom and side walls based on thermal aspect ratio (A) is investigated for a wide range of parameters (1 ? Re ? 100, 0.015 ? Pr ? 7.2, 10^?5 ? Da ? 10^?3 and 10³ ? Gr ? 10⁵). A penalty finite element method with bi-quadratic elements has been used to investigate the results in terms of streamlines, isotherms and heatlines and average Nusselt numbers. Lid driven effect is dominant at low Darcy number (Da = 10^?5), whereas buoyancy driven effect is dominant at high Darcy numbers (Da = 10^?4 and Da = 10^?3) for Re = 1. Asymmetric pattern is observed in isotherms and heatlines for Re = 100. It is found that thermal gradient is high at the center of the bottom wall for A = 0.1 due to large dense heatlines at that zone and that is low for A = 0.9 irrespective of Re, Pr and Gr. Overall heat transfer rates are higher for A = 0.1 compared to other thermal aspect ratios (A = 0.5, A = 0.9) irrespective of Darcy number, Prandtl number and Reynolds number.     

Forced convective heat transfer with impinging slot jets of meso-scale

Measurements were made to investigate the localized heat transfer behavior of submerged slot jets. The experiments were performed with kerosene jets impinging on a vertical constant-heat-flux surface from a meso-scale slot nozzle 125 μm in width with Re = 600–1200 and nozzle-to-plate spacing Z/B = 2–20. Heat transfer coefficients at the stagnation line were measured and correlated as a function of jet Reynolds numbers and Prandtl numbers. Lateral distributions of local heat transfer coefficients were also determined and correlated. Non-monotonic variations and unusual behavior of local heat transfers were observed and attributed to the possible transition from a laminar to a turbulent flow. This transition takes place within an extremely short distance of 400–500 μm.     

Heat transfer from a spinning disk during semi-confined axial impingement from a rotating nozzle

Conjugate heat transfer from a uniformly heated spinning solid disk of finite thickness and radius during a semi-confined liquid jet impingement from a rotating nozzle is studied. The model covers the entire fluid region including the impinging jet on a flat circular disk and flow spreading out downstream under the spinning confinement plate and free surface flow after exposure to the ambient gaseous medium. The model examines how the heat transfer is affected by adding a secondary rotational flow under semi-confined jet impingement. The solution is made under steady state and laminar conditions. The study considered various plate materials such as aluminum, copper, silver, constantan and silicon. Ammonia, water, flouroinert FC-77 and MIL-7808 oil were used as working fluids. The range of parameters covered included Reynolds number (220–900), Ekman number (7.08 × 10^?5–∞), nozzle-to-target spacing (β = 0.25–1.0), disk thicknesses to nozzle diameter ratio (b/d_n = 0.25–1.67), Prandtl number (1.29–124.44) and solid to fluid thermal conductivity ratio (36.91–2222). It was found that a higher Reynolds number increased local heat transfer coefficient reducing the interface temperature difference over the entire disk surface. The rotational rate also increased local heat transfer coefficient under most conditions. An engineering correlation relating the Nusselt number with other dimensionless parameters was developed for the prediction of the system performance.     

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.     

Unsteady fluid mechanics and heat transfer study in a double-tube air–combustor heat exchanger with porous medium

Fluid mechanics and heat transfer are studied in a double-tube heat exchanger that uses the combustion gases from natural gas in a porous medium located in a cylindrical tube to warm up air that flows through a cylindrical annular space. The mathematical model is constructed based on the equations of continuity, linear momentum, energy and chemical species. Unsteady fluid mechanics and heat transfer by forced gas convection in the porous media, with combustion in the inner tube, coupled to the forced convection of air in the annular cylindrical space are predicted by use of finite volumes method. Numerical simulations are made for four values of the annular air flow Reynolds number in the range 100 ? Re ? 2000, keeping constant the excess air ψ = 4.88, the porosity ε = 0.4, and the air–fuel mixture inlet speed Uo = 0.43 m/s. The results obtained allow the characterization of the velocity and temperature distributions in the inner tube and in the annular space, and at the same time to describe the displacement of the moving combustion zone and the annular porous media heat exchanger thermal efficiency. It is concluded that the temperature increase is directly related to the outer Reynolds number.     

Numerical study of the heat sink with un-uniform fin width designs

The thermal performances of the heat sink with un-uniform fin width designs with an impingement cooling were investigated numerically. The governing equations are discretized by using a control-volume-based finite-difference method with a power-law scheme on an orthogonal non-uniform staggered grid. The coupling of the velocity and the pressure terms of momentum equations are solved by the SIMPLEC algorithm. The well-known k ? ε two-equations turbulence model is employed to describe the turbulent structure and behavior. The parameters include the five Reynolds number (Re = 5000–25000), three fin heights (H = 35, 40, 45 mm), and five fin width designs (Type-1–Type-5). The objective of this study is to examine the effects of the fin shape of the heat sink on the thermal performance. The results show that the Nusselt number increases with the Reynolds number. The increment of the Nusselt number decreases gradually with the increasing Reynolds number. Furthermore, the effects of fin dimensions on the Nusselt number at high Reynolds numbers are more significant than that at low Reynolds numbers. It is also found that there is potential for optimizing the un-uniform fin width design.