| 查看: 4964 | 回复: 11 | ||||
[交流]
【求助】fluent计算中，udf编译通过，初始化出错 已有2人参与
|
||||
|
麻烦达人老兄们救急，这个程序是用来定义出口处的压力的，编译成功，但是在边界条件里选中此UDF后，在初始化时报错，不知是何缘故，请各位不吝赐教，拜谢！ 程序为：#include "udf.h" DEFINE_PROFILE(pressure_outlet,thread,position) { real x[ND_ND]; real y; real a; face_t f; cell_t c; begin_f_loop(f,thread) { F_CENTROID(x,f,thread); y=x[1]; a=F_U(f,thread); F_PROFILE(f,thread,position)=1000*9.81*(1-y/1000)-0.5*1000*a*a; } end_f_loop(f,thread) } 错误提示为： Error: FLUENT received fatal signal (ACCESS_VIOLATION) 1. Note exact events leading to error. 2. Save case/data under new name. 3. Exit program and restart to continue. 4. Report error to your distributor. Error Object: () |
回复此楼» 收录本帖的淘帖专辑推荐
 多相流UDF |
» 猜你喜欢
梦想
已经有3人回复
-
面上合作单位盖章
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有5人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有4人回复
-
基金不中，共勉
已经有7人回复
-
为什么 国际(地区)合作与交流项目 没有放榜？
已经有11人回复
-
让我中一个面上吧！
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在ANSYS14中的FLuent中，写关于流体比热物性的UDF
已经有14人回复
- udf编程通过，但是在FLUENT中运行出错，麻烦udf高手帮忙看看
已经有4人回复
- fluent 在udf编译时出'constant'错！
已经有6人回复
- FLUENT源项UDF非稳态迭代的一个疑问请教
已经有8人回复
- UDF载入fluent之后出错，求指点！
已经有7人回复
- 有关fluent的udf中的hook
已经有6人回复
- 【求助】单机多核fluent并行udf
已经有13人回复
- 【求助】linux版fluent使用udf并行问题请教
已经有4人回复
- 【求助】如何改变Fluent中的曳力模型（UDF）
已经有12人回复
onwardupward
金虫 (正式写手)
Dr.
- 应助: 8 (幼儿园)
- 金币: 589.8
- 散金: 66
- 红花: 2
- 帖子: 340
- 在线: 44.2小时
- 虫号: 1093727
- 注册: 2010-09-09
- 性别: GG
- 专业: 高分子材料的加工与成型
11楼2014-12-04 17:54:55
2楼2010-07-16 16:35:36
3楼2010-07-16 22:48:16
★
robert2020(金币+1):多谢虫友应助！鼓励继续解答！ 2010-09-06 18:56:51
尘积(金币+1): 2010-11-16 15:59:37
robert2020(金币+1):多谢虫友应助！鼓励继续解答！ 2010-09-06 18:56:51
尘积(金币+1): 2010-11-16 15:59:37
|
试了一下编辑可以通过：程序应该为 #include "udf.h" DEFINE_PROFILE(pressure_outlet,thread,position) { real x[ND_ND]; real y; real a; face_t f; begin_f_loop(f,thread) { F_CENTROID(x,f,thread); y=x[1]; a=F_U(f,thread); F_PROFILE(f,thread,position)=1000.*9.81*(1-y/1000.)-0.5*1000.*a*a; } end_f_loop(f,thread) } |
4楼2010-09-06 16:41:04
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我也遇到了这情况啊！lz大哥你就绝了没有，我是自定义的密度和粘度，编译没问题，初始化通不过，附上我的程序，求指导~ #include "udf.h" #define TC 305.25 #define PC 7380000 #define W 0.225 #define R 8.314 #define tc 305.25 #define pc 7375000 #define w 0.225 #define r 8.314 #define pc1 73.75 #define dc 468 #define m 44.01 #define v0 150 DEFINE_PROPERTY(cell_density, cell, thread) { real dd; real temp = C_T(cell, thread); real pres = C_P(cell, thread); real tr,k,a1,a2,a,b,a3,a4,a5,a6,b1,b2,b3,y1,y2,y3,x; tr = temp/TC; k = 0.37464+1.54226*W-0.26992*W*W; a1 = (1+k*(1-sqrt(tr)))*(1+k*(1-sqrt(tr))); a2 = 0.45724*R*R*TC*TC/PC; a = a1*a2; b=0.0778*R*TC/PC; a3=pres; a4=pres*b-R*temp; a5=a-3*pres*b*b-2*R*temp*b; a6=pres*b*b*b+R*temp*b*b-a*b; b1 = a4*a4-3*a3*a5; b2 = a4*a5-9*a3*a6; b3 = a5*a5-3*a4*a6; y1 = b1*a4+3*a3*(-b2+sqrt(b2*b2-4*b1*b3))/2; y2 = b1*a4+3*a3*(-b2-sqrt(b2*b2-4*b1*b3))/2; if (y1<0) {y3=-y1;} else {y3=y1;} x = (-a4-(pow(y3,0.33333)-pow(-y2,0.33333)))/(3*a3); dd =44/(1000*x); return dd; } DEFINE_PROPERTY(cell_viscosity, cell, thread) { real v; real temp = C_T(cell, thread); real pres = C_P(cell, thread); real tr,k,a1,a2,a,b,a3,a4,a5,a6,b1,b2,b3,y1,y2,y3,x,dd,b4,b5,dr; tr=temp/tc; k=0.37464+1.54226*w-0.26992*w*w; a1=(1+k*(1-sqrt(tr)))*(1+k*(1-sqrt(tr))); a2=0.45724*r*r*tc*tc/pc; a=a1*a2; b=0.0778*r*tc/pc; a3=pres; a4=pres*b-r*temp; a5=a-3*pres*b*b-2*r*temp*b; a6=pres*b*b*b+r*temp*b*b-a*b; b1=a4*a4-3*a3*a5; b2=a4*a5-9*a3*a6; b3=a5*a5-3*a4*a6; y1=b1*a4+3*a3*(-b2+sqrt(b2*b2-4*b1*b3))/2; y2=b1*a4+3*a3*(-b2-sqrt(b2*b2-4*b1*b3))/2; if (y1<0) {y3=-y1;} else {y3=y1;} x=(-a4-(pow(y3,0.333333)-pow(-y2,0.333333)))/(3*a3); dd=44/(1000*x); b4=pow(tc,0.166666)/(sqrt(m)*pow(pc1,0.666666)); dr=dd/dc; b5=1.0230+0.23364*dr+0.058533*pow(dr,2)-0.40758*(dr,3)+0.093324*pow(dr,4); v=(((pow(b5,4)-1)/b4)+v0)*0.000001; return v; } 错误提示： Error: FLUENT received fatal signal (ACCESS_VIOLATION) 1. Note exact events leading to error. 2. Save case/data under new name. 3. Exit program and restart to continue. 4. Report error to your distributor. Error Object: () |
5楼2012-07-18 10:30:31
leixiaobo713
银虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1228.5
- 散金: 148
- 红花: 4
- 帖子: 490
- 在线: 121.9小时
- 虫号: 1997539
- 注册: 2012-09-13
- 性别: GG
- 专业: 金属功能材料
6楼2013-09-04 11:23:26
gengxiaowei
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1752.8
- 帖子: 35
- 在线: 91.4小时
- 虫号: 2051553
- 注册: 2012-10-09
- 性别: GG
- 专业: 金属材料的力学行为
7楼2013-09-08 14:51:56
gengxiaowei
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1752.8
- 帖子: 35
- 在线: 91.4小时
- 虫号: 2051553
- 注册: 2012-10-09
- 性别: GG
- 专业: 金属材料的力学行为
8楼2013-09-08 14:53:49
死天
专家顾问 (著名写手)
-
专家经验: +60 - 应助: 399 (硕士)
- 金币: 3648.9
- 红花: 89
- 帖子: 1336
- 在线: 443.8小时
- 虫号: 2643057
- 注册: 2013-09-09
- 专业: 流体力学
- 管辖: 仿真模拟
9楼2013-09-09 16:14:37
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
我也遇到了一样的问题，我要把回水口的温度加上2℃，然后赋值给入口。程序如下： #include "udf.h" DEFINE_PROFILE(inlet_t_1, inlet, index) { Domain *d=Get_Domain(1); Thread *out = Lookup_Thread(d, 18); face_t face; real t_out=0.0, area_out=0.0; real NV_VEC(area); begin_f_loop(face, out) { F_AREA(area, face, out); t_out += F_T(face, out)*NV_MAG(area); area_out += NV_MAG(area); } end_f_loop(face, out) t_out = t_out/area_out+2.0; begin_f_loop(face, inlet) { F_PROFILE(face, inlet, index) = t_out; } end_f_loop(face, inlet) } 然后设置边界条件的时候选中这个udf，就出现跟楼主一样的error，你是怎么解决的啊 |
10楼2013-12-12 12:05:24
