本文以數值方法在對稱式背向階梯燃燒室中，探討噴嘴之位置及收縮比對流場之影響。利用模擬軟體FDS( Fire Dynamics Simulator )及大渦漩模擬法（Large Eddy Simulation）搭配Smagorinsky次網格黏滯應力模式( SGS )對流場的穩態、暫態流況及定量上的性質變化做一模擬分析並與實驗文獻所得之結果作一比較。當數值使用之入口假設條件所獲得之結果與實驗流場使用之入口條件所得結果比對後，顯示數值與實驗在長短邊再接觸長度上有明顯差異。由於數值之入口流發展為層流狀態，實驗之入口流則已為紊流狀態，因此形成誤差之存在。此差異亦呈現出層流狀態之再接觸長度較紊流狀態之再接觸長度來得長。為消除此誤差可利用改變入口流狀態以獲得改善。    
    穩態分析中，未加噴嘴之突張流場在短邊迴流區下游之三維壁面出現第三迴流區。加裝噴嘴則會造成短邊迴流區與第三迴流區的結合形成較對稱之流場形態。水平平均速度、垂直平均速度與紊流強度則因入口流狀態與實驗之差異，在階梯後方產生不同之分佈趨勢。
    暫態分析中，未加噴嘴之突張流場的壓力擾動頻率較為雜亂。加裝噴嘴則使的壓力擾動頻率有一明顯之的頻率出現。噴嘴之收縮比會使流場壓力擾動頻率有一明顯之主頻，噴嘴位置則會改變流場壓力擾動之頻率。利用數值可視化觀察發現渦流之配對與流逸發生在流場兩低壓區之間，壓力擾動之頻率則影響渦流之配對與流逸。透過三維流線與壓力之切面觀察得知角落渦漩之產生是來自第三維兩側壁面流體迴流至階梯角落所造成且其迴流與迴流區內之壓力變化有關。

This study presents the numerical simulation of a backward-facing step combustor with nozzle, results of influence of the contract ratio and position of nozzle in sudden expansion flow	were shown in this thesis. Utilizing FDS ( Fire Dynamics Simulator ) and LES( Large Eddy Simulation )with subgrid scale (SGS) viscosity stress of Smagorinsky model to simulate the steady flow, the transient flow and to quantitatively analyze then compared with the results of the document of the experiment. When the numerical results using the condition of the inlet flow of experiment compared with the experiment results, the results show that the short and long reattachment length of the simulation have obviously difference with the experiment. The reason is that the inlet flow of the simulation is laminar flow but  turbulent flow is in the experiment, so the error is produced by difference between conditions of the inlet flow. The difference between the simulation and the experiment also show that the reattachment length of laminar flow is more long than turbulent flow. By changing the condition of the inlet flow will reduce to the error and improve the result. 
    The analysis of the steady, sudden expansion flow without nozzle appears to the third recirculation region at the downstream zone of short recirculation region in three-dimensional side. After installing the nozzle, the third recirculation region combine into the short recirculation region, leading to the flow which is more symmetrical . The distribution of horizontal average velocity , vertical average velocity and turbulent intensity behind step in the simulation are different from that in experiment, since the conditions in the inlet flow state of experiment and simulation are not the same .
    The analysis of the transient, the fluctuant frequency of pressure of sudden expansion flow without nozzle is disorderly . After installing the nozzle , the fluctuant frequency of thepressure appears to a obvious peak. The contract ratio of nozzle can make fluctuant frequency of pressure have a fixed value . The position of nozzle can change the fluctuant frequency. The combination and flowing of the vortex appears Between two low-pressure zone by numerical visualization. the fluctuant frequency of pressure influences combination and shdding of the vortex .

目錄
中文摘要  .............................................. I
英文摘要  ........................................... III
目錄  ............................................... V
圖目錄  ............................................. VIII
表目錄  ............................................. XI
符號說明  ........................................... XIII

第一章  簡介  ……………………………………………… 1
    1-1  前言 ……………………………………………………  1
    1-2  文獻回顧 ………………………………………………  3
         1-2-1  流場描述  ……………………………………  3
         1-2-2  流場實驗研究  ………………………………  4
1-2-3  流場數值研究 ………………………………  6
1-2-4  研究目的 ……………………………………  8
第二章  理論分析   ……………………………………… 10
     2-1  理論概述   ………………………………………… 10
     2-2  統御方程式  ……………………………………… 10
     2-3  紊流模式  …………………………………………  11
     2-4  邊界條件   …………………………………………  14
第三章  數值方法   …………………………………… 16
3-1  FDS 簡介   ……………………………………………  16
3-2  格點系統   ……………………………………………  17
3-3  離散與疊代方法   ……………………………………  18
    3-4  收斂標準   ……………………………………………  19
3-5  參數設定   ……………………………………………  19
3-6  計算細節   ……………………………………………  20
第四章  結果與討論    ………………………………… 21
4-1 格點測試    …………………………………………  21
4-2 穩態流場特性分析  ……………………………………  22
     4-2-1 再接觸長度  …………………………………  22
     4-2-2 平均流場  ……………………………………  24
     4-2-3 紊流強度之探討  ……………………………  26 
4-3 暫態流場與壓力擾動頻譜    …………………………  27
     4-3-1 噴嘴對壓力擾動頻譜之影響  ………………  27
     4-3-2  壓力與渦流結構動態   ……………………  29
     4-3-3  三維角落渦漩之探討   ……………………  30
第五章  結論與建議  …………………………………… 32
    5-1 結論    …………………………………………………  32
   5-2 建議    …………………………………………………  32
參考文獻      …………………………………………… 35
  圖目錄
圖1-1 德國出產之Me-262 噴射戰機 ………………………… 38
圖1-2 德國製之 V-2 火箭 …………………………………… 38
圖1-3 典型之背向階梯流場 ………………………………… 39
圖1-4 Hyper-X 之數值模擬示意圖 …………………………… 39
圖1-5 燃燒室流場示意圖 …………………………………… 40
圖3-1 FDS 處理程序架構示意圖 …………………………… 41
圖 3-2 非均勻網格示意圖 …………………………………… 42
圖3-3 多重區間配置示意圖 ………………………………… 42
圖 3-4 FDS 之Smokeview 可視化外型圖 ………………… 43
圖 3-5 FDS 之外型網格圖 ………………………………… 43
圖 3-6 實驗外型圖 單位:mm ……………………………… 44
圖4-1 無加噴嘴之燃燒室網格 ……………………………… 45
圖4-2 Nd=60 cm 、Cr=1/2 之燃燒室網格 ………………… 46
圖4-3 無噴嘴之時間平均速度向量與速度為零之等高線分佈
圖 …………………………………………………………………… 47
圖4-4 Nd=60 cm、Cr=1/2 之時間平均速度向量與速度為零之等高
線分佈圖 ……………………………………………………… 48
IX
圖4-5 Nd=40 cm、Cr=1/3 之時間平均速度向量與速度為零之等高
線分佈圖 ………………………………………………………… 49
圖4-6 無噴嘴之水平平均速度……………………………… 50
圖4-7 Nd=30 cm、Cr=1/2 之水平平均速度 ……………… 51
圖4-8 Nd=60 cm、Cr=1/2 之水平平均速度 ……………… 52
圖4-9 Nd=30 cm、Cr=1/9 之水平平均速度 ……………… 53
圖4-10 無噴嘴之垂直平均速度 …………………………… 54
圖4-11 Nd=30 cm、Cr=1/2 之垂直平均速度 …………… 55
圖4-12 Nd=60 cm、Cr=1/2 之垂直平均速度 ……………… 56
圖4-13 Nd=30 cm、Cr=1/9 之垂直平均速度 ……………… 57
圖4-14 無噴嘴之垂直平均速度等輪廓分佈圖 …………… 5 8
圖4 - 1 5 無噴嘴之紊流強度 ……………………………… 5 9
圖4-16 Nd=30 cm、Cr=1/2 之紊流強度 …………………… 60
圖4-17 Nd=60 cm、Cr=1/2 之紊流強度 …………………… 61
圖4-18 Nd=30 cm、Cr=1/9 之紊流強度 ………………… 62
圖4-19 (a) 無噴嘴之壓力擾動頻譜 ………………………… 63
(b) Nd=40 cm、Cr=1/2 之壓力擾動頻譜 …………… 63
(c) Nd=40 cm、Cr=1/3之壓力擾動頻譜 ……………… 64
(d) Nd=60 cm、Cr=1/3之壓力擾動頻譜 ……………… 64
X
圖 4-20 無噴嘴之渦漩動態圖 …………………………………… 65
圖4-21 無噴嘴之壓力與渦流結構週期圖 ……………………… 66
圖4-22 Nd=40 cm、Cr=1/3 之壓力與渦流結構週期圖 ……… 67
圖4-23 長邊迴流區之三維暫態流線圖 ……………………… 68
圖4-24 長邊迴流區之水平切面壓力觀測圖 ………………… 69
圖4-25 短邊迴流區之三維暫態流線圖 ……………………… 70
圖4-26 短邊迴流區之水平切面壓力觀測圖 ………………… 71
  表目錄
表 3-1 不同噴嘴條件之搭配組合 …………………………… 44

參考文獻.

1.  Bradshaw, P., and Wong, F. Y. F., 1972 , ‘‘The Reattachment and Relaxation of Turbulent Shear Layer, ” Journal of Fluid Mechanics, Vol. 52, pp. 113-135.

2.  Kim, J., Kline, S. J., and Jonston, J. P., 1980 , ‘‘Investigation of a Reattaching Turbulent Shear Layer: Flow Over a Backward-Facing Step, ” Journal of Fluids Engineering, Vol. 102, pp. 302-308.

3.  Eaton, J. K., and Johnston, J. P., 1981 , ‘‘A Review of Research on Subsonic Turbulent Flow Reattachment, ” AIAA Journal, Vol. 19, pp. 1093-1100.

4.  Isomoto, K., and Honami, S., 1989 , ‘‘The Effect of Inlet Turbulence Intensity on the Reattachment Process Over a Backward-Facing Step, ” Journal of Fluids Engineering, Vol. 111, pp. 87-92.

5.  Adams, E. W., and Johnston, J. P., 1988 , ‘‘Flow Structure in the Near-Wall Zone of a Turbulent Separated Flow, ” AIAA Journal, Vol. 26, pp. 932-939.

6.  Durst,F., F. Pereira, J. C., and Tropea, C., 1993 , ‘‘The plane symmetric sudden expansion flow at low Reynolds numbers,’’ J. Fluid Mech., Vol 248, pp. 567

7.  Brown, G. L., and Roshko, A., 1974, ‘‘On Density Effects and Large Structure in Turbulent Mixing Layers, ’’ Journal of fluid Mechanics, Vol. 64, Part 4, pp. 775-816.

8.  Ganji, A. R., and Sawyer, R. F., 1980, ‘‘Experimental Study of the Flowfield of a Two-Dimensional Premixed Turbulent Flame, ” AIAA Journal, Vol. 18, No. 7,  pp. 817-824.

9.  Armaly, B. F., Durst, F., Pereira, J. C. F. & Schoenung, B., 1983, “Experimental and theoretical investigation of backward-facing step,” J. Fluid Mech ,Vol 127, pp. 473–496.

10.  Fearn, R., Mullin, T., Cliffe, K.A., 1990 , ‘‘Nonlinear flow phenomena in a symmetric sudden expansion,” J. Fluid Mech ,Vol 211, pp. 595–608.

11.  Chiang, T. P., Tony W. H. Sheu , and Fang, C. C., 1999, ‘‘Side wall effects on the structure of laminar flow over a plane-symmetric sudden expansion, ” Computers and Fluids, Vol. 29, pp. 467-492.

12.  De Zilwa, S. R. N., Khezzar, L. and Whitelaw, J. H. , 2000 , ‘‘ Flows through plane sudden-expansions,” International Journal for Numerical Methods in Fluids, Vol 32 , No 3 , pp. 313-329

13.  Lee, S. C., Yang, J. T., and Chao, Y. C., 2001, ‘‘Large Eddy Simulation of Transition from Transpiration to Film Cooling over the Backstep Flow, ” Book of Abstracts PSFVIP-3, March 18-21, Maui, Hawaii, U.S.A., pp.149-150.

14.  Manica R., De Bortoli A.L., 2004 , “Simulation of sudden expansion flows for power-law fluids,” Journal of Non-Newtonian Fluid Mechanics,Vol 121, pp 35-40

15.  Deardorff, J.W., 1970, “A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers,” J. Fluid Mech., vol. 41, pp. 453–480.

16.  Smagorinsky, J., 1963 , “General circulation experiments with the primitive equations in the basic experiment,”  Mon. Weather Rev., Vol 91, pp. 99-164  

17.  馬萬鈞，2000，「分離流場之蒸散冷卻現象研究」，國立清華大學動力機械工程學系博士論文。
18.  卓連旺，1994，「二維對稱突張燃燒室下游噴嘴對流場結構之影響」，國立清華大學動力機械工程研究所碩士論文。
19.  鄭勝文，1997，「噴嘴背壓對雙背向階梯流場之影響」，國立清華大學動力機械工程研究所碩士論文。
20.  黃興閎，2004，「背向階梯流場之剪流層非穩態特性研究」，國立清華大學動力機械工程研究所碩士論文。
21.  李世昭，2002，「背向階梯流場模態與側向噴流混合效應之大渦模擬分析」，國立清華大學動力機械工程研究所博士論文。