CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Simulating liquid droplet deformation

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

LinkBack Thread Tools Display Modes
Old   June 13, 2013, 08:50
Default Simulating liquid droplet deformation
New Member
Join Date: Apr 2013
Posts: 24
Rep Power: 4
mebinitap is on a distinguished road
Dear all,

I want to see the deformations in a liquid droplet in presence of air flow. So I used a squared domain (in blockMeshDict file) and setFieldsDict for the droplet. I stated the boundaries as atmospheric and used the freeStream BC conditions for P and U. Problem is that droplet moves away leaving the domain if I put some value to velocity. Since I want to track the shape changes I have to move along with the droplet. Can anyone please give me a hint on that? Also are the BCs correct ?

I want to model the shape changes in a static 3D droplet placed in a convective flow field and then to see how it response in presence an acoustic wave (acoustic levitation).
mebinitap is offline   Reply With Quote

Old   June 14, 2013, 05:30
Join Date: Oct 2010
Location: Delft, Netherlands
Posts: 97
Rep Power: 6
michielm is on a distinguished road
Hi mebinitap,
A note upfront: the deformation of a static and a dynamic droplet is not the same so you cannot use the deformation of the moving droplet to study your problem. Therefore, moving along with the droplet is not going to help you.

A droplet in an air flow can only remain static if the surface shows contact angle hysteresis. Experimentally, this is pretty much always the case. Numerically however, by default there will be no contact angle hysteresis which will cause the droplet to move already for the tiniest amount of air flow.
Implementing contact angle hysteresis is not that simple because you need to check whether the contact angle is in a certain range and if so, adapt the local contact angle instead of moving the contact line.

A trick you could give a try is to adjust your boundary condition with funkySetBoundaryField. You can set the contact angle to a certain value underneath the droplet and to a higher value outside of the droplet. This will cause the droplet to stick to the spot it is located and effectively generates contact angle hysteresis.
michielm is offline   Reply With Quote

Old   June 14, 2013, 08:45
New Member
Join Date: Apr 2013
Posts: 24
Rep Power: 4
mebinitap is on a distinguished road
Thank you so much for your suggestion. I will try to implement it and let you know.
mebinitap is offline   Reply With Quote

Old   June 27, 2013, 00:39
New Member
Join Date: Apr 2013
Posts: 24
Rep Power: 4
mebinitap is on a distinguished road
Hi again,

I am a bit confused about surface tension force implementation in interFoam. It seems it is handled as : fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1). But if a drop is sitting on vibrating surface and contact angle is not in equilibrium, then how are the diff sigmas at all surfaces (solid-liq, liq-gas, gas-solid) are taken care of, i mean there should be some kind of Young's type of eqn. But I dont realize how is that implemented. Should that be given as some dynamic contact angle boundary condition ? It will be great if I can get any kind of hint.

mebinitap is offline   Reply With Quote


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
Modeling Liquid Liquid droplet flow Mircro_fluidics_ STAR-CCM+ 6 July 20, 2015 21:43
Problem of simulating of small droplet with radius of 2mm liguifan OpenFOAM Running, Solving & CFD 5 June 3, 2014 02:53
UDF for droplet of immiscible liquid mixing boh_tea FLUENT 0 March 7, 2011 22:43
Liquid Droplet Impingement naimishharpal STAR-CCM+ 0 February 28, 2010 14:24
Liquid fueled Droplet Evaporation Ray FLUENT 8 January 3, 2007 01:36

All times are GMT -4. The time now is 17:24.