# Two-phase problem (Vof?)

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

 July 25, 2004, 16:15 Two-phase problem (Vof?) #1 Raj Guest   Posts: n/a I need help in solving a two phase problem. I have a closed tank(3-d, complicated geometry)that is half filled with liquid and has ullage space(with air) above it. After 4 hours, it is drained. I am using CFD-ACE to model flow and heat transfer. What is the most appropriate way to solve this? I tried to use VOF, but it requires a very small time step which is not feasible for my problem. Can anyone suggest what to do here?

 August 3, 2004, 09:31 Re: Two-phase problem (Vof?) #2 Eugene de Villiers Guest   Posts: n/a Cant you increase or add VOF subcycles to reduce the CFL limitations of the overall timestep? In a good VOF code you should be able to run the interface at a courant number of near one, provided you subcycle the phase fraction transport equation about 4 or 5 times. Let me know if you need a reference. Eugene

 August 3, 2004, 10:40 Re: Two-phase problem (Vof?) #3 Raj Guest   Posts: n/a Thanks Eugene - Yes, can you please send me the references. Also, I dont understand what do u mean by adding/increasing VOF sub-cycles. I am using CFD-ACE.

 August 3, 2004, 13:11 Re: Two-phase problem (Vof?) #4 Eugene Guest   Posts: n/a Well, I was using this code: www.nabla.co.uk for my VOF calculations of spray atomization, there are some descriptions on the website. I do not have any knowledge of CFD-ACE unfortunately. By subcycling I mean you run the momentum-pressure solver at one timestep. Then after the PISO loop has converged for that time, you solve the phase fraction transport equation at a different timestep that is an integer fraction of the main timestep. Example: main timestep = 1 second. Using the converged velocity field, solve the phase fraction transport with a timestep of say 0.25 seconds. By doing this 4 times you arrive at a phase fraction distribution at the same time as the current momentum solution. You can make it more accurate by storing old fields and using time interpolated velocity at the different subcycles. Since your fluid fraction is being transported at a much smaller courant number than your momentum you can dodge some of the stability issues involved with the sharp density gradient at the interface. Of course this isnt a fix all: you cant have too large a courant number difference between the two (at a given location) or you'll run into a different set of problems. Hope this helps