[didomenico]

Hi everybody,
I am computing a jet in still air and I want to test different discretization schemes for the convection term using simpleFoam.
I set up the geometry (axisymmetric), got the solution with default parameters (first order upwind) and then I wanted to switch to QUICK scheme (to compare the solution with Fluent). Therefore I set

div(phi,U) Gauss QUICK;

w/o changing any other parameter. After few iterations (50) I get suddenly unphysical velocities (1e3 m/s) near the outflow. I am using the standart kEpsilon model with wall function. Am I missing something?
Thanks
Massimiliano