 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1562 | 回復(fù): 15 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
幫忙看下這個udf有問題嗎？初始化沒問題，點計算出錯
|
||
|
幫忙看下這個udf有問題嗎？初始化沒問題，點計算出錯 #include "udf.h" /*包括常規(guī)宏*/ #include "sg_mphase.h" /*包括體積分?jǐn)?shù)宏*/ #define LAT_HT 2260000 /*定義水蒸汽潛熱*/ #define pi 3.1415 /*圓周率常數(shù)*/ #define R 8314 /*gas constant 8.314j/mol.K*/ #define moleweight 0.018 /*SI unit water mole weight 0.018kg/mol */ #define A 7.19621 #define B 1730.63 #define C 233.426 DEFINE_SOURCE(liq_src, cell, pri_th, dS, eqn) /*定義液相源項*/ { Thread *mix_th, *sec_th; /*混合相、第二相 定義計算區(qū)域指針*/ real m_dot_l; /*蒸發(fā)冷凝速率 定義液相質(zhì)量轉(zhuǎn)移*/ real P_SAT ; mix_th = THREAD_SUPER_THREAD(pri_th); /*指向混合區(qū)的主相即液相的指針*/ sec_th = THREAD_SUB_THREAD(mix_th, 1); /* 指向單相控制區(qū)的氣相的指針*/ P_SAT=exp(A-B/(C_T(cell,pri_th)+C-273.15))*1e3; /*定義飽和蒸汽壓與液相溫度的關(guān)系式，即antoine公式*/ if(C_P(cell, sec_th)< P_SAT) {m_dot_l = -(P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); /*定義液相向氣相轉(zhuǎn)移的速率，如果氣相的壓力小于一定溫度下的飽和蒸汽壓，液相向氣相發(fā)生質(zhì)量轉(zhuǎn)移*/ dS[eqn]= -2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th)); /*定義源項對質(zhì)量轉(zhuǎn)移偏導(dǎo)，fabs為取絕對值*/ } else {m_dot_l = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); /*定義液相向氣相轉(zhuǎn)移的速率，如果氣相壓力大于一定溫度下的飽和蒸汽壓，氣相向液相發(fā)生質(zhì)量轉(zhuǎn)移*/ dS[eqn] = 0; /*由于氣相向液相轉(zhuǎn)移，所以液相的質(zhì)量源項對質(zhì)量轉(zhuǎn)移的偏導(dǎo)為零*/ } return m_dot_l; } DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) { Thread *mix_th, *pri_th; real m_dot_v; real P_SAT ; mix_th = THREAD_SUPER_THREAD(sec_th); pri_th = THREAD_SUB_THREAD(mix_th, 0); P_SAT=exp(A-B/(C_T(cell,pri_th)+C-273.15))*1e3; if(C_P(cell, sec_th)< P_SAT) { m_dot_v = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); dS[eqn] =0; } else {m_dot_v = -(P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); dS[eqn]=-2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th)); } return m_dot_v; } DEFINE_SOURCE(enrg_src, cell, mix_th, dS, eqn) { Thread *pri_th,*sec_th; real m_dot; real P_SAT; pri_th=THREAD_SUB_THREAD(mix_th,0); sec_th=THREAD_SUB_THREAD(mix_th,1); P_SAT=exp(A-B/(C_T(cell,pri_th)+C-273.15))*1e3; if(C_P(cell, sec_th)< P_SAT) { m_dot =- (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); dS[eqn]=-2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th)); } else { m_dot = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); dS[eqn]=2*pi*R*sqrt(moleweight)*(P_SAT-P_SAT*C_T(cell,pri_th)*B/pow((C_T(cell,pri_th)+C-273.15),2)-C_P(cell, sec_th)); } return LAT_HT*m_dot; } |
回復(fù)此樓» 猜你喜歡
材料調(diào)劑
已經(jīng)有5人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有8人回復(fù)
-
請教下大家 2026年國家基金申請是雙盲審嗎？
已經(jīng)有5人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷
已經(jīng)有10人回復(fù)
-
一志愿南航材料專317分求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- FLUENT中編寫UDF時指針變量是否要初始化
已經(jīng)有4人回復(fù)
- fluent14編譯UDF時總是出錯，求高手啊�。�！難道沒人遇到過嗎？？？
已經(jīng)有11人回復(fù)
- 求助：模型初始化時出錯
已經(jīng)有5人回復(fù)
- fluent udf請教
已經(jīng)有16人回復(fù)
- 一個UDF并行報錯
已經(jīng)有17人回復(fù)
- 導(dǎo)入UDF時出這種錯誤��！
已經(jīng)有4人回復(fù)
- UDF 改為并行計算
已經(jīng)有9人回復(fù)
- FLUENT UDF 編譯通過，初始化出錯的原因？？？求助
已經(jīng)有4人回復(fù)
- 求助�。�！我想要在fluent中針對一個面進(jìn)行積分，那么UDF應(yīng)該用什么宏函數(shù)？怎么用？
已經(jīng)有5人回復(fù)
- 好心人幫幫忙吧，為什么這個UDF導(dǎo)入fluent以后，在邊界條件處沒有顯示~~~
已經(jīng)有17人回復(fù)
- 怎樣解決FLUENT received fatal signal (ACCESS_VIOLATION)的錯誤？
已經(jīng)有20人回復(fù)
- udf編譯出錯！急
已經(jīng)有10人回復(fù)
- udf編程通過，但是在FLUENT中運行出錯，麻煩udf高手幫忙看看
已經(jīng)有4人回復(fù)
- 使用UDF后，初始化錯誤
已經(jīng)有11人回復(fù)
- 【求助】fluent初始化的時候出現(xiàn)問題，請大家指點，在線等，多謝各位了
已經(jīng)有4人回復(fù)
- 【求助】fluent計算中，udf編譯通過，初始化出錯
已經(jīng)有11人回復(fù)
- 【求助】udf出錯
已經(jīng)有13人回復(fù)
轉(zhuǎn)石千仞
金蟲 (正式寫手)
擊破長空
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 941
- 散金: 164
- 紅花: 13
- 帖子: 673
- 在線: 161.6小時
- 蟲號: 1163587
- 注冊: 2010-12-06
- 性別: GG
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
4樓2013-09-05 10:52:09
2樓2013-09-05 10:39:39
轉(zhuǎn)石千仞
金蟲 (正式寫手)
擊破長空
- 應(yīng)助: 24 (小學(xué)生)
- 金幣: 941
- 散金: 164
- 紅花: 13
- 帖子: 673
- 在線: 161.6小時
- 蟲號: 1163587
- 注冊: 2010-12-06
- 性別: GG
- 專業(yè): 工程熱物理相關(guān)交叉領(lǐng)域
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
837729744: 金幣+8 2013-09-05 11:04:03
1592203609: 金幣+3, 謝謝回帖 2013-09-05 16:01:15
感謝參與，應(yīng)助指數(shù) +1
837729744: 金幣+8 2013-09-05 11:04:03
1592203609: 金幣+3, 謝謝回帖 2013-09-05 16:01:15
|
這種情況下，建議你對計算的數(shù)據(jù)進(jìn)行檢查。 你可以采用程序截斷的形式，用Message插入可能出問題的位置前后，看是否能輸出語言與數(shù)據(jù)，這樣可以檢測你的程序運行到哪里出現(xiàn)了浮點錯誤。 你還可以用Error語句來控制程序的停止，你可以讓程序運行到任意位置停止，以便不讓程序在每一個控制容積里運行。 出現(xiàn)浮點錯誤的原因，是你計算出現(xiàn)極大或極小值，比0等等。 調(diào)試程序需要耐心與經(jīng)驗，祝你好運！ |
3樓2013-09-05 10:51:03
5樓2013-09-05 11:03:37
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +8 | 迷人的哈哈 2026-03-23 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +4 | 炸呀炸呀炸薯條 2026-03-23 | 4/200 |
|
|
|
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +11 | 手機用戶 2026-03-17 | 12/600 |
|
|
|
[考研] 279分求調(diào)劑 一志愿211 +17 | chaojifeixia 2026-03-19 | 19/950 |
|
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
|
[考研] 316求調(diào)劑 +7 | 梁茜雯 2026-03-19 | 7/350 |
|
|
|
[考研] 328求調(diào)劑 +4 | LHHL66 2026-03-23 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
|
[考研] 尋找調(diào)劑 +4 | 倔強芒? 2026-03-21 | 4/200 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研]
 材料與化工 322求調(diào)劑
+4
|
然11 2026-03-19 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|