I'm running and modifying this example http://sfepy.org/doc-devel/examples/navier_stokes/navier_stokes.html. I can't see any specific mention of boundary conditions at the outflow. It it using Neumann boundary conditions with sigma_nn=sigma_nt = 0?