1 /*---------------------------------------------------------------------------*\
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
5 \\ / A nd | Copyright (C) 1991-2009 OpenCFD Ltd.
7 -------------------------------------------------------------------------------
9 This file is part of OpenFOAM.
11 OpenFOAM is free software; you can redistribute it and/or modify it
12 under the terms of the GNU General Public License as published by the
13 Free Software Foundation; either version 2 of the License, or (at your
14 option) any later version.
16 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
21 You should have received a copy of the GNU General Public License
22 along with OpenFOAM; if not, write to the Free Software Foundation,
23 Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
29 Calculate the geometric compression ratio.
30 Note that if you have valves and/or extra volumes it will not work,
31 since it calculates the volume at BDC and TCD.
33 \*---------------------------------------------------------------------------*/
36 #include "engineTime.H"
37 #include "engineMesh.H"
39 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
41 int main(int argc, char *argv[])
43 #include "setRootCase.H"
44 #include "createEngineTime.H"
45 #include "createEngineMesh.H"
47 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
50 scalar fullCycle = 360.0;
55 while (runTime.theta() > ca0)
61 while (mag(runTime.theta() - ca0) > eps)
63 scalar t0 = runTime.userTimeToTime(ca0 - runTime.theta());
64 runTime.setDeltaT(t0);
66 Info<< "CA = " << runTime.theta() << endl;
70 scalar Vmax = sum(mesh.V().field());
72 while (mag(runTime.theta()-ca1) > eps)
74 scalar t1 = runTime.userTimeToTime(ca1-runTime.theta());
75 runTime.setDeltaT(t1);
77 Info<< "CA = " << runTime.theta() << endl;
81 scalar Vmin = sum(mesh.V().field());
83 Info<< "\nVmax = " << Vmax;
84 Info<< ", Vmin = " << Vmin << endl;
85 Info<< "Vmax/Vmin = " << Vmax/Vmin << endl;
86 Info<< "\nEnd" << endl;
92 // ************************************************************************* //