In the present study,2-D large eddy simulations (LES) are conducted for flow past a porous circular army with a solid volume fraction (SVF) of 8.8％,15.4％ and 21.5％.Such simulations are relevant to understanding flow in natural streams and channels containing patches of emerged vegetation.In the simulations discussed in the paper,the porous cylinder of diameter D contains a variable number of identical solid circular cylinders (rigid plant stems) of diameter d =0.048D.Most of the simulations are conducted at a Reynolds number of 2 100 based on the diameter D and the velocity of the steady uniform incoming flow.Though in all cases wake billows are shed in the regions where the separated shear layers (SSLs) forming on the sides of the porous cylinder interact,the effect of these wake billows on the mean drag is different.While in the high SVF case (21.5％),the total drag force oscillates quasi-regularly in time,similar to the canonical case of a large solid cylinder,in the cases with a lower SVF the shedding of the wake billows takes place sufficiently far from the cylinder such that the unsteady component of the total drag force is negligible.The mean amplitude of the oscillations of the drag force on the individual cylinders is the largest in a streamwise band centered around the center of the porous cylinder,where the wake to wake interactions are the strongest.In all cases the maximum drag force on the individual cylinders is the largest for the cylinders directly exposed to the flow,but this force is always smaller than the one induced on a small isolated cylinder and the average magnitude of the force on the cylinders directly exposed to the flow decreases monotonically with the increase in the SVF.Predictions of the global drag coefficients,Strouhal numbers associated with the wake vortex shedding and individual forces on the cylinders in the array from the present LES are in very good agreement with those of 2-D direct numerical simulations conducted on finer meshes,which suggests LES is a better option to numerically investigate flow in channels containing canopy patches,given that LES is computationally much less expensive than DNS at high Reynolds number.To prove this point,the paper also discusses results of 2-D LES conducted at a much higher Reynolds number,where the near-wake flow is strongly turbulent.For the higher Reynolds number cases,where the influence of the turbulence model is important,the effect of the sub-grid scale model and the predictive capabilities of the unsteady Reynolds averaged Navier-Stokes (RANS) approach to predict flow past porous cylinders are discussed.