In this paper we use of large eddy simulation method for 2 dimensional incompressible fluid in square lid-driven cavity. We also apply smagorinsky SGS model for attach sub grid scales. In order to we have written a computer program with finite volume method and implicit for discretization unsteady governing equations and we use of simple algorithm for pressure term. Our computation get from Re=100 up to Re=10000 and we obtain mean horizontal velocity, mean vertical velocity in central line, stream line and pressure contour and also we have calculated fluctuation of velocities and in finally we accepted Reynolds stress in Re=10000. The results are in good agreement with the experimental and numerical data.