In this article, the three-dimensional unsteady multiphase flow is simulated in the whole passage of Francis hydraulic turbine. The pressure pulsation is predicted and compared with experimental data at positions in the draft tube, in front of runner, guide vanes and at the inlet of the spiral case. The relationship between pressure pulsation in the whole passage and air admission is analyzed. The computational results show: air admission from spindle hole decreases the pressure difference in the horizontal section of draft tube, which in turn decreases the amplitude of low-frequency pressure pulsation in the draft tube; the rotor-stator interaction between the air inlet and the runner increases the blade-frequency pressure pulsation in front of the runner.