[redTo]

Hello eveybody,

I'm using compressibleInterIsoFoam (OF v2006) for a nozzle spray simulation and I'm comparing the results between polyMesh and hexaMesh.

For both cases I have defined the same conditions and the results are really different (see attached file).

Any idea why is this happening?

I'm grateful for any help!

Cheers,
Liz

PS: Inlet condition: nozzle bottom 4,6 g/s
outlet condition: outlet box 1 bar

fvSolution

solvers
{
    "alpha.oil.*"
    {
        nAlphaCorr      1;
        nAlphaSubCycles 3;
        cAlpha          1;

        reconstructionScheme plicRDF;
        vof2IsoTol 	1e-8;
        surfCellTol 	1e-6; 
        nAlphaBounds    3;
        snapTol         1e-12;
        clip            true;

    }

    "pcorr.*"
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-7;
        relTol          0.01;
    }
    
    ".*(rho|rhoFinal)"
    {
        solver          diagonal;
    }

    p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-8;
        relTol          0.01;
    }

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "(U|T|k|epsilon|B).*"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-08;
        relTol          0;
    }
}

PIMPLE
{
    momentumPredictor no;    
    nOuterCorrectors 1;
    nCorrectors     3;
    nNonOrthogonalCorrectors 0;

}

relaxationFactors
{
     fields
    {
        "p.*"           0.1; 
    }
    equations
    {
        "k.*" 0.7;
	"epsilon.*" 0.7;
	"U.*" 0.8;
    }
}

fvSchemes

 ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    div(rhoPhi,U)  Gauss upwind;
    div(phi,alpha)  Gauss vanLeer;
    div(phirb,alpha) Gauss interfaceCompression;
 
    div(phi,thermo:rho.oil) Gauss upwind;
    div(phi,thermo:rho.air) Gauss upwind;

    div(phi,k)      Gauss upwind;
    div(rhoPhi,k)   Gauss upwind;
    div(phi,epsilon) Gauss upwind;
    div(rhoPhi,epsilon)  Gauss upwind;

    div(rhoPhi,T)  Gauss upwind;
    div(rhoPhi,K)  Gauss upwind;
    div(phi,p)      Gauss upwind;

    div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}