Performance of an Array of Oblate Spheroidal Heaving Wave Energy Converters in Front of a Wall

Abstract

In this paper, we investigate the performance of a linear array of five semi-immersed, oblate spheroidal heaving Wave Energy Converters (WECs) in front of a bottom-mounted, finite-length, vertical wall under perpendicular to the wall regular waves. The diffraction and radiation problems are solved in the frequency domain by utilizing the conventional boundary integral equation method. Initially, to demonstrate the enhanced absorption ability of this array, we compare results with the ones corresponding to arrays of cylindrical and hemisphere-shaped WECs. Next, we investigate the effect of the array’s distance from the wall and of the length of the wall on the physical quantities describing the array’s performance. The results illustrate that the array’s placement at successively larger distances from the wall, up to three times the WECs’ radius, induces hydrodynamic interactions that improve the array’s hydrodynamic behavior, and thus its power absorption ability. An increase in the length of the wall does not lead to any significant power absorption improvement. Compared to the isolated array, the presence of the wall affects positively the array’s power absorption ability at specific frequency ranges, depending mainly on the array’s distance from the wall. Finally, characteristic diffracted wave field patterns are presented to interpret physically the occurrence of the local minima of the heave exciting forces.

1. Introduction

Wave energy is a vast and powerful renewable energy source. Its sustainable harnessing, directly associated with the European Unions’ energy strategies and priorities for the reduction of greenhouse emissions, can contribute to climate change effects’ mitigation and energy security enhancement, and at the same time it provides the potential for long-term economic growth through the creation of jobs and the generation of export opportunities for both technology and expertise [1,2,3]. As a result, nowadays, the wave energy sector is rapidly developing with the aiming of overcoming existing technological barriers, as well as economic and social implications [3,4], and thus delivering commercially competitive solutions. In this respect, various Wave Energy Converters (WECs) characterized by a different working principle, i.e., mode of power absorption, have been proposed and designed (e.g., [5,6]). Currently, these WECs’ types have reached different technological development stages, with the oscillating water columns and the heaving WECs (point absorbers) representing the most advanced device types [7].

Although the power absorption efficiency of single-unit heaving WECs can be improved by adequately modifying its shape or its Power Take-Off (PTO) mechanism, the need to absorb significant amounts of wave power at a given location requires the deployment of multiple WECs arranged in arrays. These arrays can either be installed and operate in offshore areas, with the potential to exploit waves of higher energy content, or at near-shore locations, where specific drawbacks (e.g., survivability issues) can be overcome. Furthermore, construction, operation, and maintenance costs can be reduced. The integration of WEC technologies with other marine facilities (e.g., [8,9,10]) presents, nowadays, a potential beneficial approach for utilizing WECs, since it facilitates cost reductions through cost sharing. Following this idea, cost-efficient solutions at near-shore areas can also be realized by installing WECs arrays in front of existing coastal structures, such as vertical (wall-type) breakwaters. In these cases, the WECs’ hydrodynamic behavior and thus their power absorption efficiency can be improved by taking advantage of both the incident and the scattered from the wall boundary waves.

The diffraction and radiation problems for the case of a single, free-floating, truncated cylinder in front of a bottom-mounted wall has been examined in the frequency domain by [11] and [12], respectively, while the power absorption of a single heaving WEC in the presence of the wall boundary has been investigated and assessed, also in the frequency domain, by [13]. In the latter study, the analysis was implemented for two different WECs’ shapes (cone and hemispherical cylinder), for various sizes and drafts of the WEC, as well as for different distances of the WEC from the wall, while the action of regular and irregular waves of various wave directions was taken into account. However, it should be noted that, in all the above investigations, the method of images was utilized for taking into account the presence of the wall in the corresponding numerical model. Thus, a completely reflecting boundary of infinite extent was assumed, which although leads to positive in terms of power absorption effects [14], it does not represent adequately the real physical problem. Finally, two-dimensional approximations have been developed in [15] for describing the resonant behavior of a single, free-floating, semi-immersed cylinder in the presence of a vertical rigid wall.

With respect to WECs arrays in front of breakwaters, surge oscillating WECs placed successively in the interior of a series of U-shaped openings (harbours) embedded in a “pure” reflecting wall have been investigated in [16], while all other relevant existing studies so far deal with arrays of cylindrical-shaped heaving WECs, in contrast to the case of isolated WECs arrays, where non-cylindrical heaving devices have been also utilized (e.g., [17,18,19]). Specifically, the power absorption of a linear array of five cylindrical heaving WECs in front of a bottom-mounted vertical wall of infinite length was examined and evaluated in [20] in both frequency and time domain, by also utilizing the method of images. Recently, for an array of five truncated cylinders, Loukogeorgaki and Chatjigeorgiou [21] revealed the limitations of the “infinite wall” concept on the calculation of the cylinders’ exciting forces; moreover, they investigated the effect of the length of the wall on these quantities. The performance of nine cylindrical heaving WECs in the seaward side of a bottom-mounted vertical wall of finite length was examined in the frequency domain by [22], where emphasis was mainly given on the effect of the array’s distance from the wall on the WECs’ heave exciting forces and responses, as well as on the power absorption ability of the array. Alternatively to the existence of a bottom-mounted wall boundary, the case of a box-type truncated breakwater of infinite length behind an array of heaving WECs has been investigated analytically in [23], and the performance of nine cylindrical heaving WECs placed in front of a fixed rectangular pontoon of finite length has been studied in the frequency domain by [24]. In this study, the focus was on the effect of various geometrical parameters (e.g., dimensions of the pontoon, distance of the array from the pontoon, etc.) and the effect of the PTO characteristics on the interaction factors of both the array and the individual WECs. In addition to all the above numerical investigations, two-dimensional experiments have been conducted quite recently by [25] for an array of four cylindrical heaving WECs placed in the seaward side of a pontoon-type breakwater. The corresponding results demonstrate that by appropriately defining important design parameters of the arrangement (i.e., WECs’ draft and distance of the WECs from the breakwater), the existence of the pontoon in the leeward side of the array can significantly improve the array’s performance.

In this paper, we focus on the deployment of heaving WECs with a different shape than that of the widely deployed, so far, cylindrical one, for configuring arrays in front of a wall boundary. More specifically, we numerically investigate the performance (hydrodynamic behaviour and power absorption) of a linear array of five semi-immersed, oblate spheroidal heaving WECs placed in front of a bottom-mounted vertical wall of finite length and of negligible thickness. The WECs are distributed uniformly within the array, and each WEC is assumed to absorb power through a linear PTO mechanism, actuated from its heave motion. The analysis is implemented in the frequency domain considering the action of perpendicular to the wall regular waves. The relevant diffraction and radiation problems, taking into account the hydrodynamic interactions among the WECs and between the wall and the WECs, are solved in the frequency domain by utilizing the conventional Boundary Integral equation (BIE) method. Initially, results are compared with the ones corresponding to arrays of cylindrical and hemisphere-shaped WECs to demonstrate the enhanced absorption ability of the examined array. Next, for the array with the oblate spheroidal WECs, extended results are presented focusing on the effect of the array’s distance from the wall and of the length of the wall on the exciting forces, the WECs’ responses, and the power absorption of the array. In addition, the results are compared with the case of an isolated array to illustrate more clearly the effect of the existence of the wall on the performance of the array. Finally, characteristic diffracted wave field patterns are presented aiming at interpreting physically the occurrence of distinctive local minima of heave exciting forces at specific frequency ranges.

2. Numerical Modeling

A linear array of $M$ hydrodynamically interacting, semi-immersed heaving WECs is placed in front of a bottom-mounted vertical wall of finite length, $l_w$, and of negligible thickness, in a liquid region of constant water depth, $d$, (Figure 1). Each WECj, $j=1,\dots ,M$, of the array has an oblate spheroidal geometry with semi-major and semi-minor axes $a$ and $b$, as shown in Figure 1b. In this figure, the PTO mechanism of each WEC is schematically represented as a linear damping system, with damping coefficient $b_{PTOj}$, $j=1,\dots ,M$. All WECs are situated at a distance $c$ from the wall and they are distributed uniformly with a center-to-center spacing equal to $l_{bet}$ (see Figure 1a). Both the array and the wall are subjected to the action of monochromatic incident waves of circular frequency $\omega$ and linear amplitude $A$, which propagate at an angle $\beta$ relatively to the global horizontal $X$-axis (Figure 1a).

The hydrodynamic analysis of the array in front of the wall, taking into account the hydrodynamic interactions between the WECs and between the wall and the WECs, is performed in the frequency domain and it is based on the BIE method, which is numerically realized using the WAMIT V6.3 © software [26]. The employed methodology relies on the three-dimensional linear potential theory, where the wall is taken fixed at its position, and all WECs are assumed to oscillate freely only along their working direction, i.e., along the vertical $z$ axis (Figure 1b). Accordingly, for all WECs, all degrees of freedom, except the one corresponding to heave, are considered ideally restricted. The complex spatial part of the velocity potential is given by [26,27]:

$$\phi =\underset{{\phi }_D}{\underbrace{({\phi }_I+{\phi }_S)}}+i\omega {\displaystyle \sum }_{j=1}^M{\xi }_3^j{\phi }_j$$

(1)

where the incident wave component, ${\phi }_I$, is obtained from:

$${\phi }_I=\frac{igA}{\omega }\frac{cosh\left[k\left(Z+d\right)\right]}{cosh\left(kd\right)}{e}^{-ik\left(Xcos\beta +{\rm Y}sin\beta \right)}$$

(2)

Furthermore, ${\phi }_D$ denotes the diffraction potential, ${\phi }_S$ is the scattered potential related to the scattered disturbance of the incident waves from the WECs and the wall, and ${\phi }_j$, $j=1,\dots ,M$, correspond to the radiation potentials, associated with the waves radiated from the WECs due to their forced motion in heave. In that respect, ${\xi }_3^j$, $j=1,\dots ,M$, denote the complex amplitudes of the heave motions of the WECs. Finally, $g$ is the gravitational acceleration, $k$ is the wave number, and $i^2=-1$.

Within the realm of the linear potential theory, the velocity potentials ${\phi }_q$ ($q=D$ or $q=j$) should satisfy the Laplace equation and the following linearized boundary conditions:

$$\frac{\partial {\phi }_q}{\partial Z}-\frac{{\omega }^2}{g}{\phi }_q=0\mathrm{on}Z=0$$

(3)

$$\frac{\partial {\phi }_q}{\partial Z}=0\mathrm{on}Z=-d$$

(4)

$$\frac{\partial {\phi }_D}{\partial n}=0$$

(5)

$$\frac{\partial {\phi }_j}{\partial n}=n_3^j\mathrm{for}j=1,\dots ,M$$

(6)

Equation (3) is the combined kinematic and dynamic free-surface condition, Equation (4) is the bottom boundary condition, and Equations (5) and (6) are the Neumann conditions, which should hold on the wetted surfaces of the bodies. In Equation (6), $n_3^j$ denotes the normal unit vector of WECj in the vertical direction.

The boundary integral equations for the unknown diffraction and radiation potentials on the boundaries of all bodies (WECs and wall) and of the WECs, respectively, are formed by utilizing the Green’s theorem. Then, the corresponding boundary value problem is solved based on a three-dimensional low-order panel method [26,27]. For modeling the wetted surface of the wall, zero-thickness dipole panels [27] are deployed, given that the thickness of the wall is assumed to be negligible.

Following the solution of the first-order boundary value problem, the heave exciting force, $F_3^i$, $i=1,\dots ,M$, applied on each WECi, as well as the added mass, $A_{ij}$, and the radiation damping, $B_{ij}$, coefficients are obtained using Equations (7) and (8), respectively, where $S_b^i$ is the wetted surface of the $i$th WEC and $\rho$ is the water density:

$$F_3^i=-i\omega \rho {{\displaystyle \iint }}_{S_b^i}^{}{n}_3^i{\phi }_{D}ds,i=1,\dots ,M$$

(7)

$$A_{ij}-\frac{i}{\omega }{{\rm B}}_{ij}=\rho {{\displaystyle \iint }}_{S_b^i}^{}{n}_3^i{\phi }_{j}ds,i,j=1,\dots ,{\rm M}$$

(8)

Having calculated the first-order hydrodynamic forcing quantities, the following linear system of equations is solved to obtain the complex amplitudes of the WECs’ heave motions, ${\xi }_3^j$, $j=1,\dots ,M$:

$${\displaystyle \sum }_{j=1}^M\left[-{\omega }^2\left({{\rm M}}_{ij}+A_{ij}\right)+i\omega \left({{\rm B}}_{ij}+B_{ij}^{PTO}\right)+C_{ij}\right]{\xi }_3^j=F_3^{i}i=1,\dots ,M$$

(9)

In Equation (9), $M_{ij}$ are the mass matrix coefficients, $B_{ij}^{PTO}$ represent the damping coefficients originating from the PTO mechanism, and $C_{ij}$ correspond to the hydrostatic-gravitational stiffness coefficients. Given that each WEC is assumed to oscillate only in heave, $C_{ij}=0$ for $i\ne j$, and $C_{ij}$ for $i=j=1,\dots ,M$ are obtained as follows:

$$C_{ij}=\rho g{{\displaystyle \iint }}_{S_b^i}^{}{n}_3^{i}ds,i=j=1,\dots ,M$$

(10)

Regarding the coefficients $B_{ij}^{PTO}$, the PTO mechanism for a $j$th WEC of the array is modeled as a linear damping system (Figure 1b), with damping coefficient $b_{PTOj}$, actuated from the heave motion of the corresponding WEC. Consequently, in Equation (9), $B_{ij}^{PTO}=b_{PTOj}$ for $i=j=1,\dots ,{\rm M}$, while $B_{ij}^{PTO}=0$ for $i\ne j$.

The heave response of each WEC is expressed in terms of the Response Amplitude Operator, Equation (11), where $\left|{\xi }_3^j\right|$ denotes the amplitude of the complex quantity ${\xi }_3^j$:

$$RA{O}_3^j=\frac{\left|{\xi }_3^j\right|}{A}$$

(11)

The mean power, $p\left(\omega \right)$, absorbed by the whole array at a specific $\omega$ is, then, calculated as:

$$p\left(\omega \right)={\displaystyle \sum }_{j=1}^M{p}_j\left(\omega \right)$$

(12)

where $p_j\left(\omega \right)$,$j=1,\dots ,{\rm M}$, is the power absorbed by the $j$th WEC, calculated using Equation (13):

$$p_j\left(\omega \right)=0.5b_{PTOj}{\omega }^2{\left|{\xi }_3^j\right|}^2$$

(13)

Finally, the non-dimensional wave elevation, $\eta \left(X,Y\right)/A$, at a given field point situated at $X,Y$, is obtained from

$$\frac{\eta \left(X,Y\right)}{A}=\left(\frac{{\eta }_D\left(X,Y\right)}{A}+\frac{{\eta }_r\left(X,Y\right)}{A}\right)={\left(\overline{{\phi }_D}\left(X,Y\right)+\frac{{\omega }^2}{g}a{\displaystyle \sum }_{j=1}^{M}RA{O}_3^j\overline{{\phi }_j}\left(X,Y\right)\right)}_{Z=0}$$

(14)

where ${\eta }_D\left(X,Y\right)$ and ${\eta }_r\left(X,Y\right)$ correspond to the wave field due to diffracted and radiated waves, respectively, and $\overline{{\phi }_q}$, $q=D$ or $j$, denote non-dimensional quantities. The values of $\overline{{\phi }_D}$ and $\overline{{\phi }_j}$ are obtained by dividing the corresponding dimensional quantities with $iAg/\omega$ and $a$, respectively.

3. Results and Discussion

3.1. General Remarks

The BIE-based numerical model described in Section 2 is applied for an array of $M=5$ identical, semi-immersed oblate spheroidal WECs (Figure 1) with semi-major axis $a=2.0$ m and non-dimensional semi-minor axis (draft) $b/a=0.85$. These dimensions have been selected based on [28], who demonstrated the enhanced, as compared with other WEC geometries (i.e., cylinder, hemisphere, and prolate spheroid), absorption ability of a corresponding single isolated oblate spheroidal WEC at infinite water depths. All five WECs are considered to have the same PTO characteristics; thus, $B_{ij}^{PTO}=b_{PTO}$, $i=j=1,\dots ,{\rm M}$. Moreover, the damping coefficient of the PTO is tuned to a constant value, so that maximum energy absorption is achieved at the natural frequency of a single isolated WEC (e.g., [22,29,30,31]). Accordingly and in line with [32], the aforementioned $b_{PTO}$ value is taken equal to the heave radiation damping of a single isolated WEC at its heave natural frequency, ${\omega }_{n3}^{ISO}$, i.e., $b_{PTO}=B_{33}^{ISO}\left(\omega ={\omega }_{n3}^{ISO}\right)$. On the basis of the hydrodynamic analysis of an isolated oblate spheroidal WEC, ${\omega }_{n3}^{ISO}$ has been calculated equal to $2.4$ rad/s, resulting in $b_{PTO}=10,322.20$ Nm/s.

Initially, relevant results are compared with those corresponding to arrays of semi-immersed cylindrical and hemisphere-shaped WECs (Figure 2) in order to demonstrate the enhanced absorption ability of the aforementioned array. Aiming at comparing WECs that have the same ${\omega }_{n3}^{ISO}$, irrespectively of their geometry, the geometrical characteristics of the cylindrical and the hemisphere-shaped WECs have been defined, so that ${\omega }_{n3}^{ISO}$= 2.4 rad/s for both geometries, as in the case of the oblate spheroidal WECs. These characteristics along with the corresponding $b_{PTO}$ values are shown in

Table 1. Geometrical and $b_{PTO}$ characteristics of WEC geometries considered.

Geometry/Shape	$\mathit{a}$ (m)	$\mathit{b}$ (m)	${\mathit{b}}_{\mathit{P}\mathit{T}\mathit{O}}(\mathbf{Ns}/\mathbf{m})$
Oblate spheroid	2.0	1.7	$\mathrm{10,322.20}$
Cylinder	1.4	1.0	$2572.46$
Hemisphere	1.8	1.8	$7111.86$

, where the notations of Figure 2 are taken into account. For both the cylindrical and the hemisphere-shaped WECs, the$b_{PTO}$ values of Table 1, equal to $B_{33}^{ISO}\left(\omega ={\omega }_{n3}^{ISO}\right)$, have been obtained from the hydrodynamic analysis of the corresponding isolated body. This, in turn, leads to the consideration of different $b_{PTO}$ values among the examined geometries, due to the existence of different hydrodynamic properties (Figure 3).

For each of the three examined arrays, the WECs are placed in front of the wall at a distance $c$ of $6$ m and they are distributed uniformly within the array with non-dimensional center-to-center spacing equal to $l_{bet}/a=4$. Comparisons are made for two cases of non-dimensional length of the wall, $l_w/a$ equal to $18$ and 72 and $l_{edge}/a=1$ and $28$, respectively, with $l_{edge}$ denoting the distance between the centers of the two outer WECs in the array with the corresponding edges of the wall, as shown in Figure 1a. It is noted that for a given WEC geometry, non-dimensional quantities are defined using the corresponding value of $a$, included in Table 1.

Next, we investigate and assess the effect of the ratios $c/a$ and $l_w/a$ on the hydrodynamic behavior and the power absorption ability of the array of the five, semi-immersed oblate spheroidal WECs. For the first design parameter, six different values of $c/a$ (

Table 2. Examined cases for the array of the five semi-immersed oblate spheroidal WECs.

Examined Design Parameter	$\mathit{c}/\mathit{a}$	${\mathit{l}}_{\mathit{w}}/\mathit{a}$	${\mathit{l}}_{\mathit{e}\mathit{d}\mathit{g}\mathit{e}}/\mathit{a}$	${\mathit{l}}_{\mathit{b}\mathit{e}\mathit{t}}/\mathit{a}$
Non-dimensional distance from the wall ($c/a)$	$1.5$, $2$, $2.5$, $3$, $3.5$, $4$	$36$	$10$	$4$
Non-dimensional wall length ($l_w/a$)	$3$	$18$, $36$, $72$	$1$, $10$, $28$	$4$

) are examined for $l_w/a=36$ and $l_{bet}/a=4$, and for the second examined design parameter, three different values of $l_w/a$ are considered (Table 2) for $c/a=3$ and $l_{bet}/a=4$. Results for the isolated array (i.e., $l_w/a=0$) are also cited for illustrating more clearly the effect of the presence of the wall on the performance of the array. It is mentioned that the value of $l_{bet}/a=4$ has been chosen based on preliminary tests that revealed adequate energy absorption ability for this center-to-center spacing as compared with other $l_{bet}/a$ values examined.

In all the above cases, the array and the wall are considered to be situated in a liquid region with $d=10$ m and they are subjected to regular incident waves propagating at angle $\beta =270$ deg (Figure 1a) with $\omega$ varying between 0.01 and 4.0. The heave exciting forces are given normalized by $\rho gA\pi a^2$, while the symmetry of the array with respect to the examined incident wave direction leads to the same exiting forces and responses for WEC1 and WEC5 (outer WECs) and for WEC2 and WEC4 (Figure 1a with $M=5$).

3.2. Comparison of Arrays with Different WEC Geometries

In Figure 4, the heave exciting forces, $F_3$, applied on the semi-immersed oblate spheroidal WECs of the examined five-body array are compared with the corresponding ones applied on the WECs of the arrays that consist of semi-immersed cylindrical and hemisphere-shaped WECs (Figure 2). The results of Figure 4 correspond to the case of $l_w/a=18$, while analogous conclusions can be drawn for $l_w/a=72$. For all the oblate spheroidal WECs, $F_3$ starts its variation from the limiting value of one (at $\omega =0.01$ rad/s), as exists in the case of the isolated array (see for example Figure 7 in Section 3.3), and accordingly it increases quite rapidly up to $\omega ~0.75$ rad/s, where the $F_3$ peak (global maximum) occurs. A rapid decrease follows leading to a local minimum of $F_3$ at $1.5$ rad/s $<\omega <1.8$ rad/s, and for the remaining examined $\omega$ values, $F_3$ varies quite intensively, having multiple peaks and successive local minima, with successively decreasing values towards higher frequencies. The utilization of hemisphere-shaped WECs does not introduce significant differences in the values and the variation patterns of the heave exciting forces. However, it is worthwhile to note that for the specific WEC geometry the observed global maximum is shifted at a slightly larger $\omega$ as compared with the oblate spheroidal WECs. This shift becomes more pronounced in the case of the cylindrical WECs; moreover, $F_3$ applied on the cylindrical WECs varies more smoothly at $\omega >3.0$ rad/s as compared with the oblate spheroidal and the hemisphere-shaped WECs. The above characteristics lead to different $F_3$ values of the cylindrical WECs as compared with the other two examined geometries at specific $\omega$ values (e.g., at$0.5$ rad/s $<\omega <1.0$ rad/s for WEC2 and WEC4, Figure 4b, at$3.0$ rad/s $<\omega <3.5$ rad/s for WEC3, Figure 4c).

Figure 5 shows the comparison of $RA{O}_3$ between the three examined WEC geometries. Again, the results correspond to arrays placed in front of a wall with $l_w/a=18$, while similar conclusions can be drawn for the case of $l_w/a=72$. Irrespectively of the WECs’ geometry and the position of the WEC within the array, $RA{O}_3$ follows, in general, the variation pattern of the heave exciting forces which is characterized by the occurrence of two distinctive peaks; one at $0.5$ rad/s $<\omega \le 1.0$ rad/s as a result of the existence of the peaks in $F_3$ in the concerned frequency range (Figure 4) and a second peak at $2.0$ rad/s $<\omega <2.5$ rad/s attributed to resonance phenomena. Moreover, a local minimum is observed at $1.5$ rad/s $<\omega <2.0$ rad/s, in absolute accordance with the variation pattern of the heave exciting forces (Figure 4). Except at $2.0$ rad/s $<\omega <2.5$ rad/s, the deployment of hemisphere-shaped WECs within the array does not lead to substantial differences of $RA{O}_3$ as compared with the case of the oblate spheroidal WECs. In the concerned frequency range, where resonance phenomena occur, the consideration of a smaller $b_{PTO}$ value for the hemispheres (Table 1) has a more pronounced effect on $RA{O}_3$, leading, therefore, to slightly larger values for $RA{O}_3$ as compared with those of the oblate spheroids. In a similar manner, the utilization of an even smaller $b_{PTO}$ value for the examined cylindrical WECs (Table 1) leads to a significant increase of $RA{O}_3$ at $2.0$ rad/s $<\omega <2.5$ rad/s in the case of WECj, $j=1,2,4,$ and $5$ as compared with the corresponding oblate spheroidal and the hemisphere-shaped WECs. The same holds for WEC3 (Figure 5c) at $1.8$ rad/s $<\omega <2.7$ rad/s, except at $\omega =2.4$ rad/s, where $RA{O}_3$ of WEC3 shows a local minimum in absolute accordance with the corresponding variation of the heave exciting force (Figure 4c). Finally, for all cylindrical WECs, the existence of different $F_3$ values at specific frequency ranges (e.g., at $0.5$ rad/s $<\omega <1.0$ rad/s, Figure 4) as compared with the rest examined WEC geometries leads also to different $RA{O}_3$ values at these frequency ranges.

As far as the absorbed mean power is concerned, Figure 6 shows the comparison of $p\left(\omega \right)$ among the three examined five-body arrays. For both $l_w/a=18$ and $72$ and for all three arrays, the variation pattern of $p\left(\omega \right)$ is characterized by the existence of a first (local) peak at $0.9$ rad/s $<\omega <1.3$ rad/s and a second (global) peak at $2.0$ rad/s $<\omega <2.5$ rad/s, in accordance with the variation of $RA{O}_3$, as previously discussed for the case of $l_w/a=18$. Moreover, local minima of $p\left(\omega \right)$ occur at specific $\omega$ values as a result of the relevant substantial decrease of $RA{O}_3$. Compared to the cylindrical WECs, the deployment of either the hemisphere-shaped or the oblate spheroidal WECs greatly improves the power absorption ability of the array, since for the latter WEC geometries significantly larger $p\left(\omega \right)$ peak values occur, and, additionally, for $\omega >1.8$ rad/s the frequency range, where adequate amount of power is absorbed, becomes wider. All the above are attributed to the consideration of a quite smaller $b_{PTO}$ value in the case of the cylindrical WECs (Table 1), as a result of their intrinsic hydrodynamic characteristics. Therefore, although this $b_{PTO}$ value results in a significant increase of $RA{O}_3$ at the frequency range where resonance phenomena occur (Figure 5), it reduces remarkably the power absorption ability of the array. Following a similar rationale, the deployment of hemisphere-shaped WECs within the array leads to slightly smaller $p\left(\omega \right)$ peak values as compared with the oblate spheroidal WECs. Consequently, the oblate spheroidal WECs can be considered to have the best power absorption ability among all the three examined WEC geometries.

Finally, it is worthwhile to mention that for a given WEC’s geometry, the length of the wall affects at a small degree $p\left(\omega \right)$. Specifically, in the case of $l_w/a=18$ (Figure 6a), the existence of more intense wall-edge effects reduces the power absorption ability of the array at $2.0$ rad/s $<\omega <2.5$ rad/s, leading to slightly smaller $p\left(\omega \right)$ peak values in this frequency range as compared with $l_w/a=72$ (Figure 6b). However, at smaller frequencies, wall-edge effects seem to have a positive effect on $p\left(\omega \right)$, since at $0.8$ rad/s $<\omega <1.3$ rad/s slightly larger $p\left(\omega \right)$ values are observed for $l_w/a=18$ as compared with $l_w/a=72$.

3.3. Effect of WECs’ Distance from the Wall on the Array’s Performance

The effect of the distance of the semi-immersed oblate spheroidal WECs from the wall on the WECs’ heave exciting forces is shown in Figure 7, where the variation of $F_3$ for all WECs in the array as a function of $\omega$ is presented for all $c/a$ values examined, as well as for the isolated array. For WECj, $j=1,2,4$, and $5$ (Figure 7a,b) and for $c/a=1.5$, the variation of $F_3$ begins from the limiting value of $~1.0$ at $\omega =0.01$ rad/s and it is characterized by a rapid increase up to $\omega ~0.4$ rad/s, where the peak of $F_3$ (global maximum) occurs. Accordingly, $F_3$ decreases quite rapidly up to $\omega ~2.3$ rad/s, where it obtains its first local minimum, which approaches zero. Finally, $F_3$ is increased resulting in a second peak at $\omega ~3$ rad/s. Analogous variation pattern is observed for all the examined $c/a$ values. However, by successively increasing the array’s distance from the wall, the first local minimum of $F_3$ and, thus, its second peak are successively shifted at smaller frequencies, while the values of the latter peaks are also consecutively increased. The above trends introduce differences on the values of $F_3$ at specific frequency ranges. For example, the increase of $c/a$ leads to larger $F_3$ values at $2.0$ rad/s $<\omega <2.5$ rad/s, while the opposite holds true at $1.0$ rad/s $<\omega <1.5$ rad/s. Analogous conclusions can be drawn for the heave exciting forces applied on the middle WEC of the array (Figure 7c). However, $F_3$ of WEC3 obtains one more distinctive peak at $\omega ~1.4$ rad/s for the two smallest examined WECs’ distances from the wall (i.e., $c/a=1.5$ and $2$).

The results of Figure 7 demonstrate clearly that the presence of the wall, irrespectively of the distance of the WECs from this boundary, affects significantly the heave exciting forces of all the WECs in the entire examined frequency range, since for all $c/a$ values, $F_3$ does not show the continuous smooth decrease as in the case of the isolated array. Compared to the latter array, the existence of the wall boundary in the leeward side of the WECs at small distances from them (i.e., $c/a=1.5$ and $2.0$) leads also to a substantial increase of $F_3$ at $0.01$ rad/s $<\omega \le 1.5$ rad/s and at very high frequencies ($2.8$ rad/s $<\omega <3.2$ rad/s), while the placement of the array at further distances from the wall (i.e., $c/a=2.5$, $3.0$, $3.5$, and $4.0$) results in a significant increase of $F_3$ at $0.01$ rad/s $<\omega \le 1.0$ rad/s, as well as at $2.0$ rad/s $<\omega \le 2.5$ rad/s, where resonance phenomena of the bodies are anticipated. From a physical point of view, the above trends can be related to the formation of characteristic patterns of the diffracted wave field around the WECs, as commended in Section 3.5.

Figure 8 shows the effect of$c/a$ on $RA{O}_3$ for all the oblate spheroidal WECs. In this figure, results corresponding to the isolated array are also included to demonstrate more clearly the effect of the presence of the wall on$RA{O}_3$. In the case of $c/a=1.5$ and for WECj, $j=1,2,4$, and $5$ (Figure 8a,b), the variation pattern of $RA{O}_3$ is quite intense and it is characterized by three distinctive peaks. The first one (global maximum) and the third one occur at $\omega ~0.4$ rad/s and $\omega ~3.1$ rad/s, respectively, as a result of the relevant $F_3$ maximum values (Figure 7), while the second peak is observed at $\omega ~2.5$ rad/s and it is related to resonance phenomena. Moreover, a characteristic local minimum, following the preceding rapid decrease of $RA{O}_3$, occurs at $\omega ~2.3$ rad/s, in absolute accordance with the variation of the heave exciting forces (Figure 7). A similar variation pattern is observed for all examined $c/a$ values. However, by successively increasing $c/a$ up to $3.0$, the aforementioned local minimum and resonance-related peak of $RA{O}_3$ are successively shifted at slightly smaller frequencies as compared with the case of $c/a=1.5$, while the values of the latter peaks are increased substantially. Same conclusions can be drawn when $c/a$ is further increased up to$4.0$; nevertheless, the placement of the WECs at one of the two largest examined distances from the wall does not introduce any significant differences in the peak values of $RA{O}_3$ at $2.0$ rad/s $<\omega \le 2.5$ rad/s as compared with $c/a=3.0$. It is also worthwhile to note that at $0.4$ rad/s $<\omega \le 1.5$ rad/s the increase of $c/a$ leads to a smooth reduction of $RA{O}_3$ values, while only the curve of $c/a=2.0$ has a distinctive third $RA{O}_3$ peak at very high frequencies (i.e., at $\omega ~3.0$) as in the case of $c/a=1.5$. All the above are in absolute accordance with the results of Figure 7. Finally, with regard to WEC3, which is situated in the middle of the array (Figure 8c), conclusions similar to the case of WECj, $j=1,2,4$, and $5$ can be drawn. However, for all $c/a$ values examined, $RA{O}_3$ obtains one more distinctive peak at $\omega ~1.5$ rad/s.

Compared to the isolated array, the results of Figure 8 illustrate that the presence of the wall, has a significant impact on $RA{O}_3$ of all the WECs in almost the entire examined frequency range. Specifically, on the one hand, when the array is placed in front of the wall at a small distance from the wall (i.e., $c/a=1.5$ and $2.0$), the WECs’ heave responses increase significantly mainly at $0.01$ rad/s $<\omega \le 1.5$ rad/s, while a great reduction of $RA{O}_3$ occurs at $2.0$ rad/s $<\omega \le 2.5$ rad/s, especially for $c/a=1.5$. On the other hand, by increasing $c/a$ to values larger than 2, $RA{O}_3$, increases at $0.01$ rad/s $<\omega \le 1.0$ rad/s, but, most importantly in the frequency range, where resonance phenomena occur.

Finally, the effect of $c/a$ on the power absorbed by the whole array is shown in Figure 9, where $p\left(\omega \right)$ is plotted as a function of $\omega$ for all $c/a$ values examined, as well as for the case of the isolated array. For $c/a=1.5$, a significant amount of power is absorbed at $0.5$ rad/s $\le \omega \le 2.0$ rad/s, where $p\left(\omega \right)$ obtains its maximum value approximately equal to $140$ kW/m². Moreover, additional $p\left(\omega \right)$ peaks with smaller values occur at $\omega =2.5$ rad/s and $\omega =3.1$ rad/s as a result of the relevant $RA{O}_3$ peak values (Figure 8). By increasing $c/a$to $2.0$, a significant amount of power is absorbed at larger wave frequencies and, more specifically, at $2.1$ rad/s $<\omega <2.7$ rad/s, where resonance phenomena also occur. Moreover, the power absorption ability of the array is substantially enhanced, since the maximum value of $p\left(\omega \right)$, occurring at $\omega =2.45$ rad/s, is approximately equal to $270$ kW/m². Similar behavior is observed for the rest of the examined $c/a$ values. However, the increase of $c/a$from 2.0 to $3.0$ leads to a successive improvement of the array’s power absorption ability, as $p\left(\omega \right)$ global maxima become equal to $370$ kW/m² (at $\omega =2.4$ rad/s) and $385$ kW/m² (at $\omega =2.3$ rad/s) for $c/a=2.5$ and $c/a=3.0$, respectively. Moreover, the frequency range, where a significant amount of power is absorbed, becomes more and more wider (i.e., $1.9$ rad/s $<\omega <2.7$ rad/s for $c/a=2.5$ and $1.7$ rad/s $<\omega <2.7$ rad/s for $c/a=3.0$). By further increasing $c/a$ up to $4.0$, the power absorption ability of the array is successively reduced as compared with the case of $c/a=3.0$ ($p\left(\omega \right)$ maxima are approximately equal to $380$ kW/m² and $340$ kW/m² for $c/a=3.5$ and $c/a=4.0$, respectively) and it is realized at a slightly smaller frequency ranges (i.e., at $1.6$ rad/s $<\omega <2.5$ rad/s for $c/a=3.5$ and at $1.5$ rad/s $<\omega <2.4$ rad/s for $c/a=4.0$). For $c/a\ge 2.0$, adequate power is also absorbed at $\omega <1.5$ rad/s, where additional $p\left(\omega \right)$ peaks are observed. However, by increasing $c/a$, the corresponding power absorption ability of the array is successively reduced, and it is bounded at less wide frequency ranges. Taking all the above into consideration, it can be concluded that for a wall with $l_w/a=36$, the placement of the examined five-body array with $l_{bet}/a=4$ at a non-dimensional distance from the wall, $c/a$, equal to 3.0 leads to the best power absorption. Finally, compared to the isolated array, it is clear that the existence of the wall boundary positively affects the power absorption ability of the array leading to a significant increase of $p\left(\omega \right)$ at specific frequency ranges, as well as to more than one $p\left(\omega \right)$ peak, depending upon the value of $c/a$.

The results of Figure 9 demonstrate that $c/a$ is a critical design parameter that directly affects the power absorption ability of the examined five-body array. For the smallest examined non-dimensional distance from the wall, the power absorption ability of the array is not driven by resonance phenomena, as significant heave exciting forces and responses, resulting from the hydrodynamic interactions of the WECs with the wall boundary, exist at wave frequencies outside the frequency range ($2.1$ rad/s $<\omega \le 2.5$ rad/s), where resonance phenomena occur. However, the opposite holds true for the rest of the examined$c/a$ values, where maximum $p\left(\omega \right)$ values occur at $2.1$ rad/s $<\omega \le 2.5$ rad/s, since the presence of the wall boundary does not have a negative effect on the WECs’ heave exciting forces in this frequency range, and thus on the relevant $RA{O}_3$ amplification due to resonance. Up to $c/a=3.0,$ the placement of the array at successively larger distances from the wall induces hydrodynamic interactions between the WECs and the boundary that enhance consecutively the hydrodynamic behavior of the WECs, and thus the array’s power absorption ability. However, a further increase of $c/a$ to values larger than $3.0$ does not enable to exploit in the best possible way the disturbances induced from the wall wave, and thus it leads to arrays with consecutively reduced power absorption ability. Consequently, among the examined $c/a$ values, $c/a=3.0$ is considered to be the upper limit of this design parameter, in terms of power absorption enhancement.

3.4. Effect of the Length of the Wall on the Array’s Performance

The effect of the length of the wall on the hydrodynamic behavior of the semi-immersed oblate spheroidal WECs is shown in Figure 10 and Figure 11, where the variations of $F_3$ and $RA{O}_3$ for all WECs in the array as a function of $\omega$ are presented, respectively, for $l_w/a=18,36$, and $72$, as well as well as for the isolated array.

Starting with the heave exciting forces (Figure 10), $F_3$ for all WECs and for $l_w/a=72$ varies in a similar manner as in the case of $l_w/a=18$ and $36$, already described in Section 3.2 and Section 3.3, respectively. However, for this $l_w/a$ value, a second peak at $\omega =0.95$ rad/s occurs for WECj, $j=2,4$, and $3$ (Figure 10b,c). The increase of $l_w/a$ leads to a more rapid increase of $F_3$, from the limiting value of $1.0$ at $\omega =0.01$ rad/s, and thus to a shift of the $F_3$ first peaks (global maxima) at smaller wave frequencies. Moreover, it results to larger values of the aforementioned $F_3$ peaks. The above trends introduce differences in the values of $F_3$ at $0.1$ rad/s $<\omega <2.0$ rad/s in the case of WECj, $j=1,5$, and $3$ (Figure 10a,c) and at $0.1$ rad/s $<\omega <1.0$ rad/s in the case of WECj, $j=2$ and $4$ (Figure 10b). Compared to the isolated array, it is straightforward that the presence of the wall, irrespectively of its length, affects the WECs’ heave exciting forces in the whole examined frequency range, since for all $l_w/a$ values examined, $F_3$ does not demonstrate a continuous smooth decrease as in the case of the isolated array. This in turn leads to different $F_3$ values, especially at $0.01$ rad/s $<\omega <1.0$ rad/s, where for $l_w/a\ne 0$ a great increase of the heave exciting forces occurs as compared with the isolated array.

With regards to WECs’ heave responses (Figure 11), $RA{O}_3$ for all WECs and for $l_w/a=72$ varies in a similar manner as in the case of $l_w/a=18$ and $36$ already described in Section 3.2 and Section 3.3, respectively. Nevertheless, $RA{O}_3$ for WECj, $j=2,4$, and $3$ (Figure 11b,c) depicts a second peak at $\omega =0.95$ rad/s in absolute accordance with the corresponding heave exciting forces. On the one hand, for the outer WECs of the array (Figure 11a), slightly larger $RA{O}_3$ values compared to $l_w/a=18$ for both $l_w/a=36$ and $72$ are observed at $2.0$ rad/s $<\omega <2.5$ rad/s, where resonance phenomena occur, since the corresponding $F_3$ values also become larger at this frequency range (Figure 10a). On the other hand, for all WECs, the increase of $l_w/a$ leads to larger $RA{O}_3$ values at $0.1$ rad/s $<\omega <0.5$ rad/s, i.e., in the low frequency range, in accordance to the results of Figure 10. Finally, compared to the isolated array, the wall boundary, irrespectively of its length, significantly increases the WECs’ responses at $0.1$ rad/s $<\omega <1.0$ rad/s and at $2.0$ rad/s $<\omega <2.5$ rad/s, whereas the opposite holds true at $1.5$ rad/s $<\omega <2.0$ rad/s.

The effect of $l_w/a$ on the mean power absorbed by the whole array is shown in Figure 12, where $p\left(\omega \right)$ is plotted as a function of $\omega$ for $l_w/a=18,36$ and $72$, as well as for the case of the isolated array. Since for all $l_w/a$ cases examined, the array is placed at a distance $c/a=3.0$ from the wall, the power absorption ability of the array is mainly driven by the WECs’ resonance. Therefore, maximum $p\left(\omega \right)$ values occur at $2.0$ rad/s $<\omega <2.5$ rad/s, where resonance phenomena occur. However, the presence of the wall boundary introduces a second peak in the vicinity of $\omega =1.0$ rad/s due to the existence of significant heave exciting forces, and thus large $RA{O}_3$ values in the examined low frequency range, as previously described. Compared to $l_w/a=18$, where the maximum value of $p\left(\omega \right)$ is approximately equal to $360$ kW/m², the increase of $l_w/a$ to $36$ enhances, to a small extent, the power absorption ability of the array, since the maximum value of $p\left(\omega \right)$ becomes equal to $385$ kW/m². A further increase of $l_w/a$ does not lead to any improvement of the array’s power absorption ability as compared with $l_w/a=36$. As for $0.5$ rad/s $<\omega <1.3$ rad/s, where the second peak of $p\left(\omega \right)$ is observed, the change of $l_w/a$ has an insignificant effect on the values of $p\left(\omega \right)$. Finally, compared to the isolated array, it is clear that the existence of the wall boundary positively affects the power absorption ability of the array, as it results in a significant increase of $p\left(\omega \right)$, especially in the frequency range, where resonance phenomena occur.

3.5. Spatial Variation of the Diffracted Wave Field

For physically interpreting the existence of great differences of the heave exciting forces at specific frequency ranges for different distances of the array from the wall (Figure 7), the spatial variation of the non-dimensional wave elevation, ${\eta }_D/A$, due to only diffracted waves, i.e., first term in Equation (14), is taken into account. Specifically, for the array of the semi-immersed oblate spheroidal WECs with $l_{bet}/a=4$ placed in front of a wall of $l_w/a=36$, Figure 13 and Figure 14 show the ${\eta }_D/A$ spatial variation of the diffracted wave field indicatively for $c/a=1.5$ and $3.0$ calculated at $-12\le X/a\le 12$ and $0\le Y/a\le 8.0$ for two characteristic values of $\omega$, equal to $1.5$ rad/s and $2.3$ rad/s, respectively. At $\omega =1.5$ rad/s, the WECs’ heave exciting forces for $c/a=1.5$ show significant values as compared with $c/a=3.0$, while the opposite holds for $\omega =2.3$ rad/s, where, moreover, $F_3$ for all WECs in the case of $c/a=1.5$ obtains a characteristic local minimum (Figure 7). For the above two wave frequencies, the diffracted wave field in the seaward side of the wall in the absence of the WECs is also shown in Figure 15, which more clearly demonstrates the wave disturbances induced by the wall itself.

Starting with the case of $\omega =1.5$ rad/s, the presence of the wall boundary in the absence of the WECs (Figure 15a), at approximately $2.5\le Y/a\le 4.5$, leads to the formation of a longitudinal zone, with a sinusoidal-like variation pattern along $X$ (coined hereafter “zone A1”), where ${\eta }_D/A$ has values almost equal to zero. Outside of this zone, ${\eta }_D/A$ varies quite intensively and has significant values. By placing the WECs in front of the wall (Figure 13), zone A1 still exists; depending, however, upon the examined $c/a$ value, the zone’s width along $Y$ and its variation pattern along $X$ are modified due to the hydrodynamic interactions among the WECs and between the wall and the WECs. Accordingly, from a physical point of view, the case of $c/a=1.5$ (Figure 13a) corresponds to the placement of WECs outside of zone A1. More specifically, the WECs are placed in an area between the leeward side of zone A1 and the wall, where, the combined diffraction disturbances induced by the WECs and the wall amplify the intense variation pattern and the large values of ${\eta }_D/A$ already existing, due to the presence of the wall. As a result, a non-symmetric (with respect to the WECs’ local horizontal axes) diffracted wave field, with significant ${\eta }_D/A$ values, is formed around the WECs. On the other hand, the case of $c/a=3.0$ (Figure 13b) is physically equivalent to the placement of the WECs inside zone A1. This in turn, leads to the formation of an almost symmetric (with respect to the WECs’ local horizontal axes) diffracted wave field around the WECs, with quite small ${\eta }_D/A$ values. All the above advocate the existence of larger $F_3$ values for $c/a=1.5$ at $\omega =1.5$ rad/s as compared with the case of $c/a=3.0$, in absolute accordance with the results of Figure 7.

In the case of $\omega =2.3$ rad/s, the consideration of incident waves with shorter wave length results in the formation of three distinct longitudinal zones that are almost parallel to the wall, at the seaward side of the boundary (Figure 15b), where ${\eta }_D/A$ has almost zero values. The first zone (coined hereafter “zone B1”) is bounded at approximately $1.0\le Y/a\le 2.0$, while the second and the third zones (coined hereafter “zone B2” and “zone B3”, respectively) are bounded at approximately $4.0\le Y/a\le 5.0$ and $7.0\le Y/a\le 8.0$, respectively. These zones remain almost unaffected, when the WECs are placed close to the wall, i.e., in the case of $c/a=1.5$ (Figure 14a), whereas for $c/a=3.0$ the existence of the WECs modifies the width and the variation pattern of these zones along $X$ (Figure 14b). Following a rationale similar to the case of Figure 13, from a physical point of view, the non-dimensional distance from the wall $c/a=1.5$ (Figure 14a) corresponds to the placement of WECs inside zone B1. This, in turn, results in the existence of an almost symmetric (with respect to the WECs’ local horizontal axes) diffracted wave field around the WECs, with quite small ${\eta }_D/A$ values. However, the case of $c/a=3.0$ (Figure 14b) is physically equivalent with the placement of the WECs in the area between zones B1 and B2, where intense variation patterns and significant values of ${\eta }_D/A$ are observed. Consequently, a non-symmetric (with respect to the WECs’ local horizontal axes) diffracted wave field with significant ${\eta }_D/A$ values is formed around the WECs. All the above advocate for the existence of very small $F_3$ values for $c/a=1.5$ at $\omega =2.3$ rad/s as compared with $c/a=3.0$, in accordance with the results of Figure 7.

4. Conclusions

In this paper, we investigate in the frequency domain the performance (hydrodynamic behavior and power absorption) of a linear array of five semi-immersed oblate spheroidal heaving WECs placed in front of a bottom-mounted, vertical wall of finite-length and of negligible thickness under the action of perpendicular to the wall regular waves. A summary of the main conclusions of this investigation is the following:

• Compared to the cylindrical WECs, the deployment of either the hemisphere-shaped or the oblate spheroidal WECs improves the power absorption ability of the array to a great extent. However, among all the three WECs’ geometries examined, the oblate spheroidal WECs are considered to have the best power absorption ability due to their intrinsic hydrodynamic characteristics that also enable the consideration of larger $b_{PTO}$ values.
• Irrespectively of the values of $l_w/a$ and/or $c/a$, the presence of the wall boundary in the leeward side of the oblate spheroidal WECs, contrary to the isolated array, leads to the existence of $F_3$ global maxima at $\omega \ne 0.01$ rad/s with values larger than 1.0.
• For different $c/a$ values, the oblate spheroids’ heave exciting forces show great differences at specific frequency ranges. From a physical point of view, this fact can be related to the realization or not of an array’s placement at existing characteristic wave field zones, resulting from the presence of the wall, where the relevant wave elevation has values almost equal to zero. In cases, where $c/a$ falls within these zones, a smooth diffracted wave field around the WECs with symmetrical features is formed advocating the existence of very small $F_3$ values, while the opposite holds true when the WECs are placed outside of these zones.
• For the smallest examined distance from the wall, the power absorption ability of the oblate spheroidal WECs array is not driven by resonance phenomena, as large heave exciting forces, significant responses, and thus maximum $p\left(\omega \right)$ values, occur at wave frequencies outside the range, where WECs’ resonance occurs. The opposite holds true for the remaining examined$c/a$ values, where the presence of the wall does not impose any restrictions on the $RA{O}_3$ amplification due to resonance.
• The placement of the oblate spheroidal WECs array at successively larger distances from the wall induces hydrodynamic interactions between the spheroids and the boundary that enhance consecutively the hydrodynamic behavior of the WECs, and thus the array’s power absorption ability. However, this holds true for $c/a$ up to $3.0,$ since a further increase of $c/a$ leads to arrays with consecutively reduced power absorption ability. Consequently, among the $c/a$ values examined in this paper, $c/a=3.0$ presents an upper limit of this design parameter in terms of power absorption enhancement. Compared to the isolated array, the presence of the wall boundary positively affects the power absorption ability of the array leading to a significant increase of $p\left(\omega \right)$ at specific frequency ranges, as well as to more than one $p\left(\omega \right)$ peak, depending upon the value of $c/a$.
• The increase of $l_w/a$ from 18 to $36$ enhances the power absorption ability of the oblate spheroidal WECs array to a small extent, and a further increase of $l_w/a$ does not lead to any relevant improvement. Compared to the isolated array, the existence of the wall boundary, irrespectively of its length, affects positively the power absorption ability of the array, since it results in a significant increase of $p\left(\omega \right)$, especially at the frequency range where resonance phenomena occur.

This investigation can be further extended in order to examine the effect of other characteristic parameters (e.g., PTO characteristics, wave directionality, etc.) on the performance of the oblate spheroidal WECs array under the action of both regular and irregular waves. Moreover, in terms of power absorption maximization, the hydrodynamic model can be appropriately coupled with an optimization algorithm aiming at defining the optimum distance or location of the WECs in front of the wall.