TY - EJOU
AU - Kohno, Haruhiko
AU - Tanahashi, Takahiko
TI - Three-dimensional Numerical Simulation of Unsteady Marangoni Convection in the CZ Method using GSMAC-FEM
T2 - Computer Modeling in Engineering \& Sciences
PY - 2001
VL - 2
IS - 2
SN - 1526-1506
AB - Three-dimensional (3D) unsteady numerical simulations are carried out by means of the finite element method (FEM) with the generalized simplified marker and cell (GSMAC) method in silicon melt with a non-deformable free surface with Prandtl number Pr = 1.8534 × 10^{-2}, Marangoni number Ma = 0.0 - 6.2067 × 10^{2}, Grashof number Gr = 7.1104 × 10^{6}, and the aspect ratio As = 1.0 in the Czochralski (CZ) method. The flow state becomes unstable earlier by increasing the absolute value of the thermal coefficient of surface tension in the range of σ_{T} =0.0 - 1.5 × 10^{-5}N/mK. Although the velocity distribution in the circumferential direction is isotropy in any direction first, its magnitude becomes periodic and has the wavelength equal to 1/8 of the circumference. Then the wavelength doubles, and the flow pattern becomes finally asymmetrical. Moreover, the oscillation of the velocity distribution is observed just under the single crystal, and the amplitude is found to depend on the value of *σ*_{T} . After imposing the vertical magnetic field more than 0.05T to the melt from 50s, the flow pattern becomes restored to symmetry. But the instability remains under the single crystal and it indicates that the influence of Marangoni convection can not be neglected in the crystal growing process.
KW - Marangoni convection
KW - Czochralski method vertical magnetic field
KW - GSMAC-FEM
DO - 10.3970/cmes.2001.002.155