| 4 | 1/1 | 返回列表 |
| 查看: 1082 | 回復(fù): 3 | |||
tjj_wei鐵蟲 (初入文壇)
|
[求助]
蒸發(fā)的程序出錯,順便請大神指導(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; //定義計算區(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; //定義計算區(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; //定義計算區(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; //定義計算區(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)容,一初始化就出錯,不知道是不是我的C_YI事業(yè)有問題,還有第一相是0還是查的ID號。空埓蠹?guī)臀铱纯? |
鐵蟲 (初入文壇)
新蟲 (初入文壇)
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料科學(xué)與工程調(diào)劑 +6 | 深V宿舍吧 2026-03-30 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 310求調(diào)劑 +9 | 爭取九點睡 2026-03-30 | 9/450 |
|
|
[考研] 材料與化工272求調(diào)劑 +21 | 阿斯蒂芬2004 2026-03-28 | 21/1050 |
|
|
[考研] 295求調(diào)劑 +5 | wei-5 2026-03-26 | 5/250 |
|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|
[考研] 281求調(diào)劑 +4 | 亞克西good 2026-03-26 | 6/300 |
|
|
[考研] 289求調(diào)劑 +5 | BrightLL 2026-03-29 | 5/250 |
|
|
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
[考研] 328求調(diào)劑 +7 | 嗯滴的基本都 2026-03-27 | 7/350 |
|
|
[考研] 315分求調(diào)劑 +7 | 26考研上岸版26 2026-03-26 | 7/350 |
|
|
[考研] 322求調(diào)劑 +4 | 我真的很想學(xué)習(xí) 2026-03-23 | 4/200 |
|
|
[考研] 材料求調(diào)劑 +5 | .m.. 2026-03-25 | 5/250 |
|
|
[考研] 中國科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 07化學(xué)303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研]
|
黃粱一夢千年 2026-03-24 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|