| 5 | 1/1 | 返回列表 |
| 查看: 1580 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
加油鴨2333新蟲 (小有名氣)
|
[交流]
在壁面添加質(zhì)量源項(xiàng) 已有1人參與
|
||
|
請(qǐng)問如何在壁面添加質(zhì)量源項(xiàng)呢?我想加一個(gè)-KC,就是一級(jí)反應(yīng)動(dòng)力學(xué)里的那個(gè),使到壁面的污染物以-kc減少,我看了很久,知道需要提取臨近壁面的cell,可是,如何表示方程式和用壁面的污染物質(zhì)量分?jǐn)?shù)呢?我實(shí)在是憋不出來了,UDF最多寫成這個(gè)樣子,求助蟲友們幫我完善UDF,謝謝,謝謝。。∑渲幸患(jí)反應(yīng)動(dòng)力學(xué)常數(shù)是0.274,壁面的ID是13,我用yk表示污染物的質(zhì)量分?jǐn)?shù) #include "udf.h" #define K 0.274 DEFINE_SOURCE(face_source, cell, face, thread, yk, dS, eqn) { real source; fluid_ID = 13; Domain *domain; ?Thread *t = Lookup_Thread(domain,fluid_ID); ?Thread *t0 = THREAD_T0(t); cell_t c0 = F_C0(f,tf); tc0 = THREAD_T0(tf); source=yk*exp(-Kt) dS[eqn]=-K*yk return source } 發(fā)自小木蟲Android客戶端 |
鐵蟲 (正式寫手)
|
1 源項(xiàng)都是單位時(shí)間內(nèi)的源項(xiàng),是個(gè)速率,所以 我覺得不用加,你再研究一下; 2 根據(jù)公式 公式怎么樣就怎么樣,肯定要轉(zhuǎn)換啊,質(zhì)量分?jǐn)?shù)變成濃度; 3 你畫網(wǎng)格的時(shí)候,壁面一層單獨(dú)畫成一個(gè)zone,如果是用icem畫網(wǎng)格 就是單獨(dú)一個(gè)block 如果是gambit或者其他的 我還真不完全清楚,然后在改zone里面 加載source;如果你想在整個(gè)計(jì)算區(qū)域都加載source 然后判斷臨近邊界的區(qū)域 我覺得太麻煩了,可能還浪費(fèi)計(jì)算時(shí)間。所以,如果網(wǎng)格好畫的話,建議在壁面第一層,單獨(dú)畫個(gè)zone。如果實(shí)在不行 那你就采用判斷方式吧。 |
鐵蟲 (正式寫手)
|
首先 不知道你憋了多久寫的,不過寫的確實(shí)不好 ,建議好好看看fluent udf幫助里面的source部分,很簡(jiǎn)單的。 1 DEFINE_SOURCE (name, c, t, dS, eqn) 這個(gè)宏,你怎么在這個(gè)括號(hào)里面隨便繼續(xù)加變量?這個(gè)變量是不能變的。 2 包括cell 以及thread,通過這個(gè)宏fluent已經(jīng)自動(dòng)給你傳進(jìn)來了,你還Thread *t = Lookup_Thread(domain,fluid_ID); Thread *t0 = THREAD_T0(t); cell_t c0 = F_C0(f,tf); tc0 = THREAD_T0(tf),都是多余的。你只需要指定每個(gè)cell里面的消耗速率值即可,置于是在哪個(gè)cell 哪個(gè)thread,這個(gè)宏已經(jīng)給你傳進(jìn)來了,你不用管了。 3你的消耗速率公式是什么? 4 ds的數(shù)值 我寫udf的時(shí)候 基本不反饋; 5 消耗速率計(jì)算計(jì)算公式里面還有時(shí)間嗎,不太清楚。 5 我覺的應(yīng)該這么寫, #include "udf.h" DEFINE_SOURCE (kc_source, c, t, dS, eqn) { real source; real yk,k yk=每個(gè)cell的組分質(zhì)量百分比宏函數(shù),你找找udf幫助,我忘了; k=0.274; source=消耗速率計(jì)算公式; return source } |
新蟲 (小有名氣)
送紅花一朵 |
謝謝給予回復(fù),謝謝! 1.我的反應(yīng)公式就是一級(jí)反應(yīng)動(dòng)力學(xué)公式,dc/dt=-Kc,所以把它解析解之后,就帶了時(shí)間,請(qǐng)問我這個(gè)公式需要那么變嗎?因?yàn)榭丛错?xiàng)案例,都把他們解析解了 2.式中是濃度,但是入口條件只有寫質(zhì)量分?jǐn)?shù),我直接當(dāng)他們相等了?梢灾苯佑眠@個(gè)公式嗎?還是需要換算一下? 2.我看到別人說加壁面源項(xiàng)需要先把臨近壁面的網(wǎng)格提出來,我加的不是整體的源項(xiàng),是一個(gè)面的,所以加的那些指針,請(qǐng)問如果不加直接在壁面上的udf導(dǎo)入,不知道宏能不能找到壁面臨近的cell呢? 發(fā)自小木蟲Android客戶端 |
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +7 | Faune 2026-03-06 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +3 | 鋼鐵大炮 2026-03-04 | 3/150 |
|
|
[考研] 求調(diào)劑,學(xué)校研究所都可以,材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
[考研] 290求調(diào)劑 +5 | Bananaiy 2026-03-04 | 5/250 |
|
|
[考研] 267化工調(diào)劑求助 +5 | 聰少OZ 2026-03-04 | 5/250 |
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
[論文投稿]
EST拒稿重投
5+3
|
15102603076 2026-03-02 | 3/150 |
|
|
[考研] 0703 總分319求調(diào)劑 +4 | Xinyuu 2026-03-02 | 4/200 |
|
|
[考研] 298求調(diào)劑一志愿中海洋 +3 | lour. 2026-03-03 | 3/150 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|
|
[考研] 理學(xué),工學(xué),農(nóng)學(xué)調(diào)劑,少走彎路,這里歡迎您! +8 | likeihood 2026-03-02 | 11/550 |
|
|
[碩博家園] 2025屆雙非化工碩士畢業(yè),申博 +4 | 更多的是 2026-02-27 | 5/250 |
|
|
[考研] 288求調(diào)劑 +3 | 少71.8 2026-03-02 | 5/250 |
|
|
[考研] 求調(diào)劑 +7 | repeatt?t 2026-02-28 | 7/350 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 264求調(diào)劑 +4 | 巴拉巴拉根556 2026-02-28 | 4/200 |
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|
|
[考研] 307求調(diào)劑 +4 | 73372112 2026-02-28 | 6/300 |
|