# NS Equations for multiphase flow in porous media

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

 February 13, 2017, 11:00 NS Equations for multiphase flow in porous media #1 New Member   Quirin Grossmann Join Date: Nov 2016 Posts: 1 Rep Power: 0 Hi everyone, I've been stuck on following problem for a while: I am trying to model multiphase (gas-liquid) flow through a packed bed for which I am using twoPhaseEulerFoam in OpenFOAM 3.0.1. To my understanding, all I need to do is change following equation to where is the porosity. The alpha-equation solved in fluid.solve() is the one suggested by Weller in Rusche's thesis: However, this is only valid with equation (1). So, using equation (2) instead, an additional term appears: where . I have also added a momentum sink term to the momentum equations to account for liquid-solid and gas-solids interaction. My problem is that the s are not bounded between 0 and , but still go as high as 1. I have tried defining an alphaMax in phaseProperties which should give a bound for the MULES solver, however, this doesn't make a difference. Any ideas on what could be the problem? Thanks Last edited by quir; February 14, 2017 at 10:37.

 June 15, 2017, 08:49 #2 New Member   Bibin K.S Join Date: Oct 2014 Posts: 19 Rep Power: 8 Did u get any solution to this problem... can you share your solver... ???

May 25, 2020, 10:48
#3
New Member

roshan
Join Date: Dec 2019
Posts: 10
Rep Power: 3
Hi

Roshanak

Quote:
 Originally Posted by quir Hi everyone, I've been stuck on following problem for a while: I am trying to model multiphase (gas-liquid) flow through a packed bed for which I am using twoPhaseEulerFoam in OpenFOAM 3.0.1. To my understanding, all I need to do is change following equation to where is the porosity. The alpha-equation solved in fluid.solve() is the one suggested by Weller in Rusche's thesis: However, this is only valid with equation (1). So, using equation (2) instead, an additional term appears: where . I have also added a momentum sink term to the momentum equations to account for liquid-solid and gas-solids interaction. My problem is that the s are not bounded between 0 and , but still go as high as 1. I have tried defining an alphaMax in phaseProperties which should give a bound for the MULES solver, however, this doesn't make a difference. Any ideas on what could be the problem? Thanks

 December 26, 2020, 11:25 #4 Member   ... Join Date: May 2018 Posts: 37 Rep Power: 5 hi all can any body help me please, i want to change my alpha equation for two phase flow in porous media and i need to add effect of relative permeability and capillry pressure thanks in advance

 Tags alphaeqn, multiphase flow, porous media flow, twophaseeulerfoam