[ndtrong]

cylinder-01.jpg

cylinder-06.jpg

alpha.txt

Hi everyone,

Hope that you are fine!

I am facing problem with boundary condtion for alpha in simulation of flow with a cylinder near interface.
The simulated results of alpha is wrong and I think that is because of boundary condition.
Could someone kindly help me to fix it?

Thanks

[HPE]

what is the driving force of the simulation? how are the BCs for the other quantities?

[ndtrong]

Originally Posted by HPE

what is the driving force of the simulation? how are the BCs for the other quantities?

Thanks for your comment,

The boundary conditions for the other quantities is as below:
For p_rgh
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1906 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
#includeEtc "caseDicts/setConstraintTypes"

inlet
{
type fixedFluxPressure;
value $internalField;
}

outlet
{
type zeroGradient;
}

atmosphere
{
type totalPressure;
p0 uniform 0;
}

hull
{
type fixedFluxPressure;
value $internalField;
}

bottom
{
type fixedFluxPressure;
value $internalField;
}
}

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


For U

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1906 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Umean 0.395;
mUmean -0.395;

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

internalField uniform ($Umean 0 0);

boundaryField
{
// #includeEtc "caseDicts/setConstraintTypes"

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type outletPhaseMeanVelocity;
alpha alpha.water;
Umean $Umean;
value $internalField;
}

atmosphere
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

hull
{
type movingWallVelocity;
value uniform (0 0 0);
}

bottom
{
type noSlip;
}
}


// ************************************************** *********************** //
For k
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1906 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0.00015;

boundaryField
{
#includeEtc "caseDicts/setConstraintTypes"

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}

atmosphere
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}

hull
{
type kqRWallFunction;
value $internalField;
}
bottom
{
type kqRWallFunction;
value $internalField;
}
}

// ************************************************** *********************** //
For omega
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1906 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 2;

boundaryField
{
#includeEtc "caseDicts/setConstraintTypes"

inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}

atmosphere
{
type inletOutlet;
inletValue $internalField;
value $internalField;
}

hull
{
type omegaWallFunction;
value $internalField;
}

bottom
{
type omegaWallFunction;
value $internalField;
}
}



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