| 查看: 1688 | 回复: 1 | ||
[求助]
UDS方程中的对流项
|
 现在有一个UDS方程中的对流项如图所示，其中Wf为常数，取0.3，我编辑的UDF代码如下，但是导入FLUENT后计算不收敛，不知道是不是对流项出错了，求各位大神指导！！刚注册，没什么金币，就靠大家帮忙了！ DEFINE_UDS_FLUX(snow_uds_flux,f,t,i) { real U[ND_ND],A[ND_ND],flux; cell_t c0,c1; Thread *t0,*t1; c0=F_C0(f,t); t0=F_C0_THREAD(f,t); F_AREA(A,f,t); flux=0.0; c1=-1; t1=NULL; if(BOUNDARY_FACE_THREAD_P(t)) { U[0]=F_U(f,t); U[1]=F_V(f,t); U[2]=(F_W(f,t)+0.3); flux=U[0]*A[0]+U[1]*A[1]+U[2]*A[2]; } else { c1=F_C1(f,t); t1=F_C1_THREAD(f,t); U[0]=(C_U(c0,t0)+C_U(c1,t1))/2.0; U[1]=(C_V(c0,t0)+C_V(c1,t1))/2.0; U[2]=(C_W(c0,t0)+0.3+C_W(c1,t1)+0.3)/2.0; flux=U[0]*A[0]+U[1]*A[1]+U[2]*A[2]; } return flux; } 采用NV_宏定义的向量形式编辑如下： DEFINE_UDS_FLUX(snow_uds_flux,f,t,i) { real NV_VEC(psi_vec),NV_VEC(A),flux; real U0,V0,W0,U1,V1,W1; cell_t c0,c1; Thread *t0,*t1; c0=F_C0(f,t); t0=F_C0_THREAD(f,t); F_AREA(A,f,t); flux=0.0; c1=-1; t1=NULL; if(BOUNDARY_FACE_THREAD_P(t)) { real dens; if (NNULLP(THREAD_STORAGE(t,SV_DENSITY))) dens=F_R(f,t); else dens=C_R(c0,t0); U0=F_U(f,t); V0=F_V(f,t); W0=(F_W(f,t)+0.3); NV_DS(psi_vec,=,U0,V0,W0,*,dens); flux=NV_DOT(psi_vec,A); } else { c1=F_C1(f,t); t1=F_C1_THREAD(f,t); U0=C_U(c0,t0); V0=C_V(c0,t0); W0=(C_W(c0,t0)+0.3); U1=C_U(c1,t1); V1=C_V(c1,t1); W1=(C_W(c1,t1)+0.3); NV_DS(psi_vec,=,U0,V0,W0,*,C_R(c0,t0)); NV_DS(psi_vec,+=,U1,V1,W1,*,C_R(c1,t1)); flux=NV_DOT(psi_vec,A)/2.0; } return flux; } 这两种方法是不是一样的呢？ |
回复此楼» 猜你喜欢
关于Filecode分析方法
已经有5人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有8人回复
-
不应该看fileCode
已经有6人回复
-
好奇怪的filecode
已经有6人回复
-
FileCode能看出啥？
已经有21人回复
-
长年满屏的广告，版主太不责任了。基金也等的急
已经有4人回复
-
filecode
已经有15人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有35人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 多孔介质udf，uds编写问题
已经有33人回复
- uds方程中二阶导数边界条件，编程如下，运行时总是出错，求指导，谢谢
已经有6人回复
- uds方程中如何定义饱和浓度
已经有4人回复
- fluent求解UDS过程中，出现少数不合理的值..
已经有14人回复
- FLUENT中耦合求解两个标量方程，求指点...
已经有5人回复
- 对于uds的对流项设置问题~
已经有4人回复
- fluent-udf自定义标量方程边界条件
已经有5人回复
- fluent里对于轴边界是怎样处理的？对于uds解方程怎样体现模型是二维轴对称模型呢？
已经有10人回复
- 有没有什么办法能控制uds求解的标量在整个求解域上的范围啊
已经有6人回复
2楼2017-03-19 16:17:13
