Cavitation phase change phenomenon appears in many engineering applications, often eroding and damaging surfaces, so deteriorating the performance of devices. Therefore, it is a phenomenon of great interest for the research and industry communities. In this work, three different cavitation models, the Homogeneous Relaxation Model (HRM), the Schnerr and Sauer, and the Kunz, are implemented in a Eulerian multiphase homogeneous flow Computational Fluid Dynamics (CFD) solver previously developed for simulating fully atomized sprays. The improved solver can be used then to study not only cases with cavitation, such a hydrofoil, but also situations where cavitation occurs together with liquid atomization, such as high pressure injection systems. Validation of this solver is carried out for three different cases under diverse operating conditions: a two-dimmensional throttle, a hydrofoil and a single-hole fuel injector. The Reynolds-Averaged Navier-Stokes (RANS) approach isThe equipment and resources used in this work have been partially supported by the Universitat Politècnica de València in the framework of the PAID-06-18 program (reference SP20180170). Additionally, the Ph.D. student María Martínez has been funded by a grant from the Government of Generalitat Valenciana with reference ACIF/2018/118 and financial support from The European Union.