#include "udf.h"
#define E 2.71828
.……
DEFINE_SOURCE(my_source,c,t,dS,i)
{
real f=pow(0.1*C_T(c,t)-1.0,1.38);
dS[eqn]=
return 0.621*pow(0.1*C_T(c,t)-1.0)1.76)*pow(E,f);
}
其中dS[eqn]是S对该方程中所求解量的导数,我不知道是否为T,就没写上去。dA/dt手算后加上去的
参考udf manual的例子
参考帮助,很好写的