应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » FLUENT & CFX » 【求助】压力出口背压值udf的设定-udf贴出，各位大神快到碗里来
11/1返回列表
查看: 1511  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

feiwen90

新虫 (初入文坛)

[求助] 【求助】压力出口背压值udf的设定-udf贴出，各位大神快到碗里来

模型上面是气体，下面是堆积着一定高度的固体颗粒
想实现的功能：
1.搜索出出口处固气分界线，记边界线高度为Ymin
2.设置出口背压值：Ymin以下为固体，为防止坍塌，其背压值为计算出来的压力值，Ymin以上为气体，故其背压值为0.
同时写出来两个循环均是错误的，于是做了以下尝试
1，最简单原始模型：假设Ymin固定为2m，背压值2m以上为0,2m以下为100000。udf如下，可以正常运行：
#include"udf.h"
DEFINE_PROFILE(out_pressure,thread,position)
{
     real minY=2;                                          
     real x[ND_ND];
     real y;
     face_t f;
     begin_f_loop(f,thread)
     {                     
         F_CENTROID(x,f,thread);
         y = x[1];
         if (y<=minY)
         {
             F_PROFILE(f,thread,position)=100000;
         }
         else
         {
             F_PROFILE(f,thread,position)=0;
         }
     }
     end_f_loop(f,thread);
  }
  
2，在1的基础上加上对Ymin的判断（Ymin最大不能超过20，计算域高度为20），当空气相大于0.9时，即认为为空气相，udf编写如下。
#include"udf.h"
DEFINE_PROFILE(out2_pressure,thread,position)
{
     real minY=20;                                          
     real x[ND_ND];
     real y;
     face_t f;
     real vmax=0.9;
     real volume,minX;  
     cell_t c;
     Thread *t,*gast,*mix;
       begin_f_loop(f,thread)
     {  
        c = F_C0(f,thread);   /* returns ID for c0 */
        t = THREAD_T0(thread);
        mix=THREAD_SUPER_THREAD(t);
        gast=THREAD_SUB_THREAD(mix,0);
             C_CENTROID(x,c,t);
             volume=C_VOF(c,gast);
             if (volume > vmax )
             {
                 minX=x[1];
                 if (minX<minY)
                 {
                     minY=minX;
                 }
                 else
                 {   
                     minY=minY;
                 }
             }
             else
             {
                 minY=minY;
             }
     }
     end_f_loop(f,thread);
  
      
     begin_f_loop(f,thread)
     {                     
         F_CENTROID(x,f,thread);
         y = x[1];
         if (y<=minY)
         {
             F_PROFILE(f,thread,position)=100000;
         }
         else
         {
             F_PROFILE(f,thread,position)=0;
         }
     }
     end_f_loop(f,thread);
  }
但是build和load均成功，但一加载到fluent的压力出口边界上，就出现如下错误。
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: ()
  
3.不判断Ymin，Ymin还是固定值，但是Ymin下面的背压为计算值p，udf如下，出现了跟上面一样的错误
#include"udf.h"
DEFINE_PROFILE(out1_pressure,thread,position)
{
     real minY=2;                                          
     real x[ND_ND];
     real y,p;
     face_t f;
     real vmax=0.9;
     real vmin=0.1;
     real volume,minX;  
     cell_t c;
     Thread *t,*gast,*mix;
     begin_f_loop(f,thread)
     {                     
         c = F_C0(f,thread);   
         t = THREAD_T0(thread);
         mix=THREAD_SUPER_THREAD(t);
         gast = THREAD_SUB_THREAD(mix,0);
//        C_CENTROID(x,c,t);
         p = C_P(c,gast);
         F_CENTROID(x,f,thread);
         y = x[1];
         if (y<=minY)
         {
             F_PROFILE(f,thread,position)=100000;
         }
         else
         {
             F_PROFILE(f,thread,position)=0;
         }
//        Message("%f\n",p);
     }
     end_f_loop(f,thread);
  }
为fluent菜鸟，去网上搜索办法，基本上都是说要先初始化再加载udf（我也是这么做的），还有就是写的udf错误了。但是在找不到哪里错误，望高手指点啊···不胜感激，感激涕零！
回复此楼

» 猜你喜欢
1楼 2014-04-19 22:44:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 feiwen90 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-29 19:27 by cqwQDCxMcL3I
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-29 13:50 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
信息提示
关闭
请填处理意见
关闭