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

Interfoam Droplet under shear test case

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 25, 2009, 22:52
Default Hi Everyone, I am in the pr
  #1
Member
 
Adam Donaldson
Join Date: Mar 2009
Location: Ottawa, Ontario, Canada
Posts: 37
Rep Power: 8
adona058 is on a distinguished road
Hi Everyone,

I am in the process of validating an alternative multiphase solver, and one of my test cases has been giving me some problems.

To determine if the issue was related to my solver, or the case, I have tried solving it using interFoam (and interDyMFoam). Alas, I obtain the same incorrect results. I was hoping someone might be able to shed some light on the problem.

In the classical droplet under shear case, a Capillary number of 0.1 should result in a deformation parameter, D=(L-S)/(L+S) = 0.109375, where L and S are the longest and shortest radii in the shear plane. Interfoam and my new solver both give a value of D = 1.68 (give or take). This is significantly higher than expected. When I try the case with a higher capillary number (0.3), the droplet doesn't reach a steady state (It eventually elongates and breaks)

I am working on a 3D grid, 128x64x128, with dimensions of [-4 to 4],[-2 to 2], [-4 to 4]. A droplet with radius = 1 is placed in the center of the domain. The top and bottom boundaries are walls, where the top has an imposed velocity of +4 in the x direction, while the bottom has a velocity of -4 in the negative x direction. The remainder of the velocity field is initialized with Ux = Z (i.e. a shear rate of 1).

The two boundaries on the left and right are cyclic, and the front and back are set to zero-gradient for all variables.

Both fluids have a density of 4 kg/m3, and a viscosity of 1 kg/ms (0.25 m2/s). The surface tension is = 10 kg/s2. Under these conditions, the Capillary number for the droplet is = 0.1.


If anyone has successfully completed this test case using interFoam, I would love to hear from you with details on how you set the case up, and what you think might be causing the shear stress to be poorly accounted for.
adona058 is offline   Reply With Quote

Old   July 10, 2009, 15:54
Default
  #2
New Member
 
sesha
Join Date: Mar 2009
Posts: 23
Rep Power: 8
hamsadhwani8 is on a distinguished road
Hi,

I am trying to solve the exact problem. I am not sure if you have met with any success. I was able to get something done for a 2D case with some accuracy, i.e., 1 cell in the z direction. I also found this over estimation issue. For a 2D it is not that bad. I am trying to set up a 3D case but am getting into problems while running it. Check my other latest post today.
hamsadhwani8 is offline   Reply With Quote

Old   July 10, 2009, 15:57
Default
  #3
New Member
 
sesha
Join Date: Mar 2009
Posts: 23
Rep Power: 8
hamsadhwani8 is on a distinguished road
Adona,

You should try the inflow-outflow condition on the left and right face instead of the cyclic boundary conditions. This might be the problem.
hamsadhwani8 is offline   Reply With Quote

Old   May 3, 2010, 18:46
Default
  #4
Member
 
Join Date: Nov 2009
Posts: 48
Rep Power: 7
farhagim is on a distinguished road
Hello sir,

I am a new user in openfoam.. I want to simulate the droplet which falls on surface ( solid or liquid ). Could you please help me and tell me how do i have to define a droplet in openFoam and which solver do i have to use...

I would also be so grateful if you can help me simulate the one that u done.

thanks

Mehran

Quote:
Originally Posted by adona058 View Post
Hi Everyone,

I am in the process of validating an alternative multiphase solver, and one of my test cases has been giving me some problems.

To determine if the issue was related to my solver, or the case, I have tried solving it using interFoam (and interDyMFoam). Alas, I obtain the same incorrect results. I was hoping someone might be able to shed some light on the problem.

In the classical droplet under shear case, a Capillary number of 0.1 should result in a deformation parameter, D=(L-S)/(L+S) = 0.109375, where L and S are the longest and shortest radii in the shear plane. Interfoam and my new solver both give a value of D = 1.68 (give or take). This is significantly higher than expected. When I try the case with a higher capillary number (0.3), the droplet doesn't reach a steady state (It eventually elongates and breaks)

I am working on a 3D grid, 128x64x128, with dimensions of [-4 to 4],[-2 to 2], [-4 to 4]. A droplet with radius = 1 is placed in the center of the domain. The top and bottom boundaries are walls, where the top has an imposed velocity of +4 in the x direction, while the bottom has a velocity of -4 in the negative x direction. The remainder of the velocity field is initialized with Ux = Z (i.e. a shear rate of 1).

The two boundaries on the left and right are cyclic, and the front and back are set to zero-gradient for all variables.

Both fluids have a density of 4 kg/m3, and a viscosity of 1 kg/ms (0.25 m2/s). The surface tension is = 10 kg/s2. Under these conditions, the Capillary number for the droplet is = 0.1.


If anyone has successfully completed this test case using interFoam, I would love to hear from you with details on how you set the case up, and what you think might be causing the shear stress to be poorly accounted for.
farhagim 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
Compute shear stress interFoam hsieh OpenFOAM Post-Processing 42 June 17, 2015 07:39
3d test case Hassan Raiesi Main CFD Forum 1 August 19, 2006 12:33
Problems starting new case interFoam billy OpenFOAM Running, Solving & CFD 3 June 21, 2006 10:18
Test Case ganesh Main CFD Forum 0 March 16, 2006 13:34
My 3D interFoam test zou_mo OpenFOAM Running, Solving & CFD 4 July 19, 2005 09:20


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