Next Article in Journal
Characterization of Overtopping Waves on Sea Dikes with Gentle and Shallow Foreshores
Next Article in Special Issue
The Potential for Tidal Range Energy Systems to Provide Continuous Power: A UK Case Study
Previous Article in Journal
Comparison of Wave Energy Park Layouts by Experimental and Numerical Methods
Previous Article in Special Issue
Wave Energy Resource Assessment for Exploitation—A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Influence of the Chamber Configuration on the Hydrodynamic Efficiency of Oscillating Water Column Devices

by
Ayrton Alfonso Medina Rodríguez
1,
Jesús María Blanco Ilzarbe
2,*,
Rodolfo Silva Casarín
1 and
Urko Izquierdo Ereño
2
1
Institute of Engineering, National Autonomous University of Mexico, Circuito Escolar, Mexico City CP 04510, Mexico
2
Department of Nuclear Engineering and Fluid Mechanics, Faculty of Engineering in Bilbao-UPV/EHU, Plaza Ingeniero Torres Quevedo, 1, CP 48013 Bilbao, Spain
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2020, 8(10), 751; https://doi.org/10.3390/jmse8100751
Submission received: 31 August 2020 / Revised: 21 September 2020 / Accepted: 22 September 2020 / Published: 26 September 2020

Abstract

:
Based on the two-dimensional linear wave theory, the effects of the front wall thickness and the bottom profile of an Oscillating Water Column (OWC) device on its efficiency were analyzed. Using the potential flow approach, the solution of the associated boundary value problem was obtained via the boundary element method (BEM). Numerical results for several physical parameters and configurations were obtained. The effects of the front wall thickness on the efficiency are discussed in detail, then, various configurations of the chamber bottom are presented. A wider efficiency band was obtained with a thinner front wall. In a real scenario having a thinner front wall means that such a structure could have less capacity to withstand the impact of storm waves. Applying the model for the case of the Mutriku Wave Energy Plant (MWEP), findings showed that the proposed bottom profiles alter the efficiency curve slightly; higher periods of the incoming water waves were found. This could increase the efficiency of the device in the long-wave regime. Finally, the numerical results were compared with those available in the literature, and were found to be in good agreement.

1. Introduction

It has been suggested that wave power has the potential to provide most of the world’s electricity needs in the short term [1]. A wide variety of systems have been proposed, of which only a few have reached full-scale prototype deployments [2]. Among the deployed systems, the OWC system has been shown to be one of the most promising devices. It is probably the system that has been most studied and is one of the few to have been tested at full-scale. OWCs can be located offshore, near-shore or on the shoreline, and placed on the seabed or fixed to a rocky cliff [3]. Since the design and construction of OWCs are strongly site-dependent, their location and anchorage points are of the most critical aspects, as well as the most influential in economic terms.
In light of this, installing an OWC device into a breakwater was seen as a means to provide many benefits and thus encourage further development of OWC technology [4]. The breakwater provides shelter and contributes to coastal protection by reducing wave reflection. An OWC power plant within a breakwater has the advantage of being relatively easy to install and maintain, having no mooring systems and underwater electric cables. With construction and maintenance costs shared, and the operation of the power plant being easier, energy extraction is more cost-effective. Although the waves found near the coast are less energetic, this can, in part, be compensated by the natural concentration of wave energy due to shoaling, refraction and diffraction [5].
The first integrated breakwater-OWC system was built in the port of Sakata, Japan, in 1990 [6]. Subsequently, the Basque Energy Agency (Ente Vasco de la Energía or EVE for its acronym in Spanish) employed this concept in Mutriku, The Basque Country, Spain, opening the Mutriku Wave Energy Plant in July 2011, Figure 1. This plant consists of 16 units built onsite that are 4.5 m wide, 3.1 m depth, and 10 m high (above Maximum Equinoctial Spring Tide Low Water). For each unit, a hole of 0.75 m diameter leads to a Wells turbine and electrical generator of 18.5 kW [7], yielding the total 296 kW with a 100 m breakwater. The MWEP section of the breakwater is the first multiple OWC plant in the world and is currently the only OWC device in operation that regularly supplies power to the grid. However, regarding its performance, the initial expectations have not been met because of the poor design in some of the chambers that provide moderately different pressure at the inlet of their turbines [8]. This is because the breakwater that houses the Wells turbines was manly designed to maximizing the protection of Mutriku harbour and not for wave energy harnessing.
In this context, another factor contributing to the plant’s reduced electricity output could be the changes made to the front face of the original design. The area is regularly affected by severe storms and during the construction of the power plant, three storms hit the MWEP producing severe structural damage to a number of the OWC chambers [12]. As a consequence, the front face of the chambers was reinforced to withstand the wave loads, using prefabricated concrete slabs, so that now the thickness of the front wall has doubled the length of the chambers, Figure 2. This alteration was made to save the structure of the plant, but without considering the effect that this would have on the device performance. The main focus of this work is, therefore, the evaluation of the influence that the front wall thickness of an OWC-breakwater system has on its hydrodynamic performance as an OWC.
It is important to note that the success of the OWC system will depend on the coupling between the chamber and the power take-off (PTO) system. In this sense, a good turbine design, an effective control strategy and the matching of the turbine to the OWC collector to ensure efficient collector operation are essential [13]. Furthermore, the peak performance of most OWC systems occurs at resonance, which takes place when the incident wave frequency coincides with the natural frequency of the converter. Therefore, to operate optimally at resonance, the OWC chamber design plays a significant role to obtain higher efficiencies. Typically the chamber geometrical configuration is chosen to produce a column whose natural frequency of oscillation coincides with that of the most occurring wave at the location where the OWC will be installed [14]. In this sense, the variability of sea state conditions can influence the OWC feasibility, because, once installed, the size and shape of the structure can be hardly modified.
Over the last years, a variety of analytical, numerical and experimental techniques have been employed to study the effects of the geometrical configuration of an OWC on its hydrodynamic efficiency. Wang et al. [15] studied the hydrodynamic performance, both numerically and theoretically, of an OWC device with arbitrary topography near the shoreline. They reported that as the bottom slope increases, the peaks in capture-width ratios become lower frequency values, concluding that a change in water depth at the shoreline has a significant effect on the hydrodynamic performance of an OWC. The effect of front wall geometry on OWC hydrodynamic efficiency was analyzed by Thomas et al. [16]. Their experimental study concluded that the overall peak in hydrodynamic efficiency is not influenced greatly by the front wall geometry. Martins-rivas and Mei [17] presented a theoretical model for a cylindrical OWC installed on a cliff coast. It was found that air compressibility helps optimize the power absorption efficiency while the angle of incidence significantly affects the waves outside the chamber but not the averaged response inside or the capture length of energy absorption. Şentürk and Özdamar [18] carried out a theoretical analysis of an OWC which had a gap in its fully-submerged front wall. They showed that it is possible to increase the efficiency of an OWC with a surface piercing, barrier-type front wall when appropriate geometrical parameters are taken into consideration.
Rezanejad et al. [19] analyzed the impact of stepped bottom topography in the efficiency of a nearshore OWC device. They reported that there are significant effects when there is a stepped bottom profile outside of the chamber. Ning et al. [20] studied the performance of a fixed OWC device based on a time-domain higher-order BEM in a 2D fully nonlinear numerical wave flume. They investigated the hydrodynamic performance, with, and without, a bottom slope in the OWC chamber, and reported that the geometric parameters of the air chamber have a significant influence on hydrodynamic efficiency. The configuration of the bottom profile on the hydrodynamic performance of the OWC was investigated experimentally by Ashlin et al. [21]. Flat, circular, curved and sloped bottom profiles were tested in a wave flume. It was found that the OWC with a circular curved bottom profile was more effective in wave energy conversion, as was the wave amplification factor inside the chamber. The effects of the incident wave amplitude and geometric parameters on the hydrodynamic efficiency of a fixed OWC were investigated by Ning et al. [22]. They concluded that the incident wave amplitude and the bottom slope have a small influence on the resonant frequency, while the optimal hydrodynamic efficiency increases with an increase of bottom slope. A theoretical model based on linear potential flow theory to study the performance of a circular cylindrical OWC along a vertical coast/breakwater without the thin-wall restriction was proposed by [23]. The authors concluded that the incident wave direction and the thickness of the circular chamber wall both play an important role in the wave power captured by the OWC. By employing the eigenfunction matching method, Zheng et al. [24] developed a theoretical model to evaluate the hydrodynamic performance of multiple circular cylinder OWCs installed along a vertical straight coast. It was found that due to the effects of constructive wave interference from the OWCs array and the coast, the hydrodynamic performance of the OWC devices was enhanced significantly for a certain range of wave conditions. Zheng et al. [25] studied the effect of the radius of the entrance to the chamber and the finite wall thickness of the tubular-structure. They demonstrated that wave power extraction is greater with a thinner chamber wall thickness, mainly in terms of a broader primary band of efficiency curves. Using a coupled eigenfunction expansion—BEM—Koley and Trivedi [26] analyzed the hydrodynamic performance and efficiency of an OWC device placed on an undulated seabed. They concluded that the OWC structural design and bottom profile can significantly increase the hydrodynamic efficiency. A 2D BEM model for analyzing the OWC’s response in general bathymetry regions was carried out by [27]. They showed that the effects of the bottom slope and curvature on the OWC performance could be important, especially when the wave climate leads the site-specific optimal design to low resonance frequencies.

2. Aims and Methodology

In the specialized literature, there still remains a lot to be investigated regarding the improvement of the OWC efficiency by modifying its structural configuration. To the authors’ knowledge, a numerical study for analyzing the interaction of water waves with an OWC-breakwater system considering a wide front barrier has not been examined in the past. The fundamental hypothesis of the present work is that the hydrodynamic efficiency can be highly affected when a thick front barrier is employed. This reduction to the efficiency could be explained by the fact that the transfer of energy from the incoming wave to the internal free surface due to the orbital wave motion is reduced. A reduction in energy transmission then may lead to a decrease in the internal free surface oscillation for driving the air column, which consequently diminishes the output power.
Furthermore, by following the same idea of a wide front barrier but now considering the physical dimensions of a chamber in the MWEP, three different bottom profiles inside the chamber are then proposed to analyze their influence on the hydrodynamic efficiency. These proposed varying bottom profiles are a slope, a cycloid and an ellipse. This proposal is motivated by the fact that a curve bottom profile can exhibit better performance in terms of wave energy conversion and wave amplification inside the chamber, as it was experimentally studied by [21]. For this purpose, a numerical study for analyzing these curved profiles inside the chamber is also proposed.
Thus, this work examines the two-dimensional hydrodynamic interaction of ocean waves with an OWC device. Linear wave theory for a constant sea depth is employed and the viscous effects and the nonlinear air compressibility are neglected. The associated Boundary Value Problem (BVP) is then solved by the BEM employing three noded quadratic elements. The present formulation is novel in addressing the influence of a wide front wall and the use of BEM with a second-order discretization. The main interest of this work a) lies on the analysis of the bandwidth reduction on the efficiency curves due to an increment on the front wall thickness, and then, b) based on the geometric dimensions of the MWEP, to study alternatives for increasing the efficiency by considering a modification in the bottom profile inside the OWC chamber. Numerical estimates for the hydrodynamic efficiency and the radiation susceptance and radiation conductance coefficients are presented for a range of different parameters. Furthermore, numerical results for particular cases are validated with the previous results obtained by Evans and Porter [28] for a thin vertical surface-piercing barrier next to a vertical wall, and Şentürk and Özdamar [18] for an OWC with a gap on a fully submerged front wall.

3. The Boundary-Value Problem

For the present study, the Cartesian coordinate system was chosen, with the x axis corresponding to the opposite direction of the wave propagation and the z axis corresponding to the upward direction. The origin of the coordinate system lies on the undisturbed free surface and the left-vertical wall inside the chamber. The OWC is this rigid wall, situated at x = 0 , extending down to the sea bottom and complemented by a vertical, surface-piercing barrier, at x = b , with a thickness w and a draft h a , as shown in Figure 3. The front barrier is denoted by L b = { ( x , z ) : ( x = b , h a z 0 ) ( b < x < b + w , z = h a ) ( x = b + w , h a z 0 ) } , at the left side of the chamber entrance, the vertical length of the gap between the immersed tip of the barrier and the bottom is defined by L g = { ( x , z ) : x = b , h e z h a } , the rigid vertical wall by S w = { ( x , z ) : x = 0 , h < z < 0 } , the internal free surface inside the water column by S i = { ( x , z ) : 0 x b , z = 0 } , the external free surface by S f = { ( x , z ) : b + w x , z = 0 } and the bottom by S b = { ( x , z ) : ( 0 < x < b , z = h ) ( x = b , h < z < h e ) ( b < x < b + w , z = h e ) ( x = b + w , h < z < h e ) ( b + w < x < , z = h ) } .
The fluid is assumed to be inviscid and incompressible and linear wave theory is applied, ignoring the effect of surface tension. By assuming an irrotational flow and simple harmonic in time with angular frequency ω , there is thus a velocity potential Φ ( x , z , t ) with Φ ( x , z , t ) = Re { ϕ ( x , z ) e i ω t } , where Re { } denotes the real part of a complex expression and t is the time. The spatial velocity potential ϕ then satisfies the Laplace equation
2 ϕ x 2 + 2 ϕ z 2 = 0 ,
along with the no-flow boundary condition on the solid boundaries such as the barrier, the rigid vertical wall and the bottom described by
ϕ n = 0 for ( x , z ) S b , S w and L b ,
together with the continuity of pressure and horizontal velocity given by
ϕ = ϕ + and ϕ x = ϕ x + for ( x , z ) L g on x = b , x = b + w .
Inside the chamber, by imposing a pressure distribution over the internal free surface P ( t ) and after considering simple harmonic motions for the free surface η ¯ = Re { ζ e i ω t } and P ( t ) = Re { p e i ω t } , the dynamic and kinematic free surface boundary conditions are
ϕ + i g ω ζ = i ρ ω p on z = 0 , 0 < x < b ,
ϕ z + i ω ζ = 0 on z = 0 , 0 < x < b ,
and on the external free surface with p = 0
ϕ + i g ω ζ = 0 on z = 0 , b < x < ,
ϕ z + i ω ζ = 0 on z = 0 , b < x < .
Thus, by combining Equations (4) and (5), the internal and external linearized free surface boundary conditions are
ϕ z K ϕ = i ω p ρ g on z = 0 , 0 < x < b , 0 on z = 0 , b < x < ,
respectively, where K = ω 2 / g , with g being the gravitational constant and ρ the seawater density.
Here, as described by [28], the potential is decomposed into two parts as follows
ϕ ( x , z ) = ϕ S + i ω p ρ g ϕ R .
The scattered potential ϕ S represents the solution of the scattering of an incident wave coming from x = + in the absence of an imposed pressure on the internal free surface inside the chamber, satisfying Equations (1)–(6) with p = 0 ; while the radiated potential ϕ R represents the solution of the radiation problem due to the pressure imposed on the internal free surface and satisfies Equations (1)–(6) with Equation (6) replaced by
ϕ R z K ϕ R = 1 on z = 0 , 0 < x < b ,
which is due to an oscillating pressure distribution on the internal free surface in the absence of incoming waves.
The Sommerfeld radiation condition describes the far field boundary condition for the diffraction and radiation problems as follows:
ϕ D , R x i k ϕ D , R = 0 as x + ,
where ϕ D represents the diffracted potential that together with the incident potential ϕ I composed the scattered potential ϕ S , while k represents the wave number and is the real root of the wave dispersion relation given by
ω 2 = g k tanh k h ,
whose solution of this expression can be easily determined by a root-finding algorithm.
On the other hand, the time harmonic induced volume flux across the internal free surface, Q ( t ) = Re { q e i ω t } (see [28]) is given by
q = S i ϕ z d x = q S + i ω p ρ g q R ,
where q S and q R are the volume fluxes across S i in the scattering and radiation problems, respectively. Thus, by using the continuity of volume flux across the internal free surface and the gap between the barrier tip and the sea bottom, we obtain
q S , R = S i ϕ S , R z d x = L g ϕ S , R x d z .
Finally, the volume flux q R for the radiation problem is separated into real and imaginary parts as follows
i ω p ρ g q R = B ˜ i A ˜ p = Z p ,
where Z = B ˜ i A ˜ is a complex admittance and A ˜ and B ˜ are analogous to the added mass and the radiation damping coefficients of the forced oscillation of a rigid body system immersed in an ideal fluid and, following [28], are called the radiation susceptance and the radiation conductance parameters, respectively, described by
A ˜ = ω ρ g Re { q R } ,
B ˜ = ω ρ g Im { q R } ,
where Im { } denotes the imaginary part of a complex expression.

Efficiency Relations

Since in practice it may be easier to control the volume flux through the turbines than the pressure drop across it [29], a linear relationship between these two without a phase lag is assumed,
q = Λ i ϱ p ,
where Λ is a real control parameter, related to the damping induced to the airflow by the linear turbine and ϱ = ω V 0 / γ p a represents the effect of compressibility of air in the chamber with V 0 being the air volume inside the chamber, γ the specific heat ratio of air equal to 1.4 and p a the atmospheric air pressure [30]. The sign in Λ is taken to be positive since, in contrast to Equation (13), the pressure forces and volume fluxes are both measured vertically upwards. Equation (15) assumes that the pressure inside the chamber is uniform and the air exits to the atmosphere through the turbine, a characteristic of Wells turbines that has been widely investigated for OWC devices. Using Equations (11), (13) and (15) an expression for the pressure in the chamber can be found
p = q S Λ + Z i ϱ .
The total rate of working of the pressure forces inside the OWC is basically Q ( t ) × P ( t ) . By averaging this over one period, the mean power absorbed per unit width of pressure distribution is obtained
W = 1 2 Re { p * q } ,
where * denotes complex conjugate. Now, by using Equations (11) and (13) on this last expression, we obtain
W = 1 2 Re p * q S Z p = 1 2 Re { p * q S } B ˜ | p | 2 .
The expression (18) can be re-written in the form
W = 1 8 q S * q S B ˜ 1 1 2 B ˜ p q S 2 B ˜ * p q S 2 B ˜ = | q S | 2 8 B ˜ B ˜ 2 p q S 2 B ˜ 2 ,
where if B ˜ 1 exists, the maximum work is
W m a x = | q S | 2 8 B ˜ ,
when
p = q S 2 B ˜ ,
showing that Λ = Z i ϱ * for maximum power.
By combining previous Equations (16) and (19), it is finally obtained Equation (22):
W = | q S | 2 8 B ˜ 1 | Λ Z + i ϱ | | Λ + Z i ϱ | 2 .
Now, in order to optimize the power conversion efficiency, the last term in the square brackets must be minimized. As in Şentürk and Özdamar [18], this can be done by finding the optimum value of Λ , which can be evaluated by applying zero value to the derivative for the squared-right term inside the brackets of Equation (22) with respect to Λ , and thus obtaining
Λ o p t = | Z i ϱ | = B ˜ 2 + A ˜ + ϱ 2 1 / 2 .
Therefore, the maximum value of extracted work at this condition becomes,
W o p t = | q S | 2 8 B ˜ 1 Λ o p t B ˜ Λ o p t + B ˜ ,
where A ˜ , B ˜ and Λ are function of the angular frequency ω which means that for each wave frequency, the turbine parameter must be altered appropriately to satisfy Equation (23).
Thus, the expression for maximum efficiency is expressed as
η m a x = W o p t W m a x = 2 B ˜ Λ o p t + B ˜ ,
where the maximum hydrodynamic efficiency is bounded by 0 η m a x 1 . From expression (23), it is clear that when the radiation susceptance parameter is zero and the air compressibility term is neglected, it results in Λ o p t = B ˜ , W o p t = W m a x and η m a x = 1 , thus implying that the device has effectively absorbed all of the incident wave energy. In this situation, the free surface inside the OWC chamber is characterized by a piston-like resonant mode and the PTO damping optimization is satisfied [31]. To remain at this condition, the rate of energy extraction must equate to the rate of radiation damping while the internal free surface must remain in a state of resonance. Physically, this requires that the radiated waves, resulting from the oscillatory heave motion of the internal free surface, superpose and cancel the incident and scattered waves; in this instance, the device has thereby captured all of the incident wave energy [16].
Now, as in [28] the non-dimensionalised quantities μ and ν to represent the radiation susceptance and radiation conductance coefficients are defined as
μ = ρ g ω b A ˜ ,
ν = ρ g ω b B ˜ ,
respectively, where the radiation conductance coefficient ν is related to the transfer of energy into the system, while the radiation susceptance coefficient μ to the energy that remains un-captured [19].
Therefore, by substituting these coefficients into Equation (25), the efficiency η m a x is
η m a x = 2 1 + μ ν 2 1 / 2 + 1 ,
which is independent of the incident wave power and only depends on the radiation solution of the volume flux q.

4. Solution

In this section, the BEM is used to solve the BVP in the frequency domain. In order to solve the governing equation together with the appropriated boundary conditions, a quadratic distribution of variables along each element is considered. The integral representation of the solution for the Laplace equation Equation (1) at any point source P ¯ inside the domain Ω in terms of the boundary values of ϕ and ϕ / n is given by
α ( P ¯ ) ϕ ( P ¯ ) + Γ ϕ ( q ¯ ) ψ ( P ¯ , q ¯ ) n q ¯ d Γ q ¯ = Γ ψ ( P ¯ , q ¯ ) ϕ ( q ¯ ) n q ¯ d Γ q ¯ ,
where ϕ is the unknown flow potential; ϕ / n is the derivative of the potential relative to normal unit vector on the boundary Γ ; d Γ is the length of an infinitesimal piece of Γ ; q ¯ an arbitrary point; while ψ and ψ / n are the fundamental solution of Laplace equation and its normal derivative at point q ¯ of the boundary, respectively; and α = θ / 2 π , where θ is the internal angle of the corner in radians [32].
The fundamental solution of Laplace equation is given by
ψ = 1 2 π ln r ,
where r is the distance between the source P ¯ and the arbitrary point q ¯ .
Now, discretizing the boundary into a series of N E elements, Equation (28) can be written as
α i ϕ i + j = 1 N E Γ ϕ ψ n d Γ = j = 1 N E Γ ψ ϕ n d Γ .
In order to define the values of ϕ and ϕ / n on each element, three noded quadratic elements are employed, Figure 4. The variables ϕ and ϕ / n are thus written in terms of interpolation functions, φ ^ 1 , 2 , 3 , which are function of a homogeneous coordinate ξ as follows
ϕ ( ξ ) = φ ^ 1 ϕ 1 + φ ^ 2 ϕ 2 + φ ^ 3 ϕ 3 ,
ϕ ( ξ ) n = φ ^ 1 ϕ 1 n + φ ^ 2 ϕ 2 n + φ ^ 3 ϕ 3 n ,
where the superscript indicates the number of the node, while the interpolation functions are given by
φ ^ 1 = 1 2 ξ ξ 1 ,
φ ^ 2 = 1 2 1 ξ 1 + ξ ,
φ ^ 3 = 1 2 ξ 1 + ξ ,
with the dimensionless coordinate ξ varying from −1 to 1. Now, carrying out the integrals from Equation (30) over an element j, these can be written as
Γ j ϕ ψ n d Γ = Γ j φ ^ 1 ψ n d Γ , Γ j φ ^ 2 ψ n d Γ , Γ j φ ^ 3 ψ n d Γ · ϕ 1 ϕ 2 ϕ 3 ,
and
Γ j ψ ϕ n d Γ = Γ j φ ^ 1 ψ d Γ , Γ j φ ^ 2 ψ d Γ , Γ j φ ^ 3 ψ d Γ · ϕ 1 n ϕ 2 n ϕ 3 n .
Here, it is observed that in order to solve the integrals Equations (33) and (34), the calculation of the Jacobian is required since they are a function of the boundary Γ in the x z plane, while the interpolation functions are a function of ξ . This transformation is given by
d Γ = d x d ξ 2 + d z d ξ 2 d ξ = | J | d ξ ,
where J indicates the Jacobian and can then be substituted into Equations (33) and (34). Additionally, in order to calculate the value of Equation (35), the variation of the x and z coordinates in terms of ξ must be also known. This can be carried out in the same way as the variables ϕ and ϕ / n , with the use of the quadratic interpolation defined by
x = φ ^ 1 x 1 + φ ^ 2 x 2 + φ ^ 3 x 3 ,
z = φ ^ 1 z 1 + φ ^ 2 z 2 + φ ^ 3 z 3 ,
where again the superscript indicates the number of the node. Thus, Equation (30) can be written as
α i ϕ i + j = 1 N E h 1 i j , h 2 i j , h 3 i j · ϕ 1 ϕ 2 ϕ 3 j = j = 1 N E g 1 i j , g 2 i j , g 3 i j · ϕ n 1 ϕ n 2 ϕ n 3 j ,
where
h k i j = 1 1 φ ^ k ( ξ ) ψ n | J | d ξ ,
g k i j = 1 1 φ ^ k ( ξ ) ψ | J | d ξ ,
with k = 1 , 2 and 3 and h k i j and g k i j are estimated by using a Gauss integration method with ten points to account for the quadratic variation of the element geometry, the potential and flux. More details regarding the numerical procedure for solving these integrals can be found in [33,34].
Furthermore, as explained by [33], in order to consider the possibility of having different values of ϕ n at node 3 of an element and at node 1 of the next adjoining element, the fluxes are arranged in a 3 × N E array where a position is held for each nodal value of every element. However, in the case of the potential ϕ , its value is always considered unique in the connection between two elements. Therefore, the values of ϕ can be arranged in an N array, where N is the number of nodes equal to 2 N E for closed boundaries. Thus, Equations (37) can be written as
α i ϕ i + H ^ i 1 , H ^ i N · ϕ 1 ϕ N j = G i 1 , G i N E · ϕ n 1 ϕ n 2 ϕ n 3 1 ϕ n 1 ϕ n 2 ϕ n 3 N E j ,
with H ^ i j being equal to the h 1 term of an element plus the h 3 term of the previous element for odd nodes and equal to the h 2 term of the corresponding element for the central nodes. On the other hand, G i j are 1 × 3 matrices with the elements g 1 i j , g 2 i j , g 3 i j . Therefore, the whole system of equations can be simply written as follows
H Φ = G Φ n ,
where H is a square matrix N × N , Φ is an N × 1 vector, G is a rectangular matrix N × 3 N E and Φ n is an 3 N E × 1 vector.
Furthermore, as previously described by [19], in order to avoid the numerical errors arising from the cases where a very thin front wall of the device is considered, the method of subdomains is used to solve the BVP by applying the BEM separately to each of its regions [32,35]. The domain is then divided into three separate regions which have a common interface boundary on both lateral sides of the front wall as shown in Figure 5a. For each subdomain the following vectors are defined: in region R1
  • Φ 1 1 , Φ n 1 1 Nodal values on Γ 1 of the external boundary.
  • Φ 12 1 , Φ n 12 1 Nodal values on the interface Γ 12 ,
where the superscript denotes the region, while the subscript denotes the corresponding external boundary or interface. The number of the nodal points on Γ 1 and Γ 12 are N 1 and N 12 , respectively. In region R2
  • Φ 2 2 , Φ n 2 2 Nodal values on Γ 2 of the external boundary.
  • Φ 12 2 , Φ n 12 2 Nodal values on the interface Γ 12 .
  • Φ 23 2 , Φ n 23 2 Nodal values on the interface Γ 23 ,
where the number of the nodal points on Γ 2 , Γ 12 and Γ 23 are N 2 , N 12 and N 23 , respectively. In region R3
  • Φ 3 3 , Φ n 3 2 Nodal values on Γ 3 of the external boundary.
  • Φ 23 2 , Φ n 23 2 Nodal values on the interface Γ 23 ,
with N 3 and N 23 being the number of the nodal points on Γ 3 and Γ 23 , respectively.
Since ϕ is unknown on either side of the interfaces Γ 12 and Γ 23 , the number of boundary unknowns in each subdomain is:
  • Region 1: N 1 on Γ 1 and N 12 on Γ 12 .
  • Region 2: N 2 on Γ 2 , N 12 on Γ 12 and N 23 on Γ 23 .
  • Region 3: N 3 on Γ 3 and N 23 on Γ 23 .
On the other hand, for ϕ n , which is defined on the three nodes of each element, it is given by:
  • Region 1: M 1 = 3 N E 1 on Γ 1 and M 12 = 3 N E 12 on Γ 12 .
  • Region 2: M 2 = 3 N E 2 on Γ 2 , M 12 = 3 N E 12 on Γ 12 and M 23 = 3 N E 23 on Γ 23 .
  • Region 3: M 3 = 3 N E 3 on Γ 3 and M 23 = 3 N E 23 on Γ 23 ,
where N E j denotes the number of elements depending on the boundary or interface.
Furthermore, in order to match the regions and to obtain the same number of unknowns and equations, the physical consideration of continuity of the potential and flux at the interfaces Equation (3) should be made. Thus, assuming that the nodes in Γ 12 of R 1 and R 2 , and the nodes in Γ 23 of R 2 and R 3 are in perfect contact, Figure 5b, the following physical consideration at the interfaces can be made:
  • Continuity of the potential: The values of the potential on each side of the interface separating two subdomains must be equal
    Φ 12 1 = Φ 12 2 , Φ 23 2 = Φ 23 3 .
  • Continuity of the flux: The outcoming flux from one subdomain is equal to the incoming flux in the adjacent subdomain. Thus, the flux along the normal of the interface requires
    Φ n 12 1 = Φ n 12 2 , Φ n 23 2 = Φ n 23 3 ,
    where the minus signs in the right hand side of Equation (42) indicate that the two flux vectors at the common interface of adjacent subdomains are in opposite directions.
Therefore, the matrix equation for each boundary subdomain is as follows: for the boundary subdomain R 1
H 1 1 H 12 1 Φ 1 1 Φ 12 1 = G 1 1 G 12 1 Φ n 1 1 Φ n 12 1 ,
while for the boundary subdomain R 2
H 12 2 H 2 2 H 23 2 Φ 12 2 Φ 2 2 Φ 23 2 = G 12 2 G 2 2 G 23 2 Φ n 12 2 Φ n 2 2 Φ n 23 2 ,
and for the boundary subdomain R 3
H 3 3 H 23 3 Φ 3 3 Φ 23 3 = G 3 3 G 23 3 Φ n 3 3 Φ n 23 3 .
Equations (43)–(45) of the three subdomains may then be combined in a single matrix equation as
H 1 1 H 12 1 0 0 0 0 0 0 0 H 2 2 H 12 2 H 23 2 0 0 0 0 0 0 0 H 3 3 H 23 3 Φ 1 1 Φ 12 1 Φ 2 2 Φ 12 2 Φ 23 2 Φ 3 3 Φ 23 3 = G 1 1 G 12 1 0 0 0 0 0 0 0 G 2 2 G 12 2 G 23 2 0 0 0 0 0 0 0 G 3 3 G 23 3 Φ n 1 1 Φ n 12 1 Φ n 2 2 Φ n 12 2 Φ n 23 2 Φ n 3 3 Φ n 23 3 ,
where the left-hand side matrix and vector have dimensions of N × N and N × 1 , respectively, with N = N 1 + N 2 + N 3 + 2 N 12 + 2 N 23 , while the right-hand side matrix and vector of N × 3 N E and 3 N E × 1 , respectively, with 3 N E = M 1 + M 2 + M 3 + 2 M 12 + 2 M 23 . Now, since the nodes on Γ 12 and Γ 23 are in perfect contact, the matrices can be further arranged by combining the coefficients of the related variables as follows
H 1 1 H 12 1 0 0 0 0 0 0 H 12 2 H 2 2 0 H 23 2 0 0 0 0 0 0 H 23 3 H 3 3 0 Φ 1 1 Φ 12 1 Φ 2 2 0 Φ 23 2 Φ 3 3 0 = G 1 1 G 12 1 0 0 0 0 0 0 G 12 2 G 2 2 0 G 23 2 0 0 0 0 0 0 G 23 3 G 3 3 0 Φ n 1 1 Φ n 12 1 Φ n 2 2 0 Φ n 23 2 Φ n 3 3 0 .
Finally, after inserting the boundary conditions specified in Equations (2), (6), (8) and (9), and shifting the known variables to the right-hand side and the unknowns to the left-hand side, a matrix of the following form is obtained
A { X } = { B } ,
where { X } is a vector consisting of all the unknown values on the external boundary and on the interfaces of dimension N × 1 ; A is known square coefficient matrix of dimensions N × N whose columns are columns of H and columns of G after a change of sign or sum of two consecutive columns of G with the opposite sign when the unknown is the unique flux at the interfaces at a node connecting two elements [33]; while { B } is a known vector of dimension N × 1 .

5. Results and Discussion

In this section, numerical results based on the BEM discussed in the previous section are presented. First, by considering a thick front wall and ignoring the influence of the air compressibility (i.e., ϱ = 0 ), the effect of the chamber configuration on the hydrodynamic efficiency, radiation susceptance and radiation conductance coefficients is analyzed. Then, based on the physical dimensions of a single chamber of the MWEP, the effects of three different bottom profiles and the air compressibility on the efficiency are studied. It should be mentioned here that in order to minimize the effect of local disturbances at the far-field boundary, for the boundary discretization the distance between the front wall and the far-field boundary was considered to be 4 times the water depth h.
On the other hand, before performing the rest of the numerical calculations, a convergence analysis was carried out. In Table 1, the results of the hydrodynamic efficiency η , radiation susceptance μ and radiation conductance ν for four different K h values are given. It is observed that around 480 nodes (240 quadratic elements) are enough to ensure convergence of the numerical results within three decimal places and also to avoid numerical instabilities that arise when the front wall thickness w tends to zero. Therefore, in the present calculations, all the BVPs are discretized through 480 nodes.

5.1. Front Wall Thickness

To validate the numerical method described here, the numerical results for the limiting case of Evans and Porter [28] were used. In the case of h a / h = 0.125 , b / h = 1 , while w and h e tending to zero, the efficiency obtained by the present formulation was compared with the corresponding results of [28] for an OWC device with a horizontal topography, as shown in Figure 6a. The circles in Figure 6a depict the results of [28], whereas the line represents the results calculated by the present method. It can be seen that both results are in good agreement.
The numerical results of the efficiency η m a x versus K h for different thickness ratios w / b ( = 2.0 , 1.5 , 1.0 , 0.5 , and 0.01 ) in a flat bottom (without considering h e ) are shown in Figure 6a. In this figure, it is seen that by increasing the thickness of the front barrier, the bandwidth of the efficiency curves is reduced and their first peak frequency value is shifted to lower values of the non-dimensional frequency K h . This reduction in the efficiency is explained by the fact that the energy transfer due to the wave motion over small periods is reduced when the front barriers are greater in thickness. However, in a real scenario, during severe storm events, or during times of high water levels, the front barrier is subjected to high loads, due to direct wave action [36,37,38,39], and a slender front wall cannot offer protection to the whole system, as occurred at the MWEP [12]. Therefore, special consideration should be given to this structural aspect.
On the other hand, Figure 6b shows the effect on efficiency of different submergance ratios h a / h ( = 0.125 , 0.250 , 0.500 , and 0.750 ) , together with a front barrier of the same thickness as the OWC chamber (i.e., w / b = 1.0 ). In this figure, it is observed that the effective area of efficiency, under the curve, and the magnitude of the first natural frequency, both increase when the front wall h a / h decreases. Nevertheless, for a relatively small h a / h , and considering that changes in water depth due to tidal variations and extreme waves may take place, a small draft may mean that the trough of a wave propagates below the front wall. This should be considered at the design stage, since in this situation the pressure within the chamber would be equivalent to the atmospheric pressure, causing the power available inside the OWC device to be zero and thus decreasing the efficiency.
The variation of the radiation susceptance and radiation conductance coefficients versus K h when h a / h = 0.125 and b / h = 1.0 for different values of w / b ( = 2.0 , 1.5 , 1.0 , 0.5 , and 0.01 ) is shown in Figure 7a,b, respectively. In Figure 7a it is observed that the frequencies for which μ is zero are related to the peaks of maximum efficiency shown in Figure 6a, which is also evident from the maximum efficiency Equation (27). Figure 7a also shows that when the thickness ratio w / b increases, the range of the non-dimensional frequency K h for which μ is negative also increases. As previously mentioned, the radiation susceptance μ is related to A ˜ which is analogous to the added mass of the forced oscillation of a rigid body system. In this sense, the above-mentioned negative values in μ are may be due to the relevance of the free-surface effects of the internal free surface enclosed by the rigid wall and the surface-piercing front barrier. Negative values in the added mass are a common phenomenon in the theory of submerged and floating moving bodies in a fluid, such as when one or more elements of a structure enclose a portion of the free surface or two-dimensional cylinders are close to the free surface [40,41]. Furthermore, as demonstrated by Falnes [42], negative added mass occurs when an oscillating body produces a water motion where the associated potential energy is greater than the associated kinetic energy.
Figure 7b shows that the peaks in the radiation conductance ν are associated with those observed in Figure 6a. Together with Figure 7a, this shows that by increasing the radiation conductance coefficient ν , with respect to the radiation susceptance coefficient μ , an increase in the power extraction capacity can be attained. In this sense, the radiation conductance coefficient ν indicates the degree to which the system absorbs energy at different frequencies.
Figure 7c,d show the numerical results of the radiation susceptance and radiation conductance coefficients, respectively, versus K h with w / b = 1.0 , b / h = 1.0 and for different submergance ratios h a / h ( = 0.125 , 0.250 , 0.500 , and 0.750 ) . In Figure 7c it is observed that by increasing the front wall draft h a with respect to h, the value of K h for which μ first become zero decreases, which is associated with the first resonance frequency inside the chamber, while the range of K h for negative values of μ increases. On the other hand, Figure 7d shows that a larger draft decreases the frequency at which this resonance occurs. It can also be observed that the peak resonance become more prominent the further the barrier is submerged. With a large draft, conditions are similar to those in a closed tank with parallel sides, where a second resonance mechanism occurs when the incident wave frequency is such that the fluid inside the chamber is excited into an antisymmetric sloshing mode [28]. In this case, the sloshing frequencies occur at values of the dimensionless wave number k b = n π , with n being the sloshing mode. For the case presented in Figure 7d, it is observed that the second peaks in ν due to the first sloshing frequency take place close to K h π .
Figure 8a plots the numerical results of the efficiency η m a x versus K h for various wall to front barrier spacing ratios b / h ( = 1.0 , 0.5 , 0.25 , and 0.125 ) with h a / h = 0.125 and w / b = 1.0 with in a flat bottom. As reported by [19,28], large motions inside the chamber occur when the fluid between the front barrier and the back wall is excited by the incident wave into a resonant piston-like motion inside the OWC. An estimation of this natural frequency of oscillation can be obtained for small values of b / h and so the water contained between the walls can be regarded as a solid body. By employing simple hydrostatic modelling gives that the expected resonance occurs at K h h / h a . In the present case, this resonance would occur at K h = 8 , which is seemed to be approached for smaller values of b / h . On the other hand, for a longer chamber, the frequency at which this resonance occurs is smaller. Physically, this is due to the fact that by increasing the length of the device, the horizontal distance a typical fluid particle must travel during a period of motion increases. This can also be obtained in the vertical direction by increasing the draft of the front wall. As a consequence, a decrease in the value of K h at which resonance occurs is caused, and since an increase in b / h allows more local fluid motion inside the chamber, this leads to a breakdown in the solid-body model of resonance, and the amplitude of oscillation decreases.
Figure 8b illustrates the condition when a step, as long as the front wall, is placed below the latter with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 . First, a comparison of the present formulation with the limiting case of Şentürk and Özdamar [18] of an OWC device with a gap in its fully submerged thin front wall, together with the nondimensional parameters h a / h = 0.125 , h e / h = 0.625 , b / h = 1.0 and w / b = 0.01 , is shown in Figure 8b. It can be seen that both results agree very well. Then, it can be observed that the smaller the distance between the front wall and the step, the lower the magnitude of the non-dimensional frequency K h at which resonance occurs. This is similar to the trend observed when the draft of the front wall is increased without considering the step, Figure 6b. A larger gap leads to a wider range of frequency bandwidth as it increases the transference of energy due to wave motion.
The numerical results for the radiation susceptance and radiation conductance coefficients versus K h when h a / h = 0.125 and w / b = 1.0 for different values of b / h ( = 1.0 , 0.5 , 0.25 , and 0.125 ) are shown in Figure 9a,b, respectively. In Figure 9a, it is observed that when the wall to front barrier spacing is sufficiently small, compared to the depth h, the range of the non-dimensional frequency K h for which the radiation susceptance is negative decreases. Thus, when a small draft is considered, the b / h ratio is important for the occurrence of negative values of the radiation susceptance coefficient. On the other hand, Figure 9b shows that when the length of chamber b increases with respect to the depth h, the radiation conductance coefficient peaks are maximum for longer periods. Consequently, since ν is a measure of the transfer of energy into the system, it may be beneficial to design the chamber length of the OWC device so that the range of frequency bandwidth in the radiation conductance coincides with the most occurring wave period of a particular location and thus exploit the available wave energy as much as possible.
Figure 9c,d show the radiation susceptance and radiation conductance coefficients, respectively, versus K h for various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 . On one hand, Figure 9c shows that the variation from positive to negative in the radiation susceptance increases when the vertical spacing between the step and the front wall decreases. On the other hand, Figure 9d shows an increasing and narrowing peak for the radiation conductance coefficient, as the gap is reduced, decreasing the magnitude of the resonance frequency.
Figure 10a,b show a comparison of the present formulation with the experimental results obtained by Thomas et al. [16] and Wang et al. [43], respectively. First, in the case of a thick front wall considered by Thomas et al. [16], the OWC dimensions employed in the calculations are h = 0.92 m, b = 0.64 m h a = 0.15 m and w = 0.08 m, while in the case of comparison with Wang et al. [43], the dimensions are h = 0.80 m, b = 0.55 m h a = 0.14 m and w = 0.04 m. From Figure 10a,b, it is observed that by comparing the maximum theoretical efficiency and the experimental efficiency, the discrepancy is significant. The numerical solutions overpredict the hydrodynamic efficiency since it neglects the wave nonlinearity and the viscous dissipation, but the resonant frequencies and the shapes of the hydrodynamic efficiency curves predicted by the present numerical method agree well with each experiment [20,22]. It should be pointed out that the present formulation is based on the assumption of an ideal fluid and, therefore, viscous effects and flow separation due to the front wall are apparently the main cause of difference between the experimental and numerical results. Furthermore, another factor contributing to this discrepancy may be attributed to the rate of energy extraction modelled by the PTO system and the energy loss through it by viscous dissipation during the experiments.
In this sense, it is worth mentioning that η m a x = 1 theoretically means that the OWC device effectively captures all the incident wave energy, a condition that in practice is not feasible, due to the radiated wave generated by the oscillatory motion of the internal free surface, the scattering waves by the device and various viscous damping previously mentioned [20]. Thus, in a real scenario, a value of η m a x = 1 could never be achieved because of the energy loss through viscous dissipation as the fluid flow interacts with the chamber geometry and the PTO modelling. It should be mentioned that the effect of vortex and flow separation, which occur near the front wall, can be simulated well by introducing an artificial viscous damping term to the dynamic free surface boundary condition inside the OWC chamber and thus, to account for the energy loss due to vortex shedding and flow separation as previously reported by [20,22,43,44,45].

5.2. Bottom Profile

In this subsection, an analysis is made, based on the physical dimensions of the MWEP and on the highest and lowest tidal levels, Figure 11. The influence of different bottom profiles, the linear turbine damping and the linearized air compressibility inside the chamber on the hydrodynamic efficiency are evaluated. However, it should be noted that in this study the effect of the non-linear phenomena that occur in the interaction between the waves, the OWC device and the trapped air inside the chamber, such as viscous flow separation, turbulence, wave breaking and thermodynamic processes, are not taken into account. These aspects may play an important role in the performance of the OWC when variations in the bottom chamber are considered.
The varying bottom profiles proposed are defined in the interval 0 x b and given by:
z = h ( x ) = h a h b x h a , Sloped bottom , h h a b b 2 x 2 h , Elliptical bottom ,
while the cycloidal bottom is given by the parametric equation
x ( θ ^ ) = r ^ θ ^ + sin θ ^ π + b ,
z ( θ ^ ) = r ^ cos θ ^ + 1 h ,
with 0.56416 θ ^ π and r ^ = 1.51759 m.
In Figure 12, the numerical results for 1:1 sloped bottom profile with different linear turbine damping coefficients, Λ , are compared with the experimental results obtained by Ashlin et al. [21] for a wave steepness H w / λ varying from 0.0320 to 0.0371 . Regarding the numerical results presented here, it is worth noting that these are independent of the wave height since linear wave theory was employed for the formulation of the BVP. In order to perform the calculations, the OWC dimensions are those used by Ashlin et al. [21] in their experiments; these are h = 0.500 m, h a = 0.200 m, b = 0.300 m, w = 0.012 m, a chamber width of l = 0.471 m and a distance from the internal free surface to the top of the chamber of s = 0.400 m. In order to see better the agreement between the results, the least-squares method was applied to the experimental data presented by [21] to obtain a best-fit second-order polynomial curve. Figure 12 shows that the numerical solution and the experimental results of [21] are in good agreement for a linear turbine damping coefficient Λ = 5 × 10 4 m 4 · s/kg. However, it is observed that by comparing the maximum theoretical efficiency obtained from Equation (27) and the experimental efficiency, the discrepancy is high. Therefore, special attention should be paid to turbine damping, as well as non-linear effects, in order to make an adequate estimation of the power absorption of an OWC device.
Figure 13a–d show the numerical results of the maximum η m a x versus the incoming wave period T for the cases when the HEST and LEST take place and the linearized air compressibility is considered. The range of wave period used in these figures is related to mean wave periods reported by [8]. In Figure 13a,c the maximum efficiency was obtained from Equation (25) with ϱ = 0 , while in Figure 13b,d, the optimum value of the damping coefficient ( Λ o p t ), calculated from Equation (23), was employed to obtain the hydrodynamic efficiency by considering the linearized air compressibility inside the chamber. In all the figures it can be seen that the efficiency band shifts slightly to the right as the bottom of the chamber becomes steeper, generating a slightly higher efficiency for higher periods. Figure 13a,b show that for incoming wave periods of less than 8 s, the flat bottom inside the chamber gives maximum efficiency. As reported by [21], this is due to the higher reflection generated by the curved bottom, together with the energy reflection caused by the front and back walls of the OWC device, for shorter wavelengths. The cycloidal and elliptical bottoms decrease the efficiency because these profiles reduce the section of the entrance of fluid particles, obstructing the waves and leading to a decrease in the internal free surface oscillation which drives the air column. This consequently diminishes the output power. The natural frequency of the system is seen to alter slightly for the different bottom profile configurations. Figure 13c,d show the efficiency in the case of the lowest equinoctial spring tide. As expected, compared with the HEST, the period at which resonance takes place is reduced due to the lower front wall draft, h a = 0.60 m. Furthermore, as observed in Figure 6b, Figure 13c shows that the efficiency effective area under the curve increases with a shorter front wall draft. However, in the case of Mutriku, where significant wave heights of 4 m are common [8], this small draft increases the possibility that the wave trough propagates below the front wall, and thus reduces efficiency. Regarding the effect of the volume of air inside the chamber, compared with Figure 13b, Figure 13d shows that efficiency is significantly reduced when the air volume is larger. Comparing Figure 13d to Figure 13c shows that the period at which resonance occurs is independent of the damping condition and is mostly determined by the natural frequency of the water column.

6. Conclusions

The effects of the chamber configuration of an OWC device on efficiency were numerically analyzed, using the BEM, employing quadratic elements. Comparisons were made of these numerical results with theoretical limiting cases obtained by Evans and Porter [28] and Şentürk and Özdamar [18] for a thin front wall and very good agreement was achieved. Numerical estimates for the hydrodynamic efficiency and the radiation susceptance and radiation conductance coefficients were then obtained for different physical configurations. Experimental results reported by Thomas et al. [16] and Wang et al. [43] were compared with those obtained by the present formulation; resonance frequencies and shapes of the hydrodynamic efficiency curves were found to be in good agreement. The main conclusions drawn from this study are as follows:
  • By increasing the thickness of the front barrier, the bandwidth on the efficiency curves is reduced. This reduction in efficiency could be related to the fact that the transfer of energy from the incoming wave to the internal free surface, due to the orbital wave motion, is reduced for short wave periods when the front barrier is thicker.
  • For a thick front barrier, a further reduction in the efficiency effective area under the curve is obtained when the front wall draft is increased.
  • When the OWC chamber length-water depth ratio b / h is decreased, the period of maximum hydrodynamic efficiency is shorter. Consequently, an OWC chamber in which the range of frequency bandwidth in η m a x coincides with the predominant wave period of a particular location, will mean the available wave power will be made better use of.
  • It was observed that the incorporation of a step below the front wall reduces the bandwidth on the efficiency. This step gives a similar effect as that observed when the draft of the front wall is increased in an OWC with a completely flat bottom.
  • It was also observed that when the wall to front barrier spacing is sufficiently small, compared to the depth, the range of the non-dimensional frequency K h , for which the radiation susceptance coefficient is negative, is significantly reduced.
  • By comparing the maximum theoretical efficiency with the experimental efficiency reported by Ashlin et al. [21] for a wave steepness H / λ varying from 0.0320 to 0.0371 , the discrepancy is seen to be high. Therefore, special attention should be paid to turbine damping, as well as to non-linear effects, in order to make an adequate estimation of the power absorption of an OWC.
  • When sloped, cycloidal or elliptical bottom profiles in a chamber of the MWEP were considered, it was seen that the efficiency band slightly shifts to longer periods, as the bottom of the chamber becomes steeper, generating slightly higher efficiency for longer wavelengths.
  • For small periods, it was found that compared with the flat bottom, the sloped, cycloidal and elliptical bottoms diminish the hydrodynamic efficiency. This is due to the reduction of the part of the chamber entrance for the fluid particles, obstructing the waves and leading to a decrease in the internal free surface oscillation which drives the air column.
  • It was observed that in the case of LEST in the MWEP, the efficiency band becomes wider as the draft is reduced. However, when the air volume inside the chamber is greater, the efficiency is significantly reduced.
  • By comparing the different bottom profiles, it was found that the period in which resonance occurs is almost independent of the bottom geometrical configuration and it is mostly determined by the natural frequency of the water column.
This paper is only a numerical investigation on the hydrodynamic performance of an OWC device based on the linear wave theory. Experimental investigations that include the non-linearities on the air compressibility and the turbine damping should be carried out in the future. Finally, it is hoped that the results of this study may provide valuable information for the clean and efficient harnessing of marine renewable energy.

Author Contributions

Conceptualization, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; methodology, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; software, A.A.M.R.; validation, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; formal analysis, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; investigation, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; resources, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; data curation, A.A.M.R.; writing–original draft preparation, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; writing–review and editing, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; visualization, A.A.M.R.; supervision, J.M.B.I., U.I.E. and R.S.C.; project administration, A.A.M.R., J.M.B.I., U.I.E. and R.S.C.; funding acquisition, J.M.B.I., U.I.E. and R.S.C. All authors have read and agreed to the published version of the manuscript.

Funding

The present research has been developed under the framework of CEMIE-Océano (Mexican Centre for Innovation in Ocean Energy). Project FSE-2014-06-249795 financed by CONACYT-SENER- Sustentabilidad Energética. In addition, the authors would like also to express their gratitude for the funding provided by the UPV/EHU [PPGA20/26 research group] and the Basque Government (IT1314-19 research group).

Acknowledgments

The authors would like to thank the Energy Agency of the Basque Government for the information provided regarding the Mutriku Wave Energy Plant and also the Basque Government through the research group (GIU19/029).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BEMBoundary Element Method
BVPBoundary Value Problem
MWEPMutriku Wave Energy Plant
EVEEnte Vasco de la Energía
HESTHighest equinoctial spring tide
LESTLowest equinoctial spring tide
OWCOscillating Water Column
PTOPower take-off

Nomenclature

Awave amplitude
[ A ] square coefficient matrix
A ˜ radiation susceptance parameter
bchamber length
{ B } vector
B ˜ radiation conductance parameter
c g group velocity
Etotal energy per wave period
ggravitational acceleration
[ G ] submatrix
g k i j coefficient integrals
G ^ i j matrix coefficient
G rectangular matrix
hwater depth
[ H ] submatrix
h a front wall draft
h e step height
h k i j coefficient integrals
H ^ i j matrix coefficient
H square matrix
H w wave height
i = 1 the imaginary unit
| J | Jacobian
kwave number
L b front barrier boundary
L g submerged gap
nnormal unit vector
Nnumber of boundary nodes
N E number of boundary elements
M j number of fluxes defined at the corresponding boundary
pspatial pressure distribution
Ptime-dependent pressure distribution
P ¯ point source
p a atmospheric air pressure
qvolume flux
q ¯ arbitrary point
q R radiated volume flux
q S scattered volume flux
Qtime-dependent volume flux
rdistance between P ¯ and q ¯
r ^ radius of the cycloid curve
S b bottom boundary
S f internal free surface boundary
S i external free surface boundary
ttime
Tincident wave period
V 0 air volume inside the chamber
wfront wall thickness
Wmean work absorbed
W m a x maximum work
W o p t optimum work
xhorizontal axis
zvertical axis
Zcomplex admittance

Greek Letters

α = θ 2 π internal angle parameter
γ specific heat ratio of air
Γ boundary
ζ spatial free surface elevation
η ¯ time-dependent free surface elevation
η m a x maximum hydrodynamic efficiency
θ internal angle between two elements
θ ^ parameter
λ wavelength
Λ linear turbine damping coefficient
Λ o p t optimum linear turbine damping coefficient
μ radiation susceptance coefficient
ν radiation conductance coefficient
ξ homogeneous coordinate
ρ density of water
ϱ air compressibility term
ϕ spatial velocity potential
ϕ R radiated velocity potential
ϕ S scattered velocity potential
Φ time-dependent velocity potential
Φ vector containing the velocity potential values
φ ^ 1 , 2 , 3 interpolation functions
ψ 2D fundamental solution of Laplace equation
ω angular frequency
Ω 2D domain

References

  1. Izquierdo, U.; Adolfo Esteban, G.; Blanco, J.M.; Albaina, I.; Peña, A. Experimental validation of a CFD model using a narrow wave flume. Appl. Ocean Res. 2019, 86, 1–12. [Google Scholar] [CrossRef]
  2. Ahamed, R.; McKee, K.; Howard, I. Advancements of wave energy converters based on power take off (PTO) systems: A review. Ocean Eng. 2020, 204, 107248. [Google Scholar] [CrossRef]
  3. Falcão, A.F.O.; Henriques, J.C.C. Oscillating-water-column wave energy converters and air turbines: A review. Renew. Energy 2016, 85, 1391–1424. [Google Scholar] [CrossRef]
  4. Vicinanza, D.; Di Lauro, E.; Contestabile, P.; Gisonni, C.; Lara, J.L.; Losada, I.J. Review of Innovative Harbor Breakwaters for Wave-Energy Conversion. J. Waterway Port Coast. Ocean Eng. 2019, 145, 03119001. [Google Scholar] [CrossRef]
  5. Polinder, H.; Scuotto, M. Wave energy converters and their impact on power systems. In Proceedings of the International Conference on Future Power Systems, Amsterdam, The Netherlands, 16–18 November 2005. [Google Scholar]
  6. Torre-Enciso, Y.; Ortubia, I.; López De Aguileta, L.I.; Marqués, J. Mutriku Wave Power Plant: From the Thinking out to the Reality. In Proceedings of the 8th European Wave and Tidal Energy Conference (EWTEC), Uppsala, Sweden, 7–10 September 2009. [Google Scholar]
  7. Garrido, A.J.; Otaola, E.; Garrido, I.; Lekube, J.; Maseda, F.J.; Liria, P.; Mader, J. Mathematical Modeling of Oscillating Water Columns Wave-Structure Interaction in Ocean Energy Plants. Math. Prob. Eng. 2015, 2015, 727982. [Google Scholar] [CrossRef] [Green Version]
  8. Ibarra-Berastegi, G.; Sáenz, J.; Ulazia, A.; Serras, P.; Esnaola, G.; Garcia-Soto, C. Electricity production, capacity factor, and plant efficiency index at the Mutriku wave farm (2014–2016). Ocean Eng. 2018, 147, 20–29. [Google Scholar] [CrossRef] [Green Version]
  9. Google Maps 2020. Available online: https://www.google.es/maps (accessed on 10 July 2020).
  10. GeoEuskadi, Infraestructura de Datos Espaciales (IDE) de Euskadi. Available online: http://www.geo.euskadi.eus (accessed on 10 July 2020).
  11. Ente Vasco de la Energía. Available online: http://www.eve.es (accessed on 15 August 2020).
  12. Medina-Lopez, E.; Allsop, W.; Dimakopoulos, A.; Bruce, T. Conjectures on the Failure of the OWC Breakwater at Mutriku. In Proceedings of the Coastal Structures and Solutions to Coastal Disasters Joint Conference, Boston, MA, USA, 9–11 September 2015. [Google Scholar]
  13. Heath, T.V. A review of oscillating water columns. Philos. Trans. R. Soc. A 2012, 370, 235–245. [Google Scholar] [CrossRef] [Green Version]
  14. Delmonte, N.; Barater, D.; Giuliani, F.; Cova, P.; Buticchi, G. Review of Oscillating Water Column Converters. IEEE Trans. Ind. Appl. 2016, 52, 1698–1710. [Google Scholar] [CrossRef]
  15. Wang, D.J.; Katory, M.; Li, Y.S. Analytical and experimental investigation on the hydrodynamic performance of onshore wave-power devices. Ocean Eng. 2002, 29, 871–885. [Google Scholar] [CrossRef]
  16. Morris-Thomas, M.T.; Irvin, R.J.; Thiagarajan, K.P. An Investigation Into the Hydrodynamic Efficiency of an Oscillating Water Column. J. Offshore Mech. Arct. 2006, 129, 273–278. [Google Scholar] [CrossRef]
  17. Martins-rivas, H.; Mei, C.C. Wave power extraction from an oscillating water column along a straight coast. Ocean Eng. 2009, 36, 426–433. [Google Scholar] [CrossRef]
  18. Şentürk, U.; Özdamar, A. Wave energy extraction by an oscillating water column with a gap on the fully submerged front wall. Appl. Ocean Res. 2012, 37, 174–182. [Google Scholar] [CrossRef]
  19. Rezanejad, K.; Bhattacharjee, J.; Guedes Soares, C. Stepped sea bottom effects on the efficiency of nearshore oscillating water column device. Ocean Eng. 2013, 70, 25–38. [Google Scholar] [CrossRef]
  20. Ning, D.Z.; Shi, J.; Zou, Q.P.; Teng, B. Investigation of hydrodynamic performance of an OWC (oscillating water column) wave energy device using a fully nonlinear HOBEM (higher-order boundary element method). Energy 2015, 83, 177–188. [Google Scholar] [CrossRef]
  21. John Ashlin, S.; Sundar, V.; Sannasiraj, S.A. Effects of bottom profile of an oscillating water column device on its hydrodynamic characteristics. Renew. Energy 2016, 96, 341–353. [Google Scholar] [CrossRef]
  22. Ning, D.Z.; Wang, R.Q.; Zou, Q.P.; Teng, B. An experimental investigation of hydrodynamics of a fixed OWC Wave Energy Converter. Appl. Energy 2016, 168, 636–648. [Google Scholar] [CrossRef]
  23. Zheng, S.; Zhang, Y.; Iglesias, G. Coast/breakwater-integrated OWC: A theoretical model. Mar. Struct. 2019, 66, 121–135. [Google Scholar] [CrossRef]
  24. Zheng, S.; Alessandro, A.; Zhang, Y.; Greaves, D.; Miles, J.; Iglesias, G. Wave power extraction from multiple oscillating water columns along a straight coast. J. Fluid Mech. 2019, 878, 445–480. [Google Scholar] [CrossRef]
  25. Zheng, S.; Zhu, G.; Simmonds, D.; Greaves, D.; Iglesias, G. Wave power extraction from a tubular structure integrated oscillating water column. Renew. Energy 2020, 150, 342–355. [Google Scholar] [CrossRef]
  26. Koley, S.; Trivedi, K. Mathematical modeling of oscillating water column wave energy converter devices over the undulated sea bed. Eng. Anal. Bound. Elem. 2020, 117, 26–40. [Google Scholar] [CrossRef]
  27. Belibassakis, K.; Magkouris, A.; Rusu, E. A BEM for the Hydrodynamic Analysis of Oscillating Water Column Systems in Variable Bathymetry. Energies 2020, 13, 3403. [Google Scholar] [CrossRef]
  28. Evans, D.V.; Porter, R. Hydrodynamic characteristics of an oscillating water column device. Appl. Ocean Res. 1995, 17, 155–164. [Google Scholar] [CrossRef]
  29. Evans, D.V. Wave-power absorption by systems of oscillating surface pressure distributions. J. Fluid Mech. 1982, 114, 481–499. [Google Scholar] [CrossRef]
  30. Rezanejad, K.; Guedes Soares, C.; López, I.; Carballo, R. Experimental and numerical investigation of the hydrodynamic performance of an oscillating water column wave energy converter. Renew. Energy 2017, 106, 1–16. [Google Scholar] [CrossRef]
  31. Michele, S.; Renzi, E.; Perez-Collazo, C.; Greaves, D.; Iglesias, G. Power extraction in regular and random waves from an OWC in hybrid wind-wave energy systems. Ocean Eng. 2019, 191, 106519. [Google Scholar] [CrossRef]
  32. Katsikadelis, J. Boundary Elements. Theory and Applications; Elsevier: Amsterdam, The Netherlands, 2002. [Google Scholar]
  33. Dominguez, J. Boundary Elements in Dynamics; Computational Mechanics Publications: Southampton, UK, 1993. [Google Scholar]
  34. Brebbia, C.A.; Dominguez, J. Boundary Elements: An Introductory Course; Computational Mechanics Publications: Southampton, UK, 1992. [Google Scholar]
  35. Becker, A. The Boundary Element Method in Engineering: A Complete Course; McGraw-Hill Book Company: London, UK, 1992. [Google Scholar]
  36. Viviano, A.; Naty, S.; Foti, E.; Bruce, T.; Allsop, W.; Vicinanza, D. Large-scale experiments on the behaviour of a generalised Oscillating Water Column under random waves. Renew. Energy 2016, 99, 875–887. [Google Scholar] [CrossRef] [Green Version]
  37. Pawitan, K.A.; Dimakopoulos, A.S.; Vicinanza, D.; Allsop, W.; Bruce, T. A loading model for an OWC caisson based upon large-scale measurements. Coast. Eng. 2019, 145, 1–20. [Google Scholar] [CrossRef] [Green Version]
  38. Viviano, A.; Musumeci, R.E.; Vicinanza, D.; Foti, E. Pressures induced by regular waves on a large scale OWC. Coast. Eng. 2019, 152, 103528. [Google Scholar] [CrossRef]
  39. Pawitan, K.A.; Vicinanza, D.; Allsop, W.; Bruce, T. Front Wall and In-Chamber Impact Loads on a Breakwater-Integrated Oscillating Water Column. J. Waterway Port Coast. Ocean Eng. 2020, 146, 04020037. [Google Scholar] [CrossRef]
  40. McIver, P.; Evans, D.V. The occurrence of negative added mass in free-surface problems involving submerged oscillating bodies. J. Eng. Math. 1984, 18, 7–22. [Google Scholar] [CrossRef]
  41. McIver, M.; McIver, P. The added mass for two-dimensional floating structures. Wave Motion 2016, 64, 1–12. [Google Scholar] [CrossRef] [Green Version]
  42. Falnes, J. Added-Mass Matrix and Energy Stored in the “Near Field”; Internal Report Division of Experimental Physics, University of Trondheim: Trondheim, Norway, 1983. [Google Scholar]
  43. Wang, R.Q.; Ning, D.Z.; Zhang, C.W.; Zou, Q.P.; Liu, Z. Nonlinear and viscous effects on the hydrodynamic performance of a fixed OWC wave energy converter. Coast. Eng. 2018, 131, 42–50. [Google Scholar] [CrossRef]
  44. Ning, D.Z.; Wang, R.Q.; Gou, Y.; Zhao, M.; Teng, B. Numerical and experimental investigation of wave dynamics on a land-fixed OWC device. Energy 2016, 115, 326–337. [Google Scholar] [CrossRef]
  45. Wang, R.Q.; Ning, D.Z. Dynamic analysis of wave action on an OWC wave energy converter under the influence of viscosity. Renew. Energy 2020, 150, 578–588. [Google Scholar] [CrossRef]
Figure 1. Location and turbo-generators of the MWEP. (a) Location of the harbour at Mutriku (source: Google Maps [9]). (b) Bird’s eye view of Mutriku harbour and the OWC-breakwater system. (source: geoEuskadi [10]). (c) Mutriku Wells Turbo-generators (source: EVE [11]).
Figure 1. Location and turbo-generators of the MWEP. (a) Location of the harbour at Mutriku (source: Google Maps [9]). (b) Bird’s eye view of Mutriku harbour and the OWC-breakwater system. (source: geoEuskadi [10]). (c) Mutriku Wells Turbo-generators (source: EVE [11]).
Jmse 08 00751 g001
Figure 2. Comparison between the original and the present-day design of the integrated breakwater-OWC system in Mutriku. (a) View of the MWEP in 2008 (Reproduced with permission from [6]). (b) Mutriku OWC (source: EVE [11]).
Figure 2. Comparison between the original and the present-day design of the integrated breakwater-OWC system in Mutriku. (a) View of the MWEP in 2008 (Reproduced with permission from [6]). (b) Mutriku OWC (source: EVE [11]).
Jmse 08 00751 g002
Figure 3. Definition sketch of an OWC device with a thick front wall.
Figure 3. Definition sketch of an OWC device with a thick front wall.
Jmse 08 00751 g003
Figure 4. Schematic diagram of the modeling of the boundary with quadratic elements.
Figure 4. Schematic diagram of the modeling of the boundary with quadratic elements.
Jmse 08 00751 g004
Figure 5. Boundary element discretization of the composite domain. (a) Separation of the BVP into three regions. (b) Discretization of the three regions.
Figure 5. Boundary element discretization of the composite domain. (a) Separation of the BVP into three regions. (b) Discretization of the three regions.
Jmse 08 00751 g005
Figure 6. Hydrodynamic efficiency versus K h . (a) For various thickness ratios w / b with h a / h = 0.125 and b / h = 1.0 . (b) For different submergence ratios h a / h with w / b = 1.0 and b / h = 1.0 .
Figure 6. Hydrodynamic efficiency versus K h . (a) For various thickness ratios w / b with h a / h = 0.125 and b / h = 1.0 . (b) For different submergence ratios h a / h with w / b = 1.0 and b / h = 1.0 .
Jmse 08 00751 g006
Figure 7. The radiation susceptance and radiation conductance coefficients versus K h . (a) The radiation susceptance coefficient for various thickness ratios w / b with h a / h = 0.125 and b / h = 1.0 . (b) The radiation conductance coefficient for various thickness ratios w / b with h a / h = 0.125 and b / h = 1.0 . (c) The radiation susceptance coefficient for different submergance ratios h a / h with w / b = 1.0 and b / h = 1.0 . (d) The radiation conductance coefficient for different submergance ratios h a / h with w / b = 1.0 and b / h = 1.0 .
Figure 7. The radiation susceptance and radiation conductance coefficients versus K h . (a) The radiation susceptance coefficient for various thickness ratios w / b with h a / h = 0.125 and b / h = 1.0 . (b) The radiation conductance coefficient for various thickness ratios w / b with h a / h = 0.125 and b / h = 1.0 . (c) The radiation susceptance coefficient for different submergance ratios h a / h with w / b = 1.0 and b / h = 1.0 . (d) The radiation conductance coefficient for different submergance ratios h a / h with w / b = 1.0 and b / h = 1.0 .
Jmse 08 00751 g007
Figure 8. Hydrodynamic efficiency versus K h . (a) For various wall to front barrier spacing ratios b / h with h a / h = 0.125 and w / b = 1.0 . (b) For various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 .
Figure 8. Hydrodynamic efficiency versus K h . (a) For various wall to front barrier spacing ratios b / h with h a / h = 0.125 and w / b = 1.0 . (b) For various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 .
Jmse 08 00751 g008
Figure 9. The radiation susceptance and radiation conductance coefficients versus K h . (a) The radiation susceptance coefficient for various wall to front barrier spacing ratios b / h with h a / h = 0.125 and w / b = 1.0 . (b) The radiation conductance coefficient for various wall to front barrier spacing ratios b / h with h a / h = 0.125 and w / b = 1.0 . (c) The radiation susceptance coefficient for various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 . (d) The radiation conductance coefficient for various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 .
Figure 9. The radiation susceptance and radiation conductance coefficients versus K h . (a) The radiation susceptance coefficient for various wall to front barrier spacing ratios b / h with h a / h = 0.125 and w / b = 1.0 . (b) The radiation conductance coefficient for various wall to front barrier spacing ratios b / h with h a / h = 0.125 and w / b = 1.0 . (c) The radiation susceptance coefficient for various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 . (d) The radiation conductance coefficient for various step to bottom ratios h e / h with h a / h = 0.125 , b / h = 1.0 and w / b = 1.0 .
Jmse 08 00751 g009
Figure 10. Hydrodynamic efficiency versus (a) K h and (b) k h . (a) Comparison with experimental results obtained by Thomas et al. [16] for a thick front wall. (b) Comparison with experimental and numerical results obtained by Wang et al. [43].
Figure 10. Hydrodynamic efficiency versus (a) K h and (b) k h . (a) Comparison with experimental results obtained by Thomas et al. [16] for a thick front wall. (b) Comparison with experimental and numerical results obtained by Wang et al. [43].
Jmse 08 00751 g010
Figure 11. Definition sketch of a single chamber in the Mutriku Wave Energy Plant (Dimensions in meters. Chamber width = 4.5 m; HEST = Highest equinoctial spring tide and LEST = Lowest equinoctial spring tide).
Figure 11. Definition sketch of a single chamber in the Mutriku Wave Energy Plant (Dimensions in meters. Chamber width = 4.5 m; HEST = Highest equinoctial spring tide and LEST = Lowest equinoctial spring tide).
Jmse 08 00751 g011
Figure 12. Hydrodynamic efficiency versus period T for different values of linear turbine damping coefficient Λ , without considering the step.
Figure 12. Hydrodynamic efficiency versus period T for different values of linear turbine damping coefficient Λ , without considering the step.
Jmse 08 00751 g012
Figure 13. Hydrodynamic efficiency in the cases of HEST and LEST versus T for different bottom profiles. (a) Maximum efficiency versus T with ϱ = 0 for various bottom profiles with h = 7.90 , h a / h = 0.646 , b / h = 0.392 and w / b = 2.145 in the case of the HEST. (b) Maximum efficiency versus T for various bottom profiles with h = 7.90 , h a / h = 0.646 , b / h = 0.392 and w / b = 2.145 . in the case of the HEST. (c) Maximum efficiency versus T with ϱ = 0 for various bottom profiles with h = 3.40 , h a / h = 0.176 , b / h = 0.392 and w / b = 2.145 in the case of the LEST. (d) Maximum efficiency versus T for various bottom profiles with h = 3.40 , h a / h = 0.176 , b / h = 0.392 and w / b = 2.145 in the case of the LEST.
Figure 13. Hydrodynamic efficiency in the cases of HEST and LEST versus T for different bottom profiles. (a) Maximum efficiency versus T with ϱ = 0 for various bottom profiles with h = 7.90 , h a / h = 0.646 , b / h = 0.392 and w / b = 2.145 in the case of the HEST. (b) Maximum efficiency versus T for various bottom profiles with h = 7.90 , h a / h = 0.646 , b / h = 0.392 and w / b = 2.145 . in the case of the HEST. (c) Maximum efficiency versus T with ϱ = 0 for various bottom profiles with h = 3.40 , h a / h = 0.176 , b / h = 0.392 and w / b = 2.145 in the case of the LEST. (d) Maximum efficiency versus T for various bottom profiles with h = 3.40 , h a / h = 0.176 , b / h = 0.392 and w / b = 2.145 in the case of the LEST.
Jmse 08 00751 g013
Table 1. Values of hydrodynamic efficiency η m a x , radiation susceptance coefficient μ and radiation conductance coefficient ν computed for different number of nodes N with h a / h = 0.125 , b / h = 1.0 and w / b = 0.5 .
Table 1. Values of hydrodynamic efficiency η m a x , radiation susceptance coefficient μ and radiation conductance coefficient ν computed for different number of nodes N with h a / h = 0.125 , b / h = 1.0 and w / b = 0.5 .
N K h = 3.8329 K h = 2.2657 K h = 1.2054 K h = 0.5074
η μ ν η μ ν η μ ν η μ ν
5600.2808−0.29260.04840.4335−0.35950.10350.8621−0.62870.72990.94250.65071.2787
4800.2814−0.29400.04880.4337−0.35980.10370.8622−0.62950.73120.94250.65191.2806
4000.2822−0.29570.04920.4340−0.36020.10390.8624−0.63050.73290.94240.65341.2830
3280.2833−0.29820.04990.4343−0.36080.10420.8626−0.63180.73520.94230.65561.2861
2560.2848−0.30180.05080.4349−0.36160.10460.8629−0.63380.73860.94220.65871.2906
2000.2856−0.30710.05190.4370−0.36440.10610.8636−0.63730.74510.94180.66491.2978

Share and Cite

MDPI and ACS Style

Medina Rodríguez, A.A.; Blanco Ilzarbe, J.M.; Silva Casarín, R.; Izquierdo Ereño, U. The Influence of the Chamber Configuration on the Hydrodynamic Efficiency of Oscillating Water Column Devices. J. Mar. Sci. Eng. 2020, 8, 751. https://doi.org/10.3390/jmse8100751

AMA Style

Medina Rodríguez AA, Blanco Ilzarbe JM, Silva Casarín R, Izquierdo Ereño U. The Influence of the Chamber Configuration on the Hydrodynamic Efficiency of Oscillating Water Column Devices. Journal of Marine Science and Engineering. 2020; 8(10):751. https://doi.org/10.3390/jmse8100751

Chicago/Turabian Style

Medina Rodríguez, Ayrton Alfonso, Jesús María Blanco Ilzarbe, Rodolfo Silva Casarín, and Urko Izquierdo Ereño. 2020. "The Influence of the Chamber Configuration on the Hydrodynamic Efficiency of Oscillating Water Column Devices" Journal of Marine Science and Engineering 8, no. 10: 751. https://doi.org/10.3390/jmse8100751

APA Style

Medina Rodríguez, A. A., Blanco Ilzarbe, J. M., Silva Casarín, R., & Izquierdo Ereño, U. (2020). The Influence of the Chamber Configuration on the Hydrodynamic Efficiency of Oscillating Water Column Devices. Journal of Marine Science and Engineering, 8(10), 751. https://doi.org/10.3390/jmse8100751

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop