| 查看: 1265 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
请高手帮我看看我的关于体积力UDF中的错误
|
||
|
#include "udf.h" #define TSTART 0.0 #define d 0.00085 #define pi 3.14 #define a 3e-7 #define m0 20 #define R 0.0005 #define c 0.1 #define u 1.256e-6 DEFINE_DPM_BODY_FORCE(particle_body_force, p, i) { real VCP,cc,m,z0,r,x,y,z,bforce; if(P_TIME(p)>=TSTART) { x=P_POS(p)[0];/颗粒在X轴上的坐标 y=P_POS(p)[1];/颗粒在Y轴上的坐标 z=P_POS(p)[2];/颗粒在Z轴上的坐标 VCP=4.0*pi*a*a*a/3.0;/颗粒的体积 cc=c/(3+c); m=pow(m0,2); z0=z+R+d; r=sqrt(x*x+y*y+z0*z0); if(i==0) bforce=-9.0*u*cc*VCP*m*x*(4*z0*z0+r*r)/(16*pi*pi*(pow(r,10)));/X方向的体积力 else if(i==1) bforce=-9.0*u*cc*VCP*m*y*(4*z0*z0+r*r)/(16*pi*pi*(pow(r,10)));/Y方向的体积力 else if(i==2) bforce=-9.0*u*cc*VCP*m*(4*z0*z0*z0)/(16*pi*pi*(pow(r,10)));/Z方向的体积力 } else bforce=0.0; /* an acceleration should be returned */ return (bforce/P_MASS(p)); } [ Last edited by 王周鹏 on 2014-1-6 at 21:17 ] |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- UDF中如何编写体积分数的梯度
已经有9人回复
- 绝对吸附和超额吸附中—晶体的自由体积与气体的密度是如何确定的？请各位高手帮忙。
已经有8人回复
- udf计算两相流 获得混合区单元体积
已经有4人回复
- fluent 加载udf 失败 请高手不吝赐教
已经有8人回复
- fluent14编译UDF时总是出错，求高手啊！！！难道没人遇到过吗？？？
已经有11人回复
- Fluent关于DPM模型UDF编程的求助
已经有4人回复
- 【求助】求UDF高手帮忙看看我的问题
已经有18人回复
- 【求助】请问Fluent中如何简单地添加流动区域各坐标下的动量源项？
已经有16人回复
3楼2014-01-07 10:07:27
achilFes1990
木虫 (著名写手)
- 应助: 363 (硕士)
- 金币: 5307.5
- 红花: 21
- 帖子: 1400
- 在线: 197.1小时
- 虫号: 1115537
- 注册: 2010-10-07
- 性别: GG
- 专业: 海洋工程
2楼2014-01-07 08:08:18
