In this paper, lattice Boltzmann method and large eddy simulation are adopted to simulate the laminar and turbulent Poiseuille flow through square duct. The shear or friction Reynolds number based on the duct width are 10 for the laminar flow, whereas are 360 for the turbulent flow. We validate the governing equations and the boundary conditions by simulating laminar Poiseuille duct flow. The results of the laminar flow show good agreement compared with analytic solution. For the turbulent Poiseuille square duct flow at Re = 360, our simulation is able to capture the turbulence quantities by observing the mean streamwise velocity profile along the wall bisector compared with the direct numerical simulation data and the law of the wall. Then the simulation result of turbulence intensities and the Reynolds stress variation along the wall bisector is also shown in present research.