applications/utilities/postProcessing/miscellaneous/postChannel/calculateFields.H

   1         /*
   2         volTensorField gradU = fvc::grad(U);
   3         volSymmTensorField D = symm(fvc::grad(U));
   4         volTensorField Dprim = symm(fvc::grad(U - UMean));
   5
   6         volScalarField prod = -((U - UMean)*(U - UMean)) && D;
   7         */
   8
   9         /*
  10         volScalarField txx
  11         (
  12             IOobject
  13             (
  14                 "txx",
  15                 runTime.timeName(),
  16                 mesh,
  17                 IOobject::NO_READ,
  18                 IOobject::NO_WRITE
  19             ),
  20             mesh,
  21             dimensionSet(0, 1, -1, 0, 0)
  22         );
  23         txx =sqrt(Txx - (UMeanx*UMeanx));
  24         txx.write();
  25
  26         volScalarField tyy
  27         (
  28             IOobject
  29             (
  30                 "tyy",
  31                 runTime.timeName(),
  32                 mesh,
  33                 IOobject::NO_READ,
  34                 IOobject::NO_WRITE
  35             ),
  36             mesh,
  37             dimensionSet(0, 1, -1, 0, 0)
  38         );
  39         tyy = sqrt(Tyy - (UMeany*UMeany));
  40         tyy.write();
  41
  42         volScalarField tzz
  43         (
  44             IOobject
  45             (
  46                 "tzz",
  47                 runTime.timeName(),
  48                 mesh,
  49                 IOobject::NO_READ,
  50                 IOobject::NO_WRITE
  51             ),
  52             mesh,
  53             dimensionSet(0, 1, -1, 0, 0)
  54         );
  55         tzz = sqrt(Tzz - (UMeanz*UMeanz));
  56         tzz.write();
  57
  58         volScalarField txy
  59         (
  60             IOobject
  61             (
  62                 "txy",
  63                 runTime.timeName(),
  64                 mesh,
  65                 IOobject::NO_READ,
  66                 IOobject::NO_WRITE
  67             ),
  68             mesh,
  69             dimensionSet(0, 2, -2, 0, 0)
  70         );
  71         txy = Txy - (UMeanx*UMeany);
  72         txy.write();
  73 */