24小时热门版块排行榜

返回列表

feiwen90

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 2
帖子: 4
在线: 15.2小时
虫号: 2264473
注册: 2013-01-28
专业: 数值预报与数值模拟

[求助] 【求助】压力出口背压值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 的主题更新

返回列表