24小时热门版块排行榜    

Znn3bq.jpeg
查看: 4348  |  回复: 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: ()
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

尘积

银虫 (小有名气)

引用回帖:
Originally posted by sghudng at 2010-07-16 16:35:36:

乍一看没啥问题啊……只知道cell_t c;这一行是多余的,根本没用到,你是面循环,又没用到相邻网格,要它干嘛

编译能通过 表明语法没错  就是在初始化是出错
3楼2010-07-16 22:48:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

sghudng

铜虫 (小有名气)


robert2020(金币+1):多谢应助! 2010-07-16 18:39:29
尘积(金币+1): 2010-11-16 15:59:50
引用回帖:
Originally posted by 尘积 at 2010-07-16 14:41:36:
麻烦达人老兄们救急,这个程序是用来定义出口处的压力的,编译成功,但是在边界条件里选中此UDF后,在初始化时报错,不知是何缘故,请各位不吝赐教,拜谢!
程序为:#include "udf.h"
DEFINE_PROFILE ...

乍一看没啥问题啊……只知道cell_t c;这一行是多余的,根本没用到,你是面循环,又没用到相邻网格,要它干嘛
2楼2010-07-16 16:35:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangj25

金虫 (正式写手)


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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xingwuyan

新虫 (初入文坛)


小木虫: 金币+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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 求博导|生物质基多孔碳/超级电容方向,已有相关成果,寻能源材料/碳材料方向老师 +3 猪猪人Zzz 2026-04-12 3/150 2026-04-17 19:10 by 阳阳阳^_^
[考研] 一志愿华中农业071010,320求调剂 +16 困困困困坤坤 2026-04-14 18/900 2026-04-17 17:00 by zhanglin_84
[考研] 295分求调剂 +5 ?要上岸? 2026-04-17 5/250 2026-04-17 16:51 by fenglj492
[考研] 271求调剂 +37 2261744733 2026-04-11 39/1950 2026-04-17 10:11 by 黑科技矿业
[考研] 收到复试调剂但是去不了 +6 小蜗牛* 2026-04-16 6/300 2026-04-17 10:05 by 涵竹刘
[考研] 求调剂 +8 小聂爱学习 2026-04-16 10/500 2026-04-16 22:06 by 1shin_ichi
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 求助调剂,跨调 +19 X十甫寸Y 2026-04-11 20/1000 2026-04-15 21:18 by cuisz
[考研] 085801电气专硕272求调剂 +19 电气李 2026-04-13 21/1050 2026-04-15 13:37 by 黑科技矿业
[考研] 各位老师好,求调剂,本科211,一志愿天津大学生物与医药学硕,差两名录取。 +11 路六六jjj 2026-04-13 11/550 2026-04-14 16:01 by zs92450
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
[考研] 材料考研调剂 +29 云木达达 2026-04-11 31/1550 2026-04-13 13:32 by lyh鲁老师
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] +10 李多米lee. 2026-04-12 11/550 2026-04-12 22:58 by yuyin1233
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 一志愿厦大0856,306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 求调剂 +3 胃痉挛累了 2026-04-11 5/250 2026-04-11 14:13 by luhong1990
信息提示
请填处理意见