食品热传递工程之数值模拟.pptVIP

食品熱傳遞工程之數值模擬進度報告 指導老師：楊懷文 計畫學生：鄭湘穎 萬鴻鋼 報告日期：2010/09/17 已完成項目 閱讀相關食品熱傳遞教材 一維空間非穩定狀態 時間-從屬性計算 MatLab 程式設計實習 二維空間穩定狀態數值模擬 閱讀相關食品熱傳遞教材 熱傳導(Thermal conduction) 熱對流(Thermal convection) 非穩定狀態熱傳遞 有限差異數值分析 一維空間非穩定狀態 時間-從屬性計算 整塊容積法(Lumped Capacitance Method) 假設被加熱的物體內部無溫度梯度(差異) 熱傳導度趨近無限大 無法以傅立葉方程式評估時間重數關係 需回歸至最基礎的能量平衡方程式 整塊容積法有效性驗證 特質長度(characteristic length. Lc) 熱傳導係數(k) 表面熱傳系數(h) Lc* k/ h= 比爾數(Biot number) 當比爾數小於0.1時則整塊容積法之誤差較小，因而成立並且適用。 實例說明 以熱電偶是溫度探針為例:其探針尖點的構造可是為一圓球狀假設此溫度探針符合塊狀容積條件，可做下列二項時間從屬性質評估 設定時間常數，以評估適當的探針尖點直徑 設定測量溫度誤差允許範圍，探討所需之測量時間 非整塊容積之處理原則 食品由於熱傳導度低，加熱時內部會產生溫度梯度，整塊容積不適用，需考慮起始條件與邊界條件並以下列其中一種方法 以數學解析法評估 需以 transental 方程式之無限項次正數根(positive roots)做為確切解析 由無限項次正數根之前4~10項做近似解析 以數值分析法評估 有限差異法 MatLab 程式設計實習 處理集合數字 應用於工程問題 MatLab 程式設計實習 = ? 指定運算子 X=3 ≠3=X X=X+3 X=Y+3 MatLab 程式設計實習 * ? 乘 c=a*b / ? 除 c=a/b ; ? 不顯示運算結果 % ? 註解，指定變數資料格式 MatLab 程式設計實習 [ ]?宣告陣列元素 ` ?矩陣倒置 MatLab 程式設計實習 矩陣加減 須具有相同維度 Ex: A=[12 34 56 20]; B=[1 3 2 4]; C=A+B ? C=[13 37 58 24] MatLab 程式設計實習 矩陣乘除 Ex: A=[ 4 4 ]; B=[ 2 2 2 ]; C=A*B ? C=[ 8 8 8 8 8 8 ] 將所學習的 Matlab 相關功能套用在練習題目中，進行數值與圖形分布的模擬分析，以便於日後實際熱傳遞係數測定與液態熱對流實習。 一維空間穩定狀態數值模擬模型陣列(A) A B 0 0 0 0 0 0 … T1 E1 0 B A B 0 0 0 0 … T2 E2 0 0 B A B 0 0 0 … T3 E2 …………………… X … = … …………………… … … 0 0 0 0 0 B A B 0 Tn-1 E2 …..................B C Tn E3 一維空間非穩定狀態數值模擬模型 定義陣列=[………..] 模型陣列(A) ? MatLab 結果陣列=[………..] 結果陣列=定義陣列 結果陣列=歸零陣列 模擬結果 二維空間穩定狀態數值模擬 後續實驗規劃 * * ？ Ex: y = (5*2+0.08)/5; % 將運算結果儲存在 y，但不用顯示於螢幕 z = y*2 % 將運算結果儲存在 z，並顯示於螢幕 ? z = 4.032 Ex: [ 1 2 3 A=[ 1 2 3;4 5 6;7 8 9 ] ? A= 4 5 6