CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM Running, Solving & CFD (
-   -   contact angle behavior in micro & nano dimensions (

Jimbomet November 29, 2011 08:11

contact angle behavior in micro & nano dimensions
Hi all foamers,
I am a new entry in the world of OF and I'm trying to simulate the formation of a liquid meniscus between two flat surfaces (2D) as a process of merging of two drops.
I'm using InterFoam and I initialize two liquid half drops that intersect themselves, with their base on the two surfaces.

My domain dimension is very small (order of nanometers), and, working with these dimensions I encountered a strange behavior in the contact angle, in the sense it not respects my BC at the wall.
Dealing with micro dimensions I obtain a meniscus (but without the correct angle, that I imposed null), but when I try with nano dimensions the meniscus becomes a "cylinder" with infinite curvature (as if I Imposed alpha1=90)
Could someone help me?Any suggestion to solve this behaviour?... I can't understand what could be wrong in my case..?
Thanks in advance


duongquaphim December 2, 2011 09:06


I won't expect any difference between these two scales in the simulation if your mesh is fine enough. If possible, please post some figures and your case then we can talk more on that.



Jimbomet December 2, 2011 10:24

case description
2 Attachment(s)
Thanks for your reply Duong,
I have not solved my problem yet, I tried with a 3D simulation without changes in my behaviour...
My case consists in two flat plate at certain distance, and I want to simulate the formation of a liquid meniscus between them... : I initialise (with funckySetFields) a certain part of the domain with liquid phase (in particular two semi-sphere that intersect themselves, attached to the two surfaces), and I evaluate the equilibrium solution for each aim is to calculate the variation of pressure between wall with distance (but this is an other topic...).
My liquid is ethanol and the other phase is air. I'm using InterFoam solver. There are naturally no gravity forces in my problem.
I imposed a null contact angle at the walls and buoyant pressure.
In the 2D case I found the difference shown in the attached figures, when dealing with micro and nano scale

I thinks that behaviour is due to the very high drop of pressure at the interface of the two fluids at the interface (at nanoscale more than microscale), and perhaps the interface-compression can't correctly manage it, but I am not an expert at all... What could I do to solve that behavior, or it is the correct one and there is nothing to solve??
Thanks in advance

I forgot to give the dimension of the domain: respectively 3nm and 7.5nm for the two sides of the rectangular domain, and the same measures but in micron for the "bigger" case

duongquaphim December 2, 2011 11:27

Hi Paolo,

As I understood, in nano case, you initialized two semi-sphere and after getting to steady state, the interface become flat. It is a little bit strange. Are you using constantAlphaContactAngle or dynamicsAlphaContactAngle bc?

Also a good check might be a couple of simulations with different width from micro to nano scale (let's say 10micro, 1micro, 100nano and 10nano) to see when you observe this behavior.



Jimbomet December 2, 2011 11:53

2 Attachment(s)
I know it is a bit strange. I attach the first time step and the last one for the nano problem.
I use constant contact angle. my condition is:

type constantAlphaContactAngle;
gradient uniform 0;
limit none;
theta0 0;
value uniform 0;

Moreover, it is normal that when I initialise with funckySetFields only one of my wall have a
value nonuniform List<scalar>

I don't know where I'm wrong because performing a simulation in a bigger scale (say mm) it seems to change with respect to micro one (the contact angle seems to become smaller, more similar to null angle ).

Thanks for your suggestion, I'm going to try with dimensions between the two scales (I have only to find the correct time scales for each case...10^-9 for the nano-dimensions).
Thank you again, and good weekend

Jimbomet December 5, 2011 10:16

5 Attachment(s)
I've tried, as you suggested me, to run some simulation at different scale between nano and micro dimensions. I attach what I found in the range from mm to nm.
In the dimension of mm the BC contact angle seems to be respected while at smaller scales it seems not.
I tried to lower the viscosity in one case (100nm) case without apparent improvements for the moment...
Any suggestion?


duongquaphim December 6, 2011 11:06

Hi Paolo,

From your simulations, contact angle boundary condition worked pretty well in the scale of mm and micrometer but not in the scale of nanometer. And it is understandable since VOF is built on continuum mechanic. If you go to nano-scale, Knudsen number go to 1 and then molecular force becomes important. In that situation, you might want to use MD simulation rather than this VOF which is a continuum-based method.

I think that is the explanation for the outcome of your simulation. And I think I was wrong previously when saying that "I won't expect any difference between these two scales".



All times are GMT -4. The time now is 13:48.