|
[Sponsors] |
Calculating Dynamic Contact Angle with Slip Boundary Condition |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 7, 2017, 23:57 |
Calculating Dynamic Contact Angle with Slip Boundary Condition
|
#1 |
New Member
Faraz
Join Date: Apr 2014
Location: Toronto, Canada
Posts: 16
Rep Power: 12 |
Hi All,
I need to know exactly how the dynamic contact angle is calculated for 2phase flows. As I know this is the formulation OpenFOAM uses: θ_d=θ_s+(θ_A-θ_R)∙tanh(U_wall) Where θ_d, θ_s, θ_A, and θ_R are dynamic, static, advancing, and receding contact angles, respectively, and U_wall is the interface velocity component parallel to the wall (here): HTML Code:
Foam::dynamicAlphaContactAngleFvPatchScalarField::theta
(
const fvPatchVectorField& Up,
const fvsPatchVectorField& nHat
) const
{
if (uTheta_ < SMALL)
{
return tmp<scalarField>(new scalarField(size(), theta0_));
}
const vectorField nf(patch().nf());
// Calculated the component of the velocity parallel to the wall
vectorField Uwall(Up.patchInternalField() - Up);
Uwall -= (nf & Uwall)*nf;
// Find the direction of the interface parallel to the wall
vectorField nWall(nHat - (nf & nHat)*nf);
// Normalise nWall
nWall /= (mag(nWall) + SMALL);
// Calculate Uwall resolved normal to the interface parallel to
// the interface
scalarField uwall(nWall & Uwall);
return theta0_ + (thetaA_ - thetaR_)*tanh(uwall/uTheta_);
}
vectorField Uwall(Up.patchInternalField() - Up); Many thanks in advance,
__________________
Thermofluids for Energy and Advanced Materials (TEAM) Laboratory |
|
September 8, 2017, 04:42 |
|
#2 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 |
Up.patchInternalField is the velocity in the cell center next to the boundary face you are looking at. There is no explicit interface velocity. In case of a stationary no-slip wall, Up will be zero and Uwall is the internal field value. In case of a moving no-slip wall, you remove that contribution. In case of perfect slip, your contact angle becomes zero. In case of partial slip, you will have something in between.
Regards to Prof. Bazylak.
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
September 8, 2017, 12:41 |
|
#3 | |
New Member
Faraz
Join Date: Apr 2014
Location: Toronto, Canada
Posts: 16
Rep Power: 12 |
Quote:
Thanks for your reply. So, in case of no-slip B.C., there is still a Uwall velocity and dynamic contact angle is not equal to the static contact angle, right? Though because Uwall is regularly very small in no-slip case, the dynamic contact angle should be very close to the static contact angle, right? Therefore, is it correct to say, using dynamic contact angle only makes sense when we use a slip/partial-slip B.C.? Also, I did not understand why in the fully slip case, the contact angle becomes zero. I can say Uwall can be very large, thus its tanh will be 1, and therefore θ_d=θ_s+(θ_A-θ_R), am I wrong? Finally, as I know the advancing contact angle must be the maximum, and receding contact angle must be the minimum angle the interface can reach during motion; however, at large Uwalls, θ_d=θ_s+(θ_A-θ_R) will not give θ_A! Thanks, Faraz
__________________
Thermofluids for Energy and Advanced Materials (TEAM) Laboratory |
||
September 12, 2017, 04:14 |
|
#4 |
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,377
Rep Power: 30 |
Uwall is generally not zero, right. I would not dare to say for certain it is generally small. For one, it depends on your mesh resolution. Uwall is zero in static cases (no flow), and then the "dynamic" contact angle becomes the static contact angle (shocking, I know ).
In case of a perfect slip, Uwall is not going to be large, it will be zero parallel to the wall (since a slip boundary condition is defined by zero gradient in parallel and zero normal velocity). Thus your dynamic contact angle should be your static contact angle in this case too (so you are correct about it not being zero).
__________________
*On twitter @akidTwit *Spend as much time formulating your questions as you expect people to spend on their answer. |
|
Tags |
contact angle, dynamic contact angle, interface, multiphase, slip b.c. |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind turbine simulation | Saturn | CFX | 60 | July 17, 2024 06:45 |
sliding mesh problem in CFX | Saima | CFX | 46 | September 11, 2021 08:38 |
Radiation in semi-transparent media with surface-to-surface model? | mpeppels | CFX | 11 | August 22, 2019 08:30 |
dynamic contact angle / contact line velocity | ThomasN | CFX | 6 | September 23, 2014 07:07 |
turbineSiting tutorial: slip condition at top boundary gives unexpected result | letzel | OpenFOAM | 0 | June 6, 2014 06:25 |