粮仓干燥

咬着笔头学软件

毕业论文是关于，通风过程中粮仓内水分变化。将仓内粮食设置为多孔介质区域，只编写了一个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在后面？