applications/solvers/combustion/engineFoam/engineFoam.C

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
   6      \\/     M anipulation  |
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of OpenFOAM.
  10
  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
  13     the Free Software Foundation, either version 3 of the License, or
  14     (at your option) any later version.
  15
  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
  19     for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 Application
  25     engineFoam
  26
  27 Description
  28     Solver for internal combustion engines.
  29
  30     Combusting RANS code using the b-Xi two-equation model.
  31     Xi may be obtained by either the solution of the Xi transport
  32     equation or from an algebraic exression.  Both approaches are
  33     based on Gulder's flame speed correlation which has been shown
  34     to be appropriate by comparison with the results from the
  35     spectral model.
  36
  37     Strain effects are encorporated directly into the Xi equation
  38     but not in the algebraic approximation.  Further work need to be
  39     done on this issue, particularly regarding the enhanced removal rate
  40     caused by flame compression.  Analysis using results of the spectral
  41     model will be required.
  42
  43     For cases involving very lean Propane flames or other flames which are
  44     very strain-sensitive, a transport equation for the laminar flame
  45     speed is present.  This equation is derived using heuristic arguments
  46     involving the strain time scale and the strain-rate at extinction.
  47     the transport velocity is the same as that for the Xi equation.
  48
  49 \*---------------------------------------------------------------------------*/
  50
  51 #include "fvCFD.H"
  52 #include "engineTime.H"
  53 #include "engineMesh.H"
  54 #include "hhuCombustionThermo.H"
  55 #include "turbulenceModel.H"
  56 #include "laminarFlameSpeed.H"
  57 #include "ignition.H"
  58 #include "Switch.H"
  59 #include "OFstream.H"
  60 #include "mathematicalConstants.H"
  61 #include "pimpleControl.H"
  62
  63 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  64
  65 int main(int argc, char *argv[])
  66 {
  67     #include "setRootCase.H"
  68
  69     #include "createEngineTime.H"
  70     #include "createEngineMesh.H"
  71     #include "readCombustionProperties.H"
  72     #include "createFields.H"
  73     #include "initContinuityErrs.H"
  74     #include "readEngineTimeControls.H"
  75     #include "compressibleCourantNo.H"
  76     #include "setInitialDeltaT.H"
  77     #include "startSummary.H"
  78
  79     pimpleControl pimple(mesh);
  80
  81     // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  82
  83     Info<< "\nStarting time loop\n" << endl;
  84
  85     while (runTime.run())
  86     {
  87         #include "readEngineTimeControls.H"
  88         #include "compressibleCourantNo.H"
  89         #include "setDeltaT.H"
  90
  91         runTime++;
  92
  93         Info<< "Crank angle = " << runTime.theta() << " CA-deg" << endl;
  94
  95         mesh.move();
  96
  97         #include "rhoEqn.H"
  98
  99         // --- Pressure-velocity PIMPLE corrector loop
 100         for (pimple.start(); pimple.loop(); pimple++)
 101         {
 102             #include "UEqn.H"
 103
 104             #include "ftEqn.H"
 105             #include "bEqn.H"
 106             #include "huEqn.H"
 107             #include "hEqn.H"
 108
 109             if (!ign.ignited())
 110             {
 111                 hu == h;
 112             }
 113
 114             // --- PISO loop
 115             for (int corr=1; corr<=pimple.nCorr(); corr++)
 116             {
 117                 #include "pEqn.H"
 118             }
 119
 120             if (pimple.turbCorr())
 121             {
 122                 turbulence->correct();
 123             }
 124         }
 125
 126         #include "logSummary.H"
 127
 128         rho = thermo.rho();
 129
 130         runTime.write();
 131
 132         Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
 133             << "  ClockTime = " << runTime.elapsedClockTime() << " s"
 134             << nl << endl;
 135     }
 136
 137     Info<< "End\n" << endl;
 138
 139     return 0;
 140 }
 141
 142
 143 // ************************************************************************* //