| 5 | 1/1 | 返回列表 |
| 查看: 1252 | 回復(fù): 2 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
一只小火車新蟲 (初入文壇)
|
[求助]
氣液兩相流非均相反應(yīng)的UDF
|
|
|
參考之前的一篇帖子,自己寫的udf,計(jì)算結(jié)果始終是沒有發(fā)生反應(yīng)。 反應(yīng)式是A(l)=B(l)+C(g),VOF模型,A和B是混合相(第一相,液相),B第二項(xiàng)氣相。希望高手指點(diǎn)。 附程序: #include "udf.h" #include "mem.h" #include "sg_mphase.h" #define FLUID_ZONE_ID 7 /* zone ID of reaction */ #define R 8.31434 /* gas constant J/K-mol */ #define Pre 3.97e+12 /* pre-exponential factors units 1/s */ #define E1 228000 /* Activation energies (J/mol) */ DEFINE_HET_RXN_RATE(my_rate, c, t, r, mw, yi, rr, rr_t) { int zone_ID = 0; Thread **pt = THREAD_SUB_THREADS(t); Thread *prim_t = pt[0]; /* thread for primary phase */ Thread *sec_t = pt[1]; /* thread for secondary phase */ real T_prim = C_T(c,prim_t); /* primary phase temperature (K) */ real C2=C_YI(c, prim_t, 0); /* mass fraction of P in primary phase*/ zone_ID = THREAD_ID(t); if (zone_ID == FLUID_ZONE_ID) *rr = C2*Pre*exp(-E1/(R*T_prim)); else *rr=0.; } PS:C2應(yīng)該是摩爾分?jǐn)?shù),這里還是質(zhì)量分?jǐn)?shù)沒有改過來,依然沒有結(jié)果。 |
木蟲 (小有名氣)
新蟲 (初入文壇)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|