[Arne87]

Hi Foam community,

I have a problem with the driftFluxFoam solver. I used the damBreak setup from the interFoam toturial and edited the properties in a way that I can use the driftFluxFoam solver to compute it.
Although the simulation runs well,I have two major issues.

1. The second phase (alpha.sludge) is leaving the computational domain through the lowerwall patch although this patch is defined as wall. Maybe my boundary conditions are incorrect.

Here is a part of my blockMeshDict :

boundary
(
leftWall
{
type wall;
faces
(
(0 12 16 4)
(4 16 20 8)
);
}
rightWall
{
type wall;
faces
(
(7 19 15 3)
(11 23 19 7)
);
}
lowerWall
{
type wall;
faces
(
(0 1 13 12)
(1 5 17 13)
(5 6 18 17)
(2 14 18 6)
(2 3 15 14)
);
}
atmosphere
{
type patch;
faces
(
(8 20 21 9)
(9 21 22 10)
(10 22 23 11)
);
}
);

and a part of my alpha.sludge init file

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
leftWall
{
type zeroGradient;
}

rightWall
{
type zeroGradient;
}

lowerWall
{
type zeroGradient;
}

atmosphere
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

defaultFaces
{
type empty;
}
}

The alpha.sludge field is editedt with the setFieldsDict File of the damBreak test case (interFoam).

2. If I define a drift velocity of 0 (transportProperties v0 = 0) I still have a settling of the second phase. From my understanding of the mixture model the second phase is supposed to follow the first phase if there is no relative velocity.


I hope someone can help me Best

Arne