Dynamic contact angle calculation in multiphaseInterFoam
 

Hi,

Could anyone tell me the implementation of dynamic contact angle in multiphaseInterfoam solver?

There are four contact angle values to be input for multiphaseInterFoam, via.. equlibrium contact angle, Utheta, Advancing contact angle, Receding contact angle.

I found some details , how contact angle is being calculated in multiphaseInterFoam as follows,
θ = (θA − θR ) ∗ tanh(uwall/uθ)

my question is when i use no slip boundary condition, u wall is zero then the theta on left hand side is also becomes zero. Then how the contact angle effect is being imposed in multiphaseinterFoam under no slip BC. First of all, is this solver requires all four contact angle to be defined?

Please help me to understand this contact angle issue.

hi,
i'am also looking for experimental values for water/air ... anybody can help?
thank very much be advance LL

Could anyone reply to my query?

If you look above the calculation for theta you will see that the patchInternalField() function is used to sample the internal field values. This is then subtracted to wall velocity (you specify no-slip and so it is equal to zero). Therefore uWall will not be zero like you suggested. I cannot find this model in literature anywhere though and it's concerning.