珊瑚岛礁孤立波爬坡的平面二维数值模拟研究

为了探究岛屿周围珊瑚礁在抵御海啸灾害中的作用,采用激波捕捉类Boussinesq模型FUNWAVE-TVD,对孤立波在理想化三维岛礁地形上的传播及爬坡开展了现场尺度的平面二维数值模拟,分析了入射波高、礁坪水深、礁坪宽度、礁前斜坡坡度、礁后斜坡坡度、珊瑚礁糙率对岛屿四周孤立波爬高分布的影响。结果表明,珊瑚礁的存在总体上可有效降低岛屿四周孤立波的最大爬坡高度;入射波高、礁坪水深、礁坪宽度、珊瑚礁糙率是影响珊瑚岛礁四周孤立波爬坡分布的主要因素,岛礁四周最大爬坡高度会随入射波高和礁坪水深的增大、礁坪宽度和珊瑚礁糙率的减小而不断增大;当礁坪水深增大到一定程度时,珊瑚礁主要会对岛屿背浪面的爬高失去影响,而当礁坪宽度和珊瑚礁糙率减小至一定程度时,会出现岛礁四周最大爬高高于无珊瑚礁时爬高的现象;礁后斜坡的变缓会使岛礁周围的最大爬高有所减小,而礁前斜坡坡度对珊瑚岛礁周围的最大爬高几乎没有影响。     

孤立波对礁坪上直墙冲击试验和RANS数值模拟

在有、无礁冠两种典型岛礁剖面地形上开展孤立波冲击礁坪上直墙的二维水槽试验和数值模拟研究。通过OpenFOAM结合k-ε湍流模型进行RANS数值模拟,基于试验数据和验证后的RANS模拟结果,研究了岛礁剖面上孤立波传播的水动力特性及其对直墙的作用,并探讨了礁冠的存在对直墙所受冲击的影响。结果表明:RANS(k-ε)模型能够准确还原孤立波在岛礁地形上的主要传播特征,但在波浪破碎的情况下,RANS数值模拟不够精细,尤其对直墙所受瞬击荷载的捕捉能力欠缺;礁冠的存在,增强了波浪反射,使得波浪提前破碎,减小了直墙上的波浪爬升和最大动压。     

礁面大糙率存在下孤立波传播变形及爬高实验研究

本文采用圆柱体阵列来模拟珊瑚礁面的大糙率,通过波浪水槽实验研究礁面糙率对孤立波传播变形及岸滩爬高的影响。结果表明,粗糙礁面的存在显著削弱了礁坪上孤立波的首峰和礁后岸滩反射造成的次峰,同时降低了波浪在珊瑚礁面的传播速度;垂直于岸线方向沿礁相对波高随着入射波增大而减小,随着礁坪水深的增大而增大,粗糙礁面上波高沿礁的衰减更为显著;礁前斜坡的无量纲反射波高随无量纲入射波高的变化与礁坪水深有关,当入射波高足够大时其趋于常值,粗糙礁面略微增大了礁前斜坡的反射;无量纲透射波高和岸滩爬高随着无量纲入射波高的增大而减小,特别是礁坪水深较大时更为显著,粗糙礁面时的无量纲岸滩爬高相对于光滑礁面平均减小46%;通过回归分析得出了同时适合于光滑和粗糙礁面的预测孤立波岸滩爬高的经验关系式。     

内孤立波过陆架陆坡地形的数值模拟研究

内波是海洋中普遍存在的波动形式。内孤立波是典型的非线性内波,多发于陆架边缘海,如南海等海域,对陆架海域有重要影响。本文针对内孤立波在陆架地形上的传播问题,先基于弱非线性与全非线性数值模型,模拟了不同振幅、地形高度条件下内孤立波的演化的过程,探讨了动力系数对内孤立波演化过程的影响,对比了两模型的模拟结果在内孤立波演化过程、能量分配以及能量耗散的差异,后分析了南海的动力系数分布特征。结果表明,在内孤立波不发生破碎的情况下,弱非线性模型与全非线性模拟结果相近。当发生破碎过程时,弱非线性模型可准确模拟头波,但无法通过强非线性的破碎过程耗散能量,只能以裂变的方式辐射能量。在弱非线性模型中,随地形高度增加,频散系数减小到零,平方非线性系数由负转正,立方非线性系数绝对值增大一个量级,并主导陆架地形上内孤立波的演化过程。通过对比南海夏季与冬季非线性内波动力系数空间分布,发现内孤立波在传播过程由于夏季平方非线性效应、立方非线性效应与频散效应较强的影响,其在夏季更易发生陡化与裂变,波列发生频率高。     

基于InterDyMFoam的潜礁孤立波传播数值模拟

建立了基于OpenFOAM动边界的类fixedValuePointPatchVectorField继承的仿物理造波数值波浪水槽。对孤立波在三种不同潜礁地形上的行进、爬坡以及破碎等典型过程进行数值模拟,模拟结果与Boussinesq模型及物理模型试验所得结果进行了对比分析。结果表明,采用动边界进行仿物理造波更适合处理波高水深比较大的孤立波传播问题,可以较好地模拟孤立波在潜礁上传播引起的波浪破碎、水跃等现象。     

南沙群岛礁概化地形上拍岸浪数值模拟与统计分析

南沙群岛珊瑚岛礁众多,大多数岛礁具有向海坡陡峭、外礁坪比较平缓的特征。将南沙群岛岛礁的迎浪向地形概化为陡坡和缓坡组成的双斜坡,采用FUNWAVE-TVD模式数值模拟概化地形上的波浪,根据模拟的破碎波高分析其拍岸浪特征。对拍岸浪数值模拟结果进行比较分析,向海坡的坡度对拍岸浪影响不大,外礁坪上拍岸浪高随地形坡度增大而略有增大;向海坡和外礁坪交界位置(即坡折点)水深对拍岸浪有比较明显的影响,拍岸浪高随坡折点水深增大而减小;拍岸浪高随入射波高和波周期增大而增大。利用大量的拍岸浪数值模拟数据对国内外5种统计模型进行检验,并且基于拍岸浪数值模拟数据建立了3种南沙群岛岛礁拍岸浪统计模型,计算结果显示这些模型适用性较好。     

基于类海啸波型的岛礁水动力特性数值模拟研究

基于非静压单相流模型NHWAVE建立了高精度二维数值波浪水槽,采用日本2011年实测真实海啸波型系统研究了海啸波在岛礁上传播变形的规律,并且分析了波高、礁坪淹没水深和礁前斜坡坡度等因素对孤立波和真实海啸传播变形的影响。结果表明,相比孤立波,类海啸波的波长明显大于孤立波波长,在测点处引起的水面变化持续时间更长,同等波高情况下真实海啸波型比孤立波能够携带更多的能量,与岛礁的相互作用也更为复杂,在礁坪上形成的淹没水深约为孤立波的两倍。礁前斜坡坡度和礁坪淹没水深均对类海啸波的反射和透射系数有显著影响。随着礁前斜坡坡度的增加,反射系数和透射系数均逐渐增加。随着礁坪淹没水深的增加,反射系数逐渐减小,而透射系数逐渐增大。但是,反射系数和透射系数均随着入射波高的增加而逐渐减小。     

基于激波捕捉类Boussinesq模型对岛礁地形上波浪传播的数值模拟

为了探究激波捕捉类Boussinesq模型在模拟岛礁地形上规则波和不规则波传播的可行性,采用基于完全非线性Boussinesq方程并具有激波捕捉能力的数值模型Funwave-TVD对规则波和不规则波在岛礁地形上的传播进行了数值模拟,通过与试验数据对比,分析模型中空间步长的影响,验证模型在模拟波高、平均水位分布以及波谱空间演变的能力,结果表明:采用合适的空间步长,模型能较好地模拟规则波和不规则波在岛礁地形上的传播和演化过程。对于规则波,较小的空间步长可改善破碎点处波高峰值的预测,并能更好地预测波浪破碎后波高的空间分布,相比结合经验破碎的Boussinesq模型,Funwave-TVD能更好地模拟规则波在岛礁地形上的破碎,以及破碎以后行进涌波的再生成过程;对于不规则波,Funwave-TVD总体而言能较好地模拟涌浪有效波高、次重力波的生成及空间演化和平均水位,但会低估礁坪上次重力波波高,较粗的空间步长也会低估礁坪上涌浪有效波高。     

内孤立波过山脊地形演化实验研究

通过模型实验,研究了下沉型内孤立波通过山脊地形演化特征。实验以三角形障碍物模拟海底山脊地形,采用两种密度的分层水,对上层流体和下层流体的高度比、密度等进行了量化处理。实验研究表明:KdV理论波形可较好模拟本次实验内孤立波波形,但随着内孤立波振幅的增大,误差增加;在内孤立波与障碍物微量作用、中度作用和破波作用三种程度的相互作用中,内孤立波过障碍物具有不同的波形变化和主波能量衰减率。     

孤立波爬坡的二维数值模拟

波浪在滩地上以及遇海岸工程后传播发生变形、爬坡等现象,对其进行数值模拟具有广阔的工程应用背景.应用基于Boltzmann方程的KFVS(kinetic flux vector splitting)格式求解二维浅水方程,同时采用干底Riemann解模拟动边界问题.模型模拟了孤立波在滩地上传播变形、爬坡的过程,以及孤立波在滩地上遇圆柱后绕射、变形和爬坡的过程.计算结果与实验结果非常吻合,表明模型具有较大的推广应用价值.     

Study of irregular wave run-up over fringing reefs based on a shock-capturing Boussinesq model

This paper presents the results of a parametric study of irregular wave run-up over fringing reefs using the shock-capturing Boussinesq wave model Funwave-TVD to better understand the role of fringing reefs in the mitigation of wave-driven flooding. Laboratory experiments were newly performed with a typical fringing reef profile and typical hydrodynamic conditions to validate the model. Experimental data shows irregular wave run-ups are dominated by the low-frequency motions and confirms the run-up resonant phenomenon over the back-reef slope, which has been revealed in previous numerical studies. It is demonstrated that irregular wave evolution and run-up over fringing reefs are reasonably reproduced by the present model with a proper grid size. However, the infragravity run-up height and highest 2% run-up height over the back-reef slope are under-predicted due to the underestimation of the infragravity wave height over the reef flat. The validated model was then utilized to model irregular wave transformations and run-ups under different conditions. Through a series of numerical experiments, the effects of key hydrodynamic and reef geometry parameters, including the reef flat width, water depth over the reef flat, fore-reef slope angle and back-reef slope angle, on the irregular wave run-up were investigated. Variations of spectral components of irregular wave run-ups were examined to better understand the physical process underlying the effect of each parameter.     

珊瑚岸礁破碎带附近波浪演化实验研究

通过波浪水槽实验对珊瑚岸礁破碎带附近波浪演变规律开展研究,实验采用了概化的岸礁模型,测试了4种礁坪水深、4种礁前斜坡坡度和一系列入射波高的组合工况。对破碎带宽度和破碎带附近波浪的入射、反射、透射以及能量耗散进行了测量分析,透射波的计算考虑了礁坪上高次谐波的影响。结果表明:礁坪水深和入射深水波高的比值(即礁坪相对水深)是影响岸礁破碎带附近波浪演化的关键参数,而礁前斜坡坡度的影响在本文测量的范围内可以忽略不计。破碎带宽度与礁坪上浅水波波长为同一数量级,并与礁坪相对水深成反比;透射系数随礁坪相对水深的增大呈线性增长,而反射系数的变化却无类似规律;岸礁能够削弱超过50%入射波能,礁坪相对水深越小,波浪破碎造成的能量耗散越大。     

珊瑚礁地形上破碎波高试验研究

破碎波高是珊瑚礁地形上波浪演化的重要参数之一,对工程安全和海岸变形具有重要影响。通过二维波浪水槽,对珊瑚礁地形上破碎波高进行试验研究,分析破碎波高随波陡、礁坪水深以及礁前斜坡坡度的变化。研究表明,相对破碎波高随相对礁坪水深的增大而增大,随入射波陡的增大而减小,但礁前斜坡坡度对相对破碎波高的影响并不明显。通过引入相对礁坪水深,将经典的破碎波高计算公式拓展至珊瑚礁地形上破碎波高的计算。该公式计算值与前人的试验值进行对比验证,吻合较好。研究成果可为工程实践和数值模拟提供参考与借鉴。     

Numerical Study of the Impact of Climate Change on Irregular Wave Run-up Over Reef-Fringed Coasts

Wave hydrodynamics over fringing reefs is largely controlled by the reef surface roughness and hydrodynamic forcing. It is believed that climate change will result in a net increase in the water depth over the reef flat, a degrading of the surface roughness of coral reefs and changes in extreme incident wave heights. For an accurate assessment of how climate change affects the safety of reef-fringed coasts, a numerical study of the impact of climate change on irregular wave run-up over reef-fringed coasts was carried out based on a Boussinesq wave model,FUNWAVE-TVD. Validated with experimental data, the present model shows reasonable prediction of irregular wave evolution and run-up height over fringing reefs. Numerical experiments were then implemented based on the anticipated effects of climate change and carried out to investigate the effects of sea level rise, degrading of the reef surface roughness and increase of extreme incident wave height on the irregular wave run-up height over the backreef beach respectively. Variations of run-up components(i.e., spectral characteristics of run-up and mean water level) were examined specifically and discussed to better understand the influencing mechanism of each climate change-related effect on the run-up.     

Study on Influences of Fringing Reef on Harbor Oscillations Triggered by N-Waves

Influences of topographic variations of the offshore fringing reef on the harbor oscillations excited by incident Nwaves with different amplitudes and waveform types are studied for the first time. Both the propagation of the Nwaves over the reef and the subsequently-induced harbor oscillations are simulated by a Boussinesq-type numerical model, FUNWAVE-TVD. The present study concentrates on revealing the influences of the plane reef-face slope,the reef-face profile shape and the lagoon width on the maximum runup, the wave energy distribution and the total wave energy within the harbor. It shows that both the wave energy distribution uniformity and the total wave energy gradually increase with decreasing reef-face slope. The profile shape of the reef face suffering leading-elevation Nwaves(LEN waves) has a negligible impact on the wave energy distribution uniformity, while for leading-depression N-waves(LDN waves), the latter gradually decreases with the mean water depth over the reef face. The total wave energy always first increases and then decreases with the mean water depth over the reef face. In general, the total wave energy first sharply decreases and then slightly increases with the lagoon width, regardless of the reef-face width and the incident waveform type. The maximum runup subjected to the LEN waves decreases monotonously with the lagoon width. However, for the LDN waves, its changing trend with the lagoon width relies on the incident wave amplitude.     

Effects of reef morphology variations on wave processes over fringing reefs

Many low-lying tropical and sub-tropical atolls fringed by coral reefs are susceptible to coastal inundation during extreme wave events. Previous studies have shown that the infragravity (IG) wave is the dominant component of shoreline run-up compared to the sea and swell (SS) wave and the wave-induced setup. To better understand both the SS and IG wave dynamics over a fringing reef with various morphologies, a series of laboratory experiments were conducted in a wave flume based on an idealized fringing reef profile. The shoreline responses of waves to different reef morphologies with/without the reef crest, the lagoon and the reef surface roughness were examined. IG wave resonance on the reef flat was identified by a spectral analysis of the shoreline wave records. Subsequently, a numerical model based on the Boussinesq equations was validated by the experimental data. The model was then applied to investigate the impacts of varying reef morphologic features (fore-reef slope, reef-crest width, lagoon width, and reef roughness coefficient) on the shoreline wave motions.