Wash Waves Generated by High Speed Displacement Ships

It is difficult to compute far-field waves in a relative large area by using one wave generation model when a large calculation domain is needed because of large dimensions of the waterway and long distance of the required computing points. Variation of waterway bathymetry and nonlinearity in the far field cannot be included in a ship fixed process either. A coupled method combining a wave generation model and wave propagation model is then used in this paper to simulate the wash waves generated by the passing ship. A NURBS-based higher order panel method is adopted as the stationary wave generation model; a wave spectrum method and Boussinesq-type equation wave model are used as the wave propagation model for the constant water depth condition and variable water depth condition, respectively. The waves calculated by the NURBS-based higher order panel method in the near field are used as the input for the wave spectrum method and the Boussinesq-type equation wave model to obtain the far-field waves. With this approach it is possible to simulate the ship wash waves including the effects of water depth and waterway bathymetry. Parts of the calculated results are validated experimentally, and the agreement is demonstrated. The effects of ship wash waves on the moored ship are discussed by using a diffraction theory method. The results indicate that the prediction of the ship induced waves by coupling models is feasible.     

关于Boussinesq型水波方程理论和应用研究的综述

Boussinesq型方程是研究水波传播与演化问题的重要工具之一,本文就1967-2018年常用的Boussinesq型水波方程从理论推导和数值应用两个方面进行了回顾,以期推动该类方程在海岸(海洋)工程波浪水动力方向的深入研究和应用。此类方程推导主要从欧拉方程或Laplace方程出发。在一定的非线性和缓坡假设等条件下,国内外学者建立了多个Boussinesq型水波方程,并以Stokes波的相关理论为依据,考察了这些方程在相速度、群速度、线性变浅梯度、二阶非线性、三阶非线性、波幅离散、速度沿水深分布以及和(差)频等多方面性能的精度。将Boussinesq型水波方程分为水平二维和三维两大类,并对主要Boussinesq型水波方程的特性进行了评述。进而又对适合渗透地形和存在流体分层情况下的Boussinesq型水波方程进行了简述与评论。最后对这些方程的应用进行了总结与分析。     

Numerical Wave Flume Study on Wave Motion Around Submerged Plates

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

Boussinesq相位解析的海岸水动力学数学模型研究进展

由于在平衡计算效率和精度上具有优势,Boussinesq相位解析数学模型研究不断取得突破,已成为波浪和水流精细化模拟的较优解析方式,为海岸工程、环境、地质等问题提供了实用和高效的研究手段。本文对已有Boussinesq类模型的研究进行了评述,深入探讨其重要发展、实际应用和理论瓶颈,从高阶非静压修正、GPU准三维高性能算法编译、波浪破碎和泥沙运移沉积等4个方面提出未来可能的科学突破方向。     

应用Boussinesq方程对由航道开挖所造成水深变化对波浪传播所产生的异常波况的数值研究

应用Boussinesq方程采用有限差分法建立了非线性数值波浪模式，并应用该模式对由航道开挖所造成水深变化对波浪传播产生的异常现象进行了数值模拟研究。从结果来看这种异常波浪局部增大现象是由于入射波与航道夹角过小，使波浪无法折射入航道，而在航道两侧反射叠加的结果。因此说水深的变化对波浪产生的反射现象也是不可忽视的。     

不规则波浪数值水槽的造波和阻尼消波

利用MAC(marker and cell)法直接数值求解连续方程和N-S方程.为模拟不规则波长时段连续造波及消除波浪遇结构物后形成的二次反射,采用了源造波法.对开敞边界,采用了海绵阻尼消波和Sommerfeld条件相结合的处理方式.     

浅水中浮体在波浪作用下运动的三维时域计算模型

港口中系泊船在波浪作用下运动问题的本质是浅水波浪与浮体的相互作用。与深水情况不同,浅水问题应当考虑水底、水域边界的影响及浅水波浪自身的特性,单一模型很难实现该模拟过程。为此,建立了Boussinesq方程计算入射波和Laplace方程计算散射波的全时域组合计算模型。有限元法求解的Boussinesq方程能使入射波充分考虑到水底、水域边界的影响和浅水波浪的特性;散射波被线性化,采用边界元法求解,并以浮体运动时的物面条件为入射波和散射波求解的匹配条件。该方法为完全的时域方法,计算网格不随时间变动,计算过程较为方便。通过与实验及其他数值方法的结果进行比较,验证了本模型对非线性波面、浮体的运动都有比较理想的计算结果,显示了本模型对非线性问题具有较好的计算能力。     

Applicability of numerical models to nonlinear dispersive waves

The applicability of three different wave-propagation models in nonlinear dispersive wave fields has been investigated. The numerical models tested here are based on three different wave theories: a fully nonlinear potential theory, a Stokes second-order theory, and a Boussinesq-type theory with an improved dispersion relation. Physical experiments and computations were conducted for wave evolutions during passage over a submerged shelf under various wave conditions. As expected, the fully nonlinear solutions agree better with the measurements than do the other solutions. Although the second-order solution has sufficient accuracy for smaller-amplitude wave cases, the truncation after the third harmonics causes significant discrepancies in wave form for larger waves. In addition, the second-order model markedly overestimates the first- and second-harmonic amplitudes in transmitted waves. The Boussinesq model provides excellent predictions of wave profile over the shelf even in larger wave cases. However, this model also overestimates the magnitudes of several higher harmonics in transmitted waves. These facts may indicate that energy transfer from bound components into free waves in these higher harmonics cannot be accurately evaluated by the Boussinesq-type equations.     

一种新型二阶全非线性Boussinesq方程及其数值验证

通过改进二阶全非线性 Boussinesq 波浪方程中的色散项,得到了一组没有改变原方程的数学形式但适用于更大变化水深的新方程,其色散性能和变浅性能都比原方程有了很大改进,所适用的水深范围更大,能更好地描述从深水到近岸浅水处的波浪传播;并基于新方程建立了波浪数值模型,通过模拟波浪从浅水到深水的传播变形来验证新方程的有效性.     

A mathematical model of calculating local tsunami wave source constraints and propagation

This paper presents a local tsunami simulation, including the initial displacement field model of tsunami source and tsunami wave propagation model. We deduced the tsunami wave equation; applied the ma...     

An analytic investigation of oscillations within a harbor of constant slope

Longitudinal and transverse oscillations within a harbor of constant slope are analyzed. Based on the linear shallow water approximation, longitudinal oscillations are described with Bessel equations. Ignoring friction, oscillations are forced using the period of the incident perpendicular wave field by the method of matched asymptotics. The analytic results show that the varying depth shifts the resonant wave numbers to lower values than those for the same geometric harbor with constant depth. Furthermore, we extend the shallow water equations to a linear, weakly dispersive, Boussinesq-type equation by modifying the offshore velocity component, and then use it to investigate possible existing transverse oscillations in the harbor of constant slope. These oscillations are types of standing edge waves. Their character is quite sensitive to the boundary condition at the backwall of the harbor.     

A coupled numerical model of wave interaction with porous medium

A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a kε closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses.     

Numerical modeling of nonlinear water waves over heterogeneous porous beds

The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup.     

实验室造波条件对内孤立波发展影响的直接数值模拟

内孤立波具有振幅尺度大、能量集中的特点,其引起流场和密度场的迅速变化可能对海洋工程结构物以及水下潜体造成严重威胁。因此研究不同造波条件下生成的内孤立波运动的流场特征具有重要的学术意义和实际应用价值。采用直接数值模拟方法和给定的初始密度场密度跃迁函数,对重力塌陷激发内孤立波的运动过程进行研究,探讨了不同造波条件下,激发产生的内孤立波波型、涡度、振幅和水平速度等流场特征。结果表明:(1)直接模拟数值方法能够模拟内孤立波传播过程中的密度界面波型反转现象;(2)从定性和定量的角度,证实了不稳定内孤立波传播过程中存在能量的向后传递;(3)对于相同的台阶深度(水闸两侧初始密度界面的高度差),初始涡流保持相同,但是随着上下层水深比的减小,其强度下降显著;(4)台阶深度对初始涡流的垂直结构的影响要大于上下层水深比,且台阶深度对内孤立波的振幅、水平速度的影响显著。     

Internal wave generation in an improved two-dimensional Boussinesq model

A set of Boussinesq-type equations with improved linear frequency dispersion in deeper water is solved numerically using a fourth order accurate predictor-corrector method. The model can be used to simulate the evolution of relatively long, weakly nonlinear waves in water of constant or variable depth provided the bed slope is of the same order of magnitude as the frequency dispersion parameter. By performing a linearized stability analysis, the phase and amplitude portraits of the numerical schemes are quantified, providing important information on practical grid resolutions in time and space. In contrast to previous models of the same kind, the incident wave field is generated inside the fluid domain by considering the scattered wave field in one part of the fluid domain and the total wave field in the other. Consequently, waves leaving the fluid domain are absorbed almost perfectly in the boundary regions by employment of damping terms in the mass and momentum equations. Additionally, the form of the incident regular wave field is computed by a Fourier approximation method which satisfies the governing equations accurately in water of constant depth. Since the Fourier approximation method requires an Eulerian mean current below wave trough level or a net mass transport velocity to be specified, the method can be used to study the interaction of waves and currents in closed as well as open basins. Several computational examples are given. These illustrate the potential of the wave generation method and the capability of the developed model.     

Using RANS to simulate vortex generation and dissipation around impermeable submerged double breakwaters

In this paper, a Reynolds Averaged Navier–Stokes (RANS) model was developed to simulate the vortex generation and dissipation caused by progressive waves passing over impermeable submerged double breakwaters. The dynamics of the turbulence are described by introducing a k–ɛ model with Boussinesq closure. The Height Function (HF) is implemented to define the free-surface configuration. The governing equations are discretized by means of a finite volume method based on a staggered grid system with variable width and height. The feasibility of the numerical model was verified through a series of comparisons of numerical results with the existing analytical solutions and the experimental data. The good agreements demonstrate the satisfactory performance of the developed numerical model. The flow separation mechanism both near the upstream and the downstream edges of the obstacles demonstrates the physical and expected nature of development of the flow. The present model provides an accurate and efficient tool for the simulation of flow field and wave transformation near coastal structures without breaking.     

波浪与大孔隙多孔介质相互作用的耦合数学模型

建立了波浪与大孔隙多孔介质相互作用的耦合数学模型,波浪域的控制方程为雷诺时均方程和k-ε紊流模型。对于计算域的入射波采用推板式造波,它可以是线性波、椭圆余弦波和孤立波。采用PLIC-VOF法追踪波浪自由表面。对于多孔介质内的孔隙流场采用非线性Forchheimer方程,两区域共享连续方程,最后导出的波浪域与孔隙流域的压力修正方程具有完全相同的形式,利用这个方程能够同时而不是分别求解波浪场和孔隙流场,避免了在内部边界上给定匹配条件,实现了波浪场与孔隙流场的同步耦合。波浪与粗颗粒海床、平底床面上抛石潜堤及斜坡上抛石潜堤相互作用的验证计算结果表明该模型可用于研究波浪与大孔隙多孔介质相互作用的问题。     

完全非线性孤立波的直墙反射

报道了应用边界积分方法模拟完全非线性孤立波的传播与直墙反射,给出了波形演变过程。结果表明,本模型对计算孤立波的传播与直墙反射是有效的。三阶Boussinesq方程的孤立波解比低阶方程的孤立波解更接近完全非线性的数值解.当来波波高增大时,孤立波直墙反射的相位滞后变小。若考虑大波高孤立波的直墙反射或波——波相互作用,一阶理论预报的相位滞后往往低估实际情况。