|
|
[交流]
能量源项加载问题
坐标原点为体热源中心，加载高斯旋转体热源；程序如下，但是热源加载不上，求指点
#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 |
» 收录本帖的淘帖专辑推荐
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有5人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有7人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
» 抢金币啦！回帖就可以得到:
查看全部散金贴
|
回复此楼
+5/4015
看来今天不会放榜了？
+8