[hbulus]

Dont be modest. Thanks anyway.

Have good days!

[saidc.]

Hi all,

I'm having same trouble with a quite bad mesh. I've done the simulation with various schemes but non of them gave the good resuts. Velocity profiles are 10 times wrong from the experimental results while temperature profiles and Nusselt Numbers are close enough. It is high density sphere packing geometry and meshing whole pack is really hard to manage. I attached the properties of tetrahedral mesh.

• Should I go for finer mesh? What would you do for getting better results?

• Isn't it better treating schemes and linear solver with respect to high quantiles datas than max values of checkMesh?

application     buoyantSimpleFoam;
startFrom       startTime;
startTime       0;
stopAt          endTime;
endTime         10000;
deltaT          1;
writeControl    timeStep;
writeInterval   20;
purgeWrite      3;
writeFormat     ascii;
writePrecision  10;
writeCompression off;
timeFormat      general;
timePrecision   10;
runTimeModifiable true;

 gradSchemes
{
    default         cellMDLimited Gauss linear 0.5;
    grad(U)         cellMDLimited Gauss linear 0.5;
}
divSchemes
{
    default          none;
    div(phi,U)       Gauss linearUpwind grad(U);
    div(phi,h)       Gauss upwind;
    div(phi,K)       Gauss upwind;
    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
    default          Gauss linear limited 0.5;
}
interpolationSchemes
{
    default         linear;
}
snGradSchemes
{
    default         limited 0.5;
}

 solvers
{
    p_rgh
    {
        solver           GAMG;
        tolerance        1e-06;
        relTol           0.0;
        smoother         GaussSeidel;
        nPreSweeps       0;
        nPostSweeps      2;
        cacheAgglomeration      on;
        agglomerator            faceAreaPair;
        nCellsInCoarsestLevel   1000;
        mergeLevels      1;
    }
    "(U|h)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol          0.0;
    }
}
SIMPLE
{
    momentumPredictor no;
    nNonOrthogonalCorrectors 1;
    pRefCell        0;
    pRefValue       1e5;
    residualControl
    {
        p_rgh        1e-3;
        h              1e-3;
    }
}
relaxationFactors
{
    fields
    {
        rho             1.0;
        p_rgh           0.3;
    }
    equations
    {
        U               0.7;
        h               0.7;
    }
}

Kind regards,
Said