 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 3 | 1/1 | 返回列表 |
| 查看: 1296 | 回復(fù): 2 | |||
[求助]
相變udf求助
|
|
最近在做相變的模擬，udf使用的是論壇上大家通用的，我只是改了一下指針，結(jié)果就老是初始化報錯，把udf放在下面，希望大家給看一下啊，謝謝各位，我個人感覺是溫度的指針出了問題。 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; } |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 博士科研新路歷程 |
» 猜你喜歡
309求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿南航 335分 | 0856 | GPA 4.07 | 有科研經(jīng)歷
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
求化學調(diào)劑
已經(jīng)有5人回復(fù)
-
085600，材料與化工321分求調(diào)劑
已經(jīng)有10人回復(fù)
-
0703 化學 求調(diào)劑，一志愿山東大學 342 分
已經(jīng)有5人回復(fù)
-
一志愿南開大學0710生物學359求調(diào)劑
已經(jīng)有3人回復(fù)
-
085600，專業(yè)課化工原理，320分求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工272求調(diào)劑
已經(jīng)有16人回復(fù)
-
290求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求大神幫檢查下一個簡單的UDF
已經(jīng)有13人回復(fù)
- 有關(guān)fluent UDF的問題
已經(jīng)有15人回復(fù)
- 關(guān)于溫度的傳遞問題，能否有人幫我看一下這個UDF是哪里出了錯誤呢？
已經(jīng)有10人回復(fù)
- Ansys Fluent 編譯UDF及加載UDF出現(xiàn)錯誤
已經(jīng)有14人回復(fù)
- 急，在FLUENT中如何同時調(diào)用兩個UDF？
已經(jīng)有11人回復(fù)
- Fluent中DPM模型的UDF中使用P_POS()指針，然后給我提示undeclared variable 怎么辦?
已經(jīng)有6人回復(fù)
- 請問誰做過氣液相變的fluent模擬呀
已經(jīng)有3人回復(fù)
- fluent中相變過程怎樣設(shè)置固液密度
已經(jīng)有10人回復(fù)
- fluent UDF 修改物性
已經(jīng)有13人回復(fù)
- UDF中如何編寫體積分數(shù)的梯度
已經(jīng)有9人回復(fù)
- fluent中UDF編譯的問題！
已經(jīng)有6人回復(fù)
- 關(guān)于水蒸發(fā)udf
已經(jīng)有19人回復(fù)
- 相變模型 UDF 求助
已經(jīng)有9人回復(fù)
- Fluent中UDF的編譯問題
已經(jīng)有25人回復(fù)
- 【求助】FLUENT UDF 氣液兩相流 汽蝕 汽化潛熱(問題詳見內(nèi)容)
已經(jīng)有21人回復(fù)
- UDF應(yīng)用之后的初始化錯誤
已經(jīng)有6人回復(fù)
|
本帖內(nèi)容被屏蔽 |
2樓2020-03-12 19:06:17
zhangcx1208
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 6889.3
- 散金: 51
- 紅花: 5
- 帖子: 309
- 在線: 112.4小時
- 蟲號: 468285
- 注冊: 2007-11-27
- 專業(yè): 流體力學
3樓2020-12-26 11:19:55
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 gky227227 的主題更新
| 3 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿雙一流機械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
|
[考研] 數(shù)一英一271專碩（085401）求調(diào)劑，可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 材料277求調(diào)劑 +7 | min3 2026-03-24 | 7/350 |
|
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 339求調(diào)劑，想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
|
[考研] 0703化學求調(diào)劑，各位老師看看我！��！ +5 | 祁祺祺 2026-03-25 | 5/250 |
|
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
|
[考研] 0703化學338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關(guān)調(diào)劑 +5 | 星沉uu 2026-03-26 | 6/300 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|