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

Mass conservation problem with multiphaseInterFoam alphaContactAngle

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 14, 2010, 08:21
Default Mass conservation problem with multiphaseInterFoam alphaContactAngle
  #1
New Member
 
andres
Join Date: Aug 2009
Posts: 15
Rep Power: 8
andresbh is on a distinguished road
Hi
I am simulating with multiphaseInterFoam openfoam 1.6 a 2D sessile drop laying on a flat surface with a defined contact angle in alphaContactAngle.

During the simulation the mass of the phase laying on the surface starts do diminish until it completely disappears from the domain (the other way around for the other phase it start to take over the space left).

I have attached the test case to this posting and a plot of the volume of the phase laying on the surface.

Have I set up something incorrectly?

Thank you for your help

Andres
Attached Images
File Type: png plotMassInt.png (6.4 KB, 36 views)
Attached Files
File Type: gz sessileDrop.tar.gz (9.0 KB, 15 views)
andresbh is offline   Reply With Quote

Old   October 14, 2010, 08:24
Default Mass conservation problem with multiphaseInterFoam alphaContactAngle
  #2
New Member
 
andres
Join Date: Aug 2009
Posts: 15
Rep Power: 8
andresbh is on a distinguished road
A plot with longer time of the phase volume is attached to this post
Attached Images
File Type: png plotMassInt.png (6.3 KB, 21 views)
andresbh is offline   Reply With Quote

Old   October 25, 2010, 02:36
Default
  #3
Member
 
Richard Kenny
Join Date: Mar 2009
Posts: 59
Rep Power: 9
richpaj is on a distinguished road
Hello there,
I was intrigued by what you had to write as it overlaps with the "numerical evaporation" (!) I had observed some while ago with 3D meshes. As liquid droplets impact and spread over a surface they develop azimuthal kinks which
appear to become sites of rapid mass loss ("evaporation") when the lengthscale of the kink is similar to that of a characteristic mesh cell. Eventually the whole droplet will disappear. At that time we'd reached the limit of the resolution that could be modelled reasonably/efficiently by the computer concerned. Switching to a (high resolution) wedge/axisymmetric mesh avoided that issue owing to the suppression of the orthogonal
coordinate.

I was hopeful then (!) that your 2D problem might show similar anomalous behaviour
to that of our 3D case and which might be easier to examine (given the much smaller mesh size involved). Using OF16x's interFoam and swapping in the materials we happened to be considering over here, however, I found no "numerical evaporation" and indeed the average mass fraction remains close to the initial value.

(As for your particular setup I can only think that possibly a larger domain might help,
and/or certainly a comparison with interFoam.)

There is another factor in the above, namely that as the ratio of "interface area / initial droplet volume" increases then so does the error in the mass balance (or average gamma)
calculation. Additionally, the error increases with characteristic velocity. If these errors
are regarded as unacceptable (up to 10% in the cases we considered) then it might be
worthwhile switching codes to compressibleInterFoam.

Good luck,

Richard K.
richpaj is offline   Reply With Quote

Old   October 27, 2010, 12:15
Default
  #4
Member
 
Lars Kiewidt
Join Date: Sep 2009
Location: Germany
Posts: 54
Rep Power: 7
LarsPT is on a distinguished road
I think your problem occurs because of your pressure BC at the wall. Try this instead of zeroGradient at the walls!

<patch>
{
type buoyantPressure;
value uniform 0;
}

It should work!

Best regards, Lars!
LarsPT is offline   Reply With Quote

Old   November 2, 2010, 12:34
Default solved
  #5
New Member
 
andres
Join Date: Aug 2009
Posts: 15
Rep Power: 8
andresbh is on a distinguished road
Hi thanks for your explanations,

LarsPT's solution solved it. But I dont really get it. What does this BC do actually?

Cheers Andres
andresbh is offline   Reply With Quote

Old   November 2, 2010, 12:53
Default
  #6
Member
 
Lars Kiewidt
Join Date: Sep 2009
Location: Germany
Posts: 54
Rep Power: 7
LarsPT is on a distinguished road
Glad to hear that! It has something to do with the change from pd to p in OF-1.6 (OF-1.7 is solving for p_rgh again, which is the same as pd). There was a thread about it in this forum but I can"t find it right now. Nonetheless, zeroGradient for p results in mass conservation issues.

Lars
LarsPT 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
DPM mass problem Henrik Ström FLUENT 8 December 14, 2014 12:44
Need any WHATEVER Help,please. Heat-Flow problem Amr FLUENT 1 March 15, 2013 12:47
Star-CD transient simulation, problem of loss of mass kit STAR-CD 3 April 13, 2010 11:34
Question on Sediment Mass Conservation Scheme Wen Long Main CFD Forum 0 July 6, 2004 15:45
mass conservation in diverging flow yonghong yan Main CFD Forum 4 July 27, 2002 02:06


All times are GMT -4. The time now is 23:01.