请教：只给液相施加源项~~~~~~~~~~~~~~~~

wangxw0315

用VOF模型做气相和液相的双相流模拟时，液相为第二相，如果只给液相加源项，UDF怎样写？
我仿照别人的写了一个，蛋不能通过编译，麻烦大家帮我看看
DEFINE_SOURCE(r_momsource,c,subthread,dS,eqn)
{
    Thread *mixture_thread;
    Thread *subthread=THREAD_SUB_THREAD(mixture_thread,1);
    real xr[ND_ND];
    real I=210;
    real rmomsource;
    C_CENTROID(xr,c,subthread);
     jz1=I(3.14*pow(xr[0],2));
     rmomsource=-6.28*(1e-07)*xr[1]*pow(jz1,2);
    return rmomsource;
}

[ 本帖最后由 wangxw0315 于 2013-1-17 16:20 编辑 ]