[i3abdulla]

Hi,

I'm looking for help to find the source of error in my UDF (code below)

So when I interpret the c code into Fluent 18.2, I get an error message saying:

cpp -I"C:\PROGRA~1\ANSYSS~1\\\Desktop\\CFD Model\\Hetrogenous_UDF.c: line 13: parse error.
Error: C:\\Users\\a.alshehhi\\Desktop\\CFD Model\\Hetrogenous_UDF.c: line 20: subthread: undeclared variable

I tried to solve it but unfortunately I couldn't, I wish if someone can help me.

Regards,


- - - - - - Start of Code -------

#include "udf.h"

DEFINE_HET_RXN_RATE(rxratedecomp,c,t,r,mw,yi,rr,rr _t)

{
Domain *d = Get_Domain(1);

const real R=8.3145;
real T, P, A0, EA, mole_frac, rho, vol_frac;
real conc [2];
const int POROUS_ID=12 , LIQ_PHASE_ID=1;
const enum {gas_phase, liquid_phase}
const enum {h2o2_w, h2o_w}
const enum {air, o2_g, h2o_g}

// phase thread pointer to the liquid phase
Thread *subthread = THREAD_SUB_THREAD(t, LIQ_PHASE_ID);

//volume fraction of the liquid phase
vol_frac=C_VOF(c,subthread);

// if the liquid is in the porous zone
if (Lookup_Thread(d, POROUS_ID) == t )
{

//pre-exponential factor for catalytic decomposition
A0=2.6*pow(10,4);

// activation energy for catalytic decomposition
EA=52.5*pow(10,3);
}

else{

//pre-exponental factor for catalytic decomposition
A0=1.0*pow(10,5.8);

// activation energy for thermal decomposition
EA=71.0*pow(10,3);
}

// temperature of the liquid phase
T=C_T(c,subthread);

// dens i ty of the liquid phase
rho=C_R(c,subthread);

//mole fraction of HP
mole_frac=yi[liquid_phase][h2o2_w]/mw[liquid_phase][h2o2_w];

// concentration of HP given by ideal gas law (kmol/m3)
conc[h2o2_w]=mole_frac*rho;

// reaction rate kmol /(m3, s )
*rr = 0.02*A0*exp(-EA/(R*T))*conc[h2o2_w]*vol_frac;
}

[i3abdulla]

Is there anyone have an answer?

[i3abdulla]

I hope if someone can help me