applications/solvers/DNS/dnsFoam/globalProperties.H

   1         Info<< "k("
   2              << runTime.timeName()
   3              << ") = "
   4              << 3.0/2.0*average(magSqr(U)).value() << endl;
   5
   6         Info<< "epsilon("
   7              << runTime.timeName()
   8              << ") = "
   9              << (
  10                      0.5*nu*average
  11                      (
  12                          magSqr(fvc::grad(U) + fvc::grad(U)().T())
  13                      )
  14                 ).value() << endl;
  15
  16         Info<< "U.f("
  17              << runTime.timeName()
  18              << ") = "
  19              << 181.0*average(U & force).value() << endl;
  20