# Dynamic contact angle

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 April 8, 2009, 11:04 Dynamic contact angle #1 New Member     Roozbeh Mousavi B. T. Join Date: Mar 2009 Posts: 2 Rep Power: 0 Dear FLUENT users; I want to write a UDF for taking the dynamic contact angle into account. As you know, a contact line is the intersection of a solid wall and the interface of two fluids. The contact line is the angle between the wall surface and the tangent to the fluids interface at a point on the contact line. So, because that you have infinite points on the contact line, the contact angle is a local variable. There are several corrolation relating contact angle and the contact line speed (the fluid velocity at contact line). I can write a UDF using a Define_Profile function, but I do not know that if I speak about the cell_velocity in my UDF, does FLUENT adopt the velocity of the cell for which the value of contact line is being calculated? Many thanks for your attention, Roozbeh

 April 14, 2009, 19:05 #2 New Member   Ayyoub Join Date: Apr 2009 Posts: 2 Rep Power: 0 Hi Roozbeh. I am working on VOF model for a while. This would be really nice while in Fluent typically we are setting the static contact angle. I think that you are finding an element with in which volume fraction is about 50% on the wall (interface place) and then extracting the x-velocity of that cell as the velocity of the interface? Am I correct?

 July 22, 2009, 10:43 #3 New Member   davide Join Date: Jul 2009 Posts: 2 Rep Power: 0 Hi is there anyone around with a working UDF for setting the value of the contact angle (CA) on wall faces as a time dependent function ? or as a delta_CA increase with respect the CA value in the previous iteration Thank you in advance

 August 10, 2009, 07:43 #4 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,832 Rep Power: 100 Hi: Ayyoub's point highlights the dilemma here - with a no slip boundary condition the velocity at the wall is zero. So how then does the free surface move along the surface? If you use the velocity of the first cell in the domain then the velocity will be dependent on mesh size and your simulation will not converge as you refine the mesh. Simulating a moving contact line along a no slip wall with a specified contact angle is an area of active research and I have not yet seen a "universal" approach. What I do know is that if you refine the mesh of a simulation which uses this effect to drive the flow (eg a capilliary driven flow) it is impossible to achieve a mesh converged solution - in fact all my tests are diverging with mesh refinement. This means, in short, that the implementation of wall contact angles in Fluent (and also most other CFD codes I should add) is not physical and is not accurate. A leading researcher in this field is Shikhmurzaev (http://web.mat.bham.ac.uk/Y.D.Shikhmurzaev/) Has anyone else looked into this issue? Have you checked the mesh size independence of these flows? Glenn Agad15 likes this.

 November 1, 2010, 17:11 Dca #5 New Member   Concordia_CFD Join Date: Jul 2010 Location: Canada Posts: 24 Rep Power: 8 Hello! As Glenn pointed out there is a singularity at the contact region where 3 phases meet. Hence, most numerical models make assunptions to calculate the contact line velocity. Recently Sikalo has tackled this issue and proposed a novel method. His article is titled :"Dynamic contact angle of spreading droplets: Experiments and simulations" published in physics of fluids in 2005. You may find it interesting. Having said that, the question is now if it is possible to implement his methodology in FLUENT?!

 November 2, 2010, 03:45 dynamic contact angle ... cont #6 New Member   davide Join Date: Jul 2009 Posts: 2 Rep Power: 0 by experience ... if you like dyn. contact angle .... use immiscible lattice Boltzmann (ILB) approaches rather than traditional CFD ... ILB Methods are also 2 (o even more) order of magnitude faster than CFD on multiphase - multicomponent applications

 November 2, 2010, 10:41 #7 New Member   Concordia_CFD Join Date: Jul 2010 Location: Canada Posts: 24 Rep Power: 8 ILB is a very powerful method, the thing is that is not offered in FLUENT

 July 28, 2012, 03:24 dynamic contact angle #8 Member     Subhasish Mitra Join Date: Oct 2009 Location: Australia Posts: 56 Rep Power: 9 It's possible to capture the velocity at the interface in VOF simulation (FLUENT) and velocities thus obtained can be used to calculate the dynamic contact angle through an UDF. raj kumar saini likes this. __________________ SM

 March 5, 2014, 10:01 #9 New Member   cfd^2 Join Date: Mar 2013 Posts: 28 Rep Power: 5 Dears, I have a question on this issue. Isn't the static contact angle implemented in FLUENT suficient to simulating multiphase flow, since when one sets no slip BC at the wall the velocity is zero... I mean, if we take some empirical correlations to account for the dynamic advancing and receding contact angle we see that they depend upon the capillary number, which is a dimensionless group which takes the fluid velocity into account. So, if the velocity at the walls is zero, the capillary number also is zero and the dynamic advanced and receding contact angle is equal to the static advanced and receding contact angle. The matter now is that which we enter in FLUENT is, in fact, neither the advanced or receding contact angle, but the static contact angle. So, in my vision, which could be done to improve this implementation has to do with setting an advanced and a receding CA.

 March 5, 2014, 17:25 #10 Super Moderator   Glenn Horrocks Join Date: Mar 2009 Location: Sydney, Australia Posts: 12,832 Rep Power: 100 I do not understand your question. But if the wall/interface point is not moving you do not have the singularity and there is no problem - so of course it should work fine. It is only when the contact point is moving that there is a problem.

 March 7, 2014, 08:00 #11 New Member   cfd^2 Join Date: Mar 2013 Posts: 28 Rep Power: 5 In my case I have flow of the two phases and I have set no-slip BC at the wall and a static contact angle, since this is the only available option in fluent. Does anyone knows if the assumption of a static contact angle is feasible and which is the limit of validity of this assumption? If it is not valid, I think a dynamic contact angle should be implemented, considering a moving contact line. At this point we have the problem, right? The non-slip BC could not be used anymore...

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post hsieh Open Source Meshers: Gmsh, Netgen, CGNS, ... 32 September 13, 2011 05:50 sebastian_vogl OpenFOAM Running, Solving & CFD 3 June 22, 2009 12:25 Pulli FLUENT 0 March 1, 2007 13:31 Aireen FLUENT 1 August 10, 2006 16:01 Aireen FLUENT 2 July 5, 2006 13:14

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

 Contact Us - CFD Online - Top