第五章 結論與未來展望
5.2 未來展望
（1） 本文已成功發展出一套模組化且有系統的數值計算程式於不可 壓縮雷諾平均奈維爾－史托克方程式中，在與文獻的比較中可 以發現使用越高階的空間離散方法，在雷諾數增加的同時，其 數值模擬解則與實驗結果更接近，所以往後可將此數值模擬方 法改用較為高階的空間離散法，觀察其結果是否如文獻所記載 般精確。
（2） 本文已成功發展出一套模組化且有系統的數值計算程式於不可 壓縮雷諾平均奈維爾－史托克方程式中，但根據文獻計載其數 值方法（LU-SGS 算則）在可壓縮雷諾平均奈維爾－史托克方程 式中亦可得到相當不錯的結果，所以往後可將此數值方法推展
應用於可壓縮流場中，觀察其結果是否亦如不可壓縮流場般精 確。
（3） 對於複雜之幾何外形，應可採多區塊之方式處理，並建立多重 網格的計算功能，使其能足以應付更複雜流場的計算。
（4） 建立雙時步次疊代的計算方式，對於較複雜之非穩態流場問題 能有效解決；並因應配合電腦的快速發展，建立平行運算的功 能，以期能更加快速且有效率的分析各式各樣的
參考文獻
1. Anderson , D.A. , Tannehill , J.C. and Plether ,
R.H.(1984),”Computational Fluid Mechanics and Heat Transfer”, McGraw-Hill , New York , 1984.
2. Harlow , F.H. and Welch , J.E.(1965),”Numerical Calculation of Time-Dependent Viscous Incompressible Flow with Free Surface”, ThePhysics of Fluids , Vol.8 , No.12 , pp2182-2189.
3. Patankar , S.V. and Spalding , D.B.(1972),”A Calculation Procedurefor Heat , Mass and Momentum Transfer in
Three-Dimensional Parabolic Flows”, International Journal of Heat and Mass Transfer , Vol.15 , pp.1787-1806.
4. Patankar , S.V.(1980),”Numerical Heat Transfer and Fluid Flow”,Mcgraw-Hill, New York.
5. Chorin , A .J.(1967),”A Numerical Method for Solving
Incompressible Viscous Flow Problems”, Journal of Computational Physics , Vol.2 , pp12-26.
6. Harten , A.(1983),”High Resolution Schemes for Hyperbolic Conservation Laws “, Journal of Computational Physics , Vol.49 , pp.357-393.
7. Yee , H.C.(1989),”A Class of High-Resolution Explicit Shock-Capturing Methods”, NASA TM-101088.
8. Shu , C.W.(1987),”TVB Uniformly High-Order Schemes for Conservation Laws”, Mathematics of Computation , Vol.49 , pp.105-121.
9. Harten , A. and Osher , S. (1987),”Uniformly High-Order Accurate Non-Oscillatory Schemes I”, SIAM Journal of Numerical
Analysis ,Vol.24 , No.2 , pp.200-212.
10. Liu , X.D. , Osher , S. and Chan , T.(1994),”Weighted Essentially Non-Oscillatory Schemes”, Journal of Computational Physics , Vol.115 , pp.200-212.
11. Jiang , G.S. and Shu , C.W.(1996),”Efficient Implementation of
Weighted ENO Scheme”, Journal of Computational Physics , Vol.126 , pp.202-228.
12. Harten , A. , Engquist , B. , Osher , S. and Chakravarthy ,
S.(1987),”Uniformly High-order Accurate Non-Oscillatory Schemes Ш”, Journal of Computational Physics , Vol.71 , pp.231-303.
13. Rogers , S.E., Kwak, D.(1990),”Upwind Differencing Schemes for the Time-Accurate Incompressible Navier-Stokes Equations”, AIAA Journal , Vol.29 , pp.603-610.
14. Yoon , S. and Jameson , A.(1987),”An LU-SSOR Scheme for the Euler and Navier-Stokes Equations”, AIAA Paper 87-0600.
15. Ghia , U. , Ghia , K.N. and Shin , C.T.(1982),”High-Re Solutions for Incompressible Flow Using the Navier-Stokes Equations and a
Multigrid Method”, Journal of Computational Physics , Vol.48 , pp.387-411.
16. Kim , J. and Moin , P.(1985),”Application of a ctional-Step Method to Incompressible Navier-Stokes Equations”, Journal of Computational Physics , Vol.59 , pp.308-323.
三維空穴流部分：
[圖 5.1] 三維空穴流之流場示意圖及其邊界條件
[圖 5.2]三維空穴流計算外型及網格分布圖
[圖 5.3]三維空穴流等壓力圖(Re=100)
[圖 5.4]三維空穴流等 u 速度圖(Re=100)
[圖 5.5]三維空穴流等 v 速度圖(Re=100)
[圖 5.7]三維空穴流流線-等 U 速度圖(Y 方向切面) (RE=100)
[圖 5.8]三維空穴流流線-等 W 速度圖(Z 方向切面) (RE=100)
[圖 5.9]雷諾數為 100 時三個中剖面之向量流線圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.10]雷諾數為 100 時之等壓力圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.11]雷諾數為 100 時之等渦度圖
[圖 5.12] 雷諾數為 400 時三個中剖面之向量流線圖
Z=0.25
Z=0.50
Z=0.75
Z=0.25
Z=0.50
Z=0.75 [圖 5.14]雷諾數為 400 時之等渦度圖
[圖 5.15]三維空穴流等壓力圖(Re=1000)
[圖 5.16]三維空穴流等 u 速度圖(Re=1000)
[圖 5.17]三維空穴流等 v 速度圖(Re=1000)
[圖 5.18]三維空穴流等 w 速度圖(Re=1000)
[圖 5.19]三維空穴流流線-等 U 速度圖(Y 方向切面) (RE=1000)
[圖 5.20]三維空穴流流線-等 W 速度圖(Z 方向切面) (RE=1000)
[圖 5.21]雷諾數為 1000 時三個中剖面之向量流線圖
Z=0.25
Z=0.50
Z=0.75
Z=0.25
Z=0.50
Z=0.75 [圖 5.23]雷諾數為 1000 時之等渦度圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.25]雷諾數為 3200 時之等壓力圖
Z=0.25
Z=0.50
Z=0.75
[圖 5.27] 雷諾數為 5000 時三個中剖面之向量流線圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.28]雷諾數為 5000 時之等壓力圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.29]雷諾數為 5000 時之等渦度圖
[圖 5.30] 雷諾數為 7500 時三個中剖面之向量流線圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.31]雷諾數為 7500 時之等壓力圖
Z=0.25
Z=0.50
Z=0.75
[圖 5.33]三維空穴流等壓力圖
[圖 5.34]三維空穴流等 u 速度圖
[圖 5.35]三維空穴流等 v 速度圖
[圖 5.37]三維空穴流流線-等 U 速度圖(Y 方向切面) (RE=10000)
[圖 5.38]三維空穴流流線-等 W 速度圖(Z 方向切面) (RE=10000)
Z=0.25
Z=0.50
Z=0.75 [圖 5.40]雷諾數為 10000 時之等壓力圖
Z=0.25
Z=0.50
Z=0.75 [圖 5.41]雷諾數為 10000 時之等渦度圖
[圖 5.42]通過中心軸的速度分布比較
[圖 5.43]通過中心軸的速度分布比較
[圖 5.44]通過中心軸的速度分布比較
[圖 5.45]通過中心軸的速度分布比較
[圖 5.46]通過中心軸的速度分布比較
[圖 5.47]通過中心軸的速度分布比較
[圖 5.48]不同算則通過中心軸的速度分布比較圖
[圖 5.49] 不同算則通過中心軸的速度分布比較圖
[圖 5.50]通過中心軸的速度分布比較圖
[圖 5.51]通過中心軸的速度分布比較圖
[圖 5.52]通過中心軸的速度分布比較圖
[圖 5.54]通過中心軸的速度分布比較圖
[圖 5.55]通過中心軸的速度分布比較圖
[圖 5.56]通過中心軸的速度分布比較圖
[圖 5.58]人工壓縮因子對收斂歷線影響圖
[圖 5.59]人工壓縮因子對收斂歷線影響放大圖
[圖 5.60]時間離散方法對收斂歷線影響圖(Lusgs)
三維背向階梯流部分：
[圖 5.62] 三維背向階梯流之流場示意圖及其邊界條件
[圖 5.63]三維背向階梯流計算外型及網格分布圖 (100x20x20)
[圖 5.64]三維背向階梯流等壓力圖(Re=100)
[圖 5.65]三維背向階梯流等壓力圖(Re=100)
[圖 5.66]三維背向階梯流等 u 速度圖(Re=100)
[圖 5.67]三維背向階梯流等 u 速度圖(Re=100)
[圖 5.68]三維背向階梯流等 v 速度圖(Re=100)
[圖 5.69]三維背向階梯流等 v 速度圖(Re=100)
[圖 5.70]三維背向階梯流等 w 速度圖(Re=100)
[圖 5.71]三維背向階梯流等 w 速度圖(Re=100)
[圖 5.72]三維背向階梯流流線-等 U 速度圖(Y 方向切面) (Re=100)
[圖 5.73]三維背向階梯流流線-等 W 速度圖(Z 方向切面) (Re=100)
[圖 5.74] 雷諾數為 100 時中剖面之等速度圖
[圖 5.76] 雷諾數為 200 時中剖面之等速度圖
[圖 5.77] 雷諾數為 200 時中剖面之向量流線圖
[圖 5.78] 雷諾數為 300 時中剖面之等速度圖
[圖 5.80] 雷諾數為 400 時中剖面之等速度圖
[圖 5.81] 雷諾數為 400 時中剖面之向量流線圖
[圖 5.82]雷諾數－流離接觸點關係圖
（separated/reattached region）