源项的问题

nsdzhl

为什么，我在能量源项中加入一个常数，就可以计算下去并达到收敛
而卧改为在能量源项中加入udf后，udf不起作用。
结论udf写错了，但是错在哪儿了？
求教：大家帮忙看看这个能量源项的udf为什么不能起作用？
DEFINE_SOURCE(Jouleheating_radiationlossed_source,c,t,dS,eqn)  /*焦耳热*/
{
real x[ND_ND];
real source3;
real Q;
real temper;
real sigma;
real U0;
real U1;
temper=C_T(c,t);
U0=C_UDSI(c,t,0);
U1=C_UDSI(c,t,1);
sigma=a1*temper*temper*temper*temper*temper*temper-a2*temper*temper*temper*temper*temper+a3*temper*temper*temper*temper-a4*temper*temper*temper+a5*temper*temper-0.6169*temper+222.37;
if(temper>9500)
Q=5.6e03*(temper-9500)+181*(temper-9500)*(temper-9500);
source3=1.77e14*sigma*(U0*U0+U1*U1)-Q;
dS[eqn]=0;
return source3;
}

ltzzq

你好 楼主能伴我指点一下吗 我该用哪些宏写这个程序呢 不胜感激啊http://www.cfluid.com/bbs/viewthread.php?tid=105261&extra=page%3D1

feihu0105

看附件中的4.3.8节有详细介绍。