The flow through porous media was numerically simulated using lattice Boltzmann method (LBM). By predicting permeability

the capability and accuracy for simulating flow through porous media were evaluated for various of lattice models: lattice Bhatnagar-Gross-Krook (LBGK) or called single relaxation time (SRT)

multiple relaxation time (MRT) and entropic lattice Boltzmann method (ELBM). In order to research the influences of each free parameter to permeability in the MRT model

a total of 12 kinds of combinations were selected. Furthermore

the large eddy simulation (LES) and LBM were coupled to simulate flow through porous media at high Reynolds number and the transformation of flow patterns. The results indicate that permeability increased with viscosity predicted by LBGK and ELBM

but there are little variation for permeability predicted by MRT. In addition

the different combinations of relaxation parameters have a greater impact to the permeability

and an optimum combination is recommended for simulating flow through porous media. Simulated results are in quantitatively agreement with the empirical formula. The transformations of flow patterns are successfully predicted in porous media using LES-MRT. Higher Reynolds number

more and bigger vortexes in porous media are observed.