[求助] fluent中添加氧气源项，编译能通过，一计算就出错，附程序，求大牛解答

UDM已经定义，因为程序较多，只贴了氧气的源项，因为其他的源项加入的时候没有错误。 DEFINE_SOURCE(source_o2_g,c,t,dS,eqn)                                   /*氧气的质量源项，定义在气相中*/ {         Thread *mixture_t=THREAD_SUPER_THREAD(t);                           /*混合相的指针*/         Thread *ts=THREAD_SUB_THREAD(mixture_t,1);                                                                                                   real Ro13,Ro14,Ro15,Ro16,Ro17,Ro18,Ro19;         int n;         face_t f;         Thread *tf;         Thread *t1;         cell_t c1;         Ro17=0.0;         if(C_UDSI(c,ts,4)>1170.0)                 Ro13=1170.0;         else if(C_UDSI(c,ts,4)<730.0)                 Ro13=730.0;         else                 Ro13=C_UDSI(c,ts,4);         Ro14=2500.0*exp(-6420.0/Ro13);                                                                                                              Ro15=(Ro14+2.0)/(2.0+2.0*Ro14);                                            /*氧气的化学反应当量系数*/         Ro16=-C_UDMI(c,mixture_t,8)*16.0-C_UDMI(c,mixture_t,9)*16.0-C_UDMI(c,mixture_t,10)*16.0-C_UDMI(c,mixture_t,7)/12.0*Ro15*32.0;                                                   c_face_loop(c,t,n)                                                                                                                    {                 f=C_FACE(c,t,n);                 tf=C_FACE_THREAD(c,t,n);                 if(F_FLUX(f,tf)>0.0)                 {                         Ro17+=F_FLUX(f,tf)*C_YI(c,t,5);                 }                 else                 {                         t1=F_C1_THREAD(f,tf);                         c1=F_C1(f,tf);                         Ro17+=-F_FLUX(f,tf)*C_YI(c,t,5);                 }         }         Ro18=Ro16+Ro17/C_VOLUME(c,t);                                                                                                         return Ro18;         dS[eqn]=0; } 添加源项计算后就出现这种问题 [ Last edited by longma135 on 2014-3-13 at 16:12 ]

回复此楼