基于格子Boltzmann方法的斜坡堤越浪数值模拟研究

准确确定越浪量对于斜坡堤设计有重要意义。利用格子Boltzmann方法(LBM),并采用主动吸收式速度入口造波、出流边界消波、VOF方法追踪自由表面以及静态Smagorinsky模型模拟紊流运动,建立二维数值波浪水槽,对光滑斜坡堤上规则波与不规则波越浪进行数值模拟。模拟结果与试验值及其他数值模型结果比较表明,二维LBM数值波浪水槽具有模拟斜坡堤越浪的能力,但对于破碎较为剧烈的越浪过程模拟,该模型还存在一定的不足,未来可从提高自由表面模型精度等方面进一步改善其性能。     

Numerical study of combined overflow and wave overtopping over a smooth impermeable seawall

A numerical wave flume is used to investigate the discharge characteristics of combined overflow and wave overtopping of impermeable seawalls. The numerical procedure computes solutions to the Reynolds-averaged Navier–Stokes equations and includes the generation of an irregular train of waves, the simulation of wave breaking and interaction with a sloping, impermeable wall. The numerical model is first tested against published experimental observations, approximate analytical solutions and empirical design formulae for the cases of pure overflow and pure overtopping. A sequence of numerical experiments simulating combined overflow and overtopping are described. The results are used to determine empirical discharge formulae of the form used in current practice.     

碟形越浪式波能发电装置的三维数值模拟

基于计算流体力学软件的三维数值模拟技术,分析了碟形越浪式波能发电装置的越浪性能,通过构建基于水气两相VOF(Volume of Fluid)模型的三维数值波浪水槽对该装置进行三维数值模拟研究,数值计算结果与物理试验结果相互对比验证较为吻合,验证了所构建的三维数值波浪水槽的可靠性,通过考察装置的坡度、导流叶片个数、干舷高度对越浪性能的影响确定装置的最优结构参数。结果表明,在装置的斜坡面边缘增加回流板可减少波浪的反射,提高装置的越浪性能。在数值模拟中将装置的斜坡面边缘处安装回流板对碟形越浪式波能发电装置参数进行优化,通过分析回流板的长度对装置越浪性能的影响来探索最优回流板长度。     

孤立波低顶海堤越浪的数值模拟研究

灾害性波浪是中国沿海地区最具破坏性的自然灾害之一。采用开源程序OpenFOAM中interFoam求解器,对低顶海堤（在风暴潮和海平面上升情况下所面临的不利工况）的孤立波越浪特性开展数值模拟研究。通过孤立波冲击海堤的基准算例,验证模型在模拟波浪爬升和越浪过程中大变形波面以及剧烈波浪力方面的精度。基于验证的数值模型,对孤立波在低顶海堤上的越浪特征以及防浪墙高度对越浪的影响开展参数化研究。结果表明堤顶超高减小导致更为剧烈的越浪。针对尚无低顶海堤孤立波越浪量经验公式的问题,提出新的适用于堤顶超高小或为0的孤立波越浪量经验公式。此外,研究发现增加防浪墙高度可有效减少越浪,但防浪墙所受的波浪力也增大。综合考虑防浪墙减少越浪以及自身所受波浪力,针对文中研究采用的海堤截面和波浪条件,建议无量纲防浪墙高度取为1.00。     

不规则波与开孔沉箱相互作用数值模拟

采用三步有限元法对N-S方程进行离散,同时借助CLEAR-VOF方法追踪流体自由表面,利用主动吸收式造波等手段改进了二维不规则波浪数值水槽,使得水槽中的波浪谱与目标靶谱吻合较好。进而建立了不规则波浪与开孔沉箱作用一种新的数值模式,分析研究不规则波作用下开孔沉箱的反射率,并与现有的物模结果和数模结果进行了对比,为不规则波与开孔沉箱作用问题的研究,探求了一种新的数值手段。     

波浪与多孔介质结构相互作用SPH模拟

应用经CSPM法和黎曼解修正后的光滑粒子流体动力学(SPH)方法,建立了主动吸收无反射数值波浪水槽,研究波浪作用下多孔介质结构的水动力特性。流体运动控制方程采用N-S方程,多孔渗水结构内流体的运动控制方程考虑渗流力的影响。数值计算结果给出了水槽内不同位置测点的波面历时曲线和越浪量随时间变化曲线,并同试验结果和Philip Liu的数值计算结果进行了比较。并对一个波浪周期内斜坡堤多孔介质结构内外的速度场和压力分布进行了讨论分析。计算分析表明,数值计算波面较Philip Liu的计算结果与试验结果吻合更好。说明应用SPH方法建立的二维数学模型能够较好地模拟破碎波在多孔渗水斜坡上的爬坡和越浪。     

基于光滑质点流体动力学方法数值波浪水槽研究

基于光滑质点流体动力学(smoothed particle hydrodynamics,简称SPH)方法建立二维数值波浪水槽,利用推波板造波生成孤立波和不规则波,模拟结果跟实验值及解析解做了对比分析,结果令人满意.SPH数值水槽中生成的波浪具有明显的非线性,波浪能量沿程有所衰减,与物理水槽结果吻合良好.     

基于FLUENT的波浪溢流水动力数值模拟

波浪溢流现象使得海堤受到了越浪和溢流的联合作用,复杂的水动力过程会引起海堤后坡产生严重的侵蚀破坏。基于FLUENT软件建立了二维数值波浪溢流水槽模型,该模型运用UDF速度边界造波法分析在不同超高条件下海堤后坡流量和水流厚度的水力学特性。结果表明数值模拟结果与前人物理模型试验结果吻合,该模型可以真实地模拟出海堤波浪溢流现象。在此基础上进一步研究了波浪溢流中越浪和溢流在不同相对超高条件下的主导性作用,而后建立了十分准确的波浪溢流海堤后坡稳定水流厚度计算公式。     

Random wave runup and overtopping a steep sea wall: Shallow-water and Boussinesq modelling with generalised breaking and wall impact algorithms validated against laboratory and field measurements

A semi-implicit shallow-water and Boussinesq model has been developed to account for random wave breaking, impact and overtopping of steep sea walls including recurves. At a given time breaking is said to occur if the wave height to water depth ratio for each individual wave exceeds a critical value of 0.6 and the Boussinesq terms are simply switched off. The example is presented of waves breaking over an offshore reef and then ceasing to break as they propagate inshore into deeper water and finally break as they run up a slope. This is not possible with the conventional criterion of a single onset of breaking based on rate of change of surface elevation which was also found to be less effective generally. The runup distribution on the slope inshore of the reef was well predicted. The model is tested against field data for overtopping available for Anchorsholme, Blackpool and corresponding 1:15 scale wave flume tests. Reflection of breaking waves impacting a steep sea wall is represented as a partial reversal of momentum flux with an empirically defined coefficient. Offshore to nearshore significant wave height variation was reasonably predicted although nearshore model spectra showed distinct differences from the experiments. The breaking wave shape described by a shape parameter was also not well represented as might be expected for such a simple model. Overtopping agreement between model, field and flume was generally good although repeatability of two nominally identical flume experiments was only within 25%. Different distributions of random phase between spectral components can cause overall overtopping rates to differ by up to a factor of two. Predictions of mean discharge by EurOtop methods were within a factor of two of experimental measurements.     

基于物理造波和主动消波的黏性数值波浪水槽开发

基于推板造波理论和摇板造波理论,在Open FOAM平台上采用重叠网格技术建立黏性数值波浪水槽,并使用一种结合SIMPLE算法和PISO算法的PIMPLE算法对数值模型进行求解。利用开发的数值模型通过数值收敛性测试和网格独立性测试分别重点研究了时间步长、库朗数和网格尺寸对数值精度和计算效率的影响。并对比研究了此数值模型分别嵌入层流模型和湍流模型的计算精度和计算效率。实现的规则波和二阶有限振幅波与理论结果和试验结果吻合,验证了此黏性数值波浪水槽的造波和主动消波功能。基于二维数值波浪水槽,进一步研究了三维数值造波,数值计算结果与理论结果吻合良好。研究结果不仅验证了重叠网格在二维和三维两相流体域中求解运动物体与流场交互的可靠性和正确性,而且为使用此黏性数值波浪水槽解决更复杂的海洋工程问题提供了依据。     

Wave run-up and overtopping on smooth and rough slopes of coastal structures

A series of hydraulic model tests has been carried out in a glass wave flume to investigate the influences of wave height, wave period, wave steepness, surf similarity parameter, roughness, layer thickness and porosity on wave run-up and overtopping of 1:2 sloped impermeable and permeable breakwaters fronted by a 1:10 gentle, smooth beach slope. The analysis of results involves the correlation between the overtopping energy transfer with the relative wall height and the relationship between wave run-up and overtopping rate. Further, measured wave run-up and overtopping rates are compared with the results given in the Shore Protection Manual (1984), Automated Coastal Engineering System (1992)and results of other investigators.     

斜向和多向不规则波在斜坡堤上的单波越浪量的实验研究

通过三维物理模型实验对斜坡堤上斜向和多向不规则波的单波越浪量进行了研究.实验考察了入射方向为0°~45°的斜向波和方向分布宽度为0°~25°的多向波以及混凝土和扭工字块体两种护面形式.在混凝土护面堤上用Weibull分布函数拟合了单波越浪量的累积频率分布,在影响因素不同的条件下确定了分布函数中的系数和越浪比例,给出了计算单波越浪量的公式,同时对扭工字块体护面堤上大约100个波中最大的单波越浪量进行了估算.     

Modelling of Random Waves over Submerged Breakwaters and Its Application to Reflection Estimation

1.IntroductionGreat efforts have been devoted to the protection of coastal areas over many years by erectingdikes,seawalls,groin systems,and detached breakwaters.The sea walls,jetties,detached breakwa-ters,etc.are traditionally adopted as absorbingfacilit…     

Numerical analysis of wave overtopping of rubble mound breakwaters

The paper describes the results of a two-dimensional (2-D) numerical modelling investigation of the functionality of rubble mound breakwaters with special attention focused on wave overtopping processes. The model, COBRAS-UC, is a new version of the COBRAS (Cornell Breaking Waves and Structures) based on the Volume Averaged Reynolds Average Navier–Stokes (VARANS) equations and uses a Volume of Fluid Technique (VOF) method to capture the free surface. The nature of the model equations and solving technique provides a means to simulate wave reflection, run-up, wave breaking on the slope, transmission through rubble mounds, overtopping and agitation at the protected side due to the combined effect of wave transmission and overtopping. Also, two-dimensional experimental studies are carried out to investigate the performance of the model. The computations of the free surface and pressure time series and spectra under regular and irregular waves, are compared with the experimental data reaching a very good agreement. The model is also used to reproduce instantaneous and average wave overtopping discharge. Comparisons with existing semi-empirical formulae and experimental data show a very good performance. The present model is expected to become in the near future an excellent tool for practical applications.     

A two-dimensional numerical wave flume based on SA-MPLS method

A spatially adaptive(SA) two-dimensional(2-D) numerical wave flume is presented based on the quadtree mesh system,in which a new multiple particle level set(MPLS) method is proposed to solve the problem of interface tracking,in which common intersection may be traversed by multiple interfaces.By using the adaptive mesh technique and the MPLS method,mesh resolution is updated automatically with time according to flow characteristics in the modeling process with higher resolution around the free surface and the solid boundary and lower resolution in less important area.The model has good performance in saving computer memory and CPU time and is validated by computational examples of small amplitude wave,second-order Stokes wave and cnoidal wave.Computational results also indicate that standing wave and wave overtopping are also reasonably simulated by the model.     

不同掩护程度弧形胸墙波压力及越浪量试验研究

为了明确斜坡堤弧形胸墙越浪量及波压力的变化规律,采用波浪水槽试验测量了弧形胸墙的越浪量和波压力。试验从斜坡堤弧形胸墙前的掩护程度等因素入手对弧形胸墙的返浪效果及波压力进行初步研究,得出不同掩护程度弧形胸墙的越浪量及波压力,发现掩护程度越好,弧形胸墙所受波压力越小;半掩护情况下越浪量最小,为实际工程设计提供了依据。     

Numerical simulation of wave overtopping of coastal structures using the non-linear shallow water equations

A one-dimensional high-resolution finite volume model capable of simulating storm waves propagating in the coastal surf zone and overtopping a sea wall is presented. The model (AMAZON) is based on solving the non-linear shallow water (NLSW) equations. A modern upwind scheme of the Godunov-type using an HLL approximate Riemann solver is described which captures bore waves in both transcritical and supercritical flows. By employing a finite volume formulation, the method can be implemented on an irregular, structured, boundary-fitted computational mesh. The use of the NLSW equations to model wave overtopping is computationally efficient and practically flexible, though the detailed structure of wave breaking is of course ignored. It is shown that wave overtopping at a vertical wall may also be approximately modelled by representing the wall as a steep bed slope. The AMAZON model solutions have been compared with analytical solutions and laboratory data for wave overtopping at sloping and vertical seawalls and good agreement has been found. The model requires more verification tests for irregular waves before its application as a generic design tool.     

Probabilistic Assessment of Wave Overtopping of Nampo Dikein the West Coast of Korea

The Nampo dike, which is located at the west coast of Korea, was destroyed by wave overtopping during the storms on 30 August and 17 September in 1959. In this paper, is performed the probabilistic assessment of wave overtopping of Nampo dike by use of Owen model, Van der Meer ＆ Janssen model and Hedges ＆ Reis model for wave overtopping of seawall. Based on the available tidal and wave data for storm surges in 1989, the risk assessment of wave overtopping of the Nampo dike has been carried out by both Level Ⅱ and Level Ⅲ reliability methods. The calculated resuhs show the general agreement of failure probability between the two methods. By utilizing the rehabilitated cross section of Nampo dike, the failure probability of wave overtopping for the Nampo dike after rehabilitation will be rapidly reduced to that of initial design at crest level of 9.0 m with the improved slope from 1 ： 2 to 1 ： 4 at seaside. Since the sea level may only rise 1.0 m in the next few decades, the failure probability of Nampo dike will be still in the safe range.     

Dynamic response of a monopile wind turbine in waves: Experimental uncertainty analysis for validation of numerical tools

The responses of a monopile offshore wind turbine subjected to irregular wave loads are investigated numerically and experimentally, considering a range of sea states. An extensive experimental campaign was carried out on a fully flexible model, representative of a 5 MW offshore wind turbine, at 1:40 scale. An assessment of the experimental results for the response amplitude operator for regular waves and the 90th percentile seabed bending moment in long-crested irregular waves is carried out using two models (analytical and numerical) for uncertainty propagation, suggesting that bias errors in the model properties and in the wave elevation contribute the most to the total uncertainty. The experimental results are also compared to a numerical model using beam elements and Morison-type wave loads with second order wave kinematics. The numerical model does not capture all of the responses within the level of uncertainty of the experiments, and possible reasons for the discrepancies are discussed.