臺灣博碩士論文加值系統

本研究採垂直水平分離演算(VHS)之概念發展三維靜水壓水理模式，搭配紊流模式計算紊流黏滯係數，其中紊流模式分別採用雙方程模式之標準k-ε紊流模式及兩種零方程模式。在座標系統上，水平方向與垂直方向分別採用正交曲線座標系統與σ座標系統，使模式邊界能夠更貼近真實渠道之蜿蜒與底床起伏。水理模式之架構係先由深度平均之控制方程式求解水平二維每個計算格點之水面高程以及水深平均流速，並將之代入流速差異量方程式，求得三維每個計算格點相對於水深平均流速之流速差異量，即可求解出三維空間上每個計算格點之流速分佈。數值方法方面，水平、垂直模式與k、ε方程式皆採隱式法求解。分別採用直渠道案例以及急彎彎道案例進行模擬分析，並比較三種模式之紊流黏滯係數表現以及模擬結果之差異性，以探討明渠流之紊流流場現象。

A 3D hydrostatic model based on a vertical horizontal splitting (VHS) concept is developed in this study. The standard k-ε model, a two-equation turbulent model, and two kinds of zero-equation models are adopted to calculate eddy viscosity. The orthogonal curvilinear coordinate system and the sigma coordinate system are used to cope with the irregularity of channel geometry. The water elevation and the depth-averaged velocity will be solved by the 2D depth-averaged model, and then the velocity profile along the vertical direction will be solved by the velocity defect model. The implicit numerical schemes are used to discrete all of the equations to preserve the model stability unconditionally. Two experimental cases including the flow in straight channel and sharp bend were simulated by the model. Through the comparison between the experimental data and simulation results, the eddy viscosity computed from two-equation and zero-equation turbulent models were examined and discussed in depth.

目錄
摘要 I
ABSTRACT II
致謝 III
目錄 IV
表目錄 VI
圖目錄 VII
符號表 IX
第一章 緒論 1
1.1 研究動機與方向 1
1.2 文獻回顧 1
1.3 研究目的與方法 3
1.4 章節介紹 4
第二章 理論基礎 5
2.1 水理模式控制方程式 5
2.1.1 三維動量方程式 5
2.1.2 水平二維模式 6
2.1.3 速度差異量方程式 7
2.2 輔助關係式 9
2.2.1 層流與紊流應力 9
2.2.2 底床剪應力 9
2.2.3 垂向剪應力與流速梯度之轉換 10
2.3紊流黏滯係數 10
2.3.1 k-ε模式 10
2.3.2 零方程模式 11
2.4邊界條件設定 11
2.4.1 二維水平部分 11
2.4.2 垂直部分 12
2.4.3 k-ε模式 12
2.4.4 牆函數 12
第三章 數值架構 16
3.1水平二維模式 16
3.2 流速差異量方程式 18
3.3 k-ε模式 19
第四章 結果與討論 23
4.1 直渠道案例 23
4.1.1 數值參數敏感度測試 23
4.1.2 模擬結果與實驗之比對分析 24
4.1.3不同紊流模式之模擬比對分析 25
4.2 急彎案例 26
4.2.1 數值敏感度測試 26
4.2.2 模擬結果與實驗之比對分析 27
4.2.3 不同紊流模式之模擬比對分析 27
第五章 結論與建議 52
5.1 結論 52
5.2 建議 52
參考文獻 54
附錄A 57

顏志偉(1995)，｢複雜自由液面之模擬與分析｣，國立成功大學水利及海洋工程研究所博士論文。
謝德勇(2003)，｢二維水理、汙染傳輸及沉滓運移模式之研發與應用｣，國立交通大學土木工程研究所博士論文。
洪聖翔(2011)，｢正交曲線座標擬似三維水理模式於彎道水流之模擬研究｣，國立交通大學土木工程研究所碩士論文。
鐘浩榮(2012)，｢河道三維高含砂水流沉滓運移模式發展與應用｣，國立交通大學土木工程研究所博士論文。
陳蓉瑩(2013)，｢三維大渦模式在明渠流之應用｣，國立交通大學土木工程研究所碩士論文。
C.J. Chen & S.Y. Jaw, 1997.Fundamentals of Turbulence Modeling. Taylor & Francis.
C. W. LI & T.S. YU,1996. Numerical Investigation of Turbulent Shallow Recirculating Flows by a Quasi-three-dimensional k-ε model. International Journal for Numerical Methods in Fluids, vol. 23, 485-501.
Elder, J. W.,1959.The dispersion of marked fluid in turbulent shear flow. Journal of Fluid Mech., 5:554-560.
George L. Mellor & Alan F. Blumberg, 1985. Modeling Vertical and Horizontal Diffusivities with the Sigma Coordinate System. Mon. Wea. Rev., 113, 1379–1383.
Gosman, A.D., Khail, E. E., Whitelaw, J.H., 1974. The calculation of two-dimensional turbulent recirculating flows, Turbulent Shear Flow I, Springer Verlag, Heidelberg.
Ha Minh, H. and Chassaing, P., 1977. Pertubation of turbulent pipe flow, Proceedings Symposium on Turbulent Shear Flows, Pennsylvanis State University.
H. Tennekes& J. L. Lumley, 1972. A First Course in Turbulence. MIT Press, Cambridge, MA
I. Nezu& H. Nakagawa, 1993. Turblence in Open-Channel Flows. Balkema, Rotterdam.
Jin, X. and Kranenburg, C., 1993.Quasi-3D Numerical Modeling of Shallow Water Circulation. Journal of Hydraulic Engineering, 119:458-472.
Jobson, H. E. and Sayre, W. W., 1970.Vertical transfer in open channel flow. Journal of the Hydraulics Division 96
Keller, R.J. and Rastogi, A.K.,1975.Prediction of flow development on spillways, J. of the Hydraulics Division, ASCE, HY9,pp.1171-1184.
Lardner R.W. and Cekirge, H.M., 1988.A New Algorithm for 3-Dimensional Tidal and Storm-Surge Computations.Appl Math Model, 12:471-481.
Launder, BE & Spalding, D.B.,1974. The numerical computation of turbulent flows. Comp. Methods in Appl. Mech. and Engrg., 3, 269-289.
Leschziner, M.A. and W. Rodi, 1979. Calculation of strongly curved open channel flow. J. Hydr. Division., 105: 1297-1314.
Naot, D. and Rodi, W, 1982. Calculation of secondary currents in channel flows. Journal of the Hydraulic Division, ASCE, 108(8): 948-968.
Nezu, I., W. Rodi,1986. Open-channel flow measurements with a laser Doppler anemometer, J. Hydraul. Eng., 1125, 335–355.
Peter S. Bernard & James M. Wallace, 2001.Turbulent flow: analysis,measurement, and prediction. Hoboken, N.J.: Wiley.
Rodi, Wolfgang & International Association for Hydraulic Research, 1993. Turbulence models and their application in hydraulics: a state-of-the-art review (3rd ed). Balkema, Rotterdam
Rozovskii,I.L.,1957. Flow of Water in Bends of Open Channels.Ac. Sc. Ukr. SSR; Isr.Progr.Sc. Transl., Jerusalem.
S. B. Pope,1978. The Calculation of Turbulent Recirculating Flows in General Orthogonal Coordinates, Journal of Computational Physics 26, 197-217.
Shih, T.H., Liou, W.W., Shabbir, A., Yang, Z., and Zhu, J., 1994.A new k-epsilon Eddy-Viscosity Model for High Reynolds Number Turbulent Flows: Model development and validation. Computers and Fluids 24, no. 3, pp. 227-238.
Singhal, A.k. and Spalding, D. B., 1975.Prediction of two-dimensional boundary layers with the aid of the k-e model of turbulence, Imperial College, Mech. Eng.
S. Poroseva& G. Iaccarino, 2001. Simulating separated flows using the k-εmodel. Center for Turbulence Research Annual Research Briefs.
Stephenson, P. L., 1976. Theoretical study of heat transfer in two dimensional turbulent flow in a circular pipe and between parallel and diverging plates, Int. J. Heat Mass Transfer, Vol. 19, pp. 413-423.
Svesson, U.,1978. Mathematical model of the seasonal thermocline, University of Lund, Depth.Of Water Resource Engg., Sweden, Rept. No. 1002.
Wang, K. H., 1994. Characterization of Circulation and Salinity Change in Galveston Bay. J EngMech-Asce, 120: 557-579.
Wilcox, D. C., 2000. Turbulence Modeling for CFD, 2rd edition, DCW Industries, Inc., La Canada CA