箱式挡浪结构透空堤消浪性能数值模拟_桂劲松.pdf

箱式挡浪结构透空堤消浪性能数值模拟桂劲松，夏曦（大连海洋大学，辽宁 大连 116023）摘要：本文基于CFD方法对箱式挡浪结构透空堤的消浪性能开展数值模拟研究，对比箱式和立板式结构的透射、反射及粘性流引起的耗散系数，重点分析波浪经过箱式结构时不同水平底板位置对透射、反射及耗散的影响、及波浪经过箱式挡浪结构时流场变化情况。研究结果表明，较之传统的立板式挡浪结构，箱式结构具有较小的透射系数、较大的反射系数；分析粘性流引起的耗散发现，立板式挡浪结构主要由波浪震荡引起耗散，箱式挡浪结构主要由平板约束水质点的竖向运动引起耗散；分析流场图得到，当波浪经过箱式挡浪结构时大部分被反射，少部分透射进入的波浪受到水平底板的约束，波高减小且趋于平稳，达到更好的消浪效果。关键词：箱式挡浪结构；OpenFoam；透射系数；反射系数；流场中图分类号：TV131.2文献标识码：Adoi:10.3969/j.issn.1007-7294.2023.02.007Numerical simulation of wave dissipation performance ofbox-type wave-retaining structure permeable breakwaterGUI Jin-song,XIA Xi(Dalian Ocean University,Dalian116023,China)Abstract:In this paper,the numerical simulation of wave dissipation performance of a box-type wave-retaining structure permeable breakwater was carried out based on CFD method.The transmission,reflection anddissipation coefficients caused by viscous flow of box-type and vertical plate structures were compared.Theinfluence of different horizontal bottom plate positions on transmission,reflection and dissipation at the timewhen waves passed through box-type structure and the change of flow field at the time when waves passedthrough box-type wave retaining structure were analyzed.The results show that the box structure has a smaller transmission coefficient and a larger reflection coefficient than those of the traditional vertical plate structure,the dissipation caused by viscous flow is mainly caused by the wave oscillation in the vertical platestructure,and the dissipation is caused by the vertical movement of the water point restrained by the plate inthe box structure.Analysis of the flow field diagram shows that,when the waves pass through the box-typewave-retaining structure,most of them are reflected,and only a small part of the transmitted waves are constrained by the horizontal bottom plate,and that the wave height decreases and tends to be stable,resulting ina better wave dissipation effect.Key words:box-type wave retaining structure;OpenFoam;transmission coefficient;reflection coefficient;flow field第27卷第2期船舶力学Vol.27 No.22023年2月Journal of Ship MechanicsFeb.2023文章编号：1007-7294（2023）02-0227-11收稿日期：2022-08-29基金项目：辽宁省重点研发计划项目（2018228003）作者简介：桂劲松（1968-），男，博士，教授，E-mail:；夏曦（1997-），女，硕士研究生，通讯作者，E-mail:。0 引言当前，各国越来越重视海洋生态环境的保护。相较于传统的重力式直立堤1、抛斜坡堤2等结构对海洋生态造成的破坏3，透空式防波堤4由于其下部透水，便于港内外水体交换，在海洋环境保护方面具有明显的优势，目前已逐渐发展成为一种重要的海洋工程结构。最早对透空堤的研究主要针对单挡板透空堤。上世纪50年代，Uresll5研究了无限水深情况下单挡浪板透空堤的透射效果，该理论研究对挡板式透空堤的发展有着深远影响。单侧挡浪板结构简单，易于施工，然而对于波浪的掩护作用较为有限。近年来许多学者对不同挡浪板结构开展了研究，周效国等6根据波能分布位置研究了不同层数和不同高程的开孔多层直立式挡板的消浪特性，总结出消浪效果好的挡板布置方案；Chioukh等7考虑正向入射的规则波，分析了双挡浪板透空堤的透射和反射系数；杜沛霖等8研究了带有透浪通道的直立式透空堤的消浪特性，由于波浪在透浪通道中多次的反射及绕射消耗了大量波能，达到了良好的消浪效果；Koraim等9研究了双侧直立式挡浪板透空堤在波长以及结构尺寸变化时的水动力特性，并且建立了基于特征函数的理论模型；邵杰等10通过物理模型试验，研究了不同结构型式的双垂直挡浪板透空堤的透射系数，并指出透射系数随挡浪板入水深度的增大会显著减小；蔡丽等11研究了V形挡板式透空堤透射后的波高分布情况，结果表明，V形挡板透空堤对短周期波有很好的掩护作用；范骏等12针对双挡浪板透空堤进行了物理模型实验，不仅考虑了挡浪板入水深度和水平板板宽等因素对透射系数的影响，还对透浪系数的计算方法进行了研究；Li等13研究了水平单弧形板的波浪透射、反射、波压以及弧板周围流场变化；徐宁等14探究了双层水平板透空堤的透浪效果，发现相对板宽是影响结构透射和反射的主要因素；王丽雪等15对比了双层水平板和双层弧形板的透射效果，发现弧形板各方面性能均优于水平板。以上学者大多是对水平或垂直板透空堤消浪结构进行了研究，但对于水平板和垂直板组合的箱式挡浪结构透空堤的研究却鲜有报道。本文基于粘性流理论对比了垂直立板式透空堤和箱式透空堤的消浪特性，并且深入探究了箱式透空堤的水动力特性，揭示了箱式透空堤的消浪机制。1 数值模型建立本文基于OpenFoam求解器建立数值波浪水槽，采用有限体积法对空间离散，应用PISO算法迭代求解压力及速度，应用动边界法模拟推板造波，在数值波浪水槽的末端设置阻尼消波区。1.1 控制方程U=0(1)Ut+()UU-()U+=CK-gX-Pgh+S*(2)式中，U为平均速度矢量，为二相流的平均密度，t为时间，为二相流平均动力粘度，为雷诺应力项，其中表示由脉动值引起的雷诺应力张量。CK为表面张力项，其中C为表面张力系数，一般取为0.07 kg/s2，K为自由面的曲率，为相体积分数，g为重力加速度，X为位置矢量。Pgh是为了求解方便引入的压力项，S*为动量消波源项。1.2 VOF 法追踪自由表面对不可压气液二相流模型，VOF法可根据两相流在网格中所占比例来确定自由面流体变化，该流体的密度和动力粘度系数由以下体积分数函数表示：=1+()1-2=1+()1-2(3)式中，1为水的密度，2为空气的密度，1为水的动力粘度系数，2为空气的动力粘度系数。228船舶力学第27卷第2期1.3 边界条件本模型主要是固壁边界条件，速度场在边界满足不可滑移条件：|Uwall=0(4)压力场在壁面处满足压力沿法向方向的梯度为0的条件：pn=0(5)1.4 造波与消波本文采用动边界法来模拟实验室的推板造波，根据线性理论得出规则波波高与推板冲程的关系为Hs=HiS=4sinh2()kd2kd+sinh2kd(6)式中，Hi为入射波高，S为造波机造波板冲程，k为波数，d为水深。推板的位移X和速度U分别为X()t=S2sin()tU()t=S2cos()t(7)根据以上推板的位置和速度公式可以得到推板线性波波面表达式为()x,t=Hi2cos()kx-t(8)式中，是波浪频率，x表示推板位置。消波方法由Larsen等16提出，在动量方程中添加源项，使波浪在阻尼区中逐渐衰减。S*=u(9)式中，S*为动量消波源项，u为速度矢量，为衰减函数。本文中衰减系数取为5，阻尼区长度取为12倍波长。2 数值模型验证2.1 理论验证数值波浪水槽利用项目组前期研发的基于OpenFoam的DUTFoam求解器17模拟，尺寸为：长21.7m、宽0.7 m、高0.8 m，水深为0.3 m，波浪参数为：波高Hi=0.08 m，周期T=1.8 s，阻尼区取2倍波长。水槽中共放置四个浪高仪，位置分别在2倍波长、3倍波长、1/2阻尼区和水槽末端处放置，将其依次编号为 WG1、WG2、WG3和WG4，图1为数值水槽示意图。首先本文对数值模型的网格收敛性进行了验证，选择四种尺度的网格计算了波面变化，网格参数如表1所示。表1 网格参数Tab.1 Mesh parameters编号1234x/m0.0320.0160.0100.004z/m0.0160.0080.0050.004Nx723146421692892Nz50100150200网格总数36 150146 400325 350578 400图1 数值水槽示意图Fig.1 Schematic diagram of numerical flume第2期桂劲松等：箱式挡浪结构透空堤消浪性能 229图2为空水槽时在结构物处的四种不同网格尺度所计算的波面结果，由图可知当网格尺度为0.01 m0.005 m0.005 m（长宽高）和0.008 m0.004 m0.004 m（长宽高）时波面已平稳，表明网格已收敛。考虑到计算效率等问题，本文选取0.01 m0.005 m0.005 m作为网格尺度。图3为波面模拟结果和微幅波理论结果对比。由图可知，在x=5.8 m和x=8.7 m位置的波面和理论波面吻合较好。图3（c）（d）是阻尼区中间处水槽末端波面与理论波面的对比图，可以看出水槽末端的阻尼消波区具有较好的波浪衰减作用。综上所述，理论上数值水槽具有有效性。图2 不同网格尺度的波面历时曲线比较Fig.2 Comparison of wave surface duration curves with different grid scales（a）x=5.8 m（b）x=8.7 m（c）x=16.6 m（d）x=21.6 m图3 微幅波波面和数值模拟波面对比图Fig.3 Comparison of airy wave surface and numerical simulation wave surface2.2 物模验证在结构物前后各放置两个波高仪测量入射、反射和透射波高，计算透射系数Kt和反射系数Kr。Kt=HtHi,Kr=HrHi(10)式中，Ht为透射波高，Hi为入射波高，由Goda两点法18分离出反射波高Hr。物理模型主要由前挡浪板和后挡浪板组成，前后挡浪板的固定尺寸为0.7 m0.3 m0.025 m。试验采用波浪周期T=1.8 s，波长L=2.88 m，Hi=0.06 m，采用DS30波高水位测量仪和DJ800多功能检测系统测量和采集波面数据。图4为数模试验时的模型简图。在进行物模试验和数模试验对比验证时，采用图4中无水平底板的结构，物模试验中采用的小直径桩对流场等水动力影响较小，数模时忽略。图5为物理模型试验的布置简图，水槽长为40 m，结构物放置在水槽中后方，在结构物前方2 m和3 m处各放置两个浪高仪，在结构物后方1 m和2 m处各放置两个浪高仪。网格1网格3网格2网格4i/m微幅波数值模拟数值模拟数值模拟微幅波微幅波微幅波数值模拟i/mi/mi/mi/m230船舶力学第27卷第2期图4 数值水槽模型简图Fig.4 Numerical flume model diagram图5 物理试验布置图Fig.5 Layout of physical test图6是T=1.8 s、d=0.3 m，Hi=0.06 m时物模波面与数模波面入射波、反射波和透射波的对比结果。反射波为WG1和WG2测得的合成波采用两点法分离后的波面，透射波为WG3的波面，WG3测得的透射波和WG4测得的透射波基本一致，故在本文中只选择了WG3波面进行对比。（a）入射波波面对比（b）反射波波面对比（c）透射波波面对比图6 波面对比图Fig.6 Wave surface contrast diagram（a）透射系数对比（b）反射系数对比图7 透射系数和反射系数物模数模对比Fig.7 Modulo-digital comparison of transmission coefficient and reflection coefficient物理试验波面数值模拟波面物理试验波面数值模拟波面数值模拟波面物理试验波面/d/di/mr/mt/m第2期桂劲松等：箱式挡浪结构透空堤消浪性能 231图7为d/L=0.1、0.17、0.21（d为水深，恒为0.3 m，L为波长），波高Hi=0.06 m时，前挡浪板入水深度0.5Hi，后挡浪板相对入水深度不同时，物理试验和数值模拟的透射系数和反射系数对比图，横坐标为相对入水深度/d（为后挡板入水深度）。由图67可见，波面曲线、透射系数、反射系数的数模和物模结果基本吻合，验证了数值水槽的有效性。3 数值模型结果与讨论3.1 垂直挡板和无底板箱式结构对比图8为数值模拟的结构网格图。（a）单挡板结构（b）双垂直挡板结构（c）无底板箱式结构图8 结构网格图Fig.8 Structural grid diagram图9比较了Hi/L（Hi为入射波高）为0.0330.073，Hi/d为0.20.5时3种不同结构挡板透射系数变化趋势。结果表明：单挡板结构透射系数最大，双垂直挡板结构透射系数小于单挡板，波高变化时表现（a）不同结构透射系数随波陡的变化（b）不同结构透射系数随相对波高的变化图9 透射系数变化图Fig.9 Diagram of transmission coefficient changes（a）不同结构反射系数随波陡的变化（b）不同结构反射系数随相对波高的变化图10 反射系数变化图Fig.10 Diagram of reflection coefficient changes单挡板结构双垂直挡板结构无底板箱式结构单挡板结构双垂直挡板结构无底板箱式结构单挡板结构双垂直挡板结构无底板箱式结构单挡板结构双垂直挡板结构无底板箱式结构Hi/LHi/dHi/dHi/L232船舶力学第27卷第2期不明显，波长变化时表现较为明显，最大可减小45%，最小可减小12%。无底板箱式结构透射系数在3种结构中最小，对短周期入射波消浪效果更好，比双垂直挡板透射系数最大可消减35%。图10比较了3种不同结构挡板反射系数变化情况。结果表明：单挡板结构反射系数最小，无底板箱式结构反射系数最大，其波高不同时反射系数增大得更为明显，最大可增加70%，最小可增加60%。无底板箱式结构在半封闭区域内形成了消浪室，其消浪室内上面为空气，下面为水。图11为同一位置双垂直板内和无底板箱式结构消浪室内气压变化情况，由图可见，消浪室内波浪受到气压的影响，压强变化明显。图12为同一位置双垂直板内和消浪室内波面位移变化情况，由图可见，当波浪经过无底板箱式结构时，消浪室内波浪受到气圧影响，波面位移明显较小。图11 气压随时间变化Fig.11 Pressure varying with time（a）双垂直挡板结构内波面位移（b）无底板箱式结构内波面位移图12 结构内波面位移随时间变化Fig.12 Wave surface displacement varying with time耗散系数是衡量波能消耗的系数之一，波能耗散系数定义如下：K2d=1-K2t-K2r(11)式中，Kt为透射系数，Kr为反射系数，Kd为耗散系数图13比较Hi/L为0.0330.073，Hi/d为0.20.5时3种不同结构挡板耗散系数。结果表明：双垂直挡板耗散系数较大，无底板箱式结构耗散系数较小，在波陡较大和相对波高较大时其耗散相差较小。由图12可见，两种结构内波高相差较大，双垂直挡板和单挡板由于无水平顶板约束，不受消浪室内气压的影响，波浪上下震荡幅度较大，耗能相对较大；增加水平顶板，波浪受气压影响，上下震荡幅度变小，能量不能充分消耗，耗能相对较小。耗散系数随波陡的增加而减小，随相对波高的增大而增大。无底板箱式结构透射系数较小，但能量耗散较差，所以在结构中增加了水平底板，探究箱式透空堤的消浪特性。（a）不同结构耗散系数随波陡的变化（b）不同结构耗散系数随相对波高的变化图13 耗散系数变化图Fig.13 Diagram of dissipation coefficient changes双垂直板间气压无底板箱式内气压单挡板结构双垂直挡板结构无底板箱式结构单挡板结构双垂直挡板结构无底板箱式结构Hi/LHi/dt/mt/mP/Pa第2期桂劲松等：箱式挡浪结构透空堤消浪性能 2333.2 箱型板式透空堤水动力分析图14为数值模拟的箱型板式结构网格图。（a）平面图（b）网格图图14 箱型板式结构图Fig.14 Structural diagram of box-type structure在结构中增加水平底板，上水平顶板和水平底板中间为空心结构，水平板厚度均为0.02 m,结构宽度B=0.7 m。在水面上方1.0Hi（入射波高）、水面上、水下0.5Hi和水下1.0Hi的位置增加水平底板。图15比较了B/L为0.1750.4时水平底板在不同位置，透射系数和反射系数随相对板宽（B/L）的变化。结果表明,相对板宽对结构透射系数和反射系数影响显著，相对板宽为0.190.4时透射系数随相对板宽的增加而减小，但同时反射系数随相对板宽的增加而增大。对比不同结构，随着水平底板高度的降低透射系数随之减小，当底板放置在水面和在空气中时透射系数减小，幅度相对较小，底板在水下0.5Hi和1.0Hi时透射系数减小，幅度相对较大，而底板在水下1.0Hi时透射系数最小，与无底板结构透射系数相比最大可减小20%。反射系数随水平板高度的降低而增大，底板在水下1.0Hi时反射系数最大。综合可见，增加底板后的箱式透空堤消浪效果对长周期波和短周期波都适用。（a）箱式结构透射系数随相对板宽的变化（b）箱式结构反射系数随相对板宽的变化图15 透射系数、反射系数随相对板宽的变化Fig.15 Transmittance coefficient and reflection coefficient varying with relative plate width图16比较了Hi/d为0.20.5时水平底板不同位置处透射系数和反射系数随相对波高的变化。结果表明，相对波高对结构透射系数和反射系数的影响较为显著，透射系数和反射系数随相对波高的增（a）箱型板式透空堤透射系数随相对波高的变化（b）箱型板式透空堤反射系数随相对波高的变化图16 透射系数、反射系数随相对波高的变化Fig.16 Transmittance coefficient and reflection coefficient varying with relative water height无水平底板底板在水面上1.0Hi底板在水面底板在水下0.5Hi底板在水下1.0Hi无水平底板底板在水面上1.0Hi底板在水面底板在水下0.5Hi底板在水下1.0Hi无水平底板底板在水面上1.0Hi底板在水面底板在水下0.5Hi底板在水下1.0Hi无水平底板底板在水面上1.0Hi底板在水面底板在水下0.5Hi底板在水下1.0HiHi/dHi/dB/LB/L234船舶力学第27卷第2期加而减小。对比不同结构，随着水平底板高度的降低透射系数随之减小，底板在水下0.5Hi和1.0Hi时透射系数减小幅度相对较大，而底板在水下1.0Hi时透射系数最小，相对于无底板结构透射系数最大可减小15%。反射系数随水平板高度的降低而增大，底板在水下1.0Hi时反射系数最大。综合可见，增加底板后的箱式结构透空堤对大的波高掩护效果更好。图17比较了B/L为0.1750.4，Hi/d为0.20.5时水平底板不同位置耗散系数的变化。结果表明，增加水平底板且底板放置在水中，结构耗散系数较大。水平底板放置在水中时会约束水质点的竖向运动，消耗更多能量。（a）箱型板式透空堤耗散系数随相对堤宽的变化（b）箱型板式透空堤耗散系数随相对波高的变化图17 耗散系数变化图Fig.17 Diagram of dissipation coefficient changes水平底板放置在水中的位置对耗散系数影响较小，综合透射、反射和耗散系数，水平底板放置在水下0.5Hi和1.0Hi时透射较小，反射和耗散较大，消浪效果更为显著。3.3 流场对比图1819分别为周期T=1.8 s，无底板箱式结构和箱式结构在水下1.0Hi时的结构流场变化图。由图可见，对于无底板结构，波浪入射时，部分波浪进入消浪室内，虽有气压作用但消浪室内水质点竖向运动较为显著，使得透射波相对较大。增加水下水平板后，大部分波浪受到前挡板和水平板阻挡形成立波产生反射，透射的小部分波浪由于水平底板约束了水质点的竖向运动，从而达到良好的消浪效果。10120.00.30.6t=0.0Tz/mx/m10120.00.30.6t=1.0Tz/mx/m图18 一个周期内无底板箱式结构流场变化图Fig.18 Flow field variation diagram of box structure without bottom plate in one cycle10120.00.30.6t=0.0Tz/mx/m无水平底板底板在水面上1.0Hi底板在水面底板在水下0.5Hi底板在水下1.0Hi无水平底板底板在水面上1.0Hi底板在水面底板在水下0.5Hi底板在水下1.0Hit=T/4t=T/4t=3T/4B/LHi/dKdKd第2期桂劲松等：箱式挡浪结构透空堤消浪性能 23510120.00.30.6t=1.0Tz/mx/m图19 一个周期内箱式结构流场变化图Fig.19 Flow field diagram of box structure in one cycle4 结论本文通过数值模拟对箱式挡浪结构透空堤消浪特性进行了系统的研究，主要结论如下：（1）箱式挡浪结构透空堤透射系数随相对板宽的增加而减小，反射系数则相反。透射系数和反射系数均随相对波高的增加而减小。（2）箱式结构水平底板位于水面以上时，由于消浪室内气压的影响产生耗散，耗散系数相对较小；水平底板位于水面以下时，由于水平底板约束水质点竖向运动产生耗散，耗散系数相对较大，消浪效果更好。（3）箱式结构透射系数随水平底板高度的降低而减小，反射系数相反。参考文献：1 Madabhushi G S P,Boksmati J I,Torres S G.Numerical and centrifuge modeling of gravity wharf structures subjected to seismic loadingJ.Journal of Waterway,Port,Coastal,and Ocean Engineering,2020,146(4):04020007.2 Vijay K G,Venkateswarlu V,Karmakar D.Scattering of gravity waves by multiple submerged rubble-mound breakwatersJ.Arabian Journal for Science and Engineering,2020,45:1-22.3 Rageh O S,Koraim A S.Hydraulic performance of vertical walls with horizontal slots used as breakwaterJ.Coastal Engineering,2010,57(8):745-756.4 陈伟毅,林伊楠.透空式防波堤之特性研究综述J.中国水运（下半月）,2016,16(6):207-209.Chen Weiyi,Lin Yinan.A review on the characteristics of permeable breakwatersJ.China Water Transport(Second Half),2016,16(6):207-209.(in Chinese)5 Ursell F.The effect of a fixed vertical barrier on surface waves in deep waterC/Mathematical Proceedings of the CambridgePhilosophyical Society,1947,43:374-382.6 周效国,曹凤龙,李丁丁.多层直立挡板透空式防波堤消浪效果数值研究J.江苏科技大学学报（自然科学版）,2020,34(5):15-20.Zhou Xiaoguo,Cao Fenglong,Li Dingding.Numerical study on wave dissipation effect of multi-layer vertical baffle permeable breakwaterJ.Journal of Jiangsu University of Science and Technology(Natural Science Edition),2020,34(5):15-20.(in Chinese)7 Chioukh N,Ouazzane K,Yksel Y,Hamoudi B,evik E.Meshless method for analysis of permeable breakwaters in theproximity of a vertical wallJ.China Ocean Engineering,2019,33(2):148-159.8 杜沛霖,孙昭晨,梁书秀.新型直立式透空堤消浪性能数值研究J.海洋工程,2021,39(1):12-20.Du Peilin,Sun Zhaochen,Liang Shuxiu.Numerical study on wave dissipation performance of new vertical permeable breakwaterJ.Ocean Engineering,2021,39(1):12-20.(in Chinese)9 Koraim A S,Heikal E M,Rageh O S.Hydrodynamic characteristics of double permeable breakwater under regular wavesJ.Marine Structures,2011,24(4):503-527.10 邵杰,陈国平,严士常,等.不规则波作用下垂直挡浪板式透空堤透浪系数试验研究J.海洋工程,2016,34(1):50-57.t=3T/4236船舶力学第27卷第2期Shao Jie,Chen Guoping,Yan Shichang,et al.An experimental study of the wave-permeating coefficient of vertical bafflepermeable breakwater under irregular wavesJ.Ocean Engineering,2016,34(1):50-57.(in Chinese)11 蔡丽,王永学,王国玉,徐杰,赵东梁.V形挡板式透空堤后的波高分布数值分析J.水动力学研究与进展（A辑）,2018,33(2):188-198.Cai Li,Wang Yongxue,Wang Guoyu,Xu Jie,Zhao Dongliang.Numerical analysis of wave height distribution behind V-shaped baffle breakwaterJ.Hydrodynamic Research and Progress(Series A),2018,33(2):188-198.(in Chinese)12 范骏,王宇楠,杨斯汉,等.双挡浪板透空堤透浪与反射系数实验研究J.海洋工程,2011,29(4):60-67.Fan Jun,Wang Yunan,Yang Sihan,et al.Experimental research on reflection coefficient and transmitted coefficient of double barriers penetrated breakwaterJ.Ocean Engineering,2011,29(4):60-67.(in Chinese)13 Li Xueyan,Li Qin,Wang Qing,Hou Chengyi,et al.Numerical and experimental investigation on the hydrodynamic characteristics of an arc-shaped plate-type breakwater under the action of long-period wavesJ.Ocean Engineering,2020,219:108198.14 徐宁,杨小桦,庞丹,刘仔峰.透空式双层板防波堤消波特性研究J.中国水运（下半月）,2014,14(7):214-216.Xu Ning,Yang Xiaohua,Pang Dan,Liu Zifeng.Study on wave dissipation characteristics of permeable double deck breakwaterJ.China Shipping(Second Half),2014,14(7):214-216.(in Chinese)15 王丽雪,李雪艳,王庆,等.双弧板式透空堤消浪特性物理模型试验J.厦门大学学报（自然科学版）,2020,59(4):534-539.Wang Lixue,Li Xueyan,Wang Qing,et al.Physical model test on wave dissipation characteristics of double arc plate permeable breakwaterJ.Journal of Xiamen University(Natural Science Edition),2020,59(4):534-539.(in Chinese)16 Larsen J,Dancy H.Open boundaries in short wave simulations-A new approachJ.Coastal Engineering,1983,7:285-297.17 Wang Dongxu.Development and application of numerical wave flume based on OpenfoamD.Dalian:Dalian Ocean University,2018.18 Somervell L T,Thampi S G,Shashikala A P.Estimation of friction coefficient for double walled permeable vertical breakwaterJ.Ocean Engineering,2018,156:25-37.第2期桂劲松等：箱式挡浪结构透空堤消浪性能 237