LES transient simulation of formula car nose, Part 2 of 3 (Calculation Setup)

on
Wind speed & ground speed was set to increase from 0 m/s to 50 m/s in 0.25 sec.

        inlet
        {
            type            uniformFixedValue;
            uniformValue    table
            (
                (0    ( 0 0 0))
                (0.25 (50 0 0))
            );
        }
        ground
        {
            type            uniformFixedValue;
            uniformValue    table
            (
                (0    ( 0 0 0))
                (0.25 (50 0 0))
            );
        }

And advective BC was imposed on the outlet.

        outlet
        {
            type            advective; //zeroGradient;
            phi             phi;
        }

Other setteings are as follows:

Solver: pimpleFoam (PISO mode)
SGS model: homogeneousDynSmagorinsky
DeltaT: Adjusted at every time step so that the max. Courant number is smaller than 5

controlDict:

application     pimpleFoam;
startFrom       latestTime;
startTime       0;
stopAt          endTime;
endTime         0.375;
deltaT          0.4e-5;
adjustTimeStep  yes;
maxCo             5;
maxDeltaT       1e-3;
writeControl    adjustableRunTime;
writeInterval   1e-3;
purgeWrite      0;
writeFormat     ascii;
writePrecision  7;
writeCompression uncompressed;
timeFormat      general;
timePrecision   6;
runTimeModifiable yes;
functions
{
  fieldMinMax
  {
    type               fieldMinMax;
    functionObjectLibs ("libfieldFunctionObjects.so");
    enabled            true;
    outputControl      timeStep;
    outputInterval     1;
    log                true;
    mode               magnitude;
    fields
    (
      U p
    );
  }
}

fvSchemes:

ddtSchemes
{
    default         backward;
}

gradSchemes
{
    default         Gauss linear;

    //grad(nuTilda)   cellLimited Gauss linear 1;
    grad(U)         cellLimited Gauss linear 1;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss LUST unlimitedGrad(U);
    div((nuEff*dev(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear limited corrected 0.333;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         limited corrected 0.333;
}

fluxRequired
{
    default         no;
    p               ;
}

fvSolution:

solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-08;
        relTol          1e-04;
        smoother        GaussSeidel;
        nPreSweeps      0;
        nPostSweeps     2;
        cacheAgglomeration true;
        nCellsInCoarsestLevel 10;
        agglomerator    faceAreaPair;
        mergeLevels     1;
    }

    pFinal
    {
        $p;
        tolerance       1e-08;
        relTol          0;
    }

    "U.*"
    {
        solver          PBiCG;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol          0;
    }
}

PIMPLE
{
    momentumPreditor yes;
    nOuterCorrectors 1;
    nCorrectors 2;
    nNonOrthogonalCorrectors 2;
}

LESProperties:

LESModel        homogeneousDynSmagorinsky;

delta           cubeRootVol;

turbulence      on;

printCoeffs     on;

homogeneousDynSmagorinskyCoeffs
{
    filter simple;
}

cubeRootVolCoeffs
{
    deltaCoeff      1;
}

U:

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            uniformFixedValue;
        uniformValue    table
        (
            (0    ( 0 0 0))
            (0.25 (50 0 0))
        );
    }

    outlet
    {
        type            advective; //zeroGradient;
        phi             phi;
    }

    side
    {
        type            slip;
    }

    nose
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    ground
    {
        type            uniformFixedValue;
        uniformValue    table
        (
            (0    ( 0 0 0))
            (0.25 (50 0 0))
        );
    }
}

p:

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            fixedMean; //fixedValue;
        meanValue       0;
        value           uniform 0;
    }

    side
    {
        type            slip;
    }

    nose
    {
        type            zeroGradient;
    }

    ground
    {
        type            zeroGradient;
    }
}

k:

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 1e-10;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 1e-10;
    }

    outlet
    {
        type            inletOutlet;
        inletValue      uniform 1e-10;
        value           uniform 1e-10;
    }

    nose
    {
        type            fixedValue;
        value           uniform 1e-10;
    }

    ground
    {
        type            fixedValue;
        value           uniform 1e-10;
    }

    side
    {
        type            slip;
    }
}

nuSgs:

dimensions      [0 2 -1 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            zeroGradient;
    }

    nose
    {
        type            zeroGradient;
    }

    side
    {
        type            zeroGradient;
    }

    ground
    {
        type            zeroGradient;
    }
}

Comments

Post a Comment