This paper presents pyrolysis simulations of a multi-component structure model of Fugu sub-bituminous coal containing 23,898 atoms with ReaxFF MD. The slow heat-up and long-time isothermal simulations were performed using the GPU-enabled ReaxFF MD code of GMD-Reax and reaction analysis code of VARxMD. The simulated pyrolyzate evolution trends and pyrolysis stages proposed based on lumped pyrolyzates agree fairly with literature. Major reaction pathways for early generation of CO2, H2O, CH4 and H-2 in the activation and primary pyrolysis stages are found closely associated with carboxyl and methoxyl groups, indicating the significant role of oxygen-containing groups in the initialization of Fugu coal pyrolysis. Compared with reported experiments, the coal tar yield is over-predicted, while the gas yield is slightly under-estimated. To shorten the simulation time from 2 ns to 250 ps, an increase of 400 K in average on the basis of temperature range of 1200-2200 K is needed that will result in the over-prediction for the yields of gas and tar, as well as under prediction for the non-volatile yield. This work suggests that ReaxFF MD simulation is useful for revealing the overall scenario and detailed reactions in coal pyrolysis but cautions should be taken for proper elevated temperatures.