本研究利用ANSYS Icepak針對板式熱交換器進行熱交換流場分析，
比較K200板式熱交換器、導流型消音器流道、圓錐形消音器流道三者所產生的紊流強弱多寡層度對熱交換效率之影響。
   首先必須證明紊流強烈混合作用所產生的質量傳遞是否為影響熱交換效率的主要因素，分別列舉層流、紊流兩種流場相比較，結果顯示紊流動能、黏度比較高之流道擁有較佳之熱交換效率，接下來的研究方向為如何改善流道幾何形狀以增加紊流動能、黏度比。
   研究過程中發現單純改變流道幾何形狀無法改善紊流動能耗損，最後針對K200板式熱交換器、導流型消音器流道、圓錐形消音器流道在強制對流下之散熱效率進行比較，研究指出導流型消音器流道所產生的紊流動能、黏度比優於其他兩者。

This study used ANSYS Icepak for heat exchange flow field analysis of plate heat exchangers.Compare the influence of the turbulence intensity and the degree of turbulence generated by the K200 plate heat exchanger, the flow-guiding muffler flow path and the conical muffler flow path on the heat exchange efficiency.
    Firstly, it must be proved whether the mass transfer caused by the intense mixing of turbulent flow is the main factor affecting the heat exchange efficiency. The flow fields of laminar flow and turbulent flow are compared respectively. The results show that the flow path with higher turbulent flow energy and higher viscosity has For better heat exchange efficiency, the next research direction is how to improve the flow channel geometry to increase the turbulent flow energy and viscosity ratio.
    During the research, it was found that simply changing the geometry of the flow channel could not improve the turbulent flow energy loss. Finally, the heat dissipation efficiency of the K200 plate heat exchanger, the flow-type muffler flow channel and the conical muffler flow channel under forced convection was compared. It is pointed out that the turbulent flow energy and viscosity ratio produced by the flow-through silencer flow passage are better than the other two.

目錄
中文摘要..........................................................................................................I
英文摘要.........................................................................................................II
目錄................................................................................................................III
表目錄..........................................................................................................VII
圖目錄.........................................................................................................VIII
符號表.......................................................................................................XVII
第一章 緒論...................................................................................................1
1.1文獻回顧...................................................................................................1
1.1.1消音器................................................................................................1
1.1.2板式熱交換器....................................................................................2
1.1.3紊流熱交換........................................................................................4
1.1.4 分析工具軟體相關文獻...................................................................6
1.1.5溫度對產品壽命的影響....................................................................7
1.2研究動機...................................................................................................7
第二章 基礎理論...........................................................................................9
2.1消音器原理...............................................................................................9
2.2熱交換原理...............................................................................................9
2.2.1渦流擴散............................................................................................9
2.2.2紊流普朗特數..................................................................................10
2.2.3紊流施密特數..................................................................................11
2.3板式熱交換器原理.................................................................................11
2.3.1熱交換器之流體流動型態..............................................................11
2.3.2對數平均溫差..................................................................................12
第三章 研究方法.........................................................................................13
3.1 ANSYS Icepak軟體介紹.......................................................................13
3.1.1標準渦流黏度模型的缺點..............................................................13
3.1.2可實現k-ε模型...............................................................................13
3.1.3強化雙方程模型..............................................................................15
3.1.4強化壁面處理..................................................................................15
3.2紊流相關物理量.....................................................................................16
3.2.1紊流動能..........................................................................................16
3.2.2紊流耗散率......................................................................................16
3.2.3黏度比..............................................................................................17
3.3模型建構.................................................................................................17
3.3.1繪製流程..........................................................................................17
3.3.2模型格式轉換..................................................................................17
3.4網格建立.................................................................................................17
3.4.1網格生成方法..................................................................................17
3.4.2網格品質分析..................................................................................18
第四章 結果分析與討論.............................................................................20
4.1初步分析.................................................................................................20
4.1.1網格劃分標準..................................................................................20
4.1.2控制組_實驗組_分析結果..............................................................20
4.2排除消音器流道不利於散熱的因素.....................................................21
4.2.1控制組_實驗組_分析結果..............................................................21
4.3消音器流道分析.....................................................................................22
4.3.1控制組_實驗組_分析結果..............................................................22
4.4紊流動能耗損對熱交換效率的影響.....................................................23
4.4.1控制組_分析結果............................................................................23
第五章 實際運用與比較.............................................................................24
5.1圓錐型散熱流道與導流型散熱流道.....................................................24
5.1.1控制組_實驗組_分析結果..............................................................24
5.2板式熱交換器的實際應用.....................................................................25
5.2.1控制組_實驗組_分析結果..............................................................26
第六章 結論.................................................................................................27
參考文獻.......................................................................................................28
表目錄
表3.4.2 .........................................................................................................33
表4.1.............................................................................................................33
表4.1.1..........................................................................................................34
表4.1.2 .........................................................................................................36
表4.2.............................................................................................................36
表4.2.1..........................................................................................................36
表4.2.2 .........................................................................................................38
表4.2.3..........................................................................................................38
表4.2.4..........................................................................................................38
表4.2.5..........................................................................................................40
表5.1.............................................................................................................40
表5.1.1..........................................................................................................40
表5.1.2..........................................................................................................42
表5.2.............................................................................................................43
表5.2.1..........................................................................................................44
表5.2.2..........................................................................................................45
圖目錄
圖2-1熱邊界層厚度與速度邊界層.............................................................46
圖2-2熱邊界層厚度.....................................................................................46
圖2-3板式熱交換器組立圖.........................................................................47
圖2-4人字形波紋熱傳板.............................................................................47
圖2-5板式熱交換器內部流體運作過程.....................................................48
圖2-6熱交換流體流動型態.........................................................................48
圖2-7對數平均溫差(LMTD) ......................................................................49
圖2-8圓錐型膨脹室消音器.........................................................................49
圖2-9導流型膨脹室消音器.........................................................................49
圖3-1 ANSYS Icepak有限體積法...............................................................50
圖3-2離散為有限的控制體或元素.............................................................50
圖3-3 O形網格.............................................................................................51
圖3-4結構化網格;非結構化網格................................................................51
圖3-5面對齊率.............................................................................................51
圖4-1控制組_尺寸圖...................................................................................52
圖4-2控制組_等角視圖...............................................................................52
圖4-3實驗組_尺寸圖...................................................................................53
圖4-4實驗組_等角視圖...............................................................................53
圖4-5控制組_導流板(總表面積;容積) ......................................................54
圖4-6實驗組_導流板(總表面積;容積) ......................................................54
圖4-7控制組初步分析_網格控制...............................................................54
圖4-8實驗組初步分析_網格控制...............................................................55
圖4-9控制組初步分析_熱源物件...............................................................55
圖4-10實驗組初步分析_熱源物件.............................................................55
圖4-11控制組初步分析_無流動,熱輻射狀態,整體溫度...........................56
圖4-12控制組初步分析_無流動,熱輻射狀態_熱源溫度..........................56
圖4-13控制組初步分析_流動,熱交換狀態,熱源溫度(剖面位置)............56
圖4-14控制組初步分析_流動,熱交換狀態,整體溫度分佈.......................56
圖4-15控制組初步分析_流動,熱交換狀態,整體溫度分佈.......................57
圖4-16控制組_流動,熱交換狀態,流場壓力...............................................57
圖4-17控制組初步分析_流動,熱交換狀態,紊流動能...............................57
圖4-18控制組初步分析_流動,熱交換狀態,黏度比...................................57
圖4-19實驗組初步分析_無流動,熱輻射狀態,整體溫度...........................58
圖4-20實驗組初步分析_無流動,熱輻射狀態,熱源溫度...........................58
圖4-21實驗組初步分析_流動,熱交換狀態,熱源溫度(剖面位置)............58
圖4-22實驗組初步分析_流動,熱交換狀態,整體溫度分佈.......................58
圖4-23實驗組初步分析_流動,熱交換狀態,流場流速...............................59
圖4-24實驗組初步分析_流動,熱交換狀態,流場壓力...............................59
圖4-25實驗組初步分析_流動,熱交換狀態,紊流動能...............................59
圖4-26實驗組初步分析_流動,熱交換狀態,黏度比...................................59
圖4-27控制組_物件尺寸圖,空間尺寸圖....................................................60
圖4-28控制組_等角視圖.............................................................................60
圖4-29實驗組_物件尺寸圖,空間尺寸圖....................................................61
圖4-30實驗組_等角視圖.............................................................................61
圖4-31控制組(總表面積;容積) ..................................................................62
圖4-32實驗組(總表面積;容積) ..................................................................62
圖4-33控制組_最大網格尺寸.....................................................................62
圖4-34實驗組_最大網格尺寸.....................................................................62
圖4-35控制組_流動,熱交換狀態,出口溫度...............................................63
圖4-36控制組_流動,熱交換狀態,整體溫度分佈.......................................63
圖4-37控制組_流動,熱交換狀態,流場流速...............................................63
圖4-38控制組_流動,熱交換狀態,流場壓力..............................................63
圖4-39控制組_流動,熱交換狀態,紊流動能..............................................64
圖4-40控制組_流動,熱交換狀態,黏度比..................................................64
圖4-41實驗組_流動,熱交換狀態,出口溫度..............................................64
圖4-42實驗組_流動,熱交換狀態,整體溫度分佈......................................64
圖4-43實驗組_流動,熱交換狀態,流場流速..............................................65
圖4-44實驗組_流動,熱交換狀態,流場壓力..............................................65
圖4-45實驗組_流動,熱交換狀態,紊流動能..............................................65
圖4-46實驗組_流動,熱交換狀態,黏度比..................................................65
圖4-47控制組_尺寸圖.................................................................................66
圖4-48控制組_等角視圖.............................................................................66
圖4-49實驗組_尺寸圖.................................................................................67
圖4-50實驗組_等角視圖.............................................................................67
圖4-51控制組(總表面積;容積) ..................................................................68
圖4-52實驗組(總表面積;容積) ..................................................................68
圖4-53控制組_最大網格尺寸.....................................................................69
圖4-54實驗組_最大網格尺寸.....................................................................69
圖4-55控制組_實驗組_熱源區域...............................................................70
圖4-56控制組_無流動,熱輻射狀態............................................................70
圖4-57控制組_無流動,熱輻射狀態,熱源溫度...........................................71
圖4-58控制組_流動,熱交換狀態,熱源溫度...............................................71
圖4-59控制組_流動,熱交換狀態,整體溫度分佈.......................................72
圖4-60控制組_流動,熱交換狀態,流場流速...............................................72
圖4-61控制組_流動,熱交換狀態,流場壓力...............................................73
圖4-62控制組_流動,熱交換狀態,紊流動能...............................................73
圖4-63控制組_流動,熱交換狀態,黏度比...................................................74
圖4-64實驗組_無流動,熱輻射狀態............................................................74
圖4-65實驗組_無流動,熱輻射狀態,熱源溫度...........................................75
圖4-66實驗組_流動,熱交換狀態,熱源溫度...............................................75
圖4-67驗組_流動,熱交換狀態,整體溫度分佈...........................................76
圖4-68實驗組_流動,熱交換狀態,流場流速...............................................76
圖4-69實驗組_流動,熱交換狀態,流場壓力...............................................77
圖4-70實驗組_流動,熱交換狀態,紊流動能...............................................77
圖4-71實驗組_流動,熱交換狀態,黏度比...................................................78
圖4-72控制組_流動,熱交換狀態,紊流停滯區域剖面位置.......................78
圖4-73控制組_流動,熱交換狀態,紊流停滯區域剖面,整體溫度分佈......79
圖4-74控制組_流動,熱交換狀態,紊流停滯區域剖面,流場流速..............79
圖4-75控制組_流動,熱交換狀態,紊流停滯區域剖面,流場壓力..............79
圖4-76控制組_流動,熱交換狀態,紊流停滯區域剖面,紊流動能..............80
圖4-77控制組_流動,熱交換狀態,紊流停滯區域剖面,黏度比..................80
圖5-1控制組_尺寸圖...................................................................................81
圖5-2控制組_等角視圖...............................................................................81
圖5-3實驗組_尺寸圖...................................................................................82
圖5-4實驗組_等角視圖...............................................................................82
圖5-5控制組(總表面積;容積) ....................................................................83
圖5-6實驗組(總表面積;容積) ....................................................................83
圖5-7控制組_最大網格尺寸.......................................................................84
圖5-8實驗組_最大網格尺寸.......................................................................84
圖5-9控制組_實驗組_熱源區域.................................................................85
圖5-10控制組_無流動,熱輻射狀態............................................................85
圖5-11控制組_無流動,熱輻射狀態,熱源溫度...........................................86
圖5-12控制組_流動,熱交換狀態,熱源溫度...............................................86
圖5-13控制組_流動,熱交換狀態,整體溫度分佈.......................................87
圖5-14實驗組_流動,熱交換狀態,流場流速...............................................87
圖5-15控制組_流動,熱交換狀態,流場壓力...............................................88
圖5-16控制組_流動,熱交換狀態,紊流動能...............................................88
圖5-17控制組_流動,熱交換狀態,黏度比...................................................89
圖5-18控制組_流動,熱交換狀態,區域剖面位置.......................................89
圖5-19控制組_流動,熱交換狀態,區域剖面,整體溫度分佈......................90
圖5-20控制組_流動,熱交換狀態,區域剖面,流場流速..............................90
圖5-21控制組_流動,熱交換狀態,區域剖面,流場壓力..............................90
圖5-22控制組_流動,熱交換狀態,區域剖面,紊流動能..............................91
圖5-23控制組_流動,熱交換狀態,區域剖面,黏度比..................................91
圖5-24控制組_流動,熱交換狀態,導流缺口設計,流場軌跡線..................92
圖5-25實驗組_無流動,熱輻射狀態............................................................92
圖5-26實驗組_無流動,熱源溫度................................................................92
圖5-27實驗組_流動,熱交換狀態,熱源溫度...............................................93
圖5-28實驗組_流動,熱交換狀態,整體溫度分佈.......................................93
圖5-29實驗組_流動,熱交換狀態,流場流速...............................................94
圖5-30實驗組_流動,熱交換狀態,流場壓力...............................................94
圖5-31實驗組_流動,熱交換狀態,紊流動能...............................................95
圖5-32實驗組_流動,熱交換狀態,黏度比...................................................95
圖5-33實驗組_流動,熱交換狀態,流場軌跡線...........................................96
圖5-34實驗組_紊流分佈區域.....................................................................96
圖5-35控制組_尺寸圖.................................................................................97
圖5-36控制組_等角視圖.............................................................................97
圖5-37實驗組_尺寸圖.................................................................................97
圖5-38實驗組_等角視圖.............................................................................98
圖5-39控制組(流道容積,熱傳導金屬面積)..............................................98
圖5-40實驗組(流道容積,熱傳導金屬面積)..............................................99
圖5-41控制組_最大網格尺寸.....................................................................99
圖5-42實驗組_最大網格尺寸...................................................................100
圖5-43控制組,實驗組_人字形波紋重疊位置示意圖..............................100
圖5-44控制組_流動,熱交換狀態,整體溫度分佈(由上至下A,B)...........101
圖5-45控制組_流動,熱交換狀態,板片溫度分佈(由上至下A,B,C).......101
圖5-46控制組_流動,熱交換狀態,流場流速.............................................102
圖5-47控制組_流動,熱交換狀態,流場壓力.............................................102
圖5-48控制組_流動,熱交換狀態,紊流動能.............................................102
圖5-49控制組_流動,熱交換狀態,黏度比.................................................102
圖5-50控制組_流動,熱交換狀態,熱源層出口溫度.................................103
圖5-51控制組_流動,熱交換狀態,冷卻層出口溫度.................................103
圖5-52實驗組_流動,熱交換狀態,整體溫度分佈(由上至下A,B) ..........103
圖5-53實驗組_流動,熱交換狀態,板片溫度分佈(由上至下A,B,C) ......104
圖5-54實驗組_流動,熱交換狀態,流場流速.............................................104
圖5-55實驗組_流動,熱交換狀態,流場壓力(由上至下A,B,C)...............105
圖5-56實驗組_流動,熱交換狀態,紊流動能.............................................105
圖5-57實驗組_流動,熱交換狀態,黏度比.................................................106
圖5-58實驗組_流動,熱交換狀態,熱源層出口溫度.................................106
圖5-59實驗組_流動,熱交換狀態,冷卻層出口溫度.................................106
圖5-60實驗組_紊流分佈區域...................................................................106

[1] Bernhard, Robert J., Shape optimization of reactive mufflers,Noise Control Engineering Journal,27, 1, 10-17, 1986.

[2]Davis D. D.,Stokes J. M.,Moore D. and Steven L.,Theoretical and Experimental Investigation of Mufflers With Comments on Engine Exhaust Muffler Design, NACA Report 1192,1954.

[3]Davis D. D.,Jr. G. L. Steven,Moore Jr. D. and Stokes G. M., Theoretical and Measured Attenuation of Mufflers at Room Temperature without Flow, with Comments on Engine-Exhaust Muffler Design,NACA TN-2893, 1954.

[4]Gogate G.R. and M.L.Munjal,“Analytical and experimental aeroacoustic studies of open-ended three-duct perforated elements used in muffler,”Journal of Acoustic Society America, Vol. 97, No.5, 1995, pp. 2919–2927.

[5]Brand J. F.,Garcia P. and Wiemeler D.,“Surface Radiated Noise of Exhaust Systems–Calculation and Optimization with CAE,”Noise & Vibration Conference and Exhibition,Michigan,SAE Paper 2004-01-0407, 2004.

[6] Munjal M.L., Acoustics in Ducts and Mufflers, John Wiley and Sons, New York, USA, 1987.

[7] Fujita, K., Abe, T., Hori, Y., Simulation of Acceleration Pass-by Noise Considering the Acoustic Radiation Characteristics of a Vehicle Body, JSAE Review, Vol. 7, No. 3, pp. 52-58, 1986.

[8] Takayama,K. et al., Shock wave development and propagation in 
automobile exhaust system ,SAENO. 880082,1989. 

[9] Li Peizi,Dai Genhua,Zhu Zhichi,Noise Radiation of a Strongly Pulsation 
Tailpipe Exhaust, Journal of sound and vibration 167, 1993,pp.385-400. 

[10]Chandrashekhar Bhat,Sharma S.S.,Jagannath K,N S Mohao,Sathisha S G, “Design and Analysis of a Expansion Chamber Mufflers,” World Journal of Engineering, pp 117-118.

[11]王紅苗,人字形波紋槽道內傳熱,流阻和流型的實驗研究[D],北京：清華大學, 1989.

[12]許淑惠,周明連,板式熱交換器內壓力分佈、流阻及流型顯示實驗研究[J] ,節能,1996.

[13]Davis E. J., Exact solution for a class of heat and mass transfer problem,
Can. J. Chem. Eng.,1973，51, 562-572.

[14]Fernando V. Castellões, João N.N. Quaresma,Convective heat transfer
enhancement in low Reynolds number flows with wavy walls,Renato M. Cotta,International Journal of Heat and Mass Transfer 53 (2010)2022–2034.

[15]Marcos Vera a, Amable Liñán, Laminar counterflow parallel-plate heat
exchangers: Exact and approximate solutions, International Journal of
Heat and Mass Transfer 53 (2010) 4885–4898.

[16]Alberto E. Quintero, Marcos Vera , Bernardo Rivero-de-Aguilar, Wall
conduction effects in laminar counterflow parallel-plate heat exchangers, International Journal of Heat and Mass Transfer 70 (2014) 939–953.

[17]Gao S. X. and Hartnett J. P., Heat transfer behavior of reiner-rivlin
fluids in rectangular ducts,lnt.J.Heat Mass Transfer,1996,39,1317-1324.

[18]Lee D. L.,Thomas F. and Irvine Jr., Shear rate dependent thermal
conductivity measurements of non-newtonian fluids,Experimental
Thermal and Fluid Science,1997,15,16-24.
[19]Koh Y.H.,Ong N.S.,Chen X.Y.,Lam Y.C. and Chai J.C., Effect of
temperature and inlet velocity on the flow of a nonnewtonian fluid, Int.
Comm. Heat Mass Transfer,2004, 31,1005-1013.

[20]Harris C.,Despa M.,Kelly K.,“Design and Fabrication of a Cross Flow Micro Heat Exchanger”,Journal of Microelectromechanical Systems,Vol.9,No. 4, pp.502-508 (2000).

[21]Harris C.,Kelly K.,Wang T.,Andrew McCandless and Shariar Motakef,“Fabrication, Modeling, and Testing of Micro-Cross-Flow Heat Exchangers”,Journal of Microelectromechanical Systems,Vol.11,No. 6, pp.726-735(2002).

[22]Schubert K.,Brander J.,Fichtner M.,Linder G.,Schygulla U.,Wenka A.,
“Microstructure Devices for Applications in Thermal and Chemical Process Engineering”,Microscale Thermophysical Engineering,Vol.5,pp.17-39(2001).

[23] Brander J., Anurjew E.,Bohm L.,Hansjosten E.,Henning T.,Schygulla U.
,Wenka A.,Schubert K.,“Concepts and Realization of Microstructure Heat Exchangers for Enhanced Heat Transfer”, Experimental Thermal and Fluid Science, Vol. 30, pp. 801-809 (2006).

[24] Tunc G., Bayazitoglu Y., “Heat Transfer in Rectangular Microchannels,”International Communications in Heat and Mass Transfer,Vol.45,pp.765-773,(2002).

[25]Muley, A., and Manglik, R. M. (1999) “Experimental study of tubulent flow heat transfer and pressure drop in a plate heat exchanger with chevron plates,” J. Heat Transfer, 121:110-117.

[26] Baker E., “Liquid immersion cooling of small electronic devices”,
SDOL Microelectronic and Reliability,Vol.12, pp. 163-173,1973.

[27] 李嘉慶,"衝擊流中渦流產生器對熱傳及流場特性探討",淡江大學航空太空工程學系,碩士論文,2003.

[28] Mori Y., Uchida Y. and Ukon T.,‘‘Forced convective heat transfer in
a curved channel with a square cross section’’, Int. J. Heat MassTransfer 14, pp1787-1805, 1971.

[29]Gül H.,Evin D.,‘‘Heat transfer enhancement in circular tubes using helical swirl generator insert at the entrance’’,International Journal of Thermal Sciences 46, 1297–1303, 2007.

[30] F. J. Edwards and C.J.R.Alker, The Improvement of Forced Convection 
Surface Heat Transfer Using Surface Protrusions in the Form of(A)Cubes and 
(B)Vortex Generators,5th, Heat Transfer Conf.,Tokyo,vol.2,pp.244-248, 1974.

[31]C.M.B.Russel,T.V.Jones,and G.H.Lee,Heat Transfer En-hancement Using Vortex Generators,Seventh Int.Heat Transfer Conf.,New York,vol. 2, pp. 283-288, 1982.

[32]林岳宏"水對空氣在板式熱交換器之性能測試分析", 國立中央大學機械工程研究所碩士論文,2003.

[33]陳任元“板式熱交換器之板片數對其性能影響之分析”,私立淡江大學機械工程研究所碩士論文,1995.

[34]羅弘達“板式熱交換器之入出口壓降實驗分析”,國立中央大學機械工程研究所碩士論文,2005.

[35]Kays,W.M.(1994)."Turbulent Prandtl Number—Where Are We?".Journal of Heat Transfer.116(2):284–295.doi:10.1115/1.2911398.

[36] Nelson, R. M. and Pletcher, R. H., “An Explicit Scheme for the Calculation of Confined Turbulent Flow with Heat Transfer”, Proc. Heat Transfer and Fluid Mechanics Institute, Standford University Press, pp. 154-170, 1974.

[37] Bret A. Zahn, “Evaluating thermal characterization accuracy using
CFD codes-A Package Level Benchmark Study of IcePak and Flotherm”, IEEE , 1998.

[38]陶文銓.數值傳熱學.第2 版.西安交通大學出版社,2001,374-376.

[39]吳海玲，陳聽寬,羅毓珊.應用不同紊流模型的二維橫向射流傳熱數值模擬研究.西安交通大學學報,2001,35(9):903 - 907.

[40]王福軍.計算流體動力學分析2CFD軟件原理與應用.清華大學出版社,2004:124-126.

[41] TAO Ping Kee, LEE Hong Moon. New Physics at Work - HEAT. Oxford 
University Press. 1993. ISBN 0-19.

[42]周一塵,次世代面板框膠之可靠度預估及加速壽命試驗之研究,碩士論文,國立成功大學工程科學系,台南,2009. 

[43]周國鵬,加速壽命試驗於無線滑鼠平均壽命之研究,碩士論文,華梵大學工業工程與經營資訊學系,新北, 2009.