版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

feiwen90

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 2
帖子: 4
在線: 15.2小時
蟲號: 2264473
注冊: 2013-01-28
專業(yè): 數(shù)值預(yù)報與數(shù)值模擬

[求助] 【求助】壓力出口背壓值udf的設(shè)定-udf貼出，各位大神快到碗里來

模型上面是氣體，下面是堆積著一定高度的固體顆粒
想實(shí)現(xiàn)的功能：
1.搜索出出口處固氣分界線，記邊界線高度為Ymin
2.設(shè)置出口背壓值：Ymin以下為固體，為防止坍塌，其背壓值為計算出來的壓力值，Ymin以上為氣體，故其背壓值為0.
同時寫出來兩個循環(huán)均是錯誤的，于是做了以下嘗試
1，最簡單原始模型：假設(shè)Ymin固定為2m，背壓值2m以上為0,2m以下為100000。udf如下，可以正常運(yùn)行：
#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的基礎(chǔ)上加上對Ymin的判斷（Ymin最大不能超過20，計算域高度為20），當(dāng)空氣相大于0.9時，即認(rèn)為為空氣相，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的壓力出口邊界上，就出現(xiàn)如下錯誤。
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如下，出現(xiàn)了跟上面一樣的錯誤
#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菜鳥，去網(wǎng)上搜索辦法，基本上都是說要先初始化再加載udf（我也是這么做的），還有就是寫的udf錯誤了。但是在找不到哪里錯誤，望高手指點(diǎn)啊···不勝感激，感激涕零！

回復(fù)此樓

» 猜你喜歡

291求調(diào)劑已經(jīng)有7人回復(fù)
材料與化工272求調(diào)劑已經(jīng)有13人回復(fù)
283求調(diào)劑（080500）已經(jīng)有5人回復(fù)
材料科學(xué)與工程求調(diào)劑已經(jīng)有5人回復(fù)
一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑已經(jīng)有17人回復(fù)
考研調(diào)劑已經(jīng)有7人回復(fù)
343求調(diào)劑085601 已經(jīng)有3人回復(fù)
08工科求調(diào)劑286 已經(jīng)有3人回復(fù)
337求調(diào)劑已經(jīng)有3人回復(fù)
329求調(diào)劑，一志愿西北工業(yè)大學(xué)，材料工程（085601）已經(jīng)有10人回復(fù)

1樓 2014-04-19 22:44:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

智能機(jī)器人

Robot (super robot)

我們都愛小木蟲

找到一些相關(guān)的精華帖子，希望有用哦~

出口壓力的UDF編寫問題已經(jīng)有10人回復(fù)
全部身家求大神幫忙解決fluent壓力脈沖怎么編寫UDF的問題已經(jīng)有8人回復(fù)
Fluent 壓力出口邊界條件設(shè)置已經(jīng)有6人回復(fù)

點(diǎn)擊這里搜索更多相關(guān)資源

科研從小木蟲開始，人人為我，我為人人

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 feiwen90 的主題更新