[chaudhad]

Hi guys,
I am doing a 2-phase simulation using compressibleInterFoam with LES as my turbulence model. The geometry is quite simple. it is a cylinder of 25.9mm radius and 906 mm length. Air is initially filled inside and then water starts to fill the cylinder from the bottom to compress the air inside. Simulations have already been done before with this case but they used blockMesh to mesh the geometry. I am using snappyHexMesh to mesh my geometry but when I make my mesh coarse without adding the boundary layers, the results are a bit different, and when I try to make it fine with the addition of boundary layers like 15 in this case, the solver stops without any error. I can also see that it is taking very large iteration to converge for p_rho, but I don't know what can I do. I have attached all the files below. please any help would be highly appreciated.

o folder


alphat

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1912 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 1e-11;

boundaryField
{
"(cylinderWalls|topX|bottomX|topY|bottomY|topZ|bot tomZ)"
{
type compressible::alphatWallFunction;
value uniform 0;
}

pipeWalls
{
type compressible::alphatWallFunction;
value uniform 0;
}

top
{
type compressible::alphatWallFunction;
value uniform 0;
}

inlet
{
type calculated;
value uniform 0;
}
}


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

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1912 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 1e-11;

boundaryField
{


"(cylinderWalls|topX|bottomX|topY|bottomY|topZ|bot tomZ)"
{
type nutUWallFunction;
value uniform 0;
}

pipeWalls
{
type nutUWallFunction;
value uniform 0;
}

top
{
type nutUWallFunction;
value uniform 0;
}

inlet
{
type calculated;
value uniform 0;
}
}

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