| 4 | 1/1 | 返回列表 |
| 查看: 1075 | 回復(fù): 3 | ||
tjj_wei鐵蟲 (初入文壇)
|
[求助]
蒸發(fā)的程序出錯(cuò),順便請(qǐng)大神指導(dǎo)下C_YI的用法
|
|
#include "udf.h" //包括常規(guī)宏 #include "sg_mphase.h" //包括體積分?jǐn)?shù) #include "mem.h" //包括溫度及組分質(zhì)量分?jǐn)?shù) #include "sg.h" #include "flow.h" #define T_SAT 373 //定義蒸發(fā)溫度100℃ #define LAT_HT 1.e3 DEFINE_SOURCE(liq_src,cell,sec_th,dS,eqn) { Thread *mix_th,*pri_th; //定義計(jì)算區(qū)域指針 real m_dot_l; //定義液相質(zhì)量轉(zhuǎn)移kg/(m2*s) mix_th=THREAD_SUPER_THREAD(sec_th); pri_th=THREAD_SUB_THREAD(mix_th,2); if(C_T(cell,mix_th)>=T_SAT) { m_dot_l=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } else { m_dot_l=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } return m_dot_l; } DEFINE_SOURCE(vap_src,cell,pri_th,dS,eqn) { Thread *mix_th,*sec_th; //定義計(jì)算區(qū)域指針 real m_dot_v; //定義氣相質(zhì)量轉(zhuǎn)移kg/(m2*s) mix_th=THREAD_SUPER_THREAD(pri_th); sec_th=THREAD_SUB_THREAD(mix_th,3); if(C_T(cell,mix_th)>=T_SAT) { m_dot_v=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } else { m_dot_v=-0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_YI(cell,pri_th,2)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } return m_dot_v; } DEFINE_SOURCE(enrg_src,cell,mix_th,dS,eqn) { Thread *pri_th,*sec_th; //定義計(jì)算區(qū)域指針 real m_dot; //定義液相質(zhì)量轉(zhuǎn)移kg/(m2*s) pri_th=THREAD_SUPER_THREAD(mix_th,2); sec_th=THREAD_SUB_THREAD(mix_th,3); if(C_T(cell,mix_th)>=T_SAT) { m_dot=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)/T_SAT; } else { m_dot=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)/T_SAT; } return LAT_HT*m_dot; } DEFINE_SOURCE(h2o_src,cell,pri_th,dS,eqn) { Thread *mix_th,*sec_th; //定義計(jì)算區(qū)域指針 real m_dot_vw; //定義氣相質(zhì)量轉(zhuǎn)移kg/(m2*s) mix_th=THREAD_SUPER_THREAD(pri_th); sec_th=THREAD_SUB_THREAD(mix_th,3); if(C_T(cell,mix_th)>=T_SAT) { m_dot_vw=0.1*C_VOF(cell,sec_th)*C_R(cell,sec_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=0; } else { m_dot_vw=-0.1*C_YI(cell,pri_th,2)*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; dS[eqn]=-0.1*C_VOF(cell,pri_th)*C_R(cell,pri_th)*fabs(C_T(cell,mix_th)-T_SAT)/T_SAT; } return m_dot_vw; } 這是我的程序的內(nèi)容,一初始化就出錯(cuò),不知道是不是我的C_YI事業(yè)有問題,還有第一相是0還是查的ID號(hào)啊?請(qǐng)大家?guī)臀铱纯? |
鐵蟲 (初入文壇)
新蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|---|---|---|---|---|
|
[基金申請(qǐng)] 山東省面上項(xiàng)目限額評(píng)審 +4 | 石瑞0426 2026-03-19 | 4/200 |
|
|
[考研] 求調(diào)劑 +3 | .m.. 2026-03-21 | 4/200 |
|
|
[基金申請(qǐng)]
學(xué)校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
[考研]
|
然11 2026-03-19 | 4/200 |
|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +22 | rare12345 2026-03-18 | 22/1100 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 一志愿蘇州大學(xué)材料工程(085601)專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|