A Theoretical Study of the Hydrodynamic Performance of an Asymmetric Fixed-Detached OWC Device

: The chamber conﬁguration of an asymmetric, ﬁxed-detached Oscillating Water Column (OWC) device was investigated theoretically to analyze its effects on hydrodynamic performance. Two-dimensional linear wave theory was used, and the solutions for the associated radiation and scattering boundary value problems (BVPs) were derived through the matched eigenfunction expansion method (EEM) and the boundary element method (BEM). The results for the hydrodynamic efﬁciency and other important hydrodynamic properties were computed and analyzed for various cases. Parameters, such as the length of the chamber and the thickness and submergence of the rear and front walls, were varied. The effects on device performance of adding a step under the OWC chamber and reﬂecting wall in the downstream region were also investigated. A good agreement between the analytical and numerical results was found. Thinner walls and low submergence of the chamber were seen to increase the efﬁciency bandwidth. The inclusion of a step slightly reduced the frequency at which resonance occurs, and when a downstream reﬂecting wall is included, the hydrodynamic efﬁciency is noticeably reduced at low frequencies due to the near trapped waves in the gap between the OWC device and the rigid vertical wall.


Introduction
Ocean wave energy is a renewable and pollution-free resource with the potential to mitigate the effects of global warming and contribute to meet the world's growing demand for electricity. This globally available energy source is estimated at about 2.1 TW [1], or 18,400 TWh per year, approximately 80% of the 2018 world's demand for electricity [2]. By use of wave energy converters (WECs), this energy source can be collected and transformed into electricity. In recent decades a broad range of WEC technologies has been proposed, with the Oscillating Water Column (OWC) device emerging as one of the most successful systems for wave energy harvesting.
The OWC differs from other technologies in its simple operating principle and adaptability for a range of locations (on the coastline, in the nearshore region or offshore). Thus, OWC systems have been researched most often, with full-scale prototypes being developed a number of times [3]. The system essentially consists of two components: a partially submerged hollow structure into which water enters through a bottom opening, and a Power Take-Off (PTO) mechanism, which converts the wave energy into electrical energy. Typically, the PTO system is a Wells turbine, which is the only moving element of the mechanism and which is above the level of the seawater to avoid direct exposure. The hollow chamber has a water column and an air column. The working principle is analogous to a piston in a cylinder; the trapped water inside is forced to oscillate up and down by the wave action. This alternately compresses and rarefies the air, trapped above the water column, and drives the airflow back and forth across the turbine in a duct that is coupled to a generator [4].
In an OWC system, wave energy is absorbed by a hydrodynamic mechanism involving relatively complex diffraction and radiation wave processes [3]. In recent years, many researchers have studied the hydrodynamic performance of the OWC device using theoretical analysis, experimental tests, and numerical simulations. Pioneering studies of wave energy absorbers based on rigid body models were performed by [5][6][7][8]. Earlier, studies of the parameters affecting the hydrodynamic efficiency of land fixed OWC devices were conducted by [9][10][11][12][13][14].
Pile-supported OWC structures, with hydrodynamic behavior similar to fixed-detached devices, have been widely studied as a means of extracting wave energy for electricity generation and reducing wave transmission, mainly through the mechanisms of wave reflection [15][16][17][18][19][20][21]. The performances of pile-supported and land-fixed OWC devices were analyzed in experiments by [15]. A pile-supported OWC-type breakwater was experimentally investigated by [16]. They found that greater pneumatic damping increased energy extraction, but that with less pneumatic damping energy dissipation was more efficient in vortex-induced energy loss. He and Huang [17] investigated the hydrodynamic performance of a pile-supported OWC structure in front of a vertical wall, demonstrating that this system can serve as a wave absorber to reduce the wave reflection from vertical walls. He et al. [18] employed the matched EEM to investigate the hydrodynamics of a pile-supported OWC breakwater. Their results showed that by adjusting PTO damping for maximum power, adequate power extraction and wave transmission is possible, whereas optimizing PTO damping for minimal wave transmission significantly lowers power extraction. The combined effect of the chamber geometry and wave conditions on the performance of fixed-detached and asymmetric OWC devices was investigated by [22]. Qu et al. [23] performed experimental tests to study an integrated OWC-pile-supported permeable breakwater and found that, compared with a non-permeable scheme, wave power generation and environmental protection performance were better.
Extensive investigations have been conducted on offshore fixed OWC system performance [24][25][26][27][28]. Ref. [24] performed both experimental test and numerical simulations of a 3D offshore stationary OWC device subjected to regular waves. They showed that the total power extraction efficiency of an offshore OWC is considerably overestimated by 2D and wave flume modelling, especially for wave frequencies greater than the chamber resonant frequency. Elhanafi et al. [25] investigated the effects of the underwater front and rear lips on the hydrodynamic performance of an offshore stationary OWC using a computational fluid dynamics (CFD) model. They concluded that by selecting an appropriate submergence ratio of the asymmetric lips and lip thickness, the total hydrodynamic efficiency is considerably increased for a broad bandwidth frequency. Experimental tests were carried out by [28] on a fixed offshore OWC model to study the effect of the wave spectrum shape on its efficiency. An analytical model to analyze the effects of a skirt on an OWC device integrated into an offshore wind turbine monopile was carried out by [29]. Through both experimental tests and numerical simulations, Deng et al. [30] studied the hydrodynamic performance of an offshore-stationary OWC device with an immersed horizontal bottom plate. By means of the EEM, Deng et al. [31] made a theoretical study of an asymmetric, offshore-stationary OWC device with an outstretched bottom plate. They found that the bottom plate provides an additional resonance mechanism that can be adjusted to offer a wider bandwidth with high performance.
There is still a lot to be investigated regarding the improvement of the hydrodynamic efficiency of fixed-detached OWC devices by modifying the geometric parameters of the chamber. To the best of the authors' knowledge, a theoretical study into the interaction of water waves with a fixed-detached OWC device with asymmetric thick rear and front walls, with a step and a reflecting wall, has not been reported. Apart from absorption of wave energy, fixed-detached OWC devices can be used to protect maritime structures and human activities from wave attack, as well as to enable water exchange and sediment transport [16]. They may, however, be affected by severe storms, which can cause serious structural damage, as was the case with the Mutriku Wave Energy Plant [32]. Therefore, the thickness of the back and front barriers that comprise the OWC chamber must be designed to withstand wave loads and save the structure, while the effect that this would have on the device performance must be also considered.
Compared to experimental and CFD approaches, this theoretical investigation has the advantage of quickly evaluating different parameters of a fixed-detached OWC device. The main objectives of this analysis were to examine modifications in bandwidth on the efficiency curves and the shifting of the peak resonant frequency due to different rear and front wall configurations, chamber size, and the length of the gap between the OWC chamber and the reflecting wall.
The present study is organized as follows: Section 2 defines the physical problems corresponding to the water wave radiation and scattering by a fixed-detached OWC device, based on the potential flow theory. In Section 3 the methodologies for solving the BVPs are presented. First, the matched EEM is described, where the eigenfunction expansions along with their orthogonal properties, the continuity of velocity potential, and flux are used to reduce the BVP into a linear system of algebraic equations. Then the BEM, which is a numerical technique based on an integral equation formulation and boundary discretization, is presented. In this work, the BEM with three-noded elements is used to account for the variation of the potential and flux on the boundaries. Section 4 describes the procedure for obtaining the energy relations. In Sections 5 and 6, convergence and truncation analyses are carried out, respectively, and a comparison is made with the experimental results of [16,17]. Section 7 shows the results for various modelled cases. Finally, in Section 8, the main conclusions of this study are drawn. Figure 1 shows the fixed-detached OWC device, the fluid and the Cartesian coordinate system, with the vertical coordinate, z, measured positively upwards and the undisturbed free surface located at z = 0. The waves approach the device from the positive x direction. The device is composed of two partially immersed front and rear walls, separated by distance b. The rear wall is located at x = −w 1 , with a draft a 1 and a thickness w 1 , while the front wall is located at x = b, with draft a 2 and thickness w 2 . A step under the OWC device is considered with a height h − h e . There is assumed to be a turbine connecting the trapped air inside the chamber and the atmosphere.

Boundary Value Problem
The BVP is divided into five regions, as shown in Figure 1, and the boundaries are denoted by

•
The rear wall by The internal free surface inside the chamber by The fluid is inviscid and incompressible and the flow is irrotational; thus, the velocity potential theory is applied. The wave motion is represented by linearized wave theory, disregarding surface tension effects. A simple harmonic flow with angular frequency ω is assumed, where the velocity potential is defined by Φ(x, z, t) with Φ(x, z, t) = Re{φ(x, z)e −iωt }, with Re{ } denoting the real part of a complex expression. Under these assumptions, the spatial velocity potential φ satisfies the Laplace equation The no-flow boundary conditions are applied at the bottom and the two barriers as The continuity of pressure and horizontal velocity on the common interfaces of the five regions are given by The linearized internal and external free surface boundary conditions are respectively, where p is the harmonic pressure distribution in the internal free surface, K = ω 2 /g, where g is the gravitational constant and ρ the seawater density. The velocity potential is separated into a scattered potential φ S and a radiated potential φ R , as proposed by [33], where φ S satisfies Equations (1)-(4) with p = 0 inside the chamber; while φ R satisfies Equations (1)-(4) with Equation (4) replaced by The Sommerfeld radiation condition is imposed at the left-and right-hand ends (x → ±∞) of the domain where φ D represents the diffracted potential, which together with the incident potential φ I , composed the scattered potential φ S , while the wavenumber k = 2π/λ is the positive real root of the wave dispersion relation given by and λ is the wavelength. The time harmonic induced volume flux across F i is also decomposed into scattering and radiation volume fluxes, q S and q R , respectively, as follows

Definition of Velocity Potentials
The solution technique for the radiation and scattering problems based on the matched EEM is explained in this subsection. For the radiation problem, the method of separation of variables is used to determine the velocity potentials in the five regions represented by the eigenfunctions expansions as for regions 1 to 5, respectively, whereγ n = nπ/(h e − a 1 ) andγ n = nπ/(h e − a 2 ) for n = 0, 1,... . Additionally, the term −1/K on φ R 1 was introduced to satisfy Equation (6) and the coefficients A R and B R , A R n and B R n for n = 1, 2,... and C R n , D R n , F R n , G R n , H R n and I R n for n = 0, 1,... are the unknown constants to be determined. The eigenfunctions Ψ,χ,Ψ andχ are defined by Ψ n (z) = 1 √ N n cos k n (z + h), for n = 0, 1, 2, ..., Ψ n (z) = 1 N n cosk n (z + h e ), for n = 0, 1, 2, ..., (11b) χ n (z) = 1 Ŵ n cos nπ z + h e h e − a 1 , for n = 0, 1, 2, ..., (11c) where N n = 1 2 1 + sin 2k n h 2k n h , for n = 0, 1, 2, ..., (12a) N n = 1 2 1 + sin 2k n h e 2k n h e , for n = 0, 1, 2, ..., (12b) with k 0 = −ik andk 0 = −ik, where ik n and ik n for n ≥ 1 indicate the imaginary roots of the dispersion relation given by Equation (8), while fork this is given by The eigenfunctions in Equation (10) are orthonormal with respect to the inner product where δ nm is the Kronecker delta. As described for the radiation problem earlier, the velocity potentials for the scattering problem are also expressed by the appropriate eigenfunctions. The expansion of velocity potential φ R 2 and φ R 4 remain the same as in the case of the radiation problem, while for velocity potentials in regions 1, 3, and 5, these are expressed as follows where the associated eigenvalues and eigenfunctions are the same as defined above. The far-field behavior of the scattered and radiated potentials is given by where A R and B R represent the amplitude of the radiated wave to ±∞, respectively, while R and T are the amplitude of the reflected and transmitted waves in the scattering problem, respectively.

Matching of Regions
Applying of the second matching condition Equation (3a) at x = −w 1 on the velocity potentials described by Equation (10a,b), together with the boundary condition Equation (2) at x = −w 1 and −h ≤ z ≤ −h e and −a 1 ≤ z ≤ 0 applied to the velocity potential Equation (10a), adding their results, exploiting the orthonormality of the eigenfunctions as defined in Equation (14a), and truncating the infinite series up to N terms, we have Next, by using the first matching condition of Equation (3a) at x = −w 1 with the velocity potentials, Equation (10a,b), exploiting the orthonormality of the eigenfunctions as defined in Equation (14c) and truncating the infinite series, we obtain By applying the second matching condition of Equation (3b) at x = 0 in the velocity potentials of Regions 2 and 3, Equation (10b,c), respectively, in conjunction with the application of boundary condition Equation (2) at x = 0 and −a 1 ≤ z ≤ 0 in the velocity potential Equation (10c), adding their results, making use of the orthonormality of the eigenfunctions as given by Equation (14b) and truncating the infinite series up to N terms, gives Next, by applying the first matching condition of Equation (3b) at x = 0 on Equation (10b,c), multiplying byχ m , integrating from −h e to −a 1 , and using the orthonormal properties of the eigenfunctionχ m Equation (14c), we obtain Then, applying the second matching condition of Equation (3c) at x = b in the velocity potentials Equation (10c,d), together with the application of boundary condition Equation (2) at x = b and −a 2 ≤ z ≤ 0 in the velocity potential Equation (10c), adding their results, making use of the orthonormality of the eigenfunctions as given by Equation (14b) and truncating the infinite series up to N terms, we obtain By using the first matching condition of Equation (3c) at x = b with the velocity potentials, Equation (10c,d), exploiting the orthonormality of the eigenfunctions, as defined in Equation (14d), and truncating the infinite series, we have Then, by applying the second matching condition of Equation (3d) at x = f w on the velocity potentials of Regions 4 and 5, Equation (10d,e), respectively, in conjunction with the application of boundary condition Equation (2) at x = f w and −h ≤ z ≤ −h e and −a 2 ≤ z ≤ 0 on the velocity potential Equation (10e), adding their results, making use of the orthonormality of the eigenfunctions as given by Equation (14a) and truncating the infinite series, we obtain Finally, by using the first matching condition of Equation (3d) at x = f w with the velocity potentials, Equation (10d,e), exploiting the orthonormality of the eigenfunctions as defined in Equation (14d), and truncating the infinite series yields Thus, Equations (18)-(25b) provide a linear system of algebraic equations for solving the unknowns. The procedure for solving the scattering problem is similar to that of the radiation problem described above.

Fixed-Detached OWC Device near a Reflecting Wall
In this subsection, the interaction of water waves with a detached OWC device near a rigid, vertical wall is analyzed. For this case, the OWC device is at a finite distance, L, from the reflecting wall, as shown in Figure 2. The velocity potentials in 2, 3, 4, and 5, Equation (10b-e), and the boundary conditions, remain the same as in the previous subsection, except for the radiation condition at the left-hand end, which is now given by and the velocity potential φ 1 now has the following form Through the application of the matching conditions, as in Equation (3a-d), along with the orthogonality conditions defined in Equation (14a-d), a system of equations to determine of the unknowns is obtained. This is similar to that shown in Equations (18)-(25a), with only a difference in Equation (18), which now takes the following form:

Boundary Integral Equation
The Laplace Equation (1), in its boundary integral representation is given by where φ and ∂φ/∂n represent, respectively, the unknown velocity potential and its normal derivative with respect to the field point Y ξ ,η on the boundary Γ; X(x, z) is the source point inside the domain Ω; α = τ/2π, where τ is the angle in radians between points X and Y [34]; ψ represents the fundamental solution of Laplace equation given by where r = x −ξ 2 + (z −η) 2 is the distance between points X and Y; and ∂ψ/∂n is its The boundary integral Equation (29) is then discretized by subdividing the boundary Γ into a number of boundary elements NE as follows with i indicating the boundary node and j the boundary element. The boundary is then divided into curved quadratic elements, where φ and ∂φ/∂n are assumed to have three different values in each element, see [35]. The fluxes will then have a single nodal value for each element, as explained by [36]. Thus, the possibility of different values between neighboring elements is considered, but in the connection between two adjacent elements, the velocity potentials only have a single value. Therefore, fluxes are organized in a 3 ×NE array, whereas the velocity potentials are in an N bem array, where N bem is equivalent to the number of nodes equal to 2NE for closed boundaries [37].

Matching of subdomains
The radiation and scattering BVPs are divided into three regions, as shown in Figure 3, with common interfaces on both sides of Region 2. Since the far-field boundaries at ±∞ must be placed far enough from the OWC chamber to avoid the effect of local disturbances, the domain as a whole becomes very long, and numerical problems associated with the integration of the fundamental solution over long distances arise [34]. The subdomain method is employed to overcome these numerical issues. This also helps to avoid numerical instabilities that appear when the wall thickness is very small [13].
The subdomain method is used to match the regions at the common interfaces. To obtain the final solution matrix, each subdomain is treated separately to form matrices that are coupled together according to the relevant boundary conditions. The continuity of the velocity potential and the flux at the interfaces should, therefore, be defined to match the regions [34,38]. For this purpose, the nodes in the left and right interfaces are assumed to be in perfect contact, as shown in Figure 3. These velocity potential and flux continuity conditions at the interfaces are given by respectively.  The previously defined boundary conditions, Equations (2), (4), (6), (7), and (26), are then used, and the unknown variables are shifted to the left-hand side, while the known variables to the right-hand side. Then, the matrix below is obtained where the square matrix [A] has dimensions of N bem × N bem ; the vector of unknown velocity potentials or fluxes {X}, with dimensions N bem × 1; and the known vector {B} has dimensions of N bem × 1 [36].

Efficiency Relations
As proposed by [33], the radiation volume flux q R is divided into real and imaginary parts as iω p where Z =B − iÂ is the complex admittance, while the real coefficientsÂ andB are given byÂ which are related to the added mass and the radiation damping in a rigid body system, respectively [9]. Applying Green's theorem to the radiation potential and its conjugate yields where U is the closed boundary composed of where F +∞ and F −∞ are the imposed far-field boundaries as x → ∞ and x → −∞, respectively, at −h < z < 0. By using the boundary conditions Equations (2) and (3), the contributions from the boundaries B d , B r , B f and F e vanish. Then, with the aid of Equations (6), (16a,b) and (34), the contributions from F i , F +∞ and F −∞ are calculated, and an expression for the radiation dampingB as a function of the complex amplitudes of the right and left radiated waves A R and B R , respectively, is obtained Now, applying the Green's theorem on the radiation and the scattered potentials as in Equation (36), and using the boundary conditions as in the previous case, a relation between the right radiated amplitude A R and the induced volume flux in the scattering problem q S is obtained as follows and is non-dimensionalized by the incident volume flux as in [9], giving rise to Equations (37) and (38) can be computed by solving the radiation problem alone and are comparable to those derived by Evans [33] and used by [9,13]. Equation (38) is linked to the volume flux across F i due to the incident plus diffracted wave fields, and it is proportional to the amplitude of the radiated waves in the direction from which the incident wave comes [33]. Now, by assuming that the pressure inside the chamber is uniform and the air exits to the atmosphere through the turbine without a phase lag, we have where the real control parameter Λ is related to the linear turbine damping induced by the airflow. This is a general characteristic of Wells turbines, which, with constant rotational speed, show a linear relationship between pressure and volume flow rate [39,40]. After combining Equations (9), (34) and (40), the imposed internal pressure gives Now, the total rate of work performed by the pressure forces inside the OWC, Q(t) × P(t), is averaged over one wave period to obtain the total rate of power absorbed per unit width of pressure distribution as where the horizontal bar (−) denotes the complex conjugate. Now, by inserting Equations (9) and (34) into Equation (42) and simplifying, we have where ifB −1 exists, the maximum work gives where Λ = Z for maximum power. Thus, after combining Equations (41) and (43), we obtain Now, the optimal power conversion efficiency must be obtained. This is performed by finding the optimum value of Λ, by applying zero value to the derivative with respect to Λ for the squared-right term inside the brackets of Equation (45) [11], thus yielding and after substituting the above expression into Equation (45), we have whereÂ,B and Λ are a function of the angular frequency ω.
Hence, an expression for the maximum hydrodynamic efficiency is given as with η max in the range 0 to 1. As in [9], the following non-dimensionalized quantities are defined which represent the coefficients of radiation susceptance and radiation conductance, respectively. Finally, after inserting the above coefficients into Equation (48), the hydrodynamic efficiency η max gives

Convergence and Truncation Analyses
A convergence study for both EEM and BEM was done prior to performing the computations of the results. For the calculations, the water depth is assumed to be 7.90 m, and a wave period T bounded by 2.50 ≤ T ≤ 30 s, as in [35]. In the case of the EEM, see Table 1, a convergence analysis for the maximum hydrodynamic efficiency η max is given for different values of the non-dimensional frequency Kh(= 0.5, 1.0, 1.5, 2.0, 2.5, 3.0 and 3.5). From Table 1, it can be seen that roughly 30 evanescent modes are adequate to ensure that the analytical results converge to within three decimal places. However, it should be noted that the OWC geometry is relevant for the convergence, and the number of equations needed in the matched EEM, especially for the lower corners of the walls, since as "w 1 " and "w 2 " decrease, more terms are needed. For the case without the step, the results converge faster, to the desired accuracy even with 20 terms in the infinite series sums. However, on average, 30 terms are sufficient for the desired accuracy; therefore, in the present calculations, N = 30 is used.
In Table 2, a truncation analysis was first carried out to avoid the effect of local disturbances at the far-field boundaries. From Table 2, it is seen that at a distance of 4 times the water depth between the left and right faces of the OWC device and the left and right far-field boundaries, respectively, the results converge to five digits. In Table 3, using the BEM, the results of η max for the same Kh values are given for different numbers of nodes (N). In this case, around 800 nodes (400 quadratic elements) were determined to be sufficient to achieve numerical results convergence to three decimal places. As a result, the three-region BVP is discretized with around 800 nodes.

Comparison with Experimental Results
To validate the results obtained by the above-mentioned methods, the experimental results of [16,17] were employed. In Figure 4a, comparisons are shown of the experimental data of [16] and the present BEM results for the hydrodynamic efficiency versus b/λ with b/h = 1.0 and w 1 /h = w 2 /h = 0.025 for two different values of the walls draft to water depth ratio a 1 /h = a 2 /h. The damping coefficients for these results were Λ = 0.008 and Λ = 0.00005 m 4 ·s/kg for a 1 /h = a 2 /h = 0.375 and a 1 /h = a 2 /h = 0.50, respectively. To determine the damping coefficient Λ, the method described by [41] was employed. The ratios of slot opening to cross-sectional area of the OWC chamber considered by [16] were 0.625% and 1.875% for a 1 /h = a 2 /h equal to 0.5 and 0.375, respectively. In Figure 4b, the hydrodynamic efficiency versus b/λ with b/h = 1.0, a 1 /h = a 2 /h = 0.25 and w 1 /h = w 2 /h = 0.025 for two different values of the rear face to reflecting wall distance to chamber length ratio L/b is compared against the results obtained by [17]. In this case, the damping coefficients were Λ = 0.0018 and Λ = 0.0030 m 4 ·s/kg for L/b = 0.24 and L/b = 0.97, respectively. Here, the opening ratio considered by [17] was 1.25% in both cases.
From these figures, it can be seen that the BEM results follow a similar trend to that observed in the experimental data. However, for the case of the highest draft a 1 /h = a 2 /h = 0.5 with the smallest opening ratio, 0.625%, as is the case of the PTO mechanism in [16], the difference is more significant. This discrepancy obtained with the present results arises from the assumption of an ideal fluid, where viscous effects and flow separation due to the OWC structure are ignored. This also leads to an overestimation of the hydrodynamic performance. Other factors that also contribute to the discrepancy are the modelled PTO system, as well as the energy loss through it by viscous dissipation during the tests.   [17] for an OWC device near a reflecting vertical wall.

Results and Discussion
Analytical and numerical results based on the previously described methodologies are provided in this section. The influence of the rear and front wall thickness to water depth ratio (w 1 /h, w 2 /h), the chamber length to water depth ratio (b/h), the rear and front walls draft to water depth ratio (a 1 /h, a 2 /h), the step height to water depth ratio (h e /h) and the chamber to wall distance to water depth ratio (L/h) on the hydrodynamic performance are analyzed.

Asymmetric Fixed-Detached OWC Device over a Flat Bottom
The theoretical results of the hydrodynamic efficiency η max versus Kh for different chamber lengths b/h(= 1/4, 1/2, 3/4, and 1) are shown in Figure 5a. It can be seen in this figure that increasing the length of the chamber increases the bandwidth of the efficiency curves. In this regard, large motions occur in the internal water column when the fluid is excited into a resonant piston-like motion by the incident wave [9]. The natural frequency of oscillation can be calculated for small values of the b/h ratio, allowing the water enclosed within the chamber to be treated as a solid body. By simple hydrostatic modelling, it results that the expected resonance occurs at Kh ≈ h/a 1 . For the cases shown in Figure 5a, this resonance would occur at Kh = 2, which appears to be approached for small b/h. Furthermore, it is observed that for large b/h ratios, the peak frequency value is shifted to lower frequencies. This is because increasing the chamber length increases the horizontal distance a fluid particle must travel during a period of motion, resulting in a reduction in the value at which resonance occurs. Furthermore, since a longer chamber allows for a greater local fluid motion, the solid-body concept of resonance breaks down, and the oscillation amplitude decreases [9].
The results of the hydrodynamic efficiency versus the non-dimensional frequency Kh for different thickness ratios of the front and rear walls w 1 /h = w 2 /h(= 1/8, 1/4, and 1/2) are shown in Figure 5b. This figure shows that the bandwidth of the efficiency curves is reduced as the thickness of the walls is increased, and their peak frequency value is moved to lower values of the non-dimensional frequency Kh. This drop in efficiency at medium and high frequencies can be explained by the fact that the energy transfer due to wave motion for relatively small wavelengths decreases as wall thickness increases. On the other hand, Figure 5c,d shows the variation of η max versus Kh for different thicknesses of the rear and front walls, respectively. Similar to Figure 5b, in these figures it is observed that the bandwidth of the efficiency and the peak frequency magnitude decrease when the thickness of the front or rear walls increases. Furthermore, a second resonance peak that occurs when conditions are similar to those in a closed tank with parallel walls is observed in Figure 5d. In such conditions, the incidence wave frequency is such that the fluid in the chamber is excited in an anti-symmetric sloshing mode [9]. In this case, the sloshing frequencies are seen to occur at values of the dimensionless wavenumber kb = nπ, with n being the sloshing mode. Thus, the second peaks in η max caused by the first sloshing frequency occur at Kh ≈ π for b/h = 1.    Figure 6a shows the effect on the efficiency of different submergence ratios a 1 /h = a 2 /h(= 1/8, 1/4, 1/2, and 3/4) versus Kh. This figure shows that the effective area under the efficiency curve and the magnitude of the natural frequency both increase when the walls submergence ratio a 1 /h and a 2 /h decreases. However, while a small draft results in a better performance for a wider frequency range, in practice, tidal fluctuations and extreme waves may compromise its effectiveness since a small draft may mean that the trough of a wave propagates below the front wall. In this instance, the pressure inside the chamber would be equivalent to the ambient pressure, resulting in no power available within the OWC device. Furthermore, Figure 6b,c shows the variation of η max versus Kh for an asymmetric OWC device with different drafts of the front and rear walls, respectively. As in Figure 6a, these figures show that the bandwidth of the efficiency and the peak frequency magnitude decrease when the submergence of the front or rear walls increase. Additionally, in Figure 6b a second resonance mechanism due to the first sloshing frequency at Kh ≈ π for b/h = 1 is observed. Figure 7a-c shows the analytical results for the radiation susceptance and radiation conductance coefficients, µ and ν, respectively, and the non-dimensional induced volume flux due to the scattering potential |q S |/|q I | versus Kh for different values of b/h(= 1/4, 1/2, 3/4, and 1). Figure 7a shows that the peak and the range of positive values in the radiation susceptance coefficient increase when chamber length increases. From Figure 7b it is observed that the bandwidth in the radiation conductance curves decreases as b/h increases, similar to the trend observed in η max , see Figure 5a. Moreover, the zero values in µ are associated with the peak values in ν, |q S |/|q I |, Figure 7c, and η max , due to the fundamental resonance inside the chamber.

Asymmetric Fixed-Detached OWC Device over a
Step Figure 8a shows the analytical and numerical results of the hydrodynamic efficiency η max versus Kh for different chamber lengths to water depth ratios b/h(= 1/4, 1/2, 3/4, and 1). This figure shows that the value of the resonance peak frequency decreases when the length of the chamber increases and, compared to the case when h e /h = 1 (Figure 5a), it is observed that the step contributes to a reduction in the value of Kh at which resonance occurs.   Figure 8b shows the theoretical results of the hydrodynamic efficiency versus the nondimensional frequency Kh for different thickness ratios of the front and rear walls w 1 /h = w 2 /h(= 1/8, 1/4, and 1/2). In this case, it is seen that a reduction in the bandwidth of the efficiency curves is obtained when the thickness of the walls increases. Moreover, it is also observed that the magnitude of the resonance frequency Kh is reduced since, for a large thickness, longer waves are needed to excite the internal water column. On the other hand, the variation of η max versus Kh for different thicknesses of the rear and front walls are shown in Figure 8c,d, respectively. Again, it is observed that whether the thickness of the rear or front wall increases, the resonance frequency Kh is reduced.
The effects of different submergence ratios a 1 /h = a 2 /h(= 1/8, 1/4, and 1/2) on the efficiency versus Kh are shown in Figure 9a. It is observed that performance for higher frequencies increases when the submergence ratio a 1 /h and a 2 /h decreases in the rear and front walls, respectively. Compared to Figure 6a where a flat bottom is considered, in this case, the value of the peak frequency is slightly decreased since the gap between the lip of the walls and the bottom is reduced. On the other hand, the variations of η max versus Kh for different drafts of the front and rear walls are shown in Figure 9b,c, respectively. These figures show that the hydrodynamic performance is sensitive to variations in the draft of the walls, reducing its magnitude for high frequencies as the submergence increases. In Figure 9b it is also observed a second resonance mechanism at Kh ≈ π as in Figure 6b, while in Figure 9c this will take place at Kh ≈ 2π since b/h = 1/2. Figure 10a-c shows the analytical results for the radiation susceptance and radiation conductance coefficients, µ and ν, respectively, and the non-dimensional induced volume flux due to the scattering potential |q S |/|q I | versus Kh for different values of w 1 /h = w 2 /h(= 1/8, 1/4, and 1/2). In Figure 10a it is observed the peak in µ increases when the thickness in the walls also increases, while the peak shifts to the lower-frequency side. However, it is observed that the range of negative values increases as the thickness increases. Figure 10b shows that the curve in the radiation conductance coefficient becomes sharper as the thickness increases, which indicates a reduction in the frequency range for the energy transfer into the system. A decreasing trend in the peak of the non-dimensional induced volume flux due to the scattering potential when the thickness decreases is observed in Figure 10c.

Asymmetric Fixed-Detached OWC Device near a Reflecting Wall
The analytical and numerical results of the hydrodynamic efficiency η max versus Kh for an OWC device in the presence of a reflection wall with different chamber length to water depth ratios b/h(= 1/4, 1/2, 3/4, and 1) are shown in Figure 11a. Compared to Figure 5a, it is observed that the efficiency bandwidth is reduced when a reflecting wall is considered. In this case, it is observed that due to the near trapped waves in the gap between the OWC chamber and the rigid vertical wall, a medium-height peak at low frequencies appears in the efficiency curves. This peak increases as the chamber length to water depth ratio increases. Furthermore, a zero efficiency value appears at 0.45< Kh <0.61 that reduces the possibility of wave energy extraction at large periods.
The results of the hydrodynamic efficiency versus Kh for various thickness ratios w 1 /h = w 2 /h(= 1/8, 1/4, and 1/2) are shown in Figure 11b. It is observed that an increase in the thickness of the walls leads to both reduction of the bandwidth of the efficiency curves and the peak frequency value. Furthermore, it is seen that when w 1 /h and w 2 /h increase, the medium-height peak slightly increases, while the zero efficiency value is shifted to lower frequencies.    Figure 11c shows that with an increase in the rear wall thickness, the medium-height peak at low frequencies and the resonant frequency value of maximum efficiency both decrease. On the other hand, in Figure 11d for b/h = 1/2, it is observed that the bandwidth of the efficiency curves is highly reduced, and the medium-height peak is increased when the front wall thickness increases. Furthermore, Figure 6a-d shows a second resonance mechanism due to the first sloshing frequency at Kh ≈ π.
The effect of different draft ratios a 1 /h = a 2 /h(= 1/8, 1/4, and 1/2) on the efficiency versus Kh for a symmetric OWC device is shown in Figure 12a. It is observed that the peak resonant frequency and the efficiency bandwidth are significantly reduced when the ratios a 1 /h and a 2 /h increase. Furthermore, the frequency value of zero efficiency increases and the peak at low frequencies decreases when the submergence in the front and rear walls decreases.
The variations of maximum efficiency η max versus the non-dimensional frequency Kh for different drafts of the rear and front walls are shown in Figure 12b,c, respectively. In both of them, it is observed that by increasing the submergence of one of the walls, the values of the frequency at which resonance occurs are decreased. However, in Figure 12b, the efficiency bandwidth is increased when the submergence of the front wall decreases, while the peak at low frequencies decreases. On the other hand, in terms of the medium-height peak, Figure 12c shows an opposite trend to that observed in Figure 12b.   The results of the hydrodynamic efficiency versus the non-dimensional frequency Kh for different lengths of the gap between the rear and reflecting walls L/h(= 1, 2, 3, and 4) are shown in Figure 13a,b. In these two figures, it is observed that the bandwidth of the efficiency increases, while the medium-height peak at low frequencies reduces when the distance of separation between the fixed OWC and the reflecting wall increases. It is also observed that at L/h = 3 and 4, the shape of the curve tend to that observed in Figure 5a for b/h = 1, although a peaky trend is observed for L/h = 4. All these aspects should be considered at the design stage since the near trapped waves generated between the reflecting wall and the OWC device can reduce the efficiency bandwidth of a fixed-detached OWC device at high periods.
The analytical results for the radiation susceptance and radiation conductance coefficients, µ and ν, respectively, and the non-dimensional induced volume flux due to the scattering potential |q S |/|q I | versus the non-dimensional frequency Kh for different values of L/h(= 1, 2, 3, and 4) are shown in Figure 14a-c. Figure 14a shows a decrease in the negative values of µ when the gap between the rear and reflecting walls also decreases. As in [10], in Figure 14a,b it is observed that at resonance, the total variation from positive to negative in the radiation susceptance coefficient is related to the peak values of the radiation conductance. Finally, from Figure 14b,c it is observed that the peak in ν and |q S |/|q I |, respectively, at the fundamental resonance increase when L/h decrease.

Conclusions
The influence of the geometric parameters of an asymmetric fixed-detached OWC device on the hydrodynamic efficiency was theoretically analyzed. The BEM with quadratic elements and the matched EEM were employed to solve the radiation and scattering problems. Comparisons were made between these two methods and a good agreement was obtained. A comparison was made between the experimental data of [16,17] giving a satisfactory agreement. Results for different hydrodynamic quantities of interest were then obtained for different physical configurations. The main conclusions of this study are as follows: • An increase in the chamber length leads to an increase in the bandwidth of the efficiency curves. This is similar to the findings for a land-fixed OWC device reported by Evans and Porter [9]. • By increasing the thickness of the rear and front walls, the bandwidth of the efficiency curves and their peak frequency value are reduced. This reduction in efficiency, especially for short wavelengths, is due to the decrease in energy transfer owing to the wave motion when a large thickness is considered. • Both the effective area under the efficiency curves and the magnitude of the natural frequency increase when the submergence of the walls decreases. This is because a larger gap between the front wall lip and the bottom allows more energy to be transferred to the water column. • Regarding the presence of a step, it was observed that this reduces the frequency at which resonance occurs since the effect is similar to reducing the gap between the lip of the walls and the bottom. • When the separation distance is shorter, the near trapped and standing waves in the gap of the OWC chamber and the reflecting wall play a significant role in the hydrodynamic performance. It was observed that the efficiency bandwidth is reduced, and a medium-height peak appears at low frequencies. • At low frequencies, the presence of constructive and destructive wave interference from the OWC device and the waves reflected by the wall results in zero efficiency. This phenomenon is extremely sensitive to changes in the OWC parameters, and engineers should take it into account when constructing a fixed-detached OWC device near a vertical wall.
Finally, this work is only a theoretical investigation of the hydrodynamic performance of a fixed-detached OWC device. The construction of a fixed-detached OWC device near a rigid vertical wall may constitute an alternative for reducing the wave load on seawalls by absorbing the incident wave energy and can be extended to include an array of fixed-detached OWC devices. Numerical simulations and experimental tests of other geometrical conditions and non-linearities on the air compressibility and turbine damping due to viscous effects should be carried out in the future to improve the present results.