Abstract: The moving-particle semi-implicit method (MPS), coupledwith a large eddy simulation (LES), is presented to free surface problem forturbulent flow. The control equation for LES is derived though the filteringoperation of the Navier-Stokers equation. In the control equation,additional Reynolds stress terms are only different between LES and MPS.Therefore MPS is extended to LES by coupled with a sub-particle-scale (SPS)turbulence model and the Smagorinsky model which is able to describe theReynolds stress terms. The MPS-LES method is used in the case of simulationof resonance sloshing, and the simulation result is very close to that ofexperiment or other numerical simulation.