CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   FLUENT (https://www.cfd-online.com/Forums/fluent/)
-   -   porous turbulence model help (https://www.cfd-online.com/Forums/fluent/138978-porous-turbulence-model-help.html)

hitchenxue July 15, 2014 08:39

porous turbulence model help
 
I use ANSYS Fluent to do some work on the turbulent flow in porous media with macroscopic model. I use UDF to add the additional source caused by porous media with p-dL model. Now I would like to validate my work using the flowing papers.

Marcelo J.S. de Lemos.Turbulent kinetic energy in a moving porous bed” in International Communications in Heat and Mass Transfer 35 (2008) 1049–1052

Marcos H. J. Pedras, Marcelo J. S. de Lemos. On the Mathematical Description and Simulation of Turbulent Flow in a Porous Medium Formed by an Array of Elliptic Rods.Journal of Fluids Engineering 2001
Renato A. Silva, Marcelo J.S. de Lemos. Turbulentflow in a composite channel.International Communications in Heat and Mass Transfer 38 (2011) 1019–1023.
etc....

I did not get the same resluts.
I use physical velocity for all the simulations, and my UDF the additional terms with P-dL model are attached. please help me out, I have been doing the verification for two months, and found nothing to solve my problem.

/************************************************** *********************/
/* P-dL(,2001) */
DEFINE_SOURCE(kpa_PdL,c,t,dS,eqn)
{
real source,ck,um;
real fai,dp,k,con;

ck=0.28;

Properties(c, t, &fai, &dp, &k, &con);
um=sqrt( C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t) );
source =fai* ck*C_R(c,t)*C_K(c,t)*fai*um/sqrt(k);

dS[eqn]=fai* ck*C_R(c,t)*fai*um/sqrt(k);
return source;
}
DEFINE_SOURCE(ybsl_PdL,c,t,dS,eqn)
{
real source,ck,c2,um;
real fai,dp,k,con;

ck=0.28; c2=1.92;

Properties(c, t, &fai, &dp, &k, &con);
um=sqrt( C_U(c,t)*C_U(c,t)+C_V(c,t)*C_V(c,t) );
source =fai* ck*c2*C_R(c,t)*C_D(c,t)*fai*um/sqrt(k);

dS[eqn]=fai* ck*c2*C_R(c,t)*fai*um/sqrt(k);
return source;
}


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