This present paper proposes a two-dimensional lattice Boltzmann model coupled with a Large Eddy Simulation (LES) model and applies it to flows around a non-submerged groyne in a channel. The LES of shallow water equations is efficiently performed using the Lattice Boltzmann Method (LBM) and the turbulence can be taken into account in conjunction with the Smagorinsky Sub-Grid Stress (SGS) model. The bounce-back scheme of the non-equilibrium part of the distribution function is used to determine the unknown distribution functions at inflow boundary, the zero gradient of the distribution function is set normal to outflow boundary to obtain the unknown distribution functions here and the bounce-back scheme, which states that an incoming particle towards the boundary is bounced back into fluid, is applied to the solid wall to ensure non-slip boundary conditions. The initial flow field is defined firstly and then is used to calculate the local equilibrium distributions as initial conditions of the distribution functions. These coupled models successfully predict the flow characteristics, such as circulating flow, velocity and water depth distributions. The comparisons between the simulated results and the experimental data show that the model scheme has the capacity to solve the complex flows in shallow water with reasonable accuracy and reliability. This present paper proposes a two-dimensional lattice Boltzmann model coupled with a Large Eddy Simulation (LES) model and applies it to flows around a non-submerged groyne in a channel. The LES of shallow water equations is efficiently performed using the Lattice Boltzmann Method (LBM) and the turbulence can be taken into account in conjunction with the Smagorinsky Sub-Grid Stress (SGS) model. The bounce-back scheme of the non-equilibrium part of the distribution function is used to determine the unknown distribution functions at inflow boundary, the zero gradient of the distribution function is set normal to outflow boundary to obtain the unknown distribution functions here and the bounce-back scheme, which states that an incoming particle towards the boundary is bounced back into fluid, is applied to the solid wall to ensure non-slip boundary conditions. The initial flow field is defined first and then is used to calculate the local equilibrium distributions as initial conditions of the distribution functions. These coupled models successfully predict the flow characteristics, such as circulating flow, velocity and water depth distributions. The comparisons between the simulated results and the experimental data show that the model scheme has the capacity to solve the complex flows in Shallow water with reasonable accuracy and reliability.