Numerical simulation for the air entrainment of aerated flow with an improved multiphase SPH model