Extended mild-slope equation for random waves

A time-dependent extended mild-slope equation is derived from the elliptic equation of Chamberlain and Porter [J. Fluid Mech. 291 (1995) 393] using the Taylor series technique. Numerical tests are made on a horizontally one-dimensional case for regular waves over sloping beds and for both regular and irregular waves over a ripple patch. Numerical results prove that the proposed model gives accurate results for both regular and irregular waves over rapidly varying topography.     

Assessment of open boundary conditions on the elliptic formulation of the mild-slope equation

The performance of open boundaries in a finite differences scheme of the elliptic mild-slope equation is assessed. The wave propagation results show that lowest order parabolic radiation boundary conditions, unlike sponge layers combined with first order radiation boundary conditions, are an efficient alternative to first order radiation boundary conditions in order to improve the accuracy of the numerical solution of the problem.     

A method to determine the incident wave boundary conditions and its application

When studying the harbor water tranquility, cases are often confronted as that the verification point is not located on the generation line or that the angle between the generation line and the isobath is so large that the differences of the wave climates along the generation line can not be ignored. For these cases, the incident boundary conditions are difficult to evaluate. In order to solve this problem, a combined wave model is developed in the present paper based on the Boussinesq equation and the wave action balance equation. Instead of the one-line wave generation method, a multi-line generation method is proposed for the combined model. Application of this method is given to a case that the harbor is designed with two entrances and the angle between the generation line and the isobath is large and the results are shown reasonable. We suggest that the wave generation method on multi-lines might also be introduced to the wave physical model as the replacement for the one-line generation method.     

Simulation of small-amplitude frequency-dispersive transient waves by means of the mild-slope equation

A numerical model based on the mild-slope equation is applied to reproduce the propagation of small-amplitude transient waves. The model makes use of the Fourier Transform to convert the time-dependent hyperbolic equation into a set of elliptic equations in the frequency domain. The results of two available experimental studies on tsunamis generated by landslides are used to validate the model, which appears to be able of carefully reproducing the effects of the frequency dispersion. An example application of tsunamis propagating around the Stromboli island is also presented to show the applicability of the present approach to real life scenarios. It is finally discussed how this model could be applied as support to a tsunami early warning system.     

Modified mild slope equation and open boundary conditions

A numerical model to compute wave field is developed. It is based on the Berkhoff diffraction-refraction equation, in which an energy dissipation term is added, to take into account the breaking and the bottom friction phenomena. The energy dissipation function, by breaking and by bottom friction, is introduced in the Berkhoff equation to obtain a new equation of propagation.The resolution is done with the hybrid finite element method, where lagrangians elements are used.     

Improvement on open boundaries on a time dependent numerical model of wave propagation in the nearshore region

An improvement on the simulation of outgoing waves on a time dependent numerical model for water wave propagation in the nearshore region is presented. The governing equations consist of a system of first order partial differential equations (PDEs), the equation of continuity and the equation of motion. A comparative study of first order radiation boundary conditions (BCs) and first order radiation BCs combined with sponge layers is presented for cases where outgoing waves leave the numerical domain of calculation through the open boundary. A reduction of spurious reflections from the numerical open boundaries can be obtained with an irrelevant increase in terms of computational cost.     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。     

An energy-controlling boundary condition for partial wave reflections in the mild slope equation

An energy-controlling technique to actively manage the reflective property of waves from solid boundary is presented. As linear waves propagate through an energy-controlling area, a reduction in wave heights occurs due to energy dissipation, which can be placed under direct control through the imaginary part of the wavenumber and phase velocity. Based on this relationship, the present study investigates a new method to control reflected waves with desired heights in the mild slope equation model. The method is validated through numerical tests for various reflection coefficients and the results confirm the promising use of energy-controlling boundary condition for partial wave reflections.     

A new kind of nonlinear mild-slope equation for combined refraction-diffraction of multifrequency waves

This paper presents the numerical solution of a new nonlinear mild-slope equation governing waves with different frequency components propagating in a region of varying water depth. There are two new nonlinear equations. The linear part of the equations is the mild-slope equation, and one of the models has the same non-linearity as the Boussinesq equations. The new equations are directly applicable to the problems of nonlinear wave-wave interactions over variable depth. The equations are first simplified with the parabolic approximation, and then solved numerically with a finite difference method. The Crank-Nicolson method is used to discretize the models. The numerical models are applied to a set of published experimental cases, which are nonlinear combined refraction-diffraction with generation of higher harmonic waves. Comparison of the results shows that the present models generally predict the measurements better than other nonlinear numerical models which have been applied to the data set.     

Transient response of harbours to long waves under resonance conditions

This paper presents the results of a study aimed at quantifying the time–response of harbour basins to long waves under resonance conditions. On the basis of numerical simulations reproducing long waves in the yacht harbour of Rome (Ostia, Italy), it shows that the results valid for periodic forcing waves, acting for an infinitely long time, as those provided by models based on elliptic equations like the Helmoltz and the mild-slope equations, can be misleading with respect to the more realistic ones that can be obtained using time-varying wave equations. Taking advantage of the similarity between the processes studied here and a simple one-dimensional resonator, a method is also proposed to roughly estimate a time–response parameter of each mode of the harbour, using results from steady-state numerical model results, commonly applied for studying harbour resonance in engineering practice. On the basis of further numerical simulations, aimed at reproducing schematic harbour layouts, the effect on resonance of the position of the entrance and of an outer harbour is studied. The results indicate that the effects of design solutions to reduce resonance, by placing the entrance at the middle of the harbour, or using the outer harbour as a resonator, can be correctly evaluated only when considering the time needed for the oscillations to fully develop.     

Transparent boundary condition for simulating nonlinear water waves by a particle method

Open boundaries are important when simulating water waves. In this study, a transparent boundary condition at an open boundary was developed for simulating nonlinear water waves propagating to a distant area using the Moving Particle Semi-implicit method. The novelty of this study is that the technique of wave analysis used in the experiment was introduced into the particle simulation to absorb incident waves; the simulation cost was reduced by employing inflow and outflow regions instead of a long dissipation region. Incident waves in front of the boundary were evaluated using Fourier analysis, and the particles on the transparent boundary were forced to move at the velocity of the analytical solution for Stokes waves in order to absorb the incident waves. The analysis was restricted to periodic waves. Wave propagation was simulated for two wave periods using the developed transparent boundary condition. The results showed that this transparent boundary transmitted the incident waves with small reflection and the simulation cost was lower than that for wave damping by a conventional highly viscous region.     

Three-dimensional finite difference model for transport of conservative pollutants

A three-dimensional finite difference transport model appropriate for the coastal environment is developed for the solution of the three-dimensional convection-diffusion equation. A higher order upwind scheme is used for the convective terms of the convection-diffusion equation, to minimise the numerical diffusion. The validity of the numerical model is verified through five test problems, whose exact solutions are known.     

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.     

Generation of random waves in time-dependent extended mild-slope equations using a source function method

We develop techniques of numerical wave generation in the time-dependent extended mild-slope equations of Suh et al. [1997. Time-dependent equations for wave propagation on rapidly varying topography. Coastal Engineering 32, 91–117] and Lee et al. [2003. Extended mild-slope equation for random waves. Coastal Engineering 48, 277–287] for random waves using a source function method. Numerical results for both regular and irregular waves in one and two horizontal dimensions show that the wave heights and the frequency spectra are properly reproduced. The waves that pass through the wave generation region do not cause any numerical disturbances, showing usefulness of the source function method in avoiding re-reflection problems at the offshore boundary.     

适于模拟不规则水域波浪的缓坡方程两种数值模型比较

本文分析比较了适于不规则水域波浪模拟的椭圆型缓坡方程两种数值模型。两种数值模型均采用有限体积法离散,分别基于四叉树网格和非结构化三角形网格建立。首先结合近岸缓坡地形上波浪传播的经典物理模型实验对两种数值模型分别进行了验证,并结合计算结果对比分析了两种模型的计算精度和效率。计算结果表明,两种数值模型均可有效地模拟近岸波浪的传播变形;相对非结构化三角形网格下的模型,基于四叉树网格建立的数值模型在数值离散和求解过程中无需引入形函数、不产生复杂的交叉项,离散简单,易于程序实现,且节约计算存储空间,计算效率高。     

考虑边界波浪方向的缓坡方程自适应求解模型

缓坡方程是描述近岸波浪运动较好的数学模型之一。在发展的自适应有限元求解缓坡方程的基础上,采用迭代求解的方法,确定波浪相对于边界的入射方向,从而对边界条件进行改进,建立了求解缓坡方程的数值计算模型。典型算例表明,考虑波浪相对于边界的入射角度后,模型可以更好地模拟吸收波浪边界,同时对多向波对双突堤的绕射进行了模拟研究,与试验结果比较表明,所建立的数值计算模型能够适用于多向不规则波传播过程的模拟研究。     

基于双层Boussinesq方程的三维波浪数值模型及其验证

Liu等给出的最高导数为2的双层Boussinesq水波方程具有较好的色散性和非线性,基于该方程建立了有限差分法的三维波浪数值模型。在矩形网格上对方程进行了空间离散,采用高阶导数近似方程中的时、空项,时间积分采用混合4阶Adams-Bashforth-Moulton的预报—校正格式。模拟了深水条件下的规则波传播过程,计算波面与解析结果吻合较好,反映出数值模型能很好地刻画波面过程及波面处的速度变化;在kh=2π条件下可较为准确获得沿水深分布的水平和垂向速度,这与理论分析结果一致。最后,利用数值模型计算了规则波在三维特征地形上的传播变形,数值结果和试验数据吻合较好;高阶非线性项会对波浪数值结果产生一定的影响,当波浪非线性增强,水深减少将产生更多的高次谐波。建立的双层Boussinesq模型对强非线性波浪的演化具有较好的模拟精度。     

Numerical simulation of three-dimensional sloshing phenomena using a finite difference method with marker-density scheme

In the present study, three-dimensional sloshing phenomena occurring in liquid cargo tanks are numerically simulated. The Navier-Stokes equations and continuity equation are used for the governing equations, and solved with a finite difference method in a rectangular fixed staggered mesh system. The positions of free surface are defined by the Marker density method satisfying the free-surface boundary conditions, and the flows of the gas and liquid regions are simulated simultaneously. The irregular leg length and star method is employed on the cells near the free surface for the computations of pressure. The computation results are compared with other experimental results to verify the consistency of the present numerical method, and the agreements are reasonably good. Furthermore, the flow characteristics inside a partially filled liquid tank of a real sized ship oscillating regularly and irregularly are computed to verify the possibility of practical application of the present method.     

A finite difference solution of the shear flow over a circular cylinder

The incompressible viscous shear flow past a circular cylinder is analyzed by solving two-dimensional Navier-Stokes equations and pressure Poisson equation using a finite difference method. The shear flow is calculated for Reynolds numbers from 80 to 1000, and shear parameters up to 0.25. The numerical results indicate that the vortex shedding persists at the shear parameters up to 0.25 for the present Reynolds number range. The Strouhal number and the drag coefficient decrease as the shear parameter increases. There is a transverse force acting from the high velocity side toward the low velocity side in the shear flow.     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.