An axisymmetric three-dimensional model for temperature distribution in the electrical discharge machining process has been developed using the finite element method to estimate the surface integrity characteristics of AISI H13 tool steel as workpiece. White layer thickness, depth of heat affected zone, and arithmetical mean roughness consisting of the studied surface integrity features on which the effect of process parameters, including pulse on-time and pulse current were investigated. Additionally, the experiments were carried out under the designed full factorial procedure to validate the numerical results. Both numerical and experimental results show that increasing the pulse on-time leads to a higher white layer thickness, depth of heat affected zone, and the surface roughness. On the other hand, an increase in the pulse current results in a slight decrease of the white layer thickness and depth of heat affected zone, but a coarser surface roughness. Generally, there is a good agreement between the experimental and the numerical results.