CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

interFoam with dynamic contact angle: distorted interface

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

Reply
 
LinkBack Thread Tools Display Modes
Old   April 18, 2012, 13:10
Default interFoam with dynamic contact angle: distorted interface
  #1
Member
 
Michiel
Join Date: Oct 2010
Location: Delft, Netherlands
Posts: 97
Rep Power: 6
michielm is on a distinguished road
Hi all,
I am using interFoam to simulate a droplet rolling down a tilted plate. I started out using the static contact angle model and all seemed fine*.

Then I switched to the dynamic contact angle model, setting the parameters to get as close as possible to the Cox-Voinov relation.

This resulted in an unstable/distorted interface (near the contact line) as you can see in the attached image.

Does anyone know what might be the cause of this distortion which is clearly unphysical?

Thanks,

Michiel


* (except that the lack of dynamic CAs gave non-realistic terminal velocities)
Attached Images
File Type: jpg Static_vs_dynamicCa.jpg (59.3 KB, 77 views)
michielm is offline   Reply With Quote

Old   April 19, 2012, 13:54
Default
  #2
New Member
 
Nikhil
Join Date: Sep 2011
Posts: 9
Rep Power: 5
Nikhilcfd is on a distinguished road
Hello,

Are you using the dynamic contact angle expression (theta0+ (thetaA-thetaR)*tanh(uwall/utheta) in OpenFoam? Non realistic (or non zero) velocities might be parasitic currents (which I think originates from curvature calculation). I am not sure why a Dynamic CA implementation would cause such a distortion. Did you use the same mesh for both the simulations. What did you input for thetaA, thetaR and utheta.??
Nikhilcfd is offline   Reply With Quote

Old   April 20, 2012, 01:34
Default
  #3
Member
 
Michiel
Join Date: Oct 2010
Location: Delft, Netherlands
Posts: 97
Rep Power: 6
michielm is on a distinguished road
I used exactly the same mesh, BCs and numerical schemes. The only difference is the static vs dynamic contact angle model.

I indeed used the standard implementation in OF and adjusted the parameters (thetaA-thetaR) and utheta to approximate the cox-voinov relation (theta_dyn=(theta0^3+100*Ca).^(1/3)).

The values I used for that where
- thetaA-thetaR=30.6 (thetaA=50.6 and thetaR=20 to be precise)
- uTheta=0.0255
michielm is offline   Reply With Quote

Old   April 24, 2012, 01:40
Default
  #4
Member
 
Michiel
Join Date: Oct 2010
Location: Delft, Netherlands
Posts: 97
Rep Power: 6
michielm is on a distinguished road
I managed to solve the problem. In the end it turned out to be something rather trivial:
I had defined my thetaA and thetaR in the physically correct order (thetaA>thetaR), but had not noticed that the definition of uwall in openfoam depends on the surface normal, which switches sign depending on which phase I use as the droplet (alpha=1 or alpha=0).

I have now swapped the values of thetaA and thetaR and the system behaves as it should (at least numerically, physically I still have to check)
michielm is offline   Reply With Quote

Old   December 20, 2012, 08:54
Default
  #5
Member
 
Indrajit
Join Date: Nov 2012
Posts: 33
Rep Power: 4
IndrajitW is on a distinguished road
Hi michielm,
I wish to start with the simulation of droplet rolling on an inclined plate.However I am new to OpenFoam and am lost where to start from??? I have done the Dam break cavity and now I wish to go for droplet rolling on inclined plate.Please guide as to where I may start.
Indrajit
IndrajitW is offline   Reply With Quote

Reply

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
contact angle yan FLUENT 4 November 6, 2012 05:22
Concentric tube heat exchanger (Air-Water) Young CFX 5 October 6, 2008 23:17
Dynamic contact angle Aireen FLUENT 2 July 5, 2006 13:14
Convective Heat Transfer - Heat Exchanger Mark CFX 6 November 15, 2004 16:55
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 21:09


All times are GMT -4. The time now is 15:59.