Large eddy simulation of the passive scalar transport in a high Schmidt number turbulent confined wake flow has been performed. The results are evaluated by comparison to particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF) data, including point-wise data as well as spatial correlations. In the LES simulations, the gradient diffusion hypothesis is used to close the transport equation for the passive scalar. Different discretization schemes are investigated in order to determine the best choice for ensuring boundedness of the passive scalar and to accurately predict the mixing rate. The simulation results compare well to experimental data, demonstrating that the transport mechanisms in this high Schmidt number turbulent flow are well predicted by the LES method. Two-point spatial correlations of passive scalar with velocity predicted by the simulation show good agreement with the experimental results, indicating that the turbulent coherent structures of the flow are reproduced by the simulation. (C) 2015 Elsevier Ltd. All rights reserved.