CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum


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

LinkBack Thread Tools Display Modes
Old   June 4, 1999, 02:10
Chuck Leakeas
Posts: n/a
I have a problem where I have two stationary free surfaces. I use the continuity eqn. to compute what the velocity should be on the free surfaces every iteration. The problem is that I have alot of convergence problems with this case. Does anyone have any experience with these types of problems or any suggestions on how to rid myself of these convergence issues? I am using a SIMPLE like algorithm and solving for u,v, and p in a 2D axisymmetric geometry. Any help would be greatly appreciated.
  Reply With Quote

Old   June 8, 1999, 11:10
Neyval C. Reis
Posts: n/a
Hi Chuck

I don't know what exactly your problem is, there are certain bits and pieces of the treatment of the free-surface that are always troublesome. The major problem is usually related to the implementation of the mass conservation in the cells containg the free-surface.

First of all. What are you doing with the pressure boundary condition at the interface ? If you are using Simple I supose you are prescribing the value of P at the free surface and setting the pressure correction as zero at the cell containg the interface. Is that correct ? If so you MUST make sure that the implementation of the continuity eq. in theses cell is correct otherwise you end up loosing mass or having convergence problems.

When you are using the continuity eq. to ensure mass conservation in the cell containing the free surface you have to be be aware that somethimes a cell can have 2 empty neighbours cells. In this case when you use the continuity eq. you end up with two unknows instead of one. What usually is done is to consider rho.Ui.Ai=rho.Ui+1.Ai+1 (where Ai is the area tranversal to i) for each direction to solve the problem (the original report of the MAC method - LA3425, Los Alamos, 1966 - contains a nice explanation of that, as well as the implications of using this aproximation).

I don't know if this can help you.



  Reply With Quote

Old   June 8, 1999, 13:11
Chuck Leakeas
Posts: n/a
Neyval, thanks for the response to my problem. I use a variation of the SIMPLE method called SIMPLE-M, M means Modified for a collocated grid structure. I do not solve a pressure correction eqn., only the pressure field from the current velocity field. I use the continuity eqn to calculate the velocity on the free surface. I do this every iteration and the convergence is very slow. If I do not include the buoyancy force in the problem the convergence is quick and the solution looks just fine. It is the coupling of the momentum and energy eqns through the buoyancy term that gives me all the trouble. Any suggestions? Any help would be greatly appreciated. Thanks, Chuck
  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
Modelling free surfaces using no slip wall, and other general queries andy_mckay88 Main CFD Forum 1 November 18, 2010 06:54
Can Flow-3D plot the free surface area in Iso-surface or colour variable? therockyy FLOW-3D 1 June 20, 2010 19:36
free surfaces and two phase Arnie Phoenics 6 May 20, 2005 07:20
Modelling Free Surfaces Using Fluent 6 M Fortey FLUENT 3 October 4, 2003 09:09
CFX 4.4 New free surface option Viatcheslav Anissimov CFX 0 April 3, 2002 06:27

All times are GMT -4. The time now is 15:29.