超音速燃燒衝壓引擎的工作流場中有些地方較為稀薄，因此，巨觀的連續體假設將不適用，需改用微觀的粒子觀點來分析，本文採用數值解(Numerical Solution)內的其中一種方法，直接模擬蒙地卡羅法(DSMC)，此法是目前較精確且廣泛用於模擬分析稀薄流的方法。
    本文以DSMC法模擬超音速燃燒衝壓引擎內流場的變化，考慮燃料噴入與空氣混合情形的狀況，此時，如何達到良好的混合效率便成了主要關注研究的方向，因此，針對各種不同構型之燃燒室計算分析其混合效率。

Supersonic combustion ramjet engine work flow field are thin. Therefore, We have to use micro particle to analyze. The macro continuum hypothesis is no longer suitable. This article use one of the Numerical Solution to simulate DSMC. It is precise and used a lot of analyze recently. 
    This article adopted the variety of Scramjet inner flow field by DSMC. We consider the situation of fuel ignition and air mixing. How to achieve good mixing efficiency has become concerned about the direction of the research, therefore, for a variety of different the combustion chamber configuration calculation and analysis its mixing efficiency.

目錄
誌謝 I
目錄 IV
表目錄 VII
圖目錄 VIII
符號表 XII
第一章	緒  論 1
1-1	前言 1
1-2	文獻回顧 2
1-3	紐森數的定義 6
1-4	稀薄氣體的定義( Dilute gas) 8
1-5	波茲曼方程式( Boltzmann equation) 9
第二章	直接模擬蒙地卡羅法 11
2-1	直接模擬蒙地卡羅法介紹( DSMC Introduction) 11
2-2	網格設置與時步的計算 13
2-3	流場初始條件 13
2-4	流場邊界條件 14
2-5	氣體分子模型 17
2-5-1	硬球模型(Hard sphere ,HS) 17
2-5-2	可變軟球模型(Variable soft sphere model ,VSS) 18
2-5-3	可變硬球模型(Variable hard sphere model ,VHS) 18
2-5-4	雙原子分子模型 19
2-6	碰撞 22
第三章	超音速燃燒衝壓引擎 24
3-1	超音速燃燒衝壓引擎概述 24
3-2	超音速燃燒衝壓引擎構造 25
3-3	超音速燃燒衝壓引擎的挑戰與技術 26
第四章	模擬結果與討論 28
4-1	程式驗證 28
4-2	模擬外型 29
4-3	混合效率 30
4-4	全壓損失 30
4-5	相同外型不同馬赫數比較 31
4-5-1	自由流12、15、18馬赫算例外流場變化 31
4-5-2	自由流12、15、18馬赫算例內流場變化 32
4-6	燃燒室內部結構不同比較 34
4-6-1	自由流12馬赫算例外流場變化 34
4-6-2	自由流15、18馬赫算例外流場變化 35
4-6-3	自由流12、15、18馬赫內流場變化 36
4-7	迴流區與混合效率 38
4-8	自由流12、15、18全壓損失比 40
第五章	結論與建議 41
5-1	結論 41
5-2	未來建議與工作 42
參考文獻 43

表目錄
表4-1 自由流條件設置 47
表4-2 可變硬球模型基本參數(273K，1atm) 47

圖目錄
圖1-1 X-51進行無人駕駛飛機圖 48
圖1-2 kn值與統御方程式間的關係圖 49
圖2-1 稀薄度與密度之關係圖	49
圖2-2 直接模擬蒙地卡羅法流程圖 50
圖2-3 粒子與壁面交互做用(鏡面反射邊界) 51
圖2-4 粒子與壁面交互作用(漫反射邊界)	51
圖2-5 分子碰撞面積示意圖 51
圖3-1 渦輪引擎構造示意圖 52
圖3-2 衝壓引擎構造示意圖 52
圖3-3 超音速燃燒衝壓引擎構造圖 52
圖3-4 噴射推進系統圖 53
圖4-1 12馬赫算例與實驗數據比較圖(A)整流罩(B)中心體 54
圖4-2 15馬赫算例與實驗數據比較圖(A)整流罩(B)中心體 54
圖4-318馬赫算例與實驗數據比較圖(A)整流罩(B)中心體 54
圖4-4 12馬赫算例之馬赫數分佈比較圖 55
圖4-5 15馬赫算例之馬赫數分佈比較圖 55
圖4-6 18馬赫算例之馬赫數分佈比較圖 56
圖4-7 二維超音速燃燒衝壓引擎外型(no step) 56
圖4-8 二維超音速燃燒衝壓引擎外型(convex) 57
圖4-9 二維超音速燃燒衝壓引擎外型(cavity) 57
圖4-10二維超音速燃燒衝壓引擎外型(cavity slop) 57
圖4-11 三維超音速燃燒衝壓引擎外型 58
圖4-12 12馬赫算例外流場馬赫數分佈圖 59
圖4-13 12馬赫算例外流場壓力線分佈圖 60
圖4-14 15馬赫算例外流場馬赫數分佈圖 61
圖4-15 15馬赫算例外流場壓力分佈圖 62
圖4-16 18馬赫算例外流場馬赫數分佈圖 63
圖4-17 18馬赫算例外流場壓力線分佈圖 64
圖4-18 內流場馬赫數算例分佈圖(12MA、15MA、18MA) 65
圖4-19 內流場壓力分佈算例圖(12MA、15MA、18MA)	65
圖4-20 內流場溫度分佈算例圖(12MA、15MA、18MA)	66
圖4-21 內流場密度分佈算例圖(12MA、15MA、18MA)	66
圖4-22 12MA分佈圖 67
圖4-23 12MA壓力線分佈圖 68
圖4-24 15MA分佈圖	69
圖4-25 18MA分佈圖	70
圖4-26 12MA壓力內流場比較圖	71
圖4-27 12MA溫度內流場比較圖	72
圖4-28 12MA密度內流場比較圖	73
圖4-29 12MA中心壓力切面圖	74
圖4-30 12MA 中心溫度切面圖	74
圖4-31 15MA 壓力內流場比較圖 75
圖4-32 15MA溫度內流場比較圖	76
圖4-33 15MA密度內流場比較圖	77
圖4-34 15MA中心壓力切面圖 78
圖4-35 15MA中心溫度切面圖 78
圖4-36 18MA壓力內流場比較圖	79
圖4-37 18MA溫度內流場比較圖	80
圖4-38 18MA密度內流場比較圖	81
圖4-39 18MA中心壓力切面圖	82
圖4-40 18MA中心溫度切面圖	82
圖4-41 突張結構產生的迴流區分佈圖 83
圖4-42 凹槽結構產生的迴流區分佈 83
圖4-43 凹槽後方斜率結構產生的迴流區分佈 84
圖4-44 凹槽與凹槽後方斜率比較圖 84
圖4-45 12MA各個外型氫、氧分佈圖 85
圖4-46 12MA各個外型混合效率比較圖 85
圖4-47 15MA 各個外型氫、氧分佈圖 86
圖4-48 15MA各個外型混合效率比較圖 86
圖4-49 18MA各個外型氫、氧分佈圖 87
圖4-50 18MA各個外型混合效率比較圖 87
圖4-51 12MA算例全壓損失率 88
圖4-52 15MA算例全壓損失率 88
圖4-53 18MA算例全壓損失率 89

[1]http://www.uq.edu.au/news/index.html?article=9258 
[2]http://www.nasa.gov/missions/research/x43-main.html 
[3]http://www.airforce-technology.com/projects/x51-wave-rider/ 
[4]Bird, G. A., “Approach to Translational Equilibrium in a Rigid Sphere Gas,” Phys. Fluids Vol.6, pp. 1518-1519, 1963.
[5]Bird, G. A., “The Velocity Distribution Function Within a Shock Wave, ”Journal of  Fluid Mechanics ,Vol.30,pp .479-487,1967.
[6]Bird, G. A., “Monte-Carlo Simulation in an Engineering Context,” Proger. Astronaut. Aeronaut. , Vol. 74,pp.239-255, 1981.
[7]Bird, G. A., “Perception of Numerical Method in Rarefied Gas Dynamics,” in Rarefied Gas Dynamics: Theoretical and Computational Techniques of Progress in Aeromautics and Astronautics, AIAA, Washington, Vol.118, 1989.
[8]Muntz, E. P., “Rarefied gas dynamics,” Annu. Rev. Fluid Mech. 21, pp. 387-417, 1989.
[9]Cheng, H. K., and Emmanuel, G., “Perspectives on hypersonic non-equilibrium flow,” AIAA J., Vol.33, pp.385-400, 1995.
[10]Bird, G. A., “Recent Advance and Current Challenges for DSMC,” Comput. Math. Appl., Vol.35, pp.1-14, 1998.
[11]Guy, G. W. and Mackey E.A., “Initial Wind Tunnel Test at Mach 4 and 7 of Hydrogen Burning Airframe Integrated Scramjet Fuel Injection for Supersonic Combustion,” AIAA Paper, pp.79-7405, April 1979.
[12]雷聲遠, “高速內流場之流場分析”, 國立中央大學機械系, 桃園, 1987.
[13]P. J. Lu and K. C. Wu, “On the shock enhancement of confined supersonic mixing flows,” Phys. Fluids, December 1991.
[14]P. J. Lu and K. C. Wu, “Numerical investigation on the structure of a confined supersonic mixing layer,” Phys. Fluids, December 1991. 
[15]劉陵,張榛,唐明,劉敬華,楊雨新,王興甫, “氫燃料超音速燃燒室實驗研究”, 航空動力學報, 第6卷, 第3期, 7月, 1991.
[16]Kuh-Wen Hwang, “Effect of Fuel Inlet Location on Combustion Characteristics Within Scramjet Dump Combustor,” IAF,1993 .
[17]Minucci, M. A. S., Nagamatsu, H. T., and kim, S. C., “Investigation of a 2-D Scramjet Inlet, M∞=8.25 and To=4100K,”AIAA Paper, pp.91-5018,1991
[18]Minucci, M. A. S., and Nagamatsu, H. T., “Investigation of a Two-Dimensional Scramjet Inlet, M∞=8-18 and To =4100K,”Journal of Propulsion and Power, Vol.9, No.1, pp.139-145, 1993.
[19]黃顧文, “超音速衝壓引擎中燃料噴入位置及震波結構對燃燒室特性影響分析”, 國立中央大學機械系, 桃園, 1994.
[20]Chan H. Chung, Suk C. Kim, Kenneth J. De Witt, Henry T. Nagamatsu, “Numerical Analysis of Hypersonic Low-Density Scramjet Inlet Flow,” Journal of Spacecraft and Rockets, Vol.32, No.1 ,pp.60-66, 1995.
[21]Emami, S., Trexler, C.A., Auslender, A. H., Weidner, J. P., “Experimental Investigation of Inlet-Combustor Isolators for a Dual-Mode Scramjet at a Mach Number of 4,” NASA Technical Paper 3502, 1995.
[22]Zuu-Chang Hong and Ming-Shiann Shin, “Shock-Induced mixing enhancement of supersonic Turbulent Reacting”, 46th International Astronautical Congerss, October, 1995.
[23]Mohieldin, T. O., Abdel-salam, T. M. and Tiwari, S. N., “Numerical Study of Supersonic Mixing and Combustion Using Unstructured Gird,” 38thAerospace Sciences Meeting & Exhibit, 2000.
[24]Adela Ben-Yakar, Ronald K. Hanson, “Cavity Flame-Holders for Ignition and Flame Stabilization in Scramjet: An Overview,” Journal of Propulsion and Power, Vol.17, No.4, 2001.
[25]Situ, M., Wang, C. and Zhuang, F. G., “Investigation of Supersonic Combustion of Kerosene Jets with Hot Gas Piloted Energy and Dual-Cavity,”40th AIAA Aerospace Sciences Meeting &Exhibit, 2002.
[26]J. Philp Drummond, Glenn S Diskin, Andrew D. Cutler, “FULE-AIR MIXING AND COMBUSTION IN SCRAMJETS,”38th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit, 2002.
[27]S. N. Tiwari and A. A. Taha and T. O. Mohieldin, “Effect of Pilot Injection and Combustor Configuration on the Flame Holding and Stability of Ethylene in Scramjet Engines,”40th AIAA aerospace Sciences Meeting & Exhibit, 2002.
[28]Dean Andreadis, “Scramjets integrate air and Space”, American Institute of Physics, September 2004
[29]J. Steelant, A. Mark, K. Hannemann and A. D. Gardner, “Comparison of Supersonic Combustion Tests With Shock Tunnels, Fight and CFD” , American Institute of Aeronautics and Astronautics, 2006.
[30]劉登, “概論超音速燃燒推進載具氣動熱力特性與流場數值模擬”, 中山科學研究院第二研究所, 2007
[31]Francesco Battista, Luigi Cutrone, Stefano Amabile, Giuliano Ranuzzi, “Supersonic Combustion Models Application for Scramjet Engines”, 16th AIAA/DLR/DGLR International Space Planes and Hypersonic System and Technologies Conference, 2009.
[32]施明憲, “超音速衝壓引擎燃燒室內震波增強混合與燃燒分析模擬”, 國立中央大學機械系, 桃園, 1995.
[33]洪念慈, “以直接模擬蒙地卡羅法計算三維微管流場”, 淡江大學機械與機電工程學系, 台北, 2004.
[34]羅文彬, “以直接模擬蒙地卡羅法模擬二維微管流場”, 淡江大學機械與機電工程學系, 台北, 2004.
[35]沈青, “稀薄氣體動力學”, 國防工業出版社, 北京, 2003.
[36]李峻溪,厲復霖, “雙模式極超音速衝壓引擎進氣道之數值研究”, 中正嶺學報, 第35卷, 第1期, 2006.
[37]林雅茹, “以直接模擬蒙地卡羅法計算三維微管流場與熱傳特性探討”, 淡江大學機械與機電工程學系, 台北, 2007.
[38]林建宏, “超音速燃燒衝壓引擎流場之直接模擬蒙地卡羅法模擬探討”, 淡江大學機械與機電工程學系, 台北 2009.
[39]潘穎哲, “以直接模擬蒙地卡羅法模擬微流道之氣體流場與熱傳特性分析探討”, 淡江大學機械與機電工程學系, 台北, 2009.
[40]王振國,梁劍寒,丁猛,范曉檣,吳繼平,林志勇, “高超音速飛行器動力系統研究進展”,國防科技大學航天與材料工程學系,第39卷,第6期,湖南,2009.
[41]林紘緯, “直接模擬蒙地卡羅法於超音速燃燒衝壓引擎之增強內流場混合效率分析”, 2010.
[42]http://jetex.org/images/archive/smoky/0501/3view-x43a-nasa.gif
[43]http://upload.wikimedia.org/wikipedia/commons/archive/d/de/20080616222506!Scramjet--X-43A.jpg
[44]http://tw.knowledge.yahoo.com/question/question?qid=1011020703237