# SIMPLE algorithm for two-phase

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 April 16, 2012, 13:11 SIMPLE algorithm for two-phase #1 Member   Join Date: Dec 2011 Location: State College, PA Posts: 87 Rep Power: 14 I'm writing a computer code to model a two-phase flow. I have read Patankar's book, "Numerical Heat Transfer and Fluid Flow", which discusses the SIMPLE algorithm. It was a great read, but he doesn't talk about two-phase flows at all. So far, I have setup my momentum equations in order to obtain my guessed velocity field, but I'm having trouble on what to do to setup the pressure correction equation. For a one-phase flow, you use the momentum equation to put the velocity correction, u', in terms of the pressure correction, p'. Then you put this into the continuity equation along with your guessed velocity field results in order to form a pressure matrix for your mesh that can be solved to get the corrected pressures. But for my two-phase flow, I have 2 continuity equations instead of one, which are coupled together by the void fraction. I assume that pressure is the same in the two phases. So now how do I arrive at the pressure correction? Do I need to solve both continuity equations simultaneously, using the old-iteration void fractions for my convection terms, in order to simultaneously calculate the pressure corrections and void fractions? And what about the energy equations? They also have the void fraction terms present, so do I also need to solve them simultaneously with the continuity equations as well? I would appreciate any advice and also any information on books dealing with applying SIMPLE to two-phase flow because I've had trouble finding anything of the sort.

 April 17, 2012, 09:22 #2 Member   Join Date: Dec 2011 Location: State College, PA Posts: 87 Rep Power: 14 To form the pressure correction equation, I considered total volume continuity, so that I only have one continuity equation. In other words, I consider all mass convection into the cell, be it vapor or liquid, and the change in total mass over time (meaning change in vapor AND liquid mass). If you consider total volume, there is no source term due to evaporation or condensation. This allowed me to form the pressure correction equation. I use density and void fraction explicit from the last iteration to form the pressure correction. This allows me to correct the pressure field and the velocity field, but now I need to solve for the remaining dependent variables, which are the void fraction and the phase enthalpies. I'm stuck here because I can see that I have 2 energy equations and 2 continuity equations. However, I only have 3 dependent variables --- void, liquid enthalpy, and vapor enthalpy. So it seems like an over-defined problem to me. I picked up a paper from the 2000 Proceedings of the ASME Fluids Engineering Division, Summer Meeting, entitled "A Phase Coupled Method for Solving Multiphase Problems on Unstructured Meshes", which explains an algorithm, PC-SIMPLE (phase-coupled SIMPLE). It confirms my use of total volume continuity to form the pressure correction equation, but I'm still at a loss on how I should now get my void fraction and enthalpies out of the remaining conservation equations.

 April 17, 2012, 11:53 #3 Member   Join Date: Dec 2011 Location: State College, PA Posts: 87 Rep Power: 14 My searching has lead me to a paper by Oliveira and Issa, "Numerical aspects of an algorithm for the Eulerian Simulation of Two-Phase Flows". It discusses two methods for handling the "over-definition" that I discussed. The first method is to solve just one of the continuity equations using the updated (corrected) velocities and then simply use (1-\alpha) to get the void fraction of the other phase that the continuity equation was not solved for. The second method is to solve both equations, but for two new, introduced variables \alpha_l* and \alpha_v*, which are preliminary void fractions for liquid and vapor phases. The actual liquid and vapor voids are then obtained using a correction factor, f, as follows: \alpha_l = f x \alpha_l* \alpha_v = f x \alpha_v* The f correction factor is defined as f = 1/(\alpha_l* + \alpha_v*). This allows the physics of both phases to be captured, but also ensures that the void fractions add up to 1 at each iteration. When convergence of the solution is achieved, f should be equal to 1.

 April 20, 2012, 05:45 #4 Senior Member   Join Date: Aug 2011 Posts: 272 Rep Power: 15 Yes what you mention is a very tricky problem !!! and papers on this topic are very rare. Especially when you deal with 3 different phases.. Check the serie of papers from Darwish and Moukalled. They have proposed several approaches to deal with such problem Good luck ;-)

 April 20, 2012, 09:02 #5 Member   Join Date: Dec 2011 Location: State College, PA Posts: 87 Rep Power: 14 I looked up some of their papers and they look like they're right along the lines of what I was looking for. Thanks for the suggestion.

 Tags simple algorithm, two-phase flow

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post johnhelt Main CFD Forum 2 October 18, 2010 06:27 lost.identity Main CFD Forum 1 October 7, 2010 11:48 peck110 Main CFD Forum 1 January 30, 2008 13:45 youngan CFX 0 July 1, 2003 23:32 Jonathan Castro Main CFD Forum 3 December 10, 1999 04:59

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

 Contact Us - CFD Online - Privacy Statement - Top