CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

How to define the interface between two fluids

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 25, 2009, 08:10
Default How to define the interface between two fluids
  #1
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
Hello,

I am trying to model a tank which is half filled. Due to higher density water is is the fluid at the bottom part of the tank and air is above the water. I am simulating natural convection.

What I want to do is that at the interface the velocities of the air and water are the same in the plane of the interface but cannot penetrate one fluid in the other (i.e. normal velocity is zero). Also I want to set continuity of the heat transfer between the two fluids.

How should I model it? I did the mesh and I asig different fluids to each part, but I do not know how to model the interface. I tried with "interior" but it does not seem to work ( I get divergence and very strange patterns at the interface).

I would appreciate if somebody could help me with this,

thanks in advance
TfG is offline   Reply With Quote

Old   April 25, 2009, 10:36
Default
  #2
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,297
Rep Power: 41
-mAx- will become famous soon enough
Create two separated fluid volumes (connected), don't set any BC at the interface.
Once your model is set, then you initialize the whole domain, and you patch each one with the right vof-value.
Regarding the non-penetration, I think it should be done in the phase panel (interaction), but I am not sure
__________________
In memory of my friend Hervé: CFD engineer & freerider

Last edited by -mAx-; April 25, 2009 at 11:56.
-mAx- is offline   Reply With Quote

Old   April 25, 2009, 11:13
Default take care, that between both fluids is only one face!!
  #3
Member
 
Ralf Schmidt
Join Date: Mar 2009
Location: Austria
Posts: 67
Rep Power: 17
Ralf Schmidt is on a distinguished road
If you have two volumes, that are simply attached to each other, and you define (one) interface at interior (or whatever else) Fluent will create a wall for that other interface!!

You can assure that by using in the face command field the "connect faces" button (looks like a plug). Select (all) faces of your domain and connect. That will cause all superposed faces to be joined.

Ralf
__________________
CFD - nothing but Colourful Fluid Dynamics
Ralf Schmidt is offline   Reply With Quote

Old   April 25, 2009, 11:18
Default
  #4
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,297
Rep Power: 41
-mAx- will become famous soon enough
Quote:
Originally Posted by Ralf Schmidt View Post
If you have two volumes, that are simply attached to each other, and you define (one) interface at interior (or whatever else) Fluent will create a wall for that other interface!!
Ralf
If you split one volume into two volumes, the "interface" between the 2 volumes (connected) will be set automatically as interior (or internal).
It will set as wall, if the 2 volumes aren't connected (eg: 2 surfaces superposed)
__________________
In memory of my friend Hervé: CFD engineer & freerider

Last edited by -mAx-; April 25, 2009 at 11:55.
-mAx- is offline   Reply With Quote

Old   April 25, 2009, 12:18
Default
  #5
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
That is what I did, dont set any boundary at the interface (i.e. the same as to set interior condition). And yes, the faces were connected.

Are you suggesting me that I should use VOF model? Isn't any possibilitie to do it without it?

thanks
TfG is offline   Reply With Quote

Old   April 25, 2009, 13:01
Default
  #6
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,297
Rep Power: 41
-mAx- will become famous soon enough
I don't see other possibility else than multiphase (but I am not expert).
I did a similar calculation with a tank filled at the 2/3 with oil, and the rest with air, but without convection.
I solved it with multiphase
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   April 25, 2009, 15:23
Default
  #7
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
Tried the VOF model and patch like you said but I get and error at the beginning of the second iteration:
Updating solution at time level N...
Global Courant Number : 560.10

Error: Global courant number is greater than 250.0. The
velocity field is probably diverging. Please check the solution
and reduce the time-step if necessary.
Error Object: ()

The no penetration thing I didn`t find anywhere....
TfG is offline   Reply With Quote

Old   April 25, 2009, 16:58
Default
  #8
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,297
Rep Power: 41
-mAx- will become famous soon enough
Define/Phases...
Go to the interaction panel
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   April 26, 2009, 06:45
Default
  #9
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
I already did it, and defined my two phases (oil and air)... but, then in interaction there are only available two tabs to edit: Mass and Surface tension.... and nothing about no penetration.
I am pretty lost now... and I need this thing for my final project!
TfG is offline   Reply With Quote

Old   April 26, 2009, 06:51
Default
  #10
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,297
Rep Power: 41
-mAx- will become famous soon enough
on the 6.3 version there are more options like slip velocity between phases etc...
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   April 26, 2009, 07:55
Default
  #11
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
Yes yes, Im running 6.3.26 and I have those tabs, but it says: this page is not applicable under current settings...
TfG is offline   Reply With Quote

Old   April 26, 2009, 08:01
Default
  #12
Super Moderator
 
-mAx-'s Avatar
 
Maxime Perelli
Join Date: Mar 2009
Location: Switzerland
Posts: 3,297
Rep Power: 41
-mAx- will become famous soon enough
Then check the help, maybe those tabs aren't available with multiphase.
Try mixture model.
But as I said, I am not expert
__________________
In memory of my friend Hervé: CFD engineer & freerider
-mAx- is offline   Reply With Quote

Old   April 27, 2009, 04:12
Default
  #13
New Member
 
Nitin Suryawanshi
Join Date: Mar 2009
Location: Pune, India
Posts: 28
Rep Power: 17
suryawanshi_nitin is on a distinguished road
Dear friend
if u want to see only convection then make that interface to wall type in preprocessor and then import those two seperate volumes to fluen and to that wall & its shadow give 0 shear stress in x y & z direction so that proper convection pattern will form at interface (i.e wall) & model density with boussinesq model

regards Nitin
suryawanshi_nitin is offline   Reply With Quote

Old   April 27, 2009, 05:31
Default
  #14
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
Hi,

Thanks for the answers, I will try that one.

Only two observations, first, on thermal boundaries I shoul set coupled right?
Second, I see shear stress 0 could be a reasonable assumption for the part of water. But for the air, shouldn't it be more like a non-slip condition (the air will have it difficult to "move" the water (forgetting waves...))?

Just wondering about it, but I am not expert so I may follow your advice and set shear 0 for both sides...
TfG is offline   Reply With Quote

Old   April 29, 2009, 08:12
Default
  #15
New Member
 
Nitin Suryawanshi
Join Date: Mar 2009
Location: Pune, India
Posts: 28
Rep Power: 17
suryawanshi_nitin is on a distinguished road
yes u r right but this will give u proper convection velocity pattern for heat transfer near the inerface (avoiding no slip on both sides of interface wall.
have u tried this??

Thank you
suryawanshi_nitin is offline   Reply With Quote

Old   April 29, 2009, 14:29
Default
  #16
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
Hey,
Yes I tried with slip condition in both side and it seems to work fine, Thanks!

Still, do you think I should mantain lip conditions in the air part or no slip conditions?
TfG is offline   Reply With Quote

Old   April 30, 2009, 02:02
Default
  #17
New Member
 
Nitin Suryawanshi
Join Date: Mar 2009
Location: Pune, India
Posts: 28
Rep Power: 17
suryawanshi_nitin is on a distinguished road
with slip on air side will disturb ur convectopn (velocity pattern on air side )
if u want to see this just ceck it out..

Thanks
suryawanshi_nitin is offline   Reply With Quote

Old   May 3, 2009, 11:49
Default
  #18
TfG
New Member
 
Join Date: Mar 2009
Posts: 17
Rep Power: 16
TfG is on a distinguished road
I tried setting slip for water side and no-slip for air side and the heat transfer in the interface is smaller than if I set slip for both sides (about 20% less).

The question is which result is more "real"... ?
TfG is offline   Reply With Quote

Reply

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Missing math.h header Travis FLUENT 4 January 15, 2009 12:48
REAL GAS UDF brian FLUENT 6 September 11, 2006 09:23
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 13:24
UDF FOR UNSTEADY TIME STEP mayur FLUENT 3 August 9, 2006 11:19
Replace periodic by inlet-outlet pair lego CFX 3 November 5, 2002 21:09


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