|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?注册
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; } |
|