| 查看: 1268 | 回复: 5 | |||
[交流]
fluent udf 模拟出错 已有3人参与
|
|
#include "udf.h" #include "dpm.h" #include "surf.h" #include "math.h" DEFINE_DPM_HEAT_MASS(multiabs,p,Cp,hgas,hvap,cvap_surf,dydt,dzdt) { cell_t c=P_CELL(p); /* cell and thread*/这句有问题么 Thread *t=P_CELL_THREAD(p);/* where the particle is in*/这句有问题么 real Tp=P_T(p);/* particle temperature*/ real Tg=C_T(c,t); real mp=P_MASS(p);/* particle mass*/这句有问题么 real Dp=P_DIAM(p);/* particle diameter*/ real Ap=3.14*Dp*Dp;/*particle surface*/ real x_vel_g,x_vel_p,y_vel_g,y_vel_p,abs_v,slip_x,slip_y,reyp,rho_g,rho_p,mu_g,YI,pg,vol_c,D,Sc,Sh,kg,p_so2,pe_so2,Ns,mp_dot; x_vel_g=C_U(c,t)+1.e-8; y_vel_g=C_V(c,t)+1.e-8; x_vel_p=P_VEL(p)[0]; y_vel_p=P_VEL(p)[1]; slip_x=x_vel_g-x_vel_p; slip_y=y_vel_g-y_vel_p; rho_g=C_R(c,t)+1.e-8; rho_p=2.716; mu_g=C_MU_T(c,t)+1.e-8; YI=C_YI(c,t,0); pg=C_P(c,t)+1.e-8; vol_c=C_VOLUME(c,t)+1.e-8; /*compute slip*/ abs_v=sqrt(slip_x*slip_x+slip_y*slip_y); /*compute Reynold's number*/ reyp=rho_g*abs_v*Dp/mu_g; /* diffusion coefficient*/ D=1.e-7*pow(Tp,1.75)*pow(1.0/29+1.0/64,0.5)/(pg*pow(pow(41.8,1./3)+pow(20.1,1./3),2.)+1.e-8); /* Schmidt number*/ Sc=mu_g/(rho_g*D+1.e-8); Sh=2.+0.55*pow(reyp,0.5)*pow(Sc,1./3); kg=Sh*D/(Dp*8.314*Tp+1.e-8);/*solve kg*/ p_so2=rho_g/(rho_p+1.e-8)*YI*pg; pe_so2=103.9951*exp(-7.5778*Tg); Ns=kg/(Tg*8.314+1.e-8)*(p_so2-pe_so2); dydt[1]-=mp_dot;这里的1带表什么 dzdt->species[0]+=mp_dot;这里的0代表什么 } 气相组分一为so2，第二组分为空气  离散项为h2o 最后想加一句确定浓度的公式，公式为 dC/dt=-Ns*（单位体积内气液传质面积），dC为液相中，可以实现么  11.png |
回复此楼» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有8人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有7人回复
-
青基如何破局
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金不中，共勉
已经有10人回复
-
国自然面上复盘~欢迎讨论
已经有11人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fluent加载UDF保存时出错
已经有0人回复
- fluent模拟加载udf计算后总是出错
已经有7人回复
- FLUENT UDF求助
已经有5人回复
- Fluent双相变换热器数值模拟
已经有4人回复
- Fluent入口速度udf编写求指点。
已经有0人回复
- 求大神帮检查下一个简单的UDF
已经有13人回复
- udf读入txt作边界条件出错！求大神指点~
已经有4人回复
- fluent模拟时一调用UDF就出错，帮忙看看，急，谢谢了
已经有1人回复
- fluent空化模拟出错
已经有1人回复
- 稳定流+脉动流模拟，导入UDF函数出错
已经有1人回复
- 求问fluent大神：如何在Injection颗粒properties设置中选择UDF定义的质量流量函数
已经有2人回复
- fluent中UDF设定边界条件，导入出错
已经有0人回复
- UDF编译出错，求助？？
已经有2人回复
- fluent中UDF编译的问题！
已经有6人回复
- fluent模拟相变传热需要自定义材料物性随温度的变化吗？
已经有1人回复
- udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看
已经有4人回复
2楼2016-10-25 22:45:44
3楼2016-10-25 23:37:58
poptiger
银虫 (小有名气)
- 应助: 4 (幼儿园)
- 金币: 571.9
- 散金: 21
- 帖子: 121
- 在线: 18.8小时
- 虫号: 793516
- 注册: 2009-06-13
- 性别: GG
- 专业: 机械工程
4楼2016-10-28 20:01:08
zhangqi999
铜虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 483.1
- 红花: 1
- 帖子: 53
- 在线: 22小时
- 虫号: 4923260
- 注册: 2016-08-17
- 性别: GG
- 专业: 燃烧学
5楼2016-11-06 19:28:19
malasong9471
木虫 (正式写手)
- 应助: 58 (初中生)
- 金币: 7314
- 红花: 29
- 帖子: 639
- 在线: 117.1小时
- 虫号: 4100890
- 注册: 2015-09-25
- 性别: GG
- 专业: 工程热力学
6楼2016-11-06 22:07:31
