FreeFOAM The Cross-Platform CFD Toolkit
thermodynamicMach.H
Go to the documentation of this file.
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 
35  (
36  IOobject
37  (
38  "Ma",
39  runTime.timeName(),
40  mesh
41  ),
42  mag(U)/(sqrt(((Cv + R)/Cv)*R*T))
43  );
44  Ma.write();
45  }
46  else
47  {
48  Info<< " No U or T" << endl;
49  }
50  }
51 
52 // ************************ vim: set sw=4 sts=4 et: ************************ //