porous media in FLUENT

Natalie2210 · June 3, 2013, 03:47

Hi everybody,

I'm simulating flow through porous media and it works well with FLUENT. However, I have a question concerning the theory behind the flow:

Basically, two different formulations can be used, one is denominated the "superficial velocity formulation" and the other one is the "physical velocity formulation". In the superficial velocity formulation, a sink term is added the usual momentum equation,

S = mu/alpha * v (in my model I neglect the inertia resistance)

mu is the dynamic viscosity, v the velocity and alpha denotes the permeability. However, the solution changes (as it should) if I use different values for the porosity gamma. I now wonder where the porosity enters the model equations within the superficial velocity formulation? As far as I know, the permeability does not depend on the porosity, however, in the ANSYS FLUENT 12.0 Users-Guide - 7.2.3, I find the following statement:
" Both K and C2 are functions of (1-gamma). "
What kind of relation is that and how is it theoretically justified?

Thank you for any piece of advice,
Natalie

oj.bulmer · June 3, 2013, 10:09

Since permeability refers to viscous losses, the nature of the porous medium will influence the value of alpha. Now naturally ,the "nature" here is inherently related with porosity $\gamma$ .

I will answer your question by taking a reference from Idelchik's handbook of hydraulic resistance, which is a widely accepted standard for hydraulic loss coefficients.

Darcy's equation for a pressure drop in porous medium is given as:

$\Delta P = \frac{\mu}{a} V L +\frac{1}{2} K \rho V^2$

V here is a free stream velocity before the porous medium - or superficial velocity,

is coefficient of resistance and L is the length of porous medium.

Idelchik generalises this pressure drop with linear and quadratic terms of V:

$\Delta P = k_1 V +k_2 V^2$

where $k_1, \;k_2$ are empirical constants which can be understood easily from earlier equation.

Now, intuitively, he defines pressure drop as: $\Delta P = \frac{1}{2} \zeta \rho V^2$

where, $\zeta=f(k_1,k_2)$

Thus Idelchik combines both viscous and inertial losses into $\zeta$ .

You will see that for different types of hydraulic resistances, he formulates the value of $\zeta$ (unitless) as a function of percentage open area or porosity $\gamma$ . He adds more terms for smaller Reynolds numbers, to account for viscous losses, thus implying that value of permeability should be influenced by porosity.

Incidently, I do observe that most of the analytical formulae he gives for $\zeta$ are functions of $(1-\gamma)$

Obviously, in FLUENT, the values of C2 are per unit length, arrived by dividing $\zeta$ by the length of porous zone, and then incorporated in the momentum sink, as a pressure gradient and not as a pressure drop.

$S_i = - \left( \frac{\mu}{a} V_i +\frac{1}{2} C_2 \rho |V| V_i \right)$

Hope this is useful.

OJ

Natalie2210 · June 4, 2013, 02:33

So FLUENT does assume that the permeability depends on the porosity and modifies the value of 1/alpha accordingly?

Thanks for your answer,
Natalie

oj.bulmer · June 4, 2013, 03:53

No, you have to set the value of permeability yourself as a part of your physics setup. FLUENT is just a solver, that solves your setup.

The gist of my post above was that both the permeability and inertial coefficient of resistance depend on porosity.

OJ

Natalie2210 · June 4, 2013, 04:44

I understand, thank you!

Natalie

ftab · May 4, 2014, 09:26

Dear Natalie and OJ,
I would appreciate if either of you would answer this issue.

The original question of Natalie which led me also to this thread, is still not answered. Also defining Zeta and other coefficients would involve porosity, but ignoring the loss coefficient and making C2=0, will leave our momentum sinkterm as (miu/permeability)*Velocity. Assume that permeability is a constant number here.

Why CFX is still asking for porosity? Why changing the porosity will change the simulation results as Natalie said, although it is not in our equations?

mrh111 · July 19, 2015, 20:33

Hi every body
I don't know how determine initial velocity for the flow in porous media.
Thank you for time and consideration.
Best regards

armanpournasiri · August 15, 2019, 12:46

hello
first: i want to say thank you for your guids about porous resistance .but i dont undrestant that how can i calcculet viscous and inertial resistans.
secund :I am working on porpous media problem (heat transfer) but i cant find a good tutorial about heat transfer into porous media
also i want to use source term in my problem.
Thanks for your answer

June 3, 2013, 03:47	porous media in FLUENT	#1
Natalie2210 Member Join Date: Feb 2013 Posts: 30 Rep Power: 13	Hi everybody, I'm simulating flow through porous media and it works well with FLUENT. However, I have a question concerning the theory behind the flow: Basically, two different formulations can be used, one is denominated the "superficial velocity formulation" and the other one is the "physical velocity formulation". In the superficial velocity formulation, a sink term is added the usual momentum equation, S = mu/alpha * v (in my model I neglect the inertia resistance) mu is the dynamic viscosity, v the velocity and alpha denotes the permeability. However, the solution changes (as it should) if I use different values for the porosity gamma. I now wonder where the porosity enters the model equations within the superficial velocity formulation? As far as I know, the permeability does not depend on the porosity, however, in the ANSYS FLUENT 12.0 Users-Guide - 7.2.3, I find the following statement: " Both K and C2 are functions of (1-gamma). " What kind of relation is that and how is it theoretically justified? Thank you for any piece of advice, Natalie Takue likes this.

June 3, 2013, 10:09		#2
oj.bulmer Senior Member OJ Join Date: Apr 2012 Location: United Kindom Posts: 473 Rep Power: 20	Since permeability refers to viscous losses, the nature of the porous medium will influence the value of alpha. Now naturally ,the "nature" here is inherently related with porosity $\gamma$ . I will answer your question by taking a reference from Idelchik's handbook of hydraulic resistance, which is a widely accepted standard for hydraulic loss coefficients. Darcy's equation for a pressure drop in porous medium is given as: $\Delta P = \frac{\mu}{a} V L +\frac{1}{2} K \rho V^2$ V here is a free stream velocity before the porous medium - or superficial velocity, is coefficient of resistance and L is the length of porous medium. Idelchik generalises this pressure drop with linear and quadratic terms of V: $\Delta P = k_1 V +k_2 V^2$ where $k_1, \;k_2$ are empirical constants which can be understood easily from earlier equation. Now, intuitively, he defines pressure drop as: $\Delta P = \frac{1}{2} \zeta \rho V^2$ where, $\zeta=f(k_1,k_2)$ Thus Idelchik combines both viscous and inertial losses into $\zeta$ . You will see that for different types of hydraulic resistances, he formulates the value of $\zeta$ (unitless) as a function of percentage open area or porosity $\gamma$ . He adds more terms for smaller Reynolds numbers, to account for viscous losses, thus implying that value of permeability should be influenced by porosity. Incidently, I do observe that most of the analytical formulae he gives for $\zeta$ are functions of $(1-\gamma)$ Obviously, in FLUENT, the values of C2 are per unit length, arrived by dividing $\zeta$ by the length of porous zone, and then incorporated in the momentum sink, as a pressure gradient and not as a pressure drop. $S_i = - \left( \frac{\mu}{a} V_i +\frac{1}{2} C_2 \rho \|V\| V_i \right)$ Hope this is useful. OJ Jonathan, ghost82 and alik0 like this.

May 4, 2014, 09:26		#6
ftab Member Ftab Join Date: Sep 2011 Posts: 87 Rep Power: 14	Dear Natalie and OJ, I would appreciate if either of you would answer this issue. The original question of Natalie which led me also to this thread, is still not answered. Also defining Zeta and other coefficients would involve porosity, but ignoring the loss coefficient and making C2=0, will leave our momentum sinkterm as (miu/permeability)Velocity. Assume that permeability is a constant number here. Why CFX is still asking for porosity? Why changing the porosity will change the simulation results as Natalie said, although it is not in our equations? Last edited by ftab; May 5, 2014 at 06:36.*

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Porous Media Simulation with CFX/Fluent	Ghasem	ANSYS	1	May 14, 2013 09:07
relative velocity resistance formulation, porous media	arminheidari	FLUENT	0	January 7, 2013 04:02
compressible flow through porous media	Chirag2302	FLUENT	0	March 2, 2012 23:13
Discrete phase model in porous media	Pradeep	FLUENT	0	February 3, 2006 16:30
porous media divergence	Phil	FLUENT	0	March 13, 2002 06:03

June 4, 2013, 02:33		#3
Natalie2210 Member Join Date: Feb 2013 Posts: 30 Rep Power: 13	So FLUENT does assume that the permeability depends on the porosity and modifies the value of 1/alpha accordingly? Thanks for your answer, Natalie

June 4, 2013, 03:53		#4
oj.bulmer Senior Member OJ Join Date: Apr 2012 Location: United Kindom Posts: 473 Rep Power: 20	No, you have to set the value of permeability yourself as a part of your physics setup. FLUENT is just a solver, that solves your setup. The gist of my post above was that both the permeability and inertial coefficient of resistance depend on porosity. OJ

June 4, 2013, 04:44		#5
Natalie2210 Member Join Date: Feb 2013 Posts: 30 Rep Power: 13	I understand, thank you! Natalie

July 19, 2015, 20:33		#7
mrh111 New Member mohamadreza Join Date: Jul 2015 Posts: 16 Rep Power: 10	Hi every body I don't know how determine initial velocity for the flow in porous media. Thank you for time and consideration. Best regards

August 15, 2019, 12:46		#8
armanpournasiri New Member arman Join Date: Feb 2016 Posts: 11 Rep Power: 10	hello first: i want to say thank you for your guids about porous resistance .but i dont undrestant that how can i calcculet viscous and inertial resistans. secund :I am working on porpous media problem (heat transfer) but i cant find a good tutorial about heat transfer into porous media also i want to use source term in my problem. Thanks for your answer