Annexe B
Les programmes sources
d'introduction des équations de
transport des suies
B.1 Le fichier d'entête createFields.H
Ce fichier permet de faire la déclaration des variables et
des constantes.
Info<< "Creating reaction model\n" << endl;
autoPtr<combustionModels::psiCombustionModel> reaction
(
combustionModels::psiCombustionModel::New(mesh)
);
psiReactionThermo& thermo = reaction->thermo();
thermo.validate(args.executable(), "h", "e");
basicMultiComponentMixture& composition =
thermo.composition();
PtrList<volScalarField>& Y = composition.Y();
word inertSpecie(thermo.lookup("inertSpecie"));
volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh
),
thermo.rho()
);
Info<< "Reading field U\n" << endl; //
déclaration de U
volVectorField U
(
IOobject
(
"U",
B.1 Le fichier d'entête createFields.H VIII
Rédigé par: MBAINGUEBEM Arnaud Mémoire de
fin d'études
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Reading field fv\n" << endl; //
déclaration de fv
volScalarField fv
(
IOobject
(
"fv",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Reading field n\n" << endl; //
déclaration de n
volScalarField n
(
IOobject
(
"n",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);
Info<< "Reading transportProperties\n" << endl; //
Déclaration et définition du repertoire des constantes
IOdictionary transportProperties
(
IOobject
(
"transportProperties", runTime.constant(), mesh,
IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE
)
);
Info<< "Reading turbulent Schmidt number Sct\n" <<
endl; dimensionedScalar Sct
B.1 Le fichier d'entête createFields.H IX
(
transportProperties.lookup("Sct")
);
Info<< "Reading Avogadro number No\n" << endl;
dimensionedScalar No
(
transportProperties.lookup("No")
);
Info<< "Reading soot density rhosoot\n" <<
endl;
dimensionedScalar rhosoot
(
transportProperties.lookup("rhosoot")
);
Info<< "Reading coefficient Calphat\n" << endl;
dimensionedScalar Calphat
(
transportProperties.lookup("Calphat")
);
Info<< "Reading coefficient Cbeta\n" << endl;
dimensionedScalar Cbeta
(
transportProperties.lookup("Cbeta")
);
Info<< "Reading coefficient Cdelta\n" << endl;
dimensionedScalar Cdelta
(
transportProperties.lookup("Cdelta")
);
Info<< "Reading coefficient Cgamma\n" << endl;
dimensionedScalar Cgamma
(
transportProperties.lookup("Cgamma")
);
Info<< "Reading Temperature Ta\n" << endl;
dimensionedScalar Ta
(
transportProperties.lookup("Ta")
);
Info<< "Reading Temperature Ty\n" << endl;
dimensionedScalar Ty
(
transportProperties.lookup("Ty")
);
Info<< "Reading Avogadro inverse A\n" << endl; //
rendre l'équation homogène
dimensionedScalar A
(
transportProperties.lookup("A")
);
Rédigé par: MBAINGUEBEM Arnaud Mémoire de
fin d'études
B.1 Le fichier d'entête createFields.H X
Info<< "Reading constant B\n" << endl;
dimensionedScalar B
(
transportProperties.lookup("B")
);
Info<< "Reading molar frcation of fuel Xfu\n" <<
endl;
dimensionedScalar Xfu
(
transportProperties.lookup("Xfu")
);
volScalarField& p = thermo.p();
const volScalarField& psi = thermo.psi();
const volScalarField& T = thermo.T();
#include "compressibleCreatePhi.H"
Info << "Creating turbulence model.\n" << nl;
autoPtr<compressible::turbulenceModel> turbulence
(
compressible::turbulenceModel::New
(
rho,
U,
phi,
thermo
)
);
Info<< "Reading field alphat\n" << endl;
// Set the turbulence into the reaction model
reaction->setTurbulence(turbulence());
Info<< "Creating field dpdt\n" << endl;
volScalarField dpdt
(
IOobject
(
"dpdt",
runTime.timeName(),
mesh
),
mesh,
dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
);
Info<< "Creating field kinetic energy K\n" <<
endl;
volScalarField K("K", 0.5*magSqr(U));
multivariateSurfaceInterpolationScheme<scalar>::fieldTable
fields;
forAll(Y, i)
{
fields.add(Y[i]);
}
fields.add(thermo.he());
Rédigé par: MBAINGUEBEM Arnaud Mémoire de
fin d'études
B.2 Le fichier source programme XI
volScalarField dQ
(
IOobject
(
"dQ",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
);
| 
