Effects of viscous dissipation and fluid axial heat conduction on heat transfer for non-Newtonian fluids in ducts with uniform wall temperature: Part I: Parallel plates and circular ducts

Laminar heat transfer in parallel plates and circular ducts subject to uniform wall temperature is studied by taking into account both viscous dissipation and fluid axial heat conduction in an infinite region. Developing temperature fields are evaluated numerically by a finite-difference method for various Brinkman numbers (Br) and Peclet numbers (Pe). Nusselt numbers are presented graphically for Pe = 10 and Pe → ∞, and Br = 0, ± 0.5 and ± 1 for non-Newtonian fluids described by the power-law model with the flow index of n = 0.5, 1.0 and 1.5. It is shown that Nusselt number has a single fixed value independent of Br in the thermally developing region and its numerical value is equal to that at the fully developed region for non-zero Br, when the preheating of incoming fluid due to both viscous dissipation and fluid axial heat conduction is considered.     

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.     

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.     

A numerical study of the unsteady heat/mass transfer inside a circulating sphere

Numerical methods are used to investigate the transient, forced convection heat/mass transfer inside circulating spheres at low to moderate Reynolds numbers. The heat/mass balance equations were solved numerically in spherical coordinates system by a finite difference method. The values considered for the sphere interior Reynolds number are Re_int ? 1000. The computations were focused on the influence of the sphere Peclet number, Pe, and Re_int on heat/mass transfer rate for Pe/(1 + μ) ? 10⁴.     

Flow and heat transfer in a micro-cylindrical gas–liquid Couette flow

Analytic solutions for the gas and liquid velocity and temperature distribution are determined for steady state one-dimensional microchannel cylindrical Couette flow between a shaft and a concentric cylinder. The solution is based on the continuum model and takes into consideration the velocity slip and temperature jump in the gaseous phase defined by the Knudsen number range of 0.001 < Kn < 0.1. The two fluids are assumed immiscible. The gas layer is adjacent to the shaft which rotates with angular velocity ω_s and is thermally insulated. The outer cylinder rotates with angular velocity ω_o and is maintained at uniform temperature. The governing parameters are identified and the effects of the Knudsen number and accommodation coefficients on the velocity and temperature profiles, reduction in the overall temperature rise due to the gas layer, the Nusselt number and shear reduction are examined. It was found that the required torque to rotate the liquid in the annular space is significantly reduced by introducing a thin gas layer adjacent to the shaft. Also, reduction in shaft temperature is enhanced through a combination of high energy accommodation coefficient and low momentum accommodation coefficients. Results also indicate that the gas layer becomes more effective in reducing the shaft temperature when the housing angular velocity is much larger than the shaft angular velocity.     

Heat transfer behavior in a rotating aluminum foam heat sink with a circular impinging jet

This work experimentally studied heat transfer associated with an impinging jet onto a rotating heat sink. Air was used as the impinging coolant, and a square Al-foam heat sink was adopted. The variable parameters were the jet Reynolds number (Re), the relative nozzle-to-foam tip distance (C/d), the rotational Reynolds number (Re_r) and the relative side length of the square heat sink (L/d). The effects of Re, C/d, Re_r and L/d on the dimensionless temperature distributions and the average Nusselt number were considered. For a stationary system, the results reveal that the average Nusselt number (Nu₀) with Al-foam was two to three times that without Al-foam. Nu₀ increased with Re. A larger L/d responded to a larger Nu₀ based on the same jet flow rate. The effect of C/d on Nu₀ was negligible herein. For a rotating system, when Re and L/d were small and C/d was large, the average Nusselt number (Nu_Ω) increased considerably with Re_r. Additionally, for Nu_Ω/Nu₀ ? 1.1, the results suggest that rotation was substantial at Re_r/Re ? 1.13 when L/d = 4.615 with C/d = 0–5 and at Re_r/Re ? 1.07 when L/d = 3.0 with C/d = 0–5. For L/d = 2.222, rotation was substantial at Re_r/Re ? 1.44 when C/d = 0 and was always substantial when C/d ? 1.     

Effects of viscous dissipation during forced convection of power-law fluids in microchannels

The effect of viscous dissipation in forced convection of power-law fluids through microchannels of different cross-sectional geometries is studied numerically over the ranges of power-law index, 0.8 ≤ n ≤ 1.2 and Brinkman number, 0.001 ≤ Br ≤ 0.1 while keeping Péclet number constant at Pe = 10. Two types of thermal boundary conditions, namely, uniform wall temperature (T₂) and uniform heat flux (H₂), have been employed at the microchannel wall and the results of the temperature fields are expressed in terms of Nusselt number. The interplay between the fluid rheology and viscous dissipation effect gives rise to significant alteration in the net convective transport and thus can be beneficial in the thermal design of biofluidic devices.     

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.     

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.     

Unsteady conjugate forced convection heat/mass transfer in ensembles of spherical particles with cell models

Numerical methods are used to investigate the transient, conjugate, forced convection heat/mass transfer in multiparticle systems at low to moderate Reynolds numbers. The interparticle interactions have been accounted for by using the simple cell models. The momentum and heat/mass balance equations were solved numerically in spherical coordinates system by a finite difference method. The values considered for the sphere Reynolds number are Re < 100. The computations were focused on the influence of the voidage and physical properties ratios on the heat/mass transfer rate for sphere Peclet number, 10 ? Pe ? 1000.     

Jet impingement interaction with cross flow in horizontal porous layer under thermal non-equilibrium conditions

In the present article the jet impingement cooling of heated portion of a horizontal surface immersed in a thermally non-equilibrium porous layer is considered for investigation numerically with the presence of a cross flow. The mathematical model is derived for steady, two-dimensional laminar flow based on Darcy model and two-energy equation for fluid and solid phases. A parametric study is carried out by varying the following parameters: cross flow to jet flow velocity ratio parameter (0 ⩽ M ⩽ 1); porosity scaled thermal conductivity ratio parameter (0.1 ⩽ K_r ⩽ 1000); heat transfer coefficient parameter (0.1 ⩽ H ⩽ 1000); Péclet number (1 ⩽ Pe ⩽ 1000) and Rayleigh number (10 ⩽ Ra ⩽ 100). The total average Nusselt number is defined based on the overall thermal conductivity, which is assumed to be the arithmetic mean of the porosity scaled thermal conductivity of the fluid and solid phases. The total average Nusselt number as well as the average Nusselt number for both fluid and solid phases is presented for different governing parameters. It is found that the presence of a weak cross flow in a jet impinging jet may degrade the heat transfer. The results show that the average Nusselt number calculated from the thermal equilibrium model are the maximum possible values and these values can be reproduced by large values of H × K_r.     

Asymptotic analysis of heat transfer in turbulent Rayleigh–Bénard convection

Based on asymptotic considerations a heat transfer law for turbulent Rayleigh–Bénard convection is found that differs from existing correlations which often are of a power law type with respect to their Rayleigh number dependence. From the asymptotic temperature profile, derived by matching temperature gradients in the overlap region of the wall layer and the core layer, a Nusselt number follows which includes a logarithmic term. This correlation is in good agreement with data from highly accurate Rayleigh–Bénard experiments for Rayleigh numbers between 10⁵ and 10¹⁵ and Prandtl numbers larger than 0.5. It is an alternative to existing power laws or more complicated correlations for Nu = Nu(Ra,Pr).     

Effects of viscous dissipation on slip-flow heat transfer in a micro annulus

This study deals numerically with the laminar slip-flow forced convection in a micro annulus with constant wall temperature. The solution takes the effects of viscous dissipation, velocity-slip and temperature-jump conditions at the surface into considerations. A hybrid application of the Laplace transformation technique and the local adaptive differential quadrature method (La-DQM) is used to solve the energy equation for the developing temperature field. Of interest are the effects of the Brinkman number Br, the Knudsen number Kn, and the radius ratio β of inner to outer cylinders on the temperature distribution and the Nusselt number.The results indicate that fully developed Nusselt number increases with an increase in the radius ratio β or a decrease in the Knudsen number Kn, but it is rather insensitive to the Brinkman number. The local Nusselt number in the thermal entrance region increases as Br increases. Comparisons are also made of the dimensionless temperature gradient at the outer cylinder between the present results and published data for the circular tube case β = 0, and the agreement is found out to be generally good.     

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.     

Influence of channel aspect ratio on heat transfer in rotating rectangular ducts with skewed ribs at high rotation numbers

Centerline heat transfer measurements along two opposite ribbed walls in three rotating rectangular ducts roughened by 45° staggered ribs with channel aspect ratios (AR) of 1:1, 2:1 and 4:1 are performed at Reynolds (Re), rotation (Ro) and buoyancy (Bu) numbers in the ranges of 5000–30,000, 0–2, and 0.005–8.879, respectively. These channel geometries are in common use as the internal cooling passages of a gas turbine rotor blade and the tested Ro and Bu ranges are considerably extended from the previous experiences. This study focuses on the heat transfer characteristics in response to the change of AR under the parameter ranges examined. With zero-rotation (Ro = 0), the local Nusselt numbers (Nu₀) along the centerlines of two opposite ribbed walls increase as AR increases due to the increased rib-height to channel-height ratio. The Bu impact on heat transfer appears to be AR dependent, i.e. the increase of Bu elevates Nusselt number ratios Nu/Nu₀ in the square channel but impairs heat transfer in the rectangular channels of AR = 2 and 4. Acting by the Coriolis effect alone, all the leading edge Nu values in the present Ro range are lower than the zero-rotation references but started to recover as Ro increases from 0.1 in the channels of AR = 1, 2 and from 0.3 in the channel of AR = 4. The trailing edge Nu/Nu₀ ratios increase consistently from unity as Ro increases but their responses toward the increase of AR are less systematic than those found along the leading edge. The above findings, with the aids of extended Ro and Bu ranges achieved by this study, serve as the original contributions for this technical community. The Nu/Nu₀ ratios in the rotating channels of AR = 1, 2, and 4 fall in the ranges of 0.6–2.2, 0.5–2.7, and 0.5–2.1, respectively. A set of heat transfer correlations is derived to represent all the heat transfer data in the periodically developed flow regions of three rotating ducts.     

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.     

Flow visualizations and heat transfer measurements for a rotating pin-fin heat sink with a circular impinging jet

This work experimentally investigated the fluid flow and heat transfer behaviors of jet impingement onto the rotating heat sink. Air was used as impinging coolant, while the square heat sinks with uniformly in-line arranged 5 × 5 and 9 × 9 pin-fins were employed. The side length (L) of the heat sink equaled 60 mm and was fixed. Variable parameters were the relative length of the heat sink (L/d = 2.222 and 4.615), the relative distance of nozzle-to-fin tip (C/d = 0–11), the jet Reynolds number (Re = 5019–25,096) and the rotational Reynolds number (Re_r = 0–8114). Both flow characteristics of stationary and rotating systems were illustrated by the smoke visualization. Besides, the results of heat transfer indicate that, for a stationary system with a given air flow rate, there was a larger average Nusselt number (Nu₀) for the 9 × 9 pin-fin heat sink with L/d = 4.615 and C/d = 11. For a rotating system, a bigger Re_r meant a more obvious heat transfer enhancement (Nu_Ω/Nu₀) in the case of smaller Re, but Nu_Ω/Nu₀ decreased with increasing Re. In this work, Nu_Ω/Nu₀ in L/d = 2.222 is higher than in L/d = 4.615; among the systems in L/d = 2.222, bigger Nu_Ω/Nu₀ exists in the case of C/d = 9–11, but among the systems in L/d = 4.615, bigger Nu_Ω/Nu₀ exists in the case of C/d = 1–3. Finally, according to the base of Nu_Ω/Nu₀ ? 1.1, the criterion of the substantial rotation was suggested to be Re_r/Re ? 1.154.     

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.     

A detailed comparison between Navier–Stokes and DSMC simulations of multicomponent gaseous flow in microchannels

Numerical results of the Navier–Stokes equations and the DSMC method for H₂/N₂ and H₂/N₂/CO₂ mixtures in planar microchannels are compared in order to verify the accuracy of slip/jump boundary conditions in multicomponent flows. For the solution of the Navier–Stokes and energy equations, a finite-volume method is used together with a complete set of slip/jump boundary conditions derived from the kinetic theory of gases. In order to compare the accuracy of these two methods, different wall temperatures as well as inlet mass flow rates are examined. Isothermal flow is also considered for comparison with available analytical solutions for slip flows. The Knudsen number typically ranged between 0.015 < Kn < 0.09 with a Reynolds number of 5 < Re < 15. The two methods are generally in good agreement under the conditions studied especially at low to moderate Knudsen numbers.     

Experimental investigation of natural convection from vertical and horizontal helicoidal pipes in HVAC applications

Natural convection from uniformly heated helicoidal pipes oriented vertically and horizontally is experimentally studied. Four helicoidal pipes of different parameters are presented. The effects of pitch to pipe diameter ratio, coil diameter to pipe diameter ratio and length to pipe diameter ratio on the average heat transfer coefficient are found. The experiments covered a range of Rayleigh number based on tube diameter from 1.5 × 10³ to 1.1 × 10⁵. The results showed that the overall average Nusselt number, Nu_m, increases with the increase in pitch to pipe diameter ratio, coil diameter to pipe diameter ratio and length to pipe diameter ratio for vertical helicoidal pipes. For horizontal helicoidal pipes, the overall average Nusselt number, Nu_m, increases with the increase in pitch to pipe diameter ratio and length to pipe diameter ratio, but it decreases with the increase in coil diameter to pipe diameter ratio. New correlations are presented, and they can be used in HVAC applications.