{ Info<< "k(" << runTime.timeName() << ") = " << 3.0/2.0*average(magSqr(U)).value() << endl; Info<< "epsilon(" << runTime.timeName() << ") = " << (0.5*nu*average(magSqr(twoSymm(fvc::grad(U))))).value() << endl; Info<< "U.f(" << runTime.timeName() << ") = " << 181.0*average(U & force).value() << endl; }