IcePak是不是存在許多bug?計算熱輻射问题

我建了一個cabnet(1m*1m*1m)六個邊界面設為open,固定20度C.  
在cabnet中心放一個Block(25mm*10mm*1mm), 即面積為0.00057m^2. 固定發熱0.2654Watt.

在不計熱輻射的情況下計算得block溫度44度C; h=12.7(mean); heat flux=303(mean).
但是打開full report檢查h=1.96(mean); heat flux=47.06(mean).

在計算熱輻射的情況下計算得block溫度69度C; h=20.7(mean); heat flux=1010(mean).
但是打開full report檢查h=29.11(mean); heat flux=1427.12(mean);

第一個不合理是,在計算熱輻射的情況block溫度不應該比不計熱輻射的情況高
第二個不合理是,固定發熱0.2654Watt,散熱卻不能與之平衡. 其分別為

不計熱輻射
1a:                   h*A*dT=12.7*0.00057*24=0.174Watt;   heat flux*A=303*0.00057=0.173Watt
1b(full report):  h*A*dT=1.96*0.00057*24=0.0269Watt;  heat flux*A=47.06*0.00057=0.0268Watt

計算熱輻射
2a:                    h*A*dT=20.7*0.00057*49=0.578Watt;   heat flux*A=1010*0.00057=0.578Watt
2b(full report):  h*A*dT=29.11*0.00057*49=0.817Watt;  heat flux*A=1427.12*0.00057=0.813Watt ;radiation heat flow=0.1635Watt

第三個不合理是,full report的summary的值竟然和主視窗右下角所顯示不同

有沒有人能為我解惑?
包裝檔(.tzr)如附件

radiation.rar (28.73 KB)
.tzr files (no-radiation.tzr and with-radiation.tzr)

icepak 的輻射熱傳還有一個bug
上一帖提到icepak的三個bugs
本帖再補述另外一個bug

在計算view factor之後,我把
cabnet surface to cabnet surface 以及 cabnet surface to block surface  的 view factor "modify" 為零
其結果竟然和view factor 未經修改的結果相同---------block溫度=69度C

總結
IcePak的輻射計算似乎存在嚴重的bug

希望IcePak原廠重視此一問題

这么多bug的软件居然可以卖几十甚至上百万；这么低级的错误都有的软件居然敢在行业中称老大；这么多的业内高手居然没有发现这些bug；

不太可能把？

比较欣赏你敢于怀疑一切的精神，但是挑战权威时自己要做充分的分析阿。

你的case我看了一下，应该是网格化分不够细致造成的。

cabnet你不需要要这么大的，50*50*50cm就可以了。

没有辐射的话，温度在48.5度；

有辐射时温度在42.3度；辐射散热量在0.065W左右。

lzzmn:
萬分感謝您花費寶貴的時間看我的算例
不知能否得到您算例的tzr檔?

的確, mesh影響結果甚巨.
請教沒有輻射時您得到多少total heat flux?---我很想知道.      ----我的heat source是0.2654W
有輻射時您得到多少對流熱傳?------------是不是0.2654-0.065=0.2004W?
因為我一直很在乎熱傳的平衡...........

我是最近才開始學習icepak
尤其網格的處理我還很嫩 希望您多指教

下午我重劃網格且改用cartesian得到下列結果
没有辐射的话，温度在53.5度；
有辐射时温度在43.5度；辐射散热量在0.0759W左右。不過對流熱傳卻高達0.2786W。

0.0759+0.2786=0.3545W>>0.2654W

對於此一結果我仍不能釋懷

我把此一tzr檔上傳予你，希望你有良策

明天我會將 cabnet 改成您建議的大小，再看結果如何.......



test3.rar (16.64 KB)

我觉得ICEPAK中所谓的HEAT FLOW应该是指传热系数（包括对流，辐设）。对该算例重新计算，修改如下：热源网格加密（首层网格0.2mm），under-relaxation改为0.7－0.3，flow regime改为laminar，采用双精度模型. 计算结果如下图：

  
下载 (38.01 KB)



CFD计算结果对网格，计算模型和算法等因素比较敏感，模型计算时需要多加注意。

Storm/Chaos:
請問你們用了多少網格

大概16060个网格，附件是TZR文件。其实你应该把ICEPAK的基础培训好好学习学习，很有好处的。  

icepak.rar (16.86 KB)


screen.gif (18.62 KB)


Storm
謝謝您的.tzr file.
/我很驚訝您只用了16060個網格
我正朝減少網格數的方向努力

lzzmn:
經過三位高手指導我已漸漸入門了
對icepak的信心度急速提高
再次感謝三位


请问如何在计算前，预防这种情况的发生。因为如果模型复杂的话，会很浪费时间。