Granular pressure

Qingang Xiong · May 7, 2012, 11:37

Dear members,
I am incorporating KTGF for sovling the alphaEqn.H in twoPhaseEulerFoam. In this sovler, is the "ppMagf" to be differentiated alpha to particle granular temperature? If so, I think it is so complex. Looking forward to your reply. Thank you!

alberto · May 8, 2012, 02:24

Quote:

Originally Posted by Qingang Xiong

Dear members,
I am incorporating KTGF for sovling the alphaEqn.H in twoPhaseEulerFoam. In this sovler, is the "ppMagf" to be differentiated alpha to particle granular temperature? If so, I think it is so complex. Looking forward to your reply. Thank you!

The ppMagf term is related to the normal stress modulus of the phase $G(\alpha_s)$ . You should check the definition of this term in the relevant literature, to have some background on the model. If you consider a granular phase, with equation of state:

$p_s = p_s(\alpha_s, \Theta_s)$ ,

the normal stress modulus is

$G(\alpha_s, \Theta_s) = \frac{\partial p_s}{\partial \alpha_s}$

In twoPhaseEulerFoam ppMagf = $G(\alpha_s)/(\alpha_s \rho_s)$ , due to the treatment used in discretizing the phase momentum equation.

Best,

Qingang Xiong · May 8, 2012, 11:49

Quote:

Originally Posted by alberto

The ppMagf term is related to the normal stress modulus of the phase G(alpha). You should check the definition of this term in the relevant literature, to have some background on the model. If you consider a granular phase, with equation of state:

p_s = p_s(alpha_s, Theta_s),

the normal stress modulus is (d indicates partial derivative)

G(alpha_s, Theta_s) = d(p_s)/d(alpha_s)

In twoPhaseEulerFoam ppMagf = G(alpha_s)/(alpha_s*rho_s), due to the treatment used in discretizing the phase momentum equation.

Best,

Dear Alberto,
Thank you very much for your reply. So from your description, Theta_s can be regarded as an independent variable to alpha_s though Theta_s is implicitly obtained from alpha_s. Am I right?

alberto · May 8, 2012, 12:29

Quote:

Originally Posted by Qingang Xiong

Dear Alberto,
Thank you very much for your reply. So from your description, $\Theta_s$ can be regarded as an independent variable to $\alpha_s$ though $\Theta_s$ is implicitly obtained from $\alpha_s$ . Am I right?

You are, in short, treating the particle phase as a compressible fluid. All you need to compute it $1/\psi$ , being $\psi$ the compressibility of the particle phase, which is defined as

$\psi = \frac{\partial \alpha_s}{\partial p_s}$

This is analogous to the case of an ideal gas, where you have

$\rho = \frac{p}{RT}$

so

$\psi = \frac{\partial \rho}{\partial p} = \frac{1}{RT}$

The difference is that, in a granular flow, generally $\rho_s = const$ , and the compressibility is due to changes in the phase fraction $\alpha_s$ .

Best,

Qingang Xiong · May 8, 2012, 12:49

Quote:

Originally Posted by alberto

You are, in short, treating the particle phase as a compressible fluid. All you need to compute it $1/\psi$ , being $\psi$ the compressibility of the particle phase, which is defined as

$\psi = \frac{\partial \alpha_s}{\partial p_s}$

This is analogous to the case of an ideal gas, where you have

$\rho = \frac{p}{RT}$

so

$\psi = \frac{\partial \rho}{\partial p} = \frac{1}{RT}$

The difference is that, in a granular flow, generally $\rho_s = const$ , and the compressibility is due to changes in the phase fraction $\alpha_s$ .

Best,

Dear Alberto,
I understand what you mean. This point is clear to me. Thank you veru much!

Best

May 7, 2012, 11:37	Granular pressure	#1
Qingang Xiong New Member Qingang Xiong Join Date: May 2012 Location: Ames, IA, USA Posts: 20 Rep Power: 14	Dear members, I am incorporating KTGF for sovling the alphaEqn.H in twoPhaseEulerFoam. In this sovler, is the "ppMagf" to be differentiated alpha to particle granular temperature? If so, I think it is so complex. Looking forward to your reply. Thank you!

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Pressure instabilities with interDyMFoam for the floatingObject case	nbadano	OpenFOAM Running, Solving & CFD	15	October 15, 2021 06:35
Calculation of the Governing Equations	Mihail	CFX	7	September 7, 2014 06:27
Using pressure outlet BC	Shreyansh	FLUENT	0	November 5, 2010 14:57
Help ASAP! pressure inlet & outlet	engahmed	FLUENT	0	June 13, 2010 15:33
Setting pressure and velocity in inlet	Asghari	FLUENT	5	September 22, 2006 13:23