| 10 | 1/1 | 返回列表 |
| 查看: 4084 | 回復(fù): 9 | |||
[交流]
fluent多孔介質(zhì)孔隙率是位置的函數(shù)問題求助 已有3人參與
|
|
本帖內(nèi)容被屏蔽 |
|
改完計算之后,在cfd post 里面看壓力云圖多孔板內(nèi)部還是全部都通過流體。。。不知道是不是應(yīng)該把粘性阻力系數(shù),慣性阻力系數(shù)關(guān)于位置的函數(shù)也編入UDF里面,因為孔隙率變化這兩個系數(shù)也會隨之變化。。。 或者是因為我給的xyz的坐標(biāo)不對應(yīng)多孔板的坐標(biāo)?我用UG建的模型,ICEM畫的網(wǎng)格,導(dǎo)入到fluent里面發(fā)現(xiàn)坐標(biāo)已經(jīng)略有變化。。是不是對多孔板這個多孔域進行UDF編程有單獨的一套新的坐標(biāo)系? 我有看到有人說fluent自帶的多孔模型不考慮流體在多孔介質(zhì)內(nèi)部的蜿蜒,不考慮細觀分布,只是大體的等效,如果真是這樣那我再如何設(shè)置孔隙率是位置的函數(shù)都沒用吧 |
送紅花一朵 |
您好,我一直調(diào)試,看了一些別的例子,似乎要加上動量源項。。下面是我看到的一個二維多孔介質(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; } |
|
本帖內(nèi)容被屏蔽 |
|
謝謝您的解答。 我模擬的過程中發(fā)現(xiàn)多孔域內(nèi)部和外部的速度一樣,雖然我確實發(fā)現(xiàn)了透過多孔介質(zhì)有壓力降。本來我定義的阻力系數(shù)大的地方速度應(yīng)該小,阻力系數(shù)小的得放速度應(yīng)該大。 我在文獻中發(fā)現(xiàn) 由于數(shù)值模擬中沒有表示多孔介質(zhì)區(qū)域的實際幾何模型存在,因此 FLUENT 中認為在多孔介質(zhì)區(qū)域內(nèi)的速度與區(qū)域外氣流的速度相同。這樣無法 模擬多孔介質(zhì)區(qū)域?qū)饬魉俣仍黾拥挠绊憽!斗里L(fēng)柵氣動性能的風(fēng)洞試驗及數(shù)值模擬研究》_李兆楊 于是 在標(biāo)準動量方程中增加動量源項,我們可以通過在 FLUENT 中采用 UDF 編譯的方法直接計算多孔介質(zhì)區(qū)域的動量方程源項。由于疏透結(jié)構(gòu)自身對氣流 的體積阻塞作用,會對氣流產(chǎn)生局部加速作用,如前所述,多孔介質(zhì)模型未能 考慮這種局部加速作用,我們可以通過考慮單元平均體積疏透率,通過 UDF 編 譯程序來求解真實速度或物理速度————《防風(fēng)柵氣動性能的風(fēng)洞試驗及數(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 |
專家顧問 (正式寫手)
![]() |
專家經(jīng)驗: +78 |
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +11 | 梨花珞晚風(fēng) 2026-03-17 | 11/550 |
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
[考研] 354求調(diào)劑 +7 | Tyoumou 2026-03-18 | 10/500 |
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 4/200 |
|
|
[考研] 333求調(diào)劑 +5 | 87639 2026-03-21 | 7/350 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 265求調(diào)劑 +9 | 梁梁校校 2026-03-17 | 9/450 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 材料專碩英一數(shù)二306 +7 | z1z2z3879 2026-03-18 | 7/350 |
|
|
[考研] 材料與化工專碩調(diào)劑 +7 | heming3743 2026-03-16 | 7/350 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|