# pressure in twoPhaseEulerFoam

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

 September 8, 2010, 12:04 pressure in twoPhaseEulerFoam #1 New Member   Alessio Alexiadis Join Date: Mar 2009 Location: Stockholm, Sweden Posts: 15 Rep Power: 10 Sponsored Links Hi everybody, I have a question concerning twoPhaseEulerFoam. I have a rectangular channel where gas is introduced from one of the walls. The left wall and the bottom are just walls, the right wall is the gas-inlet and the top is open. The gas escapes from the top, which is the free liquid surface. it is laminar. As boundary conditions at the top I tried to replicate what are something called 'degassing b.c.' in other codes. The top, for the gas is an outlet. The liquid, on the other hand, cannot escape the boundary and I impose a free velocity slip as it should be for a gas-liquid surface. The pressure at the top is fixed. The results seem logical. There is a transient when the gas gradually fills the channel and when it reaches a certain value it oscillates around this value. Now, however, I have a doubt. If the liquid cannot escape the channel, how can the overall void fraction change? I inject gas, but liquid cannot exit; how it is possible in an incompressible liquid? A real system would expand, but in my case the boundary is fixed. I tried to understand it from the code, but I am not sure I really understood it. I think the pressure is not maintained fixed but it changes in order to compensate for the lack of expansion. It is correct? If it is, do you know how exactly is the pressure changed? In the b.c., I imposed fixed pressure at the top, is this overridden by the code? And, if this is the case, what exactly happens: the pressure is maintained fixed for the gas phase, but changes for the liquid phase? How is it possible? Thank you in advance and best regards Alessio

 September 9, 2010, 10:23 #2 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,910 Rep Power: 28 Hi, one question before trying to answer: is the channel filled of liquid, with phase fraction equal to 1 at the beginning of the simulation, and then you start injecting gas? Best, __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

September 9, 2010, 11:00
#3
New Member

Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 15
Rep Power: 10
Quote:
 Originally Posted by alberto Hi, one question before trying to answer: is the channel filled of liquid, with phase fraction equal to 1 at the beginning of the simulation, and then you start injecting gas? Best,

Yes, it is

 September 9, 2010, 12:05 #4 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,910 Rep Power: 28 OK. So your outlet BC actually lets the liquid go out during the expansion (most probable), or no gas could enter the domain. One trick could be to patch the initial volume fraction with a fraction lower than 1 and close to the expected value, but it is cheating a bit ;-) Best, __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

September 9, 2010, 12:29
#5
New Member

Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 15
Rep Power: 10
Quote:
 Originally Posted by alberto OK. So your outlet BC actually lets the liquid go out during the expansion (most probable), or no gas could enter the domain. One trick could be to patch the initial volume fraction with a fraction lower than 1 and close to the expected value, but it is cheating a bit ;-) Best,

Hi Alberto,
>So your outlet BC actually lets the liquid go out during the expansion

how it is possible? The b.c. I specify do not allow it. Are they, somehow overriden?

 September 10, 2010, 02:14 #6 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,910 Rep Power: 28 You fix a pressure value at the outlet and a slip condition for the velocity on the same boundary, if I understood it correctly? Best, __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. Last edited by alberto; September 10, 2010 at 18:22.

September 10, 2010, 11:42
#7
New Member

Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 15
Rep Power: 10
Quote:
 Originally Posted by alberto You fix a pressure value at the outlet and a slip condition for the velocity on the same boundary, if I understood it correctly. The slip condition does not set the velocity to zero anywhere. It states that the velocity gradient normal to the boundary is zero, and lets the flow move freely in the direction tangential to the boundary. Best,
Hi Alberto,

thanks for your answer again. Yes, the tangential velocity is not zero, but the liquid does not exit the channel because the normal velocity is zero. Moreover, I tested a case with v_liquid=0 (only the gas can exit) and it is almost the same. Imagine just a single big computational cell. At time t=0, this is completely full of liquid (alpha=0). Then, some gas is injected and let's say at time t1 alpha=0.1. My top boundary conditions does not allow the liquid to exit the computational cell, but now beta=1-alpha=0.9. Thus, if the liquid cannot exit, where is the 10% of the incompressible liquid missing?
They told me that, in this case, mass is not conserved and the code allows somehow some liquid to "disappear", but I do not understand how, numerically, this is handled. Do you have any information about this?

Alessio

 September 10, 2010, 18:47 #8 Senior Member   Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,910 Rep Power: 28 Hi, sorry for my last reply, which contained an error, since clearly the normal component of U is set to zero with the slip condition, or the "wall" would be permeable. Check the phase fluxes, and their change to see where the mass goes. The code enforces the total mass conservation. Best, __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

 Tags degassing, pressure, twophaseeulerfoam

 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 Mihail CFX 7 September 7, 2014 06:27 cdevalve FLUENT 3 January 14, 2012 01:11 Jan Main CFD Forum 1 October 30, 2006 08:54 Asghari FLUENT 5 September 22, 2006 13:23 Chris CD-adapco 2 June 21, 2005 08:20