Info<< "k("
             << runTime.timeName()
             << ") = "
             << 3.0/2.0*average(magSqr(U)).value() << endl;

        Info<< "epsilon("
             << runTime.timeName()
             << ") = "
             << (
                     0.5*nu*average
                     (
                         magSqr(fvc::grad(U) + fvc::grad(U)().T())
                     )
                ).value() << endl;

        Info<< "U.f("
             << runTime.timeName()
             << ") = "
             << 181.0*average(U & force).value() << endl;