第39卷第1期2021年1月海洋科学进展A D V A N C E S I N MA R I N E S C I E N C E V o l .39 N o .1J a n u a r y ,2021研究论文潜堤上规则波辐射应力的数值研究穆浩然1,2,马小舟1,2*,毛艳军1,2,高 翔1,2(1.大连理工大学建设工程学部,辽宁大连116024;2.海岸和近海工程国家重点实验室,辽宁大连116024)收稿日期:2020-03-22资助项目:国家重点研发计划项目 港湾低频振荡精细化预测与港口安全评估技术(2017Y F C 1404205);中央高校基本科研业务费资助 深海岛礁水动力及港湾震荡研究(D U T 18Z D 214)作者简介:穆浩然(1994 ),女,硕士研究生,主要从事低频波浪方面研究.E -m a i l :h a o r a n .m u @q q .c o m *通信作者:马小舟(1977 ),男,副研究员,博士,主要从事港湾振荡和低频波浪方面研究.E -m a i l :m a x z h @d l u t .e d u .c n (王佳实 编辑)摘 要:在近岸波浪相关研究中,辐射应力是波动在水体中引起的剩余动量流,是波浪运动的重要物理量㊂在波浪从深水逐渐传向浅水的过程中,波浪的非线性逐渐增强,甚至会发生破碎等剧烈变形,引起辐射应力的强烈变化,对次重力波生成等有重要贡献㊂应用O p e n F O AM 精细模拟波浪在潜堤上的传播,得出波浪运动的详细流场信息,计算了有波浪破碎情况下潜堤地形上波浪的辐射应力和波浪增减水情况㊂研究结果表明,在潜堤地形下,辐射应力值在堤前平底处受波幅变化影响较敏感,波浪完全破碎后其值在堤后坡面处呈增大趋势直至平底处趋于平稳㊂辐射应力沿程变化对波浪增减水的影响趋势与波浪增减水方程符合良好㊂关键词:波浪辐射应力;波浪增减水;潜堤;O pe n F O AM 中图分类号:T V 131.2 文献标识码:A 文章编号:1671-6647(2021)01-0158-09d o i :10.3969/j .i s s n .1671-6647.2021.01.016引用格式:MU H R ,MA XZ ,MA O YJ ,e t a l .N u m e r i c a l s t u d y of r a d i a t i o n s t r e s s o n t h e s u b m e rg e db a r s b y r e gu -l a rw a v e [J ].A d v a n c e s i nM a r i n e S c i e n c e ,2021,39(1):158-166.穆浩然,马小舟,毛艳军,等.潜堤上规则波辐射应力的数值研究[J ].海洋科学进展,2021,39(1):158-166.辐射应力概念最初是由L o n g u e s t -H i g gi n s 和S t e w a r t [1]提出并给出了适用于计算微幅波辐射应力值的简化公式㊂在近岸波浪的传播运动中,由于海水深度变浅导致了波浪产生严重变形甚至破碎,从而引起辐射应力的强烈改变,辐射应力变化对沿岸波生流㊁波浪增水(s e t -u p)㊁波浪减水(s e t -d o w n )有着重要影响,也是低频波浪生成的重要驱动力㊂L o n g u e s t -H i g gi n s 和S t e w a r t [2]认为波浪辐射应力的增加将促使波群中自由长波的释放;S y m o n d s 等[3]建立了控制方程中含有辐射应力受迫项的碎波拍模型,认为波群中波浪破碎点前后移动引起辐射应力在时间及空间上的改变产生了低频波浪;K o s t e n s e [4]定性地验证了短波群波浪破碎点移动导致低频波浪产生的理论㊂因此,为了更加准确地了解辐射应力这一波浪特性,众多学者通过理论分析,数值和实验模拟的手段进行辐射应力的计算㊂对于常水深G r a h a m [5]提出用有限差分法计算辐射应力张量;郑金海和严以新[6]利用线性波理论研究波浪辐射应力张量随深度变化的分布规律,并给出适用于任意波向角的辐射应力的表达式;X i a 等[7]提出了辐射应力垂向剖面的定义和计算公式,并利用斜坡上波生流的实验结果证明了模型的准确性;S t i v e 和W i n d [8]利用实验数据,通过线性外推法近似求得湍流波浪速度剖面来计算辐射应力沿程变化;温秀媛等[9]利用改进色散关系的B o u s s i n e s q 方程推导出了全新的辐射应力计算公式,探究辐射应力对波浪增减水的影响㊂然而,上述辐射应力计算方法都难以获得碎浪带内详细的流场信息,并存在着一定程度的条件假设和公式简化,因此计算所得辐射应力结果误差较大㊂基于O pe n F O AM [10]软件开源程序建立数值水槽则是模拟波浪破碎区流场信息的一个较为可行的手段,可以详细研究完整计算域内波浪辐射应力随波浪传播而产生的变化过程㊂查晶晶和万德成[11]基于1期穆浩然,等:潜堤上规则波辐射应力的数值研究159O p e n F O AM 软件进行的数值造波实验验证了此数值方法造波和消波方式的可靠性;王东旭等[12]利用此软件较好地模拟了孤立波在潜礁上的波浪破碎及水跃现象;毛艳军等[13]应用O pe n F O AM 软件求解了箱式浮式防波堤的水动力性能,并分析了其作为垂荡浮子式波能转换装置的能量转换性能;姚宇等[14]运用O pe n -F O AM 软件模拟孤立波在岛礁地形的传播,数值输出的破碎区内波形及流速与实验结果对比良好;张陈浩和郑茜[15]利用O p e n F O AM 软件模拟规则波在浅滩上的破碎变形,模拟结果与实验结果相吻合,且输出了能较好体现波浪破碎整体过程的流场信息㊂本文拟应用O pe n F O AM 软件来模拟波浪在潜堤地形下的传播破碎,在进行平底水槽线性波的数值模拟和数值模型模拟波浪的准确性及计算辐射应力方法的有效性验证后,进行规则波在潜堤地形下传播破碎的模拟研究,在此基础上对辐射应力沿程变化的特点及其对波浪增减水的影响开展了研究,并探讨了平底地形下辐射应力与波幅偏移量之间的关系㊂1 理论基础1.1 控制方程本文采用基于O p e n F O AM 软件的造波模块w a v e s 2F o a m [16]开展波浪传播变形的数值模拟㊂坐标系统为笛卡尔坐标系,x 方向为波浪的入射方向,y 方向为垂向,z 方向为沿水槽的宽度方向,基于不可压缩流体的运动特点,控制方程为:Ñ㊃u =0,(1)∂ρu ∂t +Ñ㊃(ρu u )=-Ñp +Ñ㊃μÑu -g ㊃x Ñρ+f σ,(2)f σ=σκÑα,(3)n =Ñα|Ñα|,(4)κ=-Ñ㊃n ,(5)式(1)~式(5)中,ρ为混合流体的密度;u 为流体的速度场;t 为时间;μ为动力黏性系数;g 为重力加速度;x 为位置矢量;f σ为自由表面的张力;σ为张力系数;κ为自由表面的平均曲率;n 为界面单位法向量;α为流体体积分数㊂对于自由表面的处理,根据V O F 法将2种流体看作一种混合流体,实现每个单元相界面的追踪,引入相函数:α(x ,t )=1液相 0<α<1 气液混合相0气相 ìîíïïïï,(6)在求解时,式(2)中的混合流体的密度ρ可用液体密度ρ1和气体密度ρ2表示,动力黏性系数μ可用液体黏性系数μ1和气体黏性系数μ2来表示:ρ=αρ1+(1-α)ρ2,(7)μ=αμ1+(1-α)μ2㊂(8)同时,α满足输运方程:D αD t =∂α∂t +Ñ㊃(αu )=0㊂(9)传统的V O F 法需要反复地进行自由表面的重构,降低了计算效率,因此W e l l e r [17]引入了额外的人工压缩项来实现对自由水面的精确捕捉,式(9)可变为∂α∂t +Ñ㊃(αu )+Ñ㊃[α(1-α)u c ]=0,(10)160海洋科学进展39卷式中:u c为适用于压缩自由面及调整自由面尖锐程度的速度场,新增加的人工压缩项只对界面过渡区起作用,不会影响到周围流场㊂1.2辐射应力公式将作用于单位面积水柱体的总动量流的时均值减去没有波浪作用时的静水压力定义为波浪辐射应力㊂L o n g e s t-H i g g i n s和S t e w a r t[18]给出的水平x方向辐射应力(S x x)的定义公式为S x x=ʏη-h(p+ρu2)d z-ʏη-hρg(η-z)d z,(11)式中:u为质点速度在x方向的分量;h为水深;η为水面高程;i为平均水面高程;p为总压㊂对于微幅波,辐射应力计算公式可简化为S x x=E(2n-12),(12)n=12(1+2k hs h2k h)㊂(13) 2模型验证2.1波形验证将模型作用于模拟平底水槽中线性波的传播,其模拟工况为:波高0.02m,波周期2s,水深1m,计算域长30m,高1.1m,数值水槽前端和后端各设置6m松弛区(消波),模拟计算时长30s,采用自动调整计算时间步长的方法,设置输出时间步长Δt=0.02s㊂模型应用C a s e1,C a s e2和C a s e3三种不同大小的网格设置来进行敏感性分析,具体网格尺寸见表1㊂不同网格尺寸下波面对比如图1所示,由图可见,C a s e1网格尺度表1不同工况下的网格尺寸T a b l e1 S e n s i t i v i t y a n a l y s i s o f g r i d s i z e模拟工况水平方向网格尺寸/m垂直方向网格尺寸/m总网格数/万C a s e10.0300.01011.00C a s e20.0150.00544.00C a s e30.0120.00468.75图1不同网格尺寸下波面变化与解析波面对比F i g.1 C o m p a r i s o nd i a g r a mo fw a v e s u r f a c eb e t w e e ne x p e r i m e n t a l r e s u l t s a n dn u m e r i c a l s i m u l a t i o nw i t hd i f f e r e n t g r i d s i z e1期穆浩然,等:潜堤上规则波辐射应力的数值研究161 下模拟的波浪波峰值明显低于解析结果,由此产生的误差将影响到之后的数值计算,而C a s e 2和C a s e 3网格尺度下模拟的波浪则与理论解对应很好㊂因此,为保证波浪辐射应力计算结果的准确性并且能够节约计算资源,本文采用C a s e 2网格尺寸进行模拟计算㊂得到了较好的数值模拟结果,其波面与理论波面的运动趋势及波峰波谷值均吻合较好,误差较小,由此说明本文所用数值模型可以很好地模拟波浪运动㊂2.2辐射应力计算结果验证图2 平底水槽中数值模拟辐射应力值及微辐波理论辐射应力结果对比F i g .2 C o m p a r i s o no f n u m e r i c a l s i m u l a t i o na n da i r y w a v e t h e o r y r e s u l t s o f r a d i a t i o n s t r e s s a t f l a t t a n k 将上述工况中数值模拟计算的辐射应力沿程变化结果与微幅波理论下辐射应力结果进行比对,对比结果如图2所示㊂由式12和式13计算得到微辐波理论下的辐射应力值为0.4709N /m ,数值模拟结果中辐射应力值有较小幅度的振荡,可能是由波浪沿水槽传播过程中波幅产生微小差异所致,但整体而言,本模型模拟结果与理论结果对应得较好,说明本数值模型模拟波浪运动输出的流场㊁压力场较准确,辐射应力计算方法也是可行的㊂3 潜堤地形下辐射应力计算及波浪增减水模拟潜堤地形的实验案例采用B e j i 和B a t t je s [19]的实验地形及工况,数值模拟C a s e A 和C a s e B 两种工况下波浪的传播,并计算辐射应力及波浪增减水的沿程变化情况㊂本文采用层流模型来进行数值计算,利用相对较小的网格尺度进行波浪破碎的求解,能够得到较好的波浪模拟结果,且与湍流模型R A N S 相比,层流模型计算所得的破波区内波浪运动与实际波浪运动更为相似,能更好地模拟出波浪的破碎形态和波浪破碎后波面的起伏振荡现象,由于篇幅有限,在此不作更多的说明㊂数值模型布置如图3所示,堤顶高度为0.3m ,波浪沿水槽正向传播㊂数值模拟计算域长40.0m ,前后各设置8.0m 松弛区,计算时长为50s ,C a s e A ,C a s e B 均采用自动调整计算时间步长的方法,输出时间步长分别为Δt =0.02s ,Δt =0.025s ㊂总网格量均为80万,x 方向网格尺寸为0.01m ,y 方向网格尺寸为0.0025m ㊂数值模拟实验工况如表2所示㊂注:潜堤前坡的坡度为1ʒ20,潜堤后坡的波度为1ʒ10;为显示方便,图中水平方向距离与垂直高度不成比例图3 潜堤地形布置(m )F i g .3 L a y o u t o f t h e s u b m e r ge db a r s (m )表2 潜堤地形实验工况列表T a b l e 2 L i s t o f e x p e r i m e n t a l c a s e s o f s u b m e r g e db a r s t o p o g r a p h y 工 况波高H /m 波周期T /s 水深h /m 波浪破碎类型C a s e A 0.0412.5000.4崩破波C a s e B 0.0542.5250.4卷破波162海洋科学进展39卷3.1波形对比当波浪传播到潜堤前坡坡顶位置时,非线性效应达到最大,波浪发生明显的破碎㊂将堤顶G2,G3和G4三处波面的时间变化曲线与实验结果对比,结果见图4,可以看出,在波浪初破阶段,C a s e A和C a s e B工况的G2,G3处数模波峰波谷值与实验结果对应较好;在波浪完全破碎后(G4),2个工况的数值模拟计算结果均出现了低估峰值和谷值的现象㊂但是数值模拟的两工况的整体运动趋势,特别是破碎后波浪振荡趋势与实验现象基本一致,说明该数值模型可以较好地反映实验中波浪破碎后的实际运动,模拟的波浪破碎具有较大的可信度㊂图5给出了潜堤堤顶位置(24~26m处)波浪破碎时的流场信息,可以清楚地看到,当破碎发生时,图4 C a s e A和C a s e B两种工况下G2㊁G3㊁G4处波面对比F i g.4 C o m p a r i s o no fw a v e s u r f a c e a tG2,G3a n dG4u n d e rC a s e Aa n dC a s e B图5 C a s e A工况堤顶位置流场分布F i g.5 T h e f l o wf i e l dd i s t r i b u t i o na t t h e t o p o f t h e d a mu n d e rC a s e A1期穆浩然,等:潜堤上规则波辐射应力的数值研究163 破碎点周围水质点的速度方向均指向破碎点处,越靠近破碎点的水质点速度值越大,在破碎波面处水质点运动最为剧烈,速度值达到最大㊂而距离破碎点较远的流场中,水质点速度较小,水体流动较为温和,但是受前一破碎波浪的影响,水体中水质点移动较为混乱,方向具有不确定性,这与波浪运动的实际情况较为符合,因此本数值方法计算所得流场信息较为可信㊂3.2 辐射应力与波浪增减水当波浪传播到潜堤区域时,会发生浅水变形并最终产生破碎,继而波高骤降,波浪动量减小,辐射应力将沿程发生变化,进而影响平均水平面的改变,这一波浪平均水平面的变化便称为波浪的增减水㊂当规则波沿x 方向正向入射时,依据波浪增减水方程可知:d S x x d x =-ρg (h +η)d ηd x ,(14)当辐射应力(S x x )沿程增长,即d S x x d x >0时,则d ηd x <0,波浪将产生减水现象;当辐射应力沿程降低,即d S x x d x <0时,则d ηd x >0,波浪将产生增水现象㊂波浪增减水的值即为波面升高的时均值与静水面位置的差值㊂本文根据数值模型输出的波面数据,进行周期内时间平均,得到C a s e A 和C a s e B 两个工况下的波浪增减水变化㊂图6和图7给出了潜堤地形下2种工况的波浪辐射应力和波浪增减水的沿程变化曲线㊂2个工况均能明显看出,当波浪传播在堤前平底时,辐射应力基本呈稳定值,到潜堤前坡坡底附近其值开始增大,直至堤顶处达到最大值,随后产生下降的趋势,在潜堤后坡面处由于水深的增加,辐射应力值再次增大,直至堤后平底水槽处趋于平稳㊂而波浪增减水的沿程变化趋势恰好与之相反㊂由于波浪完全破碎,波浪成分不稳定,在堤后坡面至平底处的辐射应力和平均水平面均有明显震荡,但是整体变化趋势依旧能较好地满足波浪增减水方程㊂图6 C a s e A 波浪辐射应力与增减水沿程变化情况F i g .6 V a r i a t i o no fw a v e r a d i a t i o n s t r e s s a n d w a v e s e t -u p ,w a v e s e t -d o w nu n d e rC a s eA 图7 C a s eB 波浪辐射应力与增减水沿程变化情况F i g .7 V a r i a t i o no fw a v e r a d i a t i o n s t r e s s a n d w a v e s e t -u p,w a v e s e t -d o w nu n d e rC a s e B 另外,在堤前平底处,模型稳定后的辐射应力结果较小,这是由于随着波浪的传播,堤前反射波使波幅随着时间推移产生下沉现象,从而导致辐射应力定义公式中的速度积分项和动压积分项变小,最终导致了辐射应力结果偏低㊂由图8可见,C a s e A 工况下15m 水槽处波面过程曲线及波幅偏移值随时间变化的曲线,将波幅偏移值量化为波浪上下幅值之和,并将其进行无量纲化,其中a u 表示波浪的上幅值,a d 表示波浪的下幅值㊂波幅的偏移值在30s 后基本达到稳定状态,同时由图9可见数值水槽前端平底15m 处不同周期内的波幅偏移值与其对应的辐射应力值关系,较明显地反映出潜堤工况下的平底水槽部分辐射应力结果与对164海洋科学进展39卷应周期内波浪幅值的偏移值成正相关,且受波幅变动影响较为敏感,波幅偏移值越小,辐射应力值越低㊂结合图8和图9来看,图8中当波浪波幅偏移值达到稳定状态时,其值对应于图9中辐射应力值大概为0.5 N/m处的点,较理论预估值低,而波浪传播初期的波浪偏移值对应的同周期内的辐射应力结果约为3N/m,这与理论预估值基本一致,由此可见平底水槽处波浪稳定时辐射应力结果偏低是波幅偏移值变小造成的㊂而造成波浪波幅产生振荡的原因还需进一步地分析探讨㊂图8 C a s e A工况中水槽15m处波面变化及波幅偏移量变化F i g.8 W a v e s u r f a c e c h a n g e a n da m p l i t u d e o f f s e t c h a n g e a t15mo n t h e t a n ku n d e rC a s eA图9 C a s e A工况水槽15m处不同周期内波浪波幅偏移量与其对应辐射应力值关系F i g.9 T h e r e l a t i o n f i g u r eb e t w e e nw a v e a m p l i t u d e o f f s e t a n d c o r r e s p o n d i n g v a l u e o f r a d i a t i o n s t r e s sa t15mo f t a n ku n d e rC a s e A4结语本文应用O p e n F O AM数值模型模拟有破碎情况下潜堤地形的波浪传播,详细地输出了碎浪区波浪运动的流场信息,完整地计算了波浪在堤前平底㊁堤身及堤后平底处整个沿程的辐射应力和波浪增减水变化,并探究了堤前平底处辐射应力值受波幅偏移值的影响㊂结果表明,数值模拟波浪在平底水槽中传播的波形及辐射应力结果与解析结果吻合较好;在含潜堤地形的水槽中,得到了与实验结果相吻合的破波区波面对比图,并计算输出了整个水槽内包括碎浪带及堤后平底水槽处的波浪辐射应力变化,发现波浪完全破碎后堤后坡面及平底处辐射应力呈先增大后趋于平稳的趋势,且堤后平底处较堤前平底处辐射应力值更大㊂辐射应力结果与波浪增减水变化趋势在碎浪带内也能较好地对应波浪增减水方程;并且由于波浪反射造成的波幅的上浮与下沉,会间接导致辐射应力结果的增大与减小,且辐射应力值受波幅偏移值影响较为敏感㊂参考文献(R e f e r e n c e s):[1] L O N G U E T-H I G G I N S M S,S T E WA R TR W.C h a n g e s i n t h e f o r mo f s h o r t g r a v i t y w a v e s o n l o n g w a v e s a n d t i d a l c u r r e n t s[J].J o u r n a l o fF l u i d M e c h a n i c s,1960,8(4):565-583.[2] L O N G U E T-H I G G I N S M S,S T E WA R TR W.R a d i a t i o n s t r e s s a n dm a s s t r a n s p o r t i n g r a v i t y w a v e s,w i t h a p p l i c a t i o n t o s u r f b e a t s [J].J o u r n a l o f F l u i d M e c h a n i c s,1962,13(4):481-504.[3] S YMO N D SG,HU N T L E Y D A,B OW E N AJ.T w o-d i m e n s i o n a l s u r fb e a t:L o n g w a v e g e n e r a t i o nb y at i m e-v a r y i n g b r e a k p o i n t[J].J o u r n a l o fG e o p h y s i c a lR e s e a r c hO c e a n s,1982,87(C1):492-498.[4] K O S T E N S EJK.M e a s u r e m e n t s o f s u r f b e a t a n d s e t-d o w nb e n e a t hw a v e g r o u p s[J].C o a s t a l E n g i n e e r i n g,1985,1(19):724-740.1期穆浩然,等:潜堤上规则波辐射应力的数值研究165 [5] C O P E L A N DGJM.P r a c t i c a l r a d i a t i o ns t r e s s c a l c u l a t i o n s c o n n e c t e dw i t h e q u a t i o n s o fw a v e p r o p a g a t i o n[J].C o a s t a l E n g i n e e r i n g,1985,9(3):195-219.[6] Z H E N GJH,Y A NYX.V e r t i c a l v a r i a t i o n s o f w a v e-i n d u c e d r a d i a t i o n s t r e s s t e n s o r[J].C h i n e s e J o u r n a l o fH y d r o d y n a m i c s,2001,16(2):246-253.郑金海,严以新.波浪辐射应力张量的垂向变化[J].水动力学研究与进展(A辑),2001,16(2):246-253.[7] X I A H,X I AZ,Z HUL.V e r t i c a l v a r i a t i o n i n r a d i a t i o n s t r e s s a n dw a v e-i n d u c e d c u r r e n t[J].C o a s t a l E n g i n e e r i n g,2004,51(4):309-321.[8] S T I V E MJ,W I N D H G.As t u d y o f r a d i a t i o ns t r e s s a n d s e t-u p i n t h e n e a r s h o r e r e g i o n[J].C o a s t a l E n g i n e e r i n g,1982,6(1):1-25.[9] W E N X Y,MA XZ,D O N G G H,e t a l.An u m e r i c a lm e t h o d f o r s t u d y i n g t h ew a v e s e t-u p a n ds e t-d o w n[J].C h i n a W a t e rT r a n s p o r t,2009,9(6):142-144.温秀媛,马小舟,董国海,等.一种波浪增减水的数值研究方法[J].中国水运,2009,9(6):142-144. [10] G O S MA N A D.D e v e l o p m e n t s i n i n d u s t r i a l c o m p u t a t i o n a l f l u i d d y n a m i c s[J].C h e m i c a l E n g i n e e r i n g R e s e a r c h a n dD e s i g n,1998,76(2):153-161.[11] Z HAJ J,WA NDC.N u m e r i c a lw a v e g e n e r a t i o n a n d a b s o r p t i o nb a s e d o nO p e n F O AM[J].T h eO c e a nE n g i n e e r i n g,2011,29(3):1-12.查晶晶,万德成.用O p e n F O AM实现数值水池造波和消波[J].海洋工程,2011,29(3):1-12.[12] WA N GDX,S U NJW,G U I J S,e t a l.N u m e r i c a l s i m u l a t i o n o f s o l i t a r y w a v e t r a n s f o r m a t i o n o v e r r e e f p r o f i l e b a s e d o n I n t e r D y M F o a m.王东旭,孙家文,桂劲松,等.基于I n t e r D y M F o a m的潜礁孤立波传播数值模拟[J].海洋工程,2018,36(2):47-55.[13] MA OYJ,MAZ,C H E N GY,e t a l.E f f e c t o f t h e P T Od a m p i n g f o r c e o n t h ew a v e d a m p i n g p e r f o r m a n c e a n dw a v e c a p t u r e e f f i c i e n c y o f aW E C-t y p e b r e a k w a t e r[J].C o a s t a l E n g i n e e r i n g,2019,37(4):45-53.毛艳军,马哲,程勇,等.动力输出系统(P T O)对集成波能转换装置式防波堤消波性能及波能捕获率影响研究[J].海洋工程,2019,37(4):45-53.[14] Y A O Y,H E W R,L IY,e t a l.N u m e r i c a l i n v e s t i g a t i o no f s o l i t a r y w a v e t r a n s f o r m a t i o no v e r r e e f i s l a n d s b a s e do nN a i v e r-S t o k e s e q u a-t i o n s[J].J o u r n a l o fH a r b i nE n g i n e e r i n g U n i v e r s i t y,2018,39(2):392-398.姚宇,何文润,李宇,等.基于N a v i e r-S t o k e s方程珊瑚岛礁附近孤立波传播变形数值模拟[J].哈尔滨工程大学学报,2018,39(2):392-398.[15] Z HA N GCH,Z H E N GQ.S i m u l a t i o n o f g a s-l i q u i d t w o-p h a s e f l o wc o u p l i n g i nw a v e s u r f z o n e i n s h a l l o ww a t e r b a s e d o nO p e n F O AM[J].P o r t&W a t e r w a y E n g i n e e r i n g,2018,548:38-42.张陈浩,郑茜.基于O p e n F O AM的浅滩波浪破碎区气液两相流耦合模拟[J].水运工程,2018,548:38-42.[16]J A C O B S E N N G,F UH R MA N DR,e t a l.A w a v e g e n e r a t i o n t o o l b o x f o r t h e o p e n-s o u r c eC F D l i b r a r y:O p e n F o a m [J].I n t e r n a t i o n a lJ o u r n a l f o rN u m e r i c a lM e t h o d s i nF l u i d s,2012,70(9):1073-1088.[17] W E L L E R H G,An e wa p p r o a c h t oV O F-b a s e d i n t e r f a c e c a p t u r i n g m e t h o d s f o r i n c o m p r e s s i b l e a n d c o m p r e s s i b l e f l o w[E B/O L].(2008-01)[2020-03-22].h t t p s:ʊd o i.o r g/10.1103/P h y s R e v.72.433.[18] L O N G U E T-H I G G I N S M S,S T E WA R TR W.R a d i a t i o n s t r e s s e s i nw a t e rw a v e s:a p h y s i c a l d i s c u s s i o n,w i t h a p p l i c a t i o n s[J].D e e p S e aR e s e a r c ha n dO c e a n o g r a p h i cA b s t r a c t s,1964,11(4):529-562.[19] B E J I S,B A T T J E S JA.E x p e r i m e n t a l i n v e s t i g a t i o no fw a v e p r o p a g a t i o no v e r ab a r[J].C o a s t a l E n g i n e e r i n g,1993,19(1-2):151-162.166海洋科学进展39卷N u m e r i c a l S t u d y o fR a d i a t i o nS t r e s s o n t h e S u b m e r g e dB a r sb y R e g u l a rW a v eMU H a o-r a n1,2,MA X i a o-z h o u1,2,MA O Y a n-j u n1,2,G A O X i a n g1,2(1.F a c u l t y o f I n f r a s t r u c t u r eE n g i n e e r i n g,D a l i a nU n i v e r s i t y o f T e c h n o l o g y,D a l i a n116024,C h i n a;2.T h eS t a t eK e y L a b o r a t o r y o f C o a s t a l a n dO f f s h o r eE n g i n e e r i n g,D a l i a n116024,C h i n a)A b s t r a c t:I n t h e s t u d y o f n e a r s h o r ew a v e,t h e r a d i a t i o n s t r e s s i s a n i m p o r t a n t p h y s i c a l q u a n t i t y w h i c h i s a n e x c e s sm o m e n t u mf l o wc a u s e db y w a v e s.A sw a v e s t r a n s f e r f r o md e e p w a t e r t o s h a l l o w w a t e r,t h e n o n l i n-e a r i t y o fw a v e i n c r e a s e s g r a d u a l l y a n dw a v e s a r ed e f o r m e do r e v e nb r o k e n.T h i sw i l l c a u s e s t r o n g s h i f t i n r a d i a t i o ns t r e s sa n d m a k es i g n i f i c a n t c o n t r i b u t i o n s t ot h e g e n e r a t i o no f i n f r a g r a v i t y w a v e.I nt h i s p a p e r, O p e n F O AMi s u s e d t o f i n e s i m u l a t e t h e p r o p a g a t i o n o fw a v e s o n t h e s u b m e r g e db a r s.I n t h e c a s e o f b r e a k-i n g w a v e s,t h e d e t a i l e d f l o wf i e l d i n f o r m a t i o n o fw a v em o t i o n i s o b t a i n e d.B e s i d e s,t h e r a d i a t i o n s t r e s s a n d t h ew a v e s e t-u p a n ds e t-d o w no nt h e s u b m e r g e db a r s t e r r a i na r e c a l c u l a t e d.N u m e r i c a l r e s u l t s s h o wt h a t t h e r a d i a t i o n s t r e s s v a l u e i s q u i t e s e n s i t i v e t o t h ew a v e a m p l i t u d e s h i f t a t t h e f l a tb o t t o m b e f o r e t h e s u b-m e r g e db a r s.W h e n t h ew a v e i s c o m p l e t e l y b r o k e n,t h e v a l u e o f r a d i a t i o n s t r e s s i n c r e a s e s a t t h e b a c k s u b-m e r g e d-b a r s s l o p e a n d t h e n i t t e n d s t o l e v e l o f f a t t h e f l a t b o t t o mb e h i n d t h e s u b m e r g e db a r s.T h e e f f e c t o f v a r i a t i o nt r e n d o f r a d i a t i o n s t r e s s a l o n g t h ew a y o n t h ew a v e s e t-u p a n d s e t-d o w n i s i n g o o d a g r e e m e n tw i t h t h ew a v e s e t-u p a n d s e t-d o w ne q u a t i o n.K e y w o r d s:w a v e r a d i a t i o n s t r e s s;w a v e s e t-u p a n d s e t-d o w n;s u b m e r g e db a r s;O p e n F O AMR e c e i v e d:M a r c h22,2020。
