汪 鴻，李訓強*，朱首賢，張文靜

(1.解放軍理工大學 氣象海洋學院，江蘇 南京 211101；2.河海大學 海洋學院，江蘇 南京 210098)

?

波生沿岸流的數值模擬研究*

汪 鴻1，李訓強1*，朱首賢2，張文靜1

(1.解放軍理工大學 氣象海洋學院，江蘇 南京 211101；2.河海大學 海洋學院，江蘇 南京 210098)

采用完全非線性Boussinesq方程建立的FUNWAVE模型進行波生沿岸流數值模擬研究，通過對不同斜坡地形和波浪入射條件下波生流的物理模型實驗結果進行比較，驗證了該數值模型能較準確地計算沿岸流；通過改變波浪的不同入射條件，對不同入射條件的沿岸流數值模擬得出：當其他條件不變時，僅入射波高增大時，沿岸流的流幅和幅值增大，幅值位置向深水移動；僅增大入射周期時，沿岸流的流幅顯著增加，幅值的增加較?。恍逼碌匦纹露鹊母淖兡茱@著影響波生沿岸流的流幅和幅值，但對沿岸流幅值位置的水深影響不大。采用窄縫法處理動邊界時，選擇合適的窄縫起始水深對沿岸流的準確計算是十分必要的。

沿岸流；數值模擬；Boussinesq方程；窄縫法

波生沿岸流是由于近岸波浪的波高和波向分布不均勻產生的沿岸水流，它對近岸地形的塑造、物質輸運和海岸線演變有重要影響，長期以來波生沿岸流的研究受到極大的關注[1-2]。由于波生流場的易變性，其現場觀測比較困難[3-4]，因此，采用數值模型模擬研究波生流是較為有效的手段[5]。

近年來，國內外學者對波浪數值模型的研究和應用開展了很多工作，目前用于模擬近岸波生流的數值模型可分為兩類：一類為波浪時均模型，一類為波浪時域模型。第一類模型將流體運動方程在波浪周期上進行平均，求解時均方程直接得到波生流場，而波浪則通過相應的波浪理論或模型，應用輻射應力概念進行考慮，這類模型比較多見[6-8]。第二類模型則直接針對波浪質點運動進行求解，將結果在多個波浪周期內平均得到波生流場。Boussinesq類波浪模型為第二類模型的典型代表[9]。這兩種模型各有優(yōu)缺點，前者計算速度快，但不能充分考慮波流相互作用，僅通過輻射應力不能充分考慮短波(波浪)的尺度效應，如波浪非線性以及波浪-波浪之間的非線性相互作用等；后者計算速度相對較慢，但實現了波浪運動同水流運動的完全耦合求解，更能準確地反映近岸區(qū)域波流共存這一顯著特征。

基于Boussinesq方程在近岸區(qū)域實際應用的優(yōu)勢，美國特拉華大學建立了FUNWAVE數值模式[10]，Shi等[11]和Chen等[12]為其發(fā)展做了許多貢獻。李訓強等[13]使用該模式對波浪近岸傳播進行了數值模擬研究，效果比較理想。本文以沿岸流物理模型的實驗數據為檢驗資料，利用該模式計算近岸沿岸流。同時設計了沿岸流的數值模型，針對不同的波浪入射條件，對沿岸流分布特征進行分析。

1 模式介紹

1.1 控制方程

本文采用Wei等建立的Boussinesq模型[14]，在該模型的推導中，沒有做弱非線性假設，保留了O(μ2)和所有階的ε。方程形式：

(1)

(2)

式中，η是表面高度；h是靜止水深；uα是水深為zα=-0.531h處水質點的水平速度矢量；=(?/?x,?/?y)為水平梯度算子；g為重力加速度；下標t是對時間的偏導。上述方程在實際應用中應增加破碎項和摩擦項等進行改進。

1.2 波浪破碎

借鑒Kennedy[15]利用動量混合方程描述淺水區(qū)域波浪破碎的能量損耗：

(3)

(4)

式中，下標x，y表示空間導數；υ是破碎波前部的局部湍黏性應力，在非線性波計算中，參數υ非常重要。

1.3 子網紊動混合

本文采用Smagorinsky型子網模型[16]來考慮渦動黏性對波生流流態(tài)的影響，其表達式為

(3)

(4)

1.4 邊界處理

動邊界的處理方法有2種：一種是干濕網格法，一種是窄縫法[17]，本文的數值模型采用后者。窄縫法把整個計算區(qū)域作為具有窄縫或可滲透邊界、可流入流出的活動型區(qū)域，將沒有水或只有薄層水覆蓋的固體底邊界替換成可滲透的海床，或把固體海底假設具有窄縫，使得溢出的水位不會超過海灘高程。模式中設定了κ和Λ兩個參數描述窄縫法對波浪計算的影響：

(7)

(8)

式中，κ是單位寬度海灘上的窄縫寬；δ是最小縫寬；λ是窄縫的形狀參數；h0是窄縫起始的水深值；Λ是考慮窄縫影響后的等效水深。Madsen等對Carrier和Greeenspan的理論模型進行數值模擬[18]，發(fā)現上式中z*如果取為當地陸地高程，數值模擬的最大波浪爬高與理論解存在10%的誤差。Kennedy等[15]則認為這是因為波浪上爬時首先要填滿窄縫，從而引起水體質量的損失。波浪經過窄縫時，為了保持總的水體質量和沒有窄縫時的水體質量相等，模式中z*表達式為

(9)

式中zs為陸地高程。

2 物理實驗驗證

鄒志利等[19]在大連理工大學海岸和近海工程國家重點實驗室的港池內進行波生流物理模型試驗，實驗水池部分長42.6 m，寬24.0 m。物理模型簡圖見圖1，坡度均勻的平直海岸與造波板成30°角布置，海岸坡度為1∶100和1∶40，1∶100坡前水深為18.0 cm，1∶40坡前水深為45.0 cm。從造波機到模型的坡腳處布置了波導墻，并且在波導墻內側放置了消浪網以防止波浪的反射。實驗測量了坡度分別為1∶100和1∶40兩種平直海岸上的波生沿岸流。波面升高由垂直岸線方向排列的三列電容式波高儀測量，沿岸流由應變式流速儀測量。選取如圖1中所示沿岸線方向離岸距離13.0 m觀測截面上的波高和流速實驗值對模式進行驗證。

圖1 實驗模型Fig.1 Diagram of the experiment setup

本文使用FUNWAVE模式對上述物理模型進行數值模擬，模式采用0.2 m×0.4 m的數值網格，造波源設置為數值模式的波浪入射邊界，由于造波為雙向造波，因此在造波源位置后側設置海綿邊界進行消波。斜坡上的移動邊界采用上述窄縫法進行處理，為了減小FUNWAVE模式側邊界的影響，本文設置了長45.0 m，寬34.0 m的數值水槽模型，分別去除左右側邊界附近5 m的計算結果，選取數值模型中部區(qū)域的24.0 m的計算結果進行分析。由靜止狀態(tài)開始計算波浪運動，模式運行70 s后，波面呈現很有規(guī)律的變化，波高和流場基本穩(wěn)定。選取與物理實驗條件對應的一組數值模型計算結果，圖2a,2b分別為波高H=0.09 m，周期T=1.0 s的規(guī)則波在1∶40平直海岸上波高和沿岸流模擬結果。圖2a中，波浪波高在斜坡前變化不大，在斜坡起始線處波高分布出現變化。在岸線附近，波浪由于淺化過程出現破碎并在破碎后波高迅速衰減，輻射應力變化導致波生沿岸流的產生。圖2b為數值模型計算的沿岸流分布，圖中沿岸流在沿岸方向上的分布均勻合理。此外，源函數造波(計算區(qū)域左端)和海綿吸收層(計算區(qū)域左端)的處理結果比較有效。

圖2 FUNWAVE模擬波高和沿岸流分布Fig.2 Wave height and longshore current simulated by the FUNWAVE model

本文選取與實驗觀測截面一致的斷面數值計算值與上述物理實驗結果比較。圖3為1∶100和1∶40地形下，兩組入射波的波高和波生流在垂直于岸線斷面上的實驗觀測值和數值模型計算值的比較結果。其中，圖3a，圖3b為1∶40坡度下，入射波高Hλ=0.09 m，周期T=1.0 s規(guī)則波的波高和沿岸流分布，圖3c,3d為1∶100地形坡度下，波高為H=0.07 m，周期T=1.0 s規(guī)則波的波高和沿岸流分布。圖3中，x軸為垂直岸線的水平距離，由圖3可以發(fā)現，在兩種坡度下，不同入射條件的波高和沿岸流數值模擬結果與實驗模型的數據基本一致。在近岸區(qū)域，波浪由于淺化作用發(fā)生堆積破碎，破碎后波浪波高迅速減小。波生流的最大流速位置與波浪破碎線并不重合，而是落后于波浪破碎位置。通過圖3對數值模式的驗證，證明該模式可以較好地計算近岸波生流。

圖3 數值模型和物理模型結果Fig.3 Simulated result by the FUNWAVE model versus that from experimental data

3 不同條件對波生沿岸流影響的數值分析

基于上述研究表明，FUNWAVE模式模擬近岸波生流和波高分布的效果較好。由于波生流物理模型實驗的易變形性和復雜性，現采用FUNWAVE模式對不同波高、周期、坡度條件下的沿岸波生流進行數值模擬研究，從而對不同入射條件下近岸波生流的特征進行分析。

選擇地形坡度1∶40進行數值實驗，保持坡前水深0.45 m不變。圖4a為1∶40坡度模型下，波浪入射周期T=1.0 s不變、不同入射波高條件下波生沿岸流在垂直岸線斷面上計算結果。圖4a中，從上往下依次為入射波高為0.13，0.11，0.09，0.07和0.05 m條件下沿岸流隨離岸距離的分布。由圖4a可知，當波浪的入射周期為1 s不變時，隨著波浪入射波高的增大，波浪破碎線向深水處移動，波生沿岸流的幅值明顯增大，沿岸流的最大值位置向深水處移動，且沿岸流的影響區(qū)域也相對增大。圖4a中，波浪波高的增幅為0.02 m不變，可以發(fā)現，保持波浪波高的增幅一定，隨著波高的增大，波生沿岸流的幅值及幅值的位置變化范圍減小。參照上述設置，圖4b為保持波浪入射波高H=0.09 m不變、不同周期下的波生沿岸流在垂直岸線斷面上的數值計算結果。圖4b中，3種周期條件計算的沿岸流幅值位置變化較小，隨著波浪入射周期的增大，沿岸流的流幅顯著變寬，幅值略微增大。周期的增加能顯著改變沿岸流寬度，但對沿岸流幅值大小以及位置的影響較小。

圖4 不同入射波高、周期條件下波生沿岸流計算結果Fig.4 Simulated longshore current with different incident height and period

本文通過改變1∶40坡度地形下的坡前水深至0.18 m，保持其余參數不變，研究不同坡度下波生沿岸流的性質(圖5)。圖5為波浪入射波高Hλ=0.07 m，周期T=1.0 s不變，1∶40和1∶100兩種地形坡度下，波生沿岸流在垂直岸線斷面上的計算值。圖5中1∶40較陡坡度下，波生沿岸流的流幅相對較小，且幅值較小。1∶40坡度下，沿岸流的幅值位置相對1∶100緩坡幅值位置向岸線偏移。通過兩種坡度下幅值位置的離岸距離計算兩種坡度下幅值位置的水深可以發(fā)現，1∶40和1∶100坡度下的沿岸流幅值位置附近的水深幾乎相同，為0.073 m。地形坡度的變化能顯著改變沿岸流的流幅和幅值，但對沿岸流幅值位置水深的影響不大。

圖5 不同坡度下波生沿岸流計算結果Fig.5 Simulated longshore current with different slopes

選擇1∶40坡度下，入射波高Hλ=0.09 m、周期T=1.0 s的規(guī)則波數值模型，保持其余參數不變，通過改變平直海岸與造波源的夾角，研究不同入射角條件的波生流性質(圖6)。圖6中，由上到下分別為夾角為30°，25°和20°條件下，波生沿岸流在垂直于岸線斷面上的計算結果。圖6所示，隨著入射角的增大，波生流的幅值相對增大，幅值位置和流幅區(qū)向深水處偏移。

圖6 不同入射角下沿岸流計算結果Fig.6 Simulated longshore current with different incident angles

4 窄縫法對沿岸流數值計算的影響

FUNWAVE模式采用窄縫法處理動邊界，為了保證波浪的爬高時的水體質量守恒，在使用窄縫法處理動邊界時，z*不再取值陸地高程，其取值如式(9)，從而使FUNWAVE模型計算的波浪爬高較為準確。本文在使用FUNWAVE模式計算沿岸流過程中，發(fā)現窄縫的最小縫寬(δ)和形狀參數(λ)對波高和沿岸流計算的影響不大，而窄縫的起始水深值(h0)對沿岸流數值計算的影響較大，現對其進行初步討論。

圖7 不同窄縫起始水深對數值計算的影響Fig.7 Simulated wave height and longshore current with different slot depths

本文選取1∶40坡度下，波高0.09 m的規(guī)則入射波進行分析，其余設置參照上述不變。通過設置不同的窄縫起始水深，討論不同起始水深對垂直岸線斷面上的波高和沿岸流計算的影響。圖7為窄縫起始水深0.16,0.20,0.24,0.28和0.32 m時，垂直岸線斷面上的波高和沿岸流計算結果。圖7a中，起始水深對波浪波高的計算影響較小，不同窄縫起始水深計算的波高在很小的范圍內變化。圖7b中，不同窄縫起始水深設置計算的沿岸流分布有較大的差異：不同窄縫起始水深下，沿岸流計算的幅值位置相對不變，幅值大小變化較大。保持窄縫起始水深變化幅度0.02 m時，沿岸流的幅值差異最大可達0.10 cm/s。通過上述討論，采用窄縫法計算波浪波高較為精確，但在計算沿岸流時，選擇合適的窄縫起始水深是十分重要的。

5 結 論

本文通過采用采用FUNWAVE模式對物理模型進行數值實驗，證明FUNWAVE模式能較好的模擬近岸波生沿岸流。本文利用數值模擬耗費較少，易于調整，以及可以進行重復性實驗的優(yōu)勢，通過斜坡地形下沿岸流的數值模擬研究，對改變單一入射條件情況的波生沿岸流進行數值模擬，得出：1)僅增大波浪的入射波高，其余入射條件不變時，波生沿岸流的幅值和流幅范圍增大，幅值位置向深水移動。保持波浪波高增幅一定時，隨著入射波高的增大，波生沿岸流幅值及幅值位置的變化幅度減小。2)僅增大波浪的入射周期，其余入射條件不變時，波生沿岸流的流幅顯著增大。3)僅增加斜坡地形坡度，其余入射條件不變時，波生沿岸流的流幅和幅值減小，地形坡度的改變對幅值水深影響非常小。4)僅增大波浪入射角，其余入射條件不變時，沿岸流的幅值增大，且幅值位置和流幅向深水移動。

在使用FUNWAVE模式進行沿岸流計算時，窄縫法設置的起始水深對沿岸流計算有較大影響，合適的窄縫起始水深對保證流場計算的準確度時十分必要的。FUNWAVE模型的計算網格精度要求較高，對于小區(qū)域的波浪要素計算效果較好，但在計算大區(qū)域波浪傳播過程，其計算量較大，耗時較長，在實際海區(qū)的流場計算時有一定的困難。

[1] SVENDSEN I A. Introduction to Nearshore Hydrodynamics[J].World Scientific Publishing Co Pte Ltd, 2006:744.

[2] SUN T, TAO J H. Numerical modelling and experimental verification of pollutant transport under waves in the nearshore zone[J].Haiyang Xuebao, 2003,25(3):104-112.孫濤,陶建華. 波浪作用下近岸區(qū)污染物輸移擴散的數學模型及其實驗驗證[J].海洋學報,2003,25(3):104-112.

[3] REN C P, ZOU Z L. Research development of the shear instability of longshore current[J].Marine Science Bulletin, 2010,29(5):588-593.任春平,鄒志利.沿岸流不穩(wěn)定運動研究進展[J].海洋通報,2010,29(5):588-593.

[4] BIRKEMEIER W A, DONOGHUE C, LONG C E, et al. 1990 DELILAH nearshore experiment: Summary Report[R].Delilah Nearshore Experiment Summary Report, 1997.

[5] CHOI J, KIRBY J T, YOON S B. Boussinesq modeling of longshore currents in the SandyDuck experiment under directional random wave conditions[J].Coastal Engineering, 2015,101:17-34.

[6] HAAS K A, WARNER J C. Comparing a quasi-3D to a full 3D nearshore circulation model: SHORECIRC and ROMS[J].Ocean Modelling, 2009,26(1-2):91-103.

[7] TANG J, SHEN Y, QIU D. Numerical study of pollutant movement in waves and wave-induced long-shore currents in surf zone[J].Acta Oceanologica Sinica, 2008,27(1):122-131.

[8] FANG K Z, ZOU Z L, LIU Z B. Numerical simulation of rip current generated on a barred beach[J].Chinese Journal of Hydrodynamics, 2011,26(4):479-486.房克照,鄒志利,劉忠波.沙壩海岸上裂流的數值模擬[J].水動力學研究與進展,2011,26(4):479-486.

[9] LI S W, HUANG X Y. Numerical method for calculation of longshore current by using Boussinesq equations[J].Journal of Tianjin University.2004,37(12):1059-1062.李紹武,黃筱云.用Boussinesq方程計算沿岸流的數值方法[J].天津大學學報:自然科學與工程技術版,2004,37(12):1059-1062.

[10] KIRBY J T, WEI G, CHEN Q, et al. FUNWAVE 1.0, fully nonlinear Boussinesq wave model documentation and user's manual[J].University of Delaware, 1998.

[11] SHI F, DALRYMPLE R A, KIRBY J T, et al. A fully nonlinear Boussinesq model in generalized curvilinear coordinates[J].Coastal Engineering, 1999, 42(4):337-358.

[12] QIN C, JAMES T K, DALRYMPLE R A, et al. Boussinesq modeling of wave transformation, breaking, and runup. II: 2D[J].Journal of Waterway Port Coastal & ocean Engineering, 2000,126(1):39-47.

[13] LI X Q, YANG J X, ZHU S X, et al. Numerical simulations of surfing waves on mild slopes[J].Advance in Marine Science,2015,33(2):155-163.李訓強,楊家軒,朱首賢,等.緩坡地形拍岸浪數值模擬研究[J].海洋科學進展,2015,33(2):155-163.

[14] WEI G, KIRBY J T, GRILLI S T, et al. A fully nonlinear Boussinesq model for surface waves. I. Highly nonlinear, unsteady waves.[J].Journal of Fluid Mechanics, 1995,294:71-92.

[15] KENNEDY B, CHEN Q, KIRBY J T, et al. Boussinesq modeling of wave transformation, breaking, and runup. I: 1D[J].Journal of Waterway Port Coastal & Ocean Engineering, 2002,126(1):39-47.

[16] SMAGORINSKY J. General circulation experiments with the primitive equations[J].Monthly Weather Review, 1963,91(3):99-164.

[17] MADSEN P A, MURRAY R. A new form of the Boussinesq equations with improved linear dispersion characteristics[J].Coastal Engineering, 1991,15(4):371-388.

[18] CARRIER G F, GREENSPAN H P. Water waves of finite amplitude on a sloping beach[J].Journal of Fluid Mechanics, 1958,4(1):97-109.

[19] ZOU Z L, CHANG M, QIU D H, et al. Experimental study of longshore currents[J].Journal of Hydrodynamics, 2002,17(2):174-180.鄒志利,常梅,邱大洪,等.沿岸流的實驗研究[J].水動力學研究與進展, 2002,17(2):174-180.

Received: December 2, 2015

Numerical Simulation of Longshore Currents Induced by Waves

WANG Hong1,LI Xun-qiang1,ZHU Shou-xian2,ZHANG Wen-jing1

(1.InstituteofMeteorologyandOceanography,PLAUniversityofScienceandTechnology, Nanjing, 211101, China；2.InstituteofOceanography,HohaiUniversity, Nanjing 210098, China)

The FUNWAVE model based on nonlinear Bosssinesq equations is employed to simulate the longshore current induced by waves. Comparison with the experimental results of physical model show that the FUNWAVE model can well simulate the longshore currents. The model is further employed to the simulate longshore current on different slopes, with varying heights and periods, and results show that when incident height or period increases, the magnitude and width of the longshore current will also increase, and the location of maximum longshore current moves towards deeper water as incident height increases. When terrain slope changes, the magnitude and width change significantly but the location of maximum longshore current does not change. FUNWAVE model adopts slot technique to deal with transient boundary, the initial depth of slot technique has an important impact on the longshore current simulation.

longshore current; numerical simulation; Boussinesq equation; slot technique

2015-12-02

國家自然科學基金項目——引入波浪漫灘邊界的近岸海域懸沙輸運數值模式及其應用(41076048)；中央高?；究蒲袠I(yè)務費專項項目——流和溫鹽垂向異結構的新型海洋模式研制及其應用(2011B05714)

汪 鴻(1992-)，男，碩士研究生，主要從事海洋數值模擬方面研究.E-mail：wanghong_92@126.com

*通訊作者：李訓強(1962-)，男，江蘇睢寧人，副教授，碩士研究生導師，主要從事海洋動力學與數值模擬方面研究.E-mail：lixunqiang@sina.com

P731

A

1671-6647(2016)04-0507-09

10.3969/j.issn.1671-6647.2016.04.006