 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 4083 | 回復(fù): 9 | |||
huimengliu新蟲(chóng) (初入文壇)
|
[交流]
fluent多孔介質(zhì)孔隙率是位置的函數(shù)問(wèn)題求助 已有3人參與
|
|
想要給一塊板一段流體能通過(guò)，一段流體不能通過(guò)，我想通過(guò)定義孔隙率為0和1來(lái)實(shí)現(xiàn) 定義了一個(gè)流體域，一個(gè)多孔域 對(duì)多孔域使用了UDF #include "udf.h" DEFINE_PROFILE(Bed_porosity,t,i) { real xc[ND_ND]; cell_t c; begin_c_loop(c,t) { if (xc[0]>-606 && xc[0]<-600,xc[1]>320 && xc[1]<=470,xc[2]>0 &&xc[2]<100) { F_PROFILE(c,t,i)=1; } else if (xc[0]>-606 && xc[0]<-600,xc[1]>470 && xc[1]<620,xc[2]>0 &&xc[2]<100) { F_PROFILE(c,t,i)=0; } } end_c_loop(c,t) } 在fluent里編譯也通過(guò)了，但是fluent算完后發(fā)現(xiàn)多孔域流體全都能通過(guò)，并不是一段通過(guò)，一段不通過(guò)。 求大神指導(dǎo) |
回復(fù)此樓» 猜你喜歡
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
已經(jīng)有8人回復(fù)
-
277材料科學(xué)與工程080500求調(diào)劑
已經(jīng)有7人回復(fù)
-
317求調(diào)劑
已經(jīng)有18人回復(fù)
-
293求調(diào)劑
已經(jīng)有5人回復(fù)
-
280分求調(diào)劑 一志愿085802
已經(jīng)有7人回復(fù)
-
0854電子信息求調(diào)劑
已經(jīng)有3人回復(fù)
-
263求調(diào)劑
已經(jīng)有4人回復(fù)
-
石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告
已經(jīng)有3人回復(fù)
-
一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑
已經(jīng)有13人回復(fù)
-
北科281學(xué)碩材料求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 歐拉模型+多孔介質(zhì) 求助
已經(jīng)有4人回復(fù)
- 求助關(guān)于固定床多孔介質(zhì)孔隙率的問(wèn)題
已經(jīng)有7人回復(fù)
- fluent中UDF編譯的問(wèn)題！
已經(jīng)有11人回復(fù)
- 多孔介質(zhì)孔隙率udf問(wèn)題
已經(jīng)有7人回復(fù)
- fluent中多空介質(zhì)的計(jì)算
已經(jīng)有19人回復(fù)
- fluent多孔介質(zhì)設(shè)置問(wèn)題��！求助
已經(jīng)有10人回復(fù)
- 【求助】請(qǐng)問(wèn)在fluent中怎么設(shè)置多孔介質(zhì)的內(nèi)熱源啊
已經(jīng)有13人回復(fù)
- 多孔介質(zhì)孔隙率的聲明
已經(jīng)有4人回復(fù)
- 急求各位大俠啊~fluent里多孔介質(zhì)模型參數(shù)設(shè)置的問(wèn)題~
已經(jīng)有5人回復(fù)
- 多孔介質(zhì)在gambit中建模后邊界條件設(shè)置問(wèn)題
已經(jīng)有15人回復(fù)
- 急求��！請(qǐng)問(wèn)fluent里如何利用多孔介質(zhì)模型模擬兩相流吸附過(guò)濾？
已經(jīng)有8人回復(fù)
- Fluent多孔介質(zhì)模擬問(wèn)題
已經(jīng)有15人回復(fù)
- 【求助】求助啦~雙相介質(zhì)和孔隙介質(zhì)是一回事嗎？
已經(jīng)有3人回復(fù)
- 【求助】大俠支招，另急需FLUENT多孔介質(zhì)的算例
已經(jīng)有11人回復(fù)
★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
2樓2014-09-28 12:51:21
3樓2014-09-28 16:04:03
huimengliu
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 34.5
- 帖子: 8
- 在線: 5.4小時(shí)
- 蟲(chóng)號(hào): 3425038
- 注冊(cè): 2014-09-18
4樓2014-10-09 09:18:14
huimengliu
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 34.5
- 帖子: 8
- 在線: 5.4小時(shí)
- 蟲(chóng)號(hào): 3425038
- 注冊(cè): 2014-09-18
|
改完計(jì)算之后，在cfd post 里面看壓力云圖多孔板內(nèi)部還是全部都通過(guò)流體。。。不知道是不是應(yīng)該把粘性阻力系數(shù)，慣性阻力系數(shù)關(guān)于位置的函數(shù)也編入U(xiǎn)DF里面，因?yàn)榭紫堵首兓@兩個(gè)系數(shù)也會(huì)隨之變化。。。 或者是因?yàn)槲医o的xyz的坐標(biāo)不對(duì)應(yīng)多孔板的坐標(biāo)？我用UG建的模型，ICEM畫(huà)的網(wǎng)格，導(dǎo)入到fluent里面發(fā)現(xiàn)坐標(biāo)已經(jīng)略有變化。。是不是對(duì)多孔板這個(gè)多孔域進(jìn)行UDF編程有單獨(dú)的一套新的坐標(biāo)系？ 我有看到有人說(shuō)fluent自帶的多孔模型不考慮流體在多孔介質(zhì)內(nèi)部的蜿蜒，不考慮細(xì)觀分布，只是大體的等效，如果真是這樣那我再如何設(shè)置孔隙率是位置的函數(shù)都沒(méi)用吧 |
5樓2014-10-09 17:00:20
6樓2014-10-10 14:29:07
huimengliu
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 34.5
- 帖子: 8
- 在線: 5.4小時(shí)
- 蟲(chóng)號(hào): 3425038
- 注冊(cè): 2014-09-18
送紅花一朵 |
您好，我一直調(diào)試，看了一些別的例子，似乎要加上動(dòng)量源項(xiàng)。。下面是我看到的一個(gè)二維多孔介質(zhì)用的udf。。。 #include "udf.h" DEFINE_SOURCE(xmom_source,cell,thread,dS,eqn) { const real c2=100; real x[ND_ND]; real con,source; C_CENTROID(x,c,t); con=c2*0.5*C_R(c,t);/*C_R(c,t) terms density of air*/ source=0.001-con*fabs(C_U(c,t))*C_U(c,t); dS[eqn]=-2*con*fabs(C_U(c,t)); return source; } |
7樓2014-11-25 09:26:46
★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
本帖內(nèi)容被屏蔽 |
8樓2014-11-25 10:16:01
huimengliu
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 34.5
- 帖子: 8
- 在線: 5.4小時(shí)
- 蟲(chóng)號(hào): 3425038
- 注冊(cè): 2014-09-18
|
謝謝您的解答。 我模擬的過(guò)程中發(fā)現(xiàn)多孔域內(nèi)部和外部的速度一樣，雖然我確實(shí)發(fā)現(xiàn)了透過(guò)多孔介質(zhì)有壓力降。本來(lái)我定義的阻力系數(shù)大的地方速度應(yīng)該小，阻力系數(shù)小的得放速度應(yīng)該大。 我在文獻(xiàn)中發(fā)現(xiàn) 由于數(shù)值模擬中沒(méi)有表示多孔介質(zhì)區(qū)域的實(shí)際幾何模型存在，因此 FLUENT 中認(rèn)為在多孔介質(zhì)區(qū)域內(nèi)的速度與區(qū)域外氣流的速度相同。這樣無(wú)法 模擬多孔介質(zhì)區(qū)域?qū)饬魉俣仍黾拥挠绊��！斗里L(fēng)柵氣動(dòng)性能的風(fēng)洞試驗(yàn)及數(shù)值模擬研究》_李兆楊 于是 在標(biāo)準(zhǔn)動(dòng)量方程中增加動(dòng)量源項(xiàng)，我們可以通過(guò)在 FLUENT 中采用 UDF 編譯的方法直接計(jì)算多孔介質(zhì)區(qū)域的動(dòng)量方程源項(xiàng)。由于疏透結(jié)構(gòu)自身對(duì)氣流 的體積阻塞作用，會(huì)對(duì)氣流產(chǎn)生局部加速作用，如前所述，多孔介質(zhì)模型未能 考慮這種局部加速作用，我們可以通過(guò)考慮單元平均體積疏透率，通過(guò) UDF 編 譯程序來(lái)求解真實(shí)速度或物理速度————《防風(fēng)柵氣動(dòng)性能的風(fēng)洞試驗(yàn)及數(shù)值模擬研究》_李兆楊 附上他的UDF #include "udf.h" #define C2 20112500e+001 /*porous loss coefficients*/ #define G 8.200000e-001 /*porous stands for the value per unit*/ DEFINE_ADJUST(velocity_adjust,domain) { int ID=16; /* zone ID for porous zone */ Thread *thread=Lookup_Thread[domain,16]; cell_t,c; real x[ND_ND];/*this will hold the position vector*/ begin_c_loop(c,thread) { C_CENTROID(x,c,thread); C_U(c,thread)=G*C_U(c,thread); C_V(c,thread)=G*C_V(c,thread); C_W(c,thread)=G*C_W(c,thread); } end_c_loop(c,thread) } /*source define for x-momentum*/ DEFINE_SOURCE(xmom_source,c,t,dS,eqn) { real x[ND_ND];/*this will hold the position vector*/ real con,source; C_CENTROID(x,c,t); con=c2*0.5*C_R(c,t);/*C_R(c,t) terms density of air*/ source=0.001-con*fabs(C_U(c,t))*C_U(c,t); dS[eqn]=-2*con*fabs(C_U(c,t)); return source |
9樓2014-11-25 10:50:04
xingfuww
專(zhuān)家顧問(wèn) (正式寫(xiě)手)
-
專(zhuān)家經(jīng)驗(yàn): +78 - 仿真EPI: 2
- 應(yīng)助: 46 (小學(xué)生)
- 金幣: 6630.4
- 散金: 52
- 紅花: 29
- 帖子: 542
- 在線: 283.8小時(shí)
- 蟲(chóng)號(hào): 1776059
- 注冊(cè): 2012-04-24
- 性別: GG
- 專(zhuān)業(yè): 化學(xué)工程及工業(yè)化學(xué)
- 管轄: 仿真模擬
10樓2015-02-06 15:49:35
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 huimengliu 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
