1、微機電系統分析報告-Inverse heat transfer,Introduction,在實際熱傳導問題中，若其初始條件(Initial condition)、邊界條件(Boundary condition)與材料性質(Material Properties)為已知的熱物理量時，可直接由熱傳導方程式去進行求解整個溫度場分布，此類型的熱傳問題稱為直接熱傳導問題(Direct Heat Conduction Problem DHCP)；,Introduction,然而實際工程問題常因客觀的條件限制或是量測技術的不足，無法得到完整的邊界條件，因此必須藉由量測熱導體內部之溫度數據，再經由逆向估算求出
2、未知邊界條件，此類型的熱傳問題稱為逆向熱傳導問題(Inverse Heat Conduction Problem IHCP)，經由逆向估算出未知邊界條件後，即可求解整個溫度場之溫度分布1-5。,文獻回顧,Estimation of Surface Absorptivity and Surface Temperature1(下圖擷取自此篇),文獻回顧,在1950年代後期，逆向熱傳導問題開始受到學界的注意與研究。依空氣動力學 4 理論可知，當太空設備返回地球時，設備表面將因有極嚴重的熱量產生以致於無任何感測器可殘存於其表面。因此發展出將感測器埋入物體內部而反推算出物體表面邊界條件的逆運算技術。,I
3、nverse應用-國防軍事科技5,底火瞬間點燃時所釋放出之熱通膛內彈道中槍砲管內壁熱通及溫計算藉由槍管外壁變形之測逆估算槍膛內壁壓反艦飛彈發射時推進系統之燃氣噴火焰溫對甲板之影響垂直起戰機之發動瞬間所產生之高溫噴熔裝固化過程之接觸熱傳係等皆屬逆向熱傳導問題之範疇。,概念說明,Direct ProblemKnown:a1、a2、b1、b2、Q、RUnknown:X、Y,概念說明,Inverse ProblemKnown:a1、a2、b1、b2、Y、RUnknown:X、Q,概念說明-熱傳問題6,Finite Difference Equation,擷取自6,熱傳問題-Direct,Known:T
4、1 T12Unknown:TA1 TA4利用Finite Difference Equation寫出關係式,T2+T12+TA3+TA2-4TA1=0 1 T3+T5+TA1+TA4-4TA2=0 2 T11+T9+TA1+TA4-4TA3=0 3 T6+T8+TA2+TA3-4TA4=0 4-4T1+T2+T12=0 5 T1-4T2+T3+TA1=0 6 T2-4T3+T4+TA2=0 7 T3-4T4+T5=0 8 T4-4T5+T6+TA2=0 9T5-4T6+T7+TA4=0 10T6-4T7+T8=0 11T7-4T8+T9+TA4=0 12T8-4T9+T10+TA3=0 13T
5、9-4T10+T11=0 14T10-4T11+T12+TA3=0 15T1+T11-4T12+TA1=0 16,熱傳問題-Inverse,Known:T3 T12、TA1、TA2Unknown:T1、T2、TA3、TA4利用Finite Difference Equation寫出關係式,Inverse method,敏感性(sensitivity)與最小平方法(the least square method)共軛梯度法(Conjugate Gradient Method)Modified Conjugate Gradient MethodB-Spline Method求解逆向熱傳導問題，最小
6、平方法被廣泛使用，因其精確度和穩定性均能合乎要求。,敏感性係數(Sensitivity coefficient),敏感性係數是假設系統在線性系統之下，利用疊加法則(superposition)對其邊界條件去進行展開；是一種探討已知條件對於系統貢獻度之方法。,最小平方法(the least square method),最小平方法（又稱最小二乘法）是一種數學優化技術，它通過最小化誤差的平方和找到一組數據的最佳函數匹配7 函數表示:,最小平方法(the least square method),利用最小平方法(線性組合形式)求解:其中，*為共軛轉置(associate),實例 for COMSOL
7、,由疊加法將表1之邊界條件展開為表2、3、4，其中表2為B1溫度點對於系統之貢獻度，求出敏感性係數、表3為B2溫度點對於系統之貢獻度，求出敏感性係數、表4為邊界條件對於系統之貢獻度，求出敏感性係數。其疊加的結果可以寫成式(A)。,(A),由於式(A)是三條方程式，但只有兩個未知數，因此將利用最小平方法去估算最佳的迴歸式求出近似解。,(A),最小平方法,COMSOL 驗證,參數:寬為20m長為20m統域邊界:steel AISI4340Mesh:正方形,20 x10(即每格為一平方公尺),COMSOL 驗證組(Direct),驗證組邊界設定:熱絕緣溫度:373.15K溫度:473.15K熱絕緣,
8、COMSOL 驗證組結果,P3:463.15KP4:463.15KP5:463.15K,COMSOL 實驗組(Inverse),左邊熱源結果P3:0.52307 KP4:0.073444 KP5:0.016213 K,COMSOL 實驗組(Inverse),COMSOL 實驗組(Inverse),左邊熱源結果P3:0.073444 KP4:0.49132 KP5:0.065743 K,COMSOL 實驗組(Inverse),COMSOL 實驗組(Inverse),結果P3:180.909078 KP4:195.931548 KP5:424.372615 K,COMSOL 實驗組(Inverse
9、),利用MATLAB計算,A=0.52307,0.073444;0.073444,0.49132;0.016213,0.065743;F=463.15;463.15;463.15;D=180.909078;195.931548;424.372615;x=inv(A*A)*A*(F-D)x=473.1505 473.1506=與驗證組 473.15幾乎樣(OK),COMSOL step by step,STEP1:選取模組開啟COMSOL(即自動到模型導覽視窗)點選熱傳模組點選廣義熱傳確定,Step2:繪製圖形,點選正方形圖在視窗上拉出一方形,3.在方形上按兩下，出現尺寸視窗，輸入參數，按確定4
10、.繪製完成，按 將圖形縮放至視窗大小,Step3:統御域設定,選取 物理量統御域設定出現統御域設定視窗輸入參數 STEEL AISI4340確定,Step4:邊界設定,選取 物理量邊界設定出現邊界設定視窗輸入邊界確定,Step5:設定網格,點選 網格 映射網格參數 出現映射網格參數視窗點選邊界設定網格數目 20 x10,Step6:進行求解,點選求解進行求解後處理資料顯示統御域(為了知道P3、P4、P5溫度),Step7:MATLAB,得到驗證組與對照組之數據後，利用MATLAB運算,Reference,Yann-Shouh Sun,Cheng-I Weng,Tei-Chen CHEN,and
11、 Wang-Long LI,Estimation of Surface Absorptivity and Surface Temperature,JAPANESE JOURNAL OF APPLIED PHYSICS PART 1-REGULAR PAPERS SHORT NOTES&REVIEW PAPERS 35(6A):3658-3664 JUN 1996孫彥碩，雷射表面熱處理的表面溫度和表面塗層吸收率的預測，國立成功大學機械工程研究所碩士論文，1995Tai-Sheng Chen*,Horng-Yuan Jang*,Tsung-Chien Chen*,and Pan-Chio Tuan
12、*,Adaptive Input Estimation Method Combining Finite Element Scheme for Inverse Jet Flow Heat Flux Estimating,JOURNAL OF C.C.I.T.,VOL.34,NO.2,May.,2006顏義哲，反算法於三維穩態鰭管式熱交換器鰭片熱傳系數之預測，國立成功大學造船季船舶機械工程研究所碩士論文，2002李弘毅與王景輝，逆向熱傳導問題之分析與實驗，力學系列B，第十六卷，第期，2000Introduction to Heat Transfer,THIRD EDITION,WILEYWikipadia(http:/zh.wikipedia.org/wiki/%E7%B9%81%E4%BD%93%E4%B8%AD%E6%96%87),
