Research Papers: Forced Convection

Direct Numerical Simulation of Forced Convective Heat Transfer From a Heated Rotating Sphere in Laminar Flows

[+] Author and Article Information
Zhi-Gang Feng

Department of Mechanical Engineering,
University of Texas at San Antonio,
San Antonio, TX 78249
e-mail: zhigang.feng@utsa.edu

1Corresponding author.

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received July 3, 2013; final manuscript received December 8, 2013; published online February 18, 2014. Assoc. Editor: William P. Klinzing.

J. Heat Transfer 136(4), 041707 (Feb 18, 2014) (8 pages) Paper No: HT-13-1334; doi: 10.1115/1.4026307 History: Received July 03, 2013; Revised December 08, 2013

The laminar forced convection of a heated rotating sphere in air has been studied using a three-dimensional immersed boundary based direct numerical simulation method. A regular Eulerian grid is used to solve the modified momentum and energy equations for the entire flow region simultaneously. In the region that is occupied by the rotating sphere, a moving Lagrangian grid is used, which tracks the rotational motion of the particle. A force density function or an energy density function is introduced to represent the momentum interaction or thermal interaction between the sphere and fluid. This numerical method is validated by comparing simulation results with analytical solutions of heat diffusion problem and other published experimental data. The flow structures and the mean Nusselt numbers for flow Reynolds number ranging from 0 to 1000 are obtained. We compared our simulation results of the mean Nusselt numbers with the correlations from the literature and found a good agreement for flow Reynolds number greater than 500; however, a significant discrepancy arises at flow Reynolds number below 500. This leads us to develop a new equation that correlates the mean Nusselt number of a heated rotating sphere for flows of 0Re500.

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


Lien, F. S., Chen, C. K., and Cleaver, J. W., 1986, “Mixed and Free Convection Over a Rotating Sphere With Blowing and Suction,” ASME J. Heat Transfer, 108, pp. 398–404. [CrossRef]
Palec, G. L., and Daguenet, M., 1987, “Laminar Three-Dimensional Mixed Convection About a Rotating Sphere in Stream,” Int. J. Heat Mass Transfer, 30, pp. 1511–1523. [CrossRef]
Kreith, F., Roberts, L. G., Sullivan, J. A., and Sinha, S. N., 1965, “Convection Heat Transfer Flow Phenomena of Rotating Sphere,” Int. J. Heat Mass Transfer, 6, pp. 881–895. [CrossRef]
Tieng, S. M., and Yan, A. C., 1993, “Experimental Investigation on Convective Heat Transfer of Heated Rotating Sphere,” Int. J. Heat Mass Transfer, 36, pp. 599–610. [CrossRef]
Zheng, G., and List, R., 1996, “Convective Heat Transfer of Rotating Spheres and Spheroids With Non-Uniform Surface Temperatures,” Int. J. Heat Mass Transfer, 39, pp. 1815–1826. [CrossRef]
Farouk, B., 1985, “Mixed Convective Flow Around a Slowly Rotating Isothermal Sphere,” ASME J. Heat Transfer, 107, pp. 431–438. [CrossRef]
Dorfman, L. A., and Serazetdinov, A. Z., 1973, “Laminar Flow and Heat Transfer Near Rotating Axisymmetric Surface,” Int. J. Heat Mass Transfer, 8, pp. 317–327. [CrossRef]
Banks, W. H. H., 1965, “The Thermal Laminar Boundary Layer on Rotating Sphere,” J. Appl. Math. Phys., 16, pp. 780–788. [CrossRef]
Antar, M. A., and El-Shaarawi, M. A. I., 2004, “Simultaneous Effect of Rotation and Natural Convection on the Flow About a Liquid Sphere,” Heat Mass Transfer, 40, pp. 393–399. [CrossRef]
Feng, Z.-G., and Michaelides, E. E., 2008, “Inclusion of Heat Transfer Computations for Particle Laden Flows,” Phys. Fluids, 20, pp. 675–684. [CrossRef]
Feng, Z.-G., and Michaelides, E. E., 2009, “Heat Transfer in Particulate Flows With Direct Numerical Simulation (DNS),” Int. J. Heat Mass Transfer, 52, pp. 777–786. [CrossRef]
Peskin, C., 1977, “Numerical Analysis of Blood Flow in the Heart,” J. Comput. Phys, 25, pp. 220–252. [CrossRef]
Mohd-Yusof, J., 1997, “Combined Immersed Boundaries/b-Splines Methods for Simulations of Flows in Complex Geometries,” Annual Research Briefs, Center for Turbulence Research, Stanford University, Stanford, CA, Vol. 161, pp. 35–60.
Feng, Z.-G., and Michaelides, E. E., 2005, “Proteus: A Direct Forcing Method in the Simulations of Particulate Flow,” J. Comput. Phys., 202, pp. 20–51. [CrossRef]
Uhlmann, M., 2005, “An Immersed Boundary Method With Direct Forcing for the Simulation of Particulate Flows,” J. Comput. Phys., 209, pp. 448–476. [CrossRef]
Feng, Z.-G., and Michaelides, E. E., 2010, “Robust Treatment of No-Slip Boundary Condition and Velocity Updating for the Lattice-Boltzmann Simulation of Particulate Flows,” Comput. Fluids, 38, pp. 370–381. [CrossRef]
Wang, Z., Fan, J., Luo, K., and Chen, K., 2009, “Immersed Boundary Method for the Simulation of Flows With Heat Transfer,” Int. J. Heat Mass Transfer, 52, pp. 4510–4518. [CrossRef]
Kang, S. K., and Hassan, Y. A., 2011, “A Direct-Forcing Immersed Boundary Method for the Thermal Lattice Boltzmann Method,” Comput. Fluids, 49, pp. 36–45. [CrossRef]
Deen, N., Kriebitzsch, S., van der Hoef, M., and Kuipers, J., 2012, “Direct Numerical Simulation of Flow and Heat Transfer in Dense Fluid-Particle Systems,” Chem. Eng. Sci., 81, pp. 329–344. [CrossRef]
Wachs, A., 2011, “Rising of 3D Catalyst Particles in a Natural Convection Dominated Flow by a Parallel DNS Method,” Comput. Chem. Eng., 35, pp. 2169–2185. [CrossRef]


Grahic Jump Location
Fig. 1

Conceptual DNS-IB model of a solid particle immersed in a fluid

Grahic Jump Location
Fig. 2

Setup of Cartesian coordinates for a heated sphere of radius a rotating at angular velocity ω in a viscous fluid

Grahic Jump Location
Fig. 3

Comparisons of the mean Nusselt number of the unsteady diffusion of a sphere obtained by analytical solution and numerical simulation

Grahic Jump Location
Fig. 4

Snapshots of temperature contours from both numerical (right half side of each plot) and analytical solution (left half side of each plot) at time t = 0.25 (left), t = 0.75 (middle), and t = 1.5 (right)

Grahic Jump Location
Fig. 5

Simulation results and exact solution of dimensionless temperature in the radial direction

Grahic Jump Location
Fig. 11

Simulation results of time evolution of the mean Nusselt numbers for Re≤1000

Grahic Jump Location
Fig. 6

The effect of grid spacing to the simulated mean Nusselt number at Re = 500

Grahic Jump Location
Fig. 7

The effect of time step δt to the simulated mean Nusselt number at Re = 500

Grahic Jump Location
Fig. 10

Isosurfaces for dimensionless temperature ranging from 1.0 to 0.1 at t = 50

Grahic Jump Location
Fig. 8

Visualizations of a rotating sphere in air at Re = 500. (a) Isosurface for fluid velocity magnitude of v = 0.2 colored in terms of temperature scale; (b) top view of flow streamlines; and (c) side view of flow streamlines.

Grahic Jump Location
Fig. 9

Velocity vectors on two slices passing through the center of the sphere at t = 50. (a) a horizontal slice or equatorial plane and (b) a vertical slice.

Grahic Jump Location
Fig. 12

Flow velocity magnitude contours (top row) and temperature contours (bottom row) on vertical slices (x-y plane) at t = 40. Left: Re = 100; middle: Re = 200; right: Re = 500.

Grahic Jump Location
Fig. 13

Comparison of present simulation results with correlations found in the literature and the new correlation for Re≤500



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