Numerical Simulation of Multi-Crystalline Silicon Crystal Growth Using a Macro–Micro Coupled Method during the Directional Solidification Process

In this work, the crystal growth of multi-crystalline silicon (mc-Si) during the directional solidification process was studied using the cellular automaton method. The boundary heat transfer coefficient was adjusted to get a suitable temperature field and a high-quality mc-Si ingot. Under the conditions of top adiabatic and bottom constant heat flux, the shape of the crystal-melt interface changes from concave to convex with the decrease of the heat transfer coefficient on the side boundaries. In addition, the nuclei form at the bottom boundary while columnar crystals develop into silicon melt with amzigzag-faceted interface. The higher-energy silicon grains were merged into lower energy ones. In the end, the number of silicon grains decreases with the increase of crystal length.