24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1606  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

huimengliu

新虫 (初入文坛)

[求助] fluent多孔介质阻力系数是位置的函数,但不管阻力系数大还是小,速度压力都没明显区别

一个流体域里包含了一个多孔域,编了一个udf,想要多孔域实现类似格栅的效果,隔一段阻力系数大,孔隙率接近0,空气很难通过,隔一段阻力系数小,孔隙率接近1,空气顺利通过。

采用ergun公式算的阻力系数,fluent多孔介质阻力系数是位置的函数,但不管阻力系数大还是小,速度压力都没明显区别  体积当量直径采用的是 6倍体积比表面积

采用的physical velocity 算的。

然后计算之后发现多孔域的速度压力云图,速度矢量图都没有体现出阻力系数大,阻力系数小的区别,但是却有压力降,似乎fluent自动给多孔域平均了各种阻力系数。

附上我的UDF:(有的时候孔隙率的设置能用UDF定义的,有的时候就float报错) 请大神指正,金币可以全部奉上。

#include "udf.h"
DEFINE_PROFILE(porosity,t,i)
{
   real xc[ND_ND];
   cell_t c;
   begin_c_loop(c,t)
      {
     if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=1;
    }
     else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=0;
      }
    }
end_c_loop(c,t)
}

DEFINE_PROFILE(v_resistance,t,i)
{
   real xc[ND_ND];
   cell_t c;
   begin_c_loop(c,t)
      {
     if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=0;
    }
     else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=664781714.9993;
      }
    }
end_c_loop(c,t)
}

DEFINE_PROFILE(v_resist1,t,i)
{
   real xc[ND_ND];
   cell_t c;
   begin_c_loop(c,t)
      {
     if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=0;
    }
     else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=664781714999.3;
      }
    }
end_c_loop(c,t)
}

DEFINE_PROFILE(i_resistance,t,i)
{
   real xc[ND_ND];
   cell_t c;
   begin_c_loop(c,t)
      {
     if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=0;
    }
     else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=7368209.3023;
      }
    }
end_c_loop(c,t)
}

DEFINE_PROFILE(i_resist1,t,i)
{
   real xc[ND_ND];
   cell_t c;
   begin_c_loop(c,t)
      {
     if (xc[0]>-44.1164 && xc[0]<-39.1164&&xc[0]>-54.1164 && xc[0]<-49.1164&&xc[0]>-64.1164 && xc[0]<-59.1164&&xc[0]>-74.1164 && xc[0]<-69.1164&&xc[0]>-84.1164 && xc[0]<-79.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=0;
    }
     else if (xc[0]>-39.1164 && xc[0]<-34.1164&&xc[0]>-49.1164 && xc[0]<-44.1164&&xc[0]>-59.1164 && xc[0]<-54.1164&&xc[0]>-69.1164 && xc[0]<-64.1164&&xc[0]>-79.1164 && xc[0]<-74.1164&&xc[1]>=-0.20 && xc[1]<=1.95&&xc[2]>=-3.575 &&xc[2]<=-3.4)
     {
       F_PROFILE(c,t,i)=7368209302.3;
      }
    }
end_c_loop(c,t)
}
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liukan12

银虫 (小有名气)

我看了一下 你的程序是不是用坐标来判断的 但是从书写方式上来看 你是不是没有给xc赋值呢? 是不是忘记 C_CENTROID了?

发自小木虫Android客户端
ANSYSFluent流体仿真研究人员。
3楼2018-07-15 12:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

雪呀嘎嘎

新虫 (初入文坛)

您好,我碰到和你一样的问题,请问你的问题解决了吗,可以加个扣扣交流一下吗

发自小木虫IOS客户端
2楼2018-07-13 16:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南开大学0710生物学359求调剂 +6 兔兔兔111223314 2026-03-29 8/400 2026-04-02 22:37 by louise0220
[考研] 一志愿a区211,085601-307分求调剂 +11 党嘉豪 2026-03-31 24/1200 2026-04-02 21:54 by yuq
[考研] 085600,320分求调剂 +6 大馋小子 2026-04-02 6/300 2026-04-02 21:54 by dongzh2009
[考研] 277求调剂 +3 12A3 2026-04-02 3/150 2026-04-02 21:30 by swb0218
[考研] 求调剂22408 288分 +4 new382 2026-04-02 4/200 2026-04-02 16:07 by 氮气气气
[考研] 一志愿同济大学323分(080500)求调剂 +6 yikeniu 2026-04-01 6/300 2026-04-02 14:19 by smileboy2006
[考研] 理学07化学 303求调剂 +16 睿08 2026-03-27 17/850 2026-04-02 13:29 by 七度不信任
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +4 哇呼哼呼哼 2026-04-02 4/200 2026-04-02 13:10 by yulian1987
[考研] 能源动力 调剂 +3 不破不立0 2026-04-02 3/150 2026-04-02 12:46 by ffffjjjj
[考研] 求调剂!生物与医药专硕 +4 逆转陆先生 2026-04-01 4/200 2026-04-02 11:51 by xiaoranmu
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 292求调剂 +14 是妍子也是研子 2026-03-30 14/700 2026-04-02 10:17 by sanrepian
[考研] 367求调剂 +8 芋泥啵啵… 2026-03-28 8/400 2026-04-02 09:20 by olim
[考研] 302求调剂一志愿北航070300,本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 材料科学与工程339求调剂 +11 hyz0119 2026-03-31 12/600 2026-04-01 18:40 by 伟大河北
[考研] 0703求调剂 +4 zizimo 2026-03-31 4/200 2026-04-01 16:04 by yanflower7133
[考研] 309求调剂 +19 谁不是少年 2026-03-29 19/950 2026-04-01 15:47 by jp9609
[考研] 土木304求调剂 +3 兔突突突, 2026-03-31 3/150 2026-04-01 09:42 by JourneyLucky
[考研] 105500药学求调剂,一志愿山东大学药学,348分 +3 gr哈哈哈 2026-03-28 3/150 2026-03-30 18:56 by 源_2020
信息提示
请填处理意见