Numerical studies play a major role in the understanding and prediction of plasticity induced crack closure (PICC). However, the available numerical models can be considered simplifications of reality as they consider discrete crack propagations, relatively high fatigue crack growth rates (FCGR), sharp cracks, and propagation occurring at well-defined loads. Besides, there are a great number of numerical and physical parameters affecting the predictions of PICC. The aim of this paper is to discuss the numerical study of PICC. The numerical parameters affecting the accuracy of the numerical simulations, and the dependent parameters used to characterise the plastic wake and the closure level, are identified. The influence of the radial size of crack front elements and crack propagation is analysed. An extrapolation model is proposed, with excellent results. An intrinsic uncertainty is associated with the number of load cycles between crack increments and the definition of crack closure level. Finally, the effect of the stress ratio (R) on crack closure level is analysed.