Symmetry breaking of melt flow typically encountered in a Bridgman configuration heated from below

A computational model is developed to investigate three-dimensional fluid flow instability and transition to unsteadiness. The simulated domain corresponds to a parallelepiped configuration which is a restriction to the fluid phase of commonly used configurations in crystal growth (viz. Bridgman configuration). In the present work, this problem is studied using three-dimensional simulations and focuses on symmetry breaking and transition to unsteadiness occurring for a low Prandtl number fluid (Pr = 0.01). Three-dimensional simulations presented thus overcome the limitation of the two-dimensional approach. It has been found that the initially steady symmetric flow in the three-dimensional case becomes asymmetric for lower Ra number than for the two-dimensional case. The breaking in symmetry occurs firstly in the transverse plane. For the relatively low Ra number we still have no intensification in the global heat transfer but it appears that the heat transfer increases locally on the bottom and decreases on the vertical active walls. The classical spiral flows typical characterizing the three-dimensional effect (in the third direction) are also identified.

[1]  G. D. Davis,et al.  Effects of anisotropy and solid/liquid thermal conductivity ratio on flow instabilities during inverted Bridgman growth , 2004 .

[2]  Franz Rosenberger,et al.  Analysis of the low gravity tolerance of Bridgman-Stockbarger crystal growth II. Transient and periodic accelerations , 1991 .

[3]  B. P. Leonard,et al.  A stable and accurate convective modelling procedure based on quadratic upstream interpolation , 1990 .

[4]  A. A. Wheeler,et al.  Bifurcation analysis of solutal convection during directional solidification , 1991 .

[5]  Sam R. Coriell,et al.  Thermosolutal convection during directional solidification. II. Flow transitions , 1987 .

[6]  Kwang-Tzu Yang,et al.  Rayleigh-Bénard Convection in a Small Aspect Ratio Enclosure: Part I—Bifurcation to Oscillatory Convection , 1993 .

[7]  Three-dimensional bifurcations of a two-phase Rayleigh-Benard problem in a cylinder , 2001 .

[8]  J. C. Brice,et al.  The Growth of Crystals from Liquids , 1974 .

[9]  Jalil Ouazzani,et al.  Symmetry breaking flow transitions and oscillatory flows in a 2D directional solidification model , 1994 .

[10]  G. Guj,et al.  The Rayleigh–Bénard problem in intermediate bounded domains , 1993, Journal of Fluid Mechanics.

[11]  B. P. Leonard,et al.  Why You Should Not Use 'Hybrid', 'Power-Law' or Related Exponential Schemes for Convective Modelling—There Are Much Better Alternatives , 1995 .

[12]  R. G. Rehm,et al.  Thermosolutal convection during directional solidification , 1984 .

[13]  I. Sezai,et al.  Natural convection from a discrete heat source on the bottom of a horizontal enclosure , 2000 .

[14]  B. P. Leonard,et al.  Beyond first‐order upwinding: The ultra‐sharp alternative for non‐oscillatory steady‐state simulation of convection , 1990 .