Role of marangoni tension effects on the melt convection in directional solidification process for multi-crystalline silicon ingots

We carried out global simulations to investigate the marangoni tension effect on the thermal and flow fields in the silicon melt of the directional solidification process for multi-crystalline silicon ingots. The argon flow rate was varied to provide different solidification conditions and to change the relative values between the argon shear stress and the marangoni tension at the melt free surface. We found that the marangoni tension together with the shear stress mainly influences the upper layer melt convection while the thermal buoyancy force dominates the bulk flow of the melt. At low argon flow rates, the argon shear stress can be neglected and the marangoni tension alone enhances the melt convection intensity near the gas–melt–crucible triple junction point. The marangoni tension is so weak that it cannot modify the melt flow pattern in this case. For medium flow rate, the marangoni tension can significantly weaken the shear stress effect at the outer part of the melt free surface, leading to a distinctive flow pattern in the silicon melt. With further increase in argon flow rate, the shear stress sharply increases and dominates the upper layer melt flow, limiting the marangoni tension effect to the triple point. The numerical results are helpful for better understanding and controlling of the directional solidification process for high quality multi-crystalline silicon ingots.