Free Convection in Confined Power-Law Fluids From Two Side-by-Side Cylinders in a Square Enclosure

Laminar free convection heat transfer in power-law fluids from two side-by-side cylinders (one hot and one cold) confined in a square duct has been studied numerically in the two-dimensional flow regime. For a fixed value of the ratio of cylinder radius to size of the enclosure, the effect of geometrical placement of the cylinders is studied on the resulting velocity and temperature fields in the laminar free convection regime by considering six asymmetric locations of the two cylinders. In particular, extensive results reported herein span the range of conditions of Grashof number, 10 to 10⁵; Prandtl number, 0.7 to 100, thereby yielding the range of the Rayleigh number as 7 to 10⁷; power-law index, 0.3 to 1.8; and the relative positions (dimensionless) of the cylinders with respect to the centerline, –0.25 to 0.25. The heat transfer characteristics are analyzed in terms of the local Nusselt number along the surfaces of the two cylinders and the enclosure walls. Overall, the average Nusselt number shows a positive dependence on both the Grashof number and the Prandtl number irrespective of the values of power-law index and relative positioning of the cylinders. Also, all else being equal, shear-thinning fluid behavior promotes heat transfer with reference to that in Newtonian fluids. When the two cylinders are situated close to the bottom wall, the rate of heat transfer is augmented with reference to that for the symmetric positioning of the cylinders along the horizontal mid-plane of the enclosure. Conversely, heat transfer deteriorates as the cylinders are located above the centerline of the enclosure. The present numerical results have been consolidated via the use of a modified Rayleigh number, thereby enabling the estimation of the average Nusselt number in a new application.     

Effect of Prandtl Number on Free Convection from Two Cylinders in a Square Enclosure

This study explores the effect of Prandtl number on the laminar natural convection heat transfer to Newtonian fluids in a square enclosure consisting of one hot circular cylinder and one cold circular cylinder. The walls of the square enclosure are maintained isothermal and at the same temperature as the cold cylinder and the fluid medium. The governing partial differential equations have been solved numerically over the following ranges of conditions: Grashof number, 10 to 10⁵; Prandtl number, 0.7 to 100 (or the range of Rayleigh numbers as 7 to 10⁷); and relative positioning of the cylinders, ?0.25 to 0.25. However, the ratio of the radius of the cylinder to the side of the enclosure is held fixed at 0.2. Extensive results on the streamline and isotherm contours, the local Nusselt number distribution, and the average Nusselt number are discussed to delineate the influence of Grashof and Prandtl numbers on them for a given location with respect to the horizontal center line. The surface-averaged Nusselt number shows a positive dependence on Grashof and Prandtl numbers for a fixed location of the two cylinders. The heat transfer results have been correlated as a function of the Rayleigh number and geometric parameters, thereby enabling its prediction in a new application.     

Laminar Free Convection in Power-law Fluids in a Right Angle Triangular Duct with Heated Base

ABSTRACT

Laminar free convection in power-law fluids in a triangular duct is studied numerically to delineate the effects of the height-to-base ratio of the enclosure (0.2 to 2), power-law index (0.2 to 1.8), Grashof number (10 to 10⁴) and Prandtl number (0.7 to 100). The heat transfer is analyzed for the heated base with the other two walls being cold. Detailed kinematics is characterized by the formation of multiple recirculating zones ranging from two to four cells. Shear rate contours provide additional insights about the variation of the local viscosity in the fluid. Heatlines and the values of the Bejan number over the range of conditions are calculated to delineate the contributions of the entropy generation due to thermal effects and viscous dissipation. At low Grashof and/or Prandtl numbers, conduction dominates the overall heat transfer and this transition between the conduction and convection-dominated regimes is captured in terms of a modified Rayleigh number. The effect of aspect ratio on the Nusselt number is modulated by the values of Grashof and Prandtl numbers and power-law index. The present results have been consolidated via the use of a modified Rayleigh number for estimating the value of average Nusselt number in a new application.     

Effect of Fluid Yield Stress on Natural Convection From Horizontal Cylinders in a Square Enclosure

In this study, laminar natural convection heat transfer to Bingham plastic fluids from two differentially heated isothermal cylinders confined in a square enclosure (with isothermal walls) has been investigated numerically. The governing partial differential equations have been solved over the ranges of the dimensionless parameters, namely, Rayleigh number, 10² to 10⁶, Prandtl number, 10 to 100, and Bingham number, 0.01 to 100, for seven locations of inner cylinders as ±0.25, ±0.2, ±0.1 and 0. These values correspond to the range of Grashof number varying from 10 to 10⁵. The detailed flow and temperature fields are visualized in terms of the streamlines and isotherm contours. Further insights are developed by examining the iso-shear rate contours and the yield surfaces delineating the fluid-like and solid-like regions. The corresponding heat transfer results are analyzed in terms of the distribution of the local Nusselt number along the cylinder surface together with its surface averaged value as functions of the Rayleigh number, Prandtl number, Bingham number, and positions of the cylinders. It is found that the average Nusselt number increases with the increasing values of the Rayleigh number and decreases with the increasing Bingham number. For sufficiently large values of the Bingham number, the average Nusselt number reaches its asymptotic value wherein heat transfer takes place solely by conduction. Based on the present numerical results, simple correlations for the prediction of the average Nusselt number and the limiting Bingham number have been developed. Also, a dimensionless criterion denoting the cessation of convection regime is outlined for this configuration.     

Laminar free convection from a horizontal semi-circular cylinder to power-law fluids

Extensive numerical results on the flow and thermal fields are presented for free convection from a semi-circular cylinder (flat base upward) immersed in quiescent power-law fluids for the following ranges of conditions: Grashof number, 10 ? Gr ? 10⁵, Prandtl number, 0.72 ? Pr ? 100, and power-law index, 0.2 ? n ? 1.8. The heat transfer characteristics are analyzed in terms of the isotherm patterns, local and average Nusselt number as functions of the pertinent dimensionless parameters. The flow field is visualized in terms of the streamline patterns adjacent to the surface of the cylinder for a range of values of the Grashof number, Prandtl number and power-law index. A separated flow region forms at as low values of the Prandtl number as Pr = 0.72 for n ? 1 (Newtonian and shear-thickening fluids); whereas for shear-thinning fluids (n < 1), the flow remains attached to the cylinder surface over the range of conditions encompassed here. The bubble size grows with Grashof number and it shrinks with Prandtl number. In order to quantify the deviation from the Newtonian behaviour, the normalized values of average Nusselt number are analyzed as a function of the power-law index. In addition, a correlation is proposed for average Nusselt number as a function of the Grashof number, Prandtl number and power-law index. In general terms, shear-thinning fluid behaviour enhances heat transfer whereas shear-thickening has adverse influence on it.     

Criteria for Predicting Transitions in Free Convection Heat Transfer From Isothermal Convex Bodies in Fluids With Any Prandtl Number: A New Analytical Model

This article presents criteria for predicting transitions in free convection from isothermal convex bodies in fluids with any Prandtl number (conduction–laminar and laminar–turbulent transitions). Laminar–turbulent transition results for vertical plates happening in the vicinity of Grashof number 10⁹ in fluids with any Prandtl number and other geometries show that this criterion based on Grashof number is a better choice compared to that of the Rayleigh number, as the Prandtl number yields a large change in the transition Rayleigh number. Besides, for convex geometries other than vertical plates, the laminar–turbulent transition occurs at values slightly less than Grashof number 10⁹ owing to longitudinal vortices bringing about separation over their surface and helping laminar–turbulent transition to arise sooner. On the other side, the results of the conduction–laminar transition criterion for different convex bodies disclose that the Rayleigh number is superior to the Grashof number. Likewise, the present results for five different geometries over a wide range of Prandtl numbers are compared with the available experimental data, and excellent agreement is found. This reveals that the proposed criteria is powerful enough to predict the two transition phenomena in free convection heat transfer from isothermal convex bodies for a wide range of convex geometries in fluids with any Prandtl number.     

Magnetoconvection in an enclosure with partially active vertical walls

Magnetoconvection of an electrically conducting fluid in a square cavity with partially thermally active vertical walls is investigated numerically. The active part of the left side wall is at a higher temperature than the active part of the right side wall. The top, bottom and the inactive parts of the side walls are thermally inactive. Nine different combinations of the relative positions of the active zones are considered. The governing equations are discretized by the control volume method with QUICK scheme and solved numerically by SIMPLE algorithm for the pressure–velocity coupling together with underrelaxation technique. The results are obtained for Grashof numbers between 10⁴ and 10⁶, Hartmann numbers between 0 and 100 and Prandtl numbers 0.054–2.05. The heat transfer characteristics are presented in the form of streamlines and isotherms. The heat transfer rate is maximum for the middle–middle thermally active locations while it is poor for the top–bottom thermally active locations. The average Nusselt number decreases with an increase of Hartmann number and increases with an increase of Grashof number. For sufficiently large magnetic field Ha = 100 the convective mode of heat transfer is converted into conductive mode in the low region of Grashof number than in the high region.     

Effect of aspect ratio on natural convection in power-law liquids from a heated horizontal elliptic cylinder

The influence of aspect ratio and shear-dependent viscosity on free convection heat transfer from a horizontal heated elliptic cylinder in power-law fluids has been investigated. In particular, the coupled momentum and energy equations have been solved numerically over the following ranges of conditions: Grashof number, 10 ? Gr ? 10⁵; Prandtl number, 0.72 ? Pr ? 100; power-law index, 0.3 ? n ? 1.5 and aspect ratio, 0.2 ? E ? 5. The new extensive results demonstrate the influence of the Grashof number (Gr), Prandtl number (Pr), power-law index (n) and aspect ratio (E) on the macroscopic heat and momentum transfer characteristics like local and average values of Nusselt number (Nu) and drag coefficients (C_D). Further insights are developed by examining the structure of the flow and temperature fields adjacent to the cylinder. Broadly speaking, all else being equal, shear- thinning fluid behaviour promotes heat transfer whereas shear-thickening viscosity has a deleterious effect on it with reference to that in Newtonian fluids. Also, the rate of heat transfer gradually increases as the cylinder shape passes from blunt to slender with respect to the direction of gravity. Finally, the present numerical values of the Nusselt number are correlated using a simple analytical form which facilitates interpolation of the present results for the intermediate values of the governing parameters. The paper is concluded by presenting detailed comparisons with the previous numerical and experimental results available in the literature, especially in Newtonian fluids.     

Magnetoconvection in a square cavity with partially active vertical walls: Time periodic boundary condition

Magnetoconvection of an electrically conducting fluid in a square cavity with partially thermally active sidewalls is investigated numerically. Temperature of one of the thermally active regions of the side walls is periodic in time while the opposite wall is isothermal. The horizontal walls and the remaining parts of the side walls are thermally inactive. Nine different combinations of the relative positions of the active zones are considered. The governing equations are discretized by the control volume method with QUICK scheme and solved numerically by SIMPLE algorithm for the pressure–velocity coupling together with under relaxation technique. The tests were carried out for various values of amplitude, period, Grashof number, Hartmann number and Prandtl number. The heat transfer characteristics are presented in the form of streamlines, isotherms and velocity profiles both for transient and steady state. It is observed that the flow and the heat transfer rate in the cavity are affected by the sinusoidal temperature profile and by the magnetic field at lower values of Grashof number. The rate of heat transfer oscillates for increasing periods but it is maximum for Ω = 3 and it is found to be an increasing function of amplitude but decreases for higher values of Hartmann number. The heat transfer rate is maximum for the middle–middle thermally active locations while it is poor for the top heating and bottom cooling active locations. The average Nusselt number decreases with an increase of Hartmann number and increases with increase of Prandtl number and Grashof number.     

Forced convection heat transfer between permeable rotating slender cylinders and power-law fluids

A new analysis is presented for forced convection of power-law fluids past rotating slender cylinders with fluid suction/injection at the wall. After a suitable coordinate transformation to reduce the complexity of the governing boundary-layer equations, the resulting nonlinear coupled differential equations were solved with an implicit finite difference scheme. Of interest were the effects of transverse curvature, power-law viscosity index, the type of thermal boundary condition and the generalized Prandtl number on the local skin friction coefficient and Nusselt number of rotating cylinders with porous surfaces. In general, body spin enhances convection heat transfer and increases the skin friction coefficient. Fluid injection into the boundary layer reduces and suction increases the local skin friction group. Higher Prandtl numbers increase the local heat transfer group especially in the presence of fluid withdrawal. Dilatant fluids exhibit a distinctively different behavior with respect to the heat transfer group when compared to pseudoplastics. The two thermal boundary conditions yield quite similar results.     

Numerical Analysis of the Effects of Selected Geometrical Parameters and Fluid Properties on MHD Natural Convection Flow in an Inclined Elliptic Porous Enclosure with Localized Heating

Magnetohydrodynamic (MHD) natural convection flow and associated heat convection in an oriented elliptic enclosure has been investigated with numerical simulations. A magnetic field was applied to the cylindrical wall of the configuration, the top and bottom walls of the enclosure were circumferentially cooled and heated, respectively, while the extreme ends along the cross‐section of the elliptic duct were considered adiabatic. The full governing equations in terms of continuity, momentum, and energy transport were transformed into nondimensional form and solved numerically using finite difference method adopting Gauss–Seidel iteration technique. The selected geometrical parameters and flow properties considered for the study were eccentricity (0, 0.2, 0.4, 0.6, and 0.8), angle of inclination (0°, 30°, 60°, and 90°), Hartmann number (0, 25, and 50), Grashof number (10⁴, 10⁵, and 10⁶), and Darcy number (10^?3, 10^?4, and 10^?5). The Prandtl number was held constant at 0.7. Numerical results were presented by velocity distributions as well as heat transfer characteristics in terms of local and average Nusselt numbers (i.e., rate of heat transfer). The optimum heat transfer rate was attained at e value of 0.8. Also, the heat transfer rate increased significantly between the angles of inclination 58° and 90°. In addition, Hartmann number increased with decreased heat transfer rate and flow circulation. A strong flow circulation (in terms of velocity distribution) was observed with increased Grashof and Darcy numbers. The combination of the geometric and fluid properties therefore can be used to regulate the circulation and heat transfer characteristics of the flow in the enclosure.     

Momentum and Heat Transfer from a Semi-Circular Cylinder to Power-Law Fluids in the Vortex Shedding Regime

Flow and heat transfer from a semi-circular cylinder to power-law fluids has been studied in the laminar vortex shedding regime for the range of conditions as follows: Reynolds number, 40 ≤ Re ≤ 140; Prandtl number, 0.7 ≤ Pr ≤ 50; and power-law index, 0.2 ≤ n ≤ 1.8. Extensive results are presented in terms of the streamline, vorticity and isotherm profiles, drag and lift coefficients, Strouhal number, and Nusselt number. The influence of Reynolds number, Prandtl number, and power-law index are seen to be more prominent in shear-thinning fluids (n < 1) than that in shear-thickening (n > 1) fluids. The shear-thinning fluid behavior enhances the rate of heat transfer, whereas shear-thickening fluid behavior impedes it.     

Momentum and heat transfer from a square cylinder in power-law fluids

Numerical solutions are sought, using FLUENT, to the mass, momentum and thermal energy equations for the 2-D flow of power-law fluids over a cylinder of square cross-section. The major thrust of this work is to delineate the values of the Reynolds number denoting the onset of flow separation and the limits of the steady flow regime for both shear-thinning and shear-thickening type fluids. Extensive results are reported on streamline and vorticity contours over wide ranges of power-law index (0.2–1.4) corroborating the occurrence of these two transitions. Having established the limits of the steady flow regime, drag and Nusselt number results are obtained in this regime as functions of the Reynolds number (0.1–40), of Prandtl number (0.7–100) for highly shear-thinning fluids (power-law index < 0.5) thereby extending the range of currently available results to that encountered in practical applications. The Nusselt number shows positive dependence on both the Reynolds and Prandtl numbers. Also, shear-thinning characteristics can augment the rate of heat transfer by up to 100% under appropriate conditions.     

Mixed Convection from a Heated Sphere in Bingham Plastic Fluids

In this work, the steady and laminar mixed-convection heat transfer from an isothermal sphere immersed in Bingham plastic fluids has been investigated in the aiding-buoyancy configuration. The pertinent coupled equations of motion and thermal energy have been solved numerically over the following ranges of conditions: Richardson number, 0 ≤ Ri ≤ 2, Bingham number, 0 ≤ Bn ≤ 10, Reynolds number, 0.1 ≤ Re ≤ 100 and Prandtl number, 10 ≤ Pr ≤ 100. Flow characteristics like streamlines, pressure coefficient, morphology of yielded/unyielded regions and drag coefficient are discussed extensively. Similarly, isotherms, local Nusselt number and average Nusselt number are thoroughly examined to develop an overall understanding of the corresponding heat transfer characteristics. All else being equal, in contrast to the positive role of the aiding-buoyancy free convection in Newtonian and power-law fluids, due to the fluid yield stress, heat transfer is impeded in viscoplastic fluids. While the average value of the Nusselt number is influenced by four dimensionless groups, namely, Reynolds number, Bingham number, Prandtl number and Richardson number, by using novel scaling, it has been possible to consolidate the present results via the use of the Colburn j-factor in a simple form. This is particularly suitable for predicting the value of the Nusselt number in a new application.     

Interaction Between Two Mutually Orthogonal Heat-Generating Baffles in a Square Cavity

Laminar free convection induced by two mutually orthogonal discrete heat-generating baffles in a two-dimensional square cavity is analyzed numerically. The computations were carried out for different locations and combinations of heat source strengths of the baffles for a fixed Grashof number of 10⁶and Prandtl number of 0.71. The coupled governing equations were solved bya finite-difference method using alternating direction implicit technique and successive overrelaxation methods. The obtained results clearly show that the hydrodynamic and thermal fields in the cavity depend on both the location and strength of the heat-generating baffles. Though the flow inhibition is caused by both the baffles, the baffle with higher source strength plays a decisive role in inducing the flow. The locations of baffles with unequal source strengths produce significant changes in the net heat transfer rate. This is further magnified for higher contrast in source strengths. This study provides qualitative suggestions that may improve the thermal design of sealed enclosures, which are encountered frequently in the electronics industry.     

A numerical study of mixed convection in a two-sided lid-driven tall cavity containing a heated triangular block for non-Newtonian power-law fluids

The study of hydrodynamics and thermal characteristics inside a lid-driven cavity has been one of the most captivating problems in computational fluid dynamics. In this numerical work, the mixed convection phenomenon inside a two-dimensional, tall lid-driven cavity with top and bottom lids moving in opposite directions, +x and –x, respectively, has been explored for non-Newtonian power-law fluids. The cavity contains a uniformly heated equilateral triangular obstacle at its geometric center.  Numerical experimentation is performed for a range of flow governing parameters, such as aspect ratio (0.25, 0.5, and 0.75), Prandtl number (1, 50, and 100) Richardson number (0.1, 1, and 10), power-law index (0.6–1.4) and Grashof number of 10⁴. The physical perceptions of the cavity are explained by using streamline and isotherm contours. The fluid movement is limited adjacent to the moving wall concerning the Richardson number at the lower Prandtl number. With a rise in the aspect ratio of the cavity, the flow-pattern becomes more dispersed inside the cavity. Heat transfer enhancement is observed at a lower aspect ratio equal to 0.25.     

Free Convective Flow of Electrically Conducting and Viscous Immiscible Fluid Flow in a Vertical Channel in the Presence of First‐Order Chemical Reaction

In this paper, the effects of chemical reaction on free convective flow of electrically conducting and viscous incompressible immiscible fluids are analyzed. The coupled nonlinear equations governing the heat and mass transfer are solved analytically and numerically with appropriate boundary conditions for each fluid and the solutions have been matched at the interface. The analytical solutions are solved by using regular perturbation method valid for small values of perturbation parameter and numerically by using finite difference method. The numerical results for various values of thermal Grashof number, mass Grashof number, Hartman number, viscosity ratio, width ratio, conductivity ratio, and chemical reaction parameter have been presented graphically in the presence and in the absence of electric field load parameter. In addition, the closed form expression for volumetric flow rate, Nusselt number, species concentration, and total heat rate added to the flow is also analyzed. The solutions obtained by finite difference method and perturbation method agree very well to the order of 10^?4 for small values of perturbation parameter.     

Confined Flow and Heat Transfer Phenomena of Non-Newtonian Shear-Thinning Fluids Across a Pair of Tandem Triangular Bluff Bodies

Extensive numerical computations for the confined flow and heat transfer of shear-thinning fluids over two long tandem triangular bluff bodies are performed for Reynolds number (Re) = 1–40, power-law index (n) = 0.2–1, and gap ratio (S/B) = 1–4 for a fixed blockage ratio and Prandtl number of 25% and 50, respectively. The values of critical Re are obtained for each of the two triangular cylinders for all S/B and n. Augmentation in heat transfer for the two tandem triangular cylinders is calculated with respect to both S/B and n.     

Momentum and heat transfer characteristics of a semi-circular cylinder immersed in power-law fluids in the steady flow regime

The continuity, momentum and energy equations describing the flow and heat transfer of power-law fluids over a semi-circular cylinder have been solved numerically in the two-dimensional steady flow regime. The influence of the Reynolds number (Re), Prandtl number (Pr) and power-law index (n) on the local and global flow and heat characteristics have been studied over wide ranges of conditions as follows: 0.01 ? Re ? 30, 1 ? Pr ? 100 and 0.2 ? n ? 1.8. The variation of drag coefficient and Nusselt number with the Reynolds number, Prandtl number and power-law index is shown over the aforementioned ranges of conditions. In addition, streamline and isotherm profiles along with the recirculation length and distribution of pressure coefficient and Nusselt number over the surface of the semi-circular cylinder are also presented to gain further insights into the nature of the underlying kinematics. The wake size (recirculation length) shows almost linear dependence on the Reynolds number (Re ? 1) for all values of power-law index studied herein. The drag values show the classical inverse variation with the Reynolds number, especially for shear-thinning fluids at low Reynolds numbers. The point of maximum pressure coefficient is found slightly displaced from the front stagnation point for highly shear-thinning fluids, whereas for shear-thickening and Newtonian fluids, it coincides with the front stagnation point. For fixed values of the Prandtl number and Reynolds number, the rate of heat transfer decreases with the gradual increase in power-law index; this effect is particularly striking at high Prandtl numbers due to the thinning of the thermal boundary layer. Conversely, as expected, shear-thinning behavior facilitates heat transfer and shear-thickening impedes it. The effect of power-law index on both momentum and heat-transfer characteristics is seen to be appreciable at low Reynolds numbers and it gradually diminishes with the increasing Reynolds number.     

Conjugate natural convection in an enclosure with a heat source of constant heat transfer rate

A numerical study of two-dimensional transient natural convection in a rectangular enclosure having finite thickness heat-conducting walls with a heat source of constant heat transfer rate located on the inner side of the left wall in conditions of convection–radiation heat exchange with an environment on one of the external boundaries has been performed. Mathematical simulation has been carried out in terms of the dimensionless variables such as stream function – vorticity – temperature. Stream function, vorticity and energy equations have been solved by finite difference numerical method. The relevant governing parameters were: the Grashof number from 10⁶ to 10⁸, the Prandtl number, Pr = 0.7 and the conductivity ratio. Detailed results including streamlines and temperature profiles have been obtained.