2D斷面消波數值水槽VB2026年版      2022.08.10   更新日期:2026/02/02

一、目的

    筆者於1990年代利用邊界元素法,設計出一系列消波式數值2D造波水槽如下述參考文獻。2011年退休後利用Java Applet將之可視化,花費約1年完成並上網供分享。然學生們認為這些研究屬學術對工程實務助益不大,建議編輯工程實務相關手冊。個人此時亦有同感學術成果只是自嗨,當時這些研究主要目的為計劃申請及職務上義務,對現職教師都是壓力,曾有任教職的學生跟我說,研究是件痛苦不得已的事。有沒有這些研究社會不會有實質改變,因此將手頭上的分析軟體全部丟到垃圾筒,專心編輯港灣海岸工程實務相關手冊,歷經10個歲月在2021年完成港灣海岸工程手冊,供學子及甫步入社會的菜鳥參考。2022年初Microsoft宣布中止Java Applet支援,因此辛苦完成的JAVA波浪公式集化為烏有,心有不甘於4月重新學習Visual Studio的Visual Basic,編輯波浪公式集首先將各種波浪現象製成可視化動畫,例如可能可表現瘋狗浪形成的兩波斜向交會等。今首先完成通用的消波式數值2D造波水槽,水槽內配置半潛式、潛水式結構物時的運動,或海底配置消波潛堤時的消波效果等屬客製化動畫將陸續上網。除數值2D造波水槽外,蛇型多方向不規則3D造波水槽、海灘漂砂變形、船體運勤、多本船真時操船模擬等,亦會陸續上網。

     更新2D斷面消波數值水槽,追加補充說明。2026/01/11

二、 操作說明

1. 本程式係以執行檔型式供造訪者下載,在造訪者電腦執行。
2. 本程式可為直立壁或斜坡空水槽,以事先對各種波浪設定「預設值」供參考但不保証一定可執行,必要時可合理調整各項數據至  可執行。
3. 設定步驟如下:
   1.選定造波型式:選定波型打✓後按「選定波型」,再按「確定」。
   2.設定水槽:可調整水槽長度、水槽水深、有無設置消波段(sponge zone)等後按「設定水槽」。
   3. 波浪設定:可調整造波波高、週期、時間間隔等後按「波浪設定」,造波時間間隔針對不同波型預設不同值,可適度調整。
   4.節點設定:可調整各邊界節點數後按「節點設定」。
   5.全部設定完成後按「顯示波動」

4. 所有水深均轉換為單位水深1m,其他物理項依Froude相似率配合轉換。

補充說明(2026/01/12)

1.  本軟體為專業性波浪分析工具,討論各種波浪在不同狀況的波動,可能形成重複波或形成碎波,由於達臨界碎波後流體運動不再具有速度勢,爾後計算不具實值物理意義,必要以專業終止計算。

2.  本波浪分析工具考量自由水面的強非線性運動條件,數值計算屬較不穩定狀態,影響因素除有水深、波高、週期及地形等自然因素外,尚計算元素大小及計算時間隔等數值計算因素。自然因素可控、數值計算因素雖有各種方法可調整其相互關係達穩定,但實際有因難,必要依測試覓出可適用範圍。

4.  由於水槽水深設定為單位水深1m,實際海面經常發生的波浪週期約在7~20秒範圍,依相似律換算,實驗相對週期則約在2~5秒間。

5.  對具有週期性的波,設定元素大小時為滿足週期性波的物理特性,每波長必要配置至少16個元素。由於考量各種波浪均能適用,元素大小設定不採用傳統的波浪波長,而是以水槽水深為標準,預定值為每單位水深設定4個元素,即自由水面元素數 n2 = 水槽長度 x 6 + 1,海底面元素數 n4 = (水槽長度-斜面長度) x 2 + 1,可視必要調整。

6.  對孤立波,元素大小配合波高適度調整。因孤立波造波速度快,會使造波板超越元素導致當機,造波時間間隔須較短。

7.  對不規則波,元素大小設定,對不規則波群的有義週期的波長,至少取16個元素為宜,。

8.  側壁為斜面時,元素大小設定約與單位水深相等即可。

9.  元素大小大致依上述原則設定,最後必要設定造波時間間隔,造波時間間隔理論上與波特性有關,但元素大小佔很大因素,實用上造波時間間隔設定0.01~0.02,預設值為造波週期的1/200。對週期波,起始造波時數值計算不穩定,前400次起始造波的造波間隔設定為造波週期的1/400。

10. 斜面越緩,造波時間間隔宜越短。

11. 水面元素n2>100時、間隔2點顯示,n2>200時、間隔4點顯示。

12. 目前本軟體僅顯示空水槽各種波浪的時間歷時,可發展至水面、水中或海底地形變化等狀況,邊界元素法程式僅需作部分修正,主要為結構物座標設定即可應用,其方法述於邊界元素法在海岸工程應用。由於VB的sinx檔無法上網,只能以下述壓縮檔上綢,使用者下載exe檔執行。有意開發者可連絡本人提供原始檔。

13.本軟體採自動配置,但可自行調整修正。

測試結果

   水槽造波水深設定為單位水深(1m)、水槽長度預設為50m(可調整),不規則波為100m。教學示範時為縮短時程,對垂直側壁水槽長度可調整為30m,以了解垂直側壁的影響。

 A.規則波(週期波)

1.垂直側壁

a.測試完成範圍:涵蓋實際海上6~14秒波浪。

b.換算成單位水深d(1m)時,造波波高:0.2d以下,造波週期:3秒以下,依預設值即可。

c.超出測試範圍,可變更水槽總長使其超出6個測試波長以上(使造波達穩定),或許可得期待模擬。執行1次演算,取得模型波長重新執行鍵入適當水槽總長後重新執行。模擬過程因自由水面元素增加,耗時較長。2026/01/31

   本模擬可觀測,從造波板起始造波,此時波形受造波板影響呈不安定,隨著造波的持續波前進,為取得完整波形,水槽必要有足夠長度,至少超出6個測試波長以上為宜。未達垂直側壁前以進行波形式傳播,到達側壁受側壁反射,開始從側壁反射逐一形成重複波,此時其他波尚呈進行波。如何從進行波順利轉換形成重複波是模擬上的關鍵, 一般進行波波長以8個元素就可表達出其特性,重複波則需16個元素。但是在轉換過程需要更多元素,以24個元素以上為適。

   進行波到達側壁受側壁反射後,依序逐一形成重複波,此時在重複波逐漸形成過程,某時刻該部分水面會呈靜止,表示重複波特性出現。波向造波板方向反射,當反射波到達造波板時,因整體水槽水運動呈重複波狀態。模擬大功告成。

2.陡坡側壁(坡度 1:1~1:5)

   水槽長度設定為30m,完成測試範圍為,坡度 1:1~1:5,造波週期2~3秒,造波波高0.3d以下。

3.緩坡度海底(坡度 1:10~1:30)

 

 

B.孤立波

1.  自由水面節點數=101,側壁及海底節點數適度調整。

2.  造波時間間隔,起始造波時造波時間間隔設定為0.02(可調整)。為縮短計算,起波達安定時,計算時間間隔增倍,當孤立波達水槽中段時考量波與水槽側壁(斜面)相互影響,計算時間間隔縮短為起始造波時造波時間間隔的1/2。

三、 操作方法

 1. 下載 壓縮檔解壓縮

 2. 解壓縮至適當資料夾(自行選定位置及名稱) 

 3. 執行 dissipativeChannel.exe

 4. 依序輸入波浪條件後,確定

 

四、參考文獻

1.周宗仁、石瑞祥,1996,斜坡上考量碎波效應波浪變形之研究,港灣技術第11卷第1期。

2.  Chou, Chung-Ren, Kwan Ouyang, 1999, The deformation of solitary waves on steep slopes, Journal of the Chinese Institute of Engineers, Vol.22, No 6,    pp.805-812.

3. Chou, Chung-Ren, Ruey-Syan Shih and John Z Yim, 2002, Optimizing deployment of sponge zone on numerical wave channel, Journal of the Chinese Institute of   Engineers.  Vol. 25,No.2, pp.147-156.

4. Chou, Chung-Ren, Ruey-Syan Shih, 1996, Numerical generation and propagation of periodical wave in time domain, Coastal Engineering in Japan, Vol.39, No. 2.

5. Chou, Chung-Ren, Ruey-Syan Shih, 1996, Generation and Deformation of Solitary Waves, China Ocean Engineering, Vol.10 No. 4.

6. Chou, C. R., and K. Ouyang, 1999, Breaking of Solitary Waves on Uniform Slopes, China Ocean Engineering, 13(4).(EI, SCI)

7 Chou, Chung-ren, John Z. Yim , R.S. Shih, 2003,Numerical study on the breaking criterion for solitary waves,China Ocean    Engineering .Vol.17,NO.4,pp.589-604. (SCI)

8. Shih, R.S., C.R. Chou, and John Z. Yim, 2004, Numerical investigation on the generation an propagation of irregular waves in a two-dimensional wave tank, China ocean eng., Vol. 18, No, 4, pp. 551-566. (SCI)

9. Shih Ruey-Syan ,Chungren Chou, John Z. Yim, 2006,"Modelling of random waves over submerged breakwaters and its application to reflection estimation",China Ocean Engineering, Vol. 20, No. 2, pp. 243-258. (SCI)

10.石瑞祥、周宗仁、尹 彰 ,2004,二維數值造波水槽不規則波之數值研究,中國燕山大學學報,,第28卷,第2期。

11.周宗仁、石瑞祥、韓文育 ,2004,斜坡上孤立波碎波指針之數值研究,中國燕山大學學報,,第28卷,第2期。

13. Chou, C.R., John Z. Yim, R.S. Shih and K. Ouyang, 2001, Numerical study on the breaking of solitary wave on slopes, Boundary Element Communications, Vol.12, No.4, pp.3-12.  (SCI)

14.Shih, Ruey-Syan*, Weng Wen-kai and Chou Chung-ren, 2013, “Numerical Modeling of Wave Field Around Multiple Submerged Breakwaters", Marine Science, Vol.3, No.3, pp.65-78.

15.周宗仁、石瑞祥、尹彰,2000, 消波式二維數值造波水槽的開發研究,第22屆海洋工程研討會論文集,pp.119-126

16.周宗仁、石瑞祥、尹彰,歐陽寬,2001, 孤立波碎波指標之數值研究, 23屆海洋工程研討會論文集。

17.  Chou, Chung-Ren, and Kwan Ouyang, 1998, Numerical Study of Solitary Wave on A Mild Slope, ISOPE.