Large eddy simulations on film cooling hole array on a flat plate was carried out to investigate upstream turbulence effect. Circular cylinders were configured to create a turbulent boundary layer and its diameter has been adjusted to generate 13% upstream turbulence intensity in the main flow. Due to the small pitch to diameter configuration of the cylinder, two-dimensional LES analysis was carried out in advance and the results showed that LES was an essential method to resolve flow field around and downstream circular cylinder, which was not available in RANS simulations. The three-dimensional LES results showed reasonable agreement in turbulence intensity and normalized velocity distributions along the vertical with measured data. According to the blowing ratio, the cooling flow coverage on the surface along the stream-wise direction was varied and well agreed with measured data. Additionally, upstream boundary flows were partially ingested inside the cooling hole and discharged again near along the centerline of the cooling hole. This accounted for film cooling effectiveness distribution inside the cooling hole surface and along the centerline. The current study revealed that the LES for predicting turbulent boundary layer behaviors due to upstream turbulence generation source was an effective and feasible method. Moreover, the LES effectively resolved flow fields such as film cooling flow behaviors and corresponding film cooling effectiveness distributions.