臺灣博碩士論文加值系統

本研究利用數值方法分析可壓縮流在三維垂直管道中的流動及熱傳機制。流場利用有限差分法進行計算，計算方法可分為兩部分：第一部份為非黏滯性項的尤拉方程式採用Roe方法計算通量，並且加入Preconditioning矩陣，讓程式在計算低速可壓縮流可獲得良好之收斂結果，而程式因為在使用Preconditioning時，加入Artificial time term 時，已破壞了整個統御方程式，因此需使用Dual time stepping疊代使其在Artificial domain 收斂時才能進入下一個真實時階；第二部份為黏滯性項的計算，採用二階中央插分法。在時間項方面則採用LUSGS隱式法，利用LUSGS疊代以求出下一時階物理量。出口設非反射性邊界條件避免可壓縮流中壓力波的干擾。在許多應用例子中，溫差常常大於30K，因此Boussinesq assumption不適用。本文採用OpenMP方法提升運算速度。
由數值計算的結果得知，三維垂直管道之自然對流，因浮力效應往上推升，在出口端產生最大速度，雷諾數(Re=400)之熱傳效應較雷諾數(Re=100、200)之熱傳還差，主要是在出口端雷諾數(Re=100、200)所吸進的流量，大於雷諾數(Re=400)所吸進的流量，而雷諾數(Re=100)從出口端吸入的速度衝擊到出口速度，故熱傳較雷諾數(Re=200)小。雷諾數(Re=100、200)流場之迴流接近入口處，在流場中會有速度反曲點，造成流場的不穩定。

An investigation of heat transfer in a three-dimensional tapered chimney with consideration of the flow compressibility is studied numerically.The finite difference method is adopted and the computational approaches are divided into two parts. One is the Roe scheme applied for the flux of inviscid terms and the preconditioning matrix is added for the efficiency in all speed fields. The other one is the central difference method of second order utilized to solve viscous terms. The temporal term is solved by LUSGS. Non-reflection conditions at the outlet is derived in order to resolve reflections induced by acoustic waves. In many important natural convection problems, the temperature differences are often higher than 30K. Boussinesq assumption is unreasonable. Besides, the OpenMP method is also used to promote the computing efficiency.
By numerical results, there is the greatest flow speed near the outlet in the three-dimensional vertical natural convection pushed upward by buoyancy effect. The enhancement of heat transfer of Reynolds number 400 is worse than the enhancement of Reynolds number 100, 200. It is mainly due to the more flow rate sucked from exterior near the outlet in the case of Reynolds number 100, 200. And the flow sucked from the exterior impacted the flow exiting the outlet in the case of Reynolds number 100. The heat transfer is worse than the case of Reynolds number 200. Besides, the flow field is unstable due to the backflow near the outlet in the case of Reynolds number 100, 200.

目錄
中文摘要 i
英文摘要 ii
目錄 iii
表目錄 iv
圖目錄 v
符號表 vii
一、緒論 1
二、物理模式 9
2-1物理尺寸與分析模式 9
2-2分析假設及統御方程式 10
2-3邊界條件 12
三、數值計算模式 14
3-1統御方程式 15
3-2 Roe scheme 17
3-3 MUSCL法 24
3-4 Preconditioning 26
3-5 Dual time stepping 32
3-6 LUSGS法 34
3-7非反射性邊界 36
四、結果與討論 40
五、結論 88
參考文獻 89
附錄 93

1. A.E. Bergles, 「Recent Development in Convective Heat-transfer Augmentation,」 Applied Mechanics Reviews, vol.26,pp.675-682, 1973.
2. A.E. Bergles, 「Survey and Evaluation of Techniques to Augment Convective Heat and Mass Transfer,」Heat and Mass Transfer, vol.1,pp.331-424,1969.
3. D. B. Ingham , D. J. Keen and P. I. Heggs,」 Two Dimensional Combined Convection in Vertical Parallel Plate Ducts, Including Situations of Flow Reversal」 International Journal for Numerical methods in engineering, vol. 26, pp.1645-1664, 1973
4. A. Behzadmehr, N. Galanis and A. Laneville, 「Low Reynolds numer mixed convection in vertical tubes with uniform wall heat flux」, International Journal of Heat and Mass Transfer, vol. 46, pp. 4823-4833, 2003
5. A. Barletta, 「Analysis of flow reversal for laminar mixed convection in a vertical rectangular duct with one or more isothermal walls」, International Journal of Heat and Mass Transfer, vol. 44, pp. 3481-3497, 2000
6. A. Barletta, 「Heat transfer by fully developed flow and viscous heating in a vertical channel with prescribed wall heat fluxes」, International Journal of Heat and Mass Transfer, vol. 42, pp. 3873-3885, 1998
7. A. Barletta, Laminar convection in a vertical channel with viscous dissipation and buoyancy effects, International Journal of Heat and Mass Transfer, vol. 41, , pp. 3501-3513, 1998
8. L. Martinez-Suastegui and C. Trevino, 「Transient laminar opposing mixed convection in a differentially and asymmetrically heated vertical channel of finite length」, International Journal of Heat and Mass Transfer, vol. 51, pp. 5991-6005, 2008
9. C. S. Yang, D.Z. Jeng, K. A. Yih, C. Gau and Win Aung, 「Numerical and Analytical Study of Reversed Flow and Heat Transfer in a Heated Vertical Duct」, Journal of Heat Transfer, vol. 131, Pages 072501-1~9, 2009
10. K. Boulama and N. Galanis, 「Analytical Solution for Fully Developed Mixed Convection Between Parallel Vertical Plates With Heat and Mass Transfer」, Journal of Heat Transfer, vol. 126, pp. 381-388, 2004
11. G. Desrayaud and G. Lauriat, 「Flow reversal of laminar mixed convection in the entry region of symmetrically heated, vertical plate channels」, International Journal of Thermal Sciences, vol. 48, Issue 11, pp. 2036-2045, 2009
12. M. Zghal, N. Galanis and C. T. Nguyen , 「Developing mixed convection with aiding buoyancy in vertical tubes: a numerical investigation of different flow regimes'』, International Journal of Thermal Sciences, vol. 40, Issue 9, October 2001, Pages 816-824
13. T. H. Mai, A C. V. Popa and A G. Polidori, 「 Transient mixed convection flow instabilities in a vertical pipe」, Heat Mass Transfer , vol. 41, Issue 9, 14 April 2003, Pages 216-225
14. Mayur T. Bhoite, G.S.V.L. Narasimham and M.V. Krishna Murthy,「Mixed convection in a shallow enclosure with a series of heat generating components」, International Journal for Thermal Sciences, vol. 44, pp. 121–135, 2005.
15. C. T. Nguyen , S. B. Maiga , M. Landry, N. Galanis and G. Roy, 「Numerical investigation of flow reversal and instability in mixed laminar vertical tube flow」, International Journal of Thermal, vol. 43,pp. 797–808, 2004
16. H. Koizumi, Y. Umemura, S. Hando and K. Suzuki, 「Heat transfer performance and the transition to chaos of mixed convection around an isothermally heated sphere placed in a uniform, downwardly directed flow」International Journal of Heat and Mass Transfer, vol. 53, Issues 13-14, June 2010, Pages 2602-2614.
17. Y. Stiriba, F.X. Grau, J.A. Ferre and A. Vernet , 「A numerical study of three-dimensional laminar mixed convection past an open cavity」, International Journal of Heat and Mass Transfer, vol. 53, Issues 21-22, October 2010, Pages 4797-4808
18. N. Laaroussi, G. Lauriat and G. Desrayaud , 「Effects of variable density for film evaporation on laminar mixed convection in a vertical channel」 , International Journal of Heat and Mass Transfer, vol. 52, Issues 1-2, 15 January 2009, Pages 151-164
19. Y. Azizi and B. Benhamou, 「Buoyancy effects on upward and downward laminar mixed convection heat and mass transfer in a vertical channel,」 Heat and mass transfer in a vertical channel vol.17 NO.3, pp.333-353,2007
20. M. Shahraeeni and M. Raisee, 「Investigation of turbulent mixed convection of air flow in vertical tubes using a zonal turbulence model」, International Journal of Heat and Fluid Flow, vol. 31, pp.179–190, 2010
21. C. Balaji, M. Holling and H. Herwig , 「A temperature wall function for turbulent mixed convection from vertical, parallel plate channels 「, International Journal of Thermal Sciences, vol. 47, Issue 6, June 2008, Pages 723-729
22. W.-S. Fu ,C.-C. Li , C.-P. Huang ,and J.-C. Huang ,」An Investigation of a High Temperature Difference Natural Convection in a Finite Length Channel without Boussinesq Assumpion,」 International Journal of Heat and Mass Transfer,vol.52,pp.2571-2580,2009.
23. 李崇綱「黏性流體之全域速度場數值解法研究」, 國立交通大學 碩士論文(2006)
24. E. Turkel, 「Preconditioned method for solving the incompressible and low speed compressible equations,」 Journal of Computational Physics, vol.72, pp.277-298, 1987
25. T.J. Poinsot and S.K. Lele,」Boundary Conditions for Navier–Stokes,」 Journal of Computational Physics, vol. 101,pp. 104–129,1992.
26. P.L. Roe,」Approximation Riemann Solver, Parameter Vectors, and Difference
Schemes,」 Journal of Computational Physics, vol. 43, pp. 357–372, 1981.
27. I. Abalakin, A. Dervieux and T. Kozubskaya,」 A Vertex Centered High Order MUSCL Scheme Applying to Linearised Euler Acoustics,」 INRIA, no. 4459, 2002.
28. J.M. Weiss and W.A. Simth, 「Preconditioning Applied to Variable and Constants
Density Flows,」 AIAA ,vol. 33,pp. 2050–2056,1995.
29. X.F.Xu,J.S.Lee and R.H.Pletcher,「 A compressible finite volume formulation for large eddy simulation of turbulent pipe flows at low Mach number in Cartesian coordinates.」, J.Comput.Phys.,vol.203, pp.22-48, 2005
30. S. Yoon and A. Jamesont, 「Lower-Upper Symmetric-Gauss-Seidel Method for the Euler and Navier-Stokes Equations」, AIAA,vol. 26,no. 9,1988.
31. J. Dennis, P. Thomas, and B. Pieter, Recent Enhancements to OVEERFLOW, Aerospace Sciences Meeting and Exhibit, 35th, Reno, NV, 1997