请问：如何在气粉两相流中分别对两项加源项

hahasnow

新手问：
我想在气粉两相流中分别对两项加源项。
这是我的处理(对其中一相加源项)，不知错在那里。我用的是Fluent 6.0 Version. 刚刚接触UFD，仍不知如何入手。编译已通过。但不运行。
====================
这是出错信息：
FLUENT received fatal signal (ACCESS_VIOLATION)
1. Note exact events leading to error.
2. Save case/data under new name.
3. Exit program and restart to continue.
4. Report error to your distributor.
Error Object: ()
========================
#include "udf.h"
#define C2 100.0
DEFINE_SOURCE(xmom_source, cell, cell_thread, dS, eqn)
{
  real x[ND_ND];
  real con, source;
  int phase_domain_index;
  Domain *subdomain;
  Domain *mixture_domain;
  sub_domain_loop(subdomain,mixture_domain,phase_domain_index)
{
     if (DOMAIN_ID(subdomain) == 2)
    {
      thread_loop_c (cell_thread,subdomain)
      {
        begin_c_loop_all (cell,cell_thread)
        {
         con=C2*0.5*C_R(cell,cell_thread);
         source=-con*fabs(C_U(cell,cell_thread))*C_U(cell,cell_thread);
         dS[eqn]=-2.*con*fabs(C_U(cell,cell_thread));
        
         return source;
        }
        end_c_loop_all (cell,cell_thread)
      }
    }
}
}
谢谢关注