Here are the files of the OpenFOAM case:

fvSchemes:

/*--------------------------------*- C++ -*----------------------------------*

| ========= | |

| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |

| \ / O peration | Version: 1.6 |

| \ / A nd | Web: www.OpenFOAM.org |

| \/ M anipulation | |

*---------------------------------------------------------------------------*/

FoamFile

{

version 2.0;

format ascii;

class dictionary;

location â€śsystemâ€ť;

object fvSchemes;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes

{

default steadyState;

}

gradSchemes

{

default Gauss linear;

gradÂ§ Gauss linear;

grad(U) cellLimited Gauss linear 1;

// gradU cellLimited Gauss linear 1;

// grad(nuTilda) cellLimited Gauss linear 1;

// gradNuTilda cellLimited Gauss linear 1;

}

divSchemes

{

default Gauss linear;

div(phi,U) bounded Gauss linearUpwind gradU;

div(phi,nuTilda) bounded Gauss linearUpwind gradU;

//div(phi,nuTilda) bounded Gauss upwind;

div((nuEff*dev(grad(U).T()))) Gauss linear;

}

laplacianSchemes

{

default Gauss linear corrected; // Gauss linear limited corrected 0.33;

}

interpolationSchemes

{

default linear;

interpolate(U) linear;

}

snGradSchemes

{

default corrected; //limited corrected 0.33;

}

fluxRequired

{

default no;

p ;

Phi ;

}

// ************************************************************************* //

fvSolution:

/*--------------------------------*- C++ -*----------------------------------*

| ========= | |

| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |

| \ / O peration | Version: 2.4.0 |

| \ / A nd | Web: www.OpenFOAM.org |

| \/ M anipulation | |

*---------------------------------------------------------------------------*/

FoamFile

{

version 2.0;

format ascii;

class dictionary;

object fvSolution;

}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers

{

â€śÂ§â€ť

{

solver GAMG;

agglomerator faceAreaPair;

mergeLevels 1;

cacheAgglomeration on;

nCellsInCoarsestLevel 200;

tolerance 1e-11;

relTol 0.01;

smoother GaussSeidel;

nPreSweeps 0;

nPostSweeps 2;

nFinestSweeps 2;

minIter 1;

}

```
"(U|nuTilda)"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-15;
relTol 0.01;
}
```

// Phi

// {

// $p;

// }

}

SIMPLE

{

nNonOrthogonalCorrectors 0;

}

potentialFlow

{

nNonOrthogonalCorrectors 10;

}

relaxationFactors

{

fields

{

p 0.2;

}

equations

{

U 0.5;

nuTilda 0.5;

}

}

cache

{

grad(U);

}

// ************************************************************************* //