Hydraulic Performance of an Innovative Breakwater for Overtopping Wave Energy Conversion

: The Overtopping BReakwaterfor Energy Conversion (OBREC) is an overtopping wave energy converter, totally embedded in traditional rubble mound breakwaters. The device consists of a reinforced concrete front reservoir designed with the aim of capturing the wave overtopping in order to produce electricity. The energy is extracted through low head turbines, using the difference between the water levels in the reservoir and the sea water level. This paper analyzes the OBREC hydraulic performances based on physical 2D model tests carried out at Aalborg University (DK). The analysis of the results has led to an improvement in the overall knowledge of the device behavior, completing the main observations from the complementary tests campaign carried out in 2012 in the same wave ﬂume. New prediction formula are presented for wave reﬂection, the overtopping rate inside the front reservoir and at the rear side of the structure. Such methods have been used to design the ﬁrst OBREC prototype breakwater in operation since


Introduction
Energy consumption has been one of the most salient ways of measuring progress in society.This is especially true nowadays: fossil fuel is not only a cultural phenomenon; it is an economic necessity for many developing/developed countries.However, new parameters, such as energy efficiency, are beginning to be used to estimate the well-being of individual states.In this contest, the renewable energy share of global consumption represents a world-wide index by which it is possible to assess the technological advancement of a country.
Among the various renewable energy sources, ocean energy has attracted the attention of the business and scientific community from as early as 1973.The main reason is that the resource is so vast [1,2], i.e., the theoretical potential resource of ocean energy is more than sufficient to meet present and projected global electricity demands well into the future.
In the last decade, numerous research projects aimed both at evaluating the potential energy and also designing new types of converters (called Wave Energy Converters (WECs)) were carried out.This interest is motivated by the various advantages that characterize such a source [3]: the high energy density, greater than that of solar and wind; the easy prediction of the wave characteristics through numerical models [4,5]; the reduced energy loss during wave propagation in relative water depth.However, these benefits are offset by the following drawbacks: the high variability of the wave characteristic through time [6]; WECs are exposed to large environmental forces; high production costs compared to other devices as photovoltaic and wind turbines [7].
A solution to significantly decrease the costs would be to develop hybrid devices that can be embedded within coastal or offshore infrastructures [8][9][10].This important new concept for coastal defense structures could be a realistic solution for the WEC systems to become economically competitive with other renewable energy devices, especially considering the fact that they can be integrated within existing breakwaters or by upgrades.This integration has several advantages from an economical, constructional and operational point of view.Construction costs are shared, and the access to construction, operation and maintenance of the devices would become much easier.Among the large number of WEC technologies, only very few devices have been constructed at the prototype scale, and not one is at the commercial stage.
The utilization of wave energy close to the shoreline is attractive thanks to reduced costs concerning construction, access, maintenance and grid connection.On the other hand, the available wave energy is less than at deep water sites, although recent work (e.g., [11][12][13]) demonstrated that there are sites where energy is concentrated due to wave transformation phenomena, such as wave refraction.
OBREC is a nearshore device combining rubble mound breakwaters with WECs.The device is able to extract energy through the wave overtopping phenomenon.Instead of dissipating the incoming wave energy on the breakwater armor layer, OBREC uses a concrete ramp in order to increase the overtopping discharge and a front reservoir designed to capture the wave overtopping in order to convert wave energy into potential energy.Water stored in the reservoir produces energy by flowing through low head hydraulic turbines, using the difference in water level between the reservoir and the main sea water level (Figure 1c).
To estimate the hydraulic and the structural performance of OBREC, 2D physical model tests were carried out at Aalborg University (DK) in 2012 (AAU12) [8,11,23].
The AAU12 tests were aimed at estimating the main differences between a traditional rubble mound breakwater and OBREC.Two different configurations, characterized by different ramp lengths, respectively, 7.5 cm and 12.5 cm at the model scale, were tested.Hydraulic results showed that the integration of the device in a traditional breakwater improves the overall performances [8].As regards the hydraulic performances, the main results are: (1) the device shows a similar or even reduced reflection coefficient with respect to traditional rubble mound breakwater; (2) overtopping at the rear side of the structure is reduced by adopting appropriate precautions, e.g., the realization of a parapet at the crest of the OBREC crown wall; (3) new design methods have been proposed for the estimation of the reflection coefficient, overtopping at the rear side of the structure and overtopping volume in the front reservoir.
However, the influence of some factors were not evaluated during the AAU12 tests.Indeed, the AAU12 tests did not allow for an understanding of how the hydraulic performance (reflection and wave overtopping) could be affected by the following geometrical characteristics: reservoir width, frontal ramp length and frontal ramp shape.In order to evaluate the influence of such parameters in the hydraulic performances of the device, new tests (AAU14) were carried out.The AAU14 tests were carried out at the Hydraulic and Coastal Engineering Laboratory of Aalborg University at a length scale of 1:30 (Froude scaling) compared to typical prototype dimensions.Different geometric configurations were investigated by varying the width of the reservoir, the water level and the profile of the frontal ramp.A few preliminary results on hydraulic performances have already been presented by Iuppa et al. [24].
This paper is organized as follows: Section 2 provides information on the experimental procedure and setup.In Section 3, the hydraulic performances of the device are evaluated: the reflection coefficient of the structure, overtopping at the rear side of the structure and the overtopping volume in the front reservoir.Section 4 is devoted to an overall discussion with some concluding remarks.

Wave Flume
The model tests were carried out at the Hydraulic and Coastal Engineering Laboratory of Aalborg University at a length scale of 1:30 (Froude scaling) compared to the typical prototype dimensions.The wave flume is 25 m long, 1.50 m wide and 1.20 m deep.The flume configuration is shown in Figure 2.
Moving from the paddle (a hydraulic-driven piston mode generator) to the model, the bottom was horizontal for the first 6.5 m, with a 3.5-cm step, a 1:98 slope section with a length of 9 m and, finally, a horizontal section where the model was placed.The flume was divided into two sub-flumes by a guiding wall in order to test two different device configurations.Each part was 0.73 m wide.

Tested Configurations
The tests were performed by varying: (i) the shape of the frontal ramp; (ii) the height of the frontal ramp, R r , with respect to the still water level (SWL); (iii) the width of the frontal reservoir.
Figure 3 shows the cross-section of the analyzed configurations in the AAU14 tests.In the figure: B r is the reservoir width; B s is the emerged sloping plate width; ∆B rs is the horizontal distance between the crown wall and the crest of the ramp; h r is the depth reservoir; R c is the crest free-board of the crown wall; R r is the crest free-board of the front reservoir; ∆Rc is the vertical distance between the crown wall and the crest of the ramp; d w is the height of the sloping plate; d d is the height of the submerged sloping plate; h is the water depth at the toe of the structure.As regards the shape of the frontal ramp, two different configurations were tested: the first (defined as the flat configuration) was characterized by a flat ramp with a slope angle of 34 • ; the second (the curved configuration) has a curvilinear ramp where the slope angle varies linearly between 52 • and 17 • .The slope of the flat ramp was driven by Kofoed [25] where a maximization of the amount of overtopping was observed for such a geometry.The curved configuration was tested in order to verify for the on shore condition the results of Kofoed [26].In this study, performed for the offshore wave condition, such a geometry maximize the overtopping discharge.
The height of the frontal ramp R r with respect to the s.w.l.influences the energy conversion.Indeed, a low values of R r produces high overtopping discharge in the front reservoir, but at the same time results in low hydraulic head on the turbine.The other way around, high values of R r produce small discharge in the reservoir and high hydraulic head.Typical values of R r are in the range of 1-3.5 m.Furthermore, the study of the response of the system to R r variation is relevant in order to evaluate the influence of tide.
The width of the front reservoir ∆B rs can influence the overtopping discharge trough the overall structure.Typical values of ∆B rs are in the range of 5-15 m.
The model geometrical characteristic are shown in Table 1 together with those of the AAU12 tests.A total of 9 cases were analyzed for each configuration.Hereafter, the symbols R1, R2 and R3 are used to indicate the three different values of R r .
In order to reduce the overtopping discharge through the overall structure, a parapet has been placed on the top of the upper crown wall (Figure 4).Indeed, as indicated by Vicinanza et al. [8] the presence of a "nose" causes a strong reduction of the overtopping at the rear side of the model.The parapet had the shape of an isosceles triangle with vertical and horizontal sides of 2 cm in the model scale.Such dimensions were defined on the basis of the authors' experience.The rubble mound material was chosen in order to ensure the stone stability under wave action and to reproduce the main hydraulic behavior of the structure.The equivalent cube side length exceeding by 50% the stones (D n,50 ) for the armor layer was 50 mm.In order to reproduce the turbulent flow inside the filter layer, a D n,50 equal to 20 mm was assumed in such a layer.Finally, the D n,50 was equal to 5 mm for the core, in order to prevent the washing of the core material.

Wave Characteristics
Waves were generated from a hydraulically-driven piston mode generator controlled by the software AwaSys developed by Aalborg University.Simultaneously, active absorption of reflected waves was used in all tests [27].Waves were generated based on the three parameters in the JONSWAP (JOint North Sea Wave Project) spectrum: significant wave height H m0 , peak period (T p ) and the peak enhancement factor γ (γ = 3.3 in all tests).Each test contained at least 1000 waves.
The data obtained from the eight wave gauges (four for the model) were analyzed with the software WaveLab (developed at Aalborg University).The software allowed one to estimate the wave characteristics (e.g., incident wave height, reflected wave height, energy wave period) with the method of Zelt and Skjelbreia [28].A total of 200 tests was carried out.
Table 2 shows the maximum and minimum values of the water level and of the wave characteristics estimated at the models' toe.
Table 2. Maximum and minimum values of the water level and wave characteristics at the models' toe evaluated trough the method of Zelt and Skjelbreia [28].Table 3 reports the range of some dimensionless parameters.L m−1,0 and ξ m−1,0 represent the wavelength and the breaker parameter, respectively, referenced to the spectral incident energy wave period T m−1,0 .

Instruments
The laboratory campaign was carried out in order to gather data on the following characteristics: wave reflection, wave loading, wave overtopping discharge both in the front reservoir and behind the whole structure.
In each sub-flume, the following instruments were installed: 4 resistance wave gauges in order to evaluate incident and reflected wave spectra; 2 boxes to collect the water discharge in the front reservoir and behind the whole structure; 2 depth gauges, protected by hollow cylinders in PVC in order to measure the overtopping discharge; 14 pressure transducers for the estimation of the pressures/forces induced by the waves on the structure.
The distance between the models and the wave gauges was approximately 2.80 m (Figure 2), according to the recommendations of Klopman and van der Meer [29].
In each overtopping accumulation box, a water level gauge was installed in order to measure overtopping discharge and to control evacuation pumps (Figure 5).The overtopping accumulation box was connected to the frontal reservoir by a PVC pipe, while the rear side-overtopping accumulation box received the discharge through a ramp placed on the top of the crown wall.The PVC pipe outflow was set at the same level of the reservoir bottom.Overtopping discharges at the rear side of the models and in the front reservoirs were estimated using the water level gauge measurements.Once the water level inside the box is known, it is possible to estimate the flow rate Q as: where ∆V is the overtopping volume variation in the ∆t, A(h) is the cross-sectional area of the box and h box is the water level.The area A(h) is a function of the water level due to the presence of the pipes and pump used to extract the water.A(h) was estimated by measuring the water depth in the box after a known water volume was introduced in the box.

Results
The overtopping discharge Q in through Section S 1 is the sum of three terms: where Q reservoir is the flow through Section S 2 , Q rear is the flow through Section S 3 and Q over f low is the reflected overflow outgoing from the reservoir (Figure 6).The water effectively collected in the reservoir, generating the Q turbine .In Figure 7, the non-dimensional time average wave overtopping per meter width: is plotted against the relative crest free-board: where q reservoir is the time average wave overtopping per meter width.The point data refer to tests with R r = 0.125 m.Two zones could be clearly identified: in Zone I, q * reservoir increases with the decreasing of R * r ; in Zone II, the q * reservoir decreases with the decreasing of R * r .The trend observed in Zone I is typical of an overtopping process.In this zone, Q rear and Q over f low are null, and Qin = Qreservoir = Qturbine (the overbar represent the time average values).
Zone II identifies a range where the reservoir operates in a saturated condition.In such a condition, a large volume of water is lost (i.e., Q over f low increases).For this reason, q * reservoir decreases with the decreasing of R * r .In this zone, Qin > Qreservoir and Qreservoir = Qturbine .The overtopping rate observed in Zone II is strongly affected by test-specific conditions (e.g., reservoir dimensions, localized and distributed load losses, hole geometry).For this reason, the overtopping rate observed in Zone II will not be taken into account in the analysis.Therefore, the following sections describe the results obtained in Zone I.

Overtopping Discharge in the Front Reservoir
The wave overtopping in the front reservoir can be analyzed as in the case of structures with a single slope.For these types of structure, in the case of no breaking wave, the wave overtopping q can be expressed as by the functional relationship f : where R indicates the crest free-board of the structure and γ indicates a coefficient that takes into account all of the factors that reduce the wave overtopping.In EurOtop Manual [30], the average dimensionless overtopping discharge per meter structure width is estimated as: where: c 1 and c 2 are empirical coefficients; γ f is the influence factor for the permeability and roughness of the slope; γ β is the influence factor for the oblique wave attack; γ b is the influence factor for a berm.The coefficient c 1 is equal to 0.2, while c 2 is equal to 2.6 for the probabilistic method and 2.3 for the deterministic method.Equation ( 6) was extended by Victor and Troch [31] in the case of an impermeable slope.In particular, Victor and Troch [31] have estimated the values c 1 and c 2 as a function of the slope angle of the structure, the crest free-board and the significant wave height.
The authors observed that the effect of the wave period can be neglected due to non-breaking waves on the slope.
Based on the AAU12 tests, Vicinanza et al. [8] proposed a new formula to estimate the average dimensionless overtopping discharge.The formula was developed on the basis of the tests conducted for two configurations of OBREC.In particular, two different lengths of the ramp (d w ) were tested (Section 2.2 provides a description of the geometrical and wave characteristics of the AAU12 tests).The equation proposed by Vicinanza et al. [8] is the following: where the parameter s Rr is named the wave-structure steepness: The application range of Equation ( 7) is 0.64 < d w ∆R c < 1.35 and 0.0123 < s Rr < 0.202 [8]. Figure 8 shows the comparison between the measured non-dimensional average front reservoir overtopping discharge and the prediction methods: (a) Vicinanza et al. [8]; (b) EurOtop Manual [30].The prediction methods' reliability is expressed using the root-mean square error (rmse), defined by the following relationship: where N is the number of the test, q * Obs is the dimensionless overtopping discharge observed and q * Est is the dimensionless overtopping discharge estimated using the prediction method.The Vicinanza et al. [8] formula overestimates the observed values.Indeed, the rmse values are: 8.23 for R1, 7.90 for R2 and 7.71 for R3.The reason is straightforward and addressable by the significant difference in the model setup of the AAU12 and AAU14 tests (Figures 3 and 4; Table 1).Indeed, the AAU12 model was equipped with a toe berm, and the submerged ramp was absent (i.e., new configurations are out of the range of application of Equation ( 7)).
The formula of EurOtop Manual [30] fits quite well the observed values using γ f = 1.The rmse values are: 0.47 for R1, 0.22 for R2 and 0.01 for R3.However, a more detailed analysis shows that the decrease of d d leads to a difference between the observed values and the values estimated using Equation ( 6).This behavior is due to the effect of the amour roughness.It was observed that the shorter the height of the submerged sloping plate (d d ) with respect to the wavelength, the larger the effect of the amour roughness.This aspect can be clarified observing Figure 9, which shows the ratio between the dimensionless measurement overtopping discharge and the one estimated by Equation ( 7) More in detail, the effect of the amour roughness is negligible for d d /L m−1,0 > 0.05, and Equation ( 7) with γ f = 1 predicts the observed overtopping discharge with enough accuracy (the spread is ±20%).For d d /L m−1,0 < 0.05, the effect of the amour roughness is very relevant, and Equation (7) with γ f = 1 overestimates the overtopping discharge.
Bruce et al. [32] estimated γ f for various types of homogeneous armor units.For a composite structure, as OBREC, a constant value of γ f cannot be defined, but such a parameter has been evaluated according to the incident wave characteristics and to the length of submerged ramp.In order to evaluate the effect of the amour roughness for the analyzed configurations, γ f , which minimize the differences between the values of q reservoir observed and those estimated with the EurOtop Manual [30] relationship (Equation ( 6)), was estimated (see Figure 10), and a relationship between γ f and the relative ramp depth d d L m−1,0 was determined: where s 1 = 7.47 and s 2 = 0.42.The lower limit γ f = 0.7 was detected by analyzing the AAU12 tests, which were performed with d d close to 0 m.Indeed, for such tests, it was observed that Equation ( 6) can be applied adopting values of 0.7 for γ f [8].
Figure 10a shows also the comparison between the value of γ f estimated on the basis AAU14 tests and predicted values by Equation (10). Figure 10b shows the comparison between the q reservoir observed and q reservoir estimated with Equation (6) using the γ f computed from Equation (10).10); (b) comparison between the q reservoir observed and q reservoir estimated with Equation ( 6) by using the γ f computed from Equation 10.In (b) the red line indicates the perfect matching between the two quantities.
Comparisons of q * reservoir between flat and curved configurations are presented in Figure 11.As can be seen from Table 1, the ramp crest free-board (R r ) of the flat configuration is 4 mm greater than the curved configuration.Hence, in order to make the comparisons more reliable, non-dimensional overtopping rates in Figure 11 are divided by the relative crest free-board.The comparison shows that the curved configuration gives approximately 22% less overtopping than the flat configuration.Overtopping into the front reservoir: q * reservoir /R * r observed in the presence of the flat configuration versus q * reservoir /R * r observed in the presence of the curved configuration.

Reflection
Several previous studies have been performed to analyze the reflection caused by coastal structures.Some of these studies have enabled the development of new prediction methods to estimate the reflection coefficient K r , i.e., the ratio between reflected (H m0,r ) and incident wave heights (H m0,i ).
For rubble mound breakwater, the reflection coefficient can be expressed by the following functional relationship: where α is the structure slope angle and L m−1,0 is the deep water wavelength estimated with reference to T m−1,0 at the structure toe.
Recently, Zanuttigh and van der Meer [33] have developed a new estimation formula for various types of structures.This formula allows one to estimate K r once known: incident wave characteristics (H m0 and T m−1,0 ) at the structure toe, the slope of the structure (α) and roughness coefficient γ f : where ξ m−1,0 is the breaking parameter and the coefficients a and b are defined by the following relationship: Vicinanza et al. [8] suggested the adoption of Equation ( 12) for the estimation of K r .This method, in fact, appears opportunely conservative.In more detail, the analysis showed a good agreement between the predicted and the observed K r adopting a γ f value equal to 0.55.Moreover, for the configuration in the AAU12 tests with greater q * reservoir , a reduction of K r was recognized.This behavior was explained by the fact that the device captures the incoming wave energy.Such an aspect was also observed for breakwaters by Zanuttigh and van der Meer [33], as the lower the crest, the greater the overtopping and the lower the reflection.In order to take into account such behavior, the authors introduce a reduction factor of the reflection coefficient: where R * is the relative crest free-board of the structure.Equation ( 14) was estimated for rock permeable slopes, and it can be applied to the following range: R * ≥ −1; H m0 /D n50 ≥ 1; and s m−1,0 ≥ 0.01. Figure 12 shows the comparison between the reflection coefficients (K r,obs ) measured in the AAU14 tests and those estimated with Equations ( 12) and ( 13) (K r,est ).A good correspondence can be found.In particular, the value of K r,est is estimated using three different values of γ f with the variation of the water level: 0.55 for R1; 0.8 for R2; and 0.9 for R3.
As can be seen from Figure 12, an increasing of the submerged ramp length (d d ) causes an increase of reflection.This aspect is consistent with the observation of the the higher overtopping discharges due to the decrease of the armor roughness observed in the previous section.
In order to take into account also the effect of the submerged ramp length, a new prediction method based on Equation ( 12) has been derived.Using the values of γ f estimated from Equation (10), the method of Zanuttigh and van der Meer [33] overestimates the experimental data.Such behavior is caused by the structural difference between the models used by Zanuttigh and van der Meer [33] and the OBREC.Therefore, to evaluate K r , a corrective coefficient was introduced: with γ f estimated by Equation ( 10) and c γ f evaluated through the following relationship: where s 1,r = 2.64 and s 2,r = 0.28, and X Kr is defined as: The parameters s 1,r and s 2,r were evaluated by the fitting process using the least squares method.Figure 13a shows the comparison between c γ f detected by experiments and the calculated values with Equation 16. Figure 13b shows the comparison between K r observed and K r estimated with Equations ( 12), ( 13) and (15).16; (b) comparison between K r observed and K r estimated with Equations ( 12), ( 13) and (15).
Figure 14 shows the values of K r as a function of the breaking parameter (ξ m−1,0 ) classified depending on R r and on the ramp shape (flat and curved).The comparison shows that for large values of R r (cases R1 and R2), the difference between the flat configuration and curved configuration is negligible.Such a difference becomes relevant for the case with small values of R r , that is R3.This behavior is because: the reflection coefficient tends to increase with the slope.In fact, near the ramp crest, the curved configuration has a slope smaller than the flat configuration, causing an attenuation in reflection.

Wave Overtopping at the Rear Side of the Structure
Coastal engineers design breakwaters generally with the main aim to limit the overtopping discharge at the rear side of the crown wall.
OBREC can be considered as a slope with a crown wall provided with a parapet.For this type of structure, there were no specific studies before the AAU12 tests.In order to address this shortcoming, Vicinanza et al. [8] developed a new prediction formula: where ∆R c is the difference between R c and R r .The range of application of Equation ( 18) is: 0.014 < ∆R c L m−1,0 < 0.038; 0.035 < s m−1,0 < 0.058; 1.24 < R c H m0 < 1.38 [8].Van Doorslaer et al. [34] analyzed several configurations to reduce the wave overtopping over the smooth dike slope.The authors extended the prediction method of EurOtop Manual [30] (Equation ( 6)) for structures very similar to OBREC.One of these structures is composed of a promenade and a storm wall, and another one is composed of a promenade, a storm wall and a parapet.The authors have modified the EurOtop Manual [30] formula as: where γ VD is defined as: for a smooth breakwater with a promenade and a storm wall, and as: for a smooth breakwater with a promenade, a storm wall and a parapet, where the coefficients γ v , γ par , γ prom take into account respectively the reduction effect of the storm wall, the reduction effect of the parapet and the reduction effect of the promenade.Methods for reduction factors' estimation are shown in Van Doorslaer et al. [34].
Figure 15 shows the comparison between the observed values and the Equation ( 18).The overtopping discharge for configuration ∆B rs = 0.30 m is estimated fairly well by Equation ( 18), while for configurations ∆B rs = 0.10 m and ∆B rs = 0.20 m, overtopping is underestimated.Such behavior is because of the fact that Equation (18) does not take into account the effects of the reservoir width.AAU14 measured discharges versus computed values by Equation (19) for two different configurations tested by Van Doorslaer et al. [34] are reported in Figure 16a,b.For the structure without a parapet, Equation (19) interprets quite well the experimental data for ∆B rs = 0.20 m.However, for ∆B rs = 0.10 m and ∆B rs = 0.30 m, the average overtopping discharges are respectively underestimated and overestimated.For the structure with a parapet (Figure 16a,b), the approach of Van Doorslaer et al. [34] tends to underestimate the measured values.The main reason for these discrepancies can be explained by the absence of the reservoir.Table 4 shows the performances of the two prediction methods.Based on the AAU14 tests, a new prediction method was developed.The average wave overtopping can be expressed by the following functional relationship: This expression can be arranged in terms of dimensionless parameters: On the basis of the AAU14 tests with the flat ramp, the overtopping discharge at the rear side of the breakwater could be evaluated by the following relation: where a rear = 0.0139 and b rear = −7.17,and the parameter X rear is defined by the following equation: The parameters a rear and b rear were evaluated by the fitting process using the least squares method.
Figure 17 shows the comparison between the measured overtopping discharges and calculated from Equation (24).Experimental data are excellently interpreted by Equation ( 24) with a good correlation coefficient (R 2 > 0.96) and an rmse = 0.69.
In order to make a comparison and to establish the reliability of Equation ( 24), AAU12 data are also reported in Figure 17.The new formula can also be applied to the curved configuration by adopting an appropriate correction factor.Indeed, the curved shape of the ramp causes less wave overtopping discharge than the flat configuration.A mean difference of approximately 20% was estimated.Therefore, a new correction factor γ q rear was introduced.Such a factor equals one for the flat configuration and 0.83 for the curved configuration.Equation ( 24) can be rewritten as: The comparison between the observed values (flat and curved configurations) and those estimated using Equation ( 26) is shown in Figure 18.

Conclusions
This paper presents a study of the innovative Overtopping BReakwater for wave Energy Conversion (OBREC).Moving on from the complementary tests by Vicinanza et al. [8], 2D model tests were carried out at Aalborg University.The main purpose of this new experimental campaign was to extend the knowledge on the device behavior, with particular interest in the influence of the shape and draft of the front ramp, as well as the reservoir width.This parametric study aims to bring completion to the previous one in which different ramp crest elevations were tested.
Two different device configurations are analyzed: the first is characterized by a constant slope ramp; the second presents a curved ramp profile.For both configurations, tests with three different reservoir widths and three water levels were considered.
Wave overtopping in the front reservoir could be predicted by the method of EurOtop Manual [30] with a relatively accurate estimation.However, due to the dual nature of OBREC, a high level of accuracy in overtopping discharge effectively usable for energy production is required.Thus, a new method to estimate the roughness factor as a function of the submerged ramp length and the wavelength was developed.
As regards the reflection coefficient, the two tested configurations show the same behavior.However, the reflection coefficients measured in the previous test campaign were lower than the present ones.Such a behavior is mainly caused by the different extension of the submerged ramp.The prediction method of Zanuttigh and van der Meer [33] can be used to estimate wave reflection, although it is necessary to apply a correction factor as a function of the wave characteristics and some geometrical characteristics of the ramp.
A reduction of approximately 20% in overtopping rates in the front reservoir and at the rear side of the structure was observed for the curved configuration.This could significantly affect the potential energy production of the system, but at the same time, a higher safety level at the rear side of the crown wall can be ensured in comparison to the flat configuration.
A new prediction method to estimate overtopping discharge at the rear side of the structure was proposed, which was proven to perform remarkably well with respect to the present dataset.Based on the approach proposed by Vicinanza et al. [8], the new formula allows one to take into account the effects of the reservoir width.
The analysis described in this paper represents a further step forward in the knowledge of the OBREC device.Future studies will mainly focus on the evaluation of the energy production performance and how it can be affected by the geometrical characteristics of the device (width and depth of the reservoir, height of sloping plate, etc.) and by the turbine characteristics.

Acknowledgments:
The work discussed here is a part of the National Operational Programme for the "Research and Competitiveness" 2007-2013 (NOP for R&C) founded Project PON04a3_00303 titled "DIMEMO-DIgaMarittima per l'Energia del Moto Ondoso" (Maritime Breakwater for Wave Energy Conversion), Project PON04a3_00303.The work also was partially supported by RITMARE (acronym of la Ricerca ITaliana per il MARE ) Flagship Project and by the project PON02_000153_2939551 "Development of innovative technologies for energy saving and environmental sustainability of shipyards and harbor areas" (SEAPORT, acronym of Sviluppo di tecnologie innovative per la Sostenibilità Energetica ed Ambientale di cantieri navali ed aree PORTuali) and by HYDRALAB PLUS project (Proposal Number 64110).The authors gratefully acknowledge the Italian Ministry of Education, University and Research for supporting this innovative research and the Second University of Naples and the University of Catania for encouraging the mobility of researchers.Important input from Thomas Lykke Andersen in the study design and data analyses is gratefully acknowledged.The authors would also like to thank the Aalborg University technicians, in particular Niels Drustrup, for their cheerfulness, tolerance and considerable assistance in the design evolution, model construction and testing.Assistance in the testing and supporting physical modeling from Vincenzo Ferrante (Second University of Naples) is also gratefully acknowledged.
Author Contributions: Claudio Iuppa contributed to the work described in this paper by collaborating in carrying out the laboratory campaign and analyzing the acquired data.Pasquale Contestabile contributed to the work by designing and collaborating in carrying out the laboratory campaign and analyzing the acquired data.Luca Cavallaro contributed to the work by analyzing the acquired data.Enrico Foti contributed with knowledgeable discussion and suggestion.Diego Vicinanza gave the idea, supervised the data analyses and gave the final approval.All of the co-authors participated in writing the paper.

Conflicts of Interest:
The authors declare no conflict of interest.spectral moment of order 0 m −1 (m 2 • s) spectral moment of order −1 q * rear (-) non-dimensional overtopping discharge towards the rear of the traditional rubble mound breakwater crown wall or towards the rear OBREC crown wall q * reservoir (-) non-dimensional overtopping discharge into the reservoir q rear (m 3 • m −1 • s −1 ) average overtopping discharge towards the rear of the traditional rubble mound breakwater crown wall or towards the rear of the OBREC crown wall q reservoir [m 3 • m −1 crest free-board of crown wall, i.e., the vertical distance between the crest of the vertical walland the still water level R r [m] crest free-board of front reservoir, i.e., the vertical distance between the crest of the sloping plate and the still water level rmse [-] root mean square error s m−1,0 = 2πH m0 vertical distance between the crown wall and the crest of the ramp

Figure 2 .
Figure 2. Plant and cross-section of the wave flume.OBREC, Overtopping BReakwater for Energy Conversion.

Figure 5 .
Figure 5. Wave-by-wave system for flow discharge measurement.

Figure 9 .
Figure 9.Comparison between the dimensionless measurement overtopping discharge and the one estimated by Equation (7) using γ f = 1.

Figure 10 .
Figure 10.Overtopping discharge in the front reservoir: (a) comparison between γ f estimated on the basis of the AAU14 tests and predicted values by Equation (10); (b) comparison between the q reservoir observed and q reservoir estimated with Equation (6) by using the γ f computed from Equation10.In (b) the red line indicates the perfect matching between the two quantities.

Figure 11 .
Figure 11.Overtopping into the front reservoir: q * reservoir /R * r observed in the presence of the flat configuration versus q * reservoir /R * r observed in the presence of the curved configuration.

Figure 12 .
Figure 12.Comparison between the measured K r in the AUU14 test and those estimated by Equation (12).

Figure 13 .
Figure 13.Prediction formula: (a) comparison between c γ f detected by experiments and the calculated values with Equation 16; (b) comparison between K r observed and K r estimated with Equations (12), (13) and(15).

Figure 14 .
Figure 14.Reflection coefficient: comparison between the flat and the curved configurations.

Figure 16 .
Figure 16.Comparison between the observed values and Equation (19) (γ v = 0.78 ; γ par = 0.72): (a) the case of the structure without a parapet; (b) the case of the structure with a parapet.

Figure 17 .
Figure 17.Comparison between Equation (24) and the data observed in the AAU14 tests with a flat ramp and the AAU12 tests.

Figure 18 .
Figure 18.Overtopping at the rear of the structure: comparison between Equation (26) and the observed values during the AAU14 tests (flat and curved configurations).
B r (m) reservoir width B s (m) emerged sloping plate width d d (m) height of the submerged sloping plate D n50 (m) equivalent cube side length exceed by 50% of the stones d w (m) height of sloping plate g (m • s −2 ) gravity acceleration h (m) depth at the toe of the structure h box (m) depth in the accumulation box H m0,r (m) reflected significant wave height at the toe of the structure H m0 (m) incident significant wave height at the toe of the structure h r (m) depth in the front reservoir K r (-) H m0,r H m0 reflection coefficient L m−1,0 (m) deep water wavelength referenced to T m−1,0 m 0 (m 2 ) wave steepness at the toe of the structure s Rr[-]   non-dimensional wave-structure steepnessT m−1,0 = m −1 m 0 [s] spectral incidentenergy wave period at the toe of the structure T p [s] incident peak wave period α [ • ] slope angle of the structure γ [-] peak-enhancement factor γ β [-] reduction factor for oblique wave attack γ b [-] reduction factor for berm γ f [-] reduction factor for slope roughness γ v [-] reduction factor for the storm wall γ par [-] reduction factor for the parapet γ prom [-] reduction factor for the promenade ρ [kg • m −3 ] to T m−1,0 ∆B rs = B r − B s [m] horizontal distance between the crown wall and the crest of the ramp ∆R c = R c − R r [m]

Table 1 .
Geometrical characteristics of the different configuration tested.The values refer to the Aalborg University 2014 (AAU14) and AAU12 tests. )
• s −1 ] average overtopping discharge into the reservoir R [m] crest free-board of the structure R * c = R c H m0 [-] relative crest free-board of crown wall R * r = R r H m0 [-]relative crest free-board of front reservoir R c[m]