 |
 |
 |
| 24小時(shí)熱門版塊排行榜 |
- >論壇更新日志 (1708)
- >考研 (187)
- >導(dǎo)師招生 (129)
- >蟲友互識(shí) (77)
- >考博 (27)
- >文獻(xiàn)求助 (26)
- >論文投稿 (24)
- >休閑灌水 (18)
- >博后之家 (17)
- >碩博家園 (16)
- >找工作 (15)
- >文獻(xiàn)求助完結(jié)子版 (14)
- >基金申請(qǐng) (12)
- >論文道賀祈福 (8)
- >教師之家 (8)
- >職場(chǎng)人生 (5)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 8704 | 回復(fù): 18 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
[交流]
【求助】求UDF高手幫忙看看我的問題
|
|||||
|
如圖： (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è)溶解過程瞬間即可完成），于是水區(qū)域的上部（區(qū)域3）有了甲醛，而水區(qū)域4還沒有甲醛，于是由于濃度差甲醛會(huì)從水區(qū)域的上部（區(qū)域3）往下部（區(qū)域4）擴(kuò)散。 同時(shí)空氣區(qū)域的下部（區(qū)域2）由于甲醛溶解到水中而濃度降低，也引起空氣區(qū)域中甲醛從上部（區(qū)域1）往下部（區(qū)域2）擴(kuò)散。 (3). 甲醛的溶解過程由平衡分配關(guān)系決定。假設(shè)甲醛在空氣和水中的平衡分配系數(shù)為0.1（即，上述溶解與擴(kuò)散過程經(jīng)過足夠長(zhǎng)時(shí)間后總會(huì)穩(wěn)定，設(shè)穩(wěn)定后 [甲醛在空氣中的濃度a0)：(甲醛在水中的濃度b0)=0.1] 這是一個(gè)物理定律。穩(wěn)定后的濃度比是定值，與初始狀態(tài)無關(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)定后滿足平衡分配系數(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è)平衡。于是新的平衡分配又開始了。過通過循環(huán)完成。 If (b/a<0.1) 溶解平衡分配就會(huì)發(fā)生。 比如最開始甲醛在氣相的濃度為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è)過程。 下面為根據(jù) HELP 中UDF寫的程序，但編譯有錯(cuò)。也不知道能不能正確描述我的問題。望各個(gè)高手幫忙看看。萬分感謝。 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ù)此樓» 收錄本帖的淘帖專輯推薦
 多相流UDF | fluent氣液兩相udf |
» 猜你喜歡
284求調(diào)劑
已經(jīng)有3人回復(fù)
-
085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料專業(yè)申博
已經(jīng)有4人回復(fù)
-
311求調(diào)劑一志愿合肥工業(yè)大學(xué)
已經(jīng)有5人回復(fù)
-
070300化學(xué)279求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料與化工（0856）304求B區(qū)調(diào)劑
已經(jīng)有9人回復(fù)
-
085601 329分調(diào)劑
已經(jīng)有3人回復(fù)
-
本2一志愿C9-333分，材料科學(xué)與工程，求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有5人回復(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é)果分析提出的問題。。。求專家?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ò)，求高手啊�。�！難道沒人遇到過嗎？？？
已經(jīng)有11人回復(fù)
- 高手 幫忙看看這兩數(shù)據(jù)
已經(jīng)有5人回復(fù)
- 關(guān)于udf編譯問題
已經(jīng)有6人回復(fù)
- 導(dǎo)入U(xiǎn)DF時(shí)出這種錯(cuò)誤��！
已經(jīng)有4人回復(fù)
- 求助udf的輸入輸出
已經(jīng)有16人回復(fù)
- 求助！�。∥蚁胍趂luent中針對(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模型的編譯問題
已經(jīng)有15人回復(fù)
- udf中的平方怎么寫
已經(jīng)有5人回復(fù)
- 介個(gè)UDF有啥問題
已經(jīng)有10人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
湖北大學(xué)陳苗苗教授食品安全研究團(tuán)隊(duì)誠招博士生
+1/178
-
中科院三區(qū)，審稿速度快，期刊Journal of Groundwater Science and Engineering
+1/84
-
留學(xué)--博士招生
+1/50
-
天津理工大學(xué)功能晶體研究院（晶體材料全國(guó)重點(diǎn)實(shí)驗(yàn)室）杰青團(tuán)隊(duì)招收2026年碩士研究生
+1/45
-
福建農(nóng)林大學(xué)材料院招收調(diào)劑碩士生
+1/37
-
上海大學(xué) “生物有機(jī)電子材料及器件”團(tuán)隊(duì)碩士、博士研究生招聘
+1/34
-
西京學(xué)院土木水利 2026 級(jí)研究生招生
+2/32
-
西南科技大學(xué)李永升老師接收材料和化工類研究生調(diào)劑
+1/18
-
上海中醫(yī)藥大學(xué)創(chuàng)新中藥研究院 招收審核制博士生一名
+1/16
-
西華大學(xué)材料學(xué)院表面科學(xué)與工程技術(shù)科研團(tuán)隊(duì)2026年招收研究生
+1/11
-
省雙一流重點(diǎn)學(xué)科，國(guó)家一流專業(yè)，輕工技術(shù)與工程，生物質(zhì)能源與材料
+1/10
-
南京醫(yī)科大學(xué)-譚磊課題組-誠招博士后
+1/8
-
歡迎來英國(guó)拉夫堡大學(xué)交流-結(jié)晶，蛋白質(zhì)，AI，水處理等
+1/7
-
食品調(diào)劑-總分291-有中核論文導(dǎo)師一作、本人二作-一志愿南昌大學(xué)學(xué)碩-介紹詳細(xì)
+1/7
-
華東師范大學(xué)-分析化學(xué)-招收1名博士研究生 (2026年，第二批)
+1/7
-
廣東省五邑大學(xué)MOFs催化功能材料課題組招收化學(xué)/材料背景調(diào)劑生
+1/6
-
武漢工程大學(xué)化學(xué)與環(huán)境工程學(xué)院招收調(diào)劑生
+1/4
-
天津醫(yī)科大學(xué)基礎(chǔ)醫(yī)學(xué)院張恒課題組博士后招聘
+1/3
-
南京林業(yè)大學(xué)化學(xué)工程學(xué)院柏惺峰課題組招收碩士研究生1名
+1/3
-
武漢理工大學(xué)船海與能源動(dòng)力工程學(xué)院劉鵬飛院士團(tuán)隊(duì)博士后招聘啟事
+1/1
5樓2011-01-04 16:06:50
★ ★ ★
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
|
問題是有一些： 1. domain指針沒有賦值，要加上這么一句：dimain = Get_Domain(1); 2. DEFINE_MASS_TRANSFER宏是用于每個(gè)cell的，好像不是對(duì)face賦值的，更不能用begin_f_loop了。 3. 我對(duì)這個(gè)UDF的理解還不深，或許用begin_f_loop是沒有問題的，但是在循環(huán)里面用return (c);就有問題了，你的第一次循環(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)一下單位的問題，返回值單位要求是kg/m3-s，我對(duì)你的公式不熟，不知道對(duì)不對(duì)。 |
2樓2011-01-04 10:00:54
3樓2011-01-04 10:06:00
4樓2011-01-04 10:07:28
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 22408 359分調(diào)劑 +4 | Qshers 2026-03-27 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 281求調(diào)劑 +5 | 亞克西good 2026-03-26 | 7/350 |
|
|
[考研] 085600 295分求調(diào)劑 +7 | W55j 2026-03-30 | 9/450 |
|
|
|
[考研] 一志愿華中師范化學(xué)332分求調(diào)劑 +3 | Lyy930824@ 2026-03-29 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +5 | Faune 2026-03-30 | 5/250 |
|
|
|
[考研] 085701求調(diào)劑初試286分 +5 | secret0328 2026-03-28 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
|
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
|
[考研] 調(diào)劑求院校招收 +6 | 鶴鯨鴿 2026-03-28 | 6/300 |
|
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5 | 0857求調(diào)劑. 2026-03-24 | 6/300 |
|
|
|
[考研] 11408軟件工程求調(diào)劑 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 墨墨漠 2026-03-25 | 7/350 |
|
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
|
[考研] 機(jī)械學(xué)碩310分，數(shù)一英一，一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|