Simulation of the temperature distribution and evolution during pulse electrochemical machining can be a computationally very expensive procedure. In a previous part of the work [Smets et al. J Appl Electrochem 37(11):1345, 2007] a new approach to calculate the temperature evolution was introduced: the hybrid method, which combines averaged and pulsed calculations. The averaged calculations are performed by time averaging the boundary conditions and the bulk heat sources of the system. The timesteps used during the averaged calculations are then no longer dictated by the pulse characteristics. Using this approach, computationally very cheap, yet satisfactory results can be obtained. The analysis in the previous part of the work was obtained from analytical solutions on simplified models. In this part, the more general case is solved numerically. Multiple geometries are simulated and analyzed and methods are compared. Very satisfactory, yet cheap results are obtained.