| 5 | 1/1 | 返回列表 |
| 查看: 1581 | 回復(fù): 4 | |||
加油鴨2333新蟲 (小有名氣)
|
[交流]
在壁面添加質(zhì)量源項 已有1人參與
|
|
請問如何在壁面添加質(zhì)量源項呢?我想加一個-KC,就是一級反應(yīng)動力學(xué)里的那個,使到壁面的污染物以-kc減少,我看了很久,知道需要提取臨近壁面的cell,可是,如何表示方程式和用壁面的污染物質(zhì)量分?jǐn)?shù)呢?我實在是憋不出來了,UDF最多寫成這個樣子,求助蟲友們幫我完善UDF,謝謝,謝謝。!其中一級反應(yī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客戶端 |
鐵蟲 (正式寫手)
|
首先 不知道你憋了多久寫的,不過寫的確實不好 ,建議好好看看fluent udf幫助里面的source部分,很簡單的。 1 DEFINE_SOURCE (name, c, t, dS, eqn) 這個宏,你怎么在這個括號里面隨便繼續(xù)加變量?這個變量是不能變的。 2 包括cell 以及thread,通過這個宏fluent已經(jī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),都是多余的。你只需要指定每個cell里面的消耗速率值即可,置于是在哪個cell 哪個thread,這個宏已經(jīng)給你傳進(jìn)來了,你不用管了。 3你的消耗速率公式是什么? 4 ds的數(shù)值 我寫udf的時候 基本不反饋; 5 消耗速率計算計算公式里面還有時間嗎,不太清楚。 5 我覺的應(yīng)該這么寫, #include "udf.h" DEFINE_SOURCE (kc_source, c, t, dS, eqn) { real source; real yk,k yk=每個cell的組分質(zhì)量百分比宏函數(shù),你找找udf幫助,我忘了; k=0.274; source=消耗速率計算公式; return source } |
新蟲 (小有名氣)
送紅花一朵 |
2樓: Originally posted by landykun at 2019-06-02 11:31:36 首先 不知道你憋了多久寫的,不過寫的確實不好 ,建議好好看看fluent udf幫助里面的source部分,很簡單的。 1 DEFINE_SOURCE (name, c, t, dS, eqn) 這個宏,你怎么在這個括號里面隨便繼續(xù)加變量?這個變量是不能變 ... 謝謝給予回復(fù),謝謝! 1.我的反應(yīng)公式就是一級反應(yīng)動力學(xué)公式,dc/dt=-Kc,所以把它解析解之后,就帶了時間,請問我這個公式需要那么變嗎?因為看源項案例,都把他們解析解了 2.式中是濃度,但是入口條件只有寫質(zhì)量分?jǐn)?shù),我直接當(dāng)他們相等了?梢灾苯佑眠@個公式嗎?還是需要換算一下? 2.我看到別人說加壁面源項需要先把臨近壁面的網(wǎng)格提出來,我加的不是整體的源項,是一個面的,所以加的那些指針,請問如果不加直接在壁面上的udf導(dǎo)入,不知道宏能不能找到壁面臨近的cell呢? 發(fā)自小木蟲Android客戶端 |
鐵蟲 (正式寫手)
|
3樓: Originally posted by 加油鴨2333 at 2019-06-03 08:10:35 謝謝給予回復(fù),謝謝! 1.我的反應(yīng)公式就是一級反應(yīng)動力學(xué)公式,dc/dt=-Kc,所以把它解析解之后,就帶了時間,請問我這個公式需要那么變嗎?因為看源項案例,都把他們解析解了 2.式中是濃度,但是入口條件只有寫質(zhì)量 ... 1 源項都是單位時間內(nèi)的源項,是個速率,所以 我覺得不用加,你再研究一下; 2 根據(jù)公式 公式怎么樣就怎么樣,肯定要轉(zhuǎn)換啊,質(zhì)量分?jǐn)?shù)變成濃度; 3 你畫網(wǎng)格的時候,壁面一層單獨畫成一個zone,如果是用icem畫網(wǎng)格 就是單獨一個block 如果是gambit或者其他的 我還真不完全清楚,然后在改zone里面 加載source;如果你想在整個計算區(qū)域都加載source 然后判斷臨近邊界的區(qū)域 我覺得太麻煩了,可能還浪費計算時間。所以,如果網(wǎng)格好畫的話,建議在壁面第一層,單獨畫個zone。如果實在不行 那你就采用判斷方式吧。 |
新蟲 (小有名氣)
送紅花一朵 |
4樓: Originally posted by landykun at 2019-06-03 12:38:21 1 源項都是單位時間內(nèi)的源項,是個速率,所以 我覺得不用加,你再研究一下; 2 根據(jù)公式 公式怎么樣就怎么樣,肯定要轉(zhuǎn)換啊,質(zhì)量分?jǐn)?shù)變成濃度; 3 你畫網(wǎng)格的時候,壁面一層單獨畫成一個zone,如果是用icem畫網(wǎng) ... 非常感謝。! 發(fā)自小木蟲Android客戶端 |
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請] 進(jìn)入個人成果庫好難,一下午都沒進(jìn)去 +3 | mi_dilee 2026-03-05 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +4 | 劉合華 2026-03-05 | 4/200 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-06 | 3/150 |
|
|
[考研] 274求調(diào)劑 +9 | 一個學(xué)習(xí)者 2026-03-04 | 9/450 |
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
[考研] 347求調(diào)劑 +6 | 啊歐歐歐 2026-03-03 | 8/400 |
|
|
[考研] 0703化學(xué) 學(xué)碩 理工科均可 不區(qū)分研究方向 總分279求調(diào)劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
[考研] 085600材料與化工 298 調(diào)劑 +9 | 小西笑嘻嘻 2026-03-03 | 9/450 |
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
[考研] 一志愿鄭州大學(xué),學(xué)碩,物理化學(xué), 333求調(diào)劑 +5 | 李魔女斗篷 2026-03-04 | 5/250 |
|
|
[考研] 材料化工調(diào)劑 +15 | 今夏不夏 2026-03-01 | 18/900 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 化工270求調(diào)劑 +10 | 什么名字qwq 2026-03-02 | 10/500 |
|
|
[考研] 課題組接收材料類調(diào)劑研究生 +6 | gaoxiaoniuma 2026-02-28 | 9/450 |
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
[考研] 0856求調(diào)劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 一志愿中石油(華東)本科齊魯工業(yè)大學(xué) +3 | 石能偉 2026-03-02 | 3/150 |
|
|
[考研] 311求調(diào)劑 +6 | 亭亭亭01 2026-03-01 | 6/300 |
|