Probability density function approach for multidimensional turbulentflow calculations with application to in-cylinder flows in reciprocating engines