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

Wrong interface curvature capturing and contact line calculation at low capillary nu

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

LinkBack Thread Tools Search this Thread Display Modes
Old   April 14, 2015, 14:04
Default Wrong interface curvature capturing and contact line calculation at low capillary nu
Mahmoud Aboukhedr
Join Date: Feb 2014
Location: London
Posts: 37
Rep Power: 9
Mahmoud_aboukhedr is on a distinguished road
Dear Foamers, Hopefully all are doing well with simulations,
Am facing a big problem using interFoam to simulate oil slug mobilization in capillary tube filled with water for very low capillary numbers (range of 1e-8), I will try to explain the case as clear as I could so I can get the proper help,
The experimental paper I have is describing the oil mobilization on middle section of a capillary tube of dimension (0.5mm*400mm*0.5mm), and the imposed flow rate is in the range of (1e-10 t0 1e-13 m^s), for which it represent a very low capillary number. The experimental paper shows the development of the oil slug shapes as shown int the attached file (1A).

Here is the important part, I am trying to simulate the same experimental case in my PhD, but am facing a big failure.
My simulation conditions are as following,
- Imposed volume flow rate at the inlet and atmospheric pressure at the outlet.
- Constant contact angle at the walls with value of theta=145
- The slug is initialized as box in the begging of the tube as shown in (2A)

The slug never deform like the experiment at low flow rate, at low flow rate the slug take the droplet shape at the first time steps and then start to move along the tube, while if I increase the flow rate the trailing edge of slug deform and separate see picture 3A for comparisons

The question
1) Why at low flow rate the interface curvature is wrong, and why at high flow rate it deform and does not limited by the constant contact angle.
2) Why the contact line at low flow rate is wrong? also I tried the dynamic contact angle but it does not work too.
3) How can I correct for the interface curvature in interFoam for viscous driven flow?

I Am sorry for asking to many questions, (maybe they are stupid) but I will appreciate any help on that matter.
Attached Images
File Type: jpg 1A.jpg (35.1 KB, 26 views)
File Type: jpg 3A.jpg (20.7 KB, 33 views)
File Type: jpg 2A.jpg (4.7 KB, 18 views)
Mahmoud_aboukhedr is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

All times are GMT -4. The time now is 18:27.