Numerical simulation of nonlinear dispersive waves propagating over a submerged bar by IB–VOF model

It is a good test for a numerical model to simulate progressive waves propagating over a submerged bar with a relatively high ratio of slopes. In this paper, the combined IB–VOF model is used to predict nonlinear dispersive waves propagating over a submerged bar with both slopes of 1:2. The predicted free surface elevations are compared with the experimental data and numerical results presented by other researchers. The comparison shows that the IB–VOF model is able to provide satisfactory wave profiles in the shallow water with strong nonlinear effects and in the wave transmitted region with strong wave dispersion in particular. Moreover, the wave evolution behind the submerged bar is described in detail, including the spatial wave profile modulation, spectral analysis of the time-series waves, flow velocity and pressure fields, and kinetic energy distribution. The effect of fluid viscosity on the numerical simulations is also studied, and it is found that the effect on the wave evolution considered in this paper is not significant. Finally, the hydrodynamic force acting on the bar is calculated using the IB–VOF model.     

Simulation of wave propagation over a submerged bar using the VOF method with a two-equation k– turbulence modeling

A two-equation k– turbulence model is used in this paper to simulate the propagation of cnoidal waves over a submerged bar, where the free surface is handled by the volume-of-fluid (VOF) method. Using a VOF partial-cell variable and a donor–acceptor method, the model is capable of treating irregular boundaries, including arbitrary bottom topography and internal obstacles, where the no-slip condition is satisfied. The model also allows the viscous sublayer to be modeled by a wall function approximation implemented in the grid nodes that are immediately adjacent to a wall boundary. The numerical model applied to the propagation of cnoidal waves over a submerged bar can produce results that are in general agreement with some laboratory measurements. Some remarks arising from the comparison between the computational and experimental results are presented.     

潜堤后高阶自由谐波的研究

基于高阶边界元方法的完全非线性数值水槽模型模拟潜堤地形上波浪的传播变形,通过与实验值进行比较,考察数学模型的正确性.采用两点法分离得到堤后高倍频自由波来研究入射波参数、水深对堤后高倍频自由波的影响.研究发现:基频波、二阶和三阶自由波幅值分别与入射波波幅成线性、二次和三次函数关系,基频波幅值基本不随波浪周期变化,而二阶和...     

Modelling of run up of steep non-breaking waves

A Navier–Stokes solver is used to examine steep waves as they run up a steep beach (10.54°). The volume of fluid method (VOF) is used to model the free surface. Comparison with experimental results shows reasonable overall agreement in the prediction of the free-surface, velocities and accelerations within the flow. A spurious feature at the free-surface was found which does reduce the quality of the results. For a steep wave we see the transition from a steep wave front to a smooth run-up tongue at the beach that is in qualitative agreement with experiment.     

An active wave generating–absorbing boundary condition for VOF type numerical model

The objective of the present work is to discuss the implementation of an active wave generating–absorbing boundary condition for a numerical model based on the Volume Of Fluid (VOF) method for tracking free surfaces. First an overview of the development of VOF type models with special emphasis in the field of coastal engineering is given. A new type of numerical boundary condition for combined wave generation and absorption in the numerical model VOFbreak² is presented. The numerical boundary condition is based on an active wave absorption system that was first developed in the context of physical wave flume experiments, using a wave paddle. The method applies to regular and irregular waves. Velocities are measured at one location inside the computational domain. The reflected wave train is separated from the incident wave field in front of a structure by means of digital filtering and subsequent superposition of the measured velocity signals. The incident wave signal is corrected, so that the reflected wave is effectively absorbed at the boundary. The digital filters are derived theoretically and their practical design is discussed. The practical use of this numerical boundary condition is compared to the use of the absorption system in a physical wave flume. The effectiveness of the active wave generating–absorbing boundary condition finally is proved using analytical tests and numerical simulations with VOFbreak².     

A simple two-way coupling method of BEM and VOF model for random wave calculations

A numerical method, which combines the boundary element method (BEM) and the volume of the fluid method (VOF method), has been presented to solve wave–structure interactions; the intense wave motion at the proximity of the structure is modeled by the VOF method and the rest of the fluid region is modeled by the BEM. The combined method can considerably reduce the time-consuming VOF domain, and thus practically makes it possible to apply the VOF method for random wave calculations, in which long time computations are usually required to obtain statistically meaningful results, and therefore the use of the single-VOF model often becomes prohibitive in terms of computational time and storage memories. A VOF model CADMAS-SURF, which is based on SMAC scheme and had been constructed by a number of VOF researchers in coastal engineering in Japan, is used in the combined BEM–VOF model. The two-way coupling treatment, which enables us to deal with bidirectional wave propagations, which was originally given for the SOLA-VOF model by Yan et al. (2003a) and later improved by Kim et al. (2007), was modified for the SMAC scheme. The coupling treatments are described in detail in the paper. The validity of the combined BEM–VOF model was investigated by comparing the numerical results with the theoretical results for the propagations of Stokes 5th order waves and random waves.     

A finite difference method for effective treatment of mild-slope wave equation subject to non-reflecting boundary conditions

A numerical model for coastal water wave motion that includes an effective method for treatment of non-reflecting boundaries is presented. The second-order one-way wave equation to approximate the non-reflecting boundary condition is found to be excellent and it ensures a very low level of reflection for waves approaching the boundary at a fairly wide range of the incidence angle. If the Newman approximation is adopted, the resulting boundary condition has a unique property to allow the free propagation of wave components along the boundary. The study is also based on a newly derived mild-slope wave equation system that can be easily made compatible to the one-way wave equation. The equation system is theoretically more accurate than the previous equations in terms of the mild-slope assumption. The finite difference method defined on a staggered grid is employed to solve the basic equations and to implement the non-reflecting boundary condition. For verification, the numerical model is then applied to three coastal water wave problems including the classical problem of plane wave diffraction by a vertical circular cylinder, the problem of combined wave diffraction and refraction over a submerged hump in the open sea, and the wave deformation around a detached breakwater. In all cases, the numerical results are demonstrated to agree very well with the relevant analytical solutions or with experimental data. It is thus concluded that the numerical model proposed in this study is effective and advantageous.     

Higher Harmonics Induced by Waves Propagating over A Submerged Obstacle in the Presence of Uniform Current

To investigate higher harmonics induced by a submerged obstacle in the presence of uniform current, a 2D fully nonlinear numerical wave flume（NWF） is developed by use of a time-domain higher-order boundary element method（HOBEM） based on potential flow theory. A four-point method is developed to decompose higher bound and free harmonic waves propagating upstream and downstream around the obstacle. The model predictions are in good agreement with the experimental data for free harmonics induced by a submerged horizontal cylinder in the absence of currents. This serves as a benchmark to reveal the current effects on higher harmonic waves. The peak value of non-dimensional second free harmonic amplitude is shifted upstream for the opposing current relative to that for zero current with the variation of current-free incident wave amplitude, and it is vice versa for the following current. The second-order analysis shows a resonant behavior which is related to the ratio of the cylinder diameter to the second bound mode wavelength over the cylinder. The second-order resonant position slightly downshifted for the opposing current and upshifted for the following current.     

基于高阶边界元的三维数值波浪港池--波浪破碎的模拟

在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性.     

复杂结构形式的海堤波浪力及波浪形态数值模拟

依据雷诺方程和k-ε紊流模型,按流体体积(VOF)法追踪波浪自由表面,采用源造波法,建立数值波浪水槽,数值模拟波浪对复杂结构形式海堤的作用.数值模拟结果与经验公式、物理模型试验结果基本符合,说明所建立的数值波浪水槽合理可行.揭示了不规则波作用下复杂结构形式海堤波浪力分布规律,模拟了堤前波浪形态变化,为探讨合理的海堤结构形式提供了依据.     

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied.     

多个淹没水平圆柱诱发高阶谐波特性的数值与试验研究

基于势流理论采用时域高阶边界元方法建立了模拟非线性波浪与淹没水平双、三圆柱作用的数值水槽模型,其中采用混合欧拉-拉格朗日方法更新自由水面,四阶Runga-Kutta方法进行时间步进。利用两点法分离得到潜体下游高阶自由波,进而研究潜体间距、布置形式等对高倍频自由波的影响。同时在物理水槽内开展多潜体高阶谐波物理模型试验,并将试验结果与数值结果进行对比,吻合较好。研究发现:淹没双圆柱下游高阶谐波幅值随潜体间距呈现周期性振荡变化,其振荡的重现距离约为波长一半。而淹没三圆柱下游的高阶谐波随第一个间距呈周期性振荡变化,随第二个间距的增大而减小。     

极限波浪运动特性的非线性数值模拟

利用时域高阶边界元方法建立了模拟极限波浪运动的完全非线性数值模型,其中自由水面满足完全非线性自由水面条件.采用半混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除水槽两个侧面和底面上的积分.研究中利用波浪聚焦的方法产生极限波浪,并且在水槽中开展了物理模型实验,将测点试验数据与数值结果进行了对比,两者吻合得很好.对极限波浪运动的非线性和流域内速度分布进行了研究.     

Numerical Simulation of Waves Interaction with A Submerged Horizontal Twin-Plate Breakwater

The numerical investigation of regular waves interacting with a submerged horizontal twin-plate breakwater is presented in this paper. A numerical model with an absorbing wave-maker is established based on the VOF method. The validity of the model is verified by experimental results. Comparisons between the numerical and experimental results show that both the water surface profiles and the wave-induced pressures can be modeled accurately. Wave deformation over the breakwater,water particle velocities aroun...     

An Irregular Wave Maker of Active Absorption with VOF Method

A numerical irregular wave flume with active absorption of re-reflected waves is simulated by use of volume of fluid （VOF） method. An active ＇absorbing wave-maker based on linear wave theory is set on the left boundary of the wave flume. The progressive waves and the absorbing waves are generated simultaneously at the active wave generating-absorbing boundary. The absorbing waves are generated to eliminate the waves coming back to the generating boundary due to reflection from the outflow boundary and the structures. SIRW method proposed by Frigaard and Brorsen （1995） is used to separate the incident waves and reflected waves. The digital filters are designed based on the surface elevation signals of the two wave gauges. The corrected velocity of the wave-maker paddle is the output from the digital filter in real time. The numerical results of regular and irregular waves by the active absorbing-generating boundary are compared with the numerical results by the ordinary generating boundary to verify the performance of the active absorbing-generator boundary. The differences between the initial incident waves and the estimated incident waves are analyzed.     

Numerical Wave Flume Study on Wave Motion Around Submerged Plates

1 .IntroductionApile supportedplatesubmergedatacertaindepthunderseasurfacewasdevelopedasanewtypeofunderwaterbreakwaterfortheprotectionofcoastlinesandharbors .Thisisbecauseitdoesnothinderthewaterexchangebetweentheopenseaandtheprotectedareanordoesithindertheviewovertheopensea .Besides,itislessdependentonthegeotechnicalconditionsoftheseabottomwherethestructureistobeinstalled ;however,itscostishigh ,particularlyinrelativelydeepwaters .Formanyapplicationsitispossibletoreducethewavemotionintheprotec…     

三维完全非线性波浪水槽的数值模拟

用有限元求解拉普拉斯方程,建立了三维完全非线性数值波浪水槽.跟踪流体自由表面的方法为满足完全非线性自由表面条件的半拉格朗日法,对离散单元采用20节点的六面体二次等参数单元.并把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克斯波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播.     

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…     

Large-wave simulation of three-dimensional,cross-shore and oblique,spilling breaking on constant slope beach

In this work, the large-wave simulation (LWS) method is adapted for application in spilling wave breaking over a constant slope beach. According to LWS, large scales of velocities, pressure and free-surface elevation are numerically resolved, while the corresponding unresolved scale effects are taken into consideration by a subgrid scale (SGS) model for wave and eddy stresses. The model may be not fully applicable in very shallow water, close to the shoreline, where the unresolved, turbulent, free-surface oscillation is of the same order with the water depth. Time integration of the Euler equations is achieved by a two-stage fractional scheme, combined with a hybrid scheme for spatial discretization, consisting of finite difference and pseudospectral approximation methods. Model parameters are calibrated by comparison to available experimental data of free-surface elevation and velocities in the surf zone for cross-shore incoming waves. The action of the wave SGS stresses in the outer coastal and surf zones initiates breaking and generates appropriate vorticity, in the form of an eddy structure (surface roller), at the breaking wavefront. At incipient breaking, both advection and gravity contribute to the vorticity flux at the free surface, while only after the full development of the surface roller, the effect of advection becomes stronger. The SGS model is also utilized to simulate propagation, refraction and breaking of oblique incoming waves. The gradual breaking and dissipation of wave crestlines and the surface roller structure along the breaking wavefront are automatically captured without any empirical input, such as data for the roller shape or the wave propagation angle at breaking.