# An Analytical Solution for the Interaction of Waves with Arrays of Circular Cylinders

Authors: Zhao Mi; Long Pengzhen; Wang Piguang; Zhang Chao; Du Xiuli
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

This paper presents an analytical method to investigate the multiple scattering problem within arrays of vertical bottom-mounted circular cylinders subjected to linear incident waves. Based on the Laplace equation and boundary conditions on the seabed and surface, a formulation of a two-dimensional multiple scattering problem is first obtained by using the variable separation method. Furthermore, the analytical solution of the wave forces on multiple circular cylinders is derived, which consists of the incident wave force due to the linear incident wave and the scattered wave forces considering multiple scattering waves. The presented analytical solution is validated by comparing its results with a numerical method, and the result shows that the analytical solution is in good agreement with the numerical one. Finally, the multiple scattering analysis is conducted on arrays of cylinders with different incident wave numbers, distances between cylinders, and quantities.


---

## Body

## 1. Introduction

Many sea-crossing bridges and offshore structures have been built worldwide in recent years. Pile group foundations play an important role in these offshore structures, and ocean wave is one of the major loadings that threaten the safety of offshore structures. Therefore, it is necessary to determine the wave forces on offshore structures when designing and constructing those structures. Linear diffraction theory is often used to analyze the interaction of ocean waves with cylinders based on the potential theory. Based on this theory, an analytical solution for the linear wave diffraction by a bottom-fixed vertical circular cylinder in finite water depth was proposed by MacCamy and Fuchs [1]. This paper is to extend the analytical solution for a single circular cylinder to circular cylinder arrays.Much attention has been paid to the wave forces on a single vertical cylinder. Morison et al. proposed the Morison equation to estimate the wave forces on slender cylinders, in which the force is composed of a drag force proportional to the square of the velocity and a virtual mass force proportional to the horizontal component of the accelerative force exerted on the mass of water displaced by the cylinder [2]. Except for a large circular cylinder, the diffraction wave theory was also used by researchers to obtain analytical solutions for wave-structure interaction problem of cylinders with various geometric shapes. Williams presented two different approximation methods to predict the wave forces on elliptical cylinders [3]. The first method is based on the exact expressions obtained by Chen and Mei [4], and the second is the integral equation method. Bhatta et al. presented an analytical approach to predict the wave forces due to scattering and radiation on a floating circular cylinder in the water of finite depth [5]. The inquiry into the interaction between a fixed vertical elliptical cylinder and short-crested waves was analyzed by Wang et al. [6]. Liu et al. presented analytical solutions to estimate the diffraction wave forces on a uniform vertical cylinder and a truncated cylinder whose cross sections are arbitrarily smooth [7, 8]. Zhai et al. investigated the diffraction of waves from a system that is made up of a cylinder and an arc-shaped wall and proposed an analytical solution using the eigenfunction expansion approach [9].Except for a single cylinder, the wave-structure interaction of arrays of vertical structures subjected to incident waves has drawn many scholars’ attention because of the extensive applications of cylinder arrays in ocean engineering. Linton and Evans proposed an analytical solution, as an extension of the work of MacCamy and Fuchs, for linear wave diffraction on arrays of vertical circular cylinders [10]. The simplified approach was extended to the complementary radiation problem by Kim [11] and for application to incident focused wave groups on an array of circular cylinders by Walker and Taylor [12]. Other interactions between vertical cylinders and incident waves have been studied, such as a group of dual-cylinder systems with a thin and porous outer cylinder and an impermeable inner cylinder subjected to linear waves [13], arrays of elliptical cylinders wherein wave diffraction happened [14], an array of truncated cylinders with arbitrary cross section subjected to small steepness harmonic wave [15], and collections of vertical bottom-mounted cylinders with arbitrary cross section subjected to linear waves [16].A few numerical schemes also have been proposed to investigate the wave-structure interaction. Ali et al. developed a finite difference numerical method to estimate wave forces on double vertical cylinders [17]. Wang et al. presented finite element solutions to compute the wave forces on arrays of vertical circular cylinders and circular inclined cylinders [18, 19]. Besides the finite element method (FEM), the scaled boundary finite element method (SBFEM) has been recently applied in studies on vertical cylinders with different cross-sectional shapes [20–22]. The SBFEM discretizes only the common interfaces of the subdomains with surface finite elements, and fully accurate results can be obtained by using only fewer elements.In most analytical methods, the whole diffraction analysis mainly focused on the linear incident wave on cylinders and the scattered wave generated by every cylinder, say the first-order scattered wave [10]. Nevertheless, upon impinging on each other, the first-order scattered wave is scattered by the cylinder, and then the second scattered wave is generated. Higher-order scattered wave is also generated in this manner. The term order denotes the scattering times that an incident wave has encountered. Therefore, this study aims to propose an analytical solution to investigate how multiple scattering waves influence the wave-cylinder interaction.


## 2. Mathematical Formulation

There is situated a group ofN (N ≥ 2) vertical circular cylinders in water with depth h, as shown in Figure 1. Water is assumed to be irrotational, inviscid, and incompressible. The origin of the global Cartesian coordinate system (x, y, z) is located on the seabed level, and the z-axis is directed vertically upwards. The center of the ith cylinder at (xi, yi) with radius ai is taken as the origin of a local polar coordinate system (ri, θi), where θi is measured counterclockwise from the positive x-axis.Figure 1
The interaction of water with arrays of circular cylinders.Figure2 shows the relation between a global coordinate system and a local coordinate system of the ith cylinder, where the coordinates of Q in global Cartesian coordinate system satisfy the following equations:(1)x=Xi+xi,y=Yi+yi,where Xi and Yi are the coordinates of Oi of the local coordinate system of ith cylinder and (xi, yi) are the coordinates of Q in the local coordinate system of ith cylinder.Figure 2
Coordinate systems and geometry for Graf’s addition theorem adapted for multiple scattering problem.Assuming that the direction of wave propagation of an angleβ to the positive x-axis and cylinders are exposed to the linear incident waves, the total wave pressure in the time domain can be written as pe=pe−iωt, with i=−1, time t and the wave frequency ω, and the total wave pressure in the frequency domain p. The governing equation and boundary conditions are expressed in the polar coordinate system as(2)∂2p∂r2+1r∂p∂r+1r2∂2p∂θ2+∂2p∂z2=0.(1)
Seabed boundary condition (z = 0):(3)∂p∂z=0.(2)
The free surface condition (z = h):(4)−ω2gp+∂p∂z=0.(3)
The Sommerfeld radiation conditions at infinity:(5)limr⟶∞r∂p∂r−ikp=0,where k satisfies the dispersion relation, which is expressed as ω2=gktanhkh.The total pressurep in the fluid domain consists of the incident wave pressure pI and scattered wave pressure pS, which satisfy the Laplace equation.(6)∇2p=∇2pI=∇2pS=0.


## 3. Analytical Solution

In this part, as an isolated cylinder is exposed to the monochromatic incident waves, the interaction between it and the surrounding water will be reviewed first. Then, there will be an investigation on the multiple scattering problem among an array of circular cylinders, assuming that the influences on wave force of a column caused by other columns could be superposed approximately and the influence of columnA on column B is immutable from other column.


### 3.1. Incident Wave Pressure

According to equation (6) and boundary conditions in equations (3) and (4), the incident wave pressure and the vertical mode of the fluid can be expressed as [23](7)pI=Zzeikxx+kyy,(8)Zz=ρgHcoshkz2coshkh,in which i=−1, ρ is the fluid density, H is the wave height, and k is the total wave number with kx=kcosβ and ky=ksinβ.The form ofpiI in the corresponding local polar coordinate system of the ith cylinder (ri, θi), also denoted by piI,0, can be expressed as(9)piI,0=ρgHcoshkz2coshkheiLikcosψi−β∑n=0∞βnJnkricosnθi−β,where βn=1,n=0,2in,n≥1 and Jn is the Bessel function of the first kind of order n.


### 3.2. First-Order Scattered Wave Pressure due to Incident Waves

The scattered wave pressure on theith cylinder due to the linear incident waves is denoted by piS,1 which is called the first-order scattered wave pressure. According to equation (7), after applying the variable separation method, piS,1 can be expressed as(10)piS,1x,y,z=PiS,1x,yZz.Thus, the cylinder surface boundary condition forPiS,1 can be obtained as(11)∂PiS,1∂ri|ri=ai+∂PiI,0∂ri|ri=ai=0.Substituting equations (8) and (10) into equation (6), a two-dimensional Helmholtz equation in polar coordinate system is obtained as(12)∂2PiS,1∂ri2+1ri∂PiS,1∂ri+1ri2∂2PiS,1∂θi2+k2PiS,1=0.Applying the variable separation method again,PiS,1 can be expressed as(13)PiS,1ri,θi=RriΘθi.Substituting equation (13) into equation (12) yields the following equations:(14)Θ″+n2Θ=0,(15)ri2R″+riR′+k2ri2−n2R=0.The solution of equation (14) can be written as(16)Θnθi=Ai,nS,1cosnθi+Bi,nS,1sinnθi,n=0,1,2,….According to the Sommerfeld radiation condition equation (5), the solution of equation (15) which is Hankel function is(17)Rnri=Ci,nS,1Hnkri,where Hn is the Hankel function of the first kind of order n.Hence, the solution of the scattered wave pressurePiS,1 can be expressed as(18)PiS,1=∑n=0∞HnkriEi,nS,1cosnθi+Fi,nS,1sinnθi.Substituting equations (18) and (9) into equation (11) and using the orthogonality of mode functions cosnθi or sinnθi, PiS,1 can be obtained as(19)PiS,1=eikLicosψi−β∑n=0∞−Jn′kaicosnβHn′kaiβnHnkricosnθi+∑n=0∞−Jn′kaisinnβHn′kaiβnHnkrisinnθi.Figure2 also shows geometric relations between two local polar coordinate systems. The first-order scattered wave pressure of other cylinder can be expressed in equation (19), but it is in its own local polar coordinate system. Thus, in order to calculate pressures in one coordinate, the Graf’s addition theorem for Bessel functions can be recalled as [24](20)Hnkrlcosnθl=∑v=−∞∞Hn+vkRilJvkricosvθi−nπ−v+n¯θil,Hnkrlsinnθl=−∑v=−∞∞Hn+vkRilJvkrisinvθi−nπ−v+n¯θil,where Ril=xi−xl2+yi−yl2.Thus, with first-order scattered wave from thelth cylinder as the excitation source, the incoming wave impinging ith cylinder denoted by Pidue toPlS,1I,1 can be expressed as(21)PiduetoPlS,1I,1=∑n=0∞βnJn′kalHn′kal−cosnβ∑v=−∞∞Hn+vkRilJvkricosvθi−nπ−v+n¯θil+sinnβ∑v=−∞∞Hn+vkRilJvkrisinvθi−nπ−v+n¯θileikLlcosψl−β.Finally, the total incoming wave pressure on theith cylinder due to the first-order scattered wave from other cylinders denoted by PiI,1 can be expressed as(22)PiI,1=∑l=1l≠iNPiduetoPlS,1I,1.


### 3.3. Higher-Order Scattered Wave Pressure

Higher-order scattered waves are generated in the same manner as mentioned in Section3.2. For a general qth order scattered wave pressure on the ith cylinder denoted by PiS,q (q > 1), it can be obtained as(23)PiS,q=∑n=0∞HnkriEi,nS,qcosnθi+Fi,nS,qsinnθi,where the coefficients Ei,nS,q and Fi,nS,q can be obtained as(24)Ei,nS,q=−1δπkHn′kai∫02π∂PiI,q−1∂ricosnθi|ri=aidθi,Fi,nS,q=−1δπkHn′kai∫02π∂PiI,q−1∂risinnθi|ri=aidθi,where δ=2 if n=0 or δ=1 if n≠0, and the superscript prime denotes the derivative to the Hankel function.Withq−1th order scattered wave from the lth cylinder as the excitation source, the incoming wave impinging ith cylinder denoted by PiduetoPlS,q−1I,q−1 can be expressed as(25)PiduetoPlS,q−1I,q−1=∑n=0∞∑v=−∞∞Hn+vkRilJvkriEl,nS,q−1cosvθi−nπ−v+n¯θil−Fl,nS,q−1sinvθi−nπ−v+n¯θil.The total incoming wave pressures on theith cylinder due to the q−1th order scattered wave from other cylinders denoted by PiI,q−1 can be expressed as(26)PiI,q−1=∑l=1l≠iNPiduetoPlS,q−1I,q−1.Substituting equation (26) into equation (24) and using the orthogonality, Ei,nS,q and Fi,nS,q can be simplified as(27)Ei,nS,q=−1δHn′kai∑l=1l≠iN∑j=0∞Jn′kaiHj′kalEl,jS,q−1Hj+nkRilcosjπ+j+n¯θil+−1nHj−nkRilcosjπ+j−n¯θil+Fl,jS,q−1Hj+nkRilsinjπ+j+n¯θil+−1nHj−nkRilsinjπ+j−n¯θil,Fi,nS,q=−1δHn′kai∑l=1l≠iN∑j=0∞Jn′kaiHj′kalEl,jS,q−1Hj+nkRilsinjπ+j+n¯θil−−1nHj−nkRilsinjπ+j−n¯θil−Fl,jS,q−1Hj+nkRilcosjπ+j+n¯θil−−1nHj−nkRilcosjπ+j−n¯θil.


### 3.4. Wave Forces on the Cylinders

The wave pressure on theith cylinder can be written as(28)Pi=PiI,0+PiS,1+∑q=2∞PiS,q+PiI,q−1.Furthermore, the resultant circumferential wave force alongx-axis and y-axis on the ith cylinder at height z is denoted by fixze−iωt and fiyze−iωt, where fixz and fiyz are given by(29a)fixz=FixρgHcoshkz2coshkh,(29b)fiyz=FiyρgHcoshkz2coshkh.According to the equation (28), Fix and Fiy can be expressed as(30a)Fix=−∫02πPiaicosθidθi,(30b)Fiy=−∫02πPiaisinθidθi.Finally, the total wave force on theith cylinder along the x-axis and y-axis can be calculated by(31a)Fi,xe=∫0hfixze−iωtdz,(31b)Fi,ye=∫0hfiyze−iωtdz.


## 3.1. Incident Wave Pressure

According to equation (6) and boundary conditions in equations (3) and (4), the incident wave pressure and the vertical mode of the fluid can be expressed as [23](7)pI=Zzeikxx+kyy,(8)Zz=ρgHcoshkz2coshkh,in which i=−1, ρ is the fluid density, H is the wave height, and k is the total wave number with kx=kcosβ and ky=ksinβ.The form ofpiI in the corresponding local polar coordinate system of the ith cylinder (ri, θi), also denoted by piI,0, can be expressed as(9)piI,0=ρgHcoshkz2coshkheiLikcosψi−β∑n=0∞βnJnkricosnθi−β,where βn=1,n=0,2in,n≥1 and Jn is the Bessel function of the first kind of order n.


## 3.2. First-Order Scattered Wave Pressure due to Incident Waves

The scattered wave pressure on theith cylinder due to the linear incident waves is denoted by piS,1 which is called the first-order scattered wave pressure. According to equation (7), after applying the variable separation method, piS,1 can be expressed as(10)piS,1x,y,z=PiS,1x,yZz.Thus, the cylinder surface boundary condition forPiS,1 can be obtained as(11)∂PiS,1∂ri|ri=ai+∂PiI,0∂ri|ri=ai=0.Substituting equations (8) and (10) into equation (6), a two-dimensional Helmholtz equation in polar coordinate system is obtained as(12)∂2PiS,1∂ri2+1ri∂PiS,1∂ri+1ri2∂2PiS,1∂θi2+k2PiS,1=0.Applying the variable separation method again,PiS,1 can be expressed as(13)PiS,1ri,θi=RriΘθi.Substituting equation (13) into equation (12) yields the following equations:(14)Θ″+n2Θ=0,(15)ri2R″+riR′+k2ri2−n2R=0.The solution of equation (14) can be written as(16)Θnθi=Ai,nS,1cosnθi+Bi,nS,1sinnθi,n=0,1,2,….According to the Sommerfeld radiation condition equation (5), the solution of equation (15) which is Hankel function is(17)Rnri=Ci,nS,1Hnkri,where Hn is the Hankel function of the first kind of order n.Hence, the solution of the scattered wave pressurePiS,1 can be expressed as(18)PiS,1=∑n=0∞HnkriEi,nS,1cosnθi+Fi,nS,1sinnθi.Substituting equations (18) and (9) into equation (11) and using the orthogonality of mode functions cosnθi or sinnθi, PiS,1 can be obtained as(19)PiS,1=eikLicosψi−β∑n=0∞−Jn′kaicosnβHn′kaiβnHnkricosnθi+∑n=0∞−Jn′kaisinnβHn′kaiβnHnkrisinnθi.Figure2 also shows geometric relations between two local polar coordinate systems. The first-order scattered wave pressure of other cylinder can be expressed in equation (19), but it is in its own local polar coordinate system. Thus, in order to calculate pressures in one coordinate, the Graf’s addition theorem for Bessel functions can be recalled as [24](20)Hnkrlcosnθl=∑v=−∞∞Hn+vkRilJvkricosvθi−nπ−v+n¯θil,Hnkrlsinnθl=−∑v=−∞∞Hn+vkRilJvkrisinvθi−nπ−v+n¯θil,where Ril=xi−xl2+yi−yl2.Thus, with first-order scattered wave from thelth cylinder as the excitation source, the incoming wave impinging ith cylinder denoted by Pidue toPlS,1I,1 can be expressed as(21)PiduetoPlS,1I,1=∑n=0∞βnJn′kalHn′kal−cosnβ∑v=−∞∞Hn+vkRilJvkricosvθi−nπ−v+n¯θil+sinnβ∑v=−∞∞Hn+vkRilJvkrisinvθi−nπ−v+n¯θileikLlcosψl−β.Finally, the total incoming wave pressure on theith cylinder due to the first-order scattered wave from other cylinders denoted by PiI,1 can be expressed as(22)PiI,1=∑l=1l≠iNPiduetoPlS,1I,1.


## 3.3. Higher-Order Scattered Wave Pressure

Higher-order scattered waves are generated in the same manner as mentioned in Section3.2. For a general qth order scattered wave pressure on the ith cylinder denoted by PiS,q (q > 1), it can be obtained as(23)PiS,q=∑n=0∞HnkriEi,nS,qcosnθi+Fi,nS,qsinnθi,where the coefficients Ei,nS,q and Fi,nS,q can be obtained as(24)Ei,nS,q=−1δπkHn′kai∫02π∂PiI,q−1∂ricosnθi|ri=aidθi,Fi,nS,q=−1δπkHn′kai∫02π∂PiI,q−1∂risinnθi|ri=aidθi,where δ=2 if n=0 or δ=1 if n≠0, and the superscript prime denotes the derivative to the Hankel function.Withq−1th order scattered wave from the lth cylinder as the excitation source, the incoming wave impinging ith cylinder denoted by PiduetoPlS,q−1I,q−1 can be expressed as(25)PiduetoPlS,q−1I,q−1=∑n=0∞∑v=−∞∞Hn+vkRilJvkriEl,nS,q−1cosvθi−nπ−v+n¯θil−Fl,nS,q−1sinvθi−nπ−v+n¯θil.The total incoming wave pressures on theith cylinder due to the q−1th order scattered wave from other cylinders denoted by PiI,q−1 can be expressed as(26)PiI,q−1=∑l=1l≠iNPiduetoPlS,q−1I,q−1.Substituting equation (26) into equation (24) and using the orthogonality, Ei,nS,q and Fi,nS,q can be simplified as(27)Ei,nS,q=−1δHn′kai∑l=1l≠iN∑j=0∞Jn′kaiHj′kalEl,jS,q−1Hj+nkRilcosjπ+j+n¯θil+−1nHj−nkRilcosjπ+j−n¯θil+Fl,jS,q−1Hj+nkRilsinjπ+j+n¯θil+−1nHj−nkRilsinjπ+j−n¯θil,Fi,nS,q=−1δHn′kai∑l=1l≠iN∑j=0∞Jn′kaiHj′kalEl,jS,q−1Hj+nkRilsinjπ+j+n¯θil−−1nHj−nkRilsinjπ+j−n¯θil−Fl,jS,q−1Hj+nkRilcosjπ+j+n¯θil−−1nHj−nkRilcosjπ+j−n¯θil.


## 3.4. Wave Forces on the Cylinders

The wave pressure on theith cylinder can be written as(28)Pi=PiI,0+PiS,1+∑q=2∞PiS,q+PiI,q−1.Furthermore, the resultant circumferential wave force alongx-axis and y-axis on the ith cylinder at height z is denoted by fixze−iωt and fiyze−iωt, where fixz and fiyz are given by(29a)fixz=FixρgHcoshkz2coshkh,(29b)fiyz=FiyρgHcoshkz2coshkh.According to the equation (28), Fix and Fiy can be expressed as(30a)Fix=−∫02πPiaicosθidθi,(30b)Fiy=−∫02πPiaisinθidθi.Finally, the total wave force on theith cylinder along the x-axis and y-axis can be calculated by(31a)Fi,xe=∫0hfixze−iωtdz,(31b)Fi,ye=∫0hfiyze−iωtdz.


## 4. Numerical Results

As shown in Figure3, four different configurations, two and three circular cylinders in tandem (N = 2, 3), four cylinders in square (N = 4), and five cylinders in cross star (N = 5), are used in this study to investigate the effect of higher-order scattered waves on wave forces, where D denotes the distance between two adjacent cylinders. A dimensionless parameter is defined as(32)Dr=D2a,where D denotes the distance between two adjacent cylinders and a denotes the radius of the cylinder.Figure 3
The sketch of pile group withN piles.


### 4.1. Validation

The present analytical solution consists of an infinite series including all the scattering waves in any order, which shall be truncated into a finite term in computation since the calculation speed is limited by the number of the term. As shown in Figures4 and 5, the results will converge when truncating the infinite series into 15–18 terms. Therefore, the truncation order is selected as q = 15 in the following analysis.Figure 4
The wave pressures on the first cylinderP1 versus θ for different q with β=0° and N = 2.Figure 5
The wave forcesFix versus ka for different q with β=0°. (a) N = 2. (b) N = 5.
(a)(b)The case withN = 2 is firstly used to verify the effectiveness of the proposed analytical solution. Figure 6 shows the wave force Fix versus ka on two cylinders in tandem (N = 2) with β=0°, β=45°, and β=90° of the present method and FEM [18]. The cases with N = 3, 4, and 5 are further used to verify the effectiveness of the proposed analytical solution. Figures 7 and 8 show the wave field (P) obtained by the present method and FEM, respectively. The relative errors of wave pressure on each cylinder between Figures 7 and 8 are depicted in Figure 9. It can be concluded from Figures 6–9 that the present method agrees well with the FEM.Figure 6
The wave forces on cylindersFix versus ka with N = 2 obtained by the present method and the FEM [6]. (a) β=0°. (b) β=45°. (c) β=90°.
(a)(b)(c)Figure 7
The wave field obtained by the present method withka = 1, Dr = 0.5, and a = 5 m.Figure 8
The wave field obtained by the FEM withka = 1, Dr = 0.5, and a = 5 m.Figure 9
The relative error of the wave pressure between the present method and the FEM withka = 1, Dr = 0.5, a = 5 m, and N = 3, 4, 5.


### 4.2. Effect of Higher-Order Scattered Waves

Two dimensionless parametersIix=Fi,q=15x/Fi,q=2x and Iiy=Fi,q=15y/Fi,q=2y are defined to investigate the effect of higher-order scattered waves on the wave forces on the cylinders, where Fi,q=2x and Fi,q=15x denote the total wave forces on the ith cylinder along the x-axis with q = 2 and q = 15, respectively; Fi,q=2y and Fi,q=15y denote the total wave forces on the ith cylinder along the y-axis.The case with three cylinders arranged in tandem is investigated firstly. Figure10 shows the ratios Iix and Iiy versus ka with Dr=1 for different β. Figure 11 shows the ratios Iix and Iiy versus Dr with β=45° for different ka. The effect of higher-order scattered waves on the wave forces on each cylinder along the x-axis decreases as the wave propagation angle changes from β=0° to β=60°. However, higher-order scattered waves have a limited impact on the wave forces along the y-axis.Figure 10
RatiosIix and Iiy versus ka for different β with N = 3. (a) Left pile. (b) Middle pile. (c) Right pile.
(a)(b)(c)Figure 11
RatiosI2x and I2y of the middle pile versus Dr for different ka with N = 3. (a) I2x. (b) I2y.
(a)(b)The case with four cylinders arranged in square and five cylinders arranged in cross star are investigated secondly. Note that the wave force onC3 and C4 of the four cylinders and on the C5 of the five cylinders is omitted due to the symmetry and wave propagation angle β=0°. It can be observed from Figures 12 and 13 that the effects of higher-order scattered waves on the wave forces on the cylinders can be ignored when ka is small (≤0.3), while the effects are significant when ka is large (>0.3). This means higher-order scattered waves can have more influence on the wave forces on the cylinders with larger wave number. The results also indicate that the effect of higher-order scattered waves on the wave forces on the cylinders is significant when Dr is short between cylinders and that the influence tends to drop as the relative distance increases.Figure 12
RatioIix versus ka for different Dr with β=0° and N = 4. (a) Front pile. (b) Back pile.
(a)(b)Figure 13
RatioIix versus Dr for different ka with β=0° and N = 5. (a) Left pile. (b) Middle pile. (c) Right pile. (d) Top pile.
(a)(b)(c)(d)Figure14 shows the ratio I1x of the first cylinder in all cases. It illustrates that the effect of higher-order scattered waves on the wave forces on the cylinder tends to increase as N increases.Figure 14
RatioI1x of the first pile versus Dr with β=0° and N = 2, 3, 4, 5. (a) ka=0.2. (b) ka=0.6. (c) ka=1.0.
(a)(b)(c)


## 4.1. Validation

The present analytical solution consists of an infinite series including all the scattering waves in any order, which shall be truncated into a finite term in computation since the calculation speed is limited by the number of the term. As shown in Figures4 and 5, the results will converge when truncating the infinite series into 15–18 terms. Therefore, the truncation order is selected as q = 15 in the following analysis.Figure 4
The wave pressures on the first cylinderP1 versus θ for different q with β=0° and N = 2.Figure 5
The wave forcesFix versus ka for different q with β=0°. (a) N = 2. (b) N = 5.
(a)(b)The case withN = 2 is firstly used to verify the effectiveness of the proposed analytical solution. Figure 6 shows the wave force Fix versus ka on two cylinders in tandem (N = 2) with β=0°, β=45°, and β=90° of the present method and FEM [18]. The cases with N = 3, 4, and 5 are further used to verify the effectiveness of the proposed analytical solution. Figures 7 and 8 show the wave field (P) obtained by the present method and FEM, respectively. The relative errors of wave pressure on each cylinder between Figures 7 and 8 are depicted in Figure 9. It can be concluded from Figures 6–9 that the present method agrees well with the FEM.Figure 6
The wave forces on cylindersFix versus ka with N = 2 obtained by the present method and the FEM [6]. (a) β=0°. (b) β=45°. (c) β=90°.
(a)(b)(c)Figure 7
The wave field obtained by the present method withka = 1, Dr = 0.5, and a = 5 m.Figure 8
The wave field obtained by the FEM withka = 1, Dr = 0.5, and a = 5 m.Figure 9
The relative error of the wave pressure between the present method and the FEM withka = 1, Dr = 0.5, a = 5 m, and N = 3, 4, 5.


## 4.2. Effect of Higher-Order Scattered Waves

Two dimensionless parametersIix=Fi,q=15x/Fi,q=2x and Iiy=Fi,q=15y/Fi,q=2y are defined to investigate the effect of higher-order scattered waves on the wave forces on the cylinders, where Fi,q=2x and Fi,q=15x denote the total wave forces on the ith cylinder along the x-axis with q = 2 and q = 15, respectively; Fi,q=2y and Fi,q=15y denote the total wave forces on the ith cylinder along the y-axis.The case with three cylinders arranged in tandem is investigated firstly. Figure10 shows the ratios Iix and Iiy versus ka with Dr=1 for different β. Figure 11 shows the ratios Iix and Iiy versus Dr with β=45° for different ka. The effect of higher-order scattered waves on the wave forces on each cylinder along the x-axis decreases as the wave propagation angle changes from β=0° to β=60°. However, higher-order scattered waves have a limited impact on the wave forces along the y-axis.Figure 10
RatiosIix and Iiy versus ka for different β with N = 3. (a) Left pile. (b) Middle pile. (c) Right pile.
(a)(b)(c)Figure 11
RatiosI2x and I2y of the middle pile versus Dr for different ka with N = 3. (a) I2x. (b) I2y.
(a)(b)The case with four cylinders arranged in square and five cylinders arranged in cross star are investigated secondly. Note that the wave force onC3 and C4 of the four cylinders and on the C5 of the five cylinders is omitted due to the symmetry and wave propagation angle β=0°. It can be observed from Figures 12 and 13 that the effects of higher-order scattered waves on the wave forces on the cylinders can be ignored when ka is small (≤0.3), while the effects are significant when ka is large (>0.3). This means higher-order scattered waves can have more influence on the wave forces on the cylinders with larger wave number. The results also indicate that the effect of higher-order scattered waves on the wave forces on the cylinders is significant when Dr is short between cylinders and that the influence tends to drop as the relative distance increases.Figure 12
RatioIix versus ka for different Dr with β=0° and N = 4. (a) Front pile. (b) Back pile.
(a)(b)Figure 13
RatioIix versus Dr for different ka with β=0° and N = 5. (a) Left pile. (b) Middle pile. (c) Right pile. (d) Top pile.
(a)(b)(c)(d)Figure14 shows the ratio I1x of the first cylinder in all cases. It illustrates that the effect of higher-order scattered waves on the wave forces on the cylinder tends to increase as N increases.Figure 14
RatioI1x of the first pile versus Dr with β=0° and N = 2, 3, 4, 5. (a) ka=0.2. (b) ka=0.6. (c) ka=1.0.
(a)(b)(c)


## 5. Conclusion

In this paper, an analytical method is presented to investigate wave forces on arrays of cylinders, which include the incident wave forces due to the linear incident waves, first-order scattered wave forces, and higher-order scattered wave forces due to lower-order scattered waves. The validation results indicate that the present analytical solution agrees well with the finite element method. The effect of higher-order scattered waves on the wave forces on the cylinders is investigated by using the proposed method. Based on the results, we can conclude that the effects of higher-order scattered waves on the wave forces on the cylinders tend to increase significantly as the wave number and the number of cylinders increase and as distance between cylinders decreases. When cylinders are arranged in tandem, the effects of higher-order scattered waves on the wave forces on cylinders along thex-axis decrease as the wave propagation angle increases, and the higher-order scattered waves have little influence on the wave forces on cylinders along the y-axis. It can also be concluded that the effects of higher-order scattered waves can be neglected when ka≤0.3.



---

*Source: 5710894-2021-10-11.xml*