 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 19 | 1/1 | 返回列表 |
| 查看: 8692 | 回復(fù): 18 | |||||
kylafree至尊木蟲(chóng) (知名作家)
|
[交流]
【求助】求UDF高手幫忙看看我的問(wèn)題
|
||||
|
如圖： (1).初始時(shí)刻，矩形區(qū)域內(nèi)上半部分（區(qū)域1+2）為質(zhì)量分?jǐn)?shù)為5%的甲醛+95%的空氣的混合氣體。下半部分（區(qū)域3+4）為0%的甲醛+100%的水。 (2). 空氣中的甲醛會(huì)在氣液交界面處發(fā)生溶解，即區(qū)域2中的甲醛因溶解而進(jìn)入?yún)^(qū)域3中（設(shè)溶解過(guò)程瞬間即可完成），于是水區(qū)域的上部（區(qū)域3）有了甲醛，而水區(qū)域4還沒(méi)有甲醛，于是由于濃度差甲醛會(huì)從水區(qū)域的上部（區(qū)域3）往下部（區(qū)域4）擴(kuò)散。 同時(shí)空氣區(qū)域的下部（區(qū)域2）由于甲醛溶解到水中而濃度降低，也引起空氣區(qū)域中甲醛從上部（區(qū)域1）往下部（區(qū)域2）擴(kuò)散。 (3). 甲醛的溶解過(guò)程由平衡分配關(guān)系決定。假設(shè)甲醛在空氣和水中的平衡分配系數(shù)為0.1（即，上述溶解與擴(kuò)散過(guò)程經(jīng)過(guò)足夠長(zhǎng)時(shí)間后總會(huì)穩(wěn)定，設(shè)穩(wěn)定后 [甲醛在空氣中的濃度a0)：(甲醛在水中的濃度b0)=0.1] 這是一個(gè)物理定律。穩(wěn)定后的濃度比是定值，與初始狀態(tài)無(wú)關(guān)。 (4). 選用VOF兩相流模型（主相設(shè)為甲醛（占5%）+空氣（占95%）的混合氣，第二相設(shè)為甲醛（占0%）和水（占100%）的混和物） +species transport模型。 氣液交界面為邊界條件設(shè)為wall，其它邊界也為wall 兩相交互作用Phase interacion中設(shè)置一個(gè)mass transfer mechnism （5）.編程想法： 找到氣液交界面，然后提取出交界面上方單元（區(qū)域2）甲醛的濃度，設(shè)其為a； 再提取出交界面下方單元（區(qū)域3）甲醛的濃度，設(shè)為b； 再根據(jù)甲醛在氣相和液相中的平衡分配系數(shù)0.1 由 a+b=a0+b0 （分配前后質(zhì)量守恒） a0/b0=0.1 （分配穩(wěn)定后滿(mǎn)足平衡分配系數(shù)關(guān)系） 解上面方程組得a0=（0.1/(0.1+1) ）* (a+b) ，b0=（1/(0.1+1)） * (a+b) 并把a(bǔ)更新為a0，b更新為b0。 即a=（0.1/(0.1+1)） * (a+b) b=（1/(0.1+1) ）* (a+b) (6).溶解平衡分配關(guān)系完成后，根據(jù)上面物理原理第（2）條，由于擴(kuò)散原因會(huì)打破這個(gè)平衡。于是新的平衡分配又開(kāi)始了。過(guò)通過(guò)循環(huán)完成。 If (b/a<0.1) 溶解平衡分配就會(huì)發(fā)生。 比如最開(kāi)始甲醛在氣相的濃度為a=5%，甲醛在水中的濃度為b=0， 顯然b/a<0.1，所以按平衡分配關(guān)系把a(bǔ),b更新為 a=0.1/(0.1+1) * 5% b=1/(0.1+1) *5% 這樣氣液兩相中便出現(xiàn)濃度差了，兩相中各自又會(huì)發(fā)生擴(kuò)散，擴(kuò)散又會(huì)打破平衡分配關(guān)系...... 就是這個(gè)過(guò)程。 下面為根據(jù) HELP 中UDF寫(xiě)的程序，但編譯有錯(cuò)。也不知道能不能正確描述我的問(wèn)題。望各個(gè)高手幫忙看看。萬(wàn)分感謝。 DEFINE_MASS_TRANSFER(liq_gas_transfer, cell, mix_thread, from_index, from_species_index, to_index, to_species_index) { real a,b,c; cell_t c0,c1; Thread *tc0,*tc1; Domain *domain=Get_Domain(1); face_t f; int ID=9; /* 交界處wall邊界的ID*/ Thread *zhong_thread=Lookup_Thread(domain,ID); /* 獲取交界面處wall的指針*/ begin_f_loop(f,zhong_thread) /* 循環(huán)于交界面中的所有face*/ { c0=F_CO(f,zhong_thread); /*左邊的cell標(biāo)為c0*/ tc0=THREAD_T0(zhong_thread); /*左邊的cell的指針tc0*/ c1=F_C1(f,zhong_thread); /*右邊的cell標(biāo)為c1*/ tc1=THREAD_T1(zhong_thread); /*右邊的cell的指針tc1*/ a=C_YI(c0,tc0,from_species_index); /*獲取co的質(zhì)量分?jǐn)?shù)*/ b=C_YI(c1,tc1,to_species_index); /*獲取c1的質(zhì)量分?jǐn)?shù)*/ if(a>(10.0/11.0)*(a+b)) c=(10.0/11.0)*(a+b)-a; /*按分配關(guān)系，如果a>分配后的a0，則記錄應(yīng)返回值c=a0-a*/ if(b<(1.0/11.0)*(a+b)) c=(1.0/11.0)*(a+b)-b; /*按分配關(guān)系，如果b<分配后的b0，則記錄應(yīng)返回值c=b0-b*/ } return(c); /*返回c*/ end_f_loop(f,zhong_thread) } [ Last edited by kylafree on 2011-1-5 at 16:37 ] |
回復(fù)此樓» 收錄本帖的淘帖專(zhuān)輯推薦
 多相流UDF | fluent氣液兩相udf |
» 猜你喜歡
085701環(huán)境工程，267求調(diào)劑
已經(jīng)有12人回復(fù)
-
320分，材料與化工專(zhuān)業(yè)，求調(diào)劑
已經(jīng)有7人回復(fù)
-
085701環(huán)境工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
295材料工程專(zhuān)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
265求調(diào)劑
已經(jīng)有8人回復(fù)
-
322求調(diào)劑
已經(jīng)有4人回復(fù)
-
085602 307分 求調(diào)劑
已經(jīng)有6人回復(fù)
-
0856，材料與化工321分求調(diào)劑
已經(jīng)有7人回復(fù)
-
0703化學(xué)求調(diào)劑，各位老師看看我�。�！
已經(jīng)有5人回復(fù)
-
266分，求材料冶金能源化工等調(diào)劑
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求大神，UDF導(dǎo)入之后怎么檢測(cè)UDF中的設(shè)置已經(jīng)成功應(yīng)用？
已經(jīng)有3人回復(fù)
- udf不同壁面熱流求助
已經(jīng)有14人回復(fù)
- 求助 審稿人對(duì)交互作用的結(jié)果分析提出的問(wèn)題。。。求專(zhuān)家?guī)兔纯词鞘裁匆馑迹?/a>
已經(jīng)有4人回復(fù)
- UDF并行編譯求助-附程序段
已經(jīng)有3人回復(fù)
- 求助UDF高手-----DEFINE_SOURCE中速度和面積矢量的點(diǎn)積
已經(jīng)有4人回復(fù)
- {求助，關(guān)于一個(gè)速度入口的udf}
已經(jīng)有12人回復(fù)
- fluent14編譯UDF時(shí)總是出錯(cuò)，求高手�。。�！難道沒(méi)人遇到過(guò)嗎？？？
已經(jīng)有11人回復(fù)
- 高手 幫忙看看這兩數(shù)據(jù)
已經(jīng)有5人回復(fù)
- 關(guān)于udf編譯問(wèn)題
已經(jīng)有6人回復(fù)
- 導(dǎo)入U(xiǎn)DF時(shí)出這種錯(cuò)誤��！
已經(jīng)有4人回復(fù)
- 求助udf的輸入輸出
已經(jīng)有16人回復(fù)
- 求助�。�！我想要在fluent中針對(duì)一個(gè)面進(jìn)行積分，那么UDF應(yīng)該用什么宏函數(shù)？怎么用？
已經(jīng)有5人回復(fù)
- 能量方程源項(xiàng)udf
已經(jīng)有4人回復(fù)
- 曳力函數(shù)的的UDF
已經(jīng)有15人回復(fù)
- 求助：UDF模型的編譯問(wèn)題
已經(jīng)有15人回復(fù)
- udf中的平方怎么寫(xiě)
已經(jīng)有5人回復(fù)
- 介個(gè)UDF有啥問(wèn)題
已經(jīng)有10人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
求助紫外可見(jiàn)分光光度計(jì)
+1/99
-
贛南師范大學(xué)肖興教授課題組招2026級(jí)電子信息碩士 4-6 名，歡迎調(diào)劑
+1/82
-
重慶交通大學(xué)院士團(tuán)隊(duì)招收材料、能源動(dòng)力專(zhuān)業(yè)碩士研究生
+2/50
-
天津科技大學(xué)-（院士團(tuán)隊(duì)）先進(jìn)纖維與紙基功能材料團(tuán)隊(duì) 招收碩士生
+1/46
-
化工畢業(yè)三年轉(zhuǎn)電子信息可行嗎？？各位同行前輩
+1/17
-
【博士招生】南京林業(yè)大學(xué)車(chē)輛工程課題組招收26年9月入學(xué)“申請(qǐng)-考核制”博士研究生
+1/14
-
煙臺(tái)大學(xué)精準(zhǔn)材料高等研究院26年材料方向研究生招生
+1/12
-
港科大（廣州）招收光學(xué)傳感器，集成電路設(shè)計(jì)，醫(yī)學(xué)影像設(shè)備以及可穿戴設(shè)備方向的博士
+1/9
-
西安工程大學(xué)環(huán)境與化學(xué)工程學(xué)院能源環(huán)境催化與先進(jìn)材料課題組招生
+1/9
-
江西科技師范大學(xué) 材料與能源招收調(diào)劑研究生
+1/8
-
中國(guó)科學(xué)院上海硅酸鹽研究所招聘科研助理（納米材料與生物應(yīng)用方向）
+1/7
-
南京林業(yè)大學(xué)化學(xué)工程學(xué)院柏惺峰課題組誠(chéng)聘科研助理、水杉英才E類(lèi)/F類(lèi)青年教師或博后
+1/7
-
哈工大赫曉東院士團(tuán)隊(duì)隋超教授招收材料學(xué)博士1名
+1/7
-
【MP期刊】青年編委招募開(kāi)啟！加入我們，一起推動(dòng)材料科學(xué)前沿
+1/6
-
26博士申請(qǐng)
+1/6
-
寧波東方理工大學(xué)余鵬課題組2026年上海交大聯(lián)培博士生、博士后、科研助理招聘
+1/3
-
中國(guó)科學(xué)技術(shù)大學(xué)納米科學(xué)技術(shù)學(xué)院張振教授/江雷院士誠(chéng)招2026級(jí)博士
+1/3
-
北京石油化工學(xué)院康舉老師團(tuán)招2026年研究生復(fù)試調(diào)劑考生，歡迎聯(lián)系，
+1/3
-
武漢紡織大學(xué)2026年招收碩士研究生【高性能、功能紡織復(fù)合材料方向】
+1/2
-
【博士招生】上海交通大學(xué)機(jī)動(dòng)學(xué)院-招2026年秋“申請(qǐng)-考核制”博士生1名-3D打印方向
+1/1
★ ★ ★
kylafree(金幣+30): 2011-01-04 16:07:18
robert2020(金幣+3):辛苦了！ 2011-01-11 10:48:25
kylafree(金幣+30): 2011-01-04 16:07:18
robert2020(金幣+3):辛苦了！ 2011-01-11 10:48:25
|
問(wèn)題是有一些： 1. domain指針沒(méi)有賦值，要加上這么一句：dimain = Get_Domain(1); 2. DEFINE_MASS_TRANSFER宏是用于每個(gè)cell的，好像不是對(duì)face賦值的，更不能用begin_f_loop了。 3. 我對(duì)這個(gè)UDF的理解還不深，或許用begin_f_loop是沒(méi)有問(wèn)題的，但是在循環(huán)里面用return (c);就有問(wèn)題了，你的第一次循環(huán)就會(huì)返回值并結(jié)束函數(shù)。 4. begin_f_loop(f,interior_thread) 對(duì)應(yīng)的應(yīng)該是end_f_loop(f,interior_thread) 。 5. 你確認(rèn)一下單位的問(wèn)題，返回值單位要求是kg/m3-s，我對(duì)你的公式不熟，不知道對(duì)不對(duì)。 |
2樓2011-01-04 10:00:54
★ ★
xiegangmai(金幣+2):謝謝參與 2011-01-15 08:42:43
xiegangmai(金幣+2):謝謝參與 2011-01-15 08:42:43
|
關(guān)于樓主的這個(gè)問(wèn)題，我說(shuō)說(shuō)我想的編程思路吧： 用DEFINE_MASS_TRANSFER宏，在函數(shù)里，先判斷當(dāng)前的cell是不是靠近相界面的cell；如果不是，返回0；如果是，判斷是液相還是氣相；如果是氣相，那么要找到與之相對(duì)應(yīng)的液相cell，分別獲取二者中需要的參數(shù)，求取質(zhì)量傳遞速率，并返回值；如果是液相，同樣找到對(duì)應(yīng)的氣相cell，求取傳遞速率并返回值。 |
3樓2011-01-04 10:06:00
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
8樓2011-01-04 17:01:07
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
10樓2011-01-04 17:11:01
zzt_1234
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1683.4
- 帖子: 1726
- 在線: 169.2小時(shí)
- 蟲(chóng)號(hào): 809873
13樓2011-01-14 16:35:14
14樓2012-02-19 22:32:28
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
15樓2012-02-21 11:33:12
4樓2011-01-04 10:07:28
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
5樓2011-01-04 16:06:50
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
robert2020(金幣+2):多謝應(yīng)助！ 2011-01-11 10:48:58
robert2020(金幣-2):不好意思，評(píng)帖失誤，歡迎常來(lái)仿真模擬版交流！ 2011-01-11 10:50:43
robert2020(金幣-2):不好意思，評(píng)帖失誤，歡迎常來(lái)仿真模擬版交流！ 2011-01-11 10:50:43
|
這個(gè)UDF，就是先找到氣液兩相的交界面， 然后提取出交界面上方單元甲醛的濃度，即甲醛在氣相中的濃度，設(shè)為a； 再提取出交界面下方單元甲醛的濃度，即甲醛在液相中的濃度，設(shè)為b； 再根據(jù)甲醛在氣相和液相中的平衡分配系數(shù)0.1（即若擴(kuò)散穩(wěn)定后，甲醛在氣相中的濃度：甲醛在液相中的濃度=0.1） 根據(jù) a+b=a0+b0 a0/b0=0.1 把a(bǔ)更新為0.1/(0.1+1) * (a+b) 把b更新為1/(0.1+1) * (a+b) 這實(shí)際上是空氣中的甲醛溶解于水中的過(guò)程，設(shè)它在一瞬間即完成。 只要a/b！=0.1（事實(shí)上它是>=0.1的，如初始時(shí)刻a=5%，b=0）這個(gè)過(guò)程就會(huì)發(fā)生。 比如最開(kāi)始甲醛在氣相的濃度為a=5%，甲醛在水中的濃度為b=0， 顯然a/b>0.1，所以按平衡分配關(guān)系把a(bǔ),b更新為 a=0.1/(0.1+1) * 5% b=1/(0.1+1) *5% 這樣氣液兩相中便出現(xiàn)濃度差了，兩相中各自又會(huì)發(fā)生擴(kuò)散，擴(kuò)散又會(huì)打破平衡分配關(guān)系...... 就是這個(gè)過(guò)程。 [ Last edited by kylafree on 2011-1-4 at 17:12 ] |
6樓2011-01-04 16:49:27
zzt_1234
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1683.4
- 帖子: 1726
- 在線: 169.2小時(shí)
- 蟲(chóng)號(hào): 809873
7樓2011-01-04 16:59:36
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
9樓2011-01-04 17:08:11
kylafree
至尊木蟲(chóng) (知名作家)
- 應(yīng)助: 166 (高中生)
- 金幣: 17487
- 帖子: 6212
- 在線: 1161.7小時(shí)
- 蟲(chóng)號(hào): 1061079
11樓2011-01-04 19:37:51
zzt_1234
木蟲(chóng) (著名寫(xiě)手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1683.4
- 帖子: 1726
- 在線: 169.2小時(shí)
- 蟲(chóng)號(hào): 809873
12樓2011-01-14 16:00:11
★
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
小木蟲(chóng): 金幣+0.5, 給個(gè)紅包，謝謝回帖
|
請(qǐng)問(wèn)LZ，在設(shè)置的mass transfer mechanism時(shí)采用的是什么mechanism？ fluent14的user Guide，26.2.8中有關(guān)Mass Transfer Effects有這么一句：Unidirectional constant rate mass transfer (not available for VOF calculations)，不知道LZ當(dāng)時(shí)采取的是哪一種mechanism呢？cavitation, evaporation-condensation,or boiling？還是自定義的？能夠貼出自定義的項(xiàng)目。 |
16樓2013-05-14 15:49:31
ya634984094
銅蟲(chóng) (小有名氣)
- 應(yīng)助: 39 (小學(xué)生)
- 金幣: 74.7
- 帖子: 259
- 在線: 38.5小時(shí)
- 蟲(chóng)號(hào): 1695458
17樓2014-01-20 19:45:04
18樓2014-08-08 10:37:12
19樓2015-11-07 23:15:03
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 kylafree 的主題更新
| 19 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 340求調(diào)劑 +4 | jhx777 2026-03-27 | 4/200 |
|
|
|
[考研] 07化學(xué)280分求調(diào)劑 +10 | 722865 2026-03-23 | 10/500 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-27 | 3/150 |
|
|
|
[考研] 0856材料專(zhuān)碩353求調(diào)劑 +9 | NIFFFfff 2026-03-20 | 9/450 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +6 | 小小甜面團(tuán) 2026-03-25 | 6/300 |
|
|
|
[考研] 總分322求生物學(xué)/生化與分子/生物信息學(xué)相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 086000生物與醫(yī)藥292求調(diào)劑 +6 | 小小陳小小 2026-03-22 | 9/450 |
|
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 10/500 |
|
|
|
[考研] 303求調(diào)劑 +6 | 藍(lán)山月 2026-03-25 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap－ 2026-03-22 | 5/250 |
|
|
|
[考研] 一志愿吉大化學(xué)322求調(diào)劑 +4 | 17501029541 2026-03-23 | 6/300 |
|
|
|
[考博] 26申博自薦 +3 | whh869393 2026-03-24 | 3/150 |
|
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|