<p>In this work, the phase-field theory is applied to the thermodynamic equations to effectively solve the transient crack propagation problem of ceramics under flame thermal shock. The entire dynamic cracking process recorded in experiments is reproduced numerically. The inertial effect on the crack initiation and growth is considered in the whole simulation process. Comparing the numerical results with the experimental results in the literature, it is found that the present method can accurately estimate the instantaneous speed of crack propagation. In addition, the effects of mesh size and different geometrical parameters of the prefabricated crack such as its angles, positions, and lengths on the propagation speed of the moving crack are investigated through detailed numerical simulations. The numerical results show that compared to a quasi-static model, taking the inertia effect into account in simulations can provide better provision of the cracking process due to flame thermal shock.</p>