关于UDF的问题

ryyt0525

我做壁面化学反应的，求助各位大侠，下面的情况怎么编译UDF？下面的程序对吗？在组分里面能包括site物质和固体组分嘛？
#include "udf.h"
/*ARRHENIUS CONSTANTS*/
#define PRE_EXP 79725.6
#define ACTIVE 1.25e+07
#define BETA 0.3
real arrhenius_rate(real temp)
{
return
PRE_EXP*pow(temp,BETA)*exp(-ACTIVE/)UNIVERSAL GAS CONSTANT*temp));
}
/*物种数目。在FLUENT 面板中顺序必须匹配*/
#define SO2 0
#define O2 1
#define CuO 2
#define CuSO4 3
#define NUM_SPECS 4
/*Reaction Rate Routine 被用在两个UDFs 中*/
real reaction rate(cell t c,Thread *cthread,real mw[],real yi[])
{
real concenSO2 = C R(c, cthread)*yi[SO2]/mw[SO2];
return arrhenius_rate(C T(c,cthread))*pow(concenSO2,SO2 EXP);
}
DEFINE SR RATE(arrhenius,f,fthread,r,mw,yi,rr)
{
*rr = reaction rate(F CO(f,fthread),THREAD TO(f,fthread),mw,yi);
}
real contact_area(cell t c,Thread *t,int s id, int *n)
{
int I = 0;
real area =0.0, A[ND ND];
*n = 0;
c_face_loop(c,t,I)
{
if(THREAD ID(C FACE THREAD(c,t,I)) == s id)
{
(*n)++;
F AREA(A, C FACE(c,t,I), C FACE THREAD(c,t,I));
area += NV MAG(A);
}
}
return area;
}