It is well known that the explicit time integration of heat transfer problem is conditionally stable. The very small time step leads to increase of computational time dramatically. In this work, the location of integration points in the mass matrix [1-5] using quadrilateral elements is discussed in detail in the simulation of transient heat transfer problems. From both theoretical and numerical perspectives, the stability of dynamic heat transfer model is improved significantly with adjustment of integration points in the mass matrix. The quantitative study has indicted that the numerical stability of discretised model is proportional to r value controlling the location of integration points in the mass matrix as r > 0 in the 2D problems.  Numerical experiments including heat conduction, convection and radiation with regular and irregular mesh have demonstrated the superior performance of the proposed integration points in terms of accuracy as well as stability. The successful development of robust, efficient and accurate explicit algorithms has opened a new window in the simulation of general transient heat transfer problems.