Closures for large-eddy simulations of passive scalars