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

How Can I calculate the Contact Line velocity in interface??

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 5, 2010, 12:48
Default How Can I calculate the Contact Line velocity in interface??
  #1
Member
 
Join Date: Nov 2009
Posts: 48
Rep Power: 7
farhagim is on a distinguished road
Hello,

Can anybody tell me How can I calculate the Contact line Velocity in interFoam? I want to implement new Dynamic Contact Angle, and I need Contact line velocity for Ca. I would be so grateful if the expertise in OpenFOAM help me.

Thanks
Mehran
farhagim is offline   Reply With Quote

Old   August 6, 2010, 03:24
Default
  #2
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Quote:
Originally Posted by farhagim View Post
Hello,

Can anybody tell me How can I calculate the Contact line Velocity in interFoam? I want to implement new Dynamic Contact Angle, and I need Contact line velocity for Ca. I would be so grateful if the expertise in OpenFOAM help me.

Thanks
Mehran
Quick brainstorming:

Find the face (and its coordinates) on the contact surface occupied by the interface at the current timestep (for example where VOF = 0.5)
The spatial difference to the previous timestep will be the displacement of the interface. In addition with the timestep you will get the displacement velocity - the contact line velocity.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   August 6, 2010, 12:20
Default
  #3
Member
 
Join Date: Nov 2009
Posts: 48
Rep Power: 7
farhagim is on a distinguished road
Hello Sebastion.

Thank you for your answer. but The contact line velocity is the intersection between liquid surface and the solid surface. So I think I should find the cell with alpha=0.5 near surface and find the velocity or find the displacement of the interface and so on. Which one do you think better? Actually I am kind of new in Openfoam and I have a problem to recognize the all variable defined in interfoam ie. nHat or AfnHat ...
some times its hard to find what they come from. Would you please give me some hint in code How can I define the interface?

Thanks,
MEhran

Quote:
Originally Posted by sega View Post
Quick brainstorming:

Find the face (and its coordinates) on the contact surface occupied by the interface at the current timestep (for example where VOF = 0.5)
The spatial difference to the previous timestep will be the displacement of the interface. In addition with the timestep you will get the displacement velocity - the contact line velocity.
farhagim is offline   Reply With Quote

Old   August 6, 2010, 12:47
Default
  #4
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Quote:
Originally Posted by farhagim View Post
So I think I should find the cell with alpha=0.5 near surface and find the velocity or find the displacement of the interface and so on. Which one do you think better?
Hm, interesting idea to interprete the velocity in the cell itself as the interface velocity. What do the other interFoam users think?
I think this may be appropriate at the proximity to the solid surface ... but I'm not 100% sure.

Quote:
Originally Posted by farhagim View Post
Actually I am kind of new in Openfoam and I have a problem to recognize the all variable defined in interfoam ie. nHat or AfnHat ...
some times its hard to find what they come from. Would you please give me some hint in code How can I define the interface?
In my option the interface is defined as alpha = 0.5.

The variables you mentioned have something to do with the normal direction of the interface.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   August 6, 2010, 13:14
Default
  #5
Member
 
Join Date: Nov 2009
Posts: 48
Rep Power: 7
farhagim is on a distinguished road
Thanks again for your quick answer. one very basic question: forexample: tell me what should I write in code for defining the cell filled with alpha=0.5. Sorry I know its stupid but I am new .

Quote:
Originally Posted by sega View Post
Hm, interesting idea to interprete the velocity in the cell itself as the interface velocity. What do the other interFoam users think?
I think this may be appropriate at the proximity to the solid surface ... but I'm not 100% sure.



In my option the interface is defined as alpha = 0.5.

The variables you mentioned have something to do with the normal direction of the interface.
farhagim is offline   Reply With Quote

Old   August 7, 2010, 04:01
Default
  #6
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 11
sega is on a distinguished road
Quote:
Originally Posted by farhagim View Post
Thanks again for your quick answer. one very basic question: forexample: tell me what should I write in code for defining the cell filled with alpha=0.5. Sorry I know its stupid but I am new .
Looping over all cells will be something like this:

Code:
forAll(cells, cellI){
}
cellI will then point to the current cell.

Inside the loop forAll-loop you will have to check for the interface.

Finding the interface with alpha = 0.5 may be a bit more difficult because the interface may not have an exact value of 0.5, but will be in between 0 and 1.

You should make a thread about this in the Programming section of this forum, because I'm not sure at the moment how to check for the alpha-value of the cell.
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Old   August 9, 2010, 11:55
Default
  #7
Member
 
Join Date: Nov 2009
Posts: 48
Rep Power: 7
farhagim is on a distinguished road
Thanks for your help Sebastion. Yeah I have to make a thread to ask about this too.

Quote:
Originally Posted by sega View Post
Looping over all cells will be something like this:

Code:
forAll(cells, cellI){
}
cellI will then point to the current cell.

Inside the loop forAll-loop you will have to check for the interface.

Finding the interface with alpha = 0.5 may be a bit more difficult because the interface may not have an exact value of 0.5, but will be in between 0 and 1.

You should make a thread about this in the Programming section of this forum, because I'm not sure at the moment how to check for the alpha-value of the cell.
farhagim is offline   Reply With Quote

Old   October 7, 2010, 09:49
Default
  #8
Member
 
Duong A. Hoang
Join Date: Apr 2009
Location: Delft, Netherlands
Posts: 92
Rep Power: 8
duongquaphim is on a distinguished road
Send a message via Yahoo to duongquaphim
Hi Farhagim,

I am also very interested in checking the velocity and location of the interface. I guess what we can do is getting the location and the value of alpha at every points then make the interpolation for alpha = 0.5 to get the location of the interface every time step. However, since I am a newer also, I do not actually know how to do it yet.

Did you find success in getting it?

Duong
duongquaphim is offline   Reply With Quote

Old   October 7, 2010, 19:49
Default
  #9
Member
 
Lars Kiewidt
Join Date: Sep 2009
Location: Germany
Posts: 54
Rep Power: 7
LarsPT is on a distinguished road
So, what I do to get the interface (and the positions of a moving bulk which are important for me) is just sampling the whole alpha1 field. Then I use MATLAB to find the values of alpha1 between 0.4 and 0.6 and finally get the interface. I think you can do this with a regular shell script, so you don't have to touch the OpenFOAM code.

Hope that helps!

Best regards,

Lars
LarsPT is offline   Reply With Quote

Old   October 8, 2010, 03:52
Default
  #10
Member
 
Duong A. Hoang
Join Date: Apr 2009
Location: Delft, Netherlands
Posts: 92
Rep Power: 8
duongquaphim is on a distinguished road
Send a message via Yahoo to duongquaphim
Hi Lars,

That is what I did at this moment. I wrote both Matlab script and Paraview python to get the position and velocity of the interface in post-processing. However, what I want to do now is getting the position of the interface during my simulation to implement my model.

Then I got troubles with how we define x, y, z, alpha of every grid point and the communication between cell and cell (for instance how can we know the value of the neighbor cell).

If you have any experiences about these, please let me know. Thanks.

Regards,

Duong
duongquaphim is offline   Reply With Quote

Old   May 10, 2011, 23:25
Default contact line velocity
  #11
New Member
 
pitambar randive
Join Date: May 2011
Posts: 1
Rep Power: 0
pitambar is on a distinguished road
dear

i want to know how contact line velocity is different from bulk velocity and which formula can be used

I am thinking that KWOK (2004) would be proper to use

I want your comments
pitambar 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
Conversion Error nuovodna Open Source Meshers: Gmsh, Netgen, CGNS, ... 14 October 1, 2010 11:07
Regarding FoamX running Kindly help out hariya03 OpenFOAM Pre-Processing 0 April 18, 2008 04:26
Problem with rhoSimpleFoam matteo_gautero OpenFOAM Running, Solving & CFD 0 February 28, 2008 07:51
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 18:59.