24小時熱門版塊排行榜

返回列表

加油鴨2333

新蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 121.9
散金: 5
帖子: 72
在線: 10.4小時
蟲號: 14033989
注冊: 2019-02-13
專業(yè): 宇宙學(xué)

[交流] 在壁面添加質(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客戶端

回復(fù)此樓

» 猜你喜歡

歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組已經(jīng)有24人回復(fù)
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生已經(jīng)有10人回復(fù)
招調(diào)劑已經(jīng)有11人回復(fù)
0856材料與化工求調(diào)劑！已經(jīng)有11人回復(fù)
材料與化工304求調(diào)劑已經(jīng)有10人回復(fù)
A區(qū)一本交叉課題組，低分調(diào)劑，招收機械電子信息通信等交叉方向已經(jīng)有31人回復(fù)
中科大材料299求調(diào)劑已經(jīng)有7人回復(fù)
301求調(diào)劑已經(jīng)有3人回復(fù)
材料考研339求調(diào)劑已經(jīng)有3人回復(fù)
材料專碩290求調(diào)劑已經(jīng)有5人回復(fù)

1樓 2019-05-31 23:56:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

landykun

鐵蟲 (正式寫手)

應(yīng)助: 9 (幼兒園)
金幣: 2672.6
散金: 55
紅花: 7
帖子: 643
在線: 89.9小時
蟲號: 1772570
注冊: 2012-04-23
性別: GG
專業(yè): 流體力學(xué)

★
小木蟲: 金幣+0.5, 給個紅包，謝謝回帖

首先不知道你憋了多久寫的，不過寫的確實不好，建議好好看看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
}

贊一下

回復(fù)此樓