The development of grain structures during directional solidification of multi-crystalline silicon (mc-Si) plays a crucial role in the materials quality for silicon solar cells. Three dimensional (3D) modelling of the grain boundary (GB) interaction and evolution based on phase fields by considering anisotropic GB energy and mobility for me-Si is carried out for the first time to elucidate the process. The energy and mobility of GBs are allowed to depend on misorientation and the GB plane. To examine the correctness of our method, the known the coincident site lattice (CSL) combinations such as (Sigma a+Sigma b -> Sigma ax b) or (Sigma a+Sigma b -> Sigma a/b) are verified. We frther discuss how to use the GB normal to characterize a Sigma 3 twin GB into a tilt or a twist one, and show the interaction between tilt and twist Sigma 3 twin GBs. Two experimental scenarios are considered for comparison and the results are in good agreement with the experiments as well as the theoretical predictions. (C) 2018 Elsevier B.V. All rights reserved.