A 2D1/2 model for low Prandtl number convection in an enclosure

Numerical simulations of stationary thermal convection in a differentially heated enclosure corresponding to the AFRODITE solidification benchmark experiment 1, 2, 3, 4 and 5 are presented. The cavity of relative dimensions 10:6:1 (length:height:width) is characterized by a small transverse width. T...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:International journal of thermal sciences 2013-09, Vol.71, p.53-60
Hauptverfasser: Botton, V, Boussaa, R, Debacque, R, Hachani, L, Zaidat, K, Hadid, HBen, Fautrelle, Y, Henry, D
Format: Artikel
Sprache:eng
Schlagworte:
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:Numerical simulations of stationary thermal convection in a differentially heated enclosure corresponding to the AFRODITE solidification benchmark experiment 1, 2, 3, 4 and 5 are presented. The cavity of relative dimensions 10:6:1 (length:height:width) is characterized by a small transverse width. The Prandtl number Pr is varied within the range [0.0045, 0.03], typical of liquid metals, whereas the Grashof number, defined as Gr = g beta ( Delta [thetas]/L)H4/ nu 2, is varied within the range [1.3 106, 1.6 107]. As shown by the reference 3D simulations, the temperature field in these situations is 2D (independent of the transverse direction); 2D simulations are, however, not able to catch the physics of the flow and the resulting temperature results are also erroneous. To improve these 2D simulation results while keeping reasonable computational times, a 2D1/2 model is developed, which will take into account the no-slip condition at the side walls. This model is obtained by averaging the governing equations over the width of the domain, with a transverse profile for the velocity featuring a uniform central part and two boundary layers of size delta ( delta is fixed for the whole domain). The relative deviation of the temperature field between the 2D1/2 and 3D computations is investigated as a function of the Prandtl number, the Grashof number and the chosen boundary layer thickness. It is shown that an optimum value exists for delta , which gives a mean deviation in the middle plane of less than 2%, whereas the choice of a more usual parabolic profile would lead to a twice larger deviation. Good comparisons are also obtained with the original experimental results reported at the end of the paper. The 2D1/2 model is thus able to give results which compare well with fully 3D results. It can then be used for extensive parametric studies at a reasonable cost.
ISSN:1290-0729
DOI:10.1016/j.ijthermalsci.2013.04.017