Jet impingement heat transfer is well known to enhance local heat transfer because of the local increase in the heat transfer coefficient and Nusselt number. When adding swirl, the peak heat transfer is shifted away from the stagnation point, as influenced by the jet offset distance, jet velocity or Reynolds number, and swirl number. When allowing phase change, additional enhancement is obtained through latent heat transfer. For this research, validation was conducted for single-phase axial and swirling jets, and for boiling axially impinging jets. Jet impingement boiling was modelled using CFD and the Rensselaer Polytechnic Institute (RPI) boiling model in Ansys Fluent for different jet offset distances and Reynolds numbers. Swirl caused the spreading out of the liquid phase resulting in vaporisation at the stagnation point; and the highly transient vapour formation away from the stagnation point was more spread out due to the swirling velocity profile. With the addition of swirl, the boiling curve moved to the left for low Reynolds numbers. The onset of nucleate boiling also occurred at a lower heating rate when swirl was added. A 3-D model incorporating swirl predicted a boiling curve to the left of the corresponding 2-D axisymmetric model.