| 查看: 2115 | 回复: 1 | ||
[求助]
udf能读入，但是一初始化计算就闪退！
|
|
用udf编写了如下程序，希望达到读取速度分量作为大涡模拟的入口条件，可是每次一初始化就闪退，求问各位大神问题出在哪里？ # include <udf.h> DEFINE_PROFILE(vel_X,t,index) { float V_x; FILE *fp; char ch; face_t f; int i,j; int count=1; fp=fopen("E:\LHC\LES\velocity_U","r"  ;if(fp==NULL){ printf("Can not open file!\n" ; } begin_f_loop(f,t) { int step=RP_Get_Integer("time-step" ;for(i=0;i<count;i++){ for(j=0;j<step+1;j++){ fscanf(fp,"%lf%c",&V_x,&ch); } } F_PROFILE(f,t,index)=V_x; count++; } end_f_loop(f,t); if(fclose(fp)){ printf("Can not close the file!\n" ;} } //V_y DEFINE_PROFILE(vel_y,t,index) { float V_y; FILE *fp; char ch; face_t f; int i,j; int count=1; fp=fopen("E:\LHC\LES\velocity_V","r" ;if(fp==NULL){ printf("Can not open file!\n" ; } begin_f_loop(f,t) { int step=RP_Get_Integer("time-step" ;for(i=0;i<count;i++){ for(j=0;j<step+1;j++){ fscanf(fp,"%lf%c",&V_y,&ch); } } F_PROFILE(f,t,index)=V_y; count++; } end_f_loop(f,t); if(fclose(fp)){ printf("Can not close the file!\n" ;} } //V_z DEFINE_PROFILE(vel_z,t,index) { float V_z; FILE *fp; char ch; face_t f; int i,j; int count=1; fp=fopen("E:\LHC\LES\velocity_W","r" ;if(fp==NULL){ printf("Can not open file!\n" ; } begin_f_loop(f,t) { int step=RP_Get_Integer("time-step" ;for(i=0;i<count;i++){ for(j=0;j<step+1;j++){ fscanf(fp,"%lf%c",&V_z,&ch); } } F_PROFILE(f,t,index)=V_z; count++; } end_f_loop(f,t); if(fclose(fp)){ printf("Can not close the file!\n" ;} } 非常感谢！！！ |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Fluent模拟，同样情况下 二维和三维为什么会存在很大的差别
已经有8人回复
- UDF关于两相流的燃烧的大涡模拟的相关程序
已经有8人回复
- 关于涡量-流函数法在fluent中的使用
已经有12人回复
- fluent 怎样使用柱坐标？
已经有10人回复
- 请教下可以用UDF编写第二类边界条件么？
已经有4人回复
- Fluent的UDF是否可以自定义辐射模型或湍流模型
已经有11人回复
- 请问fluent中能否实现二维的大涡模拟？
已经有15人回复
- 气液界面处两相速度不一样时，需要把气液剪切力作为动量源项吗？UDF编写气液剪切力
已经有14人回复
- Fluent 壁面传热边界条件
已经有3人回复
- Fluent中若干疑问
已经有5人回复
- 大涡模拟流场计算
已经有14人回复
- fluent时间步长的控制
已经有4人回复
- 求高手：三维FLUENT模拟到底有多难，难到哪里？
已经有25人回复
- 【求助】请问Fluent中如何简单地添加流动区域各坐标下的动量源项？
已经有16人回复
2楼2015-05-27 16:13:27
