Numerical simulation of solitary wave scattering by a circular cylinder array

The interaction of a solitary wave with an array of surface-piercing vertical circular cylinders is investigated numerically. The wave motion is modeled by a set of generalized Boussinesq equations. The governing equations are discretized using a finite element method. The numerical model is validated against the experimental data of solitary wave reflection from a vertical wall and solitary wave scattering by a vertical circular cylinder respectively. The predicted wave surface elevation and the wave forces on the cylinder agree well with the experimental data. The numerical model is then employed to study solitary wave scattering by arrays of two circular cylinders and four circular cylinders respectively. The effect of wave direction on the wave forces and the wave runup on the cylinders is quantified.     

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.     

Wave diffraction and near-trapping by a multi-column gravity-based structure

This paper is concerned with the diffraction of water waves by offshore structures, with the ultimate aim of proposing tools for guiding airgap design. The diffraction of monochromatic waves by an array of four bottom mounted cylinders and a gravity-based structure is studied in detail using linear and second order theory. The phenomenon of near-trapping is investigated, allowing guidelines for airgap design to be established. When contemplating airgap design, however, it is crucially important that consideration is given to the largest waves in a sea state. Therefore, in this study a design wave, called NewWave, is proposed as a realistic model for large ocean waves and is used as the incident wave field in the wave-structure diffraction analysis.     

Second-order resonance among an array of two rows of vertical circular cylinders

Experimental and theoretical works are conducted on the first-order and second-order wave interactions among an array of two rows of vertical circular cylinders.We show that both the first-order and the second-order resonance of dynamic free-surface displacement could in some cases be quite relevant among an array of two rows of vertical cylinders not only in theory but also in real phenomena. In some cases, on the other hand, we show that large free-surface displacement predicted theoretically is quite attenuated in reality probably due to viscous dissipation.     

Numerical implementation of the harmonic modified mild-slope equation

A numerical solver is presented of the modified time-independent mild-slope equation, which incorporates energy dissipation. Using a second-order parabolic approximation, the following external boundary conditions are modelled: open and fully transmitting to both incoming and outgoing waves; partially reflecting, and; fully absorbing. Discretisation of the governing equation and boundary conditions is by means of a second-order accurate central difference scheme. The resulting sparse-banded matrix is solved using an inexpensive banded solver with Gaussian elimination. The numerical predictions are in excellent agreement with the analytical solution for the interaction of non-breaking waves with an array of vertical surface-piercing circular cylinders on a horizontal bed. Results are compared with those for the same array on various seabed topographies. The model is robust and can be used for wave propagation in complex geometries. It has fewer restrictions associated with wave obliqueness at boundaries than traditional models based on the mild-slope equation.     

Solitary wave interactions with an array of two vertical cylinders

This paper addresses a numerical investigation of nonlinear waves interactions with an array of two surface-piercing vertical cylinders and the corresponding nonlinear hydrodynamic loads on each individual cylinder. The primary interest of this study is concentrated on the problem of three-dimensional scattering of solitary waves by cylinder arrays and the nonlinear interactions between scattered waves. The theoretical model adopted for simulation is the generalized Boussinesq two-equation model. The boundary-fitted coordinate transformation and multiple-grid technique are utilized here to simplify the computation domain and to facilitate the applications of the boundary conditions on the cylinder surfaces. The velocity potential, free-surface elevation and subsequent evolution of the scattered wave field are numerically evaluated. The hydrodynamic forces on each cylinder during wave impact are also determined. A study of the sheltering effect by the neighboring structures on wave loads is conducted. It is found that the presence of the neighboring cylinder has shown significant influence on the wave loads and the scattering of the primary incident waves. For two transversely arranged cylinders, the transverse force coefficient increases as the separation distance decreases.     

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.     

正交直墙前直立圆柱的波浪绕射的解析研究

应用镜像原理,将正交直墙前单个圆柱的波浪绕射问题,转化为4个对称布置圆柱的4向入射波浪的绕射问题。然后应用速度势的特征展开方法和Graf加法定理,建立了正交直墙前垂直圆柱的波浪绕射解析解。通过数值计算研究了圆柱与正交直墙间距离大小、波浪入射角等因素与圆柱上总波浪作用力的解析关系。     

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

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

A comparison of complete and approximate solutions for second-order diffraction loads on arrays of vertical ci cular cylinders

This paper presents a comparison between two theoretical methods for computing the second-order diffraction loads on arrays of bottom-mounted, surface-piercing vertical circular cylinders in regular waves. One method presents a complete solution for the second-order hydrodynamic loads on the cylinder array via a numerical integration over the mean fluid free-surface. The other method is based on a large spacing approximation between the array members and involves the solution of a set of equivalent isolated body problems to obtain estimates for the second-order hydrodynamic loads. Numerical results for a pair of cylinders indicate very good agreement between the two methods at center-to-center spacing of both three and five radii, indicating that the approximate method may be sufficient to compute hydrodynamic interference effects to the second-order in many practical engineering situations.     

Analytical Study on Wave Diffraction from a Vertical Circular Cylinder in Front of Orthogonal Vertical Walls

In this paper, the principle of mirror image is used to transform the problem of wave diffraction from a circular cylinder in front of orthogonal vertical walls into the problem of diffraction of four symmetric incident waves from four symmetrically arranged circular cylinders, and then the eigenfunction expansion of velocity potential and Grafs addition theorem are used to give the analytical solution to the wave diffraction problem. The relation of the total wave force on cylinder to the distance between the cylinder and orthogonal vertical walls and the incidence angle of wave is also studied by numerical computation.     

Approximation of near-resonant wave motion using a damped harmonic oscillator model

In this work, transient resonant motions excited in linearised scattering interactions are approximated using a simple damped harmonic oscillator model. The scattering interactions considered involve the diffraction of an incident wave-train with a dominant regular time-harmonic component by structures which enclose a portion of the free surface. Provided a single resonant mode only is excited significantly, the fluid oscillation in the vicinity of the structure will primarily be composed of the resonant mode and incident wave mode contributions. The forced damped harmonic oscillator equation is used to predict the fluid motion and, in particular, the elevation of the internal free surface. The predictions are compared to the results from a numerical time-domain solver based on the linearised water-wave equations. It is shown that, given a good estimate of the location of the resonance in the complex frequency plane and a priori knowledge of the form of the incident wave, the model can successfully predict the time-dependent behaviour of the motion. Both two- and three-dimensional scattering problems are considered with a variety of scatterers in each case.     

Interaction between waves and an array of floating porous circular cylinders

The present study theoretically as well as experimentally investigates the interaction between waves and an array of porous circular cylinders with or without an inner porous plate based on the linear wave theory.To design more effective floating breakwaters,the transmission rate of waves propagating through the array is evaluated.Each cylinder in the array is partly made of porous materials.Specifically,it possesses a porous sidewall and an impermeable bottom.In addition,an inner porous plate is horizontally fixed inside the cylinders.It dissipates the wave more effectively and eliminates the sloshing phenomenon.The approach suggested by Kagemoto and Yue(1986) is adopted to solve the multiple-scatter problem,while a hierarchical interaction theory is adopted to deal with hydrodynamic interactions among a great number of bodies,which efficiently saves computation time.Meanwhile,a series of model tests with an array of porous cylinders is performed in a wave basin to validate the theoretical work and the calculated results.The draft of the cylinders,the location of the inner porous plate,and the spacing between adjacent cylinders are also adjusted to investigate their effects on wave dissipation.     

大尺度海工圆柱的非线性波浪载荷

引入Stokes二阶水波理论,建立了对海洋工程中大尺度综合型圆柱式结构波浪载荷的二阶修正统一积分算式。通过对波势特征解法的推广及引入干扰波叠加法、最小二乘法,使原仅适合于轴对称单柱的二阶散射波浪载荷快速积分算法可完全适用于类型广泛的任意型实用圆柱的计算。几则实例(圆台,单、双柱)的计算检验了该方法的有效性。     

Modelling linear wave transformation induced by dissipative structures—Random waves

The relevant theory is presented and numerical results are compared with the analytical solution for the interaction of non-breaking waves with an array of vertical porous circular cylinders on a horizontal bed. The extension to the cases of unidirectional and multidirectional waves is obtained by means of a transfer function. The influence of the mechanical properties of porous structures and wave irregularity on wave transformation is analysed. Results for unidirectional and multidirectional wave spectra are compared to those obtained for regular waves. The model presented reproduces well the analytical results and provides a tool for analysing several engineering problems.     

不规则波作用下大尺度方形群墩结构的相互影响

本文用谱分析方法对方柱群墩结构在不规则波作用下的相互影响进行了计算。研究了在不同入射波向时每个方墩的群墩影响系数随墩距变化的规律并与折算圆墩的结果进行了比较。入射波谱采用Ｂｒｅｔｓｃｈｎｅｉｄｅｒ谱。计算结果表明方柱墩群与折算圆柱墩群的群墩影响系数较为接近且随着墩间距的增加趋近于孤立墩情形。     

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.     

Benchmark computations of wave run-up on single cylinder and four cylinders by naoe-FOAM-SJTU solver

The benchmark simulations of wave run-up on a fixed single truncated circular cylinder and four circular cylinders are presented in this paper. Our in-house CFD solver naoe-FOAM-SJTU is adopted which is an unsteady two-phase CFD code based on the open source package OpenFOAM. The Navier-Stokes equations are employed as the governing equations, and the volume of fluid (VOF) method is applied for capturing the free surface. Monochromatic incident waves with the specified wave period and wave height are simulated and wave run-up heights around the cylinder are computed and recorded with numerical virtual wave probes. The relationship between the wave run-up heights and the incident wave parameters are analyzed. The numerical results indicate that the presented naoe-FOAM-SJTU solver can provide accurate predictions for the wave run-up on one fixed cylinder and four cylinders, which has been proved by the comparison of simulated results with experimental data.     

Difference-frequency wave loads on a large body in multi-directional waves

The second-order difference-frequency wave forces on a large three-dimensional body in multi-directional waves are computed by the boundary integral equation method and the so-called FML formulation (assisting radiation potential method). Semi-analytic solutions for a bottom-mounted vertical circular cylinder are also developed to validate the numerical method. Difference-frequency wave loads on a bottom-mounted vertical cylinder and stationary four legs of the ISSC tension-leg platform (TLP) are presented for various combinations of incident wave frequencies and headings. These force quadratic transfer functions (QTF) can directly be used in studying slowly varying wave loads in irregular short-crested seas described by a particular directional spectrum. From our numerical results, it is seen that the slowly varying wave loads are in general very sensitive to the directional spreading function of the sea, and therefore wave directionality needs to be taken into account in relevant ocean engineering applications. It is also pointed out that the uni-directionality of the sea is not necessarily a conservative assumption when the second-order effects are concerned.