CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

inaccurate contact angle capturing in high shear/viscousity flow

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   October 23, 2015, 00:14
Default inaccurate contact angle capturing in high shear/viscousity flow
New Member
Chen Yu
Join Date: Sep 2013
Posts: 19
Rep Power: 4
chery1986 is on a distinguished road
It is a very interesting thread solved: contact angle correction in interFoam, and I think the contact angle problem is still over there after 5 years which I'm struggling with currently.

Recently, I'm working on problems with moving droplets/bubbles driven by shear flow with high viscosity (high Ca number), where the contact angle is less accurately captured than in the cases with low viscosity (low Ca number) or quasi-steady cases.

it is stated in the thread nHat and nHatf are corrected based on contact angle in interfaceProperties.C, but alpha1 and snGrad(alpha1) are still inconsistent with contact angle on the boundary when it is used to calculate pEqn and UEqn, resulting in very bad contact angle. i have several questions referring to this problem:

correctContactAngle in interfaceProperties.C has corrected nHat on patches, and set fixedGradient value of alpha1 on patches based on it. Why field of alpha1 and snGrad(alpha1) cannot reflect the correction in contact angle?

Actually, I output and checked (nHatp & np) and (gradAlphaf & nf)/mag(gradAlphaf) in correctContactAngle function in interfaceProperties.C, they are very close to each other and difference is between 1%. It shows that nHatp, nHatf and gradAlphaf are correctly set on the boundary within, but why the gradAlphaf is not used to correct the alpha1 field to match the calculated fixedGradient BC even if I call alpha1.correctBoundaryConditions() in alphaEqnSubCycle.H

BTW, I'm using constant contact angle and OF211 current, limiter of alpha1 is set as alpha and p_rgh is fixedFluxPressure

Anyone has good idea or experience on imposing contact angle BC on alpha forcely on patches?

Thanks a lot
chery1986 is offline   Reply With Quote


Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Specifying the contact angle on the interface between two zones Rui_27 Fluent Multiphase 0 July 10, 2015 05:42
Wrong interface curvature capturing and contact line calculation at low capillary nu Mahmoud_aboukhedr OpenFOAM Running, Solving & CFD 0 April 14, 2015 14:04
Slug Flow, interFoam, problems with Contact Angle PrzemekPL OpenFOAM Running, Solving & CFD 13 February 18, 2014 23:10
contact angle yan FLUENT 4 November 6, 2012 05:22
modelling inviscid 2D flow at high angle of attack Ferdinando FLUENT 2 October 30, 2007 18:26

All times are GMT -4. The time now is 11:41.