FLUENT UDS問題

and761220

我用FLUENT裡的UDS自定義了一個方程式(變數為phi)，其中phi 在我的source term裡也會用到，所以我在 DEFINE_SOURCE macro裡使用了 C_UDSI(c,t,0) 來定義我的源項，如下 :

=====================================================================
#include "udf.h"
#define R 8.314
#define E 21180
#define C 59.19
#define SAT_DEN 7259
#define A 17.61
#define B 3704.6

DEFINE_SOURCE(uds_source,c,t,dS,eqn)
{       
        real source, P_eq,p_g,tem ;       
        tem   = C_T(c,t) ;
        p_g   = C_P(c,t);
        P_eq= 1000*exp(A-B/tem);
        source = C*exp(-E/(R*tem))*log(p_g/P_eq)*(SAT_DEN-C_UDSI(c,t,0));
        dS[eqn] = 0 ;       
        return source ;
}

===========================================================

現在遇到的問題是，最後算出的結果，phi的分布會很均勻，也就是不隨位置變化，只跟時間有關，這跟我所認知的結果明顯不太一樣，我要的是phi在不同的位置上會有變化。想請問用是否為程式碼邏輯上的問題，所以沒辦法實現我的想法，謝謝!!