Wave breaking and wave runup/rundown have a major influence on nearshore hydrodynamics, morphodynamics and beach evolution. In the case of wave breaking, there is significant mixing of air and water at the wave crest, along with relatively high kinetic energy, so prediction of the free surface is complicated. Most hydrodynamic studies of surf and swash zone are derived from single-phase flow, in which the role of air is ignored. Two-phase flow modeling, consisting of both phases of water and air, may be a good alternative numerical modeling approach for simulating nearshore hydrodynamics and, consequently, sediment transport. A two-phase flow tool can compute more realistically the shape of the free surface, while the effects of air are accounted for. This paper used models based on two-dimensional, two-phase Reynolds-Averaged Navier-Stokes equations, the Volume-Of-Fluid surface capturing technique and different turbulence closure models, i.e., k-ε, k-ω and Re-Normalized Group (RNG). Our numerical results were compared with the available experimental data. Comparison of the employed method with a model not utilizing a two-phase flow modeling demonstrates that including the air phase leads to improvement in simulation of wave characteristics, especially in the vicinity of the breaking point. The numerical results revealed that the RNG turbulence model yielded better predictions of nearshore zone hydrodynamics, although the k-ε model also gave satisfactory predictions. The model provides new insights for the wave, turbulence and means flow structure in the surf and swash zones.