CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

NS Equations for multiphase flow in porous media

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 13, 2017, 11:00
Default NS Equations for multiphase flow in porous media
  #1
New Member
 
Quirin Grossmann
Join Date: Nov 2016
Posts: 1
Rep Power: 0
quir is on a distinguished road
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
(1)      \alpha_1 + \alpha_2 = 1
to
(2)      \alpha_1 + \alpha_2 = \epsilon
where \epsilon is the porosity.
The alpha-equation solved in fluid.solve() is the one suggested by Weller in Rusche's thesis:
(3)      \frac{\partial \alpha_1}{\partial t} + \nabla \cdot (\alpha_1(\alpha_1 \bar{U}_1 + \alpha_2 \bar{U}_2) + (\bar{U}_1 - \bar{U}_2) \alpha_1(1-\alpha_1))=0
However, this is only valid with equation (1). So, using equation (2) instead, an additional term appears:
(3)      \frac{\partial \alpha_1}{\partial t} + \nabla\cdot (\alpha_1(\alpha_1\bar{U}_1 + \alpha_2\bar{U}_2) + (\bar{U}_1 -  \bar{U}_2)\alpha_1(1-\alpha_1) + \alpha_1 \alpha_s U_2)=0
where \alpha_s=1-\epsilon. 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 \alphas are not bounded between 0 and \epsilon, 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.
quir is offline   Reply With Quote

Old   June 15, 2017, 08:49
Default
  #2
New Member
 
Bibin K.S
Join Date: Oct 2014
Posts: 19
Rep Power: 8
bibin.sme is on a distinguished road
Did u get any solution to this problem... can you share your solver... ???
bibin.sme is offline   Reply With Quote

Old   May 25, 2020, 10:48
Default
  #3
New Member
 
roshan
Join Date: Dec 2019
Posts: 10
Rep Power: 3
roshanak.rabiee878 is on a distinguished road
Hi
have you solved your problem?
I have questions about how you add effect of relative permeability into your solver


Roshanak

Quote:
Originally Posted by quir View Post
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
(1)      \alpha_1 + \alpha_2 = 1
to
(2)      \alpha_1 + \alpha_2 = \epsilon
where \epsilon is the porosity.
The alpha-equation solved in fluid.solve() is the one suggested by Weller in Rusche's thesis:
(3)      \frac{\partial \alpha_1}{\partial t} + \nabla \cdot (\alpha_1(\alpha_1 \bar{U}_1 + \alpha_2 \bar{U}_2) + (\bar{U}_1 - \bar{U}_2) \alpha_1(1-\alpha_1))=0
However, this is only valid with equation (1). So, using equation (2) instead, an additional term appears:
(3)      \frac{\partial \alpha_1}{\partial t} + \nabla\cdot (\alpha_1(\alpha_1\bar{U}_1 + \alpha_2\bar{U}_2) + (\bar{U}_1 -  \bar{U}_2)\alpha_1(1-\alpha_1) + \alpha_1 \alpha_s U_2)=0
where \alpha_s=1-\epsilon. 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 \alphas are not bounded between 0 and \epsilon, 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
roshanak.rabiee878 is offline   Reply With Quote

Old   December 26, 2020, 11:25
Default
  #4
Member
 
...
Join Date: May 2018
Posts: 37
Rep Power: 5
regard is on a distinguished road
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
regard is offline   Reply With Quote

Reply

Tags
alphaeqn, multiphase flow, porous media flow, twophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Multiphase flow in porous media and relative permeabilities sirpolar Fluent Multiphase 1 July 24, 2016 05:41
Flow through two combined porous media with diffrent permeability Sandee Main CFD Forum 0 March 28, 2015 10:35
Flow through porous media: permeability issue butterfly1 Main CFD Forum 1 January 2, 2014 02:34
Flow through porous media: permeability issue butterfly1 CFX 3 December 23, 2013 21:23
Two phase flow in porous media Madhavi Krishnan FLUENT 3 June 6, 2005 05:52


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