Abstract This paper uses two-dimensional boundary element method (BEM) numerical modeling to analyze the deformation and failure behavior of a coal seam and to understand the nature of gas flow into a roadway entering the Barapukuria coal mine in Bangladesh. The Barapukuria basin contains Permian-aged Gondwana coals with high volatile B bituminous rank. Three models (A, B, and C) are presented here. Model A assumes horseshoe-shaped geometry, model B assumes trapezoid-shaped geometry, and model C assumes horseshoe-shaped geometry coupled with a roof fall-induced cave generated by the break-up of rock materials along the vertical dimension of an igneous dyke. The simulation results show that there is little difference in strata deformation between models A and B. In model A, there is no horizontal tensional stress and the overall horizontal stress patterns are compressive, while the distribution and magnitude of vertical stress show higher tensional stresses on the immediate rib sides and floor. In model B, both horizontal and vertical stress distributions indicate low to medium tensional stresses on the immediate roof, floor, and rib sides, but compressive stresses are prominent toward the interior of the coal seam. Deformation vectors indicate that failure extends laterally to about 7.5 m around the excavation geometry. On the contrary, for model C, the distributions and magnitudes of horizontal and vertical stress show higher tensional stresses in both rib sides of the roof fall zone. The deformation around the dyke-induced perturbation zone affects a large volume of coal. The deformation vectors with high magnitudes are nearly horizontal and propagate laterally up to 30 m; whereas, low-magnitude deformation vectors extend about 25 m toward the roof and 20 m toward the floor. The vertical tensional displacement, which is concentrated in the floor and the left and right hand sides of the roof, propagates about 30 m on both sides and about 22 m in the floor. From these simulation results, it is thought that the extension of the dyke-induced perturbation zone toward the roof, floor, and rib sides of the entry roadway initially creates small tensional cracks that gradually grow into large-scale tensional features. These features could also be responsible for high concentrations of gas, which are emitted into the mine from fractured coals due to insufficient mine ventilation and low atmospheric pressure.