| 查看: 3343 | 回复: 19 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
能量源项加载问题
|
||||
|
坐标原点为体热源中心，加载高斯旋转体热源；程序如下，但是热源加载不上，求指点 #include "udf.h" DEFINE_SOURCE(heat_flux,cell,thread,ds,eqn) { real x[ND_ND]; real oq,t0,v0,Q,time,cs,hh; real r,rh; real source_heat; cell_t c; C_CENTROID(x,c,thread); Q=100000; /*枪功率*/ cs=3/pow(oq,2); /*形状因子*/ hh=150e-3; /*热源高度*/ oq=60e-3; /*热源半径*/ r=sqrt(pow(x[0]-0.35,2)+pow(x[1]+0.35,2)); /*坐标与热源中心的距离*/ rh=sqrt(log(hh/x[2])/3.)*oq; /*一定热源高度对应的半径rh*/ if (x[2]==0) { source_heat=3.*cs*Q/(3.14*hh*(1.-1./exp(3))); ds[eqn] =0; } else if (x[2]>0&&x[2]<=hh&&r<=rh) { source_heat=3.*cs*Q/(3.14*hh*(1.-1./exp(3)))*exp(-3.*cs/log10(hh/x[2])*pow(r,2.)); ds[eqn] =0; } else { source_heat= ds[eqn] =0; } return source_heat; } 参考信息：http://muchong.com/bbs/viewthread.php?tid=4198673 |
回复此楼» 收录本帖的淘帖专辑推荐
 多相流UDF |
» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有10人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- fluent中的体热源加载问题（实际上是热流密度，不是热生成率）
已经有19人回复
- 求助高斯能量计算问题
已经有6人回复
- 关于过渡态的计算的能量问题。
已经有7人回复
- 请教个能量均分问题
已经有6人回复
- 能量方程源项udf
已经有4人回复
- 关于过渡态和能量优化的问题
已经有7人回复
- 体系绝对能量的问题
已经有4人回复
- 【求助】求助或讨论一下关于能量传递源项中各个焓值的含义
已经有3人回复
- 【求助】关于Gaussian计算中的能量问题
已经有14人回复
- 几个合金能量学方面的问题
已经有3人回复
17楼2012-09-08 14:33:36
2楼2012-08-20 10:01:27
3楼2012-08-20 11:05:47
4楼2012-08-20 11:45:43
