Glad I could help. Just for the sake of closure, we used a parallel compliant preconditioner and made tighter tolerance on the algebraic solvers. Am I correct to say these fixed your issues with 2nd order interpolation of turbulence? Also, what does your comparison with the other variations suggest, are they comparable? Or is there a definite winner?