| 查看: 854 | 回复: 1 | ||
Torres看CFD银虫 (小有名气)
|
[求助]
粮仓模拟
|
|
各位大神，被水分源项一直困扰这。毕业论文是关于，通风过程中粮仓内水分变化，借鉴参考文献上说水分变化通过uds来实现的。现将仓内粮食设置为多孔介质区域，只编写了一个DEFINE_SOURCE(my_udf,,cell,thread,ds,eqn)行不通，具体如下： DEFINE_SOURCE(my_mass,c,t,dS,eqn) { real w,r,tem,k,we,p,psat; w=C_UDMI(c,t,0); 粮食初始水分 r=C_UDMI(c,t,1); 通风空气的初始相对湿度 x=***** 通风空气的湿含量 tem=C_T(c,t)-273.15; 读取网格温度 if(w!=0) { if(r<0.75) 相对湿度小于平衡相对湿度 { we=0.294-0.046*log(-(tem+35.703)*log(r)/log(2.7183))/log(2.7183); 粮食的平衡水分 psat=6.0e25/pow(Tabs,5) *exp(-6800/Tabs); 饱和水分蒸汽压 p=101325*x/(0.62198+); 粮堆内蒸汽压 r=p/psat; 相对湿度 k=2000*exp(-5094/(tem+273.15))*(w-we); w=w-k; x=x+k; } C_UDMI(c,t,1)=r; } else { w=w; } C_UDMI(c,t,0)=w; } 当相对湿度<0.75的时候，打算一直循环下去，但是自己查看udf源项的时候，说每次迭代都会调用一次udf，因此不知道该不该加上循环？而且这个方程采用的是unsteady模式，是不是该继续添加一个udf在后面？ |
回复此楼
王者归来.
木虫 (职业作家)
- 应助: 17 (小学生)
- 金币: 4014.3
- 红花: 4
- 帖子: 3802
- 在线: 117.6小时
- 虫号: 3542982
- 注册: 2014-11-17
- 性别: GG
- 专业: 经济学
2楼2014-12-14 13:39:26
