0
Research Papers: Natural and Mixed Convection

Laminar Natural Convection of Bingham Fluids in Inclined Differentially Heated Square Enclosures Subjected to Uniform Wall Temperatures

[+] Author and Article Information
Şahin Yİğİt

School of Mechanical and Systems Engineering,
Newcastle University,
Newcastle-Upon-Tyne NE1 7RU, UK
Department of Mechanical Engineering,
Karadeniz Technical University,
Trabzon 61080, Turkey
e-mail: syigit@ktu.edu.tr

Robert J. Poole

School of Engineering,
University of Liverpool,
Brownlow Hill,
Liverpool L69 3GH, UK
e-mail: robpoole@liv.ac.uk

Nilanjan Chakraborty

School of Mechanical and Systems Engineering,
Newcastle University,
Newcastle-Upon-Tyne NE1 7RU, UK
e-mail: nilanjan.chakraborty@ncl.ac.uk

As will be discussed in Sec. 3, the findings based on the current numerical experiment for Pr = 500 are qualitatively valid for other values of nominal Prandtl number and the resulting correlations capture the effects of Pr from 10 to 500.

R2 value is greater than 0.99 for all Ra values for the correlation given by Eq. (11).

The R2 value is greater than 0.99 for the correlation given by Eqs. (14b) and (14c). The maximum percentage error for the predictions of Eqs. (14b) and (14c) with respect to simulation data is found to be 4.8%.

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received February 7, 2014; final manuscript received January 19, 2015; published online March 3, 2015. Assoc. Editor: Oronzio Manca.

J. Heat Transfer 137(5), 052504 (May 01, 2015) (12 pages) Paper No: HT-14-1068; doi: 10.1115/1.4029763 History: Received February 07, 2014; Revised January 19, 2015; Online March 03, 2015

The effects of inclination 180degφ0deg on steady-state laminar natural convection of yield-stress fluids, modeled assuming a Bingham approach, have been numerically analyzed for nominal values of Rayleigh number Ra ranging from 103 to 105 in a square enclosure of infinite span lying horizontally at φ=0deg, then rotated about its axis for φ>0deg cases. It has been found that the mean Nusselt number Nu¯ increases with increasing values of Rayleigh number but Nu¯ values for yield-stress fluids are smaller than that obtained in the case of Newtonian fluids with the same nominal value of Rayleigh number Ra due to the weakening of convective transport. For large values of Bingham number Bn (i.e., nondimensional yield stress), the mean Nusselt number Nu¯ value settles to unity (Nu¯=1.0) as heat transfer takes place principally due to thermal conduction. The mean Nusselt number Nu¯ for both Newtonian and Bingham fluids decreases with increasing φ until reaching a local minimum at an angle φ* before rising with increasing φ until φ=90deg. For φ>90deg the mean Nusselt number Nu¯ decreases with increasing φ before assuming Nu¯=1.0 at φ=180deg for all values of Ra. The Bingham number above which Nu¯ becomes unity (denoted Bnmax) has been found to decrease with increasing φ until a local minimum is obtained at an angle φ* before rising with increasing φ until φ=90deg. However, Bnmax decreases monotonically with increasing φ for 90deg<φ<180deg. A correlation has been proposed in terms of φ, Ra, and Bn, which has been shown to satisfactorily capture Nu¯ obtained from simulation data for the range of Ra and φ considered here.

FIGURES IN THIS ARTICLE
<>
Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.

References

Catton, I., Ayyaswamy, P. S., and Clever, R. M., 1974, “Natural Convection Flow in a Finite, Rectangular Slot Arbitrarily Oriented With Respect to the Gravity Vector,” Int. J. Heat Mass Transfer, 17(2), pp. 173–184. [CrossRef]
Arnold, J. N., Catton, I., and Edwards, D. K., 1976, “Experimental Investigation of Natural Convection in Inclined Rectangular Regions of Differing Aspect Ratios,” ASME J. Heat Transfer, 98(1), pp. 67–71. [CrossRef]
Ozoe, H., Sayama, H., and Churchill, S., 1975, “Natural Convection in an Inclined Rectangular Channel at Various Aspect Ratios and Angles—Experimental Measurements,” Int. J. Heat Mass Transfer, 18(12), pp. 1425–1431. [CrossRef]
Ozoe, H., Sayama, H., and Churchill, S., 1974, “Natural Convection in an Inclined Square Channel,” Int. J. Heat Mass Transfer, 17(3), pp. 401–406. [CrossRef]
Soong, C. Y., Tzeng, P. Y., Chiang, D. C., and Sheu, T. S., 1996, “Numerical Study on Mode-Transition of Natural Convection in Differentially Heated Inclined Enclosures,” Int. J. Heat Mass Transfer, 39(14), pp. 2869–2882. [CrossRef]
Wang, H., and Hamed, M., 2006, “Flow Mode-Transition of Natural Convection in Inclined Rectangular Enclosures Subjected to Bidirectional Temperature Gradients,” Int. J. Therm. Sci., 45(8), pp. 782–795. [CrossRef]
Khezzar, L., Signiner, D., and Vinogradov, I., 2012, “Natural Convection of Power Law Fluids in Inclined Cavities,” Int. J. Therm. Sci., 53(3), pp. 8–17. [CrossRef]
Bodenschatz, E., Pesch, W., and Ahlers, G., 2000, “Recent Developments in Rayleigh–Bénard Convection,” Annu. Rev. Fluid Mech., 32(1), pp. 709–778. [CrossRef]
Vola, D., Boscardin, L., and Latché, J. C., 2003, “Laminar Unsteady Flows of Bingham Fluids: A Numerical Strategy and Some Benchmark Results,” J. Comput. Phys., 187(2), pp. 441–456. [CrossRef]
Turan, O., Chakraborty, N., and Poole, R. J., 2010, “Laminar Natural Convection of Bingham Fluids in a Square Enclosure With Differentially Heated Side Walls,” J. Non-Newtonian Fluid Mech., 165(15–16), pp. 901–913. [CrossRef]
Vikhansky, A., 2010, “On the Onset of Natural Convection of Bingham Liquid in Rectangular Enclosures,” J. Non-Newtonian Fluid Mech., 165(23–24), pp. 1713–1716. [CrossRef]
Turan, O., Poole, R. J., and Chakraborty, N., 2011, “Aspect Ratio Effects in Laminar Natural Convection of Bingham Fluids in Rectangular Enclosures With Differentially Heated Side Walls,” J. Non-Newtonian Fluid Mech., 166(3–4), pp. 208–230. [CrossRef]
Turan, O., Sachdeva, A., Poole, R. J., and Chakraborty, N., 2011, “Laminar Natural Convection of Bingham Fluids in a Square Enclosure With Vertical Walls Subjected to Constant Heat Flux,” Numer. Heat Transfer, Part A, 60(5), pp. 381–409. [CrossRef]
Turan, O., Poole, R. J., and Chakraborty, N., 2014, “Influences of Boundary Conditions on Laminar Natural Convection of Bingham Fluids in Rectangular Enclosures With Differentially Heated Side Walls,” Heat Transfer Eng., 35(9), pp. 822–849. [CrossRef]
Balmforth, N. J., and Rust, A. C., 2009, “Weakly Nonlinear Viscoplastic Convection,” J. Non-Newtonian Fluid Mech., 158(1–3), pp. 36–45. [CrossRef]
Vikhansky, A., 2009, “Thermal Convection of a Viscoplastic Liquid With High Rayleigh and Bingham Numbers,” Phys. Fluids, 21(10), p. 103103. [CrossRef]
Zhang, J., Vola, D., and Frigaard, I. A., 2006, “Yield Stress Effects on Rayleigh–Bénard Convection,” J. Fluid Mech., 566(11), pp. 389–419. [CrossRef]
Park, H. M., and Ryu, D. H., 2001, “Rayleigh–Bénard Convection of Viscoelastic Fluids in Finite Domains,” J. Non-Newtonian Fluid Mech., 98(2–3), pp. 169–184. [CrossRef]
Turan, O., Chakraborty, N., and Poole, R. J., 2012, “Laminar Rayleigh–Bénard Convection of Yield Stress Fluids in a Square Enclosure,” J. Non-Newtonian Fluid Mech., 171–172(3), pp. 83–96. [CrossRef]
Turan, O., Poole, R. J., and Chakraborty, N., 2012, “Boundary Condition Effects on Natural Convection of Bingham Fluids in a Square Enclosure With Differentially Heated Horizontal Walls,” J. Comput. Therm. Sci., 4(1), pp. 77–97. [CrossRef]
Hassan, M. A., Pathak, M., and Khan, M. K., 2013, “Natural Convection of Viscoplastic Fluids in a Square Enclosure,” ASME J. Heat Transfer, 135(12), p. 122501. [CrossRef]
Darbouli, M., Métivier, C., Piau, J. M., Magnin, A., and Abdelali, A., 2013, “Rayleigh–Bénard Convection for Viscoplastic Fluids,” Phys. Fluids., 25(2), p. 023101. [CrossRef]
Kebiche, Z., Castelain, C., and Burghelea, T., 2014, “Experimental Investigation of the Rayleigh–Bénard Convection in a Yield Stress Fluid,” J. Non-Newtonian Fluid Mech., 203, pp. 9–23. [CrossRef]
Barnes, H. A., 1999, “The Yield Stress—A Review or ‘παντα ρɛι’—Everything Flows?,” J. Non Newtonian Fluid Mech., 81(1–2), pp. 133–178. [CrossRef]
O'Donovan, E. J., and Tanner, R. I., 1984, “Numerical Study of the Bingham Squeeze Film Problem,” J. Non-Newtonian Fluid Mech., 15(1), pp. 75–83. [CrossRef]
Gebhart, B., 1962, “Effects of Viscous Dissipation in Natural Convection,” J. Fluid Mech., 14(2), pp. 225–232. [CrossRef]
Patankar, S. V., 1980, Numerical Heat Transfer and Fluid Flow, Hemisphere, Washington, DC.
de Vahl Davis, G., 1983, “Natural Convection of Air in a Square Cavity: A Bench Mark Numerical Solution,” Int. J. Numer. Methods Fluids, 3(3), pp. 249–264. [CrossRef]
Mitsoulis, E., and Zisis, T., 2001, “Flow of Bingham Plastics in a Lid-Driven Square Cavity,” J. Non-Newtonian Fluid Mech., 101(1–3), pp. 173–180. [CrossRef]
Bejan, A., 1984, Convection Heat Transfer, Wiley, New York. [CrossRef]
Chang, B.-H., 2014, “Numerical Study of Flow and Heat Transfer in Differentially Heated Enclosures,” Therm. Sci., 18(2), pp. 451–463. [CrossRef]

Figures

Grahic Jump Location
Fig. 2

Variations of nondimensional temperature θ (1st column) and U=u1L/α (2nd column) with x2/L at x1/L = 0.5 for (a) φ=45 deg, (b) φ=60 deg, (c) φ=120 deg, and (d) φ=135 deg for different values of Ra at Pr = 500 and Bn=0.05

Grahic Jump Location
Fig. 4

Contours of nondimensional stream function (Ψ = ψ/α) and unyielded zones (gray) for φ=45,60,120,135 deg at Pr = 500 and Bn=0.05 for Ra = (a) 1 × 103, (b) 5 × 103, (c) 1 × 104, (d) 5 × 104, and (e) 1 × 105

Grahic Jump Location
Fig. 8

(a) Variation of Nu¯ with φ for different values of Ra for Newtonian fluids along with the predictions of Eqs. (10) and (11). Contours of (b) nondimensional stream functions (Ψ = ψ/α) and (c) nondimensional temperature θ for Newtonian fluid (Bn=0) for different values of φ around φ=φ* for Ra = 105 at Pr = 500. The inclination φ at which the minimum Nu¯ is obtained for a given set of Ra and Pr is highlighted with an asterisk * in Figs. 8(b) and 8(c). (d) Variation of Bnmax with φ for different values of Ra along with the prediction of Eq. (13).

Grahic Jump Location
Fig. 1

Schematic diagram of the simulation domain

Grahic Jump Location
Fig. 3

Contours of θ at φ=45,60,120,135 deg at Pr = 500 and Bn=0.05 for Ra = (a) 1 × 103, (b) 5 × 103, (c) 1 × 104, (d) 5 × 104, and (e) 1 × 105

Grahic Jump Location
Fig. 5

Variations of θ (1st column) and U=u1L/α (2nd column) with x2/L at x1/L = 0.5 for different values of Bingham numbers for (a) φ=30 deg, (b) φ=45 deg, (c) φ=60 deg, (d) φ=120 deg, (e) φ=135 deg, and (f) φ=150 deg at Pr = 500 and Ra = 105

Grahic Jump Location
Fig. 6

Contours of θ and Ψ with unyielded zones (gray) for φ=45,135 deg at Ra=105 and Pr = 500 for different values of Bn

Grahic Jump Location
Fig. 7

Variation of Nu¯ with Bn for Ra = (a) 5 × 103, (b) 1 × 104, (c) 5 × 104, and (d) 1 × 105 at Pr = 500 and φ=30,60,120,150 deg along with the prediction of Eq. (16d).

Tables

Errata

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In