Research Papers: Natural and Mixed Convection

Convective Motion and Heat Transfer in a Slowly Rotating Fluid Quasi-Sphere With Uniform Heat Source and Axial Gravity

[+] Author and Article Information
Gerardo Anguiano-Orozco

Facultad de Ciencias,
UAEM, Av. Instituto Literario 100,
50000 Toluca, Mexico;
Depto. de Física,
Instituto Nacional
de Investigaciones Nucleares,
Apdo. Postal 18-1027,
Mexico D.F.
e-mail: gerardo.anguiano@gmail.com

Ruben Avila

Departamento de Termofluidos,
Facultad de Ingeniería,
Universidad Nacional Autónoma de México,
C.P. 04510,
Mexico D.F.
e-mail: ravila@unam.mx

Syed Shoaib Raza

Pakistan Institute of Engineering
and Applied Sciences (PIEAS),
P.O. Nilore,
Islamabad 45650, Pakistan
e-mail: fac246@pieas.edu.pk;

1In Memoriam, 2010.

2Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the Journal of Heat Transfer. Manuscript received November 23, 2010; final manuscript received November 29, 2012; published online March 20, 2013. Assoc. Editor: Darrell W. Pepper.

J. Heat Transfer 135(4), 042501 (Mar 20, 2013) (10 pages) Paper No: HT-10-1549; doi: 10.1115/1.4023126 History: Received November 23, 2010; Revised November 29, 2012

The laminar natural convection of a rotating fluid quasi-sphere in the presence of an axial gravity field and uniform heat source is presented. The influence of the Rayleigh number Ra and the Taylor number Ta on the flow pattern and heat transfer rate from the fluid quasi-sphere is discussed. The governing nonsteady, three-dimensional Navier–Stokes equations for an incompressible fluid, formulated in a Cartesian coordinate system, have been numerically solved by using the h/p spectral element method. It is shown that for a given Ta number, as the Ra number is increased, the heat transfer on the northern hemisphere is enhanced whereas the average Nusselt number on the southern hemisphere is reduced. On the other hand for a given Ra number, as the Ta number is increased, the heat transfer is a function of the convective motion intensity. It has been found that for low and high Ra numbers the heat transfer rate slightly depends on the rotation rate. However at intermediate Ra numbers, the net effect of an increased rotation rate is a reduction of the heat transfer through the wall, hence an increase of the maximum temperature of the fluid sphere is observed. We show that the net effect of the Coriolis force is to damp the convective motion and to allow a redistribution of the vorticity field.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

Natural convection and heat transfer in a rotating fluid sphere with uniform heat source and axial gravity. The origin of the Cartesian coordinate system is located at the center of the sphere. The dimensional constant angular velocity, axial gravity and uniform heat source are defined by Ω = Ω3i3, g = g3i3, and Qv, respectively.

Grahic Jump Location
Fig. 2

Spectral element method mesh with 256 nonregular hexahedra macroelements with polynomial expansion of order 6. The origin of the Cartesian coordinate system is located at the center of the quasi-sphere. Left panel: discretization of the bounding spherical surface. Right panel: an interior view of the mesh.

Grahic Jump Location
Fig. 3

SEM GLL points distribution in the quasi-sphere. Left panel: x2 - x3 meridional plane. Right panel: GLL points on the bounding spherical surface at three positions of the x1 axis (directed out of the page), (i) external circle at x1 = 0 (meridional plane), (ii) middle circle at x1 = 0.7, and (iii) internal circle at x1 = 0.91. Circles: SEM-GLL points. Continuous line: perfect sphere with radius 0.985. Order of the polynomial equal to 6.

Grahic Jump Location
Fig. 4

Nonsteady heat conduction within a quasi-sphere with homogeneous and constant heat source. Dimensionless temperature distribution along the radial direction. Left panel: polynomial interpolation p of order 6. Right panel: polynomial interpolation p of order 8. Continuous line: analytical solution for a perfect sphere, see Eqs. (6) and (7). Circles: SEM method results. At t = 1, the steady state condition is reached. Dimensional values: ΔT = Qv R2/6k = 161.7 K, Qv = 1000 W/m3, R = 0.985 m, k = 1 W/m-K, α = 1 m2/s.

Grahic Jump Location
Fig. 5

Natural convection in a fluid sphere for different Ra numbers without rotation Ta = 0. Meridional fields of vorticity (left column), pressure (middle column), and temperature (right column). (a): Ra = 10, (b): Ra = 1 × 103, (c): Ra = 1 × 104, (d): Ra = 1 × 105, (e): Ra = 5 × 105, (f): Ra = 1 × 106, and (g): Ra = 1 × 107.

Grahic Jump Location
Fig. 6

Dimensionless maximum temperature Tmax in the fluid sphere (left panel) and average Nusselt number Nu¯ (right panel) as functions of the Taylor number Ta and the Rayleigh number Ra. The convection coefficient h, used to obtain Nu¯, has been calculated from the energy balance (see Eqs. (12) and (13)). (i) ○ Ra = 10, (ii) * Ra = 1 × 103, (iii) □ Ra = 1 × 104, (iv) Δ Ra = 1 × 105, (v) ⋆ Ra = 1 × 106, and (vi) ⋄ Ra = 1 × 107.

Grahic Jump Location
Fig. 7

Average Nusselt number as a function of the Ra and Ta numbers. The convection coefficient h, used to obtain Nu¯, has been calculated from the energy balance (see Eqs. (12) and (13)). (i) ○ Ta = 0, (ii) □ Ta = 1600, (iii) Δ Ta = 6400, and (iv) ⋆ Ta = 14,400.

Grahic Jump Location
Fig. 8

Meridional convective vorticity fields with rotation. First column (left): Ta = 0, second column: Ta = 1600, third column: Ta = 6400, and fourth column (right): Ta = 14,400. (a): Ra = 10, (b): Ra = 1 × 103, (c): Ra = 1 × 104, (d): Ra = 1 × 105, (e): Ra = 1 × 106, and (f): Ra = 1 × 107.

Grahic Jump Location
Fig. 9

Meridional convective temperature fields with rotation. Same caption as Fig. 8.

Grahic Jump Location
Fig. 10

Dimensionless meridional (x2, x3) vorticity (i1 component) in a fluid (quasi sphere) for Ra = 10. Left panel: Ta = 1600. Middle panel: Ta = 6400. Right panel: Ta = 14,400. Bold line represents the average vorticity equal to zero, and a measure of the thickness of the Ekman boundary layer δE.

Grahic Jump Location
Fig. 11

Dimensionless Ekman boundary layer thickness δE, and thermal boundary layer thickness δT in terms of the Ra number. Left panel: Ta = 6400, δE ≈ 0.136 (continuous line). Right panel: Ta = 14,400, δE ≈ 0.113 (continuous line). ○: δT at the north pole region, □: δT at the equatorial region.




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