# multiphase flows

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

 January 13, 2003, 16:08 multiphase flows #1 Allen Badeau Guest   Posts: n/a What is the easiest way to treat d(rho)/dt in the Poisson equation for multiphase, incompressible flows, as d(rho)/dt is not equal to zero for mixtures.

 January 14, 2003, 11:32 Re: multiphase flows #2 su Guest   Posts: n/a can you give more details please!

 January 14, 2003, 11:39 Re: multiphase flows #3 Allen Badeau Jr. Guest   Posts: n/a Su, I would like to model two phase flows, specifically, liquid-liquid or liquid-gas. There will of course be two different densities, and assuming that we are only looking at the liquid-liquid mixing, with two different densities, the Poisson equation is solved for in the calculation of the pressure. In doing that, there is actually a d(rho)/dt in the equation, which is zero if there is only 1 phase and incompressibility. However, for multi-phase flows, I can't see how this is zero since it is actually the mixture density now, which does change with time. Saying that, I am not sure how to model d(rho)/dt in that equation. Once the density field is solved for, then the volume fraction of will be calculated.

 January 14, 2003, 18:07 Re: multiphase flows #4 Jim Park Guest   Posts: n/a The pressure equation is derived from the continuity equation. The details depend on the method. For instance, the formulation used in the Simple and Simple-like techniques is slightly different from the earlier work (MAC, SMAC) out of Los Alamos. In your case, you'll have two continuity equations, one for each component (or you can write one in terms of a total density and the other in terms of a mass fraction). Your best bet is to start with the momentum and continuity expressions appropriate to your problem and work through the details following the outline of one of the established methods. The d(rho)/dt term will appear as a part of that development. Good luck.

 January 14, 2003, 23:26 Re: multiphase flows #5 ptyue Guest   Posts: n/a If you mean 'mass derivative' by d/dt, d(rho)/dt=0 when fluid is incompressible, even for multiphase flow without phase change. If you want partial(rho)/partial(t), then you can calculate it from volume fraction( VOF method), level set function, etc. Good luck.

 January 15, 2003, 11:18 Re: multi-phase flows #6 su Guest   Posts: n/a Let us think in this way: first we need to write your equations then think of modeling liq-liq or gas-liq. etc. this is a simple example or just a keywork or an idea in general but let see; as i know you'll have a non-conservative form of the equations which is a MaJoR problem in two-phase flow area. once you accept that; you've 2 different densities; so you've 2 different conservation of mass, momentum and energy etc. so the frist equations may be written as (i'm writing \partial not d...) \partial_{t}(\alpha_{j} \rho_{j} + \partial_{x} (\alpha_{j} \rho_{j} u_{j}) = Sorce term (or 0). and then start solving it etc. with j = 1, 2 for phase 1 and phase 2. using mixture relations e.g. \rho = \alpha_{2} \rho_{2} + \alpha_{1} \rho_{1}, where we can write it as \partial_{t}(\rho) + \partial_{x} (\rho u) = Sorce term (or 0). so you can have a system of mixture-type two-phase flow;

 February 6, 2003, 09:38 Re: multiphase flows #7 niklas Guest   Posts: n/a As the poisson eq for pressure comes from cont.eqn I gues s one could write it as (ddt means partial ddt) [1] ddt(rho)+div(rho u)=0 [2] ddt(rho)+u div(rho)+rho div(u)=0 but if assuming that rho is proportional to eg. a convected vof or level set scalar Q ( ddt(Q)+u div(Q)=0 ) the first two terms in [2] vanishes. So no ddt(rho) left for the cont.eqn it seems.

 Thread Tools 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 On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Luk_Fiz CFX 11 April 4, 2013 05:29 lindstroem OpenFOAM Programming & Development 9 September 20, 2011 05:20 Vidar Main CFD Forum 0 September 18, 2008 05:05 Sam CFD-Wiki 1 April 26, 2007 14:50 discoganya Main CFD Forum 0 June 29, 2006 13:39

All times are GMT -4. The time now is 01:23.