In order to develop a numerical code, which can simulate a high-speed cryogenic cavitation flow in a rotating machine like a turbo-pump, a new cavitation model is considered. This model can handle whichever a pure liquid flow, a pure vapor flow, a mixture flow including bubbly flow by the bubble size distribution model. In addition, switching the flow patterns is seamless. This model is applied into a new numerical code accompanied by primitive variables, preconditioning, dual-time stepping, SD-SLAU and MFGS.