47 lines
1.1 KiB
C
47 lines
1.1 KiB
C
for (label i=startTime; i<endTime; i++)
|
|
{
|
|
runTime.setTime(Times[i], i);
|
|
|
|
Info<< "Time = " << runTime.timeName() << endl;
|
|
|
|
mesh.readUpdate();
|
|
|
|
IOobject Uheader
|
|
(
|
|
"U",
|
|
runTime.timeName(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
);
|
|
|
|
if (Uheader.headerOk())
|
|
{
|
|
volVectorField U(Uheader, mesh);
|
|
|
|
autoPtr<basicPsiThermo> thermo
|
|
(
|
|
basicPsiThermo::New(mesh)
|
|
);
|
|
|
|
volScalarField Cp = thermo->Cp();
|
|
volScalarField Cv = thermo->Cv();
|
|
|
|
volScalarField Ma
|
|
(
|
|
IOobject
|
|
(
|
|
"Ma",
|
|
runTime.timeName(),
|
|
mesh
|
|
),
|
|
mag(U)/(sqrt((Cp/Cv)*(Cp - Cv)*thermo->T()))
|
|
);
|
|
|
|
Info << "Max Mach number: " << max(Ma.internalField()) << endl;
|
|
Ma.write();
|
|
}
|
|
else
|
|
{
|
|
Info<< " No U" << endl;
|
|
}
|
|
}
|