UDF能量方程，内附程序

xi851103

做动网格计算，UDF能量源项程序如下，也设置了memory，source应该是加到能量方程中去了，但是计算后处理，C_UDMI(c,t,4)中得数为0，这是为什么呀？请高手指点，看看我的source有问题没？



DEFINE_SOURCE(energy,c,t,dS,eqn)
{
  double source;
  double A=5.0e+04;                    //电流值: 是50kA
  double time60=15e-3;            //起弧时间
  double timecur=CURRENT_TIME+15e-3;         //当前时间：
  double dt=timecur-time60;
  int    n=RP_Get_Integer("time-step");         /*当前时间步*/
  double resis;
  double TIME=CURRENT_TIMESTEP;
  double arc14ms=-(timecur-time60)*2*pi*50;
  double VOLUME;                              /*单元体积*/
  
  //Thread *t;
  //cell_t c;
    dS[eqn]=0.0;

    temp=C_T(c,t);
    if(temp>=5000)
    {source=pow(A*sin(arc14ms),2)*C_UDMI(c,t,2)-C_UDMI(c,t,1);
         return source;}
    else if((temp>=4000)&&(temp<5000))
        {resis=10000/0.005;
         source=pow(A*sin(arc14ms),2)*resis-C_UDMI(c,t,1);
         return source;}
    else if((temp>3000)&&(temp<4000))
        {resis=10000/0.002;
         source=pow(A*sin(arc14ms),2)*resis-C_UDMI(c,t,1);
         return source;}
    else
    {source=0;
         return source;}
        
     C_UDMI(c,t,4)=source;
     
     return source;
}