UDF Heat Transfer
 

I'm working with multiphase flow, particles and air, using the approach Euler-Euler. And I've had problems to specify each one of the phases>

#include "udf.h"

#define PR_NUMBER(cp,mu,k) ((cp)*(mu)/(k))
#define IP_HEAT_COEFF(vof,k,Nu,dp) ((1.-vof)*6.*(k)*(Nu)/(dp)(dp))
#define AR_NUMBER(rof,rop,dp,mu) (9.81*((rop)-(rof))*(dp)*(dp)*(dp)*(rof)/mu/mu)
#define dp 0.00313356

real heat_kmiec_1980(cell_t c, Thread *ti, Thread *tj)
{
real h;
real k = C_K_L(c,ti);
real NV_VEC(v), vel, Re, Pr, Nu, Ar, vof, H0, dc, teta, n1, n2, n3, n4;

NV_DD(v,=, C_U(c,tj),C_V(c,tj),C_W(c,tj),-,C_U(c,ti),C_V(c,ti),C_W(c,ti));
vel = NV_MAG(v);
Re = RE_NUMBER(C_R(c,ti),vel,dp,C_MU_L(c,ti));
Pr = PR_NUMBER(C_CP(c,ti),C_MU_L(c,ti),k);
Ar = AR_NUMBER(C_R(c,ti),C_R(c,tj),dp,C_MU_L(c,ti));
H0 = 0.175;
dc = 0.30;
teta = 60;
n1 = tan (teta/2);
n2 = H0/dp;
n3 = dc/dp;
n4 = 0.79;
Nu = 0.0451*pow(Re,0.644)*pow(Pr,1./3.)*pow(Ar,0.226)*pow(n1,-0.852)*pow(n2,-1.47)*pow(n3,0.947)*pow(n4,2.304);

h = IP_HEAT_COEFF(C_VOF(c,tj),k,Nu,dp);
return h;
}


DEFINE_EXCHANGE_PROPERTY(heat_udf_1,*c,*t,*i,*j)
{
Thread**ti*=*THREAD_SUB_THREAD(t,i);
Thread**tj*=*THREAD_SUB_THREAD(t,j);
real*val;

val*=*heat_kmiec_1980(c,ti,tj);
return*val;
}*

Can you help me?