请教高手:DEGASSING边界条件设置时,加载SOURCE出现错误

wust10102074

请问各位高手,我用UDF加载DEGASSING边界时,在第二相中加入MASS 时FLUENT报错,仅加反作用力项时FLUENT很正常,这该怎样处理阿,程序如下:
#include "udf.h"

DEFINE_SOURCE(degassing_source, cell, thread, dS, eqn)
{
  real source;
  Thread *tm = THREAD_SUPER_THREAD(thread);

  source =  -C_R(cell,thread)*C_VOF(cell,thread)/CURRENT_TIMESTEP ;

  C_UDMI(cell,tm,0) = source;

  dS[eqn] =  -C_R(cell,thread)/CURRENT_TIMESTEP;

  return source;
}

[ 本帖最后由 wust10102074 于 2011-12-6 10:52 编辑 ]

030301029

#include "udf.h"
#include "sg.h"
#include "sg_mphase.h"
#include "flow.h"
#include "mem.h"
#include "metric.h"