不同尺度人工岛波浪绕射系数的研究❋

于定勇, 李 龙， 张广成

(中国海洋大学工程学院，山东 青岛 266100 )

不同尺度人工岛波浪绕射系数的研究❋

于定勇, 李 龙， 张广成

(中国海洋大学工程学院，山东 青岛 266100 )

为研究人工岛沿波浪传播方向尺度变化对绕射系数的影响，本文基于MIKE21-BW模型，数值模拟研究了入射波为不规则波时人工岛尺度不同对其绕射系数的影响。通过数值结果与Briggs物理模型试验结果进行了对比，结果吻合较好，验证了模型在研究波浪绕射问题上的适用性。在未考虑结构物沿波浪传播方向尺度对绕射系数的影响时，将数值结果与《海港水文规范》(JTS145-2-2013)规范结果比对，结果差异较小；本文采用上述模型模拟了结构物沿波浪传播方向7种不同尺度时的波浪绕射情况，分别是0.5Ls,Ls, 2Ls, 3Ls, 4Ls, 5Ls及6Ls，通过对不同尺度时的波浪绕射系数变化的分析，得到结果：结构物尺度在0.5Ls～2Ls范围内时，绕射系数随着结构物尺度的增大而变大；结构物尺度在3Ls～6Ls范围内时，绕射系数随着结构物尺度的增大而减小；在人工岛岛壁附近和其掩护区内绕射系数变化较明显；人工岛结构物尺度变化对绕射系数值有影响。文中给出了不同尺度人工岛绕射系数等值线图，可为工程设计中绕射波高确定提供参考。

BW模型；人工岛；绕射系数；尺度；数值模拟

大型填海工程是我国海洋开发活动的重大举措，其中离岸式人工岛的开发及建设成为利用海洋资源的一种有效方式。人工岛从使用功能来说，可以作为修建深水泊位的海港作业区、海上机场[1-2]、海上旅游休闲区[3]等。波浪是近岸建筑物遭受的主要动力之一，人工岛在所处海域受波浪的折射、反射和绕射等作用，因此确定人工岛设计波高显得尤为重要。波浪绕射问题的研究对确定岛上装卸码头的泊稳条件、岛周围的护岸设计、挡浪墙高度及岛周围海底泥沙的局部冲淤都有着重要的现实意义。

外海深水区海浪传播到浅水区的过程，受地形、海底摩阻、近岸结构物等影响，波浪传播过程中波浪要素随之变化，研究波浪要素的变化对于海岸工程设计有着重要作用。如今对于建立波浪近岸传播过程已有多种成熟数学模型，如基于缓坡方程的数学模型[4]、基于能量平衡方程的模型[5]和基于Boussinesq型方程的模型[6]等。近年来，由于对Boussinesq方程的不断完善，不论在色散性、非线性方面都得到了很大的提高，使得它成为描述波浪近海传播变形的理想模型。

一些学者对波浪作用于人工岛的问题进行了相关研究。蔡艳君等[7]在规则波绕射数学模型的基础上，运用线性迭加不同振幅、频率和方向组成波的方法，采用文氏谱和COS2nθ型的方向分布函数建立了人工岛周围不规则波绕射的数学模型，总结了人工岛周围不规则波绕射的分布特征。陈新等[8]应用BW模型对直立式圆柱形人工岛周围的波高分布进行了数值研究，得到了圆柱形人工岛周围波高分布规律及人工岛直径变化对反射、绕射波高比的影响。竺艳蓉[9]基于孤立波理论研究了浅海波浪对人工岛的非线性绕射，提出了采用有限差分法和有限远开边界的幅射条件评算人工岛上波浪爬高的数值模型。得到的数值计算结果与整体模型试验的实测结果进行了比较，二者吻合较好。李德筠等[10]开发了基于射线理论的浅水波计算程序，提出了一种解决射线相交的简便方法。应用非线性长波方程，以波浪对圆柱形人工岛的绕射、折射及对岸坡的斜向入射为例，进行了数值计算，获得了数值结果与实际观测相吻合的结果。何国华等[11]针对计算域中存在直立岛式结构物的复连通区域，基于时间关联型缓坡方程和相应的边界条件建立了波浪传播的数值模拟模型。数值结果表明模型可以有效模拟计算域内存在直立岛式结构物的绕射和反射问题。

我国现行的《海港水文规范》(JTS145-2-2013)[12]给出了推算岛式防波堤堤后的绕射波高的规范方法，但规范方法未反映人工岛的尺度及尺度变化对其绕射系数的影响。故本文应用MIKE21-BW模型，在验证模型研究波浪绕射问题适用性的基础上，给出了不同尺度人工岛绕射系数等值线的分布图，并初步探讨了人工岛不同尺度对其掩护区绕射系数的影响。

1 数值模型介绍

MIKE21-BW模型是基于Boussinesq型方程所建立的波浪数学模型，其中所解的Boussinesq方程包含了非线性项和频率耗散项。它能够较好的模拟不同地形及平面布置条件下近岸波浪的传播与变形，包含的物理现象：折射及浅水变形、绕射、底摩阻损耗、部分反射或透射、波浪破碎等。模块中的计算方程有两种类型：经典的Boussinesq方程和改进的Boussinesq方程。计算方程选择的依据是最大水深与深水波长的比值，经典的Boussinesq方程适用于最大水深与深水波长的比值(hmax/L0)小于0.22；改进的Boussinesq因包含了深水项，并结合了改进的色散关系，从而可模拟不规则波在较深水域中或较小周期波的传播，故方程适用于最大水深与深水波长的比值(hmax/L0)小于0.5的情况。

1.1 控制方程：

模型控制方程采用Beji和Nadaoka[13]改进后的方程，连续性方程：

(1)

x方向动量方程：

(2)

y方向动量方程：

(3)

这里Ψ1、Ψ2为Boussinesq项，由以下式子定义：

(4)

式中：下标x、y和t分别表示对空间和时间的偏微分；P为x方向的流密度；Q为y方向的流密度；B为深水修正系数；Fx为x方向的水平力；Fy为y方向的水平力；d为静水水深；ξ为波面相对于静水面的高度；h为总水深(=d+ξ)；n为空隙率；C为谢才系数；α为层流阻尼系数；β为紊流阻尼系数。

(5)

式中：vt为水平方向的涡流速度；Rxx、Rxy和Ryy表示由非均匀速度引起的剩余动量：

(6)

Sorensen等[14]给出了完善的计算方法求解控制方程的数值方法，空间离散格式采用图1所示矩形网格，求解微分方程时采用交替方向隐式(Alternating Direction Implicit, ADI)算法进行求解。水面高程等标量定义在网格节点上，而流量分量等矢量定义于相应方向的网格线的中点上。

图1 BW模型离散格式Fig.1 Discretization scheme of BW model

1.2 波浪入射边界

在模型中采用沿某一直线进行流量变动的造波方法，即内波生成线法。内生波可以在生成线后设置消波海绵层，其作用是超出模型区域的波浪可以被吸收，且内波生成线与海绵层相互独立。本文模型中不规则波频谱采用JONSWAP谱[15]及方向分布函数采用与频率相关的表达式[16]，频谱表达式为：

(7)

其中：

1.094-0.019 15lnγ

(8)

式中：H1/3为有效波高；TP为谱峰周期；γ为峰高因子，取3.3；峰形参数σ=σa(当ω≤ωm时)，σ=σb(当ω>ωm时)，σa，σb分别取0.07，0.09。

方向分布函数为与频率相关的表达式：

(9)

2 数值模型的验证

2.1 试验模型简介及数值模型建立

在Briggs等[17]的试验中设置的试验区域如图2。单突堤的厚度取0.15 m，水深为0.4 m，波浪的入射边界与单突堤平行，且二者之间的距离为5 m，入射边界相邻的两侧边界取为全反射[18]，部分边界设置为吸收边界。在试验过程中选取堤后与堤轴线夹角为30°、60°和90°的三个断面为研究对象，在断面布置波高仪测量波高数据，入射波波型分别选取规则波及多向不规则波。

图2 单突堤的波浪绕射计算区域Fig.2 Calculation domain for wave diffraction around a semi-infinite breakwater

为了与Briggs等的试验结果比较，基于BW模型建立了与物理试验相同的数值计算模型。计算区域、计算步长及波浪要素采用表1数据，其中不规则波频谱采用JONSWAP谱，多向不规则波方向分布的最大偏差角度为30°，规则波的波高H=0.077 5 m，周期T=1.30 s，波长L0为2.25 m，水深波高比h/H=5.16。在数值模拟中，数值模型网格步长取0.2 m，时间步长取0.1 s。因模型在规则波及不规则波作用时最大水深与波长之比的值均为hmax/L=0.18<0.5，故在BW模型的适用范围内。波浪绕射系数采用下式计算：

(10)

式中Hi和HI分别为绕射波高和入射波高。

表1 数值模拟计算组次Table1 Calculated groups of numerical simulation

2.2 对比结果分析

对数值结果提取了与堤轴线夹角为30°，60°和90°三个断面的数据，与Briggs等[9]所建立的物理试验中对应断面测点的绕射系数值Kd进行了比较，对比结果见图3。r′代表某点绕射系数值在距堤头o点的无因次距离(r′=r/L，r代表某点距堤头o点的距离，L代表波长)。由数值结果和试验值对比可以发现数值结果与试验值大致相吻合。在规则波和多向不规则波条件下，30°和60°断面上绕射系数值基本一致。90°断面上数值结果和实验数据相差相对较大，M1组次的规则波模拟结果比试验数据值较大，数值结果与实测数据最大相差约19.6%，平均差值约11.0%，因物理试验中规则波经过单突堤，波浪在堤后90°断面仅有绕射波高，而数值模拟结果在此断面有绕射波高和反射波高的叠加，数值模拟给出的波高值大于试验值，所以根据试验数据给出的绕射系数值小于本文根据数值模拟结果给出的绕射系数值；M2组次的试验数据大于模拟值，最大相差约10.2%，平均差值为7.2%，可能原因是物理试验造波技术的影响，且试验和数值模拟过程中多向不规则波的方向分布有一定差别。总的来说，试验数据与数值结果对比结果吻合较好，得知模型在研究波浪绕射问题是可行的。

图3 数值结果与Briggs试验值对比Fig.3 Compared numerical results with Briggs’ experimental data

3 人工岛尺度的变化对绕射系数的影响

3.1 数值模型的建立

本文建立了封闭的矩形计算区域，区域范围为800 m×1 000 m。上下边界设置了海绵层吸收边界，在上边界海绵层的结束处设置了内波生成线，波浪自上而下传播，入射波为不规则波，波浪频谱采用JONSWAP谱，方向函数选用与频率相关的表达式，有效波高Hs和谱峰周期TP分别为1.5 m和6.0 s，有效波长Ls约48 m，水深取10 m，岛堤高程+10 m，水深波高比h/H=6.67，所建模型最大水深与波长之比hmax/Ls=0.21<0.5，适用于BW模型。海绵层厚度取100 m，约为有效波长的两倍，能够较好的吸收波能，减小反射波对研究区的影响。数值模型空间网格取2 m，时间步长取0.1 s。在计算区域内部布置人工岛，岛壁设为全反射边界，反射系数为1.0。入射波边界相邻两侧边界为全反射，其目的是扩大有效的计算区域。人工岛的波浪入射向的尺度定义为宽度尺度，文中所述尺度表示宽度尺度。图4给出了数值研究区域及12个特征点的位置，在岛堤的掩护区建立如图所示的直角坐标系，1#～12#特征点的空间坐标分别为(0,0)、(l/6,0)、(l/3,0)、(l/2,0)、(0,2Ls)、(l/6,2Ls)、(l/3,2Ls)、(l/2,2Ls)、(0,4Ls)、(l/6,4Ls)、(l/3,4Ls)、(l/2,4Ls)，Ls表示入射波有效波长，l表示岛堤长。

图4 研究区域及选取的特征点示意图Fig.4 Research area and the locations of the calculating points

3.2 海港水文规范相关条款

《海港水文规范》(JTS145-2-2013)中关于绕射波高求解及岛式防波堤后的不规则波绕射系数的确定方法如下所述。岛堤后某点的绕射波高可按下式计算：

Hd=KdHi，

(11)

式中：Hd为防波堤后某点的绕射波高(m)；Kd为防波堤后某点的绕射系数；Hi为防波堤入射波波高。不规则波绕射系数的确定是根据规范中9种不同工况绕射系数等值线图7.2.3确定，各图的绕射系数可通过内插取得。规范中9种绕射系数图分别为岛堤长度l/Ls=4、6、8与波浪入射方向θ0=90°、60°、30°之间的组合，θ0表示主波入射角。规范中不规则波绕射系数的确定方法是根据岛堤后规则波绕射系数数值结果，再用能量线性叠加原理得出了岛堤后的不规则波绕射系数。

3.3 绕射系数数值结果与规范值对比

对不规则波经过等水深且未考虑宽度尺度的岛式防波堤按照《海港水文规范》(JTS145-2-2013)工况进行了数值模拟。将数值计算结果分别按照主波入射方向θ0=90°、60°、30°，岛堤长度l/Ls=4、6、8绘制掩护区的绕射系数等值线图，并与规范的系数等值线进行了对比。同时，选取了岛堤掩护区12个特征点，分别给出了特征点的数值结果、规范值及差异。

图5给出了在不规则波入射条件时不同岛堤长度、入射角度岛式防波堤掩护区数值结果与规范的绕射系数等值线的对比，(a)、(b)、(c)图为岛堤长度为l/Ls=4、6、8，入射波主波向为90°时等值线对比，(d)、(e)图为岛堤长度为l/Ls=6，入射波方向分别为60°、30°时绕射系数等值线对比。图中实线为规范值，虚线为数值结果。表2给出了不规则波且主波向θ0=90°，堤长l/Ls=4、6时12个特征点的绕射系数数值结果、规范值及差异。由图4绕射系数对比得到：波浪正向入射时，开敞区1.0、0.9等直线的数值结果略大于规范值，且0.9绕射系数等值线吻合较好；在掩护区0.2～0.8等值线数值结果小于规范值；在l/Ls=4、6、8时，y/Ls分别在-4.0～0、-6.0～0、-8.0～0范围内数值结果与规范值基本吻合。在波浪入射角为60°时，0.8、0.9绕射系数数值结果略大于规范值，而0.2～0.7绕射系数数值等值线小于规范值；入射角为30°时绕射系数等值线大于规范值。由表2得：在l/Ls=4、θ0=90°时，5#、7#、8#特征点的差异小于5%；而3#、4#特征点大于20%，分别为23.8%、20.0%。在l/Ls=6、θ0=90°时，1#、5#、10#特征点差异分别为0.0%、1.5%、3.9%，差异较小；4#、7#、12#特征点规范值和数值结果的差异分别为20.0%、22.2%、27.0%，差异相对较大；总的来说，12个特征点的数值结果与规范值相差较小。由图5、表2知模型能够较好的模拟岛堤后的绕射系数值，在与《海港水文规范》相同工况条件下，规范值和数值结果的绕射系数对比结果差异较小。

3.4 人工岛不同尺度对绕射系数的影响

为研究人工岛入射波方向尺度的不同对其绕射系数Kd的影响，人工岛长度取固定值l=6Ls，宽度尺度分别取0.5Ls、Ls、2Ls、3Ls、4Ls、5Ls、6Ls，入射波取不规则波，入射主波向θ0=90°。表3给出了各特征点的规范值及0.5Ls～6Ls宽度尺度时1#～12#特征点绕射系数数值结果。图6给出了宽度尺度为Ls、2Ls、3Ls及4Ls绕射系数等值线分布。

图5 岛式防波堤后绕射系数数值结果与规范值等值线对比Fig.5 Comparison of numerical results and Code diffraction coefficients behind offshore breakwater

工况Item特征点Characteristicpoints1#2#3#4#5#6#7#8#9#10#11#12#l/Ls=4θ0=90°规范值Codevalues0．600．290．210．200．640．400．280．250．680．530．430．39数值结果Numericalresults0．640．310．260．240．630．370．280．260．640．500．380．34差异Difference/%6．76．923．820．01．67．50．04．05．95．711．612．8l/Ls=6θ0=90°规范值Codevalues0．600．260．160．150．650．390．270．230．690．510．410．37数值结果Numericalresults0．600．230．190．180．660．320．210．200．650．490．340．27差异Difference/%0．011．518．820．01．517．922．213．05．83．917．127．0

表3 不同尺度特征点绕射系数值Table 3 Diffraction coefficients of character points at different scales

图6 人工岛4种宽度尺度绕射系数等值线分布Fig.6 Distributions of diffraction coefficients of four scales of artificial islands

由表3得，人工岛宽度尺度在0.5Ls～6Ls范围，随着岛体宽度尺度的增大，选取的12个特征点绕射系数变化趋势较平缓；人工岛宽度尺度在0.5Ls～2Ls范围内，多数特征点的绕射系数值随着尺度的增大而增大，除个别点外，如10#特征点；在3Ls～6Ls范围内，特征点的绕射系数值大都随着宽度的增大而减小，但减小的幅度较小。

在0.5Ls～6Ls尺度范围内特征点1#、2#、3#、4#的绕射系数变化幅度约23.33%、16.67%、20.0%、26.09%，变化较大的原因是特征点位于岛壁边缘，受到直立岛壁宽度变化的波浪反射的影响，特征点绕射系数值有较大的变化。绕射系数值变化幅度较小的是5#、9#，变化幅度约6.35%、3.03%，原因是此特征点处于开敞区附近波浪绕射作用影响较小；在4#、7#、8#、12#特征点绕射系数的变化幅度约26.09%，31.82%，19.05%及17.85%，说明在掩护区内部的特征点随着尺度的变化，绕射系数的变化幅度较大。9#～12#特征点的绕射系数变化分别为3.03%，10.87%，11.76%，17.85%，说明随着距开敞区距离的增大，绕射系数变化幅度增大。

图6绘制了在不规则波正向入射时人工岛4种宽度尺度岛后绕射系数等值线的分布，岛后掩护区绕射系数分布在0.2～0.7。在岛堤堤头处绕射系数等值线分布较密集，说明在此绕射波高的变化较明显。各图之间绕射系数可内插，岛后的绕射系数分布为推算绕射波高提供参考。

4 结论

本文基于BW模型研究了人工岛尺度不同对其掩护区波浪绕射系数的影响，将模型的数值结果与已有物理试验数据、《海港水文规范》(JTS145-2-2013)规范值进行了对比，给出了在不规则波条件时不同尺度人工岛掩护区的绕射系数等值线分布及不同尺度对其绕射系数的影响。

(1)BW模型在研究波浪绕射问题上具有适用性。

(2)未考虑入射波方向人工岛尺度对绕射系数的影响，采用与规范中岛式防波堤相同的工况，绕射系数的数值结果与规范值对比结果基本一致。

(3)人工岛入射波方向尺度的变化对其绕射系数有一定影响。人工岛宽度尺度在0.5Ls～2Ls范围，绕射系数值随着尺度的增大而增大；在3Ls～6Ls范围，绕射系数值随着宽度的增大而减小，但减小幅度较小；在岛壁附近及掩护区域内部绕射系数变化较大，且随距开敞区距离的增大绕射系数变化幅度增大。

(4)本文给出了在不规则波正向作用时人工岛4种宽度尺度的绕射系数等值线分布图，各图之间绕射系数可内插，为推算人工岛掩护区绕射波高提供了参考。

[1] 中濑, 明男. 关西国际机场人工岛建设[J]. 海岸工程, 1994, 13(2): 197-201.

Zhong L, Ming N. Research on construction of artificial island of Kansai international airport [J]. Coastal Engineering, 1994, 13(2): 197-201.

[2] 邱毅平. 澳门国际机场海上人工岛工程设计简介[J]. 上海水利, 1993, 1(3): 28-33.

QIU Y P. Brief introduction on design of artificial island engineering of Macao international airport [J]. Shanghai Water, 1993, 1(3): 28-33.

[3] 谢浩. 迪拜奇迹之帆船酒店[J]. 绿色建筑, 2010, 2(1): 3-4.

Xie H. Research on the Burj Al Arab hotel of the Dubai miracle [J]. Green Building, 2010, 2(1): 3-4.

[4] Berkhoff J C W. Computation of combined refraction diffraction [C]. New York: Proceedings of the 13th Conference Coastal Engineering, ASCE, 1972: 471-490.

[5] Orszaghova J, Borthwick A G L, Taylor P H. From the paddle to the beach-A Boussinesq shallow water numerical wave tank based on Madsen and Sorensen’s equations [J]. Journal of Computational Physics, 2012, 231: 328-344.

[6] Rusu L, Soares C G. Evaluation of a high-resolution wave forecasting system for the approaches to ports [J]. Ocean Engineering, 2013, 58: 221-238.

[7] 蔡艳君, 谢世楞. 人工岛周围的不规则波绕射计算[J]. 港工技术, 1989, 4: 1-9.

CAI Y J, XIE S L. Diffraction calculation of random wave around artificial island [J]. Port Engineering Technology, 1989, 4: 1-9.

[8] 陈新, 刘明. 应用MIKE21对人工岛周围波高分布的数值模拟[J]. 中国水运, 2012, 12(3): 60-63.

CHEN X, LIU M. Numerical simulation on distribution of wave height around artificial island [J]. China Water Transport, 2012, 12(3): 60-63.

[9] 竺艳蓉. 浅海人工岛波浪爬高的数值计算和实验验证[J]. 水利学报, 1994, 9: 10-20.

ZHU Y R. The numerical calculation and experimental verification of wave runup around the artificial island in the Shallow Sea [J]. Journal of Hydraulic Engineering, 1994, 9: 10-20.

[10] 李德筠, 沈国光. 变深区浅水波的绕射和折射计算[J]. 海洋学报, 1995, 17(2): 106-116.

LI D J, SHEN G G. Diffraction and refraction calculation on shallow water wave at variable deep area [J]. Acta Oceanologica Sinica, 1995, 17(2): 106-116.

[11] 何国华, 陈婕. 直立岛式结构物周围波浪传播的数值模拟[J]. 浙江水利科技, 2013, 1(2): 11-14.

HE G H, CHEN J. Numerical simulation of the wave propagation around vertical island structures [J]. Zhejiang Hydrotechnics, 2013, 1(2): 11-14.

[12] JTS 145-2-2013. 海港水文规范[S]. 北京: 人民交通出版社, 2013: 30-42.

JTS 145-2-2013. Code of Hydrology for Sea Harbor [S]. Beijing: China Communications Press, 2013: 30-42.

[13] Beji S, Nadaoka K. A formal derivation and numerical modeling of the improved Boussinesq equations for varying depth [J]. Ocean Engineering, 1996, 23(8): 591-704.

[14] Sorensen O R, Hemming A, Sorensen L S. Boussinesq-type modeling using an unstructured finite element technique [J]. Coastal Engineering, 2004(50): 181-198.

[15] 俞聿修. 随机波浪及其构成应用[M]. 大连: 大连理工大学出版社, 2000: 149-154.

YU Y X. Random Wave and Its Applications to Engineering [M]. Dalian: Dalian University of Technology Press, 2000: 149-154.

[16] DHI. MIKE21 Toolbox User Guide [M]. Denmark: DHI, 2009: 148-159.

[17] Briggs M J, Thompson E F, Vincent C L. Wave diffraction around breakwater [J]. Journal of Waterway Port Coastal and Ocean Engineering, 1995, 121(1): 23-35.

[18] 柳淑学. 多向不规则波有效试验区的扩展法[J]. 海洋工程, 1996, 14(2): 40-50.

Liu S X. Extension of the effective test area of multi-directional waves by side wall reflections [J]. Ocean engineering, 1996, 14(2): 40-50.

StudyonWaveDiffractionCoefficientsofArtificialIslandwithDifferentScales

YU Ding-Yong， LI Long， ZHANG Guang-Cheng

(College of Engineering，Ocean University of China，Qingdao 266100，China)

Artificial island is a large structure in coastal engineering. The determination of design wave height has a close relation with wave diffraction. Although theCodeofHydrologyforSeaHarbor(JTS145-2-2013) proposes a method to calculate the height of diffracted wave in shield region of coastal structure, theCode’s method does not consider the effect of structure scales along wave propagation direction. In fact, the wave diffraction coefficients differ in the shield area due to the change of structure scale. In this paper, in order to study the wave diffraction coefficients of artificial island with different scales, MIKE21-BW model is employed to simulate the wave field with different structure scales. The isolines of diffraction coefficients in the shield area with different scales of artificial island are drawn based on the simulated results. The wave diffraction coefficients of artificial island with different scales are compared with experimental data as well as the results obtained by theCode’s method. The effects of structure scales to diffraction coefficients of the shield region are then analyzed.

MIKE21-BW model is an efficient model to study the combined effects of most wave phenomena, including shoaling, refraction, diffraction, partial refection and transmission and so on. The model includes two types of the Boussinesq equations, i.e., the enhanced Boussinesq equations and the classical Boussinesq equations. In this paper, the enhanced Boussinesq equation is used.

In simulating, conditions of the incident wave are supposed the same as Briggs’ physical experiment and theCode’s, respectively. The values of diffraction coefficients obtained by simulation agree well with experimental data, while, the values obtained by simulation with no consideration of structure scales along wave propagation direction differ quite small from results by theCode’s method. The comparison results can verify the applicability of BW model for this study.

In order to study the effect of structure scales along wave propagation, this paper simulates 7 cases of different structure scale, i.e., the length of artificial island along wave propagation is 0.5LsLs, 2Ls, 3Ls4Ls5Lsand 6Ls, respectively.

The following results are then obtained from the simulation:

(1) Diffraction coefficients become bigger with the increase of scales from 0.5Lsto 2Ls;

(2) Diffraction coefficients become slightly smaller with the increase of scales from 3Lsto 6Ls;

(3) Diffraction coefficients vary relatively obvious in the area near island and island-shielded region;

(4)The scales of artificial island has definite effect on values determination of diffraction coefficients.

In this paper, the isolines of diffraction coefficients with different island scales are given, which can be references for calculating the diffracted wave height.

BW model; artificial island; diffraction coefficients; scales; numerical simulation

研究简报

TV143

A

1672-5174(2018)02-096-09

10.16441/j.cnki.hdxb.20150314

于定勇, 李龙， 张广成. 不同尺度人工岛波浪绕射系数的研究[J]. 中国海洋大学学报(自然科学版), 2018, 48(2): 96-104.

YU Ding-Yong， LI Long， ZHANG Guang-Cheng. Study on wave diffraction coefficients of artificial island with different scales[J]. Periodical of Ocean University of China, 2018, 48(2): 96-104.

山东省自然科学基金项目 (ZR2013EEZ002)资助

Supported by Natural Science Foundation of Shandong Province(ZR2013EEZ002)

2015-10-12；

2016-09-20

于定勇(1964-)，男，教授。 E-mail: dyyu01@ouc.edu.cn

责任编辑 陈呈超