Forced convection heat transfer in tube banks in cross flow

The forced convection heat transfer characteristics for an incompressible, steady and Newtonian fluid flow over a bundle of circular cylinders has been investigated numerically. The inter-cylinder hydrodynamic interactions have been approximated by employing a simple cell model. The momentum and energy equations have been solved by using a finite difference based numerical solution procedure for a range of physical and kinematic conditions. Furthermore, the role of the type of thermal boundary condition, namely, a constant temperature or a constant heat flux, imposed on the surface of the cylinder has also been elucidated. Extensive results on the temperature fields, and on the variation of the Nusselt number on the surface of a typical cylinder in the assemblage have been obtained for two values of the Prandtl number (corresponding to air and water). The Reynolds number of flow was varied in the range 1-500 and the voidage of the assemblage ranged from 0.4 to 0.99 thereby covering the entire range of interest as encountered in tubular heat exchangers and in fibrous beds. The paper is concluded by presenting extensive comparisons with the limited analytical/numerical and/or experimental results available in the literature for the case of a single cylinder as well as that for tube bundles.     

Convective heat transfer for power law fluids in packed and fluidised beds of spheres

The forced convection heat transfer characteristics for an incompressible and steady flow of power law liquids in fixed and extended beds of spherical particles has been studied numerically. The sphere-sphere hydrodynamic interactions have been accounted for by using a simple cell model. Within the framework of such a cell model, the momentum and energy equations have been solved using a finite difference method to obtain the velocity and temperature fields. Extensive numerical estimates of the local and average Nusselt numbers as functions of the physical, rheological and kinematic variables have been presented and discussed for the two commonly employed thermal boundary conditions. In broad terms, the Nusselt number for power law fluids (both shear-thinning and shear-thickening conditions) normalized with respect to the corresponding value for a Newtonian fluid shows weak additional dependence on the power law flow behaviour index. The shear-thinning behaviour is seen to promote heat transfer and as expected the shear-thickening behaviour impedes heat transfer in fixed and fluidised beds. All in all, the present results encompass wide ranges of conditions as follows: Reynolds number: 1-500; Peclet number: 1-500; bed voidage: 0.4-0.8 and the flow behaviour index: 0.5-1.8 thereby covering extremely shear-thinning and shear-thickening types of fluid behaviours. The paper is concluded by presenting detailed comparisons with the limited analytical and/or experimental results available for liquid-solid mass transfer in such systems.     

Laminar natural convection characteristics in an enclosure with heated hexagonal block for non-Newtonian power law fluids

This work illustrates the steady state, two dimensional natural convective flow and heat transfer features in square enclosure containing heated hexagonal block maintained either at constant wall temperature(CWT) or uniform heat flux(UHF) thermal conditions. Governing equations(mass, momentum and energy) are solved by using finite volume method(FVM) with 3rd order accurate QUICK discretization scheme and SIMPLE algorithm for range of field pertinent parameters such as, Grashof number(10~3≤ Gr ≤ 10~6), Prandtl number(1 ≤ Pr ≤ 100) and power law index(0.5 ≤ n ≤ 1.5). The analysis of momentum and heat transfer characteristics are delineated by evolution of streamlines, isotherms, variation of average Nusselt number value and Colburn factor for natural convection(j_(nH)). A remarkable change is observed on fluid flow and thermal distribution pattern in cavity for both thermal conditions. Nusselt number shows linear variation with Grashof and Prandtl numbers; while rate of heat transfer by convection decreases for power law index value. Higher heat transfer rate can be achieved by using uniform heat flux condition. A Nusselt number correlation is developed for possible utilization in engineering/scientific design purpose.     

Using CFD to predict the behavior of power law fluids near axial-flow impellers operating in the transitional flow regime

A computational fluid dynamics (CFD) model of flow in a mixing tank with a single axial-flow impeller was developed with the Fluent^TM software. The model consists of an unstructured hexagonal mesh (158,000 total cells), dense in the region from the surface of the impeller. The flow was modeled as laminar and a multiple reference frame approach was used to solve the discretized equations of motion in one-quarter of a baffled tank. A solution of 0.1% Carbopol in water, a shear-thinning fluid, was found to be clear enough to measure impeller discharge angles using laser Doppler velocimetry. This is the first time that impeller discharge angles have been reported in the literature for a shear-thinning fluid with a hydrofoil impeller. Rheological measurements indicated that the Carbopol solution can be characterized by the power law (K=9,n=0.2) under the range of shear conditions (0.1- expected near the impeller in the mixing tank. The CFD model accurately predicted the dependence of power number and discharge angle on Reynolds number (as predicted by Metzner and Otto), for an A200 (pitched blade turbine or PBT) and an A315 (Hydrofoil) impeller operating in the transitional flow regime (Reynolds numbers: 25-400) with glycerin and 0.1% Carbopol solutions. Subsequently, the results of a systematic CFD study with power law fluids indicated that the power number and discharge angle of an axial-flow impeller in the transitional flow regime depends not only on the Reynolds number (as determined by Metzner and Otto's method) but also on the flow behavior index n. Consequently, an alternative to Metzner and Otto's method was pursued. The results of converged CFD simulations indicate that the near-impeller “average shear rate” increases not only with increasing RPM (as proposed by Metzner and Otto), but also with decreasing flow behavior index (n) and discharge angle in the transitional flow regime. Considering this result, an improved method of estimating the power number and discharge angle for power law fluids in the transitional flow regime is proposed.     

Flow of Newtonian and power law liquids in tube bundles

In this work, the annular (tangential) flow of Newtonian and non‐Newtonian fluids in tube bundles has been studied experimentally. Extensive pressure drop data has been obtained embracing wide ranges of the Reynolds number (13–6600) and for two test modules of different geometrical arrangements, but of similar overall void fraction. Preliminary experiments suggest that the pressure drop is mainly determined by the overall void fraction of the bundle and is relatively insensitive to the detailed geometrical configuration of the bundle. A simple predictive correlation has been developed which reconciles the present results for Newtonian and power law fluids with acceptable levels of reliability.     

Effect of blockage on heat transfer from a cylinder to power law liquids

The influence of planar confining walls on the steady forced convection heat transfer from a cylinder to power-law fluids has been investigated numerically by solving the field equations using FLUENT (version 6.2). Extensive results highlighting the effects of the Reynolds number (1?Re?40), power-law index (0.2?n?1.8), Prandtl number (1?Pr?100) and the blockage ratio (β=4 and 1.6) on the average Nusselt number have been presented. For a fixed value of the blockage ratio, the heat transfer is enhanced with the increasing degree of shear-thinning behaviour of the fluid, while an opposite trend was observed in shear-thickening fluids. Due to the modifications of the flow and temperature fields close to the cylinder, the closely placed walls (i.e., decreasing value of the blockage ratio) further enhance the rate of heat transfer as the fluid behaviour changes from Newtonian to shear-thickening fluids (n>1), the opposite influence is seen with the decreasing value of the flow behaviour index (n) in shear-thinning (n<1) fluids. Finally, the functional dependence of the present numerical results on the relevant dimensionless parameters has been presented in the form of closure relationships for their easy use in a new application.     

Two-dimensional unsteady flow of power-law fluids over a cylinder

The unsteady flow of incompressible power-law fluids over an unconfined circular cylinder in cross-flow arrangement has been studied numerically. The two-dimensional (2-D) field equations have been solved using a finite volume method based solver (FLUENT 6.3). In particular, the effects of the power-law index (0.4?n?1.8) and Reynolds number (40?Re?140) on the detailed kinematics of the flow (streamline, surface pressure and vorticity patterns) and on the macroscopic parameters (drag and lift coefficients, Strouhal number) are presented in detail. The periodic vortex shedding and the evolution of detailed kinematics with time are also presented to provide insights into the nature of flow. The two-dimensional flow transits from steady to unsteady behaviour at a critical value of the Reynolds number Re∼(40-50) and the von-Karman vortex street is observed beyond the critical Reynolds number (Re). Obviously, both the lift coefficient and Strouhal number values are zero for the steady flow, but their values increase with the increasing Reynolds number (Re) in the unsteady flow regime. For highly shear-thickening fluids (n=1.8), the flow becomes unsteady at Re=40 and unsteadiness in the flow appears at Re=50 for all values of power-law index (n). As expected, the evolution of the kinematics and vortex shedding show a complex dependence on the flow parameters near the transition in the flow. For a fixed value of the Reynolds number (Re), the drag coefficient increases and lift coefficient decreases with increasing value of the power-law index (n). For a fixed value of the power-law index (n), the drag coefficient gradually increases with the Reynolds number (Re). Similar to the drag coefficient, lift coefficient also shows a complex dependence on the power-law index (n) near the transition zone. The value of the Strouhal number (St) decreases with the increasing value of the power-law index (n) at a fixed value of the Reynolds number (Re).     

Steady flow of power-law fluids across an unconfined elliptical cylinder

The momentum transfer characteristics of the power-law fluid flow past an unconfined elliptic cylinder is investigated numerically by solving continuity and momentum equations using FLUENT (version 6.2) in the two-dimensional steady cross-flow regime. The influence of the power-law index (0.2?n?1.8), Reynolds number (0.01?Re?40) and the aspect ratio of the elliptic cylinder (0.2?E?5) on the local and global flow characteristics has been studied. In addition, flow patterns showing streamline and vorticity profiles, and the pressure distribution on the surface of the cylinder have also been presented to provide further physical insights into the detailed flow kinematics. For shear-thinning (n<1) behaviour and the aspect ratio E>1, flow separation is somewhat delayed and the resulting wake is also shorter; on the other hand, for shear-thickening (n>1) fluid behaviour and for E<1, the opposite behaviour is obtained. The pressure coefficient and drag coefficient show a complex dependence on the Reynolds number and power-law index. The decrease in the degree of shear-thinning behaviour increases the drag coefficient, especially at low Reynolds numbers. While the aspect ratio of the cylinder exerts significant influence on the detailed flow characteristics, the total drag coefficient is only weakly dependent on the aspect ratio in shear-thickening fluids. The effect of the flow behaviour index, however, diminishes gradually with the increasing Reynolds number. The numerical results have also been presented in terms of closure relations for easy use in a new application.     

Rising velocity of a swarm of spherical bubbles in power law fluids at high reynolds numbers

In this work, a theoretical scheme for estimating the rise velocity of a swarm of spherical bubbles through quiescent power law liquids at high Reynolds number is developed. The inter-bubble interactions have been accounted for by the use of a cell model. The effect of the power law index and the volume fraction of the gas on the rise velocity is elucidated. Depending upon the degree of shearthinning behaviour and the gas fraction, the swarm may rise slower or faster than a single bubble. This behaviour has been explained qualitatively in terms of two competing mechanisms.     

Hindered settling in non-newtonian power law liquids

New experimental results on the hindered settling of model glass bead suspensions in non-Newtonian suspending media are reported. The data presented encompass the following ranges of variables: 7.38 × 10^?4 ≤ Re_1∞ ≤ 2; 0.0083 ≤ d/D ≤ 0.0703; 0.13 ≤ C ≤ 0.43 and 1 ≥ n ≥ 0.8. In these ranges of conditions, the dependence of the hindered settling velocity on concentration is adequately represented by the corresponding Newtonian expressions available in the literature. The influence of the power law flow behaviour index is completely embodied in the modified definition of the Reynolds number used for power law liquids.     

Settling of cylinders in power law liquids

New measurements on drag coefficient and wall effects on the free settling motion of cylinders in two shearthinning polymer solutions are reported. The ranges of conditions covered in this study are: 1 < Re_∞∞ < 40; 0.25 < (LID) < 2; 0.079 < dID < 0.4, and n = 0.65 and 0.74. Both wall correction factor and drag coefficient results are in line with the Newtonian behavior provided a modified Reynolds number is used to represent the results.     

Rising velocity of a swarm of spherical bubbles in a power law non-newtonian liquid

Previous work on the rise of a swarm of bubbles in non-Newtonian media has been reviewed. Variational principles have been combined with the Happel free surface cell model to obtain upper and lower bounds on the swarm velocities of bubbles rising slowly through power law liquids. The predictions presented herein encompass wide ranges of gas holdup and shear-thinning behaviour.     

A simple correlation for gas bubbles rising in power‐law fluids

Recently, the author (Rodrigue, 2001 a, b) proposed a generalized correlation for the steady rise of gas bubbles in uncontaminated viscous Newtonian fluids of infinite extent. It is the purpose of this note to show that this model can be modified for inelastic non‐Newtonian power‐law fluids. Using data taken from six different studies, it is shown that the modified model can predict quite nicely the bubble velocity for these non‐Newtonian inelastic fluids.     

Fluid-dynamic experimental study in a slurry bubble column with a tube bank in cross flow

In this work, the fluid-dynamic of a slurry bubble column provided with a tube bank operating in cross flow was studied. The main objective was to determine the influence of the air velocity, the slurry velocity, and the solid concentration in the gas holdup and solid distribution inside a bed with internals. Pilot scale equipment with a column of square traverse section (0.144 m 2 ) of stainless steel with 2.20 m height was used. The column base has a truncated pyramidal shape with an inclination angle to the vertical of 20°. In the central part of the equipment, there is a bank of 49 tubes ( D =2.54 cm and L =38 cm) in an aligned arrangement ( S_L=S_T=4.4,rm cm ). The three-phase system used was compressed air, tap water, and sieved sand ( D_p=0.0505,rm cm ). The velocities of liquid and gas phases were varied in the range from 0.31 to 1.24 cm/s and from 1.41 to 3.15 cm/s, respectively, and the solids load concentration from 0 to 30% w/w. According to the results, when the liquid and gas phase velocities increase, the gas holdup, varepsilon_G , increases, while when the solids load concentration increases, varepsilon_G presents a local maximum at intermediate values, varepsilon_G being the largest for the two-phase system. Analogously, the axial profiles of solids concentration showed a local maximum in the central region of the column, exactly in the area of the tube bank. In the traverse direction, the solids concentration increases gradually heading to the center of the column. These profiles, axial and traverse, become more uniform as the air velocity increases and the slurry as well. In addition, according to expectations, the solids concentration increases in each point of the column as the load concentration is increased in the system.     

Unified entry length for newtonian and power–law fluids in laminar pipe flow

A numerical method based on finite differencing is used for investigating the steady–state entrance region laminar flow of incompressible Newtonian and power–law fluids in a circular pipe. The Solution method is validated by comparing the results for Newtonian fluids with those reported in the literature. For power–law fluids, the entry length results are compared with other approximate solutions in the literature. On the basis of the calculated results, a generalized entry length ξ99 = 0.056 is shown to be valid for the laminar flow at Re > 200 of both Newtonian and power–law fluids with 0.75 < n < 1.5.     

Effect of power-law index on critical parameters for power-law flow across an unconfined circular cylinder

The conditions for the formation of a wake and for the onset of wake instability for the flow of power-law fluids over an unconfined circular cylinder are investigated numerically by solving the continuity and momentum equations using FLUENT (version 6.2). The effect of power-law index on the critical Reynolds numbers, Strouhal number and drag coefficient has been presented over a wide range of power-law index (0.3?n?1.8) thereby establishing the limits of the flow without separation and the steady symmetric flow regimes, respectively. While both the shear-thinning (n<1) and the shear-thickening (n>1) seem to lower the value of the critical Reynolds number denoting the onset of wake instability as compared to that for Newtonian fluids, the effect is seen to be more prominent for shear-thickening fluids than that for shear-thinning fluids. The corresponding values of the critical Strouhal number (St_c) and drag coefficient have also been presented for the critical values of the Reynolds number. Included here are also a series of streamline plots showing the onset of asymmetry and of the time-dependent flow regime.     

Simulation of forced convection in a channel with nanofluid by the lattice Boltzmann method

This paper presents a numerical study of the thermal performance of fins mounted on the bottom wall of a horizontal channel and cooled with either pure water or an Al₂O₃-water nanofluid. The bottom wall of the channel is heated at a constant temperature and cooled by mixed convection of laminar flow at a relatively low temperature. The results of the numerical simulation indicate that the heat transfer rate of fins is significantly affected by the Reynolds number (Re) and the thermal conductivity of the fins. The influence of the solid volume fraction on the increase of heat transfer is more noticeable at higher values of the Re.     

A note on pressure drop for the cross-flow of power-law liquids and air/power law liquid mixtures past a bundle of circular rods

The steady flow of non-Newtonian polymer solutions on their own and together with air transverse to a bundle of circular rods, has been studied experimentally. In particular, the frictional pressure drop has been measured as a function of non-Newtonian power-law constants, flow rates, voidage of bundles and the input fraction of air over wide ranges of conditions as: power-law flow behaviour index, 0.5 to 1; voidage values of 0.78 and 0.87-0.88; air velocity (superficial), 0.09 to 0.28 m/s and liquid velocity (superficial) 0.018 to 0.5 m/s. While the single phase flow results for power-law liquids are in excellent agreement with an existing equation, the two-phase pressure drop results are also in accordance with the form of the well-known Lockhart-Martinelli correlation. Based on these results, it is suggested that a value of 5-10 for the Reynolds number based on rod diameter marks the limit of the laminar flow in such systems.