two-phase flow - outlet BC for pressure

michielm · November 29, 2013, 12:02

Hi,
We are trying to simulate the flow of a jet of liquid falling inside a gaseous chamber and exiting at the bottom with interFoam. We can set up the simulation such that it runs, but as soon as the liquid reaches the bottom the simulation crashes.

We believe this is caused by the following issue: we have assigned a fixedValue 0 for the pressure at the outlet and zeroGradient for the velocity.

As long as there is no liquid at the outlet this is fine, but as soon as there is both liquid and gas at the outlet you will have a capillary pressure jump across the interface. This fact, combined with the boundary condition that forces the pressure at the boundary inside BOTH fluids to 0 will mess things up.

Does anyone have advice on the appropriate pressure boundary condition that will allow both liquid and gas to flow out?

LESlie · December 11, 2013, 10:09

Hi Michiel,

Are you using VOF solvers for you two phase problem? In that case the BC of alpha can also play a role. You can try inletOutlet.

I would suggest you to use totalPressure for the pressure and pressureInletOutletVelocity for the velocity.

michielm · December 12, 2013, 05:07

Hi Leslie,
Thanks for the reply. Yes, I'm using VOF (interFoam).

I currently use 'advective' as the outflow boundary condition for alpha, which seems to do exactly what I want, namely just advect alpha along with the main flow.

I have tried the totalPressure condition for the pressure and it works somewhat (there is no crash), but it doesn't really make sense that much because totalPressure sets:
$p=p_0+1/2 \rho |U|^2$ so that would suggest that $1/2 \Delta \rho |U|^2$ equals the capillary pressure jump across the interface which is typically not the case.

What is your take on that?

November 29, 2013, 12:02	two-phase flow - outlet BC for pressure	#1
michielm Member Michiel Join Date: Oct 2010 Location: Delft, Netherlands Posts: 97 Rep Power: 16	Hi, We are trying to simulate the flow of a jet of liquid falling inside a gaseous chamber and exiting at the bottom with interFoam. We can set up the simulation such that it runs, but as soon as the liquid reaches the bottom the simulation crashes. We believe this is caused by the following issue: we have assigned a fixedValue 0 for the pressure at the outlet and zeroGradient for the velocity. As long as there is no liquid at the outlet this is fine, but as soon as there is both liquid and gas at the outlet you will have a capillary pressure jump across the interface. This fact, combined with the boundary condition that forces the pressure at the boundary inside BOTH fluids to 0 will mess things up. Does anyone have advice on the appropriate pressure boundary condition that will allow both liquid and gas to flow out? Last edited by michielm; December 12, 2013 at 05:07. Reason: Added the used solver

December 11, 2013, 10:09		#2
LESlie New Member NaiXian Leslie Lu Join Date: Jun 2009 Location: France Posts: 26 Rep Power: 17	Hi Michiel, Are you using VOF solvers for you two phase problem? In that case the BC of alpha can also play a role. You can try inletOutlet. I would suggest you to use totalPressure for the pressure and pressureInletOutletVelocity for the velocity. __________________ Cheers, Leslie LU

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Review: Reversed flow	CRT	FLUENT	1	May 7, 2018 06:36
Want Impeller Driven Fluid Flow: What Inlet and Outlet BC to use for Centrifugal Pump	Zev Xavier	FLUENT	3	May 9, 2016 07:42
Compressible flow, no data at the outlet	mireis	FLUENT	6	September 3, 2015 03:10
Disturbed flow field at outlet boundary (Multiphase flow through pipe)	Michiel	CFX	17	April 21, 2010 11:14
outlet boudary condition for a flow in the pipe	Atit	CFX	2	November 9, 2004 18:43

December 12, 2013, 05:07		#3
michielm Member Michiel Join Date: Oct 2010 Location: Delft, Netherlands Posts: 97 Rep Power: 16	Hi Leslie, Thanks for the reply. Yes, I'm using VOF (interFoam). I currently use 'advective' as the outflow boundary condition for alpha, which seems to do exactly what I want, namely just advect alpha along with the main flow. I have tried the totalPressure condition for the pressure and it works somewhat (there is no crash), but it doesn't really make sense that much because totalPressure sets: $p=p_0+1/2 \rho \|U\|^2$ so that would suggest that $1/2 \Delta \rho \|U\|^2$ equals the capillary pressure jump across the interface which is typically not the case. What is your take on that?