Hydrodynamic Investigation of a Concentric Cylindrical OWC Wave Energy Converter

A fixed, concentric, cylindrical oscillating water column (OWC) wave energy converter (WEC) is proposed for shallow offshore sites. Compared with the existing shoreline OWC device, this wave energy device is not restricted by the wave directions and coastline geography conditions. Analytical solutions are derived based on the linear potential-flow theory and eigen-function expansion technique to investigate hydrodynamic properties of the device. Three typical free-surface oscillation modes in the chamber are discussed, of which the piston-type mode makes the main contribution to the energy conversion. The effects of the geometrical parameters on the hydrodynamic properties are further investigated. The resonance frequency of the chamber, the power extraction efficiency, and the effective frequency bandwidth of the device is discussed, amongst other topics. It is found that the proposed OWC-WEC device with a lower draft and wider chamber breadth has better power extraction ability.


Introduction
Sustainable energy from ocean waves has drawn people's attention for over 40 years.So far, there have been a large number of wave energy converters proposed [1,2].The energy extraction from waves mainly relies on the hinged mechanical motions [3], individual body oscillations [4], water overtopping to reservoirs [5], or trapped air pockets in oscillating water column (OWC) chambers [6].Compared with those mechanical oscillators, the OWC wave energy converter has no moving parts in the water, which is an attractive feature in terms of survivability and maintenance convenience [7].
Up to now, the OWCs have been the most widely used wave energy converters in practice.The first OWC device deployed into the sea on a large scale was in Japan [8].After that, a number of OWC devices have been proposed and deployed.For example, a breakwater equipped with OWC plant was finished with a 40 kW Wells turbine in the Sea of Japan [9].A small shoreline OWC (75 kW) was built for academic studying at the island of Islay in Scotland.In Guangdong Province of China, a shore-fixed OWC device (100 kW) was completed in 2001 [10].Recently, a new U-type OWC has been being installed in breakwaters, which is expected to show dramatically improved absorption width for fixed installations [11].However, OWCs inserted in breakwaters can be installed in zones with low available energy [12].It should be noticed that the above-mentioned OWC-WECs are often based on shoreline or near-shore regions, which not only have the potential to harm the coastline environment but also are limited in terms of the deployment scale.To reduce the environmental impact and increase the scale of deployment, the wave energy exploitation tends to progress to offshore sites [13].Typical examples are the 110 kW Mighty Whale [14] and the OE Buoy in Galway Bay [15].In this study, a concentric cylindrical OWC-WEC is innovatively designed for offshore sites with multiple wave directions.From economic and risk considerations, a fixed OWC structure is preferred.As shown in Figure 1, the structure can be considered as an upside-down cylindrical bucket supported by a cylinder.A concentric cylindrical chamber is formed.A beam is used to connect the chamber and the cylinder as shown in Figure 1b.The energy is transformed from the wave to the air by an oscillating water column moving up and down in the chamber.As the bidirectional airflow flows through the orifice on the top of the chamber, the turbine is driven to be an electricity generator.Compared with the shoreline OWC devices that can only capture the wave energy from one direction, the application of the present-type wave energy device is not restricted by the wave direction and coastline geography condition.
increase the scale of deployment, the wave energy exploitation tends to progress to offshore sites [13].Typical examples are the 110 kW Mighty Whale [14] and the OE Buoy in Galway Bay [15].
In this study, a concentric cylindrical OWC-WEC is innovatively designed for offshore sites with multiple wave directions.From economic and risk considerations, a fixed OWC structure is preferred.As shown in Figure 1, the structure can be considered as an upside-down cylindrical bucket supported by a cylinder.A concentric cylindrical chamber is formed.A beam is used to connect the chamber and the cylinder as shown in Figure 1b.The energy is transformed from the wave to the air by an oscillating water column moving up and down in the chamber.As the bidirectional airflow flows through the orifice on the top of the chamber, the turbine is driven to be an electricity generator.Compared with the shoreline OWC devices that can only capture the wave energy from one direction, the application of the present-type wave energy device is not restricted by the wave direction and coastline geography condition.At the preliminary study stage, the efficiency and hydrodynamic properties of the structure must be identified before the optimized dimension parameters can be given.The analytic method has been used for the preliminary research of the OWC-WECs efficiently by scholars.For example, Garrett [16] derived the linear diffraction solution of a fistulous suspended cylinder in shallow water.Sarmento and Falcão [17] derived a two dimensional analytical model to investigate the effects of linear and nonlinear power take-off system on hydrodynamic properties of OWC-WEC.Zhu and Mitchell [18] presented a diffraction analytical solution of waves around a fistulous cylinder without thickness.Martin-Rivas and Mei [19] presented a linear theory of an OWC which is installed on a straight coast for studying the effects of the coastline on wave power conversion.Martin-Rivas and Mei [20] simulated a thin cylindrical OWC-WEC, which was installed in a breakwater.It is concluded that air compressibility can broaden the effective bandwidth of extraction efficiency with the specific chamber volume.Konispoliatis [21] developed an analytical model to study the hydrodynamics of an OWC-WECs array.Three-unit arrays are discussed, which can be used as a floating wind-turbine foundation.In the present study, an analytical solution based on linear potential-flow theory is derived for the present OWC-WEC.
The present paper is organized as follows.Section 2 describes the power take-off model and analytical solutions of the hydrodynamic problem.In Section 3, the accuracy and convergence of the analytical model are validated.Then, effects of environmental and geometrical parameters on the resonant wave modes and hydrodynamic properties are studied systematically.Conclusions are summarized in Section 4. At the preliminary study stage, the efficiency and hydrodynamic properties of the structure must be identified before the optimized dimension parameters can be given.The analytic method has been used for the preliminary research of the OWC-WECs efficiently by scholars.For example, Garrett [16] derived the linear diffraction solution of a fistulous suspended cylinder in shallow water.Sarmento and Falcão [17] derived a two dimensional analytical model to investigate the effects of linear and nonlinear power take-off system on hydrodynamic properties of OWC-WEC.Zhu and Mitchell [18] presented a diffraction analytical solution of waves around a fistulous cylinder without thickness.Martin-Rivas and Mei [19] presented a linear theory of an OWC which is installed on a straight coast for studying the effects of the coastline on wave power conversion.Martin-Rivas and Mei [20] simulated a thin cylindrical OWC-WEC, which was installed in a breakwater.It is concluded that air compressibility can broaden the effective bandwidth of extraction efficiency with the specific chamber volume.Konispoliatis [21] developed an analytical model to study the hydrodynamics of an OWC-WECs array.Three-unit arrays are discussed, which can be used as a floating wind-turbine foundation.In the present study, an analytical solution based on linear potential-flow theory is derived for the present OWC-WEC.
The present paper is organized as follows.Section 2 describes the power take-off model and analytical solutions of the hydrodynamic problem.In Section 3, the accuracy and convergence of the analytical model are validated.Then, effects of environmental and geometrical parameters on the resonant wave modes and hydrodynamic properties are studied systematically.Conclusions are summarized in Section 4.

Mathematical Problem
The proposed OWC-WEC structure is fixed and rigid in the sea.The submerged part can be considered as a concentric cylindrical geometry as in Figure 2. In this study, the draft of the chamber is d.The water depth is h.The bottom-mounted cylinder and the internal and external surfaces of the chamber have the radius R 1 , R 2 , and R 3 , respectively.A Cartesian coordinate system O-xyz is defined with the origin O on the undisturbed free surface, z-axis pointing vertically upwards, and x-axis coinciding with the direction of incident waves.The structure is axisymmetric about the Oz axis.

Mathematical Problem
The proposed OWC-WEC structure is fixed and rigid in the sea.The submerged part can be considered as a concentric cylindrical geometry as in Figure 2. In this study, the draft of the chamber is d.The water depth is h.The bottom-mounted cylinder and the internal and external surfaces of the chamber have the radius R1, R2, and R3, respectively.A Cartesian coordinate system O-xyz is defined with the origin O on the undisturbed free surface, z-axis pointing vertically upwards, and x-axis coinciding with the direction of incident waves.The structure is axisymmetric about the Oz axis.

Power Take-Off Model
In present study, the air in the chamber of the OWC-WEC is considered to be compressible and motion-isentropic.Additionally, all time-dependent variables in the problem are assumed to be harmonic.Then, the time-independent air pressure is introduced: in which  is the wave angular frequency, t is the time, Re is the real part of a complex variable, P0 is the complex amplitude of the pressure, and 1 i   .A Wells turbine, which is installed at the top of the air chamber, is forced to rotate by oscillating air flow.The mass flux rate of the air across the Wells turbine is proportional to the air pressure in the chamber [17,22], and the relationship between the air mass flux and the turbine characteristics can be expressed as [18] in which q0 is the amplitude of the air volume flux, V0 is the mean chamber volume, a  is the air density, N is the turbine rotational speed in r.p.m., D is diameter of the turbine rotor, c is the sound velocity in air, and K is an empirical coefficient depending on the design of the turbines.
Based on the linear theory, the volume flux in the chamber is the sum of volume fluxes due to wave radiation and diffraction.The wave radiation occurs in the situation when the wave motion is purely caused by the oscillating air pressure of the chamber.The wave diffraction is due to scattering

Power Take-Off Model
In present study, the air in the chamber of the OWC-WEC is considered to be compressible and motion-isentropic.Additionally, all time-dependent variables in the problem are assumed to be harmonic.Then, the time-independent air pressure is introduced: in which ω is the wave angular frequency, t is the time, Re is the real part of a complex variable, P 0 is the complex amplitude of the pressure, and i = √ −1.A Wells turbine, which is installed at the top of the air chamber, is forced to rotate by oscillating air flow.The mass flux rate of the air across the Wells turbine is proportional to the air pressure in the chamber [17,22], and the relationship between the air mass flux and the turbine characteristics can be expressed as [18] in which q 0 is the amplitude of the air volume flux, V 0 is the mean chamber volume, ρ a is the air density, N is the turbine rotational speed in r.p.m., D is diameter of the turbine rotor, c is the sound velocity in air, and K is an empirical coefficient depending on the design of the turbines.
Based on the linear theory, the volume flux in the chamber is the sum of volume fluxes due to wave radiation and diffraction.The wave radiation occurs in the situation when the wave motion is purely caused by the oscillating air pressure of the chamber.The wave diffraction is due to scattering of incident waves when the air pressures inside and outside the chamber are identical.Thus, complex amplitude of the volume flux has two components in which q R and q D are connected with the radiation and diffraction problems, respectively.For a radiation problem, an added mass coefficient C and a damping coefficient B can be obtained [23].
Complex amplitude of the volume flux q R , which is caused by the unit air pressure, can be expressed as By substituting Equations ( 3) and (4) into Equation ( 2), the air pressure in the chamber can be expressed as follows in which q D is solved from the wave diffraction problem, and B and C can be obtained by solving the radiation problem with unit forced pressure.The time-averaged value of the power captured by the turbine is [18] The power extraction efficiency ξ can be calculated as [18] ξ = kP out 2 ρ a are nondimensional parameters.χ characterizes the turbine, β represents the air compressibility.C g is the group velocity of the incident wave, g is the gravitational acceleration, ρ denotes the water density, and k is the incident wave number.

Solution of Boundary Value Problem
Under the assumption that the fluid is incompressible, inviscid, and flow-irrotational, the proposed problem is considered.The complex velocity potential ϕ(x, y, z, t) around the structure is introduced to describe the properties of the fluid.For small amplitude wave, the potential ϕ satisfies the following linearized boundary value problem  (11) in which φ is a spatial function.It should be highlighted that σ R is used as a switch between the radiation and diffraction problem.For the radiation problem, σ R = 1 is set and the pressure amplitude in chamber is P 0 = 1 in Equation (1).For the diffraction problem, σ R = 0 is set.Based on the axisymmetric geometry of the OWC-WEC, the radiation and diffraction problem can be solved in the cylindrical polar coordinates (r, θ, z) with x = r cos θ and y = r sin θ (12) The fluid domain can be divided into three subdomains, as shown in Figure 2. The external subdomain is Ω 1 with r ≥ R 3 and −h ≤ z ≤ 0, the lower subdomain Ω 2 with R 2 ≤ r ≤ R 3 and −h ≤ z ≤ −d, and the internal subdomain Ω 3 with R 1 ≤ r ≤ R 2 and −h ≤ z ≤ 0. Correspondingly, the potentials in Ω 1 , Ω 2 , and Ω 3 are expressed by φ (1) , φ (2) , and φ (3) , respectively.
The boundary value problem with respect to φ (1) can be founded as follows ∂φ (1)  ∂z = 0 for r ≥ R 3 and z = −h ∂φ (1)  ∂z , for r ≥ R 3 and z = 0 (15) ∂φ (1) /∂r = 0, for r = R 3 and −d ≤ z ≤ 0 The radiation condition is also needed for r → ∞ .For φ (2) , the boundary value problem corresponds to the following The φ (3) could be obtained from the following boundary value problem Energies 2018, 11, 985 In Ω 3 , P 0 is constant, and then the linear transformation can be made in which φ * satisfies the original system ( 23)-( 28) in Ω 3 , and then the inhomogeneous free surface condition ( 25) can be rewritten as

Mathematical Solution
In cylindrical polar coordinate system, the Laplace's equation of φ could be written as Based on the method of separation of variables, φ may be assumed to be φ(r, θ, z) = R(r)Θ(θ)Z(z).Then, Equation (31) can be rewritten by three ordinary differential equations as Both α 2 and β 2 are constants here.
In subdomain Ω 1 , the solution of Equation (32) has the form: in which B n denote constant coefficients and n are vertical eigen-functions.Further, taking the boundary conditions (i.e., Equations ( 14), ( 15), ( 19), ( 20), ( 24) and (30)) into account leads to the expression Z (I) n for each subdomain as and k (2) In each subdomain, the set of vertical eigen-functions is complete and satisfies the orthogonality relations as n dz = 0, and For Θ, with the continuity condition (i.e., Θ(θ) = Θ(θ + 2π)), the solution of Equation ( 33) can be written as in which B m denote constant coefficients and Θ m are circumferential eigen-functions as Equation ( 34) is the modified Bessel's equation with k 0 = −ik.The solution can be expressed as In Ω 1 , the potential of the incident wave can be written as in which A is the wave amplitude, ε m = {1, for m = 0; 2, for m ≥ 1} is the Neumann symbol, and J m is the first-kind of Bessel function of order m.The radiation condition is also satisfied to guarantee that the radiated waves propagate away from the structure.Then, the R mn term can be written as in which B mn is a constant coefficient, and H m and K m are the first-kind of Hankel function and the second-kind of modified Bessel function, respectively, both of order m.
mn has the following expression Energies 2018, 11, 985 8 of 23 in which C mn , D mn also denote constant coefficients, and I m is the first-kind of modified Bessel functions of order m.Note that the following relationship is satisfied for the case of n = 0 In Ω 3 , it has in which A mn , E mn are constant coefficients, and Y m is the second kind of Bessel functions of order m.Further, with the boundary condition at r = R 1 , the following relationship is further satisfied Then, the general solution of φ in each subdomain can be expressed as mn (r)Z ( 1) with P (1) m (kr), n = 0; K m (k n r), n ≥ 1 (55) n r , n ≥ 1 (57) Equations ( 52)-(54) have to satisfy the boundary conditions at r = R 2 and r = R 3 .Additionally, In order to take advantage of the orthogonality relations, the boundary conditions can be rewritten as Energies 2018, 11, 985 By substituting the general solutions of φ into these conditions and truncating each infinite series to M terms, it can lead to the following linear system of equations (2) n Z (2) n Z (1) (3) l dz; (67) By letting m = 0 ∼ M, l = 0 ∼ M, and n = 0 ∼ M in the eigen-function expansions terms, (M + 1) systems of equations can be obtained.For a specific m, the linear equations could be assembled as with X = [A mn , B mn , C mn , D mn ] T for n = 0, 1, . . ., M Once the coefficients A mn , B mn , C mn , and D mn are known, the velocity potential in each domain can be obtained based on Equations (52)-(54).The complex amplitude of volume flux can be calculated by the free-surface integration for the wave diffraction, and for the wave radiation.
Then, the wave loads on the OWC-WEC are given by: in which n = n x , n y , n z is the unit normal vector on the body surface, (x 0 , z 0 ) is the rotational center, f = f x , f y , f z , and p can be calculated from the linearized Bernoulli's equation as The wave loads on the beam shown in Figure 1 can be obtained if S B is replaced with the surface chamber shell in Equations ( 72) and (73).

Convergence and Validation
The convergence of the results with respect to the truncation number M in Equations ( 63)-( 66) is firstly tested.The water depth h, gravitational acceleration g, and water density ρ are used as the base of non-dimensionalization in the following sections.In the first case, the parameters are set as R 1 /h = 0.5, R 2 /h = 1.0,R 3 /h = 1.50, d/h = 0.5, and kh = 1.0.Figure 3 shows the obtained non-dimensional value of added mass B, damping coefficient C, and volume flux Q D with different M. It can be seen that the numerical convergence can be obtained as M ≥ 20.Similar convergence can also be realized for other cases with the same value of M.

Convergence and Validation
The convergence of the results with respect to the truncation number M in Equations ( 63)-( 66) is firstly tested.The water depth h, gravitational acceleration g , and water density  are used as the base of non-dimensionalization in the following sections.In the first case, the parameters are set as R1/h = 0.5, R2/h = 1.0,R3/h = 1.50, d/h = 0.5, and kh = 1.0.Figure 3 shows the obtained non-dimensional value of added mass B, damping coefficient C, and volume flux QD with different M. It can be seen that the numerical convergence can be obtained as 20 M  .Similar convergence can also be realized for other cases with the same value of M.Then, the validation of the present solution is examined.A higher-order BEM (boundary element method) [24] was applied to solve the diffraction problem with the same parameters as above.The Then, the validation of the present solution is examined.A higher-order BEM (boundary element method) [24] was applied to solve the diffraction problem with the same parameters as above.The rotation center is fixed at the bottom of the cylinder.Figure 4 shows the comparisons of wave forces in x and z directions and moment about axis y obtained by the present method and BEM, respectively.From these comparisons, the numerical results given by two methods are in good agreement, which verifies the accuracy of the present analytical model.rotation center is fixed at the bottom of the cylinder.Figure 4 shows the comparisons of wave forces in x and z directions and moment about axis y obtained by the present method and BEM, respectively.From these comparisons, the numerical results given by two methods are in good agreement, which verifies the accuracy of the present analytical model.

Determination of Turbine Rotate Speed
Effect of turbine rotate speed on the conversion efficiency is studied firstly.Martins-Rivas and Mei [18] simulated one OWC-WEC property in Pico Island, Azores, Portugal with N = 2000 r.p.m., K

Determination of Turbine Rotate Speed
Effect of turbine rotate speed on the conversion efficiency is studied firstly.Martins-Rivas and Mei [18] simulated one OWC-WEC property in Pico Island, Azores, Portugal with N = 2000 r.p.m., K = 0.45 (for one turbine).As a preliminary study, the same sea depth and power take-off system were selected to simulate the present OWC-WEC.Figure 5 shows the efficiency distribution of the OWC-WEC, with different turbine rotate speeds N = 100 (r.m.p), N = 200 (r.m.p), N = 500 (r.m.p), N = 1000 (r.m.p), and N = 2000 (r.m.p).The remaining geometrical parameters are kept as R 1 /h = 0.15, R 2 /h = 0.50, R 3 /h = 0.51, d/h = 0.2, K = 0.45, V 0 = πR 2  2 h A/h = 0.1, and h = 10 m.It is found that there is one optimal conversion frequency for different turbine rotation speeds.Additionally, the optimal frequency shifts towards a lower frequency with the decrease of the turbine rotate speed.It is clear that the case of N = 200 (r.m.p) is consistent with the results obtained in Martins-Rivas and Mei's model.This is the result of the difference between the present and Martins-Rivas and Mei's models.Additionally, based on the semi-empirical Pierson-Moscowitz (P-M) energy density spectrum in 10 m sea depth, the wave frequency is in range of 1 ≤ kh ≤ 4. Hence, the case of N = 200 (r.m.p) can be used to simulate the wave energy conversion in the present OWC-WEC.
Energies 2018, 11, x FOR PEER REVIEW 12 of 23 = 0.45 (for one turbine).As a preliminary study, the same sea depth and power take-off system were selected to simulate the present OWC-WEC.Figure 5

Wave Motion Inside the OWC Chamber
To study the wave phenomena inside the OWC chamber, an example is considered with the parameters of R1/h = 0.15, R2/h = 0.35, R3/h = 0.4, d/h = 0.2, K = 0.45, N = 200 r.p.m, and h = 10 m.The incident wave propagates in x-direction.Figure 6 shows the variation of the wave surface amplitude at two points U1(−0.20,0.0) and U2(−0.30,0.0) inside the chamber versus the dimensionless wave number kh.The vertical coordinate is normalized by the incident wave amplitude A. From the figure, it can be seen that both surface amplitudes at different positions vary in a similar trend with the wave frequency.In the low frequency region ( 1.5 kh  ), the dimensionless surface amplitudes inside the chamber are near to 1, which is due to the fact that the chamber size is much smaller than the wave length, and the relating wave effects can be ignored.Then, three peaks occur at kh = 2.83, kh = 4.68, and kh = 8.15 with the increase of the wave frequency.W1, W2, and W3 are used to denote these three situations for convenience.It is also noticed that, compared with condition W1, the dimensionless surface amplitudes at the other two conditions show sharper peaks with singular feature.

Wave Motion Inside the OWC Chamber
To study the wave phenomena inside the OWC chamber, an example is considered with the parameters of R 1 /h = 0.15, R 2 /h = 0.35, R 3 /h = 0.4, d/h = 0.2, K = 0.45, N = 200 r.p.m, V 0 = πR 2 2 h, A/h = 0.1, and h = 10 m.The incident wave propagates in x-direction.Figure 6 shows the variation of the wave surface amplitude at two points U 1 (−0.20,0.0) and U 2 (−0.30, 0.0) inside the chamber versus the dimensionless wave number kh.The vertical coordinate is normalized by the incident wave amplitude A. From the figure, it can be seen that both surface amplitudes at different positions vary in a similar trend with the wave frequency.In the low frequency region (kh < 1.5), the dimensionless surface amplitudes inside the chamber are near to 1, which is due to the fact that the chamber size is much smaller than the wave length, and the relating wave effects can be ignored.Then, three peaks occur at kh = 2.83, kh = 4.68, and kh = 8.15 with the increase of the wave frequency.W1, W2, and W3 are used to denote these three situations for convenience.It is also noticed that, compared with condition W1, the dimensionless surface amplitudes at the other two conditions show sharper peaks with singular feature.Contours of surface oscillation amplitude inside the chamber are plotted in Figure 7 at situations W1, W2, and W3, respectively.Note that all amplitudes are normalized by the incident wave height A. From Figure 7, it can be seen that the wave amplitudes of the free surface are almost the same with each another at W1, which leads to a piston-type water oscillation in the chamber.However, the distributions of the free surface amplitudes at W2 and W3 are nearly symmetrical about x = 0 or y = 0.The phenomena in Figure 7b,c are in sloshing modes occurring at circumferential direction with order m = 1 and m = 2 in Equation (43), respectively.Contours of surface oscillation amplitude inside the chamber are plotted in Figure 7 at situations W1, W2, and W3, respectively.Note that all amplitudes are normalized by the incident wave height A. From Figure 7, it can be seen that the wave amplitudes of the free surface are almost the same with each another at W1, which leads to a piston-type water oscillation in the chamber.However, the distributions of the free surface amplitudes at W2 and W3 are nearly symmetrical about x = 0 or y = 0.The phenomena in Figure 7b,c  Contours of surface oscillation amplitude inside the chamber are plotted in Figure 7 at situations W1, W2, and W3, respectively.Note that all amplitudes are normalized by the incident wave height A. From Figure 7, it can be seen that the wave amplitudes of the free surface are almost the same with each another at W1, which leads to a piston-type water oscillation in the chamber.However, the distributions of the free surface amplitudes at W2 and W3 are nearly symmetrical about x = 0 or y = 0.The phenomena in Figure 7b,c  To further analyze water motion inside the chamber, another four points U3(0.3,0),U4(0.2,0),U5(0,0.2), and U6(0,0.3)are considered as shown in Figure 7. Time series of the water elevations at these points are shown in Figure 8. From Figure 8a, it can be seen that the wave elevations at these four points in the chamber are always in phase and with the similar amplitude, as a feature of the piston-type oscillation.It further proves that the water inside the chamber at condition W1 oscillates as a piston, which can compress and expand the chamber air effectively.Strictly speaking, the wave amplitudes at these positions are not exactly same.This is due to the fact that a small-amplitude sloshing also occurs in the chamber at the same time, even though the piston motion is absolutely dominant at W1.For cases of W2 and W3, the free surface motions in the chamber are either centrosymmetric or antisymmetric about the plane x = 0 in Figure 8b,c.This means that, except for W1, there is no evident variation of the total air volume in the chamber.Figure 9 presents the distribution of power extraction efficiency  with dimensionless wave number kh.It can be seen that the extraction efficiency  reaches the peak at kh = 2.83, which just coincides with condition W1, i.e., the piston-type mode.Hence, it is the piston-type mode that makes the main contribution to the wave energy conversion.To further analyze water motion inside the chamber, another four points U 3 (0.3,0), U 4 (0.2,0), U 5 (0,0.2), and U 6 (0,0.3) are considered as shown in Figure 7. Time series of the water elevations at these points are shown in Figure 8. From Figure 8a, it can be seen that the wave elevations at these four points in the chamber are always in phase and with the similar amplitude, as a feature of the piston-type oscillation.It further proves that the water inside the chamber at condition W1 oscillates as a piston, which can compress and expand the chamber air effectively.Strictly speaking, the wave amplitudes at these positions are not exactly same.This is due to the fact that a small-amplitude sloshing also occurs in the chamber at the same time, even though the piston motion is absolutely dominant at W1.For cases of W2 and W3, the free surface motions in the chamber are either centrosymmetric or antisymmetric about the plane x = 0 in Figure 8b,c.This means that, except for W1, there is no evident variation of the total air volume in the chamber.Figure 9 presents the distribution of power extraction efficiency ξ with dimensionless wave number kh.It can be seen that the extraction efficiency ξ reaches the peak at kh = 2.83, which just coincides with condition W1, i.e., the piston-type mode.Hence, it is the piston-type mode that makes the main contribution to the wave energy conversion.

Effects of Geometry Parameters on Wave Loads and Extraction Efficiency
Both survivability and efficiency are important factors for a WEC.The effects of the OWC geometrical parameters on them are investigated in this section, which can be used as guidance for OWC optimization.To be specific, the chamber draft d, chamber breadth b = R2 − R1, and thickness of the chamber wall dwall = R3 − R2 are taken into account.As the preliminary study, the main parameters are selected based on engineering practice to be with D/h = 0.3, K = 0.45, N = 200 r.p.m, A/h = 0.1, and h = 10 m.Lateral wave force and bending moment are critical parameters for survivability and reliability in ocean engineering.Figure 10a,b shows the wave force fx and bending moment My about sea bed on the whole OWC device with four chamber drafts, i.e., d/h = 0.10, d/h = 0.20, d/h = 0.30, and d/h = 0.40.The chamber breadth b/h = 0.2 and chamber wall thickness dwall/h = 0.05 are fixed.Figure 10c shows the wave force fxbeam of the chamber shell, and Figure 10d shows the chamber shell moment Mybeam about the connection beam between the cylinder and the chamber.The beam is assumed above the water surface D/h = 0.30.In Figure 10a,b, two resonant peaks can be found within the considered range of 0 6 kh   for both the wave force and moment.It is the first peak being of engineering interest.From the figure, it can be seen that as d/h change from 0.10 to 0.40, the first peaks of both the wave force fx and moment My increase greatly due to the increase of the wet surface acting on the OWC-WEC.Hence, lower chamber draft is preferred to decrease the wave loads.In the high frequency domain, another resonant peak occurs.In Figure 10a, for the case of d/h = 0.2, the second peak occurs at kh = 4.68, which is exactly corresponding to condition W2 as shown in Figure 6.At kh = 4.68, the wave length L/h = 0.67, which is close to the chamber inner diameter 2R2 = 0.70.The incident waves are captured in m = 1 mode.Therefore, it can be concluded that the resonant sloshing mode at order m = 1 in Equation ( 43) induces an extreme wave load on the device at the second peak.The potential theory without considering the fluid viscidity leads to an over-prediction of the wave force at resonance.However, in real fluids, the wave loads at the W2 sloshing mode cannot be so large due to the viscous dissipation.Such peak wave force acting on the OWC-WEC increases sharply with the chamber draught.Additionally, the resonant sloshing frequency shifts towards a lower frequency.In other words, the sloshing resonance becomes much stronger with the increase of the chamber draught [16].In Figure 10c,d, the beam force fxbeam and bending moment Mybeam vary with the same trend in Figure 10a,b.Compared with fx and My, wave loads of the beam at the first peak decrease slightly.However, the resonant response at m = 1 mode is much higher than that in Figure 10a,b.Hence, in present OWC-WEC, lower chamber draught leads to higher survivability.

Effects of Geometry Parameters on Wave Loads and Extraction Efficiency
Both survivability and efficiency are important factors for a WEC.The effects of the OWC geometrical parameters on them are investigated in this section, which can be used as guidance for OWC optimization.To be specific, the chamber draft d, chamber breadth b = R 2 − R 1 , and thickness of the chamber wall d wall = R 3 − R 2 are taken into account.As the preliminary study, the main parameters are selected based on engineering practice to be with D/h = 0.3, K = 0.45, N = 200 r.p.m, V 0 = πR 2  2 D, A/h = 0.1, and h = 10 m.
Lateral wave force and bending moment are critical parameters for survivability and reliability in ocean engineering.Figure 10a,b shows the wave force f x and bending moment M y about sea bed on the whole OWC device with four chamber drafts, i.e., d/h = 0.10, d/h = 0.20, d/h = 0.30, and d/h = 0.40.The chamber breadth b/h = 0.2 and chamber wall thickness d wall /h = 0.05 are fixed.Figure 10c shows the wave force f xbeam of the chamber shell, and Figure 10d shows the chamber shell moment M ybeam about the connection beam between the cylinder and the chamber.The beam is assumed above the water surface D/h = 0.30.In Figure 10a,b, two resonant peaks can be found within the considered range of 0 ≤ kh ≤ 6 for both the wave force and moment.It is the first peak being of engineering interest.From the figure, it can be seen that as d/h change from 0.10 to 0.40, the first peaks of both the wave force f x and moment M y increase greatly due to the increase of the wet surface acting on the OWC-WEC.Hence, lower chamber draft is preferred to decrease the wave loads.In the high frequency domain, another resonant peak occurs.In Figure 10a, for the case of d/h = 0.2, the second peak occurs at kh = 4.68, which is exactly corresponding to condition W2 as shown in Figure 6.At kh = 4.68, the wave length L/h = 0.67, which is close to the chamber inner diameter 2R 2 = 0.70.The incident waves are captured in m = 1 mode.Therefore, it can be concluded that the resonant sloshing mode at order m = 1 in Equation ( 43) induces an extreme wave load on the device at the second peak.The potential theory without considering the fluid viscidity leads to an over-prediction of the wave force at resonance.However, in real fluids, the wave loads at the W2 sloshing mode cannot be so large due to the viscous dissipation.Such peak wave force acting on the OWC-WEC increases sharply with the chamber draught.Additionally, the resonant sloshing frequency shifts towards a lower frequency.In other words, the sloshing resonance becomes much stronger with the increase of the chamber draught [16].In Figure 10c,d, the beam force f xbeam and bending moment M ybeam vary with the same trend in Figure 10a,b.Compared with f x and M y , wave loads of the beam at the first peak decrease slightly.However, the resonant response at m = 1 mode is much higher than that in Figure 10a,b.Hence, in present OWC-WEC, lower chamber draught leads to higher survivability.Further, the effects of the chamber draught on power extraction efficiency are presented in Figure 11.A dash line at ξ = 0.3 is marked in the figure to depict the effective frequency bandwidth (i.e., the part for ξ ≥ 0.3).The results indicate that the effective frequency bandwidth decreases largely with the reflected by a larger-draught.An explanation to such a phenomenon is that, in the high-frequency zone, the corresponding short waves with lower penetrability can be easily reflected by larger-draught chamber.From the figure, it also can be observed that the resonant frequency and maximal efficiency go up with the decrease of chamber draught.In other words, the proposed OWC device with a smaller chamber draft tends to possess a better capacity for wave energy conversion.
Further, the effects of the chamber draught on power extraction efficiency are presented in Figure 11.A dash line at =0.3  is marked in the figure to depict the effective frequency bandwidth (i.e., the part for ).The results indicate that the effective frequency bandwidth decreases largely with the reflected by a larger-draught.An explanation to such a phenomenon is that, in the high-frequency zone, the corresponding short waves with lower penetrability can be easily reflected by larger-draught chamber.From the figure, it also can be observed that the resonant frequency and maximal efficiency go up with the decrease of chamber draught.In other words, the proposed OWC device with a smaller chamber draft tends to possess a better capacity for wave energy conversion.Figure 12 shows the effects of the chamber breadth on the hydrodynamic loads with four chamber breadths b/h = 0.10, b/h = 0.20, b/h = 0.30, and b/h = 0.40.The drafts and wall thickness of the chamber are set as constants d/h = 0.20 and dwall/h = 0.05.In Figure 12a,b, there are also two resonant peaks shown in the considered range of kh for both the wave force and moment.As the chamber breadth increases, the wave moment My decreases greatly at the first resonant frequency.However, the wave force fx on the device increases with the chamber breadth due to the increase of the chamber wall surface in domain Ω1 leading to a larger wave force.In the high frequency domain, as b/h increases, the second resonant peak shifts towards a lower frequency and increases its magnitude slightly.For wider chamber breadth, the resonance of the sloshing model at order m = 1 occurs at the longer wave area.Figure 12c,d shows the wave loads of the connection beam with different chamber breadth.It also can be seen that the beam force fxbeam and bending moment Mybeam at the m = 1 mode are exceptionally high.Additionally, compared with fx and My, the frequency bandwidth in W1 sloshing mode in Figure 12c,d increases with the chamber breadth; this is especially apparent for the bending moment Mybeam.Hence, for wider chamber breadth, the bending moment of the whole OWC-WEC decreases, even though the resonant loads of the connection beam increase greatly.Figure 13 shows the power extraction efficiency with different chamber breadths.It can be seen that the effective frequency bandwidth increases sharply with the chamber breadth for capturing longer waves.The reason for this is that the inertia of the water column in the chamber increases with chamber width.Thus, long wave energy captured can make a greater contribution to the piston-type oscillation in the chamber.It is also found that the resonant frequency decreases rapidly with chamber breadth.In other words, the proposed OWC with a wider chamber tends to possess better efficiency.12a,b, there are also two resonant peaks shown in the considered range of kh for both the wave force and moment.As the chamber breadth increases, the wave moment M y decreases greatly at the first resonant frequency.However, the wave force f x on the device increases with the chamber breadth due to the increase of the chamber wall surface in domain Ω 1 leading to a larger wave force.In the high frequency domain, as b/h increases, the second resonant peak shifts towards a lower frequency and increases its magnitude slightly.For wider chamber breadth, the resonance of the sloshing model at order m = 1 occurs at the longer wave area.Figure 12c,d shows the wave loads of the connection beam with different chamber breadth.It also can be seen that the beam force f xbeam and bending moment M ybeam at the m = 1 mode are exceptionally high.Additionally, compared with f x and M y , the frequency bandwidth in W1 sloshing mode in Figure 12c,d increases with the chamber breadth; this is especially apparent for the bending moment M ybeam .Hence, for wider chamber breadth, the bending moment of the whole OWC-WEC decreases, even though the resonant loads of the connection beam increase greatly.Figure 13 shows the power extraction efficiency with different chamber breadths.It can be seen that the effective frequency bandwidth increases sharply with the chamber breadth for capturing longer waves.The reason for this is that the inertia of the water column in the chamber increases with chamber width.Thus, long wave energy captured can make a greater contribution to the piston-type oscillation in the chamber.It is also found that the resonant frequency decreases rapidly with chamber breadth.In other words, the proposed OWC with a wider chamber tends to possess better efficiency.Finally, the effects of the chamber wall thickness on the hydrodynamic loads and extraction efficiency are presented in Figures 14 and 15, respectively.Four different chamber wall thicknesses, i.e., dwall/h = 0.01, dwall/h = 0.05, dwall/h = 0.10, and dwall/h = 0.15, are considered.The other parameters are set as constant d/h = b/h = 0.20.It can be seen that the chamber wall thickness has little impact on the second resonant frequency, which is mainly decided by the chamber breadth.With the increase of the chamber wall thickness, at the first resonant peak, the wave force fx increases, but the bending moment My decreases.An explanation to the latter is that, while increasing the chamber wall thickness, the wet surface of the proposed device in domain Ω2 increases, which can reduce the bending moment My about the sea bed. Figure 15 shows the variation of extraction efficiency with four different thicknesses of the chamber wall.The results indicate that the extraction efficiency  and effective frequency bandwidth decrease as the wall thickness grows.However, compared with the results shown in Figures 11 and 13, the effect of the wall thickness on the resonant frequency is weaker.Therefore, for engineering considerations, the chamber wall thickness can be designed mainly according to extreme environmental loads and material strength.the effects of the chamber wall thickness on the hydrodynamic loads and extraction efficiency are presented in Figures 14 and 15, respectively.Four different chamber wall thicknesses, i.e., d wall /h = 0.01, d wall /h = 0.05, d wall /h = 0.10, and d wall /h = 0.15, are considered.The other parameters are set as constant d/h = b/h = 0.20.It can be seen that the chamber wall thickness has little impact on the second resonant frequency, which is mainly decided by the chamber breadth.With the increase of the chamber wall thickness, at the first resonant peak, the wave force f x increases, but the bending moment M y decreases.An explanation to the latter is that, while increasing the chamber wall thickness, the wet surface of the proposed device in domain Ω 2 increases, which can reduce the bending moment M y about the sea bed. Figure 15 shows the variation of extraction efficiency with four different thicknesses of the chamber wall.The results indicate that the extraction efficiency ξ and effective frequency bandwidth decrease as the wall thickness grows.However, compared with the results shown in Figures 11 and 13, the effect of the wall thickness on the resonant frequency is weaker.Therefore, for engineering considerations, the chamber wall thickness can be designed mainly according to extreme environmental loads and material strength.Finally, the effects of the chamber wall thickness on the hydrodynamic loads and extraction efficiency are presented in Figures 14 and 15, respectively.Four different chamber wall thicknesses, i.e., dwall/h = 0.01, dwall/h = 0.05, dwall/h = 0.10, and dwall/h = 0.15, are considered.The other parameters are set as constant d/h = b/h = 0.20.It can be seen that the chamber wall thickness has little impact on the second resonant frequency, which is mainly decided by the chamber breadth.With the increase of the chamber wall thickness, at the first resonant peak, the wave force fx increases, but the bending moment My decreases.An explanation to the latter is that, while increasing the chamber wall thickness, the wet surface of the proposed device in domain Ω2 increases, which can reduce the bending moment My about the sea bed. Figure 15 shows the variation of extraction efficiency with four different thicknesses of the chamber wall.The results indicate that the extraction efficiency  and effective frequency bandwidth decrease as the wall thickness grows.However, compared with the results shown in Figures 11 and 13, the effect of the wall thickness on the resonant frequency is weaker.Therefore, for engineering considerations, the chamber wall thickness can be designed mainly according to extreme environmental loads and material strength.

Conclusions
This research has provided an easy-to-deploy OWC-WEC design with good performance.This OWC-WEC is featured with an excellent adaption to offshore sites with multiple wave directions.The overall costs are expected to be reduced when a large-scale deployment is achieved.The structure has a mushroom appearance, which can be considered as an upside down cylindrical bucket supported by a column fixed on the seabed.A concentric cylindrical chamber is formed.The energy is transformed from the wave to the air by an oscillating water column moving up and down in the chamber.As the bidirectional airflow flows through the orifice on the top of the chamber, the turbine drives an electricity generator.Compared with the existing shoreline OWC devices, this wave energy device is not restricted by the wave directions and coastline geography conditions.
Then, analytical solutions are derived based on the linear potential-flow theory to investigate hydrodynamic properties of the wave energy converter.The wave motion inside the chamber and the effect of geometry parameters on the hydrodynamic loads and extraction efficiency are investigated.Three typical free-surface oscillation modes in the chambers are found and analyzed.The piston-type mode, which can compress the chamber air efficiently, makes the largest contribution to the energy conversion efficiency.Further, it is found that, as the chamber draft decreases, the resonant frequency and peak efficiency increase apparently, and the wave loads on the whole OWC-WEC also decrease greatly.Additionally, for the design of the connection beam, the wave loads at m = 1 mode should be considered, especially the OWC-WEC with wider chamber breadth.For wider effective frequency bandwidth, the OWC shell can be selected with lower draft and wider chamber breadth.Compared with the chamber draft and breadth, the chamber wall thickness has a relatively small impact on the wave energy conversion, which can be designed by environmental loads and material strength in engineering.In the next stage, a series of model experiments (e.g., similar to those in Viviano et al. [25]) will be carried out for further investigation based on the present analytical solutions.

Figure 1 .
Figure 1.Concept of the concentric cylindrical OWC-WEC, (a) 3D effect picture; (b) the profile of the chamber.

Figure 1 .
Figure 1.Concept of the concentric cylindrical OWC-WEC, (a) 3D effect picture; (b) the profile of the chamber.

Figure 3 .
Figure 3. Convergence study with respect to parameter M, for (a) B, (b) C, and (c) QD.

Figure 3 .
Figure 3. Convergence study with respect to parameter M, for (a) B, (b) C, and (c) Q D .

Figure 4 .
Figure 4. Comparison of hydrodynamic force amplitudes by the present model and the BEM: (a) wave force in x direction, (b) wave force in z direction, and (c) wave moment about axis y.

Figure 4 .
Figure 4. Comparison of hydrodynamic force amplitudes by the present model and the BEM: (a) wave force in x direction, (b) wave force in z direction, and (c) wave moment about axis y.

Figure 5 .
Figure 5. Distribution of hydrodynamic efficiency for different turbine rotational speed N on the efficiency  .

Figure 5 .
Figure 5. Distribution of hydrodynamic efficiency for different turbine rotational speed N on the efficiency ξ.

Figure 6 .
Figure 6.Variations of dimensionless surface amplitudes at two positions with kh.

Figure 6 .
Figure 6.Variations of dimensionless surface amplitudes at two positions with kh.

Figure 6 .
Figure 6.Variations of dimensionless surface amplitudes at two positions with kh.

Figure 8 .
Figure 8.Time series of the free-surface elevation at test positions in the chamber at conditions (a) W1, (b) W2, and (c) W3.Figure 8. Time series of the free-surface elevation at test positions in the chamber at conditions (a) W1, (b) W2, and (c) W3.

Figure 8 .Figure 9 .
Figure 8.Time series of the free-surface elevation at test positions in the chamber at conditions (a) W1, (b) W2, and (c) W3.Figure 8. Time series of the free-surface elevation at test positions in the chamber at conditions (a) W1, (b) W2, and (c) W3.

Figure 12
Figure12shows the effects of the chamber breadth on the hydrodynamic loads with four chamber breadths b/h = 0.10, b/h = 0.20, b/h = 0.30, and b/h = 0.40.The drafts and wall thickness of the chamber are set as constants d/h = 0.20 and d wall /h = 0.05.In Figure12a,b, there are also two resonant peaks shown in the considered range of kh for both the wave force and moment.As the chamber breadth increases, the wave moment M y decreases greatly at the first resonant frequency.However, the wave force f x on the device increases with the chamber breadth due to the increase of the chamber wall surface in domain Ω 1 leading to a larger wave force.In the high frequency domain, as b/h increases, the second resonant peak shifts towards a lower frequency and increases its magnitude slightly.For wider chamber breadth, the resonance of the sloshing model at order m = 1 occurs at the longer wave area.Figure12c,dshows the wave loads of the connection beam with different chamber breadth.It also can be seen that the beam force f xbeam and bending moment M ybeam at the m = 1 mode are exceptionally high.Additionally, compared with f x and M y , the frequency bandwidth in W1 sloshing mode in Figure12c,d increases with the chamber breadth; this is especially apparent for the bending moment M ybeam .Hence, for wider chamber breadth, the bending moment of the whole OWC-WEC decreases, even though the resonant loads of the connection beam increase greatly.Figure13shows the power extraction efficiency with different chamber breadths.It can be seen that the effective frequency bandwidth increases sharply with the chamber breadth for capturing longer waves.The reason for this is that the inertia of the water column in the chamber increases with chamber width.Thus, long wave energy captured can make a greater contribution to the piston-type oscillation in the chamber.It is also found that the resonant frequency decreases rapidly with chamber breadth.In other words, the proposed OWC with a wider chamber tends to possess better efficiency.

Figure 12 .
Figure 12.Effect of the chamber breadth on the wave loads for b/h = 0.10, b/h = 0.20, b/h = 0.30, and b/h = 0.40: (a) wave force fx, (b) wave moment My, (c) wave force fxbeam, and (d) wave moment Mybeam.Figure 12. Effect of the chamber breadth on the wave loads for b/h = 0.10, b/h = 0.20, b/h = 0.30, and b/h = 0.40: (a) wave force f x , (b) wave moment M y , (c) wave force f xbeam , and (d) wave moment M ybeam .

Figure 13 .
Figure 13.Effect of the chamber breadth on the energy extraction efficiency for b/h = 0.10, b/h = 0.20, b/h = 0.30, and b/h = 0.40.

Figure 13 .
Figure 13.Effect of the chamber breadth on the energy extraction efficiency for b/h = 0.10, b/h = 0.20, b/h = 0.30, and b/h = 0.40.

Figure 15 .
Figure 15.Effect of the chamber wall thickness on the energy conversion efficiency, for d wall /h = 0.01, d wall /h = 0.05, d wall /h = 0.10, and d wall /h = 0.15.