| 24小時(shí)熱門(mén)版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 1295 | 回復(fù): 2 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
相變udf求助
|
|||
|
最近在做相變的模擬,udf使用的是論壇上大家通用的,我只是改了一下指針,結(jié)果就老是初始化報(bào)錯(cuò),把udf放在下面,希望大家給看一下啊,謝謝各位,我個(gè)人感覺(jué)是溫度的指針出了問(wèn)題。 DEFINE_SOURCE(liq_src, cell, thread, dS, eqn) { /*Thread *mix_th, *sec_th;*/ Thread *mix_thread, *thread_gas, *thread_liq; real m_dot_l; real T_SAT; /*定義臨界溫度*/ real tem_pressure; /*定義單元壓力*/ real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix; mix_thread= THREAD_SUPER_THREAD(thread); thread_gas= THREAD_SUB_THREAD(mix_thread, 1); thread_liq= THREAD_SUB_THREAD(mix_thread, 0); vof_gas=C_VOF(cell, thread_gas); vof_liq=C_VOF(cell, thread_liq); roul_gas=C_R(cell, thread_gas); roul_liq=C_R(cell, thread_liq); t_gas=C_T(cell, thread_gas); t_liq=C_T(cell, thread_liq); t_mix=C_T(cell, mix_thread); /*mix_th = THREAD_SUPER_THREAD(pri_th); sec_th = THREAD_SUB_THREAD(mix_th, 1);*/ tem_pressure=(C_P(cell,mix_thread)+101325)/1000000; /*讀取單元壓力*/ T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47; /*定義飽和溫度與壓力的關(guān)系*/ if(t_liq>=T_SAT) { m_dot_l = -0.1*vof_liq*roul_liq*fabs(t_liq - T_SAT)/T_SAT; dS[eqn] = -0.1*roul_liq*fabs(t_liq - T_SAT)/T_SAT; } { m_dot_l = 0.1*vof_gas*roul_gas*fabs(t_liq - T_SAT)/T_SAT; dS[eqn] =0; } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, thread, dS, eqn) { Thread *mix_thread, *thread_gas, *thread_liq; real m_dot_v; real T_SAT; /*定義臨界溫度*/ real tem_pressure; /*定義單元壓力*/ real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix; mix_thread= THREAD_SUPER_THREAD(thread); thread_gas= THREAD_SUB_THREAD(mix_thread, 1); thread_liq= THREAD_SUB_THREAD(mix_thread, 0); vof_gas=C_VOF(cell, thread_gas); vof_liq=C_VOF(cell, thread_liq); roul_gas=C_R(cell, thread_gas); roul_liq=C_R(cell, thread_liq); t_gas=C_T(cell, thread_gas); t_liq=C_T(cell, thread_liq); t_mix=C_T(cell, thread); tem_pressure=(C_P(cell,mix_thread)+101325)/1000000; /*讀取單元壓力*/ T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47; /*定義飽和溫度與壓力的關(guān)系*/ if(t_gas>=T_SAT) { m_dot_v = 0.1*vof_liq*roul_liq* fabs(t_gas - T_SAT)/T_SAT; dS[eqn] = 0; /*1*/ } else { m_dot_v = -0.1*vof_liq*roul_liq*fabs(T_SAT-t_gas)/T_SAT; dS[eqn] = -0.1*roul_liq*fabs(t_gas - T_SAT)/T_SAT; } return m_dot_v; } DEFINE_SOURCE(enrg_src, cell, thread, dS, eqn) { Thread *mix_thread, *thread_gas, *thread_liq; real m_dot_l; real vof_gas,vof_liq,roul_gas,roul_liq,t_gas,t_liq,t_mix; real T_SAT; /*定義臨界溫度*/ real tem_pressure; /*定義單元壓力*/ mix_thread= THREAD_SUPER_THREAD(thread); thread_gas= THREAD_SUB_THREAD(mix_thread, 1); thread_liq= THREAD_SUB_THREAD(mix_thread, 0); vof_gas=C_VOF(cell, thread_gas); vof_liq=C_VOF(cell, thread_liq); roul_gas=C_R(cell, thread_gas); roul_liq=C_R(cell, thread_liq); t_gas=C_R(cell, thread_gas); t_liq=C_R(cell, thread_liq); t_mix=C_T(cell, thread); tem_pressure=(C_P(cell,mix_thread)+101325)/1000000; /*讀取單元壓力*/ T_SAT=3826.36/(9.3876-log(tem_pressure))+45.47; /*定義飽和溫度與壓力的關(guān)系*/ if(t_mix>=T_SAT) { m_dot_l = -0.1*vof_liq*roul_liq* fabs(t_mix - T_SAT)/T_SAT; dS[eqn] =-0.1*vof_liq*roul_liq/T_SAT; } else { m_dot_l = 0.1*vof_gas*roul_gas* fabs(t_mix - T_SAT)/T_SAT; dS[eqn] =0.1*vof_gas*roul_gas/T_SAT; /*1*/ } return LAT_HT*m_dot_l; } |
博士科研新路歷程 |
鐵桿木蟲(chóng) (正式寫(xiě)手)
|
本帖內(nèi)容被屏蔽 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +4 | QiMing7 2026-03-25 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 化學(xué)0703 調(diào)劑 306分 一志愿211 +7 | 26要上岸 2026-03-28 | 7/350 |
|
|
[考研] 求調(diào)劑 +9 | 張zz111 2026-03-27 | 10/500 |
|
|
[考研] 327求調(diào)劑 +6 | 汲亦昊 2026-03-29 | 6/300 |
|
|
[考研] 356求調(diào)劑 +3 | gysy?s?a 2026-03-28 | 3/150 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
[考研] 調(diào)劑 +3 | 好好讀書(shū)。 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
[考研] 070300化學(xué)求調(diào)劑 +4 | 起個(gè)名咋這么難 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿華東理工大學(xué)081700,初試分?jǐn)?shù)271 +6 | kotoko_ik 2026-03-23 | 7/350 |
|
|
[考研] 343求調(diào)劑 +4 | 贈(zèng)我一本書(shū) 2026-03-23 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 327求調(diào)劑 +7 | prayer13 2026-03-23 | 7/350 |
|
|
[考研] 中國(guó)科學(xué)院深圳先進(jìn)技術(shù)研究院-光纖傳感課題組招生-中國(guó)科學(xué)院大學(xué)、深圳理工大學(xué)聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 26考研-291分-廈門(mén)大學(xué)(085601)-柔性電子學(xué)院材料工程專(zhuān)業(yè)求調(diào)劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|