Numerical Investigation of Multi-body Wave Energy Converters’ Configuration

Article information

J. Ocean Eng. Technol. 2022;36(2):132-142
Publication date (electronic) : 2022 April 12
doi : https://doi.org/10.26748/KSOE.2022.002
1Postdoctoral Researcher, School of Naval Architecture and Ocean Engineering, University of Ulsan, Ulsan, Korea
2Professor, School of Naval Architecture and Ocean Engineering, University of Ulsan, Ulsan, Korea
Corresponding author Kyeonguk Heo: +82-52-259-2158, hhggoo@ulsan.ac.kr
Received 2022 February 3; Revised 2022 March 2; Accepted 2022 March 14.

Abstract

We investigate the performance of multi-body wave energy converters (WECs). This investigation considers multiple scattering of water waves by the buoys of a WEC under the generalized mode approach. Predominantly, the effect of a WEC’s configuration on its energy extraction is studied in this research. First, single-row terminator and single-column attenuator arrays of vertical cylinders have been studied. The performance of these attenuator arrays shows that the wall effect induced by the periodic buoys influences the wave propagation and energy extraction in these WECs. Further studies show that a single-row terminator array of vertical cylinders performs better than the corresponding single-column attenuator array. Subsequently, multi-row terminator arrays of vertical cylinders are investigated by conducting a parametric study. This parametric study shows that the hydrodynamic property of three resonance phenomena makes energy extraction efficiency drop down, and the magnitude of energy extracted oscillates between the resonance points in these WECs. Finally, a 4×8 terminator array of vertical cylinders is studied to determine the effect of various dx (x-directional distance between adjacent rows) within this WEC on its performance. In particular, this study enforces at least two equal dx values within the 4×8 terminator array of vertical cylinders. It shows that a small value of this dx leads to better energy extraction efficiency in some of these various dx arrays than that of a corresponding regular array with the same dx.

1. Introduction

Renewable energy research is given the utmost importance worldwide to overcome the striking threat to the environment from non-renewable energy sources due to their conversion. In particular, energy generation using the aspects of the oceans has a huge potential as oceans are renewable sources of energy in abundance. Relatedly, wave energy converters (WECs) convert the kinematic energy of water waves into a usable form of energy. So, WECs have been researched under naval architecture and ocean engineering for a long time. WECs are of three types: oscillating body, overtopping, and oscillating water column WECs. In particular, oscillating body WECs use a buoy’s kinematic energy resulting from the water wave-body interaction and convert it into electric energy using a power take-off (PTO) device. Generally, an oscillating body WEC has several buoys. Moreover, the following oscillating body WECs result from the relative orientation between their buoys and water waves (Mei et al., 2005): terminator, attenuator, and omnidirectional oscillating body WECs. For example, the buoys of a terminator type oscillating body WEC are arranged perpendicular to the water wave direction. On the other hand, the buoys of an attenuator type oscillating body WEC are arranged along the water wave direction. Finally, an omnidirectional type oscillating body WEC performs identically for all arrangements of its buoys with respect to the water wave direction.

Several research activities have been carried out worldwide on oscillating body WECs to improve their energy conversion. Generally, it is better to use several smaller buoys instead of a larger single buoy in them (e.g., Garnaud and Mei, 2009). The oscillating body WECs with several buoys are shortly called multi-body WECs. However, the beneficial arrangement of buoys within multi-body WECs is unknown because it is hard to implement the multiple scattering of water waves by the buoys of the WECs in the analysis to find this arrangement.

Several theories have been proposed on the interaction between water waves and the buoys of multi-body WECs. Usually, a point absorber assumption is made in such theories in the long-wave region. In general, this assumption considers the incident and radiated water waves but not the diffracted water waves from the bodies (e.g., Falnes and Budal, 1982). On the other hand, a plane wave approximation considers the multiple scattering of water waves by the buoys in a multi-body WEC (Simon, 1982). However, this approximation cannot be used for an array of buoys in a multi-body WEC provided that the distance between buoys is close. So, an approximation suggested by McIver and Evans (1984) for vertical cylinders by correcting diffracted water wave components can be used in such situations. Notably, the multiple scattering of water waves by the buoys was studied by Kagemoto and Yue (1986) exactly. Their new multi-body hydrodynamic interaction theory used linear algebra of velocity potentials. In particular, they applied a characteristic matrix to the diffraction potential of the bodies by considering a generalized incident-wave potential for multi-bodies. Subsequently, another exact method for analyzing the interaction of water waves with vertical cylinders was proposed by Linton and Evans (1990) based on eigenfunction expansion. Later, increased computational capabilities led to a generalized mode approach to study multi-body interactions using the boundary element method (Newman, 1994). In particular, Newman (1994) demonstrated the applications of the generalized mode approach in studies other than multi-body interactions, such as hydroelasticity, hinge mode analysis, etc., as well. He also established a generalized boundary value problem and a generalized equation of the hydrodynamic forces in this approach. The exact methods of Kagemoto and Yue (1986) and Linton and Evans (1990) and the generalized mode approach of Newman (1994) are considered complete in analyzing the multiple scattering of water waves by the buoys under linear potential theory without special constraints. Meanwhile, these analyses become cumbersome under any numerical method when hundreds of buoys are used. Hence, Kashiwagi (2000) extended the theory of Kagemoto and Yue (1986) using a hierarchical ordering among a huge number of hydrodynamically interacting bodies and developed a new multi-body hydrodynamic interaction theory. Likewise, Tokić and Yue (2019) analyzed the multiple scattering of water waves by a large number of buoys recently. They found the hydrodynamic properties such as Bragg, Rayleigh, and Laue resonances and the existence of Rayleigh-Bloch waves in a periodic array of WECs.

This study is motivated by the research of Tokić and Yue (2019) and analyzes the effect of a multi-body WEC’s configuration on its energy extraction. We use the generalized mode approach and study the multiple scattering of water waves by the buoys (vertical cylinders) of different WECs. In addition, a parametric study has been conducted on the terminator and attenuator arrays of vertical cylinders in the absence of oblique waves.

2. Formulation

2.1 Equation of Motion of Buoys in a WEC

This study uses the generalized mode approach of Choi (2011) to consider the multiple scattering of water waves by the buoys of a multi-body WEC. In ideal fluid assumption, hydrodynamic pressure could be obtained from the velocity potential. The velocity potentials under this multiple scattering can be obtained using a boundary integral equation. However, the surface of each buoy should be discretized, and the normal vector components of the buoys’ surface should be defined to calculate the velocity potential. The schematic of a multi-body WEC used in this analysis is shown in Fig. 1. First, the velocity potential is defined in the following form:

(1) Φ=Re[{A(φI+φS)+ξjφj}e-iωt]
where subscript I means incident wave and S is scattering. A is the wave amplitude, ξj is the j-th mode of a buoy’s motion, ω is the wave frequency, and t is the time. In particular, we constrain all the modes of a buoy’s motion but its heave mode of motion in this study. Further, the body boundary conditions of the velocity potential are as follows:
(2) φSn=-φInonSBkφjn=-iωnjonSBk
where SB is the mean body surface of the k-th buoy, N is the number of buoys, and nj is the generalized directional cosine expressed by the inner product between the surface normal vector of a buoy and the vertical unit vector (Choi, 2011). Hence,
(3) nj={NjkonSBj0otherwise

Fig. 1.

Schematic of a multi-body WEC

The added mass (akj) and the damping coefficient (bkj) of the k-th buoy for the j-th mode of motion can be obtained for the generalized mode by integrating the product of the radiation velocity potential for the j-th mode of motion and the generalized directional cosine of the k-th buoy. Similarly, the unit water wave exciting force (Xk) is calculated by integrating the product of the diffraction velocity potential and the generalized directional cosine of the k-th buoy (Newman, 1994). Exactly,

(4) ω2akj+iωbkj=iωρ∫∫SBkφjnkdS
(5) Xk=iωρ∫∫SBk(φI+φS)nkdS
where ρ is the water density. In addition, the damping force induced due to the extracted energy of a power take-off (PTO) device is assumed to be linear in this analysis. So, the PTO damping force of each buoy is:
(6) FDk=-Re[-iωbPTO,kξke-iωt]
where bPTO is the damping coefficient of the PTO device. Subsequently, the equation of motion of the buoys that considers the PTO damping force is:
(7) j=1N[-ω2(mkj+akj)-iω(bkj+bPTO,jδkj+ckj]ξj=AXk
where δkjis the Kronecker delta function, mkj is the mass of the k-th buoy for the j-th mode of motion, and ckj is the hydrostatic stiffness of the k-th buoy for the j-th mode of motion.

2.2 Energy Extraction of the WEC

Generally, the time-mean extracted power of the WEC from the waves, , is written in the following form.

(8) P¯=-k=1NFDkξ˙k¯=ω22k=1NbPTO,k|ξk|2

In addition, the calculation of energy extraction efficiency of the WEC needs the capture width defined in the following form:

(9) W=P¯12ρg|A|2Cg=ω2ρgCgk=1NbPTO,k|ξkA|2
where Cg is the wave group velocity and g is the acceleration due to gravity. The capture width can be non-dimensionalized by multiplying it with the wavenumber (k). Notably, the maximum non-dimensional capture width is one at the natural frequency of the buoys when considering only one mode of motion of the buoys. Hence,
(10) (kW)max=1forbPTO(ωn)=b(ωn)

In addition, the efficiency of a multi-body WEC’s configuration can be found using the concept of array gain (Tokić and Yue, 2019) as follows:

(11) q(k,θI,C)=Pa¯(k,θI,C)NPiso¯
where Pa¯ and Piso¯ denote the time-mean extracted power of an array of buoys and an isolated buoy, respectively, θI is the angle of incidence of water waves onto an array of buoys, and C is a parameter related to the configuration of an array of buoys. In addition, all results of our multi-body WECs study are compared with those of a single buoy WEC, even though we don’t directly report the corresponding array gain. But, the better energy extraction efficiency of multi-body WECs in this comparison implies the corresponding array gain. Finally, we state that only the normal incidence of water waves onto the buoys of multi-body WECs, expressed by θI=0° or 90°, is considered in this analysis.

2.3 Resonance in Multi-Body WECs

This research investigates the energy extraction efficiency of multi- body WECs under the hydrodynamic properties of them reported by Tokić and Yue (2019). In particular, three resonance phenomena in the array of buoys of multi-body WECs are considered in this study. These resonance phenomena are explained below.

2.3.1 Rayleigh resonance

Falnes and Budal (1982) investigated the water wave power extraction and radiated water wave from a terminator array of buoys demonstrated in Fig. 2 under point absorber assumption. The radiated water waves from this terminator array of buoys can be studied under the point absorber assumption by considering the periodicity of the velocity potentials. Consequently, the wave heading angle between the incident and the respective transmitted and reflected water waves from this terminator array of buoys have the following relation:

(12) sinθm=sinθI+2πkdym,for m=0,±1,±2,
where m is the scattering mode, θm is the direction of propagation of the m-the scattering mode, and dy is the distance between two adjacent buoys oriented along the y-direction. Equation (12) shows that the scattered water waves head to the periodic buoy, i.e.,|sinθm|= 1, at the critical wave number and new scattering modes emerge. This special case is expressed by the following relation:
(13) (kdy)mcr=2π±1-sinθIm,for m=0,±1,±2,

Fig. 2.

A terminator array of buoys and the radiated water wave from it (Falnes and Budal, 1982)

Subsequently, Eq. (13) becomes the following when there are no oblique waves.

(14) (kdy)mcr=2πmfor m=1,2,

So, the wave energy is redistributed equivalently at this frequency, resulting in abrupt changes in the energy extraction. This phenomenon is called Rayleigh resonance and causes loss of the buoys’ energy extraction, especially in this terminator array of buoys (Tokić and Yue, 2019).

2.3.2 Bragg resonance

Bragg resonance is a phenomenon causing strong reflection of water waves when they pass through periodic bodies or over bottom ripples (e.g., Liu and Yue, 1998; Alam et al., 2009). Water waves passing through periodic bodies are apparently reflected when the distance between the centers of two adjacent bodies is half of the wavelength of the water wave. Hence, the condition for Bragg resonance is:

(15) kdxcosθI=nπ,for n=1,2,

Subsequently, the following equation is obtained when there are no oblique waves.

(16) kdx=nπ,forn=1,2,

Generally, Bragg resonance causes abrupt loss of the buoys’ wave energy extraction due to the strong reflection of the incident water waves occurring during it.

2.3.3 Laue resonance

The effect of Laue resonance on the energy extraction of periodic WECs was first studied by Tokić and Yue (2019). Generally, the condition for Laue resonance in an array of buoys is:

(17) (kI-km)a1=2πn,for n=1,2,
where k⃗I and k⃗m are the respective wave number vectors of the incident and the m-th scattered waves, and a⃗1 is the x-directional lattice vector of an array of buoys.

Particularly, the x-directional lattice vector of the rectangular array of buoys shown in Fig. 3 is a⃗1 = dx. Then, Eq. (17) corresponding to this rectangular array of buoys is:

(18) kdx(cosθI-cosθm)=2πn,forn=1,2,

Fig. 3.

A rectangular array of buoys and lattice vectors

Equation (18) implies that Bragg resonance is a special case of Laue resonance when θm = πθI. Finally, substituting Eq. (14) into Eq. (18) results in the following equation for the normal incidence of water waves onto this rectangular array of buoys.

(19) n(kdy)2kdx=πn2(kdy)2+πm2(kdx)2form,n=1,2,

This equation is quadratic in kdy or kdx at each (m, n) pair. Notably, Tokić and Yue (2019) reported a substantial reduction in the wave energy extraction of periodic WECs due to Laue resonance.

3. Numerical Investigation of the WECs’ Configuration

3.1 Vertical Cylinder

We consider vertical cylinders that have only the heave motion caused by their interaction with the water waves in the numerical analysis. Relatedly, the hydrodynamic properties of a vertical cylinder are well studied. For instance, the natural frequency of its heave motion due to the hydrodynamic forces is proportional to 1/T, with T denoting the draft of the vertical cylinder. In addition, the ratio of the radius of the cylinder to the draft influences the maximum response amplitude operator(RAO) of the said heave motion. Particularly, a large radius of the cylinder compared to the draft results in a small maximum RAO near the resonance frequency. On the other hand, the bandwidth of non-dimensional capture width (kW) becomes large when this ratio increases with the same maximum RAO value. This study considers vertical cylinders with a radius (a) of 1.0 m and a draft (T) of 0.67 m in the array shown in Fig. 4. The orientation of this array of cylinders with respect to the incident water waves results in different types of this array. The analysis of each of these different types is carried out as follows.

Fig. 4.

Schematic of an array of vertical cylinders

3.2 Single-column Attenuator Array

First, we compute the non-dimensional average capture width (kW) of a single-column attenuator array, i.e., θI=0°, shown in Fig. 5. The total non-dimensional capture width of buoys is divided by the number of buoys to obtain this array’s non-dimensional average capture width. Three values for the x-directional distance between adjacent cylinders (dx) are considered in this analysis: 3a, 8a, and 14a. The results of this analysis are shown in Fig. 6. In particular, the upper graphs of this figure show that the energy extraction by this array decreases abruptly near the Bragg resonance frequency, and this abrupt decrease is denoted by a red dotted line in these graphs. This effect is especially large when Nx is not small enough. In addition, this resonance effect is small in the low-frequency region, as seen from the rightmost one of these graphs. Finally, the non-dimensional average capture width of a single-column attenuator array is small compared to that of a single buoy. This reduction is large as the number of buoy increases.

Fig. 5.

Single-column attenuator array of vertical cylinders and the wall effect of periodic buoys

Fig. 6.

Non-dimensional capture width of a single-column attenuator array of vertical cylinders (Upper graphs: the average across all the buoys, Lower graphs: for each buoy in this array, and black lines: for an isolated buoy)

The lower graphs of Fig. 6 show the non-dimensional capture width of each buoy in the single-column attenuator array of vertical cylinders. The first 2 or 3 buoys achieve a large energy extraction, whereas the other buoys achieve a relatively small energy extraction except at the Bragg resonance region. This difference between the energy extractions of buoys is because the periodic buoy works like a wall, as demonstrated in the right side of Fig. 5. The water waves undergo a large reflection at the first few buoys due to the periodic buoys’ wall effect, which is concentrated at the first few buoys. It is also confirmed that there is a large reduction in the energy extraction of the first few buoys near the Bragg resonance frequency. Since buoys located at the trailing end of this array have low efficiency, this array’s non-dimensional average capture width decreases as the number of buoys increases in it.

3.3 Single-Row Terminator Array

Next, we compute the non-dimensional average capture width (kW) of a single-row terminator array of vertical cylinders shown in the left side of Fig. 7 by following a procedure similar to the one followed for the corresponding single-column attenuator array. Three values for the y-directional distance between adjacent cylinders (dy) are considered in this analysis: 3a, 8a, and 14a. The results of this analysis are shown in Fig. 8. Overall, the single-row terminator array of vertical cylinders works constructively. In particular, the energy extraction efficiency increases as the number of buoys increases, and then it converges for some number of buoys. However, energy extraction abruptly drops down near the Rayleigh resonance region, and this abrupt drop is denoted by a blue dotted line in the graphs of Fig. 8. Subsequently, this energy extraction becomes almost zero due to the emergence of new radiated water waves (Srokosz, 1980). The lower graphs of Fig. 8 show the non-dimensional capture width of each buoy in this array. These graphs indicate that each of these buoys achieves large energy extraction. In particular, the buoy in the middle of this array achieves high efficiency near the Rayleigh resonance frequency but not exactly at that frequency. Relatedly, the direction of constructive waves in this terminator array becomes the same as the array direction, i.e., |sinθm|= 1, during Rayleigh resonance. So, the wall effect of periodic buoys resulting when the constructive wave direction becomes the same as that of the array causes the water waves to concentrate between buoys due to the large reflection of them.

Fig. 7.

Single-row terminator array of vertical cylinders and the wall effect of periodic buoys

Fig. 8.

Non-dimensional capture width of a single-row terminator array of vertical cylinders (Upper graphs: the average across all the buoys, Lower graphs: for each buoy in this array, and black lines: for an isolated buoy). There is no Rayleigh resonance in 0.4–1.2 ka when dy = 3a.

Figs. 58 indicate that the wall effect induced by the periodic buoys changes the energy extraction efficiency. So, Kim and Cho (2019) used this correlation in WECs to improve their performance. In particular, they created lateral and back-side walls near only a few buoys. Hence, the lateral wall makes the buoys be in a terminator array, and the back-side wall makes buoys be similar to the first few buoys in an attenuator array, increasing the energy extraction efficiency.

3.4. Multi-row Terminator Array

The results of our analysis conducted so far indicate that the single-row terminator array of vertical cylinders shows a better energy extraction efficiency than that of the corresponding single-column attenuator array when the number of buoys is sufficiently large. Subsequently, we increase the number of rows (Nx) of this terminator array and investigate the energy extraction efficiency of the resulting arrays. Moreover, we use various dx or x-directional distance between adjacent rows in the same multi-row terminator array as the x-directional periodicity of the buoys reduces the energy extraction of this array. In particular, we study the effect of having various dx values within the same 4×8 terminator array of vertical cylinders on the energy extraction of this array under different configurations of it in detail.

3.4.1 Few rows terminator array

Generally, a multi-row terminator array of vertical cylinders can have all the already mentioned resonance phenomena occur in them. In addition, the x- and y-directional distances between adjacent rows and columns, respectively, of these arrays can result in different Bragg and Rayleigh resonance frequencies. Combining these distances also gives different Laue resonance frequencies in these arrays. The results of an analysis on the multi-row terminator arrays of vertical cylinders for dy = 3a, 8a, and 14a are given in Fig. 9. The green lines in this figure denote kdx and kdy divided by π, respectively, for different combinations of array distances (dy = dx, dy = 2dx, and dy = 0.5dx). The points of intersection between the green lines and the other blue horizontal lines, red vertical lines, and magenta curves satisfy the Bragg, Rayleigh, and Laue resonance phenomena, respectively. Subsequently, we increase the number of rows (Nx) for a fixed number of columns (Ny) of eight, as shown in Fig. 10, and analyze the resulting arrays. Notably, these Nx and Ny values are sufficient to confirm the effect of multiple scattering of water waves by the buoys in these arrays.

Fig. 9.

Resonance of a multi-row terminator array of vertical cylinders (Red line: Bragg resonance, Blue line: Rayleigh resonance, and Magenta line: Laue resonance)

Fig. 10.

Various configurations of the multi-row terminator array of vertical cylinders considered in the analysis (Number of columns (Ny): 8)

Fig. 11 shows the non-dimensional average capture width of the multi-row terminator array of vertical cylinders for all the cases considered. In particular, the increment of the number of rows in these arrays reduces their energy absorption for a small kdy, as shown in the left side graphs of this figure. This reduction is obvious when the maximum kdy is smaller than the first Rayleigh resonance frequency.

Fig. 11.

Non-dimensional average capture width of the multi-row terminator array of vertical cylinders

The maximum energy absorption of general periodic WECs is limited by the normal component of the incident energy corresponding to the WECs’ maximum average capture width,(Wa)max = d, in the absence of oblique waves. Subsequently, we compute the maximum array gain qmax for an array of N rows in the absence of oblique waves (Tokić and Yue, 2019) using the array gain given by Eq. (11) and the maximum capture width of a heaving isolated buoy(kW)max = 1:

(20) qmax=kdyN

Fig. 2 of Tokić and Yue (2019) shows that the array gain cannot always attain the maximum value. Therefore, the number of rows is crucial in the case of a small kdy to achieve array gain. Hence, the increment of the number of rows in the multi-row terminator arrays of vertical cylinders results in a decrease of the average capture width of these arrays for small kdy values.

This analysis also confirms the resonance effect in the multi-row terminator array of vertical cylinders for various combinations of Nx and Ny. The red, blue, and magenta dotted lines in Fig. 11 denote the Bragg, Rayleigh, and Laue resonances, respectively, at each frequency. In addition, the gray line of this figure denotes the simultaneous existence of two resonance phenomena, and the black vertical line in this figure denotes the occurrence of all resonances simultaneously at this frequency. First, the effect of Bragg resonance seems small in all these array configurations. Next, the magnitude of the non-dimensional average capture width oscillates between resonance points generally. In addition, a lot of resonance points occur when dx and dy are sufficiently large. Particularly, several Laue resonance points are generated under this condition. Therefore, these arrays’ non-dimensional average capture width becomes extremely oscillating if the distance between adjacent buoys of these arrays is relatively large, as shown in the upper right graph of Fig. 11. However, a large amplitude occurs only when the number of rows in these arrays is smaller than five. Further, the extraction efficiency of the multi-row terminator array of vertical cylinders decreases when the number of rows in it is larger than four generally. Consequently, a smaller value for the number of rows in these arrays gives better performance than a corresponding single-row terminator array, provided the former has a moderate kdy.

3.4.2 A 4×8 multi-row array

Tokić (2016) attempted to find the local maximum of extracted power of a multi-row terminator array by using an optimization scheme for this array’s x- and y-directional distances between its rows and columns, respectively (dx and dy), that can vary inside the array. He found that the local maximum of this extracted power could be obtained for various dx and constant dy values within the array under regular water waves. Moreover, a regular array achieves a local maximum of extracted power at least when the inclined angle of irregular water waves is small (Tokić and Yue, 2021). Tokić and Yue (2021) also concluded that a perturbed position of each buoy results in less extracted power of this array.

Our analysis of the single-row terminator array of vertical cylinders shows that a constant dy is constructive. In addition, achieving a local maximum of extracted power in a regular array under irregular water waves is possible only with an unconstrained total length of the array. On the other hand, a constrained overall length of an array together with various dx values within it results in better energy extraction efficiency of the array under some conditions. So, we investigate the extracted power of a 4×8 multi-row array of vertical cylinders with various dx values within it though under regular water waves. In addition, the distance between rows of this array is changed. In particular, we choose one regular 4×8 multi-row array of vertical cylinders and six variants of it with at least two same x-directional distances between the rows, as shown in Fig. 12. In addition, four combinations of distances are chosen for the analysis (dx = 3a, dy = 5a; dx = 3a, dy = 12a; dx = 6a, dy = 5a; and dx = 6a, dy = 12a). The ds in these arrangements is always greater than dx.

Fig. 12.

4×8 rectangular arrays (dsdx) of vertical cylinders with same and various dx values

The results of the analysis on the 4×8 multi-row arrays of vertical cylinders are summarized in Fig. 13. The black line in this figure shows the non-dimensional average capture width of a regular 4×8 multi-row array of vertical cylinders. Other colorful lines denote the results for the corresponding arrays with various dx in them. In particular, dx influences the relative magnitude of the non-dimensional average capture width between this regular array and the other corresponding arrays with various dx. In addition, a small dx is constructive only in these arrays with various dx. Notably, various small dx (2.8a < dx < 3.5a) values result in larger energy extraction efficiency of any of these arrays with various dx than that of a corresponding regular array with the same dx. However, the details of the analyses giving this finding are not reported here. The graphs corresponding to dx = 3a in Fig. 13 have the colorful lines representing the non-dimensional average capture widths of these arrays with various dx which are larger than that (represented by the black line in this figure) of the corresponding regular array. On the other hand, the black line in the graphs corresponding to dx = 6a represents the non-dimensional average capture width of a regular 4×8 multi-row array of vertical cylinders, which is greater than those (represented by the colorful lines in this figure) of the corresponding arrays with various dx values. The energy extraction efficiency of all these multi-row terminal arrays of vertical cylinders with various dx is better than that of the corresponding regular array almost in all frequency ranges. So, we expect their similar high performances under irregular water waves. Unfortunately, it is difficult to find which of the 4×8 multi-row arrays of vertical cylinders with various dx is better. However, SSX and XXS types of these arrays have relatively large and small efficiencies, respectively, in many cases.

Fig. 13.

Non-dimensional average capture width corresponding to 4×8 rectangular arrays (dsdx) of vertical cylinders with same and various dx values

5. Conclusions

The energy extraction efficiency of multi-body wave energy converters (WECs) is investigated in this study using the generalized mode approach. In particular, the effect of the WECs’ configuration on their energy extraction efficiency is parametrically studied. First, the energy extraction efficiency of the single-row terminator and single-column attenuator arrays of vertical cylinders has been computed. The wall effect induced from periodic buoys positively affected the middle buoys of this terminator array. On the other hand, the corresponding attenuator array achieves a small energy extraction due to the wall effect except at the first few buoys. Subsequently, a terminator array with a few rows of vertical cylinders has been studied, and three resonances, i.e., Bragg, Rayleigh, and Laue resonances, are checked. The increment of rows in this terminator array results in its low energy extraction efficiency when the array has a small y-directional distance between its adjacent columns. On the other hand, the energy extraction efficiency of a multi-row array of vertical cylinders oscillates between the resonance frequencies greatly when a large y-directional distance between its adjacent columns is considered, and the number of rows is below five. Finally, a 4×8 terminator array of vertical cylinders having various dx or x-directional distance between its adjacent rows is studied. This study shows that a small dx results in better energy extraction efficiency of some of these arrays with various dx than that of a corresponding regular array with the same dx.

Notes

No potential conflict of interest relevant to this article was reported.

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (2020R1I1A1A01065411)

References

Mei CC, Stiassnie MA, Yue DKP. 2005. Theory and Applications of Ocean Surface Waves: Part 1: Linear Aspects Singapore: World Scientific.
Garnaud X, Mei CC. 2009;Wave-Power Extraction by a Compact Array of Buoys. Journal of Fluid Mechanics 635:389–413. https://doi.org/10.1017/S0022112009007411 .
Falnes J, Budal K. 1982;Wave-Power Absorption by Parallel Rows of Interacting Oscillating Bodies. Applied Ocean Research 4(4):194–207. https://doi.org/10.1016/S0141-1187(82)80026-6 .
Newman JN. 1994;Wave Effects on Deformable Bodies. Applied Ocean Research 16(1):47–59. https://doi.org/10.1016/0141-1187(94)90013-2 .
Simon MJ. 1982;Multiple Scattering in Arrays of Axisymmetric Wave-Energy Devices. Part 1. A Matrix Method Using a Plane-Wave Approximation. Journal of Fluid Mechanics 120:1–25. https://doi.org/10.1017/S002211208200264X .
McIver P, Evans DV. 1984;Approximation of Wave Forces on Cylinder Arrays. Applied Ocean Research 6(2):101–107. https://doi.org/10.1016/0141-1187(84)90047-6 .
Kagemoto H, Yue DKP. 1986;Interactions Among Multiple Three-Dimensional Bodies in Water Waves: An Exact Algebraic Method. Journal of Fluid Mechanics 166:189–209. https://doi.org/10.1017/S0022112086000101 .
Linton CM, Evans DV. 1990;The Interaction of Waves with Arrays of Vertical Circular Cylinders. Journal of Fluid Mechanics 215:549–569. https://doi.org/10.1017/S0022112090002750 .
Kashiwagi M. 2000;Hydrodynamic Interactions Among a Great Number of Columns Supporting a Very Large Flexible Structure. Journal of Fluids and Structures 14(7):1013–1034. https://doi.org/10.1006/jfls.2000.0306 .
Tokić G, Yue DKP. 2019;Hydrodynamics of Periodic Wave Energy Converter Arrays. Journal of Fluid Mechanics 862:34–74. https://doi.org/10.1017/jfm.2018.911 .
Tokić G. 2016. Optimal Configuration of Large Arrays of Floating Bodies for Ocean Wave Energy Extraction (PhD Thesis). Massachusetts Institute of Technology; Cambridge, MA, U.S.A: Retrieved from http://hdl.handle.net/1721.1/104198 .
Tokić G, Yue DKP. 2021;Hydrodynamics of Large Wave Energy Converter Arrays with Random Configuration Variations. Journal of Fluid Mechanics 923:R1. https://doi.org/10.1017/jfm.2021.584 .
Choi YR. 2011;Estimation of Wave Energy Extraction Efficiency for a Compact Array System of Small Buoys. Journal of Ocean Engineering and Technology 25(1):8–13. https://doi.org/10.5574/KSOE.2011.25.1.008 .
Kim J, Cho IH. 2019;Characteristics of Wave Response in a ‘Y’ Shape Water Channel Resonator Using Resonance of Internal Fluid. Journal of Korean Society of Coastal and Ocean Engineers 31(3):170–179. https://doi.org/10.9765/KSCOE.2019.31.3.170 .
Liu Y, Yue DKP. 1998;On Generalized Bragg Scattering of Surface Waves by Bottom Ripples. Journal of Fluid Mechanics 356:297–326. https://doi.org/10.1017/S0022112097007969 .
Alam MR, Liu Y, Yue DKP. 2009;Bragg Resonance of Waves in a Two-Layer Fluid Propagating over Bottom Ripples. Part I. Perturbation Analysis. Journal of Fluid Mechanics 624:191–224. https://doi.org/10.1017/S0022112008005478 .
Srokosz MA. 1980;Some Relations for Bodies in a Canal, with an Application to Wave-Power Absorption. Journal of Fluid Mechanics 99(1):145–162. https://doi.org/10.1017/S0022112080000560 .

Article information Continued

Fig. 1.

Schematic of a multi-body WEC

Fig. 3.

A rectangular array of buoys and lattice vectors

Fig. 4.

Schematic of an array of vertical cylinders

Fig. 5.

Single-column attenuator array of vertical cylinders and the wall effect of periodic buoys

Fig. 6.

Non-dimensional capture width of a single-column attenuator array of vertical cylinders (Upper graphs: the average across all the buoys, Lower graphs: for each buoy in this array, and black lines: for an isolated buoy)

Fig. 7.

Single-row terminator array of vertical cylinders and the wall effect of periodic buoys

Fig. 8.

Non-dimensional capture width of a single-row terminator array of vertical cylinders (Upper graphs: the average across all the buoys, Lower graphs: for each buoy in this array, and black lines: for an isolated buoy). There is no Rayleigh resonance in 0.4–1.2 ka when dy = 3a.

Fig. 9.

Resonance of a multi-row terminator array of vertical cylinders (Red line: Bragg resonance, Blue line: Rayleigh resonance, and Magenta line: Laue resonance)

Fig. 10.

Various configurations of the multi-row terminator array of vertical cylinders considered in the analysis (Number of columns (Ny): 8)

Fig. 11.

Non-dimensional average capture width of the multi-row terminator array of vertical cylinders

Fig. 12.

4×8 rectangular arrays (dsdx) of vertical cylinders with same and various dx values

Fig. 13.

Non-dimensional average capture width corresponding to 4×8 rectangular arrays (dsdx) of vertical cylinders with same and various dx values