SUMMARYA large eddy simulation based on ÿltered vorticity transport equation has been coupled with ÿltered probability density function transport equation for scalar ÿeld, to predict the velocity and passive scalar ÿelds. The ÿltered vorticity transport has been formulated using di usion-velocity method and then solved using the vortex method. The methodology has been tested on a spatially growing mixing layer using the two-dimensional vortex-in-cell method in conjunction with both Smagorinsky and dynamic eddy viscosity subgrid scale models for an anisotropic ow. The transport equation for ÿltered probability density function is solved using the Lagrangian Monte-Carlo method. The unresolved subgrid scale convective term in ÿltered density function transport is modelled using the gradient di usion model. The unresolved subgrid scale mixing term is modelled using the modiÿed Curl model. The e ects of subgrid scale models on the vorticity contours, mean streamwise velocity proÿles, root-mean-square velocity and vorticity uctuations proÿles and negative cross-stream correlations are discussed. Also the characteristics of the passive scalar, i.e. mean concentration proÿles, root-mean-square concentration uctuations proÿles and ÿltered probability density function are presented and compared with previous experimental and numerical works. The sensitivity of the results to the Schmidt number, constant in mixing frequency and in ow boundary conditions are discussed.