화학공학소재연구정보센터
International Journal of Heat and Mass Transfer, Vol.44, No.10, 1983-2003, 2001
Three-dimensional numerical simulation of Marangoni instabilities in non-cylindrical liquid bridges in microgravity
Instability of Marangoni convection in non-cylindrical (convex or concave) liquid bridges of low Prandtl number fluids is investigated by direct three-dimensional and time-dependent simulation of the problem. Body-fitted curvilinear coordinates are adopted: the non-cylindrical original physical domain in the (r, z, phi) space is transformed into a cylindrical computational domain in a (xi, eta, phi) space, The geometry of the domain is transformed using a coordinate transformation method by surface fitting technique. The field equations are numerically solved explicitly in time and with a finite difference technique in a staggered grid. The numerical results are analyzed and interpleted in the general context of the bifurcation's theory. The computations show that for semiconductor melts the first bifurcation is characterized by the loss of spatial symmetry rather than by the onset of oscillatory flow and that it is hydrodynamic in nature. The flow field azimuthal organization related to the critical wave number, depends on the geometrical aspect ratio A = LID of the liquid bridge and on the shape factor S (convex S > 1, concave S < 1) of the free surface. The critical azimuthal wave number increases when the geometrical aspect ratio of the bridge is decreased and, for a fixed aspect ratio, can be shifted to higher values by increasing the volume (convex bridges) or to lower values by decreasing the volume (concave bridges). This behavior is explained on the basis of the relation between the topology of the azimuthal disturbances and the structure of the fluid-dynamic held. A generalized law is found to correlate the critical azimuthal wave number of the instability to the geometrical aspect ratio and to the shape factor. A second oscillatory (Hopf) bifurcation occurs M hen further increasing the Marangoni number. Experimental results available in literature on this second bifurcation are considered for comparison. The experimental and numerical results show a good agreement. <(c)> 2001 Elsevier Science Ltd. All rights reserved.