一个两相流UDF问题

lihongyang0

各位大侠好！我是UDF新手，在处理一个二维瞬态两相加空化问题。需要我对两相界面进行实时捕捉，并且对界面的一部分施加热流（反映在温度梯度上）。下面是我的程序：
#include "udf.h"
DEFINE_ADJUST(my_adjust,d)
{
Domain *mixture_domain;
cell_t c;
int phase_domain_index=1;
int temp1;
int temp2;
int K;
face_t f;       
Thread *mixture_thread;
Thread *subthread=THREAD_SUB_THREAD(mixture_thread,phase_domain_index);
real xc[ND_ND];
thread_loop_c(c,subthread)
{
C_CENTROID(xc,c,subthread);
if(xc[0]<0.04586)
{
temp=C_T(c,subthread);
K=C_K_L(c,subthread);

C_T_G(c,subthread)[1]=300000000/K;
}
}
}
请路过的各位大侠帮看看，指点小弟一二，感激不尽！！！