Introduction. The reduction in beach width due to erosion is a significant issue that can either be mitigated or exacerbated by coastal protection structures. Modelling breaking waves near the coast and around coastal structures can be used to determine their impact on the dynamics of the coastal zone. The objective of this study is to model and analyze the dynamics of turbulent structures around a single breakwater, obtained using two turbulence modelling schemes: RANS and LES.Materials and Methods. Turbulence induced by breaking waves was investigated. The modelling was based on bathymetric measurements conducted along the Azov Sea coast and a three-dimensional wave hydrodynamics model supplemented with various turbulence calculation configurations.Results. Modelling results of wave processes generating turbulent flows in the presence of coastal protection structures using different turbulence models were obtained. Results obtained based on Reynolds-averaged Navier-Stokes (RANS) equations are compared with the results of Large Eddy Simulation (LES) approach with Smagorinsky dynamic subgrid-scale model (DSM).Discussion and Conclusions. The results showed that wave heights simulated by LES were higher than those simulated by RANS in the front and leeward regions of the coastal protection structure and were lower in its upper part. Thus, according to LES, a greater amount of wave energy was preserved after passing over the breakwater. Velocity vectors of the water medium showed the formation of a vortex when LES was used, whereas no evidence of such turbulent vortices was detected in the case of RANS, confirming the better performance of LES for turbulence modelling in the coastal zone. According to the presented results, LES is the preferred tool for generating turbulence under incoming wave conditions in engineering practices.