| 查看: 4967 | 回复: 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 |
» 猜你喜欢
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有6人回复
-
国自然面上复盘~欢迎讨论
已经有10人回复
-
有多少人是今天查系统知道结果的？
已经有17人回复
-
基金未中，这种答复是模板吗？
已经有6人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有6人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有6人回复
-
青基如何破局
已经有4人回复
-
基金不中，共勉
已经有8人回复
-
怎么看青基中了没有啊
已经有5人回复
-
梦想
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 在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人回复
★
小木虫: 金币+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
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
10