The behavior of crack initiation, propagation and coalescence in rock is the key mechanics of rock deformation, strength and stability, and it is important to predict and analyze crack growth behavior by numerical method. Phase field(PF) method overcomes the challenge of the traditional fracture mechanics theory and numerical method to explicitly track the crack surface within the framework of strict and accurate theory, and can simulate crack initiation, propagation, branching and multi-cracking. A multi-scale model that can reflect the internal microstructure of rock materials was established at the main area of crack growth. The failure progress of crack initiation, propagation and coalescence within the basalt under three-point bending test was predicted and analyzed by phase field method. The influence of the distribution of mineral inclusions and energy release rate on the crack propagation behavior was investigated and the numerical simulation results were compared and verified with the experimental data. It is shown that the proposed multi-scale PF model considering the complex internal microstructure of basalt rock is able to simulate and predict the progress of the crack behavior. The behavior of crack initiation, propagation and coalescence in rock is the key mechanics of rock deformation, strength and stability, and it is important to predict and analyze crack growth behavior by numerical method. Phase field(PF) method overcomes the challenge of the traditional fracture mechanics theory and numerical method to explicitly track the crack surface within the framework of strict and accurate theory, and can simulate crack initiation, propagation, branching and multi-cracking. A multi-scale model that can reflect the internal microstructure of rock materials was established at the main area of crack growth. The failure progress of crack initiation, propagation and coalescence within the basalt under three-point bending test was predicted and analyzed by phase field method. The influence of the distribution of mineral inclusions and energy release rate on the crack propagation behavior was investigated and the numerical simulation results were compared and verified with the experimental data. It is shown that the proposed multi-scale PF model considering the complex internal microstructure of basalt rock is able to simulate and predict the progress of the crack behavior. [ABSTRACT FROM AUTHOR]