 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1090 | 回復(fù): 1 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者taozhenxiang將贈(zèng)送您 20 個(gè)金幣 | ||
taozhenxiang新蟲 (小有名氣)
|
[求助]
flunet 濕空氣冷凝UDF浮點(diǎn)數(shù)溢出
|
|
 各位好！小弟用網(wǎng)上的資源貼做了濕空氣在低溫環(huán)境下的冷凝，初始化后計(jì)算一直抱浮點(diǎn)數(shù)溢出的錯(cuò)誤�？紤]有這幾方面的問(wèn)題1.網(wǎng)格質(zhì)量差；2.調(diào)整松弛因子；3.調(diào)整迭代步長(zhǎng)4.邊界條件不合理。5.UDF代碼存在錯(cuò)誤。１－－４本人已經(jīng)做了多次修改，還是報(bào)同樣的錯(cuò)誤，因而推斷為UDF本身的問(wèn)題。自己剛接觸UDF，并未發(fā)現(xiàn)明顯錯(cuò)誤，還請(qǐng)各位前輩給予指點(diǎn)，具體代碼如下： #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) /*定義液相源項(xiàng)*/ { Thread *mix_th, *sec_th; /*混合相、第二相 定義計(jì)算區(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)); /*定義源項(xiàng)對(duì)質(zhì)量轉(zhuǎn)移偏導(dǎo)，fabs為取絕對(duì)值*/ } 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ì)量源項(xiàng)對(duì)質(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; } Thanks in advance |
回復(fù)此樓» 猜你喜歡
重慶交通大學(xué)光子學(xué)微結(jié)構(gòu)與器件課題組招收2026年碩士研究生信息
已經(jīng)有2人回復(fù)
-
一志愿鄭大材料學(xué)碩298分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有282人回復(fù)
-
尋合作：應(yīng)力腐蝕多尺度模擬
已經(jīng)有6人回復(fù)
-
考研交流
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語(yǔ)口語(yǔ)老師找學(xué)生
已經(jīng)有3人回復(fù)
taozhenxiang
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3487.8
- 帖子: 212
- 在線: 208.4小時(shí)
- 蟲號(hào): 2446091
- 注冊(cè): 2013-05-03
- 性別: GG
- 專業(yè): 工程熱力學(xué)
2樓2020-08-10 17:11:31
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 taozhenxiang 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +3 | 52hz~~ 2026-03-05 | 4/200 |
|
|---|---|---|---|---|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 5/250 |
|
|
|
[考研] 303求調(diào)劑 +3 | forgman95 2026-03-05 | 5/250 |
|
|
|
[考研] 283求調(diào)劑 +10 | 鹿沫笙 2026-03-02 | 11/550 |
|
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
|
[考研] 材料085600 303求調(diào)劑 +7 | 1bygone 2026-03-04 | 7/350 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 成績(jī)276，專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +8 | 醋精華了一下發(fā)?/a> 2026-03-03 | 10/500 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|
|
|
[考研] 0805總分292，求調(diào)劑 +12 | 幻想之殤 2026-03-01 | 12/600 |
|
|
|
[考研] 一志愿華中科技大學(xué)，化學(xué)專業(yè)344分，求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
|
[考研] 268求調(diào)劑 +10 | 簡(jiǎn)單點(diǎn)0 2026-03-02 | 14/700 |
|
|
|
[考研] 271求調(diào)劑 +4 | Ricardo1113 2026-03-02 | 4/200 |
|
|
|
[考研] 298求調(diào)劑 +7 | axyz3 2026-02-28 | 8/400 |
|
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。，。，。，。i 2026-03-02 | 3/150 |
|
|
|
[考研] 284求調(diào)劑 +10 | 天下熯 2026-02-28 | 11/550 |
|
|
|
[考研] 一志愿鄭大材料學(xué)碩298分，求調(diào)劑 +6 | wsl111 2026-03-01 | 6/300 |
|