| 查看: 578 | 回复: 0 | ||
[求助]
udf中SOx代码
|
|
if (POLLUT_EQN(Pollut_Par) == EQ_SO2) 这是fluent的udf中做SOx，这个表示什么意思？没看懂。下面是完整的 static void so2_so3_rate(cell_t c, Thread* t, Pollut_Cell *Pollut, Pollut_Parameter *Pollut_Par, SOx_Parameter *SOx) {real kf,kr,rf=0,rr=0; real o_eq; real r_volatile,Ys_volatile,fuels_so2_frac; Rate_Const K_F = {1.2e6, 0.0, 39765.575}; Rate_Const K_R = {1.0e4, -1.0, 10464.625}; Rate_Const K_O = {36.64, 0.5, 27123.0}; /* SO3 + O <-> SO2 + O2 */ kf = ARRH(Pollut, K_F); kr = ARRH(Pollut, K_R); o_eq = ARRH(Pollut, K_O)*sqrt(MOLECON(Pollut, O2)); if (POLLUT_EQN(Pollut_Par) == EQ_SO2) { int ifstream; Ys_volatile = 1.e-04; fuels_so2_frac = 1.; for(ifstream=0; ifstream<Pollut_Par->nfstreams; ifstream++) { rf += Pollut->r_fuel_gls[ifstream]*SOx->Ys_fuelvolat[ifstream] *SOx->fuels_so2_frac[ifstream]*1000./Pollut_Par->sp[S].mw; } rf += kf*o_eq*MOLECON(Pollut, IDX(SO3)); rr = -kr*MOLECON(Pollut, O2)*MOLECON(Pollut, IDX(SO2)); } else if (POLLUT_EQN(Pollut_Par) == EQ_SO3) { rf = kr*MOLECON(Pollut, O2)*MOLECON(Pollut, IDX(SO2)); rr = -kf*o_eq*MOLECON(Pollut, IDX(SO3)); } POLLUT_FRATE(Pollut) += rf; POLLUT_RRATE(Pollut) += rr; } |
回复此楼» 猜你喜欢
长年满屏的广告，版主太不责任了。基金也等的急
已经有3人回复
-
filecode
已经有12人回复
-
不应该看fileCode
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
