24小时热门版块排行榜

返回列表

凯仂

铜虫 (小有名气)

应助: 3 (幼儿园)
金币: 79.8
帖子: 66
在线: 36.2小时
虫号: 1983930

[交流] FLUENT DUF DEFINE_INIT 导入不能初始化

fluent模拟流化床传热传质，编程设定初始的颗粒湿度，并计算求得空气的相对湿度，但是一导进去一点初始化就出现错误，DEFINE_INIT去掉以后可以初始化，求高手帮忙解决这个问题...
#include "udf.h"
/* Heat and moisture transfer in bulks of grains */
/* Graham Thorpe */
/* Victoria University, Melbourne, Australia */
/* Chung-Pfost constants */
#define ACP 921.69
#define BCP 18.077
#define CCP 112.35
#define RHOB 639.2
#define PATM 101325
/* Define ambient temperature conditions */
#define TMEAN 288.15
/* Define a global variable that acts as a comparitor in DEFINE_ADJUST */
int last_ts = -1;
DEFINE_INIT(moisure_init,d)
/* The initial solids moisture content, W, is set */
/* and the corresponding humidity, w, of the intergranular air is calculated */
{
cell_t c;
Thread *t;
real W,Tabs,TC,r,psat,p,w;
/* Loop over all of the cell threads in the bulk of grain */
thread_loop_c(t,d)
{
begin_c_loop(c,t)
{
C_UDMI(c,t,0) = 0.1364;
W = C_UDMI(c,t,0);
Tabs = C_T(c,t);
TC = Tabs-273.15;
/* Calculate relative humidity of intergranular air */
r = exp(-ACP/(TC+CCP)*exp(-BCP*W));
/* Ensure the air does not become supersaturated */
if (r>0.99)
{r = 0.99;}
psat = 6.0e25/pow(Tabs,5) *exp(-6800/Tabs);
p = r*psat;
w = 0.622*p/(PATM-p); /* Eqn 12 */
C_UDSI(c,t,0) = w;
}
end_c_loop(c,t);
}
}

_ISCN{O85`GZ`[81Y69R~}M.jpg

回复此楼