Provided are a sub-particle scale turbulence simulation method for SPH fluid, and a system and recording medium for the method. In the present disclosure, by combining a multi Eulerian grid with a SPH system, various Eulerian grids are combined with the SPH system while the vorticity of particles is efficiently calculated, which allows firm detection of a deformation region. For this reason, along with the flexibility and simplicity of the multiple grids system, the present disclosure may be easily expanded to a broad spectrum in aspects of time and space. Moreover, the present disclosure may express multi-level vorticity, which could not be expressed by an existing SPH system, and give a stable and visually satisfactory result.