
[Sponsors] 
January 30, 2009, 15:09 
multiphase flow, quick divergence of contuinity eq

#1 
Guest
Posts: n/a

Hi all,
I am modeling a two phase flow, air and molten glass. temperature ~ 1600K. the viscosity ratio of the two phases are high ~ 10^5 and the density ration is also high ~ 10^3. The domain consists of one horizontal velocity inlet of pure molten glass, one vertical outlet (pointing downward) of molten glass. There is another upward vertical vent where the molten glass is in contact with air. I am using VOF model so that the interface of air and molten glass can be tracked. The flow domain is mainly molten glass. I've tried all the possible ways suggested in the User's Guide, like using "node based" gradient option, CICSAM scheme for the Volume Fraction Discretization, enabling "Specified operating density", reducing underrelaxation factor for density, using very small time steps, etc. However, the continuity equation diverges after several iterations. I am desperate now and any suggestions about making it converge would be helpful. Thanks a million. BTW, I used unstructured mesh. There are a few highly skewed meshes (>0.97) on the volume connecting the horizontal domain and the upwardfacing vent. The geometry are originally constructed in Solidworks and that is the best mesh I could get. I wonder if this affects the divergence greatly. 

January 31, 2009, 09:10 
Re: multiphase flow, quick divergence of contuinit

#2 
Guest
Posts: n/a

hello do you patch the initial condition?


January 31, 2009, 14:08 
Re: multiphase flow, quick divergence of contuinit

#3 
Guest
Posts: n/a

yes, I did. I patched a volume fraction of 1 (for air) to the area that is originally occupied by air.


February 1, 2009, 07:16 
Re: multiphase flow, quick divergence of contuinit

#4 
Guest
Posts: n/a

r u using presto scheme or not. if not first thing is to switch it to presto. do it pronto.


February 1, 2009, 17:25 
Re: multiphase flow, quick divergence of contuinit

#5 
Guest
Posts: n/a

just tried using "Presto!" for pressure, diverged in 7 iterations (Error: > (greaterthan): invalid argument [2]: wrong type [not a number] Error Object: nan), though it does seem like it is the "velocity" fields that are diverging.


February 2, 2009, 22:42 
Re: multiphase flow, quick divergence of contuinit

#6 
Guest
Posts: n/a

Try putting some boundary layer mesh at the interface between the two fluids. this will make it easier for the interface tracking algorithm. Having highly skewed mesh will in fact affect convergence, but it shouldn't be that significant.
For the specifying of operating density, please specify the lighter of the two fluids, and also please specify a reference point that is explicitly inside the lighter phased fluid. To be honest, I faced exactly the same problem with my simulation. All the momentum and volume fraction equations were converging, but the continuity equation diverged, and very rapidly I have to say. The funny thing is, (since I have access to both FLUENT and CFX) I tried the VOF simulation in CFX, and it converges very nicely indeed (I have now ditched FLUENT for CFX for all my VOF simulation). So I'm suspecting that the segregated solver inside FLUENT may not be up to the task for VOF, as CFX is inherently a coupled solver (if there are any multiphase flow developers from FLUENT reading this, please verify this please!). Although I haven't tried this myself, try using one of the density based coupled solver and see if that works. Do let me know how it goes. Hope that helps. 

February 2, 2009, 23:39 
Re: multiphase flow, quick divergence of contuinit

#7 
Guest
Posts: n/a

Thanks very much for you detailed response.
1. I am not sure how to do this. The interface is not known a prior, at least I do not know. It should be somewhere in the upward pointing pipe, but I do not know how much the molten glass will fill the pipe until I do have a sucessful simulation. I will try to ask someone else to see if they have information regarding this. 2. Yes I used the lighter fluid density (air) for the operating density. I have pressure boundary conditions at the vent, so according to the User's guide, the reference point for pressure should be useless. 3. If this is the case, I'll have to stay desperate since I do not have CFX. Thanks again! 

February 3, 2009, 18:54 
Re: multiphase flow, quick divergence of contuinit

#8 
Guest
Posts: n/a

could you try a simple test for me on cfx, it is useful for us.
along with air try a fluid with very high viscosity say : 8E4 and density say 1200 and let me know if CFX can handle it. Fluent can only handle this with presto scheme, and it is really tough to do. I wish to see how cfx does when viscosity ratio is very high. Thanks. 

February 16, 2016, 06:32 

#9 
New Member
Join Date: Jul 2015
Posts: 3
Rep Power: 3 
@violet :Were you able to figure it out later? I'm stuck with a similar case right now.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Nozzle Flow  Divergence Problem  Ijaz  FLUENT  9  January 11, 2014 05:36 
flow past a missle: how to solve divergence?  xiaofish  FLUENT  0  September 9, 2007 22:53 
strange divergence when solving multiphase problem  tanghao  FLUENT  2  July 27, 2006 19:47 
Divergence of Laminar Flow With Eddy ???  Peter Kostka  Main CFD Forum  3  May 9, 2002 17:02 
How to avoid Divergence of multiphase (gas/liquid system) in Fluent :???  mounir  Main CFD Forum  1  September 1, 1999 16:30 