區(qū)域如圖所示,白色部分是固體,下面是氣體。其實想模擬的是固體受熱后發(fā)生燃燒。思路是取與固體接觸的那一層氣體單元,判斷那一側(cè)氣體單元,如果單元內(nèi)的溫度、壓力以及與它接觸的固體單元內(nèi)的溫度達到一定值后就從這層氣體單元內(nèi)加質(zhì)量源項和能量源項等。源項的代碼(示意代碼)如下:
其中的C_T(cell,thread和C_P(cell,thread)都是判斷氣體層單元的溫度和壓力,而F_T(f,t)本來是想著無法獲得旁邊固體單元的溫度,所以用這層邊界面的溫度代替,但是計算后發(fā)現(xiàn)這個F_T(f,t)變量根本沒起到作用。雖然源項能計算,但是明明表面溫度還很低源項就加進去了。請問應該如果處理?或者能不能直接得到與氣體相鄰的固體單元內(nèi)的溫度?
DEFINE_SOURCE(c_m, cell, thread, dS, eqn)
{
real source;
if (C_T(cell,thread) >= 320.&&C_P(cell,thread) >= 120000.&&F_T(f,t)>=1000.)
{
source =1
dS[eqn]=0;
}
else
{
source =0;
dS[eqn] = 0;
}
return source;
}
DEFINE_SOURCE(c_e, cell, thread, dS, eqn)
{
real source;
if (C_T(cell,thread) >= 320.&&C_P(cell,thread) >= 120000.&&F_T(f,t)>=1000.)
{
source =1;
dS[eqn] = 0;
}
else
{
source = 0;
dS[eqn] = 0;
}
return source;
}
![近壁區(qū)域的源項(壁面溫度)]()
123321.png |