Lithium-ion batteries with high energy density, long cycle life and other advantages, have been widely used to energy storage systems(ESS). But as ESS fires frequently occur, the safety concern has become the main obstacle that hinders the large-scale applications of lithium-ion batteries. Especially, thermal runaway is the key scientific problem in battery safety research. Therefore, in this study, we performed a numerical analysis on the thermal runaway phenomenon of NCM111, NCM523 and NCM622 batteries using a two-dimensional analysis model. The results show that the two-dimensional simulation results are generally matched with three-dimensional simulation. Also, In the case of NCM111 with a low Ni content in the temperature range used in this study, thermal runaway phenomenon does occurred very slowly, but as the Ni content is increased, the thermal runaway phenomenon occurs rapidly and the thermal stability tends to be decreased. And, in NCM523 and NCM622 batteries, chain reactions occur almost simultaneously, but in the case of NCM111 battery, it is found that after the SEI(Solid Electrolyte Interface) layer decomposition reaction, the cathode-electrolyte reaction is appeared sequentially. After that, the anodic decomposition reaction is increased and leads to the thermal runaway reaction.