This study performed the numerical analysis of the internal nozzle flows including cavitation phenomena by using the automated body-fitted grid generator and the multi-fluid model. The effect of grid refinement and the validation of multifluid model were investigated using four computational meshes under two test conditions. The mesh #3 was chosen as the optimum which can reduce the computational time and have good prediction ability to identify the cavitation region simultaneously. In addition, the computed results using multi-fluid model were compared with the reference experimental observations and numerical simulation results using homogeneous equilibrium model. From the distribution of volume fraction and velocity field, the multi-fluid model predicted the internal nozzle flows well when the liquid quality parameters were selected as $1.0{ imes}10^{12}$ for initial number density and 25 ${mu}m$ for bubble diameter.