In previous papers, we proposed finite element schemes based on the Petrov-
Galerkin weak formulation using exponential weighting functions for solving
accurately, and in a stable manner, the flow field of an incompressible vi
scous fluid. In this paper, we present the Petrov-Galerkin finite element s
cheme for turbulent flow fields based on large eddy simulation using the st
andard Smagorinsky model with the Van Driest damping function. The filtered
incompressible Navier-Stokes equations are numerically integrated in time
by using a fractional step strategy with second-order accurate Adams-Bashfo
rth explicit differencing for both convection an diffusion terms. In order
to evaluate more accurately a mass matrix, the well-known multi-pass algori
thm was also adopted in this study. Numerical results obtained are compared
through flow around a rectangular cylinder at Re = 22,000 with the experim
ental data and other existing numerical data.