large eddy simulation cooperated with a physical fractional-step method was applied to simulate steady flow around two parallel circular cylinders. The total velocity vectors, pressure contours and vorticity magnitude are obtained. The modeling results conform to physical law, and show that the large eddy simulation theory has powerful capacity in simulation of microstructures of turbulent flows, and can be widely applied to the solution of real engineering problems.