This is what I did basically and the values you see there is the exact one I used. And the mesh has non-orhtogonality < 60 and skewness < 4. I'm in a good range imo. The only things where it could have an issue, like I mentionned, are the schemes, BC or IC. But if I stick to what you said, even if the residuals don't reach the target, it make sense that the results at that point can be a relatively good IC in transient mode.