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

InterFoam with lid driven cavity

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 10, 2005, 21:50
Default I tried to produce a lid drive
  #1
Senior Member
 
Billy
Join Date: Mar 2009
Posts: 167
Rep Power: 8
billy is on a distinguished road
I tried to produce a lid driven cavity flow using interFoam. I set initial values of gamma in some cells to 1.0 and the rest to 0.0. I wanted to simulate a mixing of two fluids, however I get always the same message even using smaller time steps:

...
Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Calculating field g.h

time step continuity errors : sum local = 0, global = 0, cumulative = 0

Starting time loop


Max Courant Number = -0
Time = 0.001

Liquid phase volume fraction = 0 Min(gamma) = 0 Max(gamma) = 0
Liquid phase volume fraction = 0 Min(gamma) = 0 Max(gamma) = 0
Liquid phase volume fraction = 0 Min(gamma) = 0 Max(gamma) = 0
Liquid phase volume fraction = 0 Min(gamma) = 0 Max(gamma) = 0
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ExecutionTime = 0.12 s

...

Does anyone know what might be wrong?
billy is offline   Reply With Quote

Old   July 10, 2005, 22:22
Default Actually the error looks like
  #2
Senior Member
 
Billy
Join Date: Mar 2009
Posts: 167
Rep Power: 8
billy is on a distinguished road
Actually the error looks like this:

...

Max Courant Number = -0
deltaT = 0.001
Time = 0.001

Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ExecutionTime = 0.1 s



Max Courant Number = -0
deltaT = 0.001
Time = 0.002

Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
Liquid phase volume fraction = 0.0625 Min(gamma) = 0 Max(gamma) = 1
time step continuity errors : sum local = nan, global = nan, cumulative = nan
ExecutionTime = 0.22 s

...
billy is offline   Reply With Quote

Old   July 11, 2005, 06:40
Default You seem to get a division by
  #3
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
You seem to get a division by 0.

Also you have no flow? Nothing is solved for it seems? Do you have valid boundary conditions? Valid initial conditions?

If you are interested in where it goes wrong: set the environment variable FOAM_SIGFPE to 1 which will make the code coredump and then look at the core in a debugger.
mattijs is offline   Reply With Quote

Old   July 11, 2005, 14:37
Default I don't know what happened, bu
  #4
Senior Member
 
Billy
Join Date: Mar 2009
Posts: 167
Rep Power: 8
billy is on a distinguished road
I don't know what happened, but I know it was due to the mesh file. I created one using blockMesh and it worked OK. The original mesh file seemed OK when viewing it in paraview though. Does OpenFOAM automatically orient the face normals?
billy is offline   Reply With Quote

Old   July 11, 2005, 17:50
Default Yes, all boundary faces point
  #5
Super Moderator
 
Mattijs Janssens
Join Date: Mar 2009
Posts: 1,416
Rep Power: 16
mattijs is on a distinguished road
Yes, all boundary faces point away from the owner cell. There is a section in the user guide about the mesh structure.
mattijs is offline   Reply With Quote

Old   July 13, 2005, 15:05
Default Hi, After I got this case w
  #6
Senior Member
 
Billy
Join Date: Mar 2009
Posts: 167
Rep Power: 8
billy is on a distinguished road
Hi,

After I got this case working (Lid-driven cavity flow with two phases) I noticed that altering the surface tension produces little effect on the output. I tried positive, negative and zero values of the surface tension and there seems to be no mix between the two fluids. Also, both fluids have same material properties such as density and viscosity.

I assume that:

A postive surface tension -> immiscible
A negative surface tension -> miscible

Is this right?

PS: I wish to thank the OpenFOAM team for making available such a great software.
billy is offline   Reply With Quote

Old   July 13, 2005, 15:44
Default interFoam uses a VOF-type meth
  #7
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 13
henry is on a distinguished road
interFoam uses a VOF-type method to capture the interface which of course does not allow mixing of the phases, that is the point. Negative surface tension has no physical meaning and may very well destabilise the solution algorithm.
henry 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
driven cavity neelam Main CFD Forum 0 June 13, 2007 00:36
Lid Driven Cavity 3D Mathias Krause Main CFD Forum 5 October 17, 2006 06:34
3D Lid Driven Cavity Danny Main CFD Forum 0 October 22, 2002 03:09
2D Lid-Driven cavity Sophie Main CFD Forum 8 April 28, 2002 18:34
Lid-Driven Cavity Paul Safier Main CFD Forum 4 February 19, 2002 07:20


All times are GMT -4. The time now is 03:33.