Numerical Calculation of Three-Dimensional Turbulent Natural Convection in a Cubical Enclosure Using a Two-Equation Model for Turbulence

[+] Author and Article Information
H. Ozoe, A. Mouri, M. Hiramitsu

Department of Industrial and Mechanical Engineering, School of Engineering, Okayama University, Okayama 700 Japan

S. W. Churchill

Department of Chemical Engineering, University of Pennsylvania, Philadelphia, PA 19104 U.S.A.

N. Lior

Department of Mechanical Engineering and Applied Mechanics, University of Pennsylvania, Philadelphia, PA 19104 U.S.A.

J. Heat Transfer 108(4), 806-813 (Nov 01, 1986) (8 pages) doi:10.1115/1.3247016 History: Received February 07, 1985; Online October 20, 2009


This paper presents a model and numerical results for turbulent natural convection in a cubical enclosure heated from below, cooled on a portion of one vertical side wall and insulated on all other surfaces. Three-dimensional balances were derived for material, energy, and the three components of momentum, as well as for the turbulent kinetic energy k and the rate of dissipation of turbulent kinetic energy ε. The constants used in the model were the same as those used by Fraikin et al. for two-dimensional convection in a channel. Illustrative transient calculations were carried out for Ra = 106 and 107 and Pr = 0.7. Both the dominant component of the vector potential and the Nusselt number were found to converge to a steady state. Isothermal lines and velocity vectors for vertical cross sections normal to the cooled wall indicated three-dimensional effects near the side walls. A top view of the velocity vectors revealed a downward spiral flow near the side walls along the cooled vertical wall. A weak spiral flow was also found along the side walls near the wall opposing the partially cooled one. The highest values of the eddy diffusivity were 2.6 and 5.8 times the molecular kinematic viscosity for Ra = 106 and 107 , respectively. A coaxial double spiral movement, similar to that previously reported for laminar natural convection, was found for the time-averaged flow field. This computing scheme is expected to be applicable to other thermal boundary conditions.

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





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