9a536b02a7
fvOptions are transferred to the database on construction using fv::options::New which returns a reference. The same function can be use for construction and lookup so that fvOptions are now entirely demand-driven. The abstract base-classes for fvOptions now reside in the finiteVolume library simplifying compilation and linkage. The concrete implementations of fvOptions are still in the single monolithic fvOptions library but in the future this will be separated into smaller libraries based on application area which may be linked at run-time in the same manner as functionObjects.
150 lines
4.5 KiB
C
150 lines
4.5 KiB
C
/*---------------------------------------------------------------------------*\
|
|
========= |
|
|
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
|
\\ / O peration |
|
|
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
|
|
\\/ M anipulation |
|
|
-------------------------------------------------------------------------------
|
|
License
|
|
This file is part of OpenFOAM.
|
|
|
|
OpenFOAM is free software: you can redistribute it and/or modify it
|
|
under the terms of the GNU General Public License as published by
|
|
the Free Software Foundation, either version 3 of the License, or
|
|
(at your option) any later version.
|
|
|
|
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
|
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
for more details.
|
|
|
|
You should have received a copy of the GNU General Public License
|
|
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
|
|
|
Application
|
|
engineFoam
|
|
|
|
Description
|
|
Solver for internal combustion engines.
|
|
|
|
Combusting RANS code using the b-Xi two-equation model.
|
|
Xi may be obtained by either the solution of the Xi transport
|
|
equation or from an algebraic exression. Both approaches are
|
|
based on Gulder's flame speed correlation which has been shown
|
|
to be appropriate by comparison with the results from the
|
|
spectral model.
|
|
|
|
Strain effects are encorporated directly into the Xi equation
|
|
but not in the algebraic approximation. Further work need to be
|
|
done on this issue, particularly regarding the enhanced removal rate
|
|
caused by flame compression. Analysis using results of the spectral
|
|
model will be required.
|
|
|
|
For cases involving very lean Propane flames or other flames which are
|
|
very strain-sensitive, a transport equation for the laminar flame
|
|
speed is present. This equation is derived using heuristic arguments
|
|
involving the strain time scale and the strain-rate at extinction.
|
|
the transport velocity is the same as that for the Xi equation.
|
|
|
|
\*---------------------------------------------------------------------------*/
|
|
|
|
#include "fvCFD.H"
|
|
#include "engineTime.H"
|
|
#include "engineMesh.H"
|
|
#include "psiuReactionThermo.H"
|
|
#include "turbulentFluidThermoModel.H"
|
|
#include "laminarFlameSpeed.H"
|
|
#include "ignition.H"
|
|
#include "Switch.H"
|
|
#include "OFstream.H"
|
|
#include "mathematicalConstants.H"
|
|
#include "pimpleControl.H"
|
|
#include "fvOptions.H"
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
int main(int argc, char *argv[])
|
|
{
|
|
#include "setRootCase.H"
|
|
|
|
#include "createEngineTime.H"
|
|
#include "createEngineMesh.H"
|
|
|
|
pimpleControl pimple(mesh);
|
|
|
|
#include "readCombustionProperties.H"
|
|
#include "createFields.H"
|
|
#include "createMRF.H"
|
|
#include "createFvOptions.H"
|
|
#include "createRhoUf.H"
|
|
#include "initContinuityErrs.H"
|
|
#include "readEngineTimeControls.H"
|
|
#include "compressibleCourantNo.H"
|
|
#include "setInitialDeltaT.H"
|
|
#include "startSummary.H"
|
|
|
|
turbulence->validate();
|
|
|
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
|
|
|
Info<< "\nStarting time loop\n" << endl;
|
|
|
|
while (runTime.run())
|
|
{
|
|
#include "readEngineTimeControls.H"
|
|
#include "compressibleCourantNo.H"
|
|
#include "setDeltaT.H"
|
|
|
|
runTime++;
|
|
|
|
Info<< "Crank angle = " << runTime.theta() << " CA-deg" << endl;
|
|
|
|
mesh.move();
|
|
|
|
#include "rhoEqn.H"
|
|
|
|
// --- Pressure-velocity PIMPLE corrector loop
|
|
while (pimple.loop())
|
|
{
|
|
#include "UEqn.H"
|
|
|
|
#include "ftEqn.H"
|
|
#include "bEqn.H"
|
|
#include "EauEqn.H"
|
|
#include "EaEqn.H"
|
|
|
|
if (!ign.ignited())
|
|
{
|
|
thermo.heu() == thermo.he();
|
|
}
|
|
|
|
// --- Pressure corrector loop
|
|
while (pimple.correct())
|
|
{
|
|
#include "pEqn.H"
|
|
}
|
|
|
|
if (pimple.turbCorr())
|
|
{
|
|
turbulence->correct();
|
|
}
|
|
}
|
|
|
|
#include "logSummary.H"
|
|
|
|
rho = thermo.rho();
|
|
|
|
runTime.write();
|
|
|
|
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
|
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
|
<< nl << endl;
|
|
}
|
|
|
|
Info<< "End\n" << endl;
|
|
|
|
return 0;
|
|
}
|
|
|
|
|
|
// ************************************************************************* //
|