Some approximation procedures are presented for the system of equations arising from the large eddy simulation of turbulent flows. Existence of solutions to the approximate problems is proved. Discrete solutions generate a strongly convergent subsequence whose limit is a weak solution of the original problem. To prove the convergence theorem we use Young measures and related tools. We do not limit ourselves to divergence-free functions and our results are in particular valid for finite element approximations where one usually does not use functions with divergence equal to zero.