 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1095 | 回復: 1 | |||
| 【懸賞金幣】回答本帖問題，作者taozhenxiang將贈送您 20 個金幣 | |||
taozhenxiang新蟲 (小有名氣)
|
[求助]
flunet 濕空氣冷凝UDF浮點數(shù)溢出
|
||
 各位好！小弟用網(wǎng)上的資源貼做了濕空氣在低溫環(huán)境下的冷凝，初始化后計算一直抱浮點數(shù)溢出的錯誤�？紤]有這幾方面的問題1.網(wǎng)格質(zhì)量差；2.調(diào)整松弛因子；3.調(diào)整迭代步長4.邊界條件不合理。5.UDF代碼存在錯誤。１－－４本人已經(jīng)做了多次修改，還是報同樣的錯誤，因而推斷為UDF本身的問題。自己剛接觸UDF，并未發(fā)現(xiàn)明顯錯誤，還請各位前輩給予指點，具體代碼如下： #include "udf.h" /*包括常規(guī)宏*/ #include "sg_mphase.h" /*包括體積分數(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ì)量轉移*/ 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; /*定義飽和蒸汽壓與液相溫度的關系式，即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)); /*定義液相向氣相轉移的速率，如果氣相的壓力小于一定溫度下的飽和蒸汽壓，液相向氣相發(fā)生質(zhì)量轉移*/ 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ì)量轉移偏導，fabs為取絕對值*/ } else {m_dot_l = (P_SAT-C_P(cell,sec_th))*sqrt(moleweight)*(2*pi*R*C_T(cell,pri_th)); /*定義液相向氣相轉移的速率，如果氣相壓力大于一定溫度下的飽和蒸汽壓，氣相向液相發(fā)生質(zhì)量轉移*/ dS[eqn] = 0; /*由于氣相向液相轉移，所以液相的質(zhì)量源項對質(zhì)量轉移的偏導為零*/ } 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; } Thanks in advance |
回復此樓» 猜你喜歡
幫我的英語口語老師找學生
已經(jīng)有4人回復
-
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有132人回復
-
2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復
-
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有2人回復
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復
-
求調(diào)劑
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復
taozhenxiang
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 3487.8
- 帖子: 212
- 在線: 208.4小時
- 蟲號: 2446091
- 注冊: 2013-05-03
- 性別: GG
- 專業(yè): 工程熱力學
2樓2020-08-10 17:11:31
相關版塊跳轉
我要訂閱樓主 taozhenxiang 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 283求調(diào)劑，材料、化工皆可 +5 | 蘇打水7777 2026-03-11 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 材料工程,326分，求調(diào)劑 +5 | KRSLSR 2026-03-10 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +4 | 步川酷紫123 2026-03-11 | 4/200 |
|
|
[考研] 求調(diào)劑 +5 | yfihxh 2026-03-09 | 5/250 |
|
|
|
[考研] 一志愿安徽大學材料工程專碩313分，求調(diào)劑的學校 +7 | Yu先生 2026-03-10 | 9/450 |
|
|
|
[考研] 083000環(huán)境科學與工程調(diào)劑 +8 | mingmingry 2026-03-09 | 9/450 |
|
|
|
[基金申請] 提交后的基金本子，已讓學校撤回了，可否換口子提交 +3 | dut_pfx 2026-03-10 | 3/150 |
|
|
|
[考研] 一志愿天津大學，英一數(shù)二305分求調(diào)劑，四六級已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 牛乳糖的卡卡 2026-03-10 | 3/150 |
|
|
|
[考研] 0703化學調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
|
[考博] 2026博士申請 +6 | 起泡酒 2026-03-08 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 +7 | 小奶蒼 2026-03-07 | 7/350 |
|
|
|
[基金申請] 面上項目還需要AI說明嗎？ +3 | liyundong 2026-03-08 | 3/150 |
|
|
|
[考研] 0832食品科學與工程293調(diào)劑 +3 | 東東不東 2026-03-07 | 3/150 |
|
|
|
[考研] 293一志愿華東理工 0817化學工程與技術 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
|
[考研] 理學287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
|
[考研] 085701環(huán)境工程專業(yè)，初試305，均過國家A區(qū)線 +7 | 卡卡來了@ 2026-03-07 | 8/400 |
|
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 材料與化工354調(diào)劑 +4 | Lucy-xiao 2026-03-06 | 7/350 |
|