For a single particle phase the divDevTau term in fvSchemes is now specified
thus:
divDevTau(U.particles) Gauss linear;
previously:
div((((thermo:rho.particles*nut.particles)*dev2(T(grad(U.particles))))+(((thermo:rho.particles*lambda.particles)*div(phi.particles))*I))) Gauss linear;