本研究透過數值分析之方式，模擬以理想氣體為工作流體之gamma型態史特靈引擎內部震盪流場與整體溫度分布情形，係利用計算流體動力學分析軟體ESI-CFD來執行運算；於本研究中，將以2-D之數學模型進行整體gamma型態史特靈引擎結構之建模計算，並針對單、雙動力汽缸做比較。時間與體積之變化將表現出對應之暫態的壓力變化及溫度場之分布；模擬之溫度場變化與實際實驗結果呈現出相當好的比對。最後，利用porous-media模組進行初步之復熱再生器的模擬研究。

This research study a numerical analysis on the oscillatory flow and temperature distribution inside a gamma type Stirling engine with ideal gas as working fluid. The CFD code used in the investigation for modeling is ESI-CFD. In this study, a 2-D mathematical model is formulated for a full scale gamma-configuration Stirling engine with single and dual power pistons. Variation of time, piston volume with pressure and temperature are presented. The calculated temperature distribution results compare well with experimental results from the other research in our group. Also, a preliminary study for an initial thermal equilibrium porous-media model used in CFD codes for simulation of Stirling regenerators is presented.

目  錄
誌謝 I
中文摘要 II
英文摘要 III
目錄 IV
圖目錄 VI
表目錄 VIII
符號說明 IX
第一章 緒論 01
1.1 研究背景與動機 01
1.2 史特靈引擎簡介 03
1.2.1 作動原理 03
1.2.2 引擎形式 03
1.3 文獻回顧 07
1.4 研究目的 12
第二章 理論分析 13
2.1 模擬分析軟體簡介 13
2.2 基本假設 15
2.3 統馭方程式（Governing Equations） 15
2.3.1 質量守恆方程式 16
2.3.2 動量守恆方程式 16
   2.3.3 能量守恆方程式 18
   2.3.4 紊流傳遞方程式 18
2.4 數值方法 19
2.4.1 有限體積法 19
2.4.2 內插法 23
2.4.3 有限差分方程式 24
2.5 邊界條件 25
第三章 數值求解 27
3.1 SIMPLEC演算法 27
3.2 低鬆弛係數（Under-Relaxation） 30
    3.2.1 虛擬時間鬆弛（false time step relaxation） 30
    3.2.2 線性鬆弛（linear relaxation） 31
3.3 建模及網格化設置 31
3.4 參數分析 35
3.5 求解流程 39
第四章 結果與討論 41
4.1 模擬驗證結果 41
4.2 再生器的設計 51
4.3 參數設定與結果 51
第五章 結論與未來建議 53
5.1 結論 55
5.2 未來建議 56
參考文獻 57
圖目錄
圖1-1 史特靈引擎構造與作動方式 06
圖1-2 史特靈引擎類型 06
圖1-3 高階矩陣之階數越高數值收斂越快 10
圖1-4 Free piston之2D建模 11
圖1-5 無重力情況下的流線 11
圖1-6 有重力情況下的流線 11
圖2-1 三維計算單元（控制體積） 20
圖2-2 二維計算單元（控制體積） 20
圖2-3 計算邊界單元 26
圖3-1 史特靈引擎尺寸圖（單缸） 33
圖3-2 2D建模示意圖（單缸） 33
圖3-3 史特靈引擎尺寸圖（雙缸） 33
圖3-4 2D建模示意圖（雙缸） 33
圖3-5 史特靈引擎設計組合圖a 34
圖3-6 史特靈引擎設計組合圖b 34
圖3-7 結構性網格架構 34
圖3-8 非結構性網格架構 35
圖3-9 各邊界示意圖 37
圖3-10 連桿機構示意圖 38
圖3-11 活塞與移氣器位移軌跡圖 38
圖3-12 CFD-ACE之求解流程圖 40
圖4-1 汽缸內壓力穩態示意圖（單缸） 42
圖4-2 出水口溫度變化圖（單缸） 43
圖4-3 動力活塞於上死點整體溫度分布圖（單缸） 43
圖4-4 動力活塞於下死點整體溫度分布圖（單缸） 44
圖4-5 實驗與模擬出水口溫度變化圖（單缸） 44
圖4-6 汽缸內壓力穩態示意圖（雙缸） 45
圖4-7 出水口溫度變化圖（雙缸） 45
圖4-8 實驗與模擬出水口溫度變化圖（雙缸） 46
圖4-9 動力活塞於上死點整體溫度分布圖（雙缸） 46
圖4-10 動力活塞於下死點整體溫度分布圖（雙缸） 47
圖4-11 穩態前五個週期之P-V圖（單缸） 48
圖4-12 穩態前五個週期之P-V圖（雙缸） 49
圖4-13 模擬與實驗P-V圖（單缸） 50
圖4-14 模擬與實驗P-V圖（雙缸） 50
圖4-15 模擬單缸與雙缸P-V圖 51
圖4-16 復熱再生器位置示意圖 52
圖4-17 孔隙率70﹪之P-V圖 53
圖4-18 孔隙率75﹪之P-V圖 53
圖4-19 孔隙率80﹪之P-V圖 54
表目錄
表3-1 模擬計算參數設定 37
表4-1 復熱再生器設計參數 52

參考文獻
[1]	士福金，「漫談史特靈機械」，船舶科技期刊  第14期  第38-48頁。
[2]	林貴添、林祥輝、簡國祥，2003，「史特靈引擎的性能數值分析」，能源季刊  第三十三卷  第三期　第78-93頁。
[3]	林貴添、簡國祥，2003，「史特靈引擎之復熱再生器的設計規劃分析」，能源季刊　第三十三卷　第一期　第114-125頁。
[4]	施長江，2004，「史特靈引擎菱形驅動結構之機構設計與熱流分析」，大同大學碩士論文。
[5]	劉賽、李青、李正宇、李強，2006，「微型整体式史特靈致冷機系統的計算機仿真優化研究」，低溫工程期刊 第三期 第43-47頁。
[6]	Ibrahim, M.B., Tew, R.C., Zhang, Z., Gedeon, D., and Simon, T.W., “CFD Modeling of Free-Piston Stirling Engines,” Proceedings of the 36th Intersociety Energy Conversion Engineering Conference, Savannah,GA,2001(IECEC2001-CT-38), NASA/TM –2001-211132, 2001.
[7]	Kornhauser, A.A., “Gas-Wall Heat Transfer During Compression and Expansion, ” Massachusetts Institute of Technology, S.D.Thesis, 1989.
[8]	Dyson, R. W., “Technique for Very-High-Order-Nonlinear Simulation and Validation,” J. Comp. Acoustics, Vol. 10, No. 2, pp.211-229, 2002.
[9]	Dyson, R. W., Wilson, S. D., and Tew, R. C., “Review of Computational Stirling Analysis Methods,” in proceedings of 2nd International Energy Conversion Engineering Conference, Paper No. AIAA-2004-5582, Providence, RI, August 15-19, 2004a.
[10]	Willems, D., “High-Power Cryocooling,”2007
[11]	Tew, R.C., Simon, T.W.,Gedeon, D., Ibrahim, M.B., “An Initial Non-Equilibrium Porous-Media Model for CFD Simulation of Stirling Regenerators.” NASA/TM—2006-214391.
[12]	Ibrahim, M.B., Zhang, Z., Kembhavi, S., Simon, T.W., and Gedeon, D., “A 2-D Axisymmetric CFD Model of Oscillatory Flow with Separation,” Proceedings of the 37th Intersociety Energy Conversion Engineering Conference, Washington, D.C.,IECEC 2002 Paper No. 20121, (2002a).
[13]	Dyson, R. W., Wilson, S. D., and Tew, R. C., “Review of Computational Stirling Analysis Methods,” in proceedings of 2nd International Energy Conversion Engineering Conference, Paper No. AIAA-2004-5582, Providence, RI, August 15-19, 2004a.
[14]	Tew, R.C., “Study of Two-Dimensional Compressible Non-Acoustic Modeling of Stirling Machine Type Components,”Proceedings of the 36th Intersociety Energy Conversion Engineering Conference, Savannah,GA,2001 (IECEC2001-CT-27),NASA/TM –2001-211066, 2001.
[15]	Ebiana, A., Savadekar, R., and Vallury, A., “2nd Law Analysis of Sage and CFD-ACE+ Models of MIT Gas Spring and “Two-Space” Test Rigs,” in proceedings of 2nd International Energy Conversion Engineering Conference, Paper No. AIAA-2004-5778, Providence, RI, August 15-19, 2004.
[16]	Tew, R.C., Cairelli, J.E., “Overview of NASA Multi-Dimensional Stirling Convertor Code Development and Validation Effort. ” NASA/TM—2002-211997.
[17]	Tew, R.C., Dyson, R.W., “Overview 2004 of NASA-Stirling Convertor CFD Model Development and Regenerator R&D Efforts. ” NASA/TM—2004-213404.
[18]	Versteeg, H.K., Malalasekera, W.,“An introduction to computational fluid dynamics The finite volume method. ”1995.
[19]	CFD Research Corp., “CFDRC User Manual, ”2004.