The objective of the current paper is to study the propagating and breaking of solitary waves on stepped slopes beaches, to simulate the shoaling and breaking, specifically the location of breaking point Xb, and solitary wave height at breaking Hb of solitary waves on the different stepped slopes. Ansys Fluent is used to implement the simulation, a two-dimensional volume of fluid (VOF) which is based on the Reynolds-Averaged Navier–Stokes (RANS) equations and the k–ε turbulence closure solver. The obtained results were firstly validated with existing empirical formulas for solitary wave run-up on the slope without stepped structure and are compared with the experimental and numerical results. The numerical computation has been carried out for several, configurations of beach slopes with tan ß= 1:15, 1:20, 1:25, wave height H0= 0.04, 0.06, 0.08m, water depth h0= 0.15, 0.2, 0.25m, and step height Sh= 0.025, 0.05, 0.075m. A set of numerical simulations were implemented to analyze shoaling and breaking of solitary waves, wave reflection, wave transmission, and wave run-up with various parameters wave heights, water depth, beach slopes, and Sh step height.