contact angle / multiphaseInterFoam
 

Hallo everybody!

I want to use the contact angle definition for walls with the multiphaseInterFoam solver.

As basis for my own modell I chose the dambreak4phase tutorial but I don't know how to modify the boundary conditions to get the alphaContactAngle working.

When I searched the Forum I only found out that I have to modify the alpha... files in 0/ directory before running setFields but what is the exactly code for it?

Greatings

Andreas

One step further
 

Hello, it's me again!

I think that I reached the next level with my model ;)

In the damBreak4Phase-tutorial one can see the correct definition for the alphaContactAngle at the end of the file:
type alphaContactAngle;
thetaProperties ( ( water air ) 90 0 0 0 ... );
value uniform 0;

I would assume that the values of the thetaProps are defined in the user guide of OpenFoam.
But the next question is why the value must be 0 and not 1, what makes the difference?

P.S. Obviously there is a bug in setFields, because running setFields with a "fresh" alphaPhase file results in a mixed up List of phase entries in this file.
Here funkySetFields seems to work OK (I use the -keepPatches option), so thank you very much for this little program!

Greetings

Andreas

Hey cookie!

i read your post carefully. however, it is not clear to me which file you have edited. Could please post the exact procedure to set up the alphaContactAngle boundary condition.

kinda regrards

alphaContactAngle
 

Hi Claus and thanks for your reply!

You're right, I missed to post the name of the file.
But now, I think, it works:
I chose the 0/alphaAir.org file to set the contact angles for all phases in my own model:

boundaryField
{
boundary_1
{
type alphaContactAngle;
thetaProperties ( ( Water Air ) 90 0 0 0 );
value uniform 0;
}

boundary_2
{
type alphaContactAngle;
thetaProperties ( ( Water Air ) 90 0 0 0 );
value uniform 0;
}
boundary_3
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
defaultFaces
{
type empty;
}
}

There are two phases, Water and Air and in this case, the static contact angle is 90°.
Please be aware that this is only an example, but you can find it also at the end of damBreak4Phase/0/alphaair.
Than you can run funkySetFields to set the positions of the phases and afterwards you will find this text-structure at the end of the alphaAir file.

When you use setFields, this may be a problem and you have to correct the position manually.

Please post, if you have any further questions.

Kind regards

Andreas

Hallo cookie

thanks for your post, it helped me a lot. However, I have just started to get involved into multiphaseInterFoam solver and I am not very far. Therefore, I want to ask you if you have done a simulation with channal flows?

cheers,

Claus

Channel flow
 

Hi Claus,

thank you very much!
Unfortunately I must say, that I'm quite new to openFoam, so I do not have any greater experience. That also means, that I haven't worked with channel flow so far, sorry for that.

Greetings Andreas

FunkySetFields and alphaContactAngle
 

Hi cookie

I hope after 11 years you still look at this message.

I have very similar problem with "alphaContactAngle" in alpha.air. unfortunately "funkySetFields" does not recognise "alphaContactAngle" as a known patchField. can you please clarify more about "-keepPatches option" that you mentioned earlier. how I can apply it in funkySetFields.

My Erorr:

--> FOAM FATAL IO ERROR:
Unknown patchField type alphaContactAngle for patch type patch

Valid patchField types are :

118