Coal and gas outburst is a common hazard in coalmines and due to the complexity of its nature, efficaciously predicting the possibility of its occurrence is intricate. The mechanical and gas flow properties of coal, in situ and mining-induced stresses and geological structure are involved in outburst occurrence, whereas suggested thresholds for its prediction do not consider them all. Numerical modelling is known as an effective method that can not only consider couple the effect of different physics but also analyse the effect of associated parameters. In this paper, a coupled fluid flow and geomechanics model has been developed in COMSOL Multiphysics to simulate outbursts by considering the effect of gas pressure and composition, damage mechanics on related coal properties, and progressive excavation and mining-induced stress. For the proposed fully coupled outburst model, the effect of pore pressure and sorption are applied to the stress while the influence of stress changes and damage are considered in updating coal porosity and permeability in each time step. In addition, to study the effect of mining-induced stress, the excavation process is simulated by changing the material properties over time according to the schedule of heading development. Then, outburst initiation is analysed from the energy point of view by using the concept of energy conservation law. The simulation results are used in calculating energy sources of the outburst initiation, also, the effect of particle size and ejection velocity are studied, and the significance of particle size effect on the energy analysis and deformed coal on outburst initiation is discussed.