On waves propagating over poro-elastic seabed

In this study, an analytical solution is developed for the problem of periodic waves propagating over a poro-elastic seabed of infinite depth. Water waves above the seabed are described using the linear wave theory. The poro-elastic seabed is modelled based on the Biot theory in which the inertia effect and Darcy's friction are added. Continuity of dynamic pressure and flow flux at the interfacial seabed surface are considered. Adopting an approach similar to Hsu et al. (1993), the governing equations for the pore pressure and displacements of the poro-elastic medium are derived. The present analytic solution compares favorably well with experimental results by Yamamoto et al. (1978), and analytical results by Song (1993) for the case of fine sand. Using the present theory, variations of the wavelength and fluid pressure caused by coupling of waves and the poro-elastic seabed are discussed. Results show that higher elasticity of the poro-elastic seabed induces larger interface pressure, but higher permeability causes smaller pressure on the seabed interface. The wave length is affected by the poro-elastic seabed and becomes shorter for softer seabed and shallower water depth.     

Numerical calculation of forces induced by short-crested waves on a vertical cylinder of arbitrary cross-section

Most off-shore oil platforms are supported by vertical cylinders extending to the ocean floor. An important problem in off-shore engineering is the calculation of the wave loading exerted on these vertical cylinders. Analytical solutions have been found for the case of plane incident waves incident on a circular cylinder by MacCamy and Fuchs [(1954), Wave forces on piles: a diffraction theory. U.S. Army Corps of Engineering, Beach Erosion Board, Technical Memorandum No. 69] and also for short-crested waves incident on a circular cylinder by Zhu [(1993), Diffraction of short-crested waves around a circular cylinder. Ocean Engng 20, 389–407]. However, for a cylinder of arbitrary cross-section, no analytic solutions currently exist. Au and Brebbia [(1983), Diffraction of water waves for vertical cylinders using boundary elements. Appl. Math. Modelling 7, 106–114] proposed an efficient numerical approach to calculate the wave loads induced by plane waves on vertical cylinders by using the boundary element method. However, wind-generated waves are better modelled by short-crested waves. Whether or not these short-crested waves can induce larger wave forces on a structure is of great concern to ocean engineers. In this paper wave loads, induced by short-crested incident waves, on a vertical cylinder of arbitrary cross-section are discussed. For a cylinder of certain cross-section, the wave loads induced by short-crested waves can be larger than those induced by plane waves with the same total wave number.     

An analytical solution for wave-induced seabed response in a multi-layered poro-elastic seabed

In this study, a new analytical solution for the wave-induced seabed response in a multi-layered poro-elastic seabed is developed. The seabed is treated as a multi-layered porous medium and characterized by Biot’s theory. The displacements of the solid skeleton and the pore pressure are expressed in terms of two scalar potentials and one vector. Then, the Biot’s dynamic equation can be solved using Fourier transformation and reducing to Helmholtz equations. To obtain the general solutions for the multi-layered poro-elastic seabed in the frequency-wave-number domain, the transmission and reflection matrices (TRM) method is used to form the equivalent stiffness. Using the boundary conditions and continuous conditions, the frequency-wave-number domain solutions are obtained. Finally, the time-space domain solutions for the multi-layered poro-elastic seabed are obtained by means of the inverse Fourier transformation with respect to the horizontal coordinate. Based on the new solution, a parametric study is carried out to examine the effects of soil characteristics (number of layers, permeability and shear modulus) and wave characteristics (water depth and wave steepness) on seabed responses. The results indicate that the seabed response is affected significantly by permeability, shear modulus and relative water depth.     

Scattering of oblique waves by permeable vertical flexible membrane wave barriers

The interaction of obliquely incident surface gravity waves with a vertical flexible permeable submerged membrane wave barrier is investigated in the context of three-dimensional linear water wave theory. From the general formulation of the submerged membrane barrier, the performance of bottom-standing, surface-piercing and fully extended membrane wave barriers are analyzed for various values of wave and structural parameters. The analytic solution of the physical problem is obtained using eigenfunction expansion method and a coupled boundary element-finite difference method has been used to get the numerical solution. In the boundary element method, since the boundary condition on the membrane barrier is not known a priori, the membrane response and velocity potentials are solved simultaneously using appropriate discretization with the help of finite difference scheme. The convergence of the analytic and numerical solution techniques is discussed. The study reveals that for suitable combination of wave and structural parameters, approximately (45–50)% incident wave energy can be dissipated irrespective of membrane barrier configurations. Further, in certain situations, nearly full wave reflection and zero transmission occur for all barrier configurations. The study will be useful in the design of flexible permeable membrane to act as an effective wave barrier for creation of tranquility zone in the marine environment.     

淹没矩形防波堤透反射系数特性研究

采用解析方法研究了斜向入射波作用下淹没矩形防波堤的透反射系数特性.首先利用特征函数展开法导出了绕射势函数的分析解和透反射系数的计算公式,然后利用边界元方法验证了解析解,在此基础上利用解析解分析了若干工况下的防波堤透反射特性.计算结果表明,淹没矩形防波堤截面的宽度、高度和相对位置以及入射角的改变都不同程度影响反射系数和透射系数.在中等深度条件下,对于一定频率的波浪,位置和尺寸适当的淹没矩形堤可以反射大部分斜向入射波.研究结果对设计淹没的矩形防波堤具有重要的参考价值.     

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

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

The prediction of the dynamic and structural motions of a floating-pier system in waves

In this study, a two-dimensional floating pier consists of single rectangular impermeable pontoon with side supporting pile-columns is studied. The purpose of this study is to present a theoretical solution for the linearized problem of incident waves exerting on a floating pier with pile-restrained. All boundary conditions are linearized in the problem, which is incorporated into a scattering problem and radiation problem with unit displacement. The method of separation of variables is used to solve for velocity potentials. For the radiation problem with unit heave and pitch amplitude, the boundary value problem with non-homogeneous boundary condition beneath the structure is solved by using a solution scheme. By calculating the wave force from velocity potential and solving the equation of motion of the floating structure simultaneously a close form theoretical solution for the problem is developed. The finite element method was also applied to calculate the dynamic responses on the supporting piles subjected to the pontoon motions and incident waves.     

Trapping mechanism of submerged ridge on trans-oceanic tsunami propagation

Based on the linear shallow water equations,an analytic solution of trapped waves over a symmetric parabolicprofile submerged ridge is derived.The trapped waves act as propagating waves along the ridge and as standing waves across the ridge.The amplitude gets the maximum at the ridge top and decays gradually towards both sides.The decaying rate gets more gently with higher modes.Besides,an explicit first-order approximate dispersion relation is derived to simplify transcendental functions in the exact solution,which is useful to describe trapped waves over shallowly submerged ridges in reality.Furthermore,the trapping mechanism of the submerged ridge waveguides on the trans-oceanic tsunami propagation can be explained by the ray theory.A critical incident angle exists as a criterion to determine whether the wave is trapped.Besides,a trapped parameter γ is proposed to estimate the ratio of the energy trapped by the oceanic ridge if a tsunami is generated at its top.     

Wave Radiation by A Submerged Ring Plate in Water of Finite Depth

A plate submerged at a certain depth underneath the sea surface has been proposed as a structure type for different purposes, including motion response reduction, wave control, and wave energy harvesting. In the present study, the three-dimensional wave radiation problem is investigated in the context of the linear potential theory for a submerged ring plate in isolation or attached to a floating column as an appendage. In the latter case, the ring plate is attached at a certain distance above the column bottom. The structure is assumed to undergo a heave motion. An analytical model is developed to solve the wave radiation problem via the eigenfunction expansion method in association with the region-matching technique. With the velocity potential being available, the hydrodynamic coefficients, such as added mass and radiation damping, are obtained through the direct pressure integration. An alternative solution of radiation damping has also been developed in this study, in which the radiation damping is related to the Kochin function in the wave radiation problem. After validating the present model, numerical analysis is performed in detail to assess the influence of various plate parameters, such as the plate size and submergence depth. It is noted that the additional added mass due to the attached ring plate is larger than that when the plate is in isolation. Meanwhile, the radiation damping of the column for the heave motion can vanish at a specific wave frequency by attaching a ring plate, corresponding to a condition that there exist no progressive waves in the exterior region.     

Scattering of obliquely incident water waves by partially reflecting non-transmitting breakwaters

In the present paper, analytic solutions are derived for scattering of water waves obliquely incident to a partially reflecting semi-infinite breakwater or breakwater gap. In order to examine the correctness of the derived solutions, they are compared with the solutions derived by McIver (1999) and Bowen and McIver (2002) for a semi-infinite breakwater and a breakwater gap, respectively, in the case of perfect reflection. The derived analytic solutions are used to investigate the effect of reflection coefficient of the breakwater and wave incident angle upon the tranquility at harbor entrance. The tranquility is deteriorated by the reflected waves as the reflection coefficient increases and as waves are incident more obliquely.     

Wave diffraction by elliptical breakwaters in shallow water

The scattering of waves by both floating and submerged stationary elliptical breakwaters is investigated by means of linearised shallow water wave theory. This formulation leads to solutions for the fluid velocity potential in terms of Mathieu functions of real argument. Expressions are derived for the wave-induced forces and moments on the structures and their total and differential scattering cross-sections. Numerical results are presented for a range of wave and structural parameters.The present analysis serves as a prelude to a more comprehensive study of the problem without the shallow water restriction.     

A theory for waves interacting with porous structures with multiple regions

This study presents an analytical solution for the problem of waves passing a submerged porous structure, using a multi-region method in the solution scheme considering the characteristics of geometry and composing materials of the porous structure. Using the flux and pressure conditions on horizontal boundaries and interfaces, the orthogonal property of wave motion within the porous layers through water depth is derived, and applied in the solution process. The flux and pressure conditions on vertical boundaries and interfaces are integrated to give a set of linear matrix equations, through which the unknown coefficients are solved. Comparisons of the present method with previous studies are preceded in verification, which suggests the validity and practicability of the present study, with a further expectation of extending our work to build a mild-slope equation over multiple-layer porous medium in the future.     

Wave interaction with cylindrical porous piles

The interaction of waves with arrays of porous circular cylinders is studied theoretically and, under the assumption of potential flow and linear wave theory, an analytical solution is derived. The solution is valid for either submerged or emerged structures. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. For specific conditions the model gives the same solution as those previously presented by other authors. Numerical results are presented which exemplify diverse wave and mechanical parameters on the wave transformation due to the presence of a system of circular cylinders.     

Piezoelectric Energy Analysis on Diverse Buoy Coupling with Hydrodynamic Parameters

This paper mainly describes the influence factors of the captured energy power by huge wave energy harvesters, in which the vertical motion of buoy can transform ocean’s potential energy into piezoelectric energy power by undulating waves. Firstly, related environmental coefficients are analyzed by means of the incident wave theory. Besides, the geometric structural parameters are also analyzed and compared under optimal environmental coefficients with semi-analytical solutions. Thirdly, the numerical results also show the impact trend of hydrodynamic parameters and geometric volume on motion, voltage and power with qualitative agreement. The numerical simulation confirms that the improved structure parameters could markedly deliver sufficient power under the same conditions with long-time stability.     

Wave motion over a breakwater system of a horizontal plate and a vertical porous wall

A two-dimensional analytical solution is presented to study the reflection and transmission of linear water waves propagating past a submerged horizontal plate and through a vertical porous wall. The velocity potential in each fluid domain is formulated using three sets of orthogonal eigenfunctions and the unknown coefficients are determined from the matching conditions. Wave elevations and hydrodynamic forces acting on the porous wall are computed. Reflection and transmission coefficients are presented to examine the performance of the breakwater system. The present analytical solutions are found in fairly good agreement with the available laboratory data. The results indicate that the plate length, the porous-effect, the gap between plate and porous wall, and the submerged depth of the plate all show a significant influence on the reflected and transmitted wave fields. It is also interesting to note that the submerged plate plays an important role in reducing the transmitted wave height, especially for long incident waves.     

Interaction of oblique waves with infinite number of perforated caissons

An analytic solution based on the division of the fluid domain is developed for the interaction of obliquely incident waves with infinite number of perforated caissons. The whole fluid domain is firstly divided into infinite sub-domains according to the division of structures, and subsequently eigenfunction expansion is employed to represent the velocity potential in each domain. A phase relation is utilized for the analysis of wave oscillation in each caisson, and the character of structure geometry is considered in setting up the mathematical model of reflection waves. The reflection waves from the present analysis include many propagation waves traveling in different directions when the incident wave frequency is high. Benchmark examinations show that the continuous condition of water particle velocity is satisfied at the front walls of caissons, and the reflection coefficients keep agreement with the energy conservation relation very well when porous effect parameter is infinite. Numerical results show that the reflection coefficients of obliquely incident waves are smaller when the length of caissons is shorter at low frequency. The wave reflection coefficients and the wave forces normal to caissons decrease and the wave forces along caissons increase with the increase of the wave incident angle.     

Numerical Analysis on the Effects of Submerged Depth of the Grid and Direction of Incident Wave on Gravity Cage

In this paper,the numerical model of the net cage with the grid mooring system in waves is set up by the lumped mass method and rigid kinematics theory,and then the motion equations of floating system,net system,mooring system,and floaters are solved by the Runge-Kutta fifth-order method.For the verification of the numerical model,a series of physical model tests have been carried out.According to the comparisons between the simulated and experimental results,it can be found that the simulated and experimental results agree well in each condition.Then,the effects of submerged depth of grid and direction of incident wave propagation on hydrodynamic behaviors of the net cage are analyzed.According to the simulated results,it can be found that with the increase of submerged depth of grid,the forces acting on mooring lines and bridle lines increase,while the forces on grid lines decrease;the horizontal motion amplitudes of floating collar decrease obviously,while the vertical motion amplitudes of floating collar change little.When the direction of incident wave propagation changes,forces on mooting lines and motion of net cage also change accordingly.When the propagation direction of incident wave changes from 0° to 45°,forces on the main ropes and bridle ropes increase,while the forces on the grid ropes decrease.With the increasing propagation direction of incident wave,the horizontal amplitude of the forces collar decreases,while the vertical amplitude of the floating collar has little variation.     

Numerical simulation of non-linear regular and focused waves in an infinite water-depth

Inviscid three-dimensional free surface wave motions are simulated using a novel quadratic higher order boundary element model (HOBEM) based on potential theory for irrotational, incompressible fluid flow in an infinite water-depth. The free surface boundary conditions are fully non-linear. Based on the use of images, a channel Green function is developed and applied to the present model so that two lateral surfaces of an infinite-depth wave tank can be excluded from the calculation domain. In order to generate incident waves and dissipate outgoing waves, a non-reflective wave generator, composed of a series of vertically aligned point sources in the computational domain, is used in conjunction with upstream and downstream damping layers. Numerical experiments are carried out, with linear and fully non-linear, regular and focused waves. It can be seen from the results that the present approach is effective in generating a specified wave profile in an infinite water-depth without reflection at the open boundaries, and fully non-linear numerical simulations compare well with theoretical solutions. The present numerical technique is aimed at efficient modelling of the non-linear wave interactions with ocean structures in deep water.     

Decomposing damped incident and reflected waves using correlation and quasi-linearization methods

Water waves propagating over a layer of soft mud or submerged aquatic vegetation can drastically attenuate over distances comparable to several wave lengths. The attenuation in the case of mud has been found previously to be reasonably described by an exponential decay. Waves reflect from beaches and any structures that they impact. The reflected waves affect wave heights measured in the field or laboratory wave basins.Decomposition of small amplitude waves into incident and reflected waves is a linear problem. However, the presence of the exponential damping introduces nonlinearity to the decomposition problem and requires an iterative process for solving the problem. Despite considerable experimental research on attenuation of waves over mud, none of the existing methods for decomposition of incident and reflected waves have accounted for this case.Here, the Newton Algorithm was used to account for the effect of wave decay over mud by quasi-linearizing the nonlinear equations. Also, a second method using a new error function and a commercial nonlinear solver was proposed in both time and frequency domain. The performance of both methods has been verified using artificial as well as laboratory data.     

Long waves dissipation and harmonic generation by coastal vegetation

In this paper, we study the harmonic generation and energy dissipation as water waves propagating through coastal vegetation. Applying the homogenization theory, linear wave models have been developed for a heterogeneous coastal forest in previous works (e.g. [17], [10], [11]). In this study, the weakly nonlinear effects are investigated. The coastal forest is modeled by an array of rigid and vertically surface-piercing cylinders. Assuming monochromatic waves with weak nonlinearity incident upon the forest, higher harmonic waves are expected to be generated and radiated into open water. Using the multi-scale perturbation theory, micro-scale flows in the vicinity of cylinders and macro-scale wave dynamics are separated. Expressing the unknown variables (e.g. velocity, free surface elevation) as a superposition of different harmonic components, the governing equations for each mode are derived while different harmonics are interacting with each other because of nonlinearity in the cell problem. Different from the linear models, the leading-order cell problem for micro-scale flow motion, driven by the macro-scale pressure gradient, is now a nonlinear boundary-value-problem, while the wavelength-scale problem for wave dynamics remains linear. A modified pressure correction method is employed to solve the nonlinear cell problem. An iterative scheme is introduced to connect the micro-scale and macro-scale problems. To demonstrate the theoretical results, we consider incident waves scattered by a homogeneous forest belt in a constant shallow depth. Higher harmonic waves are generated within the cylinder array and radiated out to the open water region. The comparisons of numerical results obtained by linear and nonlinear models are presented and the behavior of different harmonic components is discussed. The effects of different physical parameters on wave solutions are discussed as well.