找回密码
 注册
查看: 1270|回复: 1

请大神帮我看一下这个多相流的DEFINE_SOURCE宏有什么问题

[复制链接]
发表于 2013-3-28 17:21:13 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?注册

x
大家好我想做MIxture模型溶解氧的数值模拟,这是我写的UDF,不能通过编译,C语言苦手死活找不出错误,望大神帮我
#include <udf.h>
#define Dl 2.5e-9/*定义分子扩散系数*/
#define H 4.40e6/*定义氧气亨利常数*/
#defune r 0.005/*假设气泡大小保持不变*/
DdEFINE_SOURCE(DO,c,t,ds,eqn)
{
        real area,Kl,urelx,urely,Vr,Cs;/*分别定义单位体积气液表面积,传质系数,相对速度,当地氧气饱和浓度*/
        Thread **pt = THREAD_SUB_THREADS(t);
Thread *tl = pt[0];
Thread *tg = pt[1];/*分别定义气液相的指针*/
                urelx = C_U(c,tg) - C_U(c,tl);
urely = C_V(c,tp) - C_V(c,tl);
Vr = sqrt(urelx*urelx + urely*urely);/*求出气液相相对速度*/
Cs=0.21C_P(c,t)/H;/*由亨利定律求出当地饱和氧气浓度*/
area=3C_VOF(c,tg)/r/*表面积与体积分数的关系*/
Kl=2sqrt(Dl*Vr/(2*3.14*r));
source=Kl*area*(Cs-DO);
dS=0
return source;        }
 楼主| 发表于 2013-3-28 17:21:45 | 显示全部楼层

回复 1# jiandaolaoda 的帖子

是不是气液相指针这样定义不对啊
您需要登录后才可以回帖 登录 | 注册

本版积分规则

快速回复 返回顶部 返回列表