臺灣博碩士論文加值系統

本研究利用分區拆解計算方法，並應用緩衝區(Buffer zone)的概念作為分區銜接的方式，以解決在地形劇烈變化處產生的計算不穩定或分區尺度差異太大的問題。緩衝區的概念是將地形劇烈變化處視為灰色地帶且不予計算，或者將緩衝區應用於不同大小尺度區域以達到分區同時計算之效果。此銜接方法雖無法得到相當準確的數值解，但對於地形變化過度劇烈、空間尺度差異過大而導致數值模擬無法計算的問題時，為相當值得採用的方法。

為探討緩衝區銜接方式的可行性，故將緩衝區概念應用於擬三維地下水模式和水平二維水理模式中，並設計單連通區域和複連通區域的案例進行測試。利用克利金法(Kriging)作為內外插工具，並藉由反覆疊代的過程直到緩衝區兩側邊界值滿足收斂條件為止。其中緩衝區在波潮流模式需滿足質量守恆而地下水模式則需滿足能量守恆的方式來銜接。

利用實測資料與緩衝區概念進行模擬，達到不同尺度的網格同時計算的目的，例如河海相接、濁水溪平原設置集水廊道和基樁等案例。

經過無數的測試與案例模擬，結果顯示緩衝區銜接分區的計算能得到相當不錯的成果，顯現緩衝區概念的可行性。緩衝區不規則性的切割法展現其威力與彈性。

In the numerical modeling of flows in oceans, estuaries, rivers or groundwater, the effect of abrupt changes in terrain usually causes numerical instability. In this work,
a numerical treatment using “hydraulic buffer zone” for solving such computational difficulty is proposed. Hydraulic buffer zone is set to be located at the area with
abrupt changes in terrain, and decomposes the computational domain into two subregions. Direct simulation of the flow in the hydraulic buffer zone is ignored, and the connection of the regions divided by the buffer zone is performed using the condition of mass or energy conservation. Hydraulic buffer zone also can be applied to the connection of two computational regions with different spatial resolutions.

In this research, the connecting schemes using hydraulic buffer zone in the modeling of groundwater flow and estuary hydrodynamics are developed. The kriging method is employed to interpolate/extrapolate the water head and water surface elevation across the buffer zone in the groundwater flow modeling and the estuary hydrodynamics modeling, respectively. At every time step, the iterative computations across the buffer zone are performed until the results of two successive iterations converge.

To test the capability of hydraulic buffer zone, several calculations for simplified problems in groundwater and estuary hydrodynamics modeling are performed and reveal reasonable results, including irregular boundary shapes.


Applications to the groundwater modeling in Zhuoshui river and the hydrodynamics modeling in the estuary of the Mei-Luan Creek are given. The results show that the treatment of hydraulic buffer zone is capability of handling the numerical instability caused by the abrupt changes in terrain and connecting two computational regions with different resolutions.

口試委員會審定書 II
摘 要 II
Abstract III
目錄 V
表目錄 VII
圖目錄 VIII
符號表 XV

第一章 緒論 1
1.1研究動機與目的 1
1.2 文獻回顧 2
1.3 研究內容與方法 3
1.4 章節介紹 3

第二章 緩衝區的建立與銜接方法 4
2.1地下水模式 4
2.1.1緩衝區之設置與計算步驟 5
2.1.2 地下水模式垂直虛擬分層的改善 12
2.2 波潮流模式 15
2.2.1緩衝區之設置與計算步驟 17

第三章 緩衝區的測試和驗證 22
3.1地下水模式之測試 22
3.1.1 緩衝區之測試 24
3.1.2 單連通區域測試 39
3.1.3 複連通區域測試 46
3.2波潮流模式之測試 52
3.2.1斜向緩衝區應用於迴水模式 55
3.2.2緩衝區應用於地形突變之測試 58
3.2.3 單連通區域測試 60
3.2.4 複連通區域測試 66

第四章 緩衝區的實際應用 76
4.1 地下水濁水溪案例應用 76
4.1.1 基樁設置於濁水溪平原之模擬 88
4.1.2 集水廊道設置於濁水溪平原之模擬 92
4.2 波潮流模式在美崙溪河口海域之應用 97
4.2.1 美崙溪河海相接之模擬(單連通區域) 101
4.2.2 美崙溪河口底床大型坑洞之模擬(複連通區域) 108

第五章 結論與建議 116
5.1 結論 116
5.2 建議 117

參考文獻 118

1. Chu, W.S. , Barker, B.L. and Akbar, A.M. ”Modeling tidal transport in the Arabian Gulf ”, ASCE, J. of Waterway, Port, Coastal, and Ocean Engineering,114, 4, pp. 37-53, 1988.

2. Goede, E.D. , Groeneweg , J. Tan, K.H., Borsboom, M.J.A. and Stelling, G.S., “A Domain Decomposition Method for the Three-Dimensional Shallow Water Equations”, Simulation Practice and Theory,3,307-325, 1995.

3. Lin, M-Y. and Huang, L.-H., “Computation of nonlinear vertical variation of depth integrated velocity profile”, In Proc. Of Int. Conf. on Hydrodynamics, Tainan, Taiwan,395-400, 2002.

4. Lin, M-Y. and Huang, L.-H., “Velocity profiles of nonlinear shallow-water flows”, Journal of the Chinese Institute of Engineers, 31 (1), 105-120, 2008.

5. Leendertse, J.J. and S.K. Liu , ”A Three-Dimensional Model for Estuaries and Coastal Seas, Vol II, Aspects of Computation ”, R-17640-OWRT, Rand Corporation, Santa Monica,CA, pp. 1-29., 1975.

6. Liu, S.K. and J.J. Leendertse, ”Multidimensional numerical modeling of esturies and coastal seas”, Adv. Hdroscienses II,pp.95-164, 1978.

7. Paglieri, L., Ambrosi, D., Formaggia, L. Quarteroni, A. and Scheinine, A.L. ”Parallel computation for shallow water flow : a domain decomposition approach“, Parallel computing,23,1261-1277, 1997.

8.Richtmeyer, R. D. and Morton, K. W. ”Difference methods for initial value problems ”, J. Hydraul. Dev., Am. Soc. Civ.,1967.

9.黃良雄，朱文生，”台灣沿海地區港灣與河口水理水質模式之發展”，國立台灣大學土木工程研究所水利組研究報告第7960號，1990.

10.吳仁友，”擬似三維海岸水動力計算模式發展”，國立台灣大學土木工程研究所碩士論文，1997.

11.林孟郁，黃良雄，”花蓮市排水及美崙溪河口淤砂之研究及改善計畫”，經濟部水利處，1999.

12.蔡東霖，”大區域地下水超抽導致地層下陷模式之發展與應用”，國立交通大學博士論文，2001.

13.吳仁友、黃良雄、曾鈞敏，”感潮河口之水理計算-以鹽水溪為例”，中國土木工程學刊，第十三卷，第三期，533頁-541頁.，2001.

14.林孟郁，黃良雄，”分區演算在感潮河口水理模擬之應用”，第25屆海洋工程研討會論文集，603-610，2003.

15.林孟郁，徐浩仁，黃良雄，”建立波潮流與海岸變遷模式”，經濟部水利署水利規劃試驗所，2003.

16.曹之獻，”擬似三維感潮河口水流及質量傳輸計算”，國立台灣大學土木工程研究所碩士論文，2005.

17.張正緯，”水平分區、垂直分層之三維地下水計算”，國立台灣大學土木工程研究所碩士論文，2007.