Duct Attachment on Improving Breaking Wave Zone Energy Extractor Device Performance

A challenging wave energy converter design that utilized the denser energy part of the nearshore breaking wave zone to generate electricity was introduced in 2016 by Shintake. The Okinawa Institute of Science and Technology Graduate University’s project aims to take advantage of breaking wave energy to harness electricity. The 2016 version of the device consisted only of a bare turbine and power generator. Early exploration of the design recorded short periods and high impact wave pressures were experienced by the structure, with the turbine unable to harvest energy effectively. Additional structure to not only reduce incoming impact pressure but also increase the duration of water flow through the turbine was needed. These are the main reasons behind incorporating the duct attachment into the design. This paper show that the duct is capable of halving the impact pressure experienced by the turbine and can increase the energy exposure by up to 1.6 times the bare turbine configuration. Furthermore, it is also said that wave angle (β) = 40° is the critical angle, although the duct still increases wave energy exposure to the power take-off up to β = 60°.


Introduction
In 2020, the Japanese government pledged for Japan to be carbon neutral by 2050. This ambitious goal requires huge efforts from research and scientific communities to find more sustainable and renewable sources of electricity. It is estimated that Japan has up to 260 TWh/year of ocean energy within its reachable offshore (30 km off the coastline) [1]. That is roughly 40% of the country's total electricity consumption per year and thus could provide a major contribution to this carbon-neutral goal. Takahashi and Adachi [2] estimate that around 7 kW/m of power hits Japan's coastline every year. Sasaki [3] further analyzed long-term wave data across Japan at 25 different locations from 1980 to 2009-and found that there was an increasing trend of wave energy, especially around Okinawa Prefecture, by about 0.4 kW/m over the next 30 years. This data increases the appeal of using wave energy as a source of electricity.
There have been many attempts to harvest energy delivered by waves over the years with main focuses on oscillating water column (OWC) type and point absorber type wave energy converters (WEC), which deployed in the shallow water part of the ocean. In principle an OWC device takes advantage of the hydrodynamics changes induced by passing waves to oscillate the water column inside the OWC structure. This movement then forces the air to inhale into and exhale from the chamber. The power is generated by placing a power take-off (PTO) mechanism on the airways. In 2001, a groundbreaking OWC type WEC called Land Installed Marine Power Energy Transmitter (LIMPET) was installed in the Isle of Isla, Scotland, marking the first WEC device to be connected to the main power grid in the UK [1][2][3]. Together with the development of LIMPET, another near shore stand-alone OWC project was also started at Pico in the Azores, Portugal. The be mitigated by utilizing simpler design, smaller size, and closer installation location. Based on these considerations, the Okinawa Institute of Science and Technology Graduate University (OIST) WEC design was design to be placed in the shallower part of the ocean where the waves sometime break. Because of this, it is better for the device not to be designed for very harsh sea conditions, but rather moderate sea conditions with typical significant wave heights of 0.75 m to 1.0 m.
When designing an effective wave energy extractor, one must consider the wave height, wave period, and direction of the wave [17], as an unfavorable wave angle could reduce the device's performance by as much as 75% [18]. This condition is important as [19] argues that by 2042 the long wave, where the peak wave period T p > 8 s, will become more frequent, so much so that the wave energy converter's design made under current wave conditions may become less efficient in the future. This is critical as other coastal structures, such as seawalls and breakwaters, are designed to last centuries.
Unlike deep-water waves, once a wave travels to a shallower area, the wave diffracts due to the interaction between the wave and the seabed. This action forces the direction of the wave to be parallel with the shoreline. Furthermore, around the breaking wave zone area, the wavelength shortens, and the wave height increases due to shoaling, before it ultimately breaks. Following the linear wave theory, in shallow water, wave phase velocity also becomes independent of the wave period (T) and both the group velocity and phase velocity become similar. Thus, wave characteristics become much more uniform in shallow water surf zones, as they are less affected by deep-water sea characteristics. Due to this consideration, one can speculate that a device designed for a narrower band of wave characteristics when placed in this area may enable the device to work more efficiently during operational hours. The shape of the wave, furthermore, is slightly elongated as the water around the still water level travels faster than water near the seafloor. This condition means that the horizontal water flow is faster compared to the vertical motion of the water; therefore, a widely available horizontal axis-type turbine would be suitable to capture wave energy. This is the idea behind the breaking wave zone Wave Energy Converter (WEC) developed by the Okinawa Institute of Science and Technology Graduate University (OIST) [20]. Furthermore, since it is less likely to be affected by deep-sea waves, the same design may work well in many locations, as long as there are enough waves (energy) to justify the production and installation cost.
To take advantage of the breaking wave vortex, the diameter of the OIST WEC's turbine is limited to the wave height. This relatively small dimension means that the OIST WEC is relatively cheaper to fabricate, easier to install, and can be located closer to the shore, thus reducing the cost further in both installation and maintenance, and in line with the lesson learned from the failures of WEC projects in the past. The power generator itself consists of a five-blade turbine which is directly connected to a permanent magnet generator. It is expected that the turbine will operate in harsh environments, so this configuration is expected to be able to withstand the breaking wave impact. To minimize the possibility of sea water entering the generator, the inside of the generator is filled with pressurized oil. More details of the turbine design and considerations can be seen in [21].
To select a suitable design for the intended testing location, an onsite wave measurement was done at the Port of Seragaki in October 2015. The measurement was recorded using five Teledyne RD Instrument Acoustic Doppler Current Profilers (ADCP) called workhorse sentinels. The recording frequency was 2 Hz for one month. It was found that the Port of Seragaki, Okinawa, has a horizontal velocity between 1.8 m/s to 3.5 m/s at a depth of 0.4 m. Furthermore, the significant wave height (H s ) of the location was also measured to be 1.7 m. For more details on this research, please see [22].
Based on the sea characteristics, the selected turbine was 700 mm in diameter and rated to generate 12 kW at peak power. The turbine characteristics are summarized in [23] and a preliminary field test on a 1:2 scale physical model was done in Maldives. The Maldives location was chosen due to its relatively uniform wave characteristics throughout the year (T p = 9.8 s; 1.0 m < H s < 2.0 m for 68% of the year) [24]. Two half-scale physical models of the turbine tested here were located at the same distance from the shore, about 10 m distance between each model. Each model, however, was mounted at a different height. One was about 1 m from the seabed and roughly located just above the still water level during low tide. The other was mounted about 1.5 m from the seabed, located just above the still water level during high tide. Due to the impulsive nature of breaking waves, the installation work of the WECs brought a significant challenge. First, a 75 cc two-cycle gasoline engine core drill was utilized to drill four 110 mm diameter holes on the seabed, at the designated location. Each hole had the depth of about 400 mm. Then the mount sockets were inserted into the holes with the gap between the hole and the socket filled with 2-3 mm diameter ceramic balls and epoxy cement. Once the cement was dry and the socket was strongly attached to the ocean floor, the half scale turbine generator was towed to the location by boat and fastened to the socket. To store the generated electricity momentarily, a set of 24 electronic double layer capacitors was used and connected to a set of LED lights to lighten the monitoring room. The experimental results showed that the turbine installed around the low tide still water level produced, on average, more power than the one located in a higher position. The peak power measurement, however, showed that the higher location had higher peaks, albeit with longer non-productive periods. The details of the experiments, along with the novel monitoring system of the half-scale physical model testing in Maldives, can be read in [25,26].
Full-scale turbine testing of the turbine generator was done in November 2018. Two full-scale 0.6 m diameter OIST WEC turbines were installed at about 10 m landward from the half-scale turbines with the same 10 m distance between them, using the same installation method of the half-scale experimental campaign. This location was selected based on video observation of the breaking wave location. The center of the turbines was located 0.6 m above the seafloor. This allowed the bottom of the turbines to reach the low tide still water level. For reference, the difference between low tide and high tide in Maldives during this experiment was around 0.8 m, so it was expected that the turbines would be fully submerged during high tide. It is also important to note that half-scale data revealed that energy was still generated while submerged. The full-scale generator showed a lower turbine RPM when compared to the half-scale WEC's model, while the power produced in total was larger. Figure 1 shows the relationship between the turbine rotation speed in revolutions per minute (RPM), against the voltage output of the generator. As can be seen from the figure, the full-scale (blue triangle) line has a much steeper gradient when compared to the half-scale (orange circle) result. For more detailed analysis, please see [27]. Based on the Maldives experiment, it was found that current devices are still unable to harness the breaking wave energy optimally [28]. The preliminary observation showed that most of the time, the very fast and turbulent water and air mixture of a breaking wave hit the turbine when it was not rotating (still) and passed before it could absorb all of its energy. Because of this, a duct was incorporated into the turbine design to act as a concentrator of the water flow toward the turbine and slow down the incoming water flow. Furthermore, these breaking waveloads on a still turbine led to a high impact pressure experienced by the turbine's blade. Bagnold [29] speculates that high impact pressure occurs due to air trapped and compressed inside a body of water. Reference [30] further shows that a small bubble may exacerbate this impact pressure peak. Due to this, removing the bubble may be beneficial in improving the survivability of the turbine, so the duct is expected to also suppress the incoming near-breaking wave and reduce the aeration of the breaking wave.
There are two parameters of the duct design that need to be taken into consideration: The inlet dimension and the duct length. In order to make sure that most, if not all, of the water will enter the duct, it needs to be at least twice the breaking wave height. The breaking wave height can be calculated using as suggested by [31], where L pi denotes the peak period wavelength in the water depth h, and can be calculated using where g denotes the gravity acceleration and T p denotes the peak period of the significant wave. Furthermore, C * is the reflection constant depending on the overall reflection coefficient of the beach C r and can be calculated using Unfortunately, getting the necessary parameters in the real ocean was difficult. Because of that, early exploration on the design of the duct was done numerically by recreating Maldives' Kandomaa Island bathymetry in a numerical wave flume. The simulations were calculated by Free-Surface Flow Solver for Gravity-Driven Hydraulic Events, or SURGE, which is a computational fluid dynamics solver for free-surface flow, modelled using the volume-of-fluid method (VOF) based on the finite difference method. The results show a breaking wave height (H br ) equal to 1.0 m when the incident wave height at water depth (h) equal to 5 m is H inc = 0.75 m, which are the typical wave condition at the designated location. More details of the experiment can be seen in [32]. Based on this result, the inlet dimension of the duct is selected to be two times the H br , or a 2.0 m by 2.0 m rectangle, where the still water level is located at the center of the duct. This should guarantee that most, if not all, of the water of the incoming wave enters. A contraction factor of 9 is selected to increase the flow velocity and increases the performance of the device, leading to a 0.67 m by 0.67 m rectangular outlet dimension. To decide the length of the duct, one needs to consider the distance a breaking wave travels from the moment it reached the breaking wave height (H br ) until the wave is broken. One can imagine that this distance is similar to the diameter of the circular air pockets, which form when a wave is breaking, and can be approximated by as proposed by [33], where D denotes the characteristic linear dimension of the circular air pockets formed by breaking wave with breaking wave height H br . Assuming H br equals 1, as predicted by the numerical simulation, the minimum duct length should be D = 0.26 m.
To be sure that the duct will be able to capture the breaking wave, the early duct design was chosen to be about 20 times that number, to be 4.0. In addition to the breaking wave height analysis, two design configurations were tested in the numerical wave flume, using the same numerical method, as reported in [34]. The experiments compared the total pressure (P tot ) and the effective pressure (P e f f ) at the duct outlet, but only if the breaking wave occurred in the middle of the duct. Effective pressure is the total pressure with the gravitational potential power removed because the gravitational potential power cannot be captured by the electric generator.  Based on the numerical results shown in [34], the first configuration was selected for a preliminary half scale field testing which was done in Maldives (see Figure 3). Early prototypes of the duct design were made from aluminum plates with a rectangular crosssection as shown in Figure 3a. The aerial view of the initial testing using the painted prototype (corrosion prevention) in Maldives can be seen in Figure 3b, and most of the sea water enters the duct as expected. The half-scale turbine was then mounted at the outlet of the tapered duct to measure the power generated and to compare it to the bare turbine configuration. The measured voltage generated by the half-scale turbine and duct configuration showed an increase by a factor of 2.5 when compared to the bare turbine as reported in [26]. Theoretically, this may lead to six times more power being harvested just from adding the duct. Unfortunately, due to the very high breaking wave impact, the corner wall of the duct failed with the water leaked. Furthermore, the turbine mount attached at the outlet bent. Because of this, both designs shown in Figure 2 were combined for the final design. While the outlet pressure is expected to be lower, it is expected that the device will be able to hold the incoming breaking wave pressure. Figure 4a,b show the duct design based on the lesson learned from preliminary testing in Maldives with dimensions shown in mm; it is the one being considered in current study. The length is shortened to 2.4 m to simplify the fabrication, thus further reducing the cost associated with the installation. This length is about 10 times the minimum length calculated using Equation (4). Furthermore, data also show that a previous design's contraction factor may lead to a significant amount of water being reflected seaward, so the contraction factor is lowered to 2.8. Furthermore, a rectangular cross-section, although easier to fabricate, may lead to a pressure concentration at the corner, so a circular cross-section is considered for the current study. In summary, the current study aims to characterize the influence of the most recent duct design to the overall performance of the OIST WEC device, by means of Figure 4a quantifying the increase in energy exposure to the turbine by adding the duct into the WEC design and Figure 4b observing the impact of the directionality of the incoming wave to the performance of the ducted WEC design when compared to a bare turbine design.

Methodology and Experimental Apparatus
To further observe the effects of the current duct design to the overall performance of the device, small-scale physical model experiments were done at the OIST's wave flume facility. The wave flume is a 1 m × 0.5 m × 0.2 m and is made of a 10 mm thick perplex board to enable multi-angle observations. A paddle-type wavemaker is installed at one end and wave-absorbing sponges are placed at the other end of the tank. The paddle is connected to a step-motor which is operated by an Arduino ® board using a computer. The wave generator is tuned to a single monochromatic wave characteristic with a nominal wave height (H) equal to 0.07 m and an incident wave period of 0.57 s. The wave height is chosen based on real wave conditions in Kandooma Island, Maldives, and scaled down following Froude similitude law. For the wave period, however, the flume was unable to produce the scaled down version of Maldives' significant wave period. Nevertheless, because this study compared the pressure measurement with and without the duct, it is safe to assume that by keeping similar wave conditions for both the bare turbine configuration and the ducted turbine configuration, the results should be comparable with the one that used the appropriate scaled down version of the wave period of the real sea in Maldives. The wave period is selected as 0.57 s, which is the largest wave period the wave generator can produce at H = 0.07 m.
The simple schematics of the flume with the beach installed can be seen in (a) with the dimensions shown in mm. Here h denotes the water depth at various bare turbine and ducted turbine configurations and the value of h can be seen in Table 1. The beach bathymetry installed is loosely based on Kandooma Island beach bathymetry. The breaking wave's profile at the breaking location generated by the monochromatic wave can be seen in Figure 5b with the profile highlighted with a dashed red line. The wave height generated during the experiment was measured using a high-speed camera located about 230 mm from the duct model, as can be seen in Figure 5c. The camera position was maintained to be perpendicular to the inlet of the duct model in order to minimize error due to parallax. The flume and camera set up for the experiment can be viewed in Figure 5d.
The 1:20 scale of the duct design can be seen in Figure 6a, with the dimensions shown in mm. The small-scale physical model made from acrylic tube is then mounted onto the extruded aluminum holder and fixed into the position using C-clamp, see Figure 6b. This small-scale physical model design consists of a straight tube, followed by a nozzle-like chamber following the one illustrated in Figure 4. The small-scale physical model of the duct is then placed at various positions relative to the location of the breaking waves.    To measure the wave pressure as experienced by the turbine, a P-8300-501G-10 series made with a Nidec-Copal pressure transducer (PT) was placed at the outlet of the duct. The pressure transducer can measure up to 49 kPa with a response time of about 2 ms. The sensor has a diameter of 20 mm which is suitable to simulate the swept area of the scaled-down version of the turbine (about 29 mm). The pressure sensor is placed by itself to simulate a bare turbine configuration and placed at a 10 mm distance landward from the outlet of the duct to allow water to flow past the PT with minimum water reflection. This is also the scaled down distance between the duct and the turbine at full scale.
To get an accurate wave height measurement, a 240-fps high-speed camera was utilized to record a slow-motion image of the wave during experiment. A previously taken measurement grid image placed on the wall of the flume from the same camera angle and location is then superimposed to the recorded video to measure the incident wave height. Because the pressure measured in the current study is only focused on the landward wave pressure, the reflected wave is omitted in the current study and the reflection analysis was, consequently, not done.
The second part of the experiment is to explore the influence of wave angle on the performance of the ducted OIST WEC. Wave angle is defined as the angle between the direction of the incident wave to the normal direction of the WEC device. Naturally, a significant wave angle will reduce the wave energy exposure to the turbine. Contrary to a bare turbine, a ducted turbine should be more sensitive to a big wave angle, as the angle between the incoming wave and the inlet of the duct could potentially prevent the water from entering the duct and flow through the power take-off. Due to this consideration, it is important to extend the characterization of the duct to include the directionality characteristics of the device.
Due to the wave flume and wave generator limitation, generating an angled wave is not possible. So, to do the directional exploration, 3D printed holders were utilized allowing a multiple angle movement to the physical model at the same location. Figure 7a illustrates the holder design using Solidworks ® with each gear representing 10 • of rotation. The printed product made from tough polylactic acid (PLA) is shown in Figure 7b. This rotating mount is placed on the extruded aluminum holder of the duct as shown in Figure 6. To record the pressure measured by the pressure transducer, a Graphtech ® midi logger HV GL2000 series is deployed to record the pressure measured by the pressure sensor. The data recording frequency is 1000 Hz for a duration of about 15 s to 20 s with at least 10 clear breaking waves recorded. Furthermore, the location of the duct and the pressure transducer relative to the breaking wave position (x/L) and the breaking wave height (y/H br ), the corresponding water depth h, and the incident wave angle β used for the current study are shown in Table 1. For the bare turbine case (ND), the x/L and y/H br location indicates the location of the pressure sensor, while for the ducted case (DT) the x/L and y/H br location indicates the inlet of the duct with the pressure sensor located 10 mm landward from the duct's outlet. Here, the x/L = 0 value indicates the location of the breaking wave, and the y/H br = 0 indicates the still water level.

Landward Pressure at Various Horizontal and Vertical Positions at β = 0 •
The first case considered here is the wave angle (β) = 0 • (wave coming from normal direction of the device) case. For this, the measurements were done at multiple horizontal positions relative to the shallow water wavelength and at multiple vertical elevations, relative to the wave height, with zero indicating the still water level. Measurements were done twice, once with only the pressure sensor to measure the pressure experienced by a bare turbine, and then another one with the duct attached in front of the pressure sensor. It is important to note that for the bare pressure sensor case, the location is stated at the PT location, while for the case where the duct is deployed, the location stated is at the inlet of the duct, not the PT location. Figure 8 shows the comparison of pressure measured at the still water level of the bare pressure transducer (blue circle) and the ducted case (orange triangle). The y-axis shows the impact pressure measured, normalized using the hydrostatic pressure (ρgH) where ρ is water density at 25 • C, g is gravitational acceleration, and H is the wave height of 0.07 m. The x-axis, on the other hand, represents the location of the measurement relative to the breaking position (zero value) and normalized with the shallow water wavelength in accordance with the linear wave theory of period (T) = 0.057 s. As demonstrated in Figure 8, the swl impact pressure is at its highest at or near the breaking wave for the bare turbine case (blue circle), up to about 5.5 times the hydrostatic pressure. This result agrees with the impact pressure measured by [35][36][37] on a vertical wall breakwater. When placed on the broken wave zone (positive x-axis value), the pressure is drastically reduced until x/L = 0.1 and then rises again. This increase in pressure is observed to be the result of the breaking trajectory, where the lowest landward pressure (x/L = 0.1) occurs, as the tip of the breaking wave curves downward due to gravity and then "plunges" forward after touching the water surface (x/L > 0.1).
The ducted case (orange triangle) on the other hand shows a very steady pressure measurement along the breaking position. Both the non-breaking wave zone and the broken wave zone have a higher pressure. This indicates the duct reduced the impulsive pressure due to breaking waves from around 5.5 times to 2.6 times the ρgH, a reduction by a factor of about two. This may reduce the possibility of failure associated with breaking wave impact, which is generally considered as one of the main sources of a coastal structure's failure [38]. As pressure is expected to be lower towards a further seaward direction, the exploration for the ducted case is stopped at x/L equal to 0.1. In addition, while water will flow past the turbine (unlike a breakwater), the duration of the impact is also crucial in determining how much energy can be harvested. Based on these results and to simplify the analysis, the observation will be divided into three regions: Near-breaking wave (x/L < 0), broken wave (x/L > 0), and breaking wave (x/L = 0) zones. Figure 9 shows the impact measurement time-series data of Figure 8 at x/L = 0 of the no duct case. The rise time (t r ), defined as the time from minimum pressure to maximum pressure, of 0.003 s, can be considered as between severe to moderate impact [39], with the quasi-static loads following behind the impact. Given the reaction time of the turbine is much larger than t r [40], the turbine may not be able to utilize the momentum before water flow passes the turbine. This is another important function of the duct, as it is expected to enable the reduction of impact, which may lead to a shorter rise time, giving the turbine more time to react to incoming wave flow. It is important to identify if the reduction of impact pressure leads to an increase in energy exposure to the power take-off mechanism. For this comparison, the energy (Joule) per m 3 of water per wave is calculated based on the area under the pressure time series for each wave cycle. The pressure time series plotted in Figure 10 appears to support this statement, as the figure is arranged similarly to Figure 9, with the bare turbine case (blue circle) and the ducted case (orange triangle). The ducted case measurements indicate that the turbine is exposed to almost twice the energy of a bare turbine. There is a significant increase after x/L = 0.15 in the bare turbine case, as the same location could not test the ducted case due to the pressure sensor touching the seabed in the small-scale model. It is, however, expected to have a similar increasing trend in practice as the duct should act like a concentrating channel. It is also observed that the significant increase in impact pressure shown in Figure 9 does not lead to a higher energy exposure to the turbine. Cases observed so far are around the still water level. However, because of the ocean tide, the location of a deployed device relative to the swl will change over time. Further exploration on multiple vertical elevations is important and can be seen in Figure 11 (at the breaking wave location) for the no duct case. The y-axis in the figure represents the normalized vertical elevation divided by wave height (H), while the x-axis shows the pressure measured, normalized by ρgH, where ρ denotes water density at room temperature and g represents gravitational acceleration. The circular blue line represents the highest impact pressure measured in each case, with the highest pressure measured around y/H br = 0.4. It is important to note the zero-valued wave amplitude represents the still water level with a maximum pressure measured up to 8.5 times the hydrodynamic pressure. The exploration is then done for multiple locations for both near-breaking wave regions (negative x/L) and broken wave regions (positive x/L). The results can be seen in Figure 12, although, unlike the previous figure, the graphs are shown together with the same range for both the x-axis (p meas /pgH) and y-axis (vertical elevation/H). Furthermore, in addition to the maximum pressure measured (blue circle), the average of the 10 pressures (green rectangular) is added. The locations of each measurement (x/L) are represented below the y-axis of each figure. A simplified breaking wave illustration is added to the figure and intended as a visual aid and not as a representation of the actual breaking wave measured during the experiment. Both the x-axis and the y-axis range are kept the same for all of the cases. As demonstrated in Figure 12, in the near breaking wave zone, the maximum wave pressure (blue circle) occurred around the wave amplitude of y/H br = 0.5, with the highest pressure appearing to be up to seven times the hydrodynamic pressure at around x/L = 0.25. The average of 10 waves of pressure (dashed green square), however, shows a consistent maximum pressure of about 5ρgH for both cases. The pressure profile, furthermore, is quite different for the near-braking wave when compared to the broken wave zone. At breaking position, the pressure has a wider higher-pressure profile starting from the swl up to y/H br = 0.6, while in the near-breaking zone, the significant landward pressure is mainly measured only around the wave amplitude. This near breaking condition may create a twisting torque in the wave energy device during operation, which may lead to the device overturning or fatigue failure. In the broken wave zone, on the other hand, the energy dissipates and spreads relatively more uniformly along the vertical elevation. Figure 13 shows pressure measurement once the duct is incorporated and is displayed in a similar way. Three measurement locations were carried out, to represent: Nearbreaking wave zone, breaking wave zone, and broken wave zone for both the maximum pressure (solid orange line diamond) and the average of 10 pressure measurements (dashed grey asterisk). The pressure measured shows a more uniform profile above the still water level (swl) for every location. The x-axis, which shows the normalized pressure measurement (p meas / ρgH), is kept in the same range as Figure 11 to make the comparison clearer. This shows that the duct managed to reduce the incoming impact pressure at all positions and also managed to redistribute the pressure in the near-breaking wave zone more evenly. In the near-breaking zone, the duct reduced the impact pressure received by the turbine by a factor of 1.9 from about 7 to 3.6 times the hydrostatic pressure. In the breaking wave position, the duct is even more effective in suppressing the impact pressure by a factor of 2.2, from 8.6 to about 3.8 times the hydrostatics pressure. In the broken wave zone, the duct had minimal impact on the pressure as both the no duct case and the ducted case show similar pressure measurements of around 5.5 ρgH. This is expected in this zone because most of the wave's energy has been dissipated. Furthermore, looking at the difference between the maximum pressure measured and the average 10 highest waves of the same case, one can see that the duct could reduce the unpredictability of the impact pressure experienced by the power take-off mechanism. Moving from the impact pressure, Figure 14 shows the wave energy expected to be experienced by the turbine per m 3 of water (Joule/m 3 ) (x-axis) at the breaking wave position at various vertical elevations (y-axis). Here, the duct appears to increase the turbine's exposure to the wave energy when compared to the bare turbine case. If the 0 (still water level) ≤ y/H br ≤ 0.5 (wave amplitude) region is considered the main operational elevation of the device, then it can be shown that incorporating the duct into the design will increase the exposure of the turbine to 1.6 times more wave energy when compared to a bare turbine. In addition to the breaking wave position, Figure 15a,b show the results of the experiments for the near-breaking wave zone and the broken wave zone. The duct appears to increase and focus the energy so that the energy exposure per m 3 of water is higher near the still water level but may reduce quickly as it goes passed the y/L = 0.5 elevation. Nevertheless, in operational elevation, the duct may expose 1.3 time more energy to the turbine in both the near-breaking wave zone and the broken wave zone. This additional energy exposure can occur due to the duct not only reducing the impact pressure, but also lengthening the period where the water flows past the turbine as shown in the pressure time series data. Figure 16a shows the pressure time series for the no duct case where the impact rise time was very high and dropped immediately, followed by the quasi-static pressure from the wave crest before the pressure became flat due to water falling below the turbine location during the wave's trough. Once the duct was added, see Figure 16b, the pressure flattens for longer durations (red arrows), thus increasing the potential of energy absorption by the power take-off.

Wave Angle (β) Effect to the Performance of the Ducted OIST WEC
The measurements result of the average of 10 highest horizontal wave pressures measured can be seen in Figure 17. As expected, the pressure significantly reduces once the wave reaches the device at an angle. Interestingly, both the ducted case (orange triangle) and the no duct case (blue circle) exhibit similar results after just 10 • of wave angle. Relatively low landward wave pressure was measured at β > 60 • .  Figure 18a shows the average wave pressure measurements for both the ducted case and the no duct case with the average pressure normalized by the maximum average pressure of each case. A cosine function line (blue) is added as theoretical guidance as the pressure experienced by the turbine at an angle should be a cosine function and therefore should follow this line. The measurements show that the no duct case appears to consistently follow this line. For the ducted case, on the other hand, it follows the line up to around β = 40 • , then it begins to decline, meaning the reduction in wave energy exposure to the turbine is more than anticipated when the duct is added to the overall device after this critical angle. However, it is important to note that less impact pressure does not equate to less energy exposure. Figure 18b shows that in terms of energy exposure per m 3 of water the duct still provides an advantage compared to the no duct case by a factor of two up to about β = 40 • and declining up to β = 60 • . Any angle larger than this is expected for the duct to have a negative impact on the device's performance because there may not be enough water to enter the duct's inlet, thus limiting the waterflow through the turbine. Figure 18. (a) The normalized pressure measurement for the no duct case (blue circle) and ducted case (orange triangle) at various wave angles (β). A theoretical cosine line (blue line) is added as a guidance; (b) the energy per m 3 of water measurements for the no duct case (blue circle) and the ducted case (orange triangle) at various wave angles (β).

Conclusions
A challenging wave energy converter design placed in a breaking wave zone was introduced in 2016. This WEC project, run by the Okinawa Institute of Science and Technology Graduate University (OIST), aims to take advantage of the breaking wave energy to harness electricity. The directionality of the incident wave plays a big role in a device's performance, leading to the placement of an energy converter in a place where the wave typically travels in a narrow band of directional wave (toward the coast). This may lead to a more specifically designed device, enabling higher efficiency in the future. The initial testing of the device, including a full-scale prototype deployment in Maldives has been described in Section 2.
A campaign of experimental exploration was carried out to comprehensively observe the effect of adding a duct structure to the early design of the OIST Wave Energy Converter device. The exploration was divided into three regions of near breaking, breaking, and broken wave zones.
For the near-breaking zone, it was observed that adding the duct into early designs of the WEC not only reduces the impact pressure by a factor of 1.9, but also increases the wave energy exposure to the turbine by 1.3 times at the main operational elevation of 0 (still water level) ≤ y/H br ≤ 0.5 (wave amplitude). In the broken wave zone, the duct had less effect on the reduction of the impact pressure; however, it managed to increase the energy exposure by 1.3 times to the power take-off mechanism, similar to the near-breaking location results. In the breaking wave zone, a 1.6 higher increase in wave energy exposure was recorded with the ducted case compared to that of the no duct case, while the duct also managed to reduce the impact pressure by a factor of 2.2.
Continuing the wave angle exploration, both the no duct case and the ducted case were exposed to multiple directional waves. The test was done for the breaking wave position at still water level (y/H br = 0). The results show that the no duct case follows a cosine function as expected with the impact pressure significantly reduced once the wave hits at an angle, even as small as 10 • , and becomes similar to the ducted case results. For the ducted case, however, it is found that 40 • wave angle (β) was the critical angle where the performance becomes significantly worse at any angle greater than this. The duct, however, still provides a beneficial influence on the device performance up to β = 60 • . From this, one can expect the duct may prevent the waterflow to the turbine at any degree greater than that.