PrantikDas
Subscriber
#include"udf.h"
DEFINE_VECTOR_EXCHANGE_PROPERTY(burnetal,cell,mixture,liquid,solid,Ftdvector)
{
Thread*thread_l,*thread_s;
realx_vel_l,y_vel_l,z_vel_l,x_vel_s,y_vel_s,z_vel_s,slip_x,slip_y,slip_z,slip,rho_l,Cd,Sigmapq,diam,alpha_l;
real alpha_s,term,rey;
realpx,py,pz,sx,sy,sz,tau_p,D,K,mu_l,Kpq;
realCtd=1;
thread_l=THREAD_SUB_THREAD(mixture,liquid);
thread_s=THREAD_SUB_THREAD(mixture,solid);
/*Definingtheslipvelocity*/
x_vel_l=C_U(cell,thread_l);
y_vel_l=C_V(cell,thread_l);
z_vel_l=C_W(cell,thread_l);
x_vel_s=C_U(cell,thread_s);
y_vel_s=C_V(cell,thread_s);
z_vel_s=C_W(cell,thread_s);
slip_x=x_vel_l-x_vel_s;
slip_y=y_vel_l-y_vel_s;
slip_z=z_vel_l-z_vel_s;
slip=sqrt(pow(slip_x,2)+pow(slip_y,2)+pow(slip_z,2));
/*Densityoftheliquid*/
rho_l=C_R(cell,thread_l);
/*DefiningsigmausedinBurnsetalformulainANSYSTHEORYGUIDE*/
Sigmapq=0.9;
/*Defininglaminarviscosity,diameterandvolumefraction*/
mu_l=C_MU_L(cell,thread_l);
diam=C_PHASE_DIAMETER(cell,thread_s);
alpha_l=C_VOF(cell,thread_l);
alpha_s=C_VOF(cell,thread_s);
/*DefiningReynoldsNumber*/
rey=(rho_l*slip*diam)/mu_l;
/*CalculationofCdbyusinggidaspowmodel*/
/*DeterminationofmomentumexchangecoefficientbygidaspowmodelforsolidliquidmomentumexchangeandtheformulahasbeentakenfromANSYSTHEORYGUIDE*/
if(alpha_l>.8)
{
if(rey==0)
{
Cd=0;
}
else
{
Cd=(24*pow((rey*alpha_l),-1)*(1+(0.15*pow((alpha_l*rey),0.687))));
}
Kpq=(.75*Cd*alpha_l*alpha_s*rho_l*slip*pow(alpha_l,-2.65))*pow((diam*4),-1);
}
else
{
Kpq=(150*(alpha_s*(1-alpha_l)*mu_l)*pow(alpha_l*pow(diam,2),-1))+(1.75*rho_l*alpha_s*slip*pow(diam,-1));
}
/*Definingthedispersionscalar*/
D=C_MU_T(cell,thread_l)*pow(rho_l,-1);
/*Terminfrontofgradient*/
term=Ctd*Kpq*D*pow(Sigmapq,-1);
/*CalculationofGradients*/
px=C_VOF_G(cell,thread_l)[0]/alpha_l;
py=C_VOF_G(cell,thread_l)[1]/alpha_l;
pz=C_VOF_G(cell,thread_l)[2]/alpha_l;
sx=C_VOF_G(cell,thread_s)[0]/alpha_s;
sy=C_VOF_G(cell,thread_s)[1]/alpha_s;
sz=C_VOF_G(cell,thread_s)[2]/alpha_s;
/*CalculationoftheFinalExpressionandtheformulahasbeentakenfromANSYSTHEORYGUIDE*/
Ftdvector[0]=term*(sx-px);
Ftdvector[1]=term*(sy-py);
Ftdvector[2]=term*(sz-pz);
}
I have commented on the code though I am getting the same error.