applications/utilities/postProcessing/velocityField/Mach/thermodynamicMach.H

   1 #   include "readThermodynamicProperties.H"
   2
   3     for (label i=startTime; i<endTime; i++)
   4     {
   5         runTime.setTime(Times[i], i);
   6
   7         Info<< "Time = " << runTime.timeName() << endl;
   8
   9         IOobject Uheader
  10         (
  11             "U",
  12             runTime.timeName(),
  13             mesh,
  14             IOobject::MUST_READ
  15         );
  16
  17         IOobject Theader
  18         (
  19             "T",
  20             runTime.timeName(),
  21             mesh,
  22             IOobject::MUST_READ
  23         );
  24
  25
  26         // Check U exists
  27         if (Uheader.headerOk() && Theader.headerOk())
  28         {
  29             mesh.readUpdate();
  30
  31             volVectorField U(Uheader, mesh);
  32             volScalarField T(Theader, mesh);
  33
  34             volScalarField Ma
  35             (
  36                 IOobject
  37                 (
  38                     "Ma",
  39                     runTime.timeName(),
  40                     mesh
  41                 ),
  42                 mag(U)/(sqrt(((Cv + R)/Cv)*R*T))
  43             );
  44
  45             Info << "Max Mach number: " << max(Ma.internalField()) << endl;
  46             Ma.write();
  47         }
  48         else
  49         {
  50             Info<< "    No U or T" << endl;
  51         }
  52     }