[saidc.]

Hi,

I've a really complex geometry and I can't create better mesh than I've already done. Because of this I need to force the foam to get accurate results but it diverges with high nNonOrthogonalCorrectors. I calculated 0.25, 0.5, 0.75, 0.99, 0.999, 0.9999 quantiles of mesh data to describe it as much as possible. I believe 0.9999 is the one that carrying more information. The simulation is laminar, closed natural convection case. All boundries are wall and has adiabatic bc. Solver is buoyantPimpleFoam.

I want to know your ideas how you would set the fvSolution and fvSchemes for this mesh. I attached meshDescription.png, please check it.

solvers
{
    "rho.*"
    {
    solver             PCG;
    preconditioner  DIC;
    tolerance         0;
    relTol              0;
    }
    p_rgh
    {
        solver             GAMG;
        tolerance        1e-06;
        relTol              0.01;
        smoother         GaussSeidel;
    }
    p_rghFinal
    {
        solver             GAMG;
        tolerance        1e-07;
        relTol              0.0;
        smoother         GaussSeidel;
    }
    "(U|h)"
    {
        solver             PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-08;
        relTol             0.0;
    }
    "(U|h)Final"
    {
        $U;
        tolerance       1e-08;
        relTol             0.0;
    }
}
PIMPLE
{
    momentumPredictor no;
    nNonOrthogonalCorrectors 2;
    nOuterCorrectors 1;
    nCorrectors         3;
    pRefCell              0;
    pRefValue           1e5;
}

ddtSchemes
{
    default         Euler;
}
gradSchemes
{
    default          Gauss linear;
    grad(U)         cellLimited Gauss linear 0.75;
}
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.75;
}
interpolationSchemes
{
    default          linear;
}
snGradSchemes
{
    default          limited 0.75;
}