A computer code, ELANEX, including several Homogenous-Equilibrium-Model (HEM) type cavitation models, were developed, to numerically simulate natural cavitation phenomena. The effectiveness of the code was checked by cavitation flows around the disk and cylinder body for a wide range of different cavitation numbers. Cavity profiles were compared with the analytic solution of disk and empirical formulae fitted from the experiment data, and contrast between different cavitation models were fulfilled as well. The cavity length and maximal cavity diameter were found to agree well with the analytic solutions, and detailed cavity profiles were in accordance with the experimental formula. Comparison with the hemisphere headed cylinder body presented a good agreement of the pressure coefficient with the experiment data. Reasonable drag-force coefficient variation and drag-force reduction effect were obtained.