Conditions aux limites

Paramètres de simulation

Fichier controlDict

application     scalarTurbulentTransportFoam;

startFrom       startTime;

startTime       0;

stopAt          endTime;

endTime         60;

deltaT          0.01;

writeControl    runTime;

writeInterval   5;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

Le pas de temps deltaT a été calculé de manière à avoir un nombre de Courant égal à 0,05 :

$$ \Delta T = \frac{C_o . \Delta x}{u} = \frac{0,05 . 1,0}{5,0} = 0,01$$

Avec une taille de maille minimale de 1 m et une vitesse de 5 m/s.

Fichier fvSolution

solvers
{
    ALPHA
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-06;
        relTol          0;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;

    residualControl
    {
          ALPHA          1e-5;
    }
}

relaxationFactors
{
    ALPHA               0.5;
}

Fichier fvSchemes

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,ALPHA)  Gauss limitedLinear 1;
}

laplacianSchemes
{
    default         none;
    laplacian((DALPHA+(((0.09*k)*k)|(epsilon*0.7))),ALPHA) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
}