Application of desingularized approach to water wave propagation over three-dimensional topography

A numerical approach based on desingularized boundary element method and mixed Eulerian–Lagrangian formulation [Zhang et al., 2006. Wave propagation in a fully nonlinear numerical wave tank: a desingularized method. Ocean Engineering 33, 2310–2331] is extended to solve the water wave propagation over arbitrary topography in a three-dimensional wave tank. A robust damping layer applicable for regular and irregular incident waves is employed to minimize the outgoing wave reflection back into the wave tank. Numerical results on the propagation of regular and irregular incident waves over the flat bottom and linear incident waves over an elliptical shoal show good concurrence with the corresponding analytical solutions and experimental data.     

Linear and nonlinear irregular waves and forces in a numerical wave tank

A time-domain higher-order boundary element scheme was utilized to simulate the linear and nonlinear irregular waves and diffractions due to a structure. Upon the second-order irregular waves with four Airy wave components being fed through the inflow boundary, the fully nonlinear boundary problem was solved in a time-marching scheme. The open boundary was modeled by combining an absorbing beach and the stretching technique. The proposed numerical scheme was verified by simulating the linear regular and irregular waves. The scheme was further applied to compute the linear and nonlinear irregular wave diffraction forces acting on a vertical truncated circular cylinder. The nonlinear results were also verified by checking the accuracy of the nonlinear simulation.     

Bragg reflection in a fully nonlinear numerical wave tank based on boundary integral equation method

We investigated the phenomenon of Bragg reflection of submerged structures in a 2D fully nonlinear numerical wave tank (NWT) based on the boundary integral equation method (BIEM). This model was validated by comparing not only the free surface elevations with that of the analytic solution of Stokes’ second-order wave theory, but also the reflection coefficients of submerged bars with that from other sources. The results of the present model show that the free surface nonlinear effect on the reflection coefficient of the primary resonance reduces significantly for all of the submerged bars considered. Finally, a case study is presented to demonstrate the reflecting capacity and overall performance of various submerged bars. Results indicate that sinusoidal bar has the maximum reflection capacity at the primary resonance, but the trapezoidal submerged bar is suggested as the better option for the practical convenience of coastal underwater construction.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

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

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

数值波浪水槽中完全非线性浅水波仿真特性研究

应用基于势流理论的时域高阶边界元方法,建立一个完全非线性的三维数值波浪水槽,通过实时模拟推板造波运动的方式产生波浪。通过混合欧拉-拉格朗日方法和四阶Runge-Kutta方法更新自由水面和造波板的瞬时位置。利用所建模型分别模拟了有限水深波和浅水波,与试验结果、相关文献结果和浅水理论结果吻合较好,且波浪能够稳定传播。系统地讨论造波板的运动圆频率、振幅和水深等对波浪传播和波浪特性的影响,并对波浪的非线性特性进行分析,研究发现造波板运动频率、运动振幅以及水深均将对波浪形态和波浪非线性产生显著影响。结果为真实水槽造波机的运动控制以及波浪生成试验提供了依据,便于实验室设置更合理的参数来准确模拟不同条件下的波浪。     

一个二维数值波浪水槽的改进

对Drimer,Agnon及Segre研发的数值波浪水槽(简称DAS),首先通过修正DAS中某些积分计算的错误,建立了第一个修正模型(MDAS1),使得模型更加准确稳定。其次,将DAS模型自由表面线性元近似替换为更为合理准确的三阶元近似,建立了第二个修正模型(MDAS2),使得模型计算的波面更合理可靠。平底水槽中波群的生成及传播算例表明:与DAS结果相比,MDAS1对波面没有显著影响,而采用三阶元的MDAS2对波面有显著影响。与Hansen和Svendsen的实测资料相比较证实了两个修正模型的有效性。     

数值水槽内的阻尼消波和波浪变形计算

采用时域内对波面运动位置追踪的边界元方法,建立了一种非线性波浪变形计算的三维数值模式,并借助Spongelayer阻尼消波和Sommerfeld放射条件相匹配的处理方式,开发了一条三维非线性波的数值造波水槽,进而对水槽内的Stokes波进行了波浪变形计算。     

完全非线性深水波的数值模拟

基于势流理论,并结合深水波质点运动从水面向下呈e指数衰减的特性,建立了完全非线性数值变深水槽模型,通过实时模拟活塞式造波机运动来产生波浪.采用时域高阶边界元法进行模拟,利用混合欧拉-拉格朗日方法和四阶Runge-Kutta方法追踪流体瞬时水面,应用镜像格林函数消除了水槽两个侧面的积分,在水槽末端布置人工阻尼层来消除反射...     

基于高阶边界元的三维数值波浪港池

初步建立了一个基于高阶边界元的三维数值波浪港池,港池具有造波和消波功能。采用高阶边界元16节点四边形单元和基于二阶显式泰勒展开的混合欧拉-拉格朗日时间步进求解带自由表面的完全非线性势流方程。模型中对于影响数值精度的问题作了细致的处理。数值计算结果表明本港池可以用来模拟非线性波浪的传播,具有很高的数值精度和稳定性。     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

A numerical study on local characteristics of predetermined irregular wave trains

In actual sea states, damage to offshore floating structures is usually caused by a few extreme waves or wave groups in an irregular wave train. Accurate simulation of the irregular wave trains can lay a solid foundation for understanding the local flow field and impact loads that would potentially cause such damage. This paper describes how the generation of a single extreme wave was investigated. Determination of the wave-maker motion for generating specified irregular wave trains is the key to this work. First, an experimental irregular wave train was decomposed into a certain number of small-amplitude waves. Fourier series expansion was performed to determine the amplitude and the initial phase angle of each wave component. Then a hydrodynamic transfer function was used to calculate the amplitude of the wave-maker motion associated with each wave component. Superposition was made on all the wave components to get the final wave-maker motion. During the numerical simulation, calculated horizontal velocity profiles of the extreme wave at different moments were analyzed and compared with experimental results, and a satisfactory agreement was obtained. In the simulation, VOF method was employed to capture the free surface, and a dissipation zone was used to deal with wave reflection.     

Direct time domain analysis of floating structures with linear and nonlinear mooring stiffness in a 3D numerical wave tank

In this paper, motion response of a moored floating structure interacting with a large amplitude and steep incident wave field is studied using a coupled time domain solution scheme. Solution of the hydrodynamic boundary value problem is achieved using a three-dimensional numerical wave tank (3D NWT) approach based upon a form of Mixed-Eulerian–Lagrangian (MEL) scheme. In the developed method, nonlinearity arising due to incident wave as well as nonlinear hydrostatics is completely captured while the hydrodynamic interactions of radiation and diffraction are determined at every time step based on certain simplifying approximations. Mooring lines are modelled as linear as well as nonlinear springs. The horizontal tension for each individual mooring line is obtained from the nonlinear load-excursion plot of the lines computed using catenary theory, from which the linear and nonlinear line stiffness are determined. Motions of three realistic floating structures with different mooring systems are analyzed considering various combinations of linear and approximate nonlinear hydrodynamic load computations and linear/nonlinear mooring line stiffness. Results are discussed to bring out the influence and need for consideration of nonlinearities in the hydrodynamics and hydrostatics as well as the nonlinear modelling of the line stiffness.     

Interactions between fully nonlinear water waves and cylinder arrays in a wave tank

Fully nonlinear interactions between water waves and vertical cylinder arrays in a numerical tank are studied based on a finite element method (FEM). The three-dimensional (3D) mesh is constructed through an extension of a 2D Delaunay surface grid along the vertical line. The velocity potential is obtained by solving a linear matrix system of FEM, and a difference scheme is then used to calculate the velocity on the free surface to track its movement. Waves and hydrodynamic forces are obtained for both bottom mounted and truncated cylinders. The simulations have provided many results to show the nature of mutual interference between cylinders in arrays and its effects on waves and forces at the nearly trapped mode frequency. The effect of the tank wall on waves and forces has been investigated, and the nonlinear features of waves and forces have also been discussed.     

圆形数值波浪水池斜向和多向入射波浪模拟研究

为了研究斜向入射波浪,基于三维不可压缩两相流模型,开发了一套圆形数值波浪水池数值模型。在圆形波浪水池中,通过源项造波法成功生成了任意入射方向的波浪,并且利用人工摩擦项模拟阻尼区以数值耗散反射波浪。模型基于嵌入式多块网格体系,采用FVM法(finite volume method)离散Navier-Stokes方程,VOF法(volume of fluid)追踪自由水面。试验结果表明,斜向入射波浪的模拟结果与理论值基本一致,圆形波浪水池在模拟斜向入射波浪时,有效区域的面积较传统波浪水池显著增大,而且有效区域受波浪入射角度的影响也较小。同时,通过叠加多列斜向入射波浪,模拟出了多向交叉波列,并通过与理论结果对比,发现其具有较高的精度。     

Wave prediction in a port using a fully nonlinear Boussinesq wave model

A wave forecasting system using FUNWAVE-TVD which is based on the fully nonlinear Boussinesq equations by Chen(2006) was developed to provide an accurate wave prediction in the Port of Busan, South Korea. This system is linked to the Korea Operational Oceanographic System(KOOS) developed by Park et al.(2015). The computational domain covers a region of 9.6 km×7.0 km with a grid size of 2 m in both directions, which is sufficient to resolve short waves and dominant sea states. The total number of grid points exceeds 16 millions,making the model computational expensive. To provide real-time forecasting, an interpolation method, which is based on pre-calculated results of FUNWAVE-TVD and SWAN forecasting results at the FUNWAVE-TVD offshore boundary, was used. A total of 45 cases were pre-calculated, which took 71 days on 924 computational cores of a Linux cluster system. Wind wave generation and propagation from the deep water were computed using the SWAN in KOOS. SWAN results provided a boundary condition for the FUNWAVE-TVD forecasting system. To verify the model, wave observations were conducted at three locations inside the port in a time period of more than 7 months. A model/model comparison between FUNWAVE-TVD and SWAN was also carried out. It is found that, FUNWAVE-TVD improves the forecasting results significantly compared to SWAN which underestimates wave heights in sheltered areas due to incorrect physical mechanism of wave diffraction, as well as large wave heights caused by wave reflections inside the port.     

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

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

Determining the onset and strength of unforced wave breaking in a numerical wave tank

A numerical wave tank is used to investigate the onset and strength of unforced wave breaking, and the waves have three types of initial spectra： constant amplitude spectrum, constant steepness spectrum and Pierson-Moscowitz spectrum. Numerical tests are performed to validate the model results. Then, the onset of wave breaking is discussed with geometric, kinematic, and dynamic breaking criteria. The strength of wave breaking, which is always characterized by the fractional energy loss and breaking strength coefficient, is studied for different spectra. The results show how the energy growth rate is better than the initial wave steepness on estimating the fractional energy losses as well as breaking strength coefficient.     

一维全非线性Green-Naghdi水波方程数值模型

建立了求解一维全非线性Green-Naghdi水波方程的中心有限体积/有限差分混合数值格式。采用结构化网格对守恒形式的控制方程进行离散和积分,界面数值通量采用有限体积法计算,剩余项则采用中心有限差分格式求解。其中,采用中心迎风有限体积格式计算控制体界面数值通量,并结合界面变量的线性重构方法,使其在空间上具有四阶精度,通过引入静压重构技术和波浪破碎指标使模型具备处理海岸水-陆动边界及波浪破碎的能力。时间积分则采用具有总时间变差减小(Total Variation Diminishing,TVD)性质的三阶龙格-库塔法进行。应用该模型对孤立波在常水深和斜坡海岸上的传播过程及规则波跨越潜堤传播的实验进行了数值模型研究,数值计算同解析解及实验数据吻合良好。