| 查看: 605 | 回复: 0 | ||
[求助]
编写了热源项S的UDF，计算结果和理论值相差很多
|
|
编写一个流体内部内热源的udf，没有报错，fluent也能正常运行，但是得出的结果和计算的不一样。 一个长方体流体区域（0<x<2m，801个节点，-0.005m<z<0，21个节点，y0.01m51个节点，单双精度没差多少），流体匀速流动，四面绝热，加这个内热源，五个热源项从小到大，但是出口温度却不是规律性变化的。 ufd如下 DEFINE_SOURCE(heat_gen,cell,thread,dS,eqn) { real x[ND_ND]; real a,b,l; real source; cell_t c; begin_c_loop(c,thread) { C_CENTROID(x,c,thread); a=x[2]; l=x[0]; b=-a; if (0.<l<2.0) { source=1450180/(2.82941+1000*b)/(2.82941+1000*b);/*1 /*1553160/(2.35057+1000*b)/(2.35057+1000*b);/*2 /*1446253/(2.1629+1000*b)/(2.1629+1000*b);/*3 /*1140643/(1.7078+1000*b)/(1.7078+1000*b);/*4 /*3*890782/(1.3403+1000*b)/(1.3403+1000*b);/*5*/ dS[eqn]=0; } else { source=0; dS[eqn]=0; } return source; } end_c_loop(c,thread) } 实在找不到问题在哪，请各位大神帮帮忙。 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
10