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.     

Numerical simulation of natural convection heat transfer in a porous cavity heated from below using a non-Darcian and thermal non-equilibrium model

The present paper investigates the numerical simulation of steady laminar incompressible natural convection heat transfer in an enclosed cavity that is filled with a fluid-saturated porous medium. The bottom wall is subjected to a relatively higher temperature than the top wall while the vertical walls are considered to be insulated. The flow field is modeled upon incorporating different non-Darcian effects, such as the convective term, Brinkman effect and Forchhiemer quadratic inertial effect. Moreover the two-equation model is used to separately account for the local fluid and solid temperatures. The numerical solution is obtained through the application of the finite volume method. The appraisals of the sought objectives are performed upon identifying key dimensionless groups of parameters. These dimensionless groups along with their operating domains are: Rayleigh number 1 ⩽ Ra ⩽ 400, Darcy number 10⁻⁴ ⩽ Da ⩽ 10⁻³, effective fluid-to-solid thermal conductivity ratio 0.1 ⩽ κ ⩽ 1.0, and the modified Biot number 1 ⩽ χ ⩽ 100. The non-Darcian effects are first examined over a broad range of Rayleigh number. Next, the implications of the group of parameters on the flow circulation intensity, local thermal non-equilibrium (LTNE) and average Nusselt number are highlighted and pertinent observations are documented.     

Jet impingement cooling of a horizontal surface in an unconfined porous medium: Mixed convection regime

Two-dimensional slot jet impingement cooling of an isothermal horizontal surface immersed in an unconfined porous medium is simulated numerically to gain insight into thermal characteristics under mixed convection conditions with the limitation of the Darcy model. The jet direction is considered to be perpendicular from the top to the horizontal heated element; therefore, the jet flow and the buoyancy driven flow are in opposite directions. The results are presented in the mixed convection regime with wide ranges of the governing parameters: Péclet number (1 ? Pe ? 1000), Rayleigh number (10 ? Ra ? 100), half jet width (0.1 ? D ? 0.5), and the distance between the jet and the heated portion (0.1 ? H ? 1.0). It is found that the average Nusselt number increases with increase in either Rayleigh number or jet width for high values of Péclet number. The average Nusselt number also increases with decrease in the distance between the jet and the heated portion. It is shown that mixed convection mode can cause minimum average Nusselt number at two values of Péclet number and a maximum average Nusselt number occurs in between theses two Péclet numbers at higher Rayleigh number due to counteraction of jet flow against buoyancy driven flow. Hence careful consideration must be given while designing a system of jet impingement cooling through porous medium.     

Jet impingement cooling of a constant heat flux horizontal surface in a confined porous medium: Mixed convection regime

In the present study, numerical investigation of jet impingement cooling of a constant heat flux horizontal surface immersed in a confined porous channel is performed under mixed convection conditions with the limitation of the Darcy model. The results are presented in the mixed convection regime with wide ranges of the governing parameters: Péclet number (1 ? Pe ? 1000), Rayleigh number (10 ? Ra ? 100), half jet width (0.1 ? D ? 1.0), and the distance between the jet and the heated portion (0.1 ? H ? 1.0). It is found that the average Nusselt number increases with increase in either Rayleigh number or jet width for high values of Péclet number. The average Nusselt number also increases with decrease in the distance between the jet and the heated portion. The correlation for Nu_avg in the forced convection regime is suggested. It is shown that mixed convection mode can cause minimum average Nusselt number unfavorably due to counteraction of jet flow against buoyancy driven flow. Hence, careful consideration must be given while designing a system of jet impingement cooling through porous medium.     

Flow structures and heat transfer of swirling jet impinging on a flat surface with micro-vibrations

This study focuses on the changes in the flow characteristics of a round jet issuing from a straight tube inserted with longitudinal swirling strips and impinging on a constant-heat-flux flat surface undergoing forced vibrations in the vertical plane. Smoke flow visualization is used to investigate the nature of the complicated flow phenomena under the swirling-flow jet for this impingement cooling. Effects of flow Reynolds number (440 ⩽ Re ⩽ 27 000), the geometries of the nozzle (BR, LSS and CSS), jet-to-test plate placement (3 ⩽ H/d ⩽ 16), and surface vibration frequencies, f [0.3–10.19 Hz (the relative amplitude of the flat surface ranged from 0.5 to 8.1 mm)] are examined. In addition, correlations were developed to predict the Nusselt number for the vibration using the results of Wen and Jang [An impingement cooling on a flat surface by using circular jet with longitudinal swirling strips, Int. J. Heat Mass Transfer 46 (2003) 4657–4667] for the no-vibration case of the present study.     

A numerical study on the forced convection heat transfer from an isothermal and isoflux sphere in the steady symmetric flow regime

The effects of Reynolds and Prandtl numbers on the heat transfer characteristics of an unconfined sphere for different thermal boundary conditions (isothermal and isoflux) on the sphere surface have been investigated numerically by using a finite volume method for the range of conditions as 5 ⩽ Re ⩽ 200 and 0.7 ⩽ Pr ⩽ 400 (the maximum value of Peclet number being 2000). Based on the numerical results obtained herein, heat transfer correlations are developed for the constant temperature and the constant heat flux boundary conditions on the solid sphere surface in the steady symmetric flow regime. The variation of local Nusselt number on the sphere surface shows the effect of Prandtl number on heat transfer from a sphere in this flow regime. In addition, this work also demonstrates an approach to solve such flow problems using the Cartesian form of the field equations.     

Steady natural convection flows in a square cavity with linearly heated side wall(s)

The present numerical study deals with natural convection flow in a closed square cavity when the bottom wall is uniformly heated and vertical wall(s) are linearly heated whereas the top wall is well insulated. Non-linear coupled PDEs governing the flow have been solved by penalty finite element method with bi-quadratic rectangular elements. Numerical results are obtained for various values of Rayleigh number (Ra) (10³ ⩽ Ra ⩽ 10⁵) and Prandtl number (Pr) (0.7 ⩽ Pr ⩽ 10). Results are presented in the form of streamlines, isotherm contours, local Nusselt number and the average Nusselt as a function of Rayleigh number.     

Heat transfer of impinging jet-array over convex-dimpled surface

A detailed heat transfer measurement over a convex-dimpled surface of impinging jet-array with three eccentricities (E/H) between jet-centre and dimple-centre is performed. These surface dimples considerably modify heat transfers from smooth-walled scenarios due to different impinging topologies for jet array with modified inter-jet reactions. Heat transfer variations caused by adjusting jet Reynolds number (Re) and separation distance (S/D_j) over the ranges of 5000 ⩽ Re ⩽ 15,000 and 0.5 ⩽ S/D_j ⩽ 11 with three eccentricities of E/H = 0, 1/4 and 1/2 are examined. A selection of experimental data illustrates the isolated and interactive influences of Re, S/D_j and E/H on local and spatially averaged heat transfers. In conformity with the experimentally revealed heat transfer physics, a regression-type analysis is performed to generate a set of heat transfer correlations, which permit the evaluations of spatially averaged Nusselt numbers over central jet region of dimpled impinging surface.     

Natural convection in a square cavity filled with a porous medium: Effects of various thermal boundary conditions

Natural convection flows in a square cavity filled with a porous matrix has been studied numerically using penalty finite element method for uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls. Darcy–Forchheimer model is used to simulate the momentum transfer in the porous medium. The numerical procedure is adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁶, Darcy number Da, 10⁻⁵ ⩽ Da ⩽ 10⁻³, and Prandtl number Pr, 0.71 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous thermal boundary conditions. Numerical results are presented in terms of stream functions, temperature profiles and Nusselt numbers. Non-uniform heating of the bottom wall produces greater heat transfer rate at the center of the bottom wall than uniform heating case for all Rayleigh numbers but average Nusselt number shows overall lower heat transfer rate for non-uniform heating case. It has been found that the heat transfer is primarily due to conduction for Da ⩽ 10⁻⁵ irrespective of Ra and Pr. The conductive heat transfer regime as a function of Ra has also been reported for Da ⩾ 10⁻⁴. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes the power law correlations between average Nusselt number and Rayleigh numbers are presented.     

Estimating effective thermal conductivity of two-phase materials

In this paper, correlations are proposed to estimate the effective thermal conductivity of two-phase materials. For any α, Maxwell equation for 0.0 < c ⩽ 0.10 and phase inverted Maxwell for 0.9 ⩽ c ⩽ 1 are considered. For concentrations between 10% and 90%, and low α (<20), an equation based on the unit-cell approach (constant isotherms) is proposed. For α > 20, three correlations are proposed based on field solution approach which includes three α ranges viz. medium (20 ⩽ α ⩽ 100), high (100 ⩽ α ⩽ 1000) and very high (1000 < α). The predicted effective thermal conductivity of two-phase system is compared with well-established models. Comparison of the predicted values of the correlations with experimental results is also made. The predictions of effective thermal conductivity of two-phase materials match well with the experimental values.     

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.     

Effects of thermal boundary conditions on natural convection flows within a square cavity

A numerical study to investigate the steady laminar natural convection flow in a square cavity with uniformly and non-uniformly heated bottom wall, and adiabatic top wall maintaining constant temperature of cold vertical walls has been performed. A penalty finite element method with bi-quadratic rectangular elements has been used to solve the governing mass, momentum and energy equations. The numerical procedure adopted in the present study yields consistent performance over a wide range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁵ and Prandtl number Pr, 0.7 ⩽ Pr ⩽ 10) with respect to continuous and discontinuous Dirichlet boundary conditions. Non-uniform heating of the bottom wall produces greater heat transfer rates at the center of the bottom wall than the uniform heating case for all Rayleigh numbers; however, average Nusselt numbers show overall lower heat transfer rates for the non-uniform heating case. Critical Rayleigh numbers for conduction dominant heat transfer cases have been obtained and for convection dominated regimes, power law correlations between average Nusselt number and Rayleigh numbers are presented.     

The effect of inclination on impinging jets at small nozzle-to-plate spacing

The effect of inclination on heat transfer characteristics of an impinging slot air jet is experimentally investigated. The effects of inclination angle (0° ? θ ? 40°) and dimensionless pumping power on the Nusselt number are considered. The focus is on cases where the nozzle-to-plate spacing is equal to or less than one nozzle diameter (H/d_h ? 1.0). The results show that the heat transfer characteristics of small nozzle-to-plate spacings are significantly different from those of large nozzle-to-plate spacings. In the cases of fixed flow rate conditions, the impingement point and average Nusselt numbers at small nozzle-to-plate spacing (H/d_h ? 1.0) increase as the inclination angle increases due to an increase in the pumping power, while the impingement point and average Nusselt numbers at large nozzle-to-plate spacing (H/d_h > 1.0) decrease as the inclination angle increases due to momentum loss of the wall jet. In the cases of fixed pumping power conditions, the impingement point and average Nusselt numbers at both of small and large nozzle-to-plate spacings are independent of the inclination angle. Based on the experimental results, correlations for the impingement point and average Nusselt numbers of the impinging jet are suggested as a function of the pumping power alone.     

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.     

Heat transfer of rectangular narrow channel with two opposite scale-roughened walls

An experimental study of heat transfer and pressure drop in a rectangular channel roughened by scaled surfaces on two opposite walls with flows directed in the forward and downward directions for Reynolds numbers (Re) in the range of 1500 ⩽ Re ⩽ 15,000 was performed. Nusselt number ratios between the scale-roughened and smooth-walled ducted flows (Nu/Nu_∞) were in the range of 7.4–9.2 and 6.2–7.4 for laminar forward and downward flows respectively. The Nu/Nu_∞ values for turbulent developed flows in the scale-roughened channel with forward and downward flows were about 4.5 and 3 respectively. A comparison of present data with reported results using different types of surface roughness demonstrated the better thermal performances of present scale-roughened channel with forward flow at conditions of Re > 10,000. Experimental correlations of heat transfer and friction coefficient were derived for the present scale-roughened rectangular channel.     

Analysis of free convection about a horizontal cylinder in a porous media using a thermal non-equilibrium model

The thermally non-equilibrium model is used to study the free convection from a horizontal cylinder immersed in porous media. The governing equations are transformed to dimensionless form by introducing the boundary layer dimensionless variables. The resultant parabolic system of differential equations is solved by using an implicit finite difference method based on Keller box algorithm. The results of the developed code are validated with different mesh sizes, which can be used as benchmark results for thermally equilibrium condition. Numerical results are obtained for non-equilibrium model to analyze the effect of the governing parameters, which are the heat transfer coefficient between the solid and fluid phases H and the porosity scaled thermal conductivity ratio K_r. The results show that increasing H or K_r leads to increase in the total average Nusselt number. The value of the average Nusselt number for both fluid and solid phases as well as the total average Nusselt number have approached the corresponding values of the thermally equilibrium model at high value of H × K_r.     

Effect of nanofluid variable properties on natural convection in enclosures filled with a CuO–EG–Water nanofluid

This work focuses on the study of natural convection heat transfer characteristics in a differentially-heated enclosure filled with a CuO–EG–Water nanofluid for different published variable thermal conductivity and variable viscosity models. The problem is given in terms of the vorticity–stream function formulation and the resulting governing equations are solved numerically using an efficient finite-volume method. Comparisons with previously published work are performed and the results are found to be in good agreement. Various results for the streamline and isotherm contours as well as the local and average Nusselt numbers are presented for a wide range of Rayleigh numbers (Ra = 10³–10⁵), volume fractions of nanoparticles (0 ≤ φ ≤ 6%), and enclosure aspect ratios (½ ≤ A ≤ 2). Different behaviors (enhancement or deterioration) are predicted in the average Nusselt number as the volume fraction of nanoparticles increases depending on the combination of CuO–EG–Water variable thermal conductivity and viscosity models employed. In general, the effects the viscosity models are predicted to be more predominant on the behavior of the average Nusselt number than the influence of the thermal conductivity models. The enclosure aspect ratio is predicted to have significant effects on the behavior of the average Nusselt number which decreases as the enclosure aspect ratio increases.     

Steady natural convection flow in a square cavity filled with a porous medium for linearly heated side wall(s)

In this paper natural convection flows in a square cavity filled with a porous matrix has been investigated numerically when the bottom wall is uniformly heated and vertical wall(s) are linearly heated whereas the top wall is well insulated. Darcy–Forchheimer model without the inertia term is used to simulate the momentum transfer in the porous medium. Penalty finite element method with bi-quadratic rectangular elements is used to solve the non-dimensional governing equations. Numerical results are presented for a range of parameters (Rayleigh number Ra, 10³ ⩽ Ra ⩽ 10⁶, Darcy number Da, 10⁻⁵ ⩽ Da ⩽ 10⁻³, and Prandtl number Pr, 0.2 ⩽ Pr ⩽ 100) in terms of stream functions and isotherm contours, and local and average Nusselt numbers.     

Numerical research of nature convective heat transfer enhancement filled with nanofluids in rectangular enclosures

Heat transfer enhancement utilizing nanofluids in a two-dimensional enclosure is investigated for various pertinent parameters. The Khanafer's model is used to analyze heat transfer performance of nanofluids inside an enclosure taking into account the solid particle dispersion. Transport equations are model by a stream function-vorticity formulation and are solved numerically by finite-difference approach. Based upon the numerical predictions, the effects of Rayleigh number (Ra) and aspect ratio (AR) on the flow pattern and energy transport within the thermal boundary layer are presented. The diameter of the nanoparticle d_p is taken as 10 nm in nanofluids. The buoyancy parameter is 10³ ≤ Ra ≤ 10⁶ and aspect ratios (AR) of two-dimensional enclosure are 1/2, 1, 2. Results show that increasing the buoyancy parameter and volume fraction of nanofluids cause an increase in the average heat transfer coefficient. Finally, the empirical equation was built between average Nusselt number and volume fraction.     

Effects of jet plate size and plate spacing on the stagnation Nusselt number for a confined circular air jet impinging on a flat surface

This work deals with the effects of jet plate size and plate spacing (jet height) on the heat transfer characteristics for a confined circular air jet vertically impinging on a flat plate. The jet after impingement was restricted to flow in two opposite directions. A constant surface heat flux of 1000 W/m² was arranged. Totally 88 experiments were performed. Jet orifices individually with diameter of 1.5, 3, 6 and 9 mm were adopted. Jet Reynolds number (Re) was in the range 10,000–30,000 and plate spacing-to-jet diameter ratio (H/d) was in the range 1–6. Eleven jet plate width-to-jet diameter ratios (W/d = 4.17–41.7) and seven jet plate length-to-jet diameter ratios (L/d = 5.5–166.7) were individually considered. The measured data were correlated into a simple equation. It was found that the stagnation Nusselt number is proportional to the 0.638 power of the Re and inversely proportional to the 0.3 power of the H/d. The stagnation Nusselt number was also found to be a function of exp[−0.044(W/d) − 0.011(L/d)]. Through comparisons among the present obtained data and documented results, it may infer that, for a jet impingement, the impingement-plate heating condition and flow arrangement of the jet after impingement are two important factors affecting the dependence of the stagnation Nusselt number on H/d.