本文將三維NS-VOF模式結合賓漢雙黏性流模式，以探討海嘯湧潮與颱風洪水造成結構物基礎沖刷之問題，模式以兩黏性參數以及降伏應力表示泥沙受水流沖刷之運動行為，其中 為牛頓流體運動行為與 為泥沙固態之特性， 表示泥沙被水流帶起之臨界應力值。在進行研究之前，本文先進行模式驗證。自由液面驗證方面，本研究設計潰壩湧潮撞擊結構物陣列之實驗，以驗證模式預測自由液面發展之預測準確度，其模擬結果與實驗高速攝影機捕捉之水位十分吻合。在賓漢雙黏性流方面，與 Bird et al.（1983）所推導之解析解驗證，模擬求得之流速與其求得之解析解流速相當一致，確認模式發展無誤。在沖刷坑發展方面，採用Dey and Barbhuiya (2005) 所發表水流對於圓形橋墩之沖刷坑實驗作為驗證對象，其結果表示模式能適切描述底床沖淤之過程與範圍，且能預測最大之沖刷深度，然於下游堆積處有較明顯之差異。本模式對於水流衝擊結構物之墩前壅水現象、向下射流與沖刷坑周圍近似馬蹄型渦流之效應產生皆能確切描述。
模式驗證並瞭解其準確性後，本文以2009莫拉克颱洪之雙園大橋斷橋事件為三維真實尺寸為研究對象，分別探討未補樁之初始橋墩與補樁擴座後之兩種橋墩型態對於沖刷坑變化之影響，並進一步分析基樁之水流壓力與公路橋梁設計規範（交通部，2009）之經驗公式值之適用性。後續並探討洪水挾帶漂流物掛淤於結構物間對於底床沖刷變化之影響，其中掛淤位置分別為浮木掛淤於兩座橋墩間、浮木掛淤於橋墩上游側、浮木掛淤於群基樁間及浮木掛淤於橋墩跨距間之極端掛淤案例四類逐一探討，研究結果發現浮木掛淤於兩座橋墩間、掛淤於上游與掛淤於群基樁間對沖刷坑並無顯著之影響，然而增加浮木掛淤數量則會加速沖刷坑發展且使沖刷坑深度增大。
由以上驗證與案例分析可知，賓漢雙黏性流模式能適切的探討洪流與結構物交互作用之沖刷坑發展過程，且能避免使用過多經驗公式以及經驗參數。研究亦證明本模式能廣泛運用於湧潮或海嘯對於沖刷結構物基礎之沖刷深度預測。

We integrate the 3D NS-VOF model with Bingham Bi-viscous fluid model to investigate the scouring problems around the structures induced by tsunami bores or typhoon floods. In this model, the scour mechanism is modeled by two viscosities and a yield stress. In which is used to describe Newtonian viscosity of liquefied sediment, and is for the solid-state properties of the sediment, and is for presenting the yield stress induced by the current. Before implementing this model, benchmark tests are performed. For validating the model accuracy on predicting the free-surface kinematics, we designed a case in which a dam-break bore is impinging with a structure array. A high-speed camera is used to record the free-surface movement. Very good agreement between the model prediction and experimental data can be seen. For validating the Bingham Bi-viscous fluid model, we compare the result with the analytical solution derived by Bird et al.（1983）. The comparison shows that the numerical result is almost identical to the analytical solution in terms of the velocity distribution. For validating the profile of scour hold, we simulate the case proposed by Dey and Barbhuiya (2005) in which the current is flowing around a circular pier. The result shows that our model is able to predict the maximum scouring hole and properly describe the scour mechanism and also the scour area around the cylinder. However, larger discrepancy can be found in the accumulation zone at the downstream direction. In this case, the present model can accurately present several phenomena, such as the bow wave, downword jet current, and horseshoe vortex around the cylinder.
After validating the model, we study the event of Shuangyuan bridge failure caused by 2009 Moarkot typhoon flood. We focus the discussion on the shape difference of the scour holes around the origional of bridge plies and the extended bridge piles. We also explore the suitability of using the empirical formula issued by the government (Ministry of transportation and communications R. O. C., 2009) We further study the the blocking effect from the drifting obstacles. Four scenarios are studied for understanding the effect from different stocking locations. The result shows no significant difference for difference stocking locations. However, adding the quantity of the drifting obstacles will increase the scouring rate and increase scour depth.
Overall speaking, Bingham Bi-viscous fluid model can be used to explore the development of scour hold, and avoid using too many empirical formulations and parameters. We also demonstrate that the present model can be widely used to predict the local scour depth induced by floods or tsunami bores.

[1] Assier-Rzadkiewicz, S., Mariotti, C., and Heinrich, P., "Numerical Simulation of Submarine Landslides and their hydraulic effects,” J. waterway, port, Coastal and Ocean Engineering, Vol. 123, No. 4, pp. 149-157, 1997.
[2] Balmforth, N. J. and Liu, J. J., “Roll waves in mud,” J. Fluid Mech., Vol. 519, pp. 33-54, 2004.
[3] Bird, R. B., Dai, G. C., and Yarusso, B. J., “The rheology and flow of viscoplastic materials,” Rev of Chemical Engrg., Vol. 1, No. 1, pp. 1-70, 1983.
[4] Borrero, J., Yeh H., Peterson, C., Chadha, R. K., Latha, G. and Katada, T., “Learning from earthquakes: The great Sumatra earthquake and Indian Ocean tsunami of December 26, 2004,” EERI Special Earthquake Report, March 2005.
[5] Brørs, B., “Numerical modeling of flow and scour at pipelines,” J. Hydr. Eng. Vol. 125, pp. 511–523, 1999.
[6] Chuang, M. H., “Developing a Two-way Coupled of Moving Solid Method for Solving Landslide Generated Tsunamis,” Master dissertation, National Central University, 2009.
[7] Clowes and Sons, Ltd., “Numerical simulation of pipeline local scour with lee-wake effects,” Int. J. Offshore Polar Eng., Vol. 10, pp. 195–199, 2000.
[8] Colicchio, G., Colagrossi, A., Greco, M. and Landrini, M., “Free-surface flow after a dam Break: a comparative study,” Ship Technol, Res. Vol. 49, No. 3, pp. 95–104, 2002.
[9] Danielsen, F., Sorensen, M. K., Olwig, M. F., Selvam, V., Parish, F., Burgess, N. D., Hiraishi, T., Karunagaran, V. M., Rasmussen, M. S., Hansen, L. B., Quarto A., Suryadiputra, N., “The Asian tsunami: a protective role for coastal vegetation,” Science Vol. 310, pp. 643, 2005.
[10] Dey, S., Bose, S.K. and Sastry, G.L.N. 1, “Clear water scour at circular piers: a model,” Journal of Hydraulic Engineering, Vol. 121, No. 12, pp. 869-876, 1995.
[11] Dey, S., and Barbhuiya, A. K., “Turbulent flow field in a scour hold at a semicircular abutment,” Can. J. Eng., Vol. 32, pp. 213-232, 2005.
[12] Ettema, R., Kirkil, G. and Muste, M., “Similitude of large-scale turbulence in experiments on local scour at cylinders,” Journal of Hydraulic Engineering, ASCE Vol. 132, No. 1, pp. 33–40, 2006.
[13] Fukui, Y., Nakamura, M., Shiraishi, H., and Sasaki, Y., “Hydraulic study on tsunami,” Coastal Engineering, Japan, Vol. 6, pp. 67–82, 1963a.
[14] Harlow, F. H. and Welch, J. E., 1965, “Numerical calculation of time-dependent viscous incompressible flow of fluid with a free surface,” Phys. Fluids, Vol. 8, pp. 2182-2189, 1965.
[15] Hirt, C. W., Nichols, B. D. and Romero, N. C., “SOLA-a numerical solution algorithm for transient fluid flows,” Los Alamos Scientific Laboratory, LA-582, pp. 1-50, 1975.
[16] Hirt, C. W. and Nichols, B. D., “Volume of Fluid (VOF) method for the dynamics of free surface boundaries,” J. Comput. Phys., pp.201-225, 1981.
[17] Huang, X. and Garcia, M. H., “A HerschalYHulkley model for mud flow down a slope,” J. Fluid Mech., Vol. 374, pp. 305-333, 1998.
[18] International Tsunami Survey Team (ITST) “The 26 December 2004 Indian Ocean tsunami:Initial findings from Sumatra,” Western Coastal and Marine Geology, U.S. Geological Survey, 2005, Internet Resource, http://walrus.wr.usgs.gov/tsunami/sumatra05/.
[19] Julien, P. Y. and Lan, Y., “Rheology of hyperconcentrations,” J. Hydraul. Eng. ASCE, Vol. 117, pp. 346-353, 1991.
[20] Kobayashi, N. and Lawrence, A. R., “Cross-shore sediment transport under breaking solitary waves,” J. Geophys. R. Vol. 109, C03047, 2004.
[21] Lacey, G., “Stable channels in alluvium,” Paper 4736, Minutes of the Proc., Institution of civil Engineers, William, 1930.
[22] Liu, P. L. F., Wu, T. R., Raichlen, Synolakis, C. E. and Borrero, J. C., “Runup and rundown generated by three-dimensional sliding masses,” J. Fluid Mech, pp. 107-144, 2005.
[23] Liu, K.F. and Mei, C.C., “Slow spreading of a sheet of Bingham fluid on an inclined plane,” J. Fluid Mech. Vol. 207, pp. 505-529, 1989.
[24] Liu, X. and García, M. H., “A three-dimensional numerical model with free water surface and mesh deformation for local sediment scour,” Journal of Waterway, Port, Coastal, and Ocean Engineering, Vol. 134, No. 4, pp. 203-217, 2008.
[25] Mao, Y., “The interaction between a pipeline and an erodible bed,” PhD dissertation, Institute of Hydrodynamics and Hydraulic Engineering, Technical University of Denmark, Lyngby, Denmark, 1986.
[26] Mascarenhas, A., Jayakumar, S., “An environmental perspective of the post-tsunami scenario along the coast of Tamil Nadu, India: Role of sand dunes and forests,” Journal of Environmental Management, Vol. 89, pp. 24–34, 2008.
[27] Melville, B. W., and Raudkivi, A. J., “Flow characteristics in local scour at bridge piers,” Journal of Hydraulic Research, Vol. 15, No. 4, pp. 373-380, 1977.
[28] Mei, C. C. and Yuhi, M., “Slow flow of a Bingham fluid in a shallow channel of finite width,” J. Fluid Mech., Vol. 431, pp. 135-159, 2001.
[29] Monaghan, J. J. and Kos, A., “Scott Russell’s wave generator,” Phys. Fluids, Vol. 12, No. 3, pp. 622–630, 2000.
[30] Monaghan, J. J., “Simulating free surface flows with SPH,” J. Comput. Phys, Vol. 110, pp. 399–406, 1994.
[31] Monaghan, J. J., “Smoothed particle hydrodynamics,” Rep. Prog. Phys, Vol. 68, pp. 1703–1759, 2005.
[32] Ng, C. and Mei, C.C., “Roll waves on a shallow layer of mud modeled as a power-law fluid,” J. Fluid Mech. Vol. 263, pp. 151-183, 1994.
[33] Noguchi, K., Sato, S. and Tanaka, S. “Large-scale experiments on tsunami overtopping and bed scour around coastal revetment,” Proc. Coastal Eng. Vol. 44, pp. 296–300, 1997.
[34] O’Brien, J.S. and Julien, P.Y., “On the importance of mudflow routing,” Proceedings of the 2nd International Conference on Debris Flow Hazards Mitigation, Taipei, Taiwan, Aug. 16-18, pp. 677-686, 2000.
[35] Olsen, N. R. B., and Melaaen, M. C., “Three-dimensional calculation of scour around cylinder,” J. Hydraul. Eng., Vol. 119, No. 9, pp. 1048-1054, 1993.
[36] Olsen, N. R. B., and Kjellesvig, H. M., “Three-dimensional numerical flow modeling for estimation of maximum local scour depth,” J. Hydraul. Res., Vol. 36, No. 4, pp. 579-590, 1998.
[37] Osher, S. and Sethian, J. A., “Fronts Propagating with Curvature Dependent Speed: Algorithms Based on Hamilton-Jacobi Formulations,” J. Comput. Phys. 79, 12-49, 1988.
[38] Richardson, J. E. and Panchang, V.G., “Three-dimensional simulation of scour-inducing flow at bridge piers,” J. Hydraul. Eng., Vol. 124, No. 5, pp. 530-540, 1998.
[39] Roulund, A., Sumer, B. M., Fredsoe, J. and Michelsen, J. “3D mathematical modeling of scour around a circular pile in current,” Proc., 7th Int. Symposium on River Sedimentation and 2nd Int. Symposium on Environmental Hydraulics 98, Hong Kong, 1999.
[40] Shao, S. and Lo, E. Y. M., “Incompressible SPH method for simulating Newtonian and non-Newtonian flows with a free surface,” Advances in Water Resources, Vol. 26, pp. 787-800, 2003.
[41] Sheppard, D.M., Odeh, M. and Glasser, T., “Large scale clear-water local pier scour experiments,” Journal of Hydraulic Engineering, ASCE Vol. 130, No. 10, pp. 957–963, 2004.
[42] Tanaka, N., Nandasena, N. A. K., Jinadasa, K. B. S. N., Sasaki, Y., Tanimoto, K., Mowjood, M. I. M., “Developing effective vegetation bioshield for tsunami protection,” Civil Engineering and Environmental Systems, Vol. 26, No. 2, pp. 163—180, 2009.
[43] Uda, T., Omata, A. and Yokoyama, Y., “Experimental study on tsunami run-up – the effects of coastal topography and structures against tsunami run-up,” Technical Note of Public Works Research Institute 2486, 1987.
[44] Wu, T. R., “A numerical study of three-dimensional breaking waves and turbulence effects,” PhD dissertation, Cornell University, 2004.
[45] Zhao, M., Cheng, L. and Zang, Z., “Experimental and numerical investigation of local scour around a submerged vertical circular cylinder in steady currents,” Coastal Engineering, Vol. 57, pp. 709–721, 2010.
[46] Zhao, Z. and Fernando, H. J. S., “Numerical simulation of scour around pipelines using an Euler-Euler coupled two phase model,” Env. Fluid Dyn., Vol. 7, No. 2, pp. 121-142, 2007.
[47] 交通部，「橋梁設計規範」，民國98年12月頒布。
[48] 吳祚任，「2011日本大海嘯之研究與省思」，土水會刊，第三十八卷第二期，2011年3月。
[49] 財團法人中華顧問工程司，「莫拉克颱風雙園大橋災害致災原因分析研究委託服務工作─期初報告書」，交通部公路總局第三區養護工程處，2010年4月。
[50] 國道45號線、5本の橋が落下，TBS系JNN，2011年3月12日。取自http://news.tbs.co.jp/newsi_sp/shinsai2011/tbs_newseye4671736.html
[51] 國道45號の橋流失6カ所に＝國交省，時事通信，2011年3月12日。取自http://news.nicovideo.jp/watch/nw40990