In the present study, sand behavior on the transformation of sandy beaches is simulated by a Lagrangian approach using a solver based on smooth particle hydrodynamics (SPH) method. There are interpolation points that are assumed as fluid particles with arbitrary distribution. These particles carry mass, velocity, density, and other material properties depending on the given problem. The multi-phase flows can be modeled by SPH and each particle is assigned to a different phase. The WCSPH method is used to solve Newtonian and non-Newtonian fluid flow which causes smaller computation power in comparison with the fluid dynamics. The present method is validated by solving a dam break, submerged hatch and sedimentation problems. In addition, the computational results are evaluated with the experimental and numerical data. The formation of bed ripples in sandy beaches is investigated in details and an intended model is accomplished.