Performance of an Array of Oscillating Water Column Devices in Front of a Fixed Vertical Breakwater

: The present study explores the performance of an array of cylindrical oscillating water column (OWC) devices, having a vertical symmetry axis, placed in front of a bottom seated, surface piercing, vertical breakwater. The main goal of this study is the investigation of a possible increase in the power efficiency of an OWC array by applying, in the vicinity of the array, a barrier to the wave propagation, aiming at amplifying the scattered and reflected waves originating from the presence of the devices and the wall. To cope with the set goal, a theoretical analysis is presented in the framework of linear potential theory, based on the solution of the proper diffraction, and pressure-radiation problems in the frequency domain, using the image theory, the matched axisymmetric eigenfunction expansion formulation, and the multiple scattering approach. Numerical results are presented and discussed in terms of the expected power absorption by the OWCs comparing different array’s characteristics i.e.,: (a) angle of incidence of the incoming wave train; (b) distances from the breakwater; and (c) geometric characteristics of the different arrangements. The results show that compared to the isolated OWC array (i.e., no presence of the wall), the power efficiency of the OWCs in front of a breakwater is amplified at specific frequency ranges.


Introduction
Sea waves have enormous power, therefore the construction of structures for mitigating such power is not easily accomplished. Breakwaters are widely used in coastal and offshore engineering, frequently applied in coastal protection and restoration schemes. They are barriers, either watersurface piercing or submerged, that are frequently displaced perpendicularly to the dominant direction of the incoming waves, which absorb, diffract, and reflect part of the wave energy, reducing the amount of energy that reaches the shoreline.
Breakwaters are primarily classified according to their structural features as restrained to the wave impact (fixed) or floating structures, whereas under these broad classifications, they can be further subdivided regarding their construction materials, shape, etc. [1]. Specifically, concerning the fixed breakwaters, they operate by reflecting the incoming wave train, as a bottom mounted rigid structure. Rubble mound, seawall, and barrier types of breakwaters fall in this category. Rubble mound breakwaters have probably existed for around 3000 years [2] and they are still applied for sheltering coastlines from wave action. Construction of breakwaters during the ancient time around the Mediterranean Sea was done by blocky stones, sometimes cementitious infill [3], while in recent times, numerous design methods for hydraulic performance and structural stability of rubble mound breakwaters have been developed i.e., [4][5][6][7]. A seawall breakwater is a most common defensive structure, acting as one large solid block at locations where the ocean environments are dominating at the coast. It is usually composed of prefabricated reinforced concrete caissons, representing a better alternative in terms of performance, construction rapidity, standardization, environmental implications, and construction and maintenance costs when compared with the rubble mound type [8]. Many analytical and laboratory studies and field observations have been undertaken concerning the design and construction of vertical breakwaters (seawalls). Indicative, recent studies are [9][10][11]. On the other hand, in situations where the complete protection from the waves is not required, thin barriers, impermeable or permeable, supported with piles can be used as breakwaters. Pile breakwaters comprise a series of piles which partially attenuate the wave energy due to turbulences and eddies created around the solids, preventing also effectively the soar from sediment siltation [12][13][14]. In addition, barrier type of breakwater consists of perforated and slotted structures, operating as permeable breakwaters, increasing wave reflection as their porosity increases [15,16].
Although a breakwater structure is constructed to minimize the wave action in areas behind of it, in front of the structure the incoming wave energy is amplified due to the scattered and reflected waves originating from the presence of the vertical wall. This phenomenon has triggered increased interest on wave energy conversion systems operating near and/or on a breakwater. Moreover, the installation of wave energy converters (WEC) in front and/or on breakwaters is facilitated by easier electricity transmission to the mainland, allowing for common usage of infrastructure (i.e., electrical cable, power transfer equipment, etc.).
In the context of a breakwater-WEC system, numerous projects and studies have been presented globally, emphasizing the wall's positive effect on the converter's efficiency. Although several different types of WECs are under development, only few typologies have been commonly used in conjunction with coastal protection structures, namely (a) the overtopping devices (OTD); (b) the oscillating wave surge converter (OWSG); (c) the point absorbers (i.e., heaving devices); and (d) the oscillating water column devices (OWC).
Rubble mound breakwater-OTDs utilize a frontal sloping plate that leads the incident waves to overtop into one or more storage basins placed at a higher level than the seawater level. In its natural way back to the sea, the water passes through turbines, generating electricity [17]. Indicative studies concerning the design optimization of an OTD for efficiency maximization are [18][19][20][21]. Regarding the OWSC operation, the converter typically has one end fixed to the sea bed while the other end is free to move. The wave energy is absorbed by the relative motion of the body (i.e., the converter comes in the form of floats, flaps, or membranes) compared to the fixed point. Analytical formulations on the effect of an onshore barrier (i.e., a straight coast and a vertical breakwater) on the device's efficiency are presented in [22,23].
As far as point absorber devices are concerned, they gain the energy from their oscillation in heave direction, which is driven by means of their interaction with the wave field. The arrangement of a breakwater-heaving device system involves partial-floating bodies placed parallel to the predominant wave direction, in front of a vertical seawall. The bodies' movements due to the scattered and radiated waves from the array's members and their interaction with the vertical wall, allows the conversion of the floaters' kinetic energy to electricity using hydraulic or mechanical transmission [24]. Most recent indicative studies concerning hydrodynamic analysis and efficiency estimation of arrays of heaving devices placed in front of a linear vertical wall are presented in [25,26].
An OWC device is a partially submerged, hollow structure open to the seabed below the water line. The vertical motion of the sea surface alternatively pressurizes and depressurizes the air inside the OWC's chambers, generating a reciprocating flow through a self-rectifying turbine which is installed beneath the roof of the device. Looking towards the multi abilities provided by the breakwater-OWC system, several different designs have been presented in the literature, while most of them concern OWCs integrated at a vertical breakwater. Specifically, in [27,28], a theoretical study of an OWC standing at the tip of a breakwater and along a straight coast is presented, respectively, whereas in [29][30][31], an OWC integrated at a flat breakwater is theoretically and experimentally investigated. In [32,33], a detailed analysis concerning the structural and economic feasibility of integrated OWC within a Mediterranean port is presented and in [34], a linearized theory of an array OWC installed on a straight coast is described. In addition, the effect of a breakwater on OWC performance using CFD analysis under the action of regular and irregular waves is presented in [35]. Finally, in [36], a modified integrated breakwater-OWC system is investigated using numerical and experimental simulations in terms of its power performance. As far as the performance of an OWC device placed in front of a vertical breakwater is concerned, only few studies are presented in the literature. In [37], a theoretical analysis of a vertical OWC device placed in front of a vertical wall is presented, whereas in [38], the efficiency of an array of five OWCs placed in parallel direction to a vertical breakwater is examined, for installation at the port of Heraklion (i.e., Crete island).
The scope of this work is to examine the effect of a bottom mounted, surface piercing, breakwater of infinite length, on the power efficiency of an array of OWCs placed in a random location in front of the wall. The examined converters consist of an exterior partially immersed toroidal body supplemented by a coaxial interior, bottom mounted, free-surface piercing vertical cylinder. In the annulus between the internal cylinder and the external torus, a finite volume air chamber is formed in which the oscillating air pressure is developed. An analogous OWC type has been examined to operate alone in the open sea, combined with a wind turbine supported on the converter's interior concentric cylindrical body, and as part of a Jacket platform [39][40][41][42]. A theoretical model is presented, taking into account the wave hydrodynamic interactions among the OWCs and the fluid flow in front of the breakwater. Several distances between the bodies and the breakwater are examined, along with different array configurations (i.e., OWCs in a rectangular, parallel, and perpendicular arrangement to the front wall) and wave heading angles, to assess the array's efficiency towards its optimization. The presented results show that the power efficiency of an array of OWCs in front of a breakwater is amplified compared to the one in unbounded waters (i.e., without the presence of the vertical wall).
This work is organized as follows: Section 2 describes the solution of the corresponding diffraction and pressure-radiation problems, while in Section 3, the OWCs' hydrodynamic characteristics and absorbed power are presented. Section 4 provides and discusses the numerical results and finally the conclusions are drawn in Section 5.

Formulation of the Hydrodynamic Problem
The diffraction and the pressure-radiation problems under consideration are examined within the context of the arrangement shown in Figure 1. An array of N similar OWCs is assumed, situated in the vicinity of a vertical, bottom-mounted, surface piercing, breakwater of infinite length. The water depth is denoted by h, assuming the sea bottom flat and horizontal. The outer and inner radii of each device's chamber are denoted by α, b, respectively, whereas the distance between the bottom of the external torus and the seabed is denoted by hc. The radius of the interior, bottom seated, coaxial cylindrical body is denoted by c. In addition, the distance between the center of the closest to the wall converter and the breakwater is denoted by Lw, whereas the distance between adjacent OWCs is by Lb. Small amplitude harmonic waves (with angular frequency ω, wave height H, and wave length λ) are incident to the breakwater at an angle β. A global, right-handed Cartesian co-ordinate system Oxyz is introduced with origin located at the bottom of the breakwater, with its vertical axis Oz directed upwards, while N local cylindrical co-ordinate systems , , , q = 1, 2, …, N are defined with origins at the intersection , of the sea bottom with the vertical axis of symmetry of each converter. The geometric layout of the breakwater-OWC system is illustrated in Figure 1. In the present analysis, the fluid is assumed non viscous and incompressible and the water flow irrotational so that linear water wave theory may be employed. Under this assumption, the fluid flow around each device q, q = 1, 2, …, N, described by the potential function: , , ; = , , can be decomposed on the basis of linear modeling as [43]: , , = , , + , , In Equation (1), stands for the velocity potential of the undisturbed incident harmonic wave; is the scattered potential around the q OWC, assuming atmospheric air pressure inside the chamber; denotes the pressure-dependent radiation potential around the q OWC when it is considered as an open-duct body (i.e., atmospheric air pressure inside the chamber) due to unit air pressure , , in the chamber of the p device. The term = + denotes the diffracted component of the corresponding total wave potential around the q body.
The potentials , , are solutions of the Laplace equation in the entire fluid domain and satisfy the proper boundary conditions on the sea bed and the water free surface; the kinematic conditions on the mean body's wetted surface and the no-flux boundary condition on the breakwater's surface [44]: In the present study, the method of images is applied together with the assumption of a fully reflecting breakwater. Specifically, the presence of the breakwater in front of the OWC array is represented by taking into consideration the image "virtual" devices with respect to the wall, without the presence of the wall. The equivalent array of 2N devices is exposed to the action of two-directional surface waves (i.e., one propagating at angle β and a second at angle 180-β) [23,37] (see Figure 2). Herein, based on the method of images, 2N local cylindrical co-ordinate systems , , , = 1, … ,2 are defined with origins at the intersection , of the sea bottom with the vertical axis of symmetry of each OWC. Thus, Equation (1) can be written as: In Equation (3), , , correspond to the velocity potential of the undisturbed incident harmonic wave; the scattered potential around the q, OWC, q = 1, 2, …, 2N, and the pressuredependent radiation potential around the q OWC q = 1, 2, …, 2N, respectively. In addition, , corresponds to the air pressure in the chamber of the p device, p = 1, …, 2N.
(c) kinematic condition on the mean device's wetted surface: Finally, a radiation condition stating that propagating disturbances must be outgoing is imposed.
In Equation (4), is the gravitational acceleration; is the water density; , is the Kronecker's symbol. In addition, in Equation (6), the () term denotes the derivative in the direction of the outward unit normal vector , to the mean wetted surface on the q OWC. The wave interaction phenomena among the OWCs (i.e., initial and image solids) have been taken into consideration through the physical idea of multiple scattering. Specifically, by properly superposing the incident wave potential and the propagating and evanescent modes that are scattered and radiated by the OWC's, exact representations of the fluid's velocity potentials around each device can be obtained, based on the single device hydrodynamic characteristics. The latter are derived through the use of matched axisymmetric eigenfunction expansions [45] for the velocity potential around each single OWC, considered alone in the wave field. Based on this method, the flow field around each OWC is subdivided in coaxial ring-shaped fluid regions, defined as I, II, III, (see Figure 1b) in each of which appropriate series representations of velocity potential can be established. These series representations are solutions of Laplace equation; satisfy the Equations (4)-(6); the radiation condition at infinity and the continuity relation, of the velocity potentials and their radial derivatives, at the vertical boundaries of neighborhood fluid regions.
The method for the solution of the diffraction and the pressure-radiation problems of a single OWC device along with the implementation of the multiple scattering approach on array of OWCs has been thoroughly described in [46]. Nevertheless, by the way of example, the velocity potentials around an isolated OWC device are presented in the Appendix A.

Array's Efficiency
Due to the water oscillation inside the OWCs' chambers, the dry air above the free surface is pushed through an air turbine, located at the top of the chamber, producing an air volume flow inside each chamber, denoted as , , ; = q , , . Here, q equals to: The term denotes the vertical velocity of the water surface in the q OWC, q = 1, .., N; is the cross-sectional area of the inner water surface inside the q device; and is the velocity potential inside the OWCs' chambers (i.e., III fluid domain).
Since the velocity potential around each device of the array can be described as a superposition of the diffraction and the pressure-radiation velocity potentials (see Equation (1)), similarly the air volume flow q inside the q OWC, q = 1, .., N, can be written in the form: In Equation (8), the term q stands for the diffraction air volume flow, while q denotes the pressure-dependent volume flow, known as radiation admittance [47].
The radiation admittance of the q OWC q = 1, …, N due to unit air pressure head inside the p device, p = 1, …, N can be also written as a function of the radiation conductance and the radiation susceptance coefficients, i.e.,: Based on the method of images, the considered breakwater-OWC system interacting with an incident wave of angle β, is equivalent to an array of 2N devices, mirror between each other with respect to the breakwater, which are exposed to the action of two wave trains at angles β and 180−β, without the presence of the vertical wall. Thus, the diffraction volume flow inside the q OWC q = 1, …, N, equals to the sum of the corresponding diffraction volume flows inside the q device for wave angle propagation β and 180-β, denoted as , and , , respectively, when the device is assumed part of an array of 2N bodies, i.e.,: Here, the terms , , = β, 180 − β are derived from Equation (7) for every examined wave heading angle.
The radiation admittance, q , of the q OWC device, q = 1, 2, ..., N, placed in front of a vertical wall can be also derived through the method of images. More specific, the radiation conductance and the radiation susceptance coefficients of the q OWC, q = 1, …, N due to unit air pressure head inside the p OWC, p = 1, …, N can be derived by summing properly the radiation conductance and susceptance coefficients of the q OWC q = 1, ..., N, due to unit air pressure head inside the p OWC, p = 1, …, N with the corresponding coefficients of the q OWC q = 1, ..., N, due to unit air pressure head inside the p OWC's image device, denoted as p', p' = N+1, …, 2N [37,38].
In the present paper, an air turbine is assumed to be placed in the OWCs' ducts, between the chamber and the outer atmosphere that exhibits an approximately linear relationship between the inner air pressure and the volume flow (e.g., a Wells type air turbine) i.e.,: q = Λ (11) Here, Λ represents the complex pneumatic admittance of the air turbine [48]. The real part of Λ is related to the pressure drop through the turbine, whereas the imaginary part represents the effect of air compressibility inside the chamber. In the presented numerical results (see Section 4), the pneumatic admittance of the air turbine is assumed to attain an optimum value, Λ , as presented in [49], which maximizes the power efficiency of a similar OWC device, considered alone in the wave field, without the presence of the breakwater.
Having determined all the pressure coefficients of an array of N OWCs in front of a vertical breakwater, each devices' wave absorbed power, , = 1, … , can be written as: Here, ω stands for the wave frequency. In Equation (12), the term denotes the efficiency of the q OWC device since in the present paper, the losses that occur in the energy conversion chain are not taken into consideration.
To evaluate the constructive or destructive effect of a vertical seawall on the WECs' power efficiency, a "q-factor", , is introduced which is defined as the ratio of the total wave power absorbed by the OWC array to N (the number of OWCs in the array) times the absorbed power by the same converter in isolation condition [26,50]. Thus, it holds: For > 1, the scattered and reflected waves due to the presence of the solids and the vertical wall have a constructive effect on the array's efficiency, increasing the value of the total absorbed wave power compared to the corresponding efficiency of N OWCs placed in isolation (alone in the wave field). On the other hand, when < 1, the wave interaction phenomena between the devices and the breakwater have a destructive effect on the amount of the array's efficiency compared to the wave absorbed power by N number of isolated OWCs.

Test Cases
The theoretical method described in the present paper is applied to an array of OWCs placed in front of a breakwater. The examined converter's external and inner radii equal to α and b = 0.9α, respectively. The distance between the bottom of the external torus and the seabed equals to hc and the water depth equals to h = 7.14α. Concerning the inner cylindrical body, it is assumed to be sea bottom seated with a radius of c = 0.4α The OWCs' air turbine characteristics are considered to be equal to the Λopt value of a similar OWC in isolation condition at its pumping resonance wave frequency [51], whereas the distance between the center of the closest to the wall converter and the breakwater is Lw, and the distance between adjacent OWCs is Lb = 4α (see Figures 1 and 2).
The calculation of the Fourier coefficients of the velocity potentials (see Appendix A) around each OWC device in the array is the most significant part of the presented theoretical analysis because of their influence on the accuracy solution. Here, for the I and III ring elements n = i = 60, whereas for the II fluid domain n = 100. The considered number of modes are m = ±7 and the number of wave interactions were taken equal to 7. The presented numerical results were obtained using the in-house developed computer software HAMVAB [52]. The software being relied on analytical representations of the velocity potential around each cylinder-type OWC device of the array was preferred in the present contribution against other available numerical tools applicable to general 3D geometries, since by keeping the same accuracy with them, it is usually less CPU time-consuming (i.e., for the solution of the diffraction and radiation problem in each wave frequency a time of 40 s is required).
In the following subsections, the absorbed power by the OWC array is presented for several examined parameters, namely: (a) number of converters; (b) array orientation to the breakwater; (c) wave heading angles; (d) distances between the device and the breakwater; (e) device's draught; and (f) distances between the devices.

Effect of the OWCs Orientation to the Breakwater
In this subsection, the effect of the OWC-array orientation with respect to the breakwater on the array's efficiency is presented. The examined converters are placed in front of a vertical breakwater in: (a) parallel to the wall arrangement; (b) perpendicular to the wall arrangement; and (c) in a rectangular arrangement in front of the wall. Furthermore, in each of these arrangements, several numbers of OWC devices are considered, i.e., configurations 1, 2, 3, 4, 5. In Figure 3, the examined orientation of the OWC's arrays arrangements and the breakwater are depicted. Herein, the distance between the center of the closest to the wall converter and the breakwater equals to Lw = 3α, whereas the distance between the bottom of the external torus and the seabed equals to hc = 6.14α.  According to the presented analysis, the pumping resonance of the water column inside the oscillating chamber occurs at wave frequency equal to 2.62 [rad/s], for α = 1 m, thus the air turbine coefficient inside the OWCs equals to 10.60 [m 5 /(kN.s)]. Figure 4 depicts the modulus of the inner air pressure head inside the 1st OWC of the array, i.e., /( /2), q = 1, normalized by the wave amplitude along with the factor (as described in Equation (13)) for each examined OWC configuration for the parallel to the wall array arrangement. The results are plotted against the non-dimensional wave numbers, kα, in the range of kα ∈ [0.05, 1.5] and for wave heading angles β = 0, π/6, π/4 (i.e., here k stands for the wave number). Similarly, in Figures 5 and 6, the corresponding inner air pressure of the 1st OWC device and the factor for the perpendicular and rectangular to the wall array arrangements, respectively, are presented.   It can be seen from Figure 4 that the breakwater has a significant effect on the inner air pressure head of the 1st OWC device. The values of the air pressure inside the OWC in front of a breakwater do not follow a similar variation pattern to the corresponding values for the unbounded water case. Specifically, the values of the air pressure inside the OWC when placed in front of a breakwater attain almost two times larger values than those of the air pressure inside the same OWC, without, however, the presence of the breakwater, for kα tending to zero, at all the examined wave heading angles and OWC configurations (i.e., 1, …, 5 number of bodies in front of the breakwater). Furthermore, it can be seen that the wave interaction phenomena between the OWCs and the vertical wall affect the values of the inner air pressure of the 1st converter also at higher values of kα. It is depicted that the inner air pressure tends to zero for kα values in the range of [0.5, 0.6] for β = 0; π/6; π/4. This behavior does not appear in the case of the isolated OWC without the presence of the wall. The zeroing of the air pressure head can be attributed to the standing wave due to the presence of the vertical breakwater and in particular it appears when the distance between the initial and the image converter equals to half wave length [53]. Furthermore, it is also notable that for kα ≈ 0.7, the inner air pressure of the converter, when it is considered part of an array or isolated in the wave field, exhibits a resonant peak regardless the wave heading angle. This resonance is associated with the pumping modes of the interior basin due to the existence of the moonpool [54]. As far as the number of the OWCs in each examined configuration is concerned, it can be seen that as the number of the OWCs in the array increases, the values of the inner air pressure inside the 1st OWC device attain an oscillatory behavior around the corresponding values of the single breakwater-OWC case, for large values of kα.
Concerning the effect of the breakwater on the array's power efficiency, it can be seen from the factor figures, that this effect is constructive or destructive depends on the examined wave number. Specifically, at small wave numbers, the efficiency of the array is four times higher than the absorbed power by an isolated OWC (without the presence of the vertical wall). On the other hand, at kα ∈ [0.5, 0.6], the factor attains values lower than one, thus the breakwater has a destructive effect on the array's power efficiency. Furthermore, at higher wave numbers, i.e., kα > 0.6 the presence of the breakwater increases the array's absorbed power compared to the isolated case, since > 1. Nevertheless, as the wave heading angle increases, the wave numbers in which > 1 are sifted to higher values.
As far as the perpendicular arrangement of the OWCs in front of the breakwater is concerned, a similar conclusion to Figure 4 can be derived concerning the effect of the vertical wall on the values of the inner air pressure of the 1st OWC (i.e., the one closest to the breakwater). Specifically, as in the parallel arrangement case, for kα tending to zero, the doubling of the air pressure values for the breakwater-OWC case is also notable, compared to the no-wall case. In addition, the zeroing of the values at the wave number, which corresponds to a wavelength, equals to two times the distance between the initial and the image converter is presented here. Furthermore, the resonances at the pumping wave frequency, i.e., kα ≈ 0.7 are also depicted regardless of the number of the converters in the array. On the other hand, the factor graphs do not follow the same pattern as the corresponding graphs in the parallel arrangement. It can be reported that the values of the oscillate around the corresponding values of the single OWC-breakwater case. These oscillations become more pronounced as the number of the OWCs in the array increases and can be attributed explicitly to the scattered and reflected waves that seem to strongly affect the wave field around the converters with respect to the incident wave direction. Concerning the examined wave heading angles, it can be seen that the constructive effect of the breakwater on the array's efficiency, compared to the isolated OWC case, is increased as the values of the wave heading angle and the wave number increase. Furthermore, compared to the parallel OWC arrangement, here the values do not zero at kα ∈ [0.5, 0.6], since the distances between the OWCs and the breakwater are not the same for each OWC (i.e., the distance between the initial and image converters does not correspond to half wave length for each OWC at this wave number band). As far as the remaining examined kα values are concerned, the does not seem to attain significantly higher values than the corresponding factor of the parallel arrangement (see Figure 4), despite that the wave interactions between the bodies and the breakwater are amplified in this arrangement. Hence, it can be concluded that the increased interaction phenomena between the OWCs and the breakwater are not always beneficial for the array's efficiency at every wave frequency.
Continuing with the results of the rectangular arrangement, presented in Figure 6, it can be noted that the graphs of the air pressure inside the 1st OWC device follow a similar pattern as in the aforementioned arrangements (i.e., parallel and perpendicular). It is depicted that due to the presence of the breakwater, the air pressure attains double values than the no-wall case, at wave numbers tending to zero, regardless of the wave heading angle and the number of bodies in the array. Additionally, the resonance of the inner air pressure, at the pumping wave frequency i.e., kα ≈ 0.7, is notable. However, this resonance is dictated by the wave heading angle and the wave interaction phenomena between the converters and the vertical wall, since as the wave angle increases, the air pressure values at kα ≈ 0.7 decrease. Furthermore, the zeroing of the air pressure values at kα ∈ [0.5, 0.6] is also notable in the rectangular arrangement. This was also the case for the parallel and the perpendicular examined arrangements. Thus, it can be obtained that the wave number in which the OWC's inner air pressure zeros is affected by the device's distance from its image device and not by the examined array configuration and orientation to the wave impact. As far as the array's efficiency is concerned, it can be derived that for small wave numbers, the breakwater has a constructive effect on the array's power performance regardless of the examined body configuration, or array arrangement, attaining four times higher values of absorbed power than the no-wall cases. On the other hand, as the wave number increases, the distances between the devices and the vertical wall, as well as the wave heading angles and the number of bodies in the array do not always have a constructive effect on the array's power efficiency.
In conclusion, when comparing the results from Figures 4-6, it can be obtained that for wave numbers tending to zero, the three examined configurations attain similar results. For wave numbers in the range of kα ∈ [0.1, 0.3], the parallel to the breakwater OWC arrangement performs better than the rectangular or perpendicular arrangements, attaining higher values of . On the other hand, the breakwater has a destructive effect on the efficiency of the array arranged parallel to the vertical wall for wave numbers in the range of kα ∈ [0.5, 0.6], whereas for the rectangular and perpendicular array, the presence of the vertical wall has a positive impact on the array's power performance at the same wave number range. At wave numbers around the pumping resonance wave frequency, i.e., kα ≈ 0.7, the three examined arrangements attain similar results. Nevertheless, the array's efficiency is increased as the number of bodies in the array increases. Finally, at large wave numbers i.e., kα > 1, the array's power efficiency is mainly affected by the wave heading angle, i.e., the values of increase as the wave heading angle increases.

Effect of the OWCs' Distance from the Breakwater
In this subsection, the effect of the distance between the OWCs' and the breakwater on their wave power efficiency is examined. Here, a parallel array arrangement of five similar OWCs in front of a breakwater is considered (see Configuration 5 in Figure 3a). The characteristics of the examined OWCs are presented in Section 4.1, whereas here, hc = 6.14α. The examined distances from the wall equal to Lw = 3α, 6α, 9α, 12α, and the wave heading angles to β = 0, π/6, π/4. Additionally, the air turbine coefficient inside the OWCs equals to 10.60 [m 5 /(kN.s)]. Figure 7 depicts the modulus of the air pressure head inside the 1st OWC of the array, as well as the factor for each examined distance between the devices and the breakwater and wave heading angles. It can be noted that the distance between the center of the devices and the breakwater significantly affects the air pressure head inside the 1st OWC device. Specifically, the values of the air pressure when the array is placed away from the wall (i.e., Lw = 6α, 9α, 12α) oscillate around those of the Lw = 3α case. It can be also seen that the larger the distance of the converters from the wall is, the stronger the oscillatory behavior of the air pressure RAOs. As far as the examined wave heading angles are concerned, it can be seen from the left side of Figure 7, that the inner air pressure graphs follow a similar pattern, regardless of the wave heading angle. Nevertheless, it should be mentioned that the wave numbers in which oscillations in the air pressure occur, are shifted to higher values as the wave heading angle increases. Concerning the efficiency of the array at several distances from the breakwater, it can be seen (right side of Figure 7) that the examined distances between the devices and the vertical wall cause the array's absorbed power to increase at some kα values and to decrease at another range of wave number compared to the efficiency of five OWCs in isolation condition. This is due to the reflected waves from the wall, which cause a dominant oscillating behavior on the array's absorbed power values. Moreover, as the distance of the OWC from the vertical wall increases, the earlier the oscillatory behavior of the wave absorbed power occurs. As far as the incident wave angles are concerned, it can be obtained that the wave numbers in which oscillations in the power efficiency occur are shifted to higher values as β increases. Finally, it can concluded that the efficiency results for the close to wall distances (i.e., Lw = 3α, 6α) follow a more steady pattern, at every examined wave number, compared to the results of the large-distance cases (i.e., Lw = 9α, 12α), of which they attain a large oscillatory behavior.

Effect of the OWCs' Draught
Next, the effect of the OWCs' draught on their wave power efficiency is examined. Here, also, a parallel arrangement of an array of five similar OWCs in front of a breakwater is considered (see Configuration 5 in Figure 3a). The characteristics of the examined OWCs are presented in Section 4.1, whereas Lw = 3α. The examined OWCs' draughts equal to hc = 6.14α, 5.64α, 5.14α, 4.64α, and the wave heading angles to β = 0, π/6, π/4. Due to the variation of the converter's draught, and consequently to its pumping resonance wave frequency value the air turbine characteristics for the OWC devices equal to: Λ = 10.6; 9.04; 10.47; 18.49 [m 5 /(kN.s)] for each examined draught hc = 6.14α, 5.64α, 5.14α, 4.64α, respectively, for α = 1 m. Figure 8 depicts the inner air pressure head of the 1st OWC device (left column) and the array's total absorbed power (as presented in Equation (12)), i.e., P(ω)/(H/2) 2 , (right column) for various examined OWC draughts and wave heading angles. From the depicted results, it can be seen that the draught of the OWCs affects the inner air pressure head as well as the array's efficiency. Specifically, as the draught of the converter increases, the pumping resonance wave frequency is shifted at lower values, thus the corresponding wave numbers in which the inner air pressure values resonance are transferred to lower values too. Furthermore, it can be noted that as the draught of the OWC increases, the air pressure head inside the 1st OWC decreases. As far as the examined wave heading angles is concerned, it is depicted that β has a minor effect on the inner air pressure head for kα < 0.5. On the other hand, for kα > 0.5, the inner air pressure head inside the 1st OWC decreases as the wave heading angle increases. Concerning the absorbed wave power, it can be seen from the results that the wave numbers in which the absorbed power maximizes (sharp peaks) are not remained constant regardless of the floater's draught, since the wave frequencies at pumping resonance are shifted to lower values. Furthermore, it is noted that the increase of the OWCs' draught has a destructive effect on the array's efficiency as it results in a significant decrease of P(ω). Concluding, it can be obtained that the increase of the OWCs' draught does not increase the array's efficiency. Contrary, as the draught decreases, the array's absorbed power increases at various wave numbers.

Effect of the Distance between the OWCs
This subsection examines the effect of the distance between the devices on their wave power efficiency. Here, a parallel arrangement of an array of five similar OWCs in front of a breakwater is considered (see Configuration 5 in Figure 3a). The characteristics of the examined OWCs are presented in Section 4.1. Additionally, the distance between the center of the closest to the wall converter and the breakwater equals to Lw = 3α, whereas the distance between the bottom of the external torus and the seabed equals to hc = 6.14α. The air turbine coefficient inside the OWCs equals to 10.60 [m 5 /(kN.s)]. Several distances between the converters of the array are examined, i.e., Lb = 4α; 8α; 16α; 32α. Figure 9 depicts the modulus of the air pressure head inside the 1st OWC of the array, as well as the factor for each of the examined distance cases for zero wave heading angle. It can be seen that the inner air pressure of the 1st device seems not to be affected by the distance between the OWCs since similar results are depicted for each examined distance. Furthermore, the presented results follow the same variation pattern regardless the distance between the OWCs in most of the examined wave numbers. It should be noted that the values of the air pressure inside the 1st OWC for the Lb = 32α case are similar to the corresponding values of an isolated OWC in front of a breakwater (see Figure 4a). On the other hand, it can be seen that the qf factor is significantly affected by the distance between the devices, especially at higher wave frequencies, in which the qf values appear as an oscillatory behavior (i.e., kα > 0.7). This behavior becomes more pronounced as the distance between the adjacent OWCs increases. Contrary, for kα < 0.7, the efficiency of the array seems not to be affected by the distance between the devices.

Conclusions
In this study, the efficiency of an OWC array placed in front of a bottom-mounted, surface piercing, fully reflecting, vertical breakwater is investigated in the frequency domain. The examined OWC-array consists of identical OWC devices containing an exterior partially immersed toroidal body supplemented by a coaxial interior, bottom mounted, free-surface piercing cylinder. In the annulus between the internal cylinder and the external torus, a finite volume air chamber is formed in which the oscillating air pressure is developed. A theoretical formulation based on the image method has been applied to simulate the effect of the vertical wall on the array's power absorption, whereas the wave interaction phenomena between the vertical wall and the converters have been taken into account using the multiple scattering approach.
Three different types of array configurations in front of the vertical wall have been studied, namely a parallel, a perpendicular, and a rectangular arrangement. Furthermore, five types of OWC's arrays configurations have been also examined (i.e., arrays consisting of 1, 2, …., 5 bodies) for various cases of distances between the devices and the breakwater, wave heading angles, devices' draughts, and distances between adjacent devices. Based on the theoretical computations shown and discussed in the dedicated sections, the main findings of the present research contribution concern the significant effect of the breakwater on the array's efficiency, which is amplified compared to the performance of isolated OWCs (without the presence of the vertical walls). The power performance amplification is strongly affected by the array's orientation to the breakwater; the number of the bodies that consist the array; the wave heading angle; the devices' distance from the breakwater and the devices' draught. Therefore, it can be derived that the installation of an OWC array in front of a vertical breakwater can be an effective way to improve its power absorption efficiency.
However, the infinite wall assumption realized theoretically by the image method should be further examined concerning probable limitations attained at specific wave frequencies. Following the remarks from [55] concerning an array of cylindrical bodies in front of a vertical breakwater, discrepancies between the results of a finite-and an infinite-length breakwater occur. These discrepancies are more pronounced at wave numbers tending to zero, whereas at higher wave numbers, the results from both methods convergence.
The present research will be continued further by determining in the effect of a V-shaped breakwater of a random angle on the power efficiency of an array of OWCs and comparing it with the results of the present study.
Funding: This research received no external funding. Wave length rk,θk,zk Local co-ordinate system of the k OWC Φ Time harmonic complex velocity potential Velocity potential of the undisturbed incident harmonic wave Scattered velocity potential of the q OWC Diffraction velocity potential of the q OWC Radiation velocity potential resulting from the inner air pressure in p OWC Amplitude of the oscillating pressure head in the chamber of the p OWC g

Conflicts of
Gravitational acceleration ρ Water density Unit normal vector δq,p Kronecker's symbol Mean wetted surface of the q OWC I The infinite ring element around the q OWC II The ring element below the q OWC III The ring element inside the chamber of the q OWC ( ) Time dependent air volume flow Vertical velocity of the water surface in the OWC Cross-sectional area of the inner water surface inside the OWC q Diffraction volume flow of the q OWC q Pressure-dependent volume flow of the q OWC Radiation conductance of the q OWC Radiation susceptance of the q OWC Λ Complex pneumatic admittance (air turbine coefficient) Λ Air turbine coefficient optimum value ( ) Absorbed wave power by each OWC of the array q-factor term 1