Energy dissipation mechamism is the key to study tsunami hazard mitigation. Numerical method is adopted to study the interaction between bores and square cylinders. The model solves the three-dimensional Navier-Stokes equations with Large-Eddy Simulation turbulence model. The Volume-of-fluid (VOF) method is used to track the complex free surface. We focus the investigation on the effect of cylinder height on the flow field. The results show that the turbulence diffusion is the main mechanism for energy dissipation. The flow patterns are significantly different within and beyond the cylinder array. The taller cylinders cause smaller velocity magnitude in the downstream area. In addition, a larger value of velocity magnitude and vorticity near the bottom is identified in the tall-cylinder case. These unique featuers make different dissipation rates.