Immersed boundary method and wall models are adopted to simulate the laminar and turbulent flows over periodic hills with large eddy simulation. The Reynolds numbers based on the hill hight are 25∼100 for the laminar flow, whereas are 2800, 5600 and 10595 for the turbulent flows. The results show that the flow over periodic hills consists of separation, recirculation and reattachment for the laminar and turbulent flows. The locations of separation and reattachment points are linear with the Reynolds number for laminar flows, but exhibit slightly difference with the Reynolds number for turbulent flows. For the simulation of turbulent flow over periodic hills, the results reveal that the present numerical treatment with grid resolutions (200 × 95 × 96) are sufficient for capturing the turbulence quantities compared with the DNS data of Breuer et al. [1] at lower Reynolds number Re = 2800. However, due to the inappropriate implement of wall model, the near wall turbulence is fail to be predicted as Reynolds number increased. It indicates that the adopted wall model with velocity correction could not accurately capture the near wall Reynolds stresses on the curvature wall boundary.