| 5 | 1/1 | 返回列表 |
| 查看: 1584 | 回復(fù): 4 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
加油鴨2333新蟲 (小有名氣)
|
[交流]
在壁面添加質(zhì)量源項(xiàng) 已有1人參與
|
||
|
請問如何在壁面添加質(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部分,很簡單的。 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í)間,請問我這個(gè)公式需要那么變嗎?因?yàn)榭丛错?xiàng)案例,都把他們解析解了 2.式中是濃度,但是入口條件只有寫質(zhì)量分?jǐn)?shù),我直接當(dāng)他們相等了?梢灾苯佑眠@個(gè)公式嗎?還是需要換算一下? 2.我看到別人說加壁面源項(xiàng)需要先把臨近壁面的網(wǎng)格提出來,我加的不是整體的源項(xiàng),是一個(gè)面的,所以加的那些指針,請問如果不加直接在壁面上的udf導(dǎo)入,不知道宏能不能找到壁面臨近的cell呢? 發(fā)自小木蟲Android客戶端 |
新蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
likeihood 2026-03-06 | 16/800 |
|
|---|---|---|---|---|
|
[考研] 一志愿清華深研院材料專碩294分,專業(yè)課111分,本科中南大學(xué)材料,有六級(jí),有工作經(jīng)驗(yàn) +4 | H14528 2026-03-04 | 4/200 |
|
|
[考研] 一志愿北京化工大學(xué)材料與化工296分求調(diào)劑 +9 | 稻妻小編 2026-03-09 | 9/450 |
|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
10+4
|
穿只靴子 2026-03-07 | 15/750 |
|
|
[考研] 085600材料與化工,一志愿廣州985,求調(diào)劑 +15 | qqyyaill 2026-03-05 | 15/750 |
|
|
[考研] 264求調(diào)劑 +5 | thext 2026-03-03 | 5/250 |
|
|
[考研] 320材料與化工,求調(diào)劑 +11 | 鶴遨予卿 2026-03-04 | 13/650 |
|
|
[考研] 考研調(diào)劑,一志愿山東大學(xué)材料與化工,328分,政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
[考研] 材料調(diào)劑 +4 | xxxcm 2026-03-08 | 7/350 |
|
|
[考研] 一志愿東北大學(xué)080500,英一數(shù)二343分求調(diào)劑 +7 | Seele_v5 2026-03-04 | 7/350 |
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
[考研] 考研求調(diào)劑,工科,總分304 +4 | 矮子不想長高了 2026-03-08 | 7/350 |
|
|
[考研] 083000,總分284,求調(diào)劑 +5 | 徐yr 2026-03-04 | 5/250 |
|
|
[考研] 化學(xué)工程求調(diào)劑 +12 | 化工人999 2026-03-04 | 12/600 |
|
|
[考研] 070300化學(xué)求調(diào)劑292分 +3 | 打烊eee 2026-03-07 | 3/150 |
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
[考研] 材料學(xué)碩080500復(fù)試調(diào)劑294 +3 | 四葉zjz 2026-03-04 | 3/150 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 0703化學(xué)求調(diào)劑 +3 | 將未有期 2026-03-04 | 3/150 |
|