離岸風電場比陸域風電場更具優勢與潛能，但海上環境條件如颱風、波浪及海流等影響亦比陸域環境嚴苛。台灣海峽平均觀測風速高於國外已開發之離岸風電場，具有良好發展潛力，但台灣位於地震帶且多颱風侵襲，風險亦較其他離岸風電場高。為確保國內離岸風機的運轉安全及可靠性，分析與評估風機支撐結構是否得以承受場址特定環境條件(Site Specific Conditions)包括極端風況、海況及地震等衝擊，勢必成為未來離岸風電場建置前之重要工作項目。
本研究透過Splash3D模式進行三維數值模擬與分析，求解離岸風機在已知極端事件下受波浪及海流衝擊之承載力。Splash3D模式以美國國家實驗室（Los Alamos National Laboratory, LANL）所發展之Truchas模式為基礎進行研究與相關功能開發，模式可同時求解複雜且破碎之自由液面與結構物受力之動態過程，其以部份網格法（Partial-cell Treatment）於計算域中建構風機基礎，並用不規則網格包覆結構、配置數值壓力計以求結構物整體受力分佈；以投影法（Projection Method）求解大渦模擬（Large Eddy Simulation）方程式；以流體體積法（Volume of Fluid, VOF），描述表面碎波行為。
模式中為能夠適切模擬季風以至於颱風之波浪特性，以頻散關係式為基礎設計之內造波機模組（Internal Wave Maker Module）可依據自氣象局取得的單一波浪波譜觀測資料生成JONSWAP理論頻譜，並以此做為輸入條件產生更符合實際海況的亂數波浪。研究中考量波浪與海流之加總效應，比較相同示性波高條件下，規則波與不規則波衝擊結構時其受力分布，以期提供未來離岸風電場建置之參考。

In this research, the three-dimensional numerical simulation and analysis for solving the dynamic loads from waves and currents on the offshore wind turbines is performed. Scenarios focused on the extreme weather conditions. During the typhoon event, the wind-driven storm waves and currents have to be considered while solving the dynamic load on the structures. The Splash3D model was adopted to perform the simulation of the interaction between breaking waves and structures. The core of the Splash3D model is the Truchas model which was developed by Los Alamos National Laboratory (LANL) and featured as high accuracy. Splash3D is capable of solving the dynamic process for the interaction between the structure and fluids with complex breaking free-surface. In order to adequately simulate the waves under monsoon or typhoon, a new wave generation module based on the dispersion relationship is developed. This wave-maker module was used to generate regular waves, irregular waves, and breaking waves under the extreme weather condition. The module was used to simulate the synthetic effect under the effects of waves and currents for obtaining the force distribution on the foundation of the offshore wind turbine.

論文指導教授推薦書 III
論文口試委員審定書 IV
中文摘要 V
英文摘要 VI
目錄 VII
圖目錄 IX
表目錄 XVI
第1章 緒論 1
1-1 研究動機 1
第2章 模式介紹與研究方法 7
2-1 模式介紹 7
2-1-1 控制方程式（Governing Equations） 7
2-1-2 有限體積法（Finite Volume Method） 9
2-1-3 流體體積法與PLIC法（Piecewise Linear Interface Calculation） 10
2-1-4 部分網格法（Partial-Cell Method） 13
2-1-5 大渦模擬（Large Eddy Simulation , LES） 14
2-2 依數值情境模擬結果輸出結構物表面受力 17
2-2-1 物體浮力驗證 17
2-2-2 潰壩波衝擊結構物之受力驗證 20
2-3 依據研究區域觀測資料決定三維數值模擬造波輸入條件 24
2-3-1 自氣象局新竹浮標 (46757B) 觀測資料選用適當之波譜 24
2-3-2 依據選用之波譜觀測資料產生與其相符之JONSWAP理論波譜 28
2-3-3 以波譜能量密度函數曲線求得自由水面隨時間變化之序列 31
第3章 模擬結果與討論 37
3-1 規則波對於5MW參考離岸風機桁架式(Lattice)支撐結構，在不同波向衝擊下之模擬及受力分析 37
3-1-1 情境模擬之波浪條件及計算域設置 37
3-1-2 規則波與流加總效應之二維流況情境模擬 39
3-1-3 三維規則波與定速海流對單圓柱之模擬及受力分析 53
3-1-4 三維規則波與定速海流對結構物在攻角0度衝擊下之模擬及受力分析 64
3-1-5 三維規則波與定速海流對結構物在攻角45度衝擊下之情境模擬及受力分析 75
3-2 不規則波對於5MW參考離岸風機桁架式(Lattice)支撐結構，在不同波向衝擊下之模擬及受力分析 86
3-2-1 情境模擬之波浪條件及計算域設置 86
3-2-2 不規則波與定速海流加總效應之二維流況情境模擬 89
3-2-3 三維不規則波與定速海流對單圓柱之模擬及受力分析 101
3-2-4 三維不規則波與定速海流對5MW參考離岸風機桁架式(Lattice)支撐結構在攻角0度衝擊下之模擬及受力分析 112
3-2-5 三維不規則波與定速海流對5MW參考離岸風機桁架式(Lattice)支撐結構在攻角45度衝擊下之模擬及受力分析 123
第4章 結論與建議未來之研究方向 134
附錄：SPLASH3D沖刷模式應用 136
參考文獻 148
口試書面答覆表 152
 
圖目錄
圖 2-1 1 流體體積法之體積分率示意圖 12
圖 2-1 2 自由液面重建示意圖（Wu, 2004） 13
圖 2-2 1模擬設置圖 18
圖 2-2 2壓力積分點布置全貌與局部放大圖 18
圖 2-2 3擷取受力數據示意圖 19
圖 2-2 4 Wu (2004)、 Hu, et al. (2012) 潰壩波衝擊方柱之數值模擬設置圖 21
圖 2-2 5 Wu (2004)、Hu, et al. (2012) 潰壩波衝擊方柱之數值模擬設置全貌 21
圖 2-2 6 壓力積分點於方柱上之分佈情形 22
圖 2-2 7（上）Wu (2004)潰壩波衝擊方柱之文獻中實驗（虛線）與模擬（藍色實線）比對結果與（下）壓力積分點擷取受力之結果（綠色實線）比較 22
圖 2-2 8於Splash 3D中建構之風機基礎 23
圖 2-2 9數值壓力計分布示意圖 23
圖 2-3 1 2013年1月3日上午7點之波譜觀測資料 26
圖 2-3 2 2013年2月8日上午4點之波譜觀測資料 27
圖 2-3 3 2014年5月5日上午7點之波譜觀測資料 27
圖 2-3 4 2013年7月13日上午6點之波譜觀測資料 28
圖 2-3 5 2013年7月13日上午6點之波譜觀測資料 (角頻率域) 30
圖 2-3 6 選用波譜觀測資料(2013/07/13/06)與JONSWAP理論波譜比對結果 31
圖 2-3 7 JONSWAP理論頻譜能量密度分佈曲線下各頻率成分所占曲線下面積示意圖 34
圖 2-3 8 依據JONSWAP理論頻譜生成之自由水面隨時間變化之序列圖。 34
圖 3-1 1 結構物與水流方向之試驗條件關係圖 38
圖 3-1 2 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖 41
圖 3-1 3 數值水位計所在位置之波高歷線圖。 41
圖 3-1 4 數值水位計所在位置之波高歷線圖。 42
圖 3-1 5 數值水位計所在位置之波高歷線圖。 43
圖 3-1 6 數值水位計所在位置之波高歷線圖。 44
圖 3-1 7 數值水位計所在位置之波高歷線圖。 45
圖 3-1 8 數值水位計所在位置之波高歷線圖。 46
圖 3-1 9 數值水位計所在位置之波高歷線圖。 47
圖 3-1 10 二維模擬之規則波流場沿水流方向流速分布圖（第90至第99秒）。 48
圖 3-1 11 二維模擬之規則波流場垂直水流方向流速分布圖（第90至第99秒）。 49
圖 3-1 12 二維模擬之規則波流場流速絕對值分布圖（第90至第99秒）。 50
圖 3-1 13 二維模擬之規則波流場總壓分布圖（第90至第99秒）。 51
圖 3-1 14 二維模擬之規則波流場動壓分布圖（第90至第99秒）。 52
圖 3-1 15 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖。 54
圖 3-1 16 數值壓力計分布示意圖。藍點為數值壓力計位置。 54
圖 3-1 17 三維規則波與定速海流對單一圓柱衝擊之水平方向流速分布剖面圖。 57
圖 3-1 18 三維規則波與定速海流對單一圓柱衝擊之垂直方向流速分布剖面圖。 58
圖 3-1 19 三維規則波與定速海流對單一圓柱衝擊之流速絕對值分布剖面圖。 59
圖 3-1 20 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之 三維流線分布圖（103至113秒）。 60
圖 3-1 21 三維規則波與定速海流對單一圓柱衝擊之流場總壓分布剖面圖。 61
圖 3-1 22 三維規則波與定速海流對單一圓柱衝擊之流場動壓分布剖面圖。 62
圖 3-1 23三維規則波與定速海流對單一圓柱衝擊之波高歷線圖（上）、 平行水流方向受力歷線圖（中）及所受總力歷線圖（下） 63
圖 3-1 24 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖。 65
圖 3-1 25 數值壓力計分布示意圖。藍點為數值壓力計位置。 65
圖 3-1 26 桁架式支撐結構受攻角0度之規則波與流衝擊之水平方向流速分布剖面圖。 68
圖 3-1 27 桁架式支撐結構受攻角0度之規則波與流衝擊之垂直方向流速分布剖面圖。 69
圖 3-1 28 桁架式支撐結構受攻角0度之規則波與流衝擊之流速絕對值分布剖面圖。 70
圖 3-1 29 桁架式支撐結構受攻角0度之規則波與流衝擊之三維流線分布圖。 71
圖 3-1 30 桁架式支撐結構受攻角0度之規則波與流衝擊之流場總壓分布剖面圖。 72
圖 3-1 31 桁架式支撐結構受攻角0度之規則波與流衝擊之流場動壓分布剖面圖。 73
圖 3-1 32桁架式支撐結構受攻角0度之規則波與流加總效應衝擊之波高歷線圖（上）、平行水流方向受力歷線圖（中）及結構物所受總力歷線圖（下）。 74
圖 3-1 33 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖。 76
圖 3-1 34 數值壓力計分布示意圖。藍點為數值壓力計位置。 76
圖 3-1 35 桁架式支撐結構受攻角45度之規則波與定速海流衝擊之 水平方向流速分布剖面圖。 79
圖 3-1 36 桁架式支撐結構受攻角45度之規則波與定速海流衝擊之 垂直方向流速分布剖面圖。 80
圖 3-1 37 桁架式支撐結構受攻角45度之規則波與定速海流衝擊之 流速絕對值分布剖面圖。 81
圖 3-1 38 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之 三維流線分布圖（103至113秒）。 82
圖 3-1 39 桁架式支撐結構受攻角45度之規則波與定速海流衝擊之 流場總壓分布剖面圖。 83
圖 3-1 40 桁架式支撐結構受攻角45度之規則波與定速海流衝擊之 流場動壓分布剖面圖。 84
圖 3-1 41 桁架式支撐結構受攻角45度之規則波與定速海流衝擊之 波高歷線圖（上）、平行水流方向受力歷線圖（中）及結構物所受總力歷線圖（下） 85
圖 3-2 1 2013年7月13日上午6點之波譜觀測資料與JONSWAP理論波譜比對結果 87
圖 3-2 2 依據JONSWAP理論頻譜生成之自由水面隨時間變化之序列圖。 88
圖 3-2 3數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖 92
圖 3-2 4 0-1號數值水位計所在位置之波高歷線圖。 92
圖 3-2 5 2-4號數值水位計所在位置之波高歷線圖。 93
圖 3-2 6 5-6號數值水位計所在位置之波高歷線圖。 94
圖 3-2 7 7-9號數值水位計與輸入條件之波高歷線圖比較結果。 95
圖 3-2 8二維模擬之不規則波與定速海流流場之沿水流方向流速分布圖。 96
圖 3-2 9 二維模擬之不規則波與定速海流流場之垂直水流方向流速分布圖。 97
圖 3-2 10 二維模擬之不規則波與定速海流流場之流速絕對值分布圖。 98
圖 3-2 11 二維模擬之不規則波與定速海流流場之總壓分布圖。 99
圖 3-2 12二維模擬之不規則波與定速海流流場之動壓分布圖。 100
圖 3-2 13 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖。 102
圖 3-2 14 數值壓力計分布示意圖。藍點為數值壓力計位置。 102
圖 3-2 15 圓柱受不規則波與定速海流加總效應衝擊之水平方向流速分布剖面圖。 105
圖 3-2 16 圓柱受不規則波與定速海流加總效應衝擊之 垂直方向流速分布剖面圖。 106
圖 3-2 17 圓柱受不規則波與定速海流加總效應衝擊之流速絕對值分布剖面圖。 107
圖 3-2 18 圓柱受不規則波與定速海流加總效應衝擊之三維流線分布圖。 108
圖 3-2 19 圓柱受不規則波與定速海流加總效應衝擊之流場總壓分布剖面圖。 109
圖 3-2 20 圓柱受不規則波與定速海流加總效應衝擊之流場動壓分布剖面圖。 110
圖 3-2 21 圓柱受不規則波與定速海流加總效應衝擊之波高歷線圖（上）、平行水流方向受力歷線圖（中）及結構物所受總力歷線圖（下） 111
圖 3-2 22 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖。 113
圖 3-2 23 數值壓力計分布示意圖。藍點為數值壓力計位置。 113
圖 3-2 24 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之水平方向流速分布剖面圖。 116
圖 3-2 25 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之垂直方向流速分布剖面圖。 117
圖 3-2 26 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之流速絕對值分布剖面圖。 118
圖 3-2 27 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之 三維流線分布圖。 119
圖 3-2 28 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之流場總壓分布剖面圖。 120
圖 3-2 29 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之流場動壓分布剖面圖。 121
圖 3-2 30 桁架式支撐結構受攻角0度之不規則波與定速海流加總效應衝擊之波高歷線圖（上）、平行水流方向受力歷線圖（中）及結構物所受總力歷線圖（下） 122
圖 3-2 31 數值計算域之設置圖。上圖為鳥瞰圖，下圖為側視圖。 124
圖 3-2 32 數值壓力計分布示意圖。藍點為數值壓力計位置。 124
圖 3-2 33 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之水平方向流速分布剖面圖。 127
圖 3-2 34 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之垂直方向流速分布剖面圖。 128
圖 3-2 35 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之流速絕對值分布剖面圖。 129
圖 3-2 36 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之 三維流線分布圖。 130
圖 3-2 37 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之流場總壓分布剖面圖。 131
圖 3-2 38 桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之流場動壓分布剖面圖。 132
圖 3-2 39桁架式支撐結構受攻角45度之不規則波與定速海流加總效應衝擊之 波高歷線圖（上）、平行水流方向受力歷線圖（中）及結構物所受總力歷線圖（下） 133
 
表目錄
表 2 1 新竹浮標每月波高統計表 (2000年～2014年) 25
表 2 2 JONSWAP理論波譜個頻率分量參數 (1) 35
表 2 3 JONSWAP理論波譜個頻率分量參數 (2) 36
表 3 1模式中採用之波浪條件 38
表 3 2 模式計算條件 87

[1] Morison, J. R., Johnson, J. W., & Schaaf, S. A. (1950). The force exerted by surface waves on piles. Journal of Petroleum Technology, 2(05), 149-154.
[2] Wienke, J. and Oumeraci, H., (2005). Breaking Wave Impact Force on a Vertical and Inclined Slender Pile - Theoretical and Large-Scale Model Investigations. Coastal Engineering, 52:435 462
[3] Aune, L. (2011). Forces from plunging breaking waves on a truss structure (Doctoral dissertation, Master thesis submitted June 2011, Norwegian University of Science and Technology, Department of Civil and Transport Engineering, Trondheim, Norway).
[4] Aashamar, M. (2012). Wave slamming forces on truss support structures for wind turbines. Master thesis, submitted June 12, 2012, Norwegian University of Science and Technology, Department of Civil and Transport Engineering, Trondheim, Norway.
[5] Goda, Y., Haranaka, S. and Kitahata, M. (1966). Study of impulsive breaking wave forces on piles. Report of Port and Harbor Research Institute, Japan, Vol. 5. No. 6, pp. 1 – 30 (in Japanese). Concept also in English language in: Watanabe, A. and Horikawa, K. (1974). Breaking wave forces on large diameter cell. Proc. of 14th International Conference on Coastal Engineering, Ch. 102, pp. 1741 – 1760.
[6] Navaratnam, C. U., Tørum, A., & Arntsen, Ø. A. (2013). Preliminary analysis of wave slamming force response data from tests on a truss structure in large wave flume, Hannover, Germany. Department of Civil and Transport Engineering, NTNU, Trondheim, Norway.
[7] Arntsen, Ø. A., & Gudmestad, O. T. (2014). Wave slamming forces on truss structures in shallow water. In Proceedings of the HYDRALAB IV Joint User Meeting, Lisbon.
http://www.hydralab.eu/flashdrive/documents/TA-FZK-05.pdf
[8] Christensen, E. D., Bredmose, H., & Hansen, E. A. (2005). Extreme wave forces and wave run-up on offshore wind turbine foundations. Proceedings of Copenhagen Offshore Wind, 1-10.U.S. Army, 1984. Shore Protection Manual, Vol. 1. Coastal Engineering Research Center.
[9] Bredmose, H., & Jacobsen, N. G. (2010). Breaking wave impacts on offshore wind turbine foundations: focused wave groups and CFD. In ASME 2010 29th International Conference on Ocean, Offshore and Arctic Engineering (pp. 397-404). American Society of Mechanical Engineers.
[10] Marino, E., Borri, C., & Peil, U. (2011). A fully nonlinear wave model to account for breaking wave impact loads on offshore wind turbines. Journal of Wind Engineering and Industrial Aerodynamics, 99(4), 483-490.
[11] Rausa, I. E., Muskulus, M., Arntsen, Ø. A., & Wåsjø, K. (2015). Characterization of wave slamming forces for a truss structure within the framework of the WaveSlam project. Energy Procedia, 80, 276-283.
[12] Wu, T. R. (2004). A numerical study of three-dimensional breaking waves and turbulence effects. PhD dissertation, Cornell University
[13] Liu, P. F., Wu, T. R., Raichlen, F., Synolakis, C. E., & Borrero, J. C. (2005). Runup and rundown generated by three-dimensional sliding masses. Journal of fluid Mechanics, 536, 107-144.
[14] Hirt, C. W., & Nichols, B. D. (1981). Volume of fluid (VOF) method for the dynamics of free boundaries. Journal of computational physics, 39(1), 201-225.
[15] Kothe, D. B., Rider W. J., Mosso, S. J., Brock, J. S., and Hochstein, J. I. (1996). Volume tracking of interfaces having surface tension in two and three dimensions. Technical Report, AIAA 96-0859.
[16] Rider, W. J. and Kothe, D. B. (1998). Reconstructing Volume Tracking. J. Comp.Phys., 141, 112-152.
[17] Deardorff, J. W. (1970). A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers. Journal of Fluid Mechanics, 41(02), 453-480.
[18] Leonard, A. (1975). Energy cascade in large-eddy simulations of turbulent fluid flows. Advances in geophysics, 18, 237-248.
[19] Smagorinsky, J. (1963). General circulation experiments with the primitive equations: I. the basic experiment*. Monthly weather review, 91(3), 99-164.
[20] Lin, P., & Li, C. W. (2003). Wave–current interaction with a vertical square cylinder. Ocean Engineering, 30(7), 855-876.
[21] Cabot, W. & Moin, P. (2000). Approximate wall boundary conditions in the large-eddy simulation of high Reynolds number flow. Flow Turb. Combust. 63, 269-291.
[22] Hu, K. C., Hsiao, S. C., Hwung, H. H., & Wu, T. R. (2012). Three-dimensional numerical modeling of the interaction of dam-break waves and porous media. Advances in Water Resources, 47, 14-30.
[23] Vorpahl, F., Popko, W., & Kaufer, D. (2011). Description of a basic model of the” UpWind reference jacket” for code comparison in the OC4 project under IEA Wind Annex XXX. Fraunhofer Institute for Wind Energy and Energy System Technology (IWES), 4, 1-14.
[24] Hasselmann, K., Barnett, T. P., Bouws, E., Carlson, H., Cartwright, D. E., Enke, K., ... & Meerburg, A. (1973). Measurements of wind-wave growth and swell decay during the Joint North Sea Wave Project (JONSWAP). Deutches Hydrographisches Institut.
[25] WAFO Group. (2000). A Matlab toolbox for analysis of random waves and loads. Lund University, Lund Institute of Technology, Centre for Mathematic Sciences, Mathematical Statistics.
[26] Lin, P., & Liu, P. L. F. (1999). Internal wave-maker for Navier-Stokes equations models. Journal of waterway, port, coastal, and ocean engineering, 125(4), 207-215.
[27] Sumer, B. M., & Fredsøe, J. (2002). The mechanics of scour in the marine environment. World Scientific.
[28] Herbich, J. B. (1981). Scour around pipelines and other objects. Offshore pipeline design elements.
[29] Herbich, J. B. (1984). Seafloor scour: Design guidelines for ocean-founded structures (Vol. 4). Marcel Dekker Inc.
[30] Whitehouse, R. (1998). Scour at marine structures: A manual for practical applications. Thomas Telford.
[31] Dey, S. (1999). Time-variation of scour in the vicinity of circular piers. Proceedings of the Institution of Civil Engineers-Water Maritime and Energy, 136(2), 67-75.
[32] Miller Jr, W., & Sheppard, D. M. (2002). Time rate of local scour at a circular pile. In First International Conference on Scour of Foundations.
[33] Bingham, E. C. (1922). Fluidity and plasticity (Vol. 2). McGraw-Hill Book Compny, Incorporated.
[34] Herschel, W. H., & Bulkley, R. (1926). Konsistenzmessungen von gummi-benzollösungen. Colloid & Polymer Science, 39(4), 291-300.
[35] Bagnold, R. A. (1954, August). Experiments on a gravity-free dispersion of large solid spheres in a Newtonian fluid under shear. In Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences (Vol. 225, No. 1160, pp. 49-63). The Royal Society.
[36] Rzadkiewicz, S. A., Mariotti, C., & Heinrich, P. (1997). Numerical simulation of submarine landslides and their hydraulic effects. Journal of Waterway, Port, Coastal, and Ocean Engineering, 123(4), 149-157.
[37] Liu, K. F., & Mei, C. C. (1989). Slow spreading of a sheet of Bingham fluid on an inclined plane. Journal of fluid mechanics, 207, 505-529.
[38] 李宇弘，2011，非連續賓漢流變模式之參數訂定與應用The Parameter Determination and Model Application of Discontinuous Bi-viscous Model，國立中央大學，碩士論文。
[39] Hadush, S., Yashima, A., Uzuoka, R., Moriguchi, S., & Sawada, K. (2001). Liquefaction induced lateral spread analysis using the CIP method. Computers and Geotechnics, 28(8), 549-574.
[40] Chatterjee, S. S., Ghosh, S. N., & Chatterjee, M. (1994). Local scour due to submerged horizontal jet. Journal of Hydraulic Engineering, 120(8), 973-992.