CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Fluent UDF and Scheme Programming (http://www.cfd-online.com/Forums/fluent-udf/)
-   -   Heat Transfer coefficient udf (http://www.cfd-online.com/Forums/fluent-udf/108041-heat-transfer-coefficient-udf.html)

Musa October 12, 2012 22:52

Heat Transfer coefficient udf
 
I am attempting to write a udf to calculate the htc from the gunn correlation but I keep getting a parse error in line 6. I have tried reviewing the code several times but I cannot find the source of the error. Any suggestions, recommendations?
#include "udf.h"

#define PR_NUMBER(cp,mu,k) ((cp)*(mu)/(k))
#define IP_HEAT_COEFF(vof,k,nu,d) ((vof)*6.*(k)*(Nu)/(d)/(d))

static real heat_gunn(cell_t c, Thread *ti, Thread *tj)
{
real h;
real vof_g=C_VOF(c,ti);
real d = C_PHASE_DIAMETER(c,tj);
real k = C_K_L(c,ti);
real NV_VEC(v), vel, Re, Pr, Nu;

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,d,C_MU_L(c,ti));
Pr = PR_NUMBER (C_CP(c,ti),C_MU_L(c,ti),k);
Nu = (7.-10.*vof_g+5.*pow(vof_g,2.))*(1+0.7*pow(Re,0.7)*pow (Pr,1./3.))+(1.33-2.40*vof_g+1.20*pow(vof_g,2))*pow(Re,0.7)*pow(Pr,1 ./3.);

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

sadjad.s October 13, 2012 04:29

hi mate,
use "THREAD_SUB_THREADS()" macro to separate ti and tj.
and also are you sure lines 2 and 3 ('#define' ) are right? i don't know either.

junianatalia May 18, 2016 13:38

Kmiec correlation
 
Hello,

I'm try write a Kmiec correlation, but I don't get sucess.
#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(rhof,rhop,dp,mu) (9.81*((rhop)-(rhof))*(dp)*(dp)*(dp)*(rhof)/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, 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;
n1 = 0.5753;
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,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?


All times are GMT -4. The time now is 03:26.