By extending LVOF to a two-phase model, we study impulsive waves created by landslides, based on our design of the mass-force coupling model (Li et al., 2007a; 2008). LVOF is a highly efficient Navier-Stokes solver, constructed by a novel VOF finite volume cut-cell approach that incorporates surface tension plus a dynamic subgrid-scale model (Li et al., 2004; 2007). According to the fictitious domain method, the mass-force model represents the coupling of a moving body on the flow, by handling the rigid body distributed over the particle domain as a fluid. Additionally, a ghost cell method (GFM) is used to treat the boundary conditions near the particle-fluid interface. Grid refinement studies are performed for test problems involving aerial landslides. Very encouragingly, the results agree well with measurements available. It is demonstrated that our results can provide valuable base for the analysis of near-field characteristics of landslide generated impulsive waves, by LVOF.