Set mean velocity (flux) instead of mean pressure gradient. Got a turbulent flow solution with Re=10000. Note that with this change, Re is the bulk Re not Re_tau.