风作用下孤立波在珊瑚礁上传播变形机制数值模拟研究

基于二维不可压缩两相流模型建立了数值风浪水槽,采用SST k-ω雷诺时均湍流模型,研究了风作用下孤立波在珊瑚礁上的传播变形规律。将计算结果与实验数据对比,证明了该两相流模型计算孤立波在珊瑚礁上传播的准确性,并进一步分析了不同风速对珊瑚礁上孤立波传播变形的影响。结果表明：风的作用会使波面发生随机脉动特征。当地波高随风速的增大而增大;当地波高关于风速的变化梯度随入射波高的增大而增大。风的作用会加快孤立波的传播并且使孤立波提前发生破碎;孤立波开始破碎的位置随风速的增大向远离礁坪的方向移动。反射系数随风速的增大而增大;反射系数关于风速的变化梯度随入射波高的增大而减小;透射系数随风速的增大呈增大趋势。平底区波峰剖面同一水深处的水平流速随风速的增大而增大;且一定的风速不改变水平流速沿水深的变化梯度。有风时波面上方的矢量密度和大小均明显高于无风时且与风速呈正相关,并且波峰上方气流不再循环。随着风速的增大,水气交界面附近的正涡量和负湍流剪应力减小,负涡量和正湍流剪应力增大。水体动能、势能和总能达到高值的时间随风速的增大而减少;水体动能、势能和总能随风速的增大而增大,并且风速对水体动能的相对影响大于势能。     

Numerical calculation of drift current in wind waves

-Drift current induced by wind and waves is investigated with phase-averaged Navier-Stokes equation in which the Reynolds stress is closed by k-ε model. The governing equations are solved by the finite volume method in a system of nonorthogonal coordinates which is fitted to the phase-averaged wave surface. The predicted drift current is fairly reasonable and the drag coefficient of sea-surface predicted with the newly developed interface conditions shows good agreement with previous measurements when breaking waves do not exist.     

风压作用下非线性极端波浪的数值模拟

本文在时域非线性数值波浪水槽中,研究了不同风速条件下极端波浪的特性。采用推板造波的方式生成非线性波浪,基于Jeffrey遮蔽理论将风压项引入自由面动力学边界条件来模拟风压作用,通过高阶边界元法和混合欧拉-拉格朗日时间步进法来求解初边值问题。通过与已发表的聚焦波群实验结果对比验证了该数值模型的准确性,并研究了风压对极端波浪的最大波高、聚焦位置的偏移和波浪谱的演变等波浪性质的影响。本文进一步在数值波浪水槽中引入均匀水流,来模拟风生流对波浪演变的影响。结果表明,风压的存在会少量增大极端波浪的最大波高,波浪的聚焦和解焦过程伴随着明显的能量传递,并且风生流进一步导致了波浪聚焦位置的偏移。     

利用波浪频谱数值模拟畸形波的方法探究

畸形波是一种特殊的波浪,它的波高极高,波峰尖瘦,能量很集中[1].图1是一个比较典型的畸形波[1],这个波列的有效波高是5.65 m,而其中的畸形波波高是有效波高的3.19倍,达到了18.04 m,波峰的高度是13.90 m,波的周期是9.8 s,畸形程度很大.最有名的畸形波是图2给出的1995-01-01-15:20在北海Draupner石油平台记录的波高25.6m的“新年波”,1 200 s长的观测记录显示当时的有效波高是11.92 m,波峰的高度达到18.4 m.由水深70 m,畸形波的特征周期12 s,根据线性色散关系得出波长为220 m.     

Numerical simulation and mechanism analysis of freak waves

A numerical wave model based on the modified fourth-order nonlinear Schroe dinger equation(mNLSE) in deep water was developed to simulate the formation of freak waves and a standard split-step,pseudo-spectral method was used to solve the equation.The validation of the model is firstly verified,then the simulation of freak waves was performed by changing sideband conditions,and the variation of wave energy was also analyzed in the evolution.The results indicate that Benjamin-Feir instability(sideband instability) is an important mechanism for freak wave formation.     

西北太平洋台风浪数值模拟

台风浪的研究对于船舶航行、港口、海洋和近岸建筑物的安全有着重要的现实意义。应用经验公式(MEF)和大风风圈反推(MGCR)两种方法计算最大风速半径,并基于Jelesnianski2型风场模式、藤田气压模式导出的风场模式、藤田气压场和第三代深海波浪模型(WAVEWATCHⅢ)进行了台风Winnie引起的台风浪的后报,将后报值与实测值进行了详细的比较和分析。结果表明:台风期间,MEF模拟出的西北太平洋海域风浪场变化合理,数值模拟结果与实测资料吻合良好,3个波浪测站的有效波高对比误差约为10%,而MGCR模拟结果误差偏大。因此,MEF的模拟结果能够为船舶、海洋与海岸工程的设计提供科学的依据。     

深水波中进动共振的数值模拟研究

进动(precession)共振是一种非线性共振相互作用,2016年才有学者对这一现象进行研究。采用非静压二维自由表面流模型模拟了深水条件下重力波的进动共振现象。通过边界造波的方法产生双色波,分析了触发进动共振的初始条件;探讨了进动共振在小振幅前提条件下发生的简化初始条件。数值模拟分析两组对称测点,对不同测点的波面、能量谱进行对比分析。数值结果表明:非静压二维自由表面流模型可以模拟进动共振现象,并且可以采用双色波作为条件来研究深水五波进动共振现象,进动共振需要一定的能量转化时间,进动共振发生的条件是三波组合的进动频率等于一个系统存在的非线性频率。     

非均匀风影响下风力机三维气动粘性流场的数值模拟

基于面向对象的开源软件OpenFOAM,选择美国国家新能源实验室(NREL)Phase VI风力机为对象,对以往研究较少的非均匀来流风速作用下风力机三维气动粘性流场进行数值模拟。采用较为接近于真实情况的指数型风剖面,计算了轮毂处风速分别为5、10、15和25 m/s四种工况下的叶片表面压力分布、叶片的推力、尾涡等气动力数据,并与均匀来流风速下的风力机气动力学性能进行详细的对比,探讨非均匀风剖面对风力机流场结构和流动特性影响的物理现象和规律。     

Numerical simulation for the two-dimensional nonlinear shallow water waves

This study deals with the general numerical model to simulate the two-dimensional tidal flow, flooding wave (long wave) and shallow water waves (short wave). The foundational model is based on nonlinear Boussinesq equations. Numerical method for modelling the short waves is investigated in detail. The forces, such as Coriolis forces, wind stress, atmosphere and bottom friction, are considered. A two-dimensional implicit difference scheme of Boussinesq equations is proposed. The low-reflection outflow open boundary is suggested. By means of this model,both velocity fields of circulation current in a channel with step expansion and the wave diffraction behind a semi-infinite breakwater are computed, and the results are satisfactory.     

Numerical modeling of wind waves generated by tropical cyclones using moving grids

A version of the WAVEWATCH III wave model featuring a continuously moving spatial grid is presented. The new model option/version is intended for research into wind waves generated by tropical cyclones in deep water away from the coast. The main advantage of such an approach is that the cyclones can be modeled with spatial grids that cover much smaller areas than conventional fixed grids, making model runs with high spatial resolution more economically feasible. The model modifications necessary are fairly trivial. Most complications occur due to the Garden Sprinkler effect (GSE) and methods used to mitigate it. The basic testing of the model is performed using idealized wind fields consisting of a Rankine vortex. The model is also applied to hurricane Lili in the Gulf of Mexico in October 2002. The latter application shows that the moving grid approach provides a natural way to deal with hurricane wind fields that have a high-resolution in space, but a low resolution in time. Although the new model version is originally intended for tropical cyclones, it is suitable for high-resolution modeling of waves due to any moving weather pattern.     

近岸波浪破碎区不规则波浪的数值模拟

基于近岸不规则波浪传播的抛物型缓坡方程和两类波浪破碎能量损耗因子,对近岸波浪破碎区不规则波浪的波高分布进行了数值模拟,并结合实验结果对数值模拟结果进行了验证分析,结果表明采用两类波浪破碎能量损耗因子所模拟的破碎区波高与实测值均吻合良好,波浪破碎能量损耗因子及波浪破碎指标对破碎区波浪波高分布影响较明显。     

基于OpenFOAM三维多向畸形波数值模拟研究

基于二次开发开源计算流体力学(CFD)软件包Open FOAM中的非稳态不可压缩两相流求解器inter Foam建立数值模型,增加了基于方向谱的三维造波边界和吸收边界,成功模拟了三维多向畸形波过程。通过与模型试验波面和目标谱对比,验证了数值模型模拟三维多向畸形波的有效性。另外,分析了网格尺度和柯朗数对模拟结果的影响,并使用连长统计和SIWEH两种方法分析了包含畸形波波列的群性。研究结果表明:在本研究范围内,网格尺度设置0.02 m×0.01 m×0.02 m,最大柯朗数选择0.25,模拟出的三维畸形波效果最好;从能量角度描述畸形波的群性更为合理。     

Numerical modeling of wind and waves for Typhoon Betty (8710)

Ｎｕｍｅｒｉｃａｌ　ｍｏｄｅｌｉｎｇ　ｏｆ　ｗｉｎｄ　ａｎｄ　ｗａｖｅｓ　ｆｏｒ　ＴｙｐｈｏｏｎＢｅｔｔｙ（８７１０）TXＮｕｍｅｒｉｃａｌｍｏｄｅｌｉｎｇｏｆｗｉｎｄａｎｄｗａｖｅｓｆｏｒＴｙｐｈｏｏｎＢｅｔｔｙ（８７１０）ＹｕＷｅｉｄｏｎｇ，Ｑｉａ...     

Numerical simulation of breaking waves by a multi-scale turbulence model

In this paper, a two-dimensional multi-scale turbulence model is proposed to study breaking waves. The purpose of developing this model is to produce a relatively accurate model with moderate computer requirements. The free surface is tracked by the VOF technique, the log-law profile for the mean velocity is applied at the bottom. Comparing with the Reynolds-Averaged Navier-Stokes models (RANS), the present model shows improving agreement with experimental measurements in terms of surface elevations, particle velocities, wave height distributions and undertow profiles. The subgrid scale (SGS) turbulent transport mechanism is also discussed in the paper. It is found that turbulent production and dissipation are of the same order, but turbulent production is primarily located at the wavefront and above the wave trough, whereas turbulent dissipation is primarily located at the back face of a wave, indicating that in these regions, the assumption of equilibrium is not correct. Below the trough level, the local equilibrium assumption is reasonable. Turbulent convection and diffusion are of the same order at the trough level. Above the trough level, turbulent convection dominates. Under the spilling breaking wave, turbulent kinetic energy is continue to dissipate in the bore region, whereas under the plunging breaking wave, the turbulent kinetic energy is dissipated very rapidly within one wave period.     

Numerical simulation of tsunami waves generated by deformable submarine landslides

This paper presents a new submarine landslide model based on the non-hydrostatic wave model NHWAVE of Ma et al. (2012). The landslide is modeled as a water–sediment mixture. The dense plume is driven by baroclinic pressure forcing introduced by spatial density variations. The model is validated using laboratory measurements of turbidity currents and of water wave generation by a granular landslide. The model is then utilized to study the dependence of landslide motion and associated tsunami wave generation on parameters including sediment settling velocity, initial depth of the landslide and slide density. Model results show that the slide motion and water waves which it generates are both sensitive to these parameters. The relative tsunamigenic response to rigid and deformable landslides of equal initial geometry and density is also examined. It is found that the wave energy is mostly concentrated on a narrow band of the dominant slide direction for the waves generated by rigid landslides, while directional spreading is more significant for waves generated by deformable landslides. The deformable landslide has larger speed and acceleration at the early stage of landslide, resulting in larger surface waves. The numerical results indicate that the model is capable of reasonably simulating tsunami wave generation by submarine landslides.     

浮式风机气动—水动—气弹性耦合响应数值模拟

随着海上风电产业的快速发展,大型浮式风机逐渐从概念设计走向工程应用,但仍面临较大的挑战。一方面,在风、浪等环境载荷的作用下,浮式风机的气动载荷和水动力响应之间存在明显的相互干扰作用;另一方面,风力机大型化使得叶片细、长、薄的特点愈发突出,叶片柔性变形十分显著,这会影响到浮式风机的耦合性能。基于两相流CFD求解器naoe-FOAM-SJTU,结合弹性致动线模型和等效梁理论,建立了浮式风机气动—水动—气弹性耦合响应计算模型,并对规则波和剪切风作用下Spar型浮式风机的气动—水动—气弹性耦合响应进行了数值模拟分析。结果表明,风力机气动载荷使得叶片挥舞变形十分显著,而叶片的扭转变形会明显降低风力机的气动载荷。此外,风力机气动载荷会增大浮式平台的纵荡位移和纵摇角,同时,浮式平台运动响应会导致风力机气动载荷产生大幅度周期性变化。进一步地,叶片结构变形响应会使得浮式风机尾流场的速度损失和湍动能有所降低。     

内孤立波与深海立管相互作用数值模拟

提出了一种对内孤立波与深海立管相互作用耦合数值模拟方法。流场采用内孤立波数值水槽方法进行模拟,结构响应采用基于薄壳理论的有限元方法进行计算,采用一种将流场和结构响应数据进行实时传输的方法,实现了流体与固体之间的耦合数值模拟。对内孤立波作用下某长径比为1 200的深海立管载荷及其动力响应特性进行了数值模拟与分析。结果表明内孤立波不仅会对深海立管产生突发性剪切载荷作用,而且还会使立管产生大幅度变形响应现象,因此在深海立管设计与应用中,内孤立波的影响是不可忽视的。研究表明,该方法为研究内孤立波作用下深海立管动力特性及其工程预报相关问题提供了一种有效的手段。     

利用被动示踪物模拟对黑潮入侵南海的数值研究

由于缺少观测数据和对黑潮水准确定义,很难识别出从太平洋入侵到南海的黑潮水团。本文基于一个经过观测验证的三维模式MITgcm,利用被动示踪物标记黑潮水,研究了入侵南海的黑潮水的时空变化。研究表明,在冬季,黑潮水入侵的范围最广,几乎占据了18°N-23°N和114°E-121°E的区域;并有一个分支进入台湾海峡;黑潮入侵的范围随深度增加逐渐减小。在夏季,黑潮水被限制在118°E以东,且没有分支进入台湾海峡;入侵的范围从海面到约205米是增大的,之后随深度增加逐渐减小。通过分析从2003年到2012年黑潮入侵的年际变化,与厄尔尼诺年和正常年相比,冬季黑潮入侵后向台湾海峡的分支在拉尼娜年是最弱的,这可能与中国大陆东南方向的风应力旋度有关。通过吕宋海峡的黑潮入侵通量（KIT）是西向的,其年平均值约为-3.86×10⁶ m³/s,大于吕宋海峡通量（LST,约-3.15×10⁶ m³/s）。250米以上的KIT约占了全深度通量的60-80%。此外,从2003年到2012年KIT与Niño 3.4指数的相关系数到达0.41,小于LST与Niño 3.4指数的相关系数0.78。     

椭圆余弦内波及其对墩柱作用数值模拟

基于FLUENT软件建立数值波浪水槽,研究椭圆余弦内波及其对墩柱的作用。椭圆余弦内波的生成采用推板造波方法,自由面捕获采用VOF方法。首先,模拟了椭圆余弦内波的生成,比较了不同周期和波高下椭圆余弦内波的波浪形态。然后,在特定的周期和波高条件下模拟了三维椭圆余弦内波对单个墩柱及多个敦柱的作用;分析了各墩柱上受到的惯性力和粘性力的变化趋势,并对不同墩柱下的总波浪力进行了对比;得出了墩柱上受到的波浪力也呈现周期性变化并且在波浪前进方向上后墩柱的受力小于靠前墩柱的结论。