 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3379 | 回復(fù): 12 | ||
BeObjective銀蟲 (正式寫手)
|
[求助]
關(guān)于UDF源項的編寫遇到了一下問題，有大神做過類似的東西的望指點 已有1人參與
|
|
|
首先介紹一下我的問題吧，做混合氣體（空氣和水蒸氣）的冷凝，網(wǎng)格如下圖所示： 二維網(wǎng)格，上部是進(jìn)口，下面是出口，左邊是恒溫壁面，右邊是絕熱壁面，蒸汽會在左邊的壁面上進(jìn)行冷凝，使用的是species模型，冷凝的過程是通過UDF添加質(zhì)量、能量源項的方式來實現(xiàn)的。公式如下圖： 這幾個是原論文中的質(zhì)量源項公式： 我自己的想法是這樣的：質(zhì)量源項出現(xiàn)的區(qū)域肯定是在靠近壁面的地方，而在這里的模型中忽略了液膜的厚度，那么靠近壁面的那一層網(wǎng)格就是源項需要處理的區(qū)域，于是通過FLUENT的Adapt>mark，標(biāo)記那一層網(wǎng)格，然后separate zone將整個流體域分成了兩個zone，這樣udf源項只在靠近壁面的那一層網(wǎng)格中去做。為了簡化問題，我先只做了質(zhì)量源項的部分，能量源項暫時不考慮，寫完程序如下： #include "udf.h" DEFINE_SOURCE(wall_msource,c,t,dS,eqn) { real source,temp,pres,dens; real m_f,diff; int i=0; temp=C_T(c,t); pres=C_P(c,t); dens=C_R(c,t); m_f=C_YI(c,t,i); C_UDMI(c,t,0)=C_YI_G(c,t,i)[0]; diff=(2.775e-6+4.479e-8*temp+1.656e-10*temp*temp)*pres/100; source=dens*diff*C_UDMI(c,t,0)/(1-m_f); return source; } 關(guān)于這個程序，有一點說明，論文中公式后面幾個是計算質(zhì)量分?jǐn)?shù)的，但是，udf能夠存儲水蒸汽質(zhì)量分?jǐn)?shù)，為什么論文中還多此一舉去求質(zhì)量分?jǐn)?shù)？這是我的一點疑問，于是我這里就直接獲得求解器的質(zhì)量分?jǐn)?shù)，并且獲得了質(zhì)量分?jǐn)?shù)的梯度，這樣問題就來了，在這個區(qū)域只有一層網(wǎng)格，獲得質(zhì)量分?jǐn)?shù)沒問題，但是這個質(zhì)量分?jǐn)?shù)梯度的獲得有沒有問題，按理說質(zhì)量分?jǐn)?shù)的梯度至少是兩個網(wǎng)格單元做的差分，一層網(wǎng)格獲得梯度貌似存在問題。這個程序?qū)懲晡覍?dǎo)入FLUENT進(jìn)行計算，就出現(xiàn)問題了，在殘差收斂后，我監(jiān)控了進(jìn)出口的質(zhì)量流量，發(fā)現(xiàn)這個源項產(chǎn)生的作用幾乎為0，進(jìn)出口完全持平，于是我為了排除是不是源項太小產(chǎn)生的誤差直接把質(zhì)量源項部分改寫成-0.4，得到非常顯著的進(jìn)出口流量差異。以下兩個圖就是UDF源項和constant-0.4源項所產(chǎn)生的水蒸氣質(zhì)量分?jǐn)?shù)分布云圖： 圖一是UDF的云圖，而且放大了很多： 圖二是常數(shù)源項的云圖： 兩者對比非常的明顯。圖一貌似沒有任何變化，圖二的分布很明顯。這兩個計算的CASE都已經(jīng)上傳上去了。 最近也是做這個東東做的有點麻木了，其中也做了很多交叉對比，沒有實質(zhì)性的進(jìn)展，看看有沒有有經(jīng)驗的大神一眼看穿問題出在哪里？只要確實能點到問題所在，不會吝嗇金幣，謝了！ 2.jpg 3.jpg 4.jpg 5.jpg 無標(biāo)題.jpg 1.jpg 2.jpg 3.jpg 4.jpg 5.jpg |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : fluent111.cas
- 附件 2 : fluent112.cas
2015-09-21 22:26:13, 694.51 K
2015-09-21 22:26:14, 694.61 K
» 猜你喜歡
一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分
已經(jīng)有5人回復(fù)
-
0856材料化工調(diào)劑 總分330
已經(jīng)有9人回復(fù)
-
085602 化工專碩 338分 求調(diào)劑
已經(jīng)有9人回復(fù)
-
352分 化工與材料
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有8人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有3人回復(fù)
-
材料求調(diào)劑一志愿哈工大324
已經(jīng)有4人回復(fù)
-
085600 286分 材料求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
275求調(diào)劑
已經(jīng)有10人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 氣液兩相流問題求助啊
已經(jīng)有1人回復(fù)
- 關(guān)于組分傳輸傳質(zhì)迭代出現(xiàn)致命錯誤，求指點如何調(diào)試udf
已經(jīng)有2人回復(fù)
- 糧倉模擬
已經(jīng)有1人回復(fù)
- 關(guān)于C_YI_G（c,t，i）和F_AREA(A,f,t)的使用
已經(jīng)有1人回復(fù)
- 關(guān)于udf中多相流的指針設(shè)置問題
已經(jīng)有0人回復(fù)
- 氣相在液相中的擴(kuò)散
已經(jīng)有0人回復(fù)
- 關(guān)于UDM的使用順序問題
已經(jīng)有9人回復(fù)
- 關(guān)于燃燒反應(yīng)放出的熱量的計算
已經(jīng)有0人回復(fù)
- 關(guān)于編寫兩相流UDF中使用UDM的一些疑問——請大家多交流
已經(jīng)有2人回復(fù)
- fluent UDF的編寫 一個發(fā)熱球體的熱量隨溫度變化，與空氣對流散熱，請問怎么編寫啊
已經(jīng)有0人回復(fù)
mopsite
木蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 8901.7
- 散金: 424
- 紅花: 2
- 帖子: 2588
- 在線: 615小時
- 蟲號: 172296
- 注冊: 2006-01-19
- 專業(yè): 傳熱傳質(zhì)學(xué)
6樓2015-09-24 13:44:54
唐小虎
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 92.5
- 帖子: 7
- 在線: 5.4小時
- 蟲號: 3393590
- 注冊: 2014-09-02
- 性別: GG
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
9樓2015-10-29 14:25:20
BeObjective
銀蟲 (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1475.1
- 紅花: 6
- 帖子: 733
- 在線: 82.5小時
- 蟲號: 3311402
- 注冊: 2014-07-07
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
2樓2015-09-22 12:04:47
BeObjective
銀蟲 (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1475.1
- 紅花: 6
- 帖子: 733
- 在線: 82.5小時
- 蟲號: 3311402
- 注冊: 2014-07-07
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
3樓2015-09-22 12:07:32
830706ljx
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 273.6
- 散金: 20
- 帖子: 39
- 在線: 9.8小時
- 蟲號: 4092655
- 注冊: 2015-09-22
- 專業(yè): 流體力學(xué)
4樓2015-09-22 15:23:31
BeObjective
銀蟲 (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1475.1
- 紅花: 6
- 帖子: 733
- 在線: 82.5小時
- 蟲號: 3311402
- 注冊: 2014-07-07
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
5樓2015-09-23 10:45:49
mopsite
木蟲 (著名寫手)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 8901.7
- 散金: 424
- 紅花: 2
- 帖子: 2588
- 在線: 615小時
- 蟲號: 172296
- 注冊: 2006-01-19
- 專業(yè): 傳熱傳質(zhì)學(xué)
7樓2015-10-10 12:59:26
BeObjective
銀蟲 (正式寫手)
- 應(yīng)助: 7 (幼兒園)
- 金幣: 1475.1
- 紅花: 6
- 帖子: 733
- 在線: 82.5小時
- 蟲號: 3311402
- 注冊: 2014-07-07
- 性別: GG
- 專業(yè): 傳熱傳質(zhì)學(xué)
|
http://www.gaoyang168.com/bbs/viewthread.php?tid=9465108這是我另外發(fā)的一個帖子，這幾天也在琢磨這個事，自己想法是問題可能出在兩個地方，一個是梯度的處理，另外一個是邊界源相位置的處理，我整個流體域分成了兩個部分，靠近壁面那層網(wǎng)格被我處理成另外一個流體域，這個zone就成了源相的處理區(qū)域，今天在論文看到有人用了其他的辦法。謝謝 有興趣的話，你可以做一做 |
8樓2015-10-11 22:31:06
10樓2015-12-28 17:07:44
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 BeObjective 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +9 | 路癡小琪 2026-03-27 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 295材料工程專碩求調(diào)劑 +5 | 1428151015 2026-03-27 | 5/250 |
|
|
|
[考研] 322求調(diào)劑 +4 | 宋明欣 2026-03-27 | 4/200 |
|
|
|
[考研] 086000調(diào)劑 +3 | 7901117076 2026-03-26 | 3/150 |
|
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
|
[考研] 求調(diào)劑推薦 材料 304 +14 | 荷包蛋hyj 2026-03-26 | 14/700 |
|
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢 2026-03-23 | 16/800 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 4/200 |
|
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 材料科學(xué)與工程 317求調(diào)劑 +4 | JKSOIID 2026-03-26 | 4/200 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。。� +4 | Acaciad 2026-03-25 | 4/200 |
|
|
|
[考研] 各位老師您好：本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +8 | 皮皮皮123456 2026-03-22 | 8/400 |
|
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|