| 5 | 1/1 | 返回列表 |
| 查看: 1582 | 回復(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í)在是憋不出來(lái)了,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客戶端 |
新蟲 (小有名氣)
送紅花一朵 |
謝謝給予回復(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.我看到別人說(shuō)加壁面源項(xiàng)需要先把臨近壁面的網(wǎng)格提出來(lái),我加的不是整體的源項(xiàng),是一個(gè)面的,所以加的那些指針,請(qǐng)問如果不加直接在壁面上的udf導(dǎo)入,不知道宏能不能找到壁面臨近的cell呢? 發(fā)自小木蟲Android客戶端 |
鐵蟲 (正式寫手)
|
首先 不知道你憋了多久寫的,不過(guò)寫的確實(shí)不好 ,建議好好看看fluent udf幫助里面的source部分,很簡(jiǎn)單的。 1 DEFINE_SOURCE (name, c, t, dS, eqn) 這個(gè)宏,你怎么在這個(gè)括號(hào)里面隨便繼續(xù)加變量?這個(gè)變量是不能變的。 2 包括cell 以及thread,通過(guò)這個(gè)宏fluent已經(jīng)自動(dòng)給你傳進(jìn)來(lái)了,你還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)來(lái)了,你不用管了。 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 } |
鐵蟲 (正式寫手)
|
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í)在不行 那你就采用判斷方式吧。 |
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 家人們 調(diào)劑不迷路 看這里 +8 | likeihood 2026-03-09 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑 +4 | kk扛 2026-03-07 | 4/200 |
|
|
[考博] 2026博士申請(qǐng) +5 | 起泡酒 2026-03-08 | 5/250 |
|
|
[考研] 材料工程085601調(diào)劑求老師收留 +5 | 強(qiáng)木木木 2026-03-07 | 5/250 |
|
|
[考研] 0703求調(diào)劑 +10 | c297914 2026-03-06 | 11/550 |
|
|
[考研] 307求調(diào)劑 +3 | 辛仁豆腐 2026-03-08 | 5/250 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +9 | 1一11 2026-03-03 | 9/450 |
|
|
[考研] 一志愿武理314求調(diào)劑 +4 | ( ̄~ ̄;) 2026-03-08 | 5/250 |
|
|
[考研] 269求調(diào)劑 +3 | 朔朔話 2026-03-08 | 4/200 |
|
|
[考研] 求調(diào)劑,不管什么專業(yè),我是可塑造的人才一枚,希望遇到知己老師撈撈我 +4 | 13102137290 2026-03-06 | 5/250 |
|
|
[考研] 085701環(huán)境工程295求調(diào)劑,一志愿A區(qū)211 +7 | 等春來(lái), 2026-03-06 | 8/400 |
|
|
[考研] 085600材料調(diào)劑 總分330 +6 | 池池丶 2026-03-03 | 6/300 |
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|
|
[考研] 材料與化工,291,求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
[考研] 293求調(diào)劑 +3 | 是樂渝哇 2026-03-04 | 3/150 |
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|