Edited UEqn to calculate the stress term as  fvm::laplacian(mu, U) + fvc::div(mu*fvc::grad(U)().T())