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

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