資源描述：
《基于緩坡方程在島礁地形上波浪破碎的模擬研究_方亞冰》由會員上傳分享，免費在線閱讀，更多相關內容在教育資源-天天文庫。
1、基于緩坡方程在島礁地形上波浪破碎的模擬研究_方亞冰第36卷第4期JournalofWaterwayandHarbor水道港口Vol.36No.4基于緩坡方程在島礁地形上波浪破碎的模擬研究方亞冰，柳淑學，李金宣，劉思（大連理工大學海岸和近海工程國家重點實驗室，大連116024）摘要：波浪從深海傳至近岸島礁時，一般需要經(jīng)過落差較大的礁坪邊緣，水深急劇變化，導致波浪在傳播過程中發(fā)生破碎，因此準確模擬波浪的破碎過程以及破碎后的波高大小，對于島礁海岸工程建設具有重要的意義。緩坡方程是描述近岸波浪傳播變形較好的數(shù)值模型之一，文章在采用自適
2、應有限元求解緩坡方程所建立的數(shù)值計算模型的基礎上，引入描述波浪破碎的模型，建立可以描述波浪破碎影響的近岸波浪數(shù)值模型。基于二維島礁地形上的波浪實驗，比較分析了4種不同的波浪破碎能量損失因子，給出了適合于島礁地形條件下波浪傳播破碎模擬的模型。關鍵詞：二維島礁地形；緩坡方程；波浪破碎中圖分類號：TV142；TV131.6文獻標識碼：A文章編號：1005-8443（2015）04-0290-07中國海域幅員遼闊、珊瑚島礁眾多，波浪從深海傳至近岸島礁時，一般需經(jīng)過落差較大的礁坪邊緣。而由于水深的急劇變化，導致與在———————————
3、————————————————————————————————————大多數(shù)緩變地形上波浪傳播相比，其規(guī)律有所不同。因此建立起能夠模擬波浪在島礁地形上傳播的數(shù)值計算模型，對波浪變形進行較為準確的預報是十分有必要的?？紤]到在波浪傳播的過程中，會受到折射、繞射、反射等因素的影響，很多學者在二維橢圓型緩坡方程的基礎上，建立起了多個數(shù)值模型（例如：Mei［1］；TsayandLiu［2］；ChenandHouston［3］；Mattioli［4］），并運用到了實際工程當中。但是，大部分模型都沒有考慮到如島礁地形這種水深差別較大情況下
4、波浪破碎的影響。本文在Berkhoff［5］發(fā)展的緩坡方程的基礎上，基于自適應有限元方法來求解緩坡方程［6］，著重考慮波浪破碎因素的影響，針對不同學者提出的4種不同的波浪破碎能量損失因子，對二維島礁地形條件下波浪的傳播進行模擬，通過與實驗結果進行對比、分析，建立起適用于島礁地形下波浪傳播的數(shù)值計算模型。11.1數(shù)值計算模型控制方程考慮到波浪在傳播過程中波浪破碎的效應，有些學者認為可以把波浪破碎作為一項耗散項加入到Berk?hoff的原始緩坡方程當中［7-9］。?·(F·??)+(ω2G+iCgωγ)?=0F=CCg，G=（1）
5、Cg（2）Cg(x,y)=?ω?k為群速度，γ為波浪破碎能量損失因子,可有多種計算模型，將其中C（x,y）=ω/k為相速度，式中在下文給出。k為波數(shù)，與波浪頻率滿足如下色散關系ω2=gktanh(kh)（3）收稿日期：2014-12-05；修回日期：2015-01-04———————————————————————————————————————————————基金項目：國家重點基礎研究發(fā)展（973）計劃資助項目（2013CB036101，2011CB013703）；國家自然科學基金創(chuàng)新研究群體基金（51221961）作者簡介
6、：方亞冰（1989-），男，湖南省郴州市桂東縣人，碩士研究生，主要從事近岸波浪傳播特性的模擬研究。Biography：FANGYa?bing(1989-),male,masterstudent.方亞冰，等基于緩坡方程在島礁地形上波浪破碎的模擬研究291式中：h為水深。1.2邊界條件在近岸波浪傳播模擬過程中，一般需考慮如下邊界條件。=0（4）（1）全反射邊界。+αΦ=0（2）部分反射邊界。（5）式中：α=α1+iα2為復常數(shù)，ε為入射波和反射波之間相與復反射系數(shù)有關，復反射系數(shù)可表示為ρ=Reiε；R為振幅衰減；α與ρ的關系為α
7、1=k（6）位差；2α2=k（7）（3）入射邊界條件。對于一般波浪傳播模擬，可認為波浪是從一直線上開始向計算區(qū)域內入射傳播的，而在這條直線上有和散射勢ΦS，且不考慮入射邊界對于散射勢的反射，因此有入射勢（ΦI=a0eikx）（8）=-ikΦSikx———————————————————————————————————————————————e=（9）+=+(-ikΦS)=ikΦI-ik(Φ-ΦI)=-ikΦ+2ikΦI1.3方程求解假設將要求解的區(qū)域劃分為有限個單元，取其中一個具有代表性的單元，假定權函數(shù)為N，采用加權余量法、
8、分部積分和格林公式等，可以建立求解單元上的有限元方程，最后將單元有限元方程在全區(qū)域上進行疊加，即可得到如下的有限元方程組[K]{Φ}={Q}用雙共軛梯度法，模型的建立及具體求解過程參見LiuSX等［6］。（10）計算中采用索引存儲法，只儲存上三角的非零元素。方程的求解則采其中
