在固體wall邊界附近的一層單元內(nèi)添加源項(xiàng)，本意是當(dāng)單元壓力大于某一值、單元上wall的那條面溫度大于某一值時(shí)就釋放源項(xiàng)，代碼如下：
可是運(yùn)行的時(shí)候總是整個(gè)wall的溫度大于設(shè)定值時(shí)才釋放，而不是單元內(nèi)的與wall對(duì)應(yīng)的那個(gè)面的溫度大于值就在該單元釋放，問(wèn)題在哪？
C_T到底是整個(gè)wall的溫度還是單元內(nèi)那個(gè)面的溫度？
#include "udf.h"
DEFINE_SOURCE(c_e, cell, thread, dS, eqn)
{
    real source;
  int n;
   double st;
   face_t f;
   Thread *t;
  Domain *domain=Get_Domain(1);
   t=Lookup_Thread(domain,15);  

begin_f_loop(f,t)
{
       st=F_T(f,t);
    if ( C_P(cell,thread) >= 120000. && st>=310.)
    {
      source =1.10921e11;
      dS[eqn] = 0;
     }
    else
    {
      source = 0;
      dS[eqn] = 0;
     }
  return source;
}   
end_f_loop(f,thread)
} 返回小木蟲查看更多