à‰tude numérique du feu généralisé avec sortie de flamme de diffusion en situation d'incendie.

B.2 Le fichier source programme

FvEqn.H

{

volScalarField alphat("alphat", turbulence->mut()/Sct );

//Concentration des suies

fvScalarMatrix NEqn

(

fvm::ddt(rho, n)

+ fvm::div(phi, n)

- fvm::laplacian(alphat, n)

(((Calpha*No*pow(rho,3)*sqrt(T)*Xc*exp(-Ta / T)) / A) - ((A*(rho*Cbeta*sqrt(T)*sqr(n)) / No)))

+ fvOptions(rho, n)

);

NEqn.relax();

fvOptions.constrain(NEqn);

NEqn.solve();

fvOptions.correct(n);

//fraction volumique des suies fvScalarMatrix FEqn

(

fvm::ddt(rho, fv)

+ fvm::div(phi, fv)

- fvm::laplacian(alphat, fv)

R*(((Cdelta / rhosoot)*Calpha*pow(rho,3)*sqrt(T)*Xc*exp(-Ta / T)) + ((Cgamma/pow(rhosoot,1 / 3))*sqr(rho)*pow(n,1 / 3)*pow(fv,2 / 3)

Rédigé par: MBAINGUEBEM Arnaud Mémoire de fin d'études

B.3 Le programme principal XII

*sqrt(T)*exp(-Tgamma/ T)))

+ fvOptions(rho, fv) );

FEqn.relax();

fvOptions.constrain(FEqn);

FEqn.solve();

fvOptions.correct(fv);

}

B.3 Le programme principal

mareactingFoam.C

#include "fvCFD.H"

#include "turbulenceModel.H"

#include "psiCombustionModel.H"

#include "multivariateScheme.H"

#include "pimpleControl.H"

#include "fvIOoptionList.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

int main(int argc, char *argv[])

{

#include "setRootCase.H"

#include "createTime.H"

#include "createMesh.H"

#include "readGravitationalAcceleration.H"

#include "createFields.H"

#include "createFvOptions.H"

#include "initContinuityErrs.H"

#include "readTimeControls.H"

#include "compressibleCourantNo.H"

#include "setInitialDeltaT.H"

pimpleControl pimple(mesh);

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nStarting time loop\n" << endl;

while (runTime.run())

{

#include "readTimeControls.H"

#include "compressibleCourantNo.H"

#include "setDeltaT.H"

runTime++;

Info<< "Time = " << runTime.timeName() << nl << endl;

#include "rhoEqn.H"

while (pimple.loop())

{

#include "UEqn.H"

#include "YEqn.H"

Rédigé par: MBAINGUEBEM Arnaud Mémoire de fin d'études

B.4 Constantes du modèle des suies XIII

Rédigé par: MBAINGUEBEM Arnaud Mémoire de fin d'études

#include "FvEqn.H" // inclusion du programme source #include "EEqn.H"

while (pimple.correct())

{

#include "pEqn.H"

}

if (pimple.turbCorr())

{

turbulence->correct();

}

}

runTime.write();

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"

<< " ClockTime = " << runTime.elapsedClockTime() << " s"

<< nl << endl;

}

Info<< "End\n" << endl;

return 0;

}

// ************************************************************************* //

B.4 Constantes du modèle des suies

Les constantes A et R ont été définies pour homogenéiser les dimensions des équations telque : A = 1mol^-1 et R = 1kg

TABLE B.1 - Les constantes des équations des transports des suies

	Constante Valeur Dimension
	Cá 6, 54.10^4kg^-2.m^3.s^-1.K^-1/^2 Câ 1,3.10^7m^3.s^-1.K^-1/^2 Cã 1, 0.10^-4kg^-2.m^3.s^-1.K^-1/^2 NA 6, 022.10^23at.mol^-1 Psuies 1800 kg.m^-3 C5 144 - Tá 46100 K Tã 12600 K Sct 0,7 - XFu 0,07