求助UDF编译可变体积源程序

学习蜗牛

各位大哥拉我一把，纠结了很久了 ，物理问题见附件    下面是我编的程序，加载没错，但结果温度基本没变，郁闷啊，一天一夜了 。。。
#include "udf.h"
DEFINE_SOURCE(cell_q_source,cell,thread,dS,eqn)
{
   real x[ND_ND];
   real source;
   real temp=C_T(cell,thread);
   real ti = RP_Get_Real("flow-time");
C_CENTROID(x, cell,thread);
   if(ti<=1)
         {source=502400*exp(-59512*x[0]*x[0])-2*1260*(temp-298);
dS[eqn]=0;
}
else
   {source=-2*1260*(temp-298);
       dS[eqn]=0;
}
return source;
}

swimming02

udf本身我感觉没什么问题，看看计算过程中的设定是否正确，包括边界条件，计算时间步长等。

学习蜗牛

出来了 呵呵 原来书上的数据给错了 谢谢哈