Capture width of the three-float multi-mode multi-resonance broadband wave energy line absorber M4 from laboratory studies with irregular waves of different spectral shape and directional spread

A moored multi-body line absorber is an attractive option for offshore wave energy conversion. Laboratory studies have beenundertaken to determine capturewidthwith multi-mode excitation and heave resonance for the three-float system M4 where the adjacent float spacing is about half a typical wavelength giving anti-phase forcing. The floats increase in diameter and draft from bow to stern and the bow and mid float are rigidly connected by a beam. A hinge with a damper above the mid float absorbs power from the relative rotation between the bow/mid float and the stern float. The resonant heave frequency for each float is different. Anti-phase surge forcing between mid and stern floats is substantial, while there is no hydrostatic stiffness producing resonance. This represents a hydrodynamically complex system and the laboratory experiments indicate high overall capture widths in irregular waves across a range of peak periods without damping optimisation. With different spectral peakedness and directional spread, the capture width is greater than 20 % of a wavelength (based on the energy period) across a range of peak periods typical of an offshore site for floats with a rounded base. The maximum capture width was about 37 % of a wavelength with rounded base floats; having rounded rather than flat bases increased energy capture by up to 60 % by reducing energy losses due to drag. For floats with flat bases comparisons with a geometrically scaled device five times larger and with similar magnitudes of equivalent damping showed similar capture widths as a proportion of wavelength. B Peter Stansby p.k.stansby@manchester.ac.uk 1 School of Mechanical, Aerospace and Civil Engineering, University of Manchester, Manchester M13 9PL, UK


Introduction
Wave energy is a massive renewable energy resource, but effective conversion has proved difficult largely due to its irregular nature. First, we give some general background. Point absorbers are single devices that move in heave, pitch or surge or some combination. Resonance amplifies power generation such that the theoretical maxima in terms of capture width (length of wave crest converted) are 1/2π , 1/π and 1/π wavelengths, respectively, e.g. Falnes (2002). Examples of heaving point absorbers are Archimedes Wave Swing (Polinder et al. 2004) and CETO (Caljouw et al. 2011); examples of pitching devices are PS Frog (McCabe et al. 2006) and Oyster (Cameron et al. 2010). Response to resonance is generally narrow band, although this is normally broader for pitching devices and may be broadened through latching control, e.g. Babarit and Clément (2006). To be effective most devices are designed to function as arrays, either separately tethered to the bed or from a fixed platform, e.g. Manchester Bobber (Lok et al. 2014;Stallard et al. 2009), WaveStar (Hansen and Kramer 2011) or from a floating platform, e.g. Langlee (Lavelle and Kofoed 2011). In another form of point absorption, wave motion is transmitted to an air column, driving oscillatory air motion through a turbine, usually a Wells turbine rotating in one direction, e.g. Masuda et al. (1995) and Ocean energy (2014). Another concept is based on line absorption with Pelamis (Yemm et al. 2012) a well-known example. Here, the device consists of a number of longitudinal cylindrical segments, aligned with the wave direction, connected by hinges at which power is taken off. A segment is typically half a wavelength long, so forcing in pitch is maximised. This device is driven principally by buoyancy. The device is floating with a mooring and is usually about two wavelengths long. This has the potential to exceed the capture widths of a single point absorber. A different form of line absorber known as Anaconda (Farley et al. 2012) has the form of a flexible submerged tube designed so that a bulge of water in the tube forms due to the wave pressure and travels at the wave speed, effectively in resonance. There are several good reviews available, e.g. Cruz (2008), Drew et al. (2009), Falcão andde (2010), Falnes (2002), and general hydrodynamic theory for wave-body interaction is developed in Mei (1989).
Some basic principles are apparent. Wave energy resource is greatest offshore which implies relatively deep water deployment (greater than 20 m say) and hence floating wave energy converters which require mooring. Single devices are of limited value for large-scale generation. Resonance is desirable to optimise power generation, but this is a narrow band process for a single mode and geometry. In terms of engineering practicality, floating moored systems are relatively easy to deploy and maintain relative to fixed supporting structures. Power take-off systems may take various forms with hydraulic systems quite mature, although accessibility above water level is desirable for maintenance.
The present study aims to evaluate a line absorber concept satisfying the following criteria: 1. Automatic alignment with wave direction 2. Forcing in several modes with power contributions combined 3. Multiple resonance with significant power capture over a range of wave frequencies 4. Low damping due to drag and sloshing (due to a body submerging and emerging in water) 5. Single hinge for power take-off (PTO) accessible above deck 6. High capture widths.
An overview of the design process for the system and the range of experimental tests conducted are described in the following section. The experimental results at small scale (geometric scale of approximately 1:40) with devices comprising both flat and rounded base floats are then presented. Experimental results at a five times larger scale with flat bases follow allowing comparison. The results and design are then discussed before conclusions are drawn.

Design
Three cylindrical vertically axisymmetric floats with circular cross section of different sizes are used with the smallest float for the bow and the largest for the stern as shown in Fig. 1. With the bow float moored, the largest drift force is applied to the stern float and this is beneficial for alignment with the wave direction. Drift force is a second-order inviscid effect, with a mean value in regular waves and slowly varying in irregular waves. For a given float motion, kinetic energy variation and hence average power increase with mass of float. The float diameter should however be less than about 20 % of a wavelength to be in the inertia regime capturing the phase of a wave and negligible diffraction will be desirable for deployment in an array. The stern float should have heave resonance around a prominent wave frequency; this is mainly dependent on draft. The mid and bow float are rigidly connected by a beam to generate a moment relative to the large stern float. This hinge may be above the mid float or at some point between the bow and stern float. Tests have shown that the position above the mid float centreline is desirable and results will only be presented for this case. A single hinge is advantageous in that the PTO will be concentrated in a single hydraulic unit. The mid float should have a draft which gives resonance within the band of prominent wave frequencies; this should be different from that for the stern float to generate power across the range of wave frequencies for a given site. However, the mid float should not be too big for diffraction to reduce wave action at the stern float reducing its response. The idealised case of two floats acting solely in heave is studied to guide a suitable diameter for mid float. The analytical formulation is shown in the Appendix. A natural period difference between stern and mid float of 20 % is chosen giving drafts in the ratio 1.58:1. Results of this formulation show that maximum power is achieved with corresponding diameters in the ratio 1.33:1. Further optimisation with three floats was not undertaken. Anti-phase heave forcing on the combined bow and mid float will generate pitch. Pitch moment on individual floats will also be generated, but will be relatively small. Surge forces will be substantial, but without a hydrostatic restoring force and at least in a direct linear sense they will not provide a resonant response. The resulting moment about the hinge will be dependent on the hinge position above water level. The moment will be generated by the forces on the combined bow and mid float and the stern float. The force on the bow float will be in anti-phase with the mid float, thus reducing the combined force. The surge force on the bow float should thus be as small as possible to maximise surge moment; this is consistent with what has been already proposed.
Some linearised mathematical modelling has been undertaken in Stansby et al. (2015) for regular waves using diffraction coefficients for heave, pitch and surge from WAMIT (2004) in a time stepping formulation. The floats had flat bases, and imposing drag forces as a function of the velocity of each float with a drag coefficient of about 1.8 gave power output corresponding with the larger-scale experiments (about 1:8 scale). It was shown that removing heave forcing reduced power out by 46 %, removing pitch forcing (on individual floats) reduced power by 12 % and, significantly, removing surge forcing reduced power by 57 % for a typical case near resonance. This reinforces the importance of surge forcing and heave resonance and the contribution to power output from different modes being combined. With zero drag coefficient, the maximum capture width was about 75 % of a wavelength which is clearly a high value and may be compared with three floats operating optimally in heave (48 %) and in surge or pitch (96 %), all hypothetical cases. With further geometry optimisation, the maximum could almost certainly be increased and minimising drag is clearly beneficial. Some experimental results were presented in Stansby et al. (2015). Capture widths were greater than 25 % of a wavelength in regular waves and greater than 20 % of a wavelength in irregular waves across a broad range of wave periods This paper presents the experimental results for this system obtained with two different base shapes and at two scales in irregular waves. Capture width is considered to be the vital hydrodynamic criterion for the performance. To optimise the geometric configuration preliminary tests were undertaken in regular waves and the effect of relative float spacing, hinge position and bow float size on power have been investigated. Power was then measured with irregular waves with a JONSWAP spectrum with two spectral peakedness factors for unidirectional waves and with two directional spreading factors. Comparisons are made with tests at a five times larger scale with a similar equivalent mechanical damping. Survivability in extreme waves will be considered separately.

Experimental arrangement
The smaller-scale experiments were undertaken in a wide laboratory flume with 0.45 m depth, 5 m width and 18 m length. Waves are generated by an Edinburgh Designs wavemaker with eight segments enabling regular and irregular waves with directional spreading to be generated using the OCEAN software (Rogers and Bolton King 1997). For each wave condition waves were measured at the device location without the device in position. Wave periods were in the range 0.8 to 1.65 s and the geometric scale was about 1:40.
Two device configurations were tested. The first has a flat base with slightly rounded corners as already shown in Fig The hinges comprised simple plastic bearings which will be shown to have very low friction. This is important since small forces were measured. The force on the actuator was measured by a strain-gauged load cell measuring axial force in the range ±20 N, and rotation of the beam about the column was measured by an angular encoder HEDS-9000, code wheel W2000M with 2000 counts per revolution (accurate to 0.18 • ). The force in the load cell is converted to moment about the hinge point and power is simply the actuator force times the piston velocity or equivalently the moment about the hinge times the angular velocity. Ballast is placed in the floats to provide the desired drafts. The mass distribution for each float and components is shown in Table 1 with inertias about the hinge point. This information is also necessary for mathematical modelling. The heave natural periods for each float are shown in Table 2 with added mass from linear diffraction theory provided by WAMIT (2004). The roll period of the combined system acting as a rigid body is close to 1.05 s.
The flat bases were originally chosen because heave forcing is due to oscillating wave pressure on a float base which  increases with proximity to the water surface; the flat base enables closest proximity for a given submerged volume. However, drag forces dissipate energy and floats with flat bases have higher drag coefficients than floats with more rounded bases. For a given submerged volume, the draft with more rounded base will be greater and the exciting force less.
To investigate their relative significance the mid float was given a hemispherical base of radius 0.15 m and the stern float was given a rounded corner also of 0.15 m radius, so that there was a flat central base of 0.1 m diameter as shown in Fig. 2. The mass distribution for each float and components is shown in Table 3 with inertias about the hinge. The heave natural periods for each float are shown in Table 4 and are only different from the flat base case due to different added mass. The roll period of the combined system was again close to 1.05 s.

Results at smaller scale
Results are mainly presented as average power generated for different wave height and period with average power converted into capture width as the average power generated divided by the wave power per metre width of crest; this width divided by a wavelength is referred to as capture width ratio. The formula for wave power per metre for regular waves based on linear theory is given by where H is wave height, h is depth, k is wave number, c is wave velocity and c g is group velocity, and for irregular waves is given by summation of the power of regular wave components with amplitude defined by a JONSWAP spectrum. To determine the capture width ratio in irregular waves, the wavelength associated with the energy period T e is most appropriate, equal to 0.84 T p for γ = 3.3 and 0.78 T p for γ = 1 where T p is the peak period and γ is the spectral peakedness factor of the JONSWAP spectrum.

Preliminary regular wave tests
In the smaller-scale experiments, reflections were evident with regular waves making absolute magnitudes difficult to assess. Reflection was relatively low in the larger-scale tests, approximately 6 % without the body present in the basin and are reported elsewhere (Stansby et al. 2015). Regular wave results at smaller scale are thus not presented, although relative power magnitudes were used to assess the geometric configuration. The power with a smaller bow float diameter of 0.16 m was compared with that due to 0.2 m. Both diameters gave very similar power output overall, but the larger diameter had greater buoyancy to support the connecting beam which would be beneficial at full scale. The effect of reducing the spacing between the bow and mid float was tested and 0.8 m gave higher power outputs than 0.7 and 0.6 m. The effect of hinge elevations on average power at 0.15, 0.21 and 0.25 m above mean water level was tested and the 0.21 m elevation overall gave marginally larger overall power output and was subsequently used.
The variations of average power output with H 2 showed a closely linear dependence for different periods, which may be expected since wave power is proportional to H 2 . The straight line fits for different wave frequencies passing close to the origin, indicating that power losses due to friction in the bearings are very small. Such small values are difficult to   quantify accurately and will be ignored; the power estimates are thus slightly conservative.

Irregular waves
Capture width ratios (with wavelength defined by T e ) are shown in Fig. 3 for the rounded base with a flat base comparison for the narrow band JONSWAP spectrum with γ = 3.3. Note the mechanical damping was not adjusted and will be shown to be suboptimal. There is a single peak close to T p /T r3 = 0.9 where T r3 is the heave resonance period on the stern float with added mass defined by WAMIT (2004) and T p = 1.0 s (T e = 0.84 s). The capture width ratios are always greater with the rounded bases with a maximum of 0.32 rather than 0.20 with significant wave height H s ≈ 0.04 m, a 60 % increase. Note with H s ≈ 0.027 m the maximum ratio is 0.37, indicating a smaller drag effect with the smaller wave height. The H s values given in the figure captions are those measured without the device in position and are an approximate value across a frequency range, as the target value was not produced directly. The measured values for each test are used to specify the non-dimensional capture width ratio. An example of time variation of moment M d , angular displacement θ and power output is shown in Fig. 4 for H s = 0.043 m, T p = 1.11 s and γ = 3.3. In a perfect damper, the coefficient B d would be constant but this was not the case. The M d time variation is fitted by the least squares method to a model defined as M fit = B 0 + B dθ + B aθ , where B a allows for inertia effects and B 0 accounts for the small mean value observed. The least squares goodness-of-fit R 2 of 0.9 is a typical value. The B d value was 1.65 Nms in this case and was in the range 1.4-2.3 Nms for all tests.
The effect of reducing the spectral peakedness factor γ from 3.3 to 1.0 is shown in Fig. 5 for the rounded base floats. The capture width ratio is reduced for the smaller γ and is again smaller for the larger H s , although this is a small effect with γ = 1.
The effect of directional spreading is shown in Fig. 6 for γ = 3.3 also for the rounded base floats. The capture width ratio is largest for uni-directional waves (s = ∞) and is reduced with s = 30 and 5 which are similar. Again, capture width ratio is greater for smaller H s . The directional spreading parameter s defines a spreading function of cosine shape, cos s (θ ), proposed by Borgman (1969) following Longuet-Higgins et al. (1963), and is the standard in OCEAN software.  Surprisingly with γ = 1.0 as shown in Fig. 7, capture width ratio increases as directional spreading increases (s reduces) and again magnitudes are larger for smaller H s ; the values for s = 5 and 30 are very similar. For both γ values, the capture width ratio is greater than 0.2 across a broad range of periods.

Results at larger scale with flat base floats
Experiments at five times larger scale were undertaken in the University of Plymouth wave basin using floats with flat bases, shown in Fig. 8. Note the radius of curvature of the rounded base corners was 10 cm, double the relative radius of the smaller scale. The wavemakers were again from Edinburgh Designs driven by the OCEAN software; the depth is 2.9 m and mass and inertia distributions are given in Table 5. With a geometric scale factor of five applied, the mean  (Stansby et al. 2015) indicate that the maximum power capture at larger scale is achieved by increasing B d by a factor of between 2 and 4 increasing power capture by 24 %. Since B d values do not exactly scale the comparison between scales can only be approximate.
In the experiments with regular waves, the variation of average power output is again close to proportional to H 2 for different periods. Best-fit lines pass very close to the origin indicating negligible friction losses in the bearings.
For the flat base floats, Fig. 9 shows that the variation in capture width ratio is very similar at small and large scale with γ = 3.3 for the smaller equivalent H s values. For the larger equivalent H s values, the capture width ratio is somewhat smaller at small scale. At the larger scale, curves are quite similar for γ = 1.0 and 3.3.
An example of time variation of moment M d , angular displacement θ and power output is shown in Fig. 10 for H s = 0.23 m, T p = 2.8 s and γ = 3.3. The M d time variation is again fitted to a model defined as The least squares goodness-of-fit R 2 is 0.94 in this case; the B d value was 2536 Nms.
Capture width ratio from tests at both scales with directional spreading s = 30 are shown in Fig. 11 for γ = 3.3 and in Fig. 12 for s = 5. It can be seen that values with equivalent larger H s are somewhat larger at the larger scale with s = 30 although quite similar with s = 5. The values for the smaller wave height, tested at smaller geometric scale only, are generally slightly larger.

Discussion
Experimental analysis is undertaken with geometrically scaled models for which time and force scales are based on Froude scaling. The scale of the systems studied differs by a factor of five. The moments generated by the mechanical actuators were fitted reasonably by terms proportional to angular velocity and acceleration and a mean term, enabling the damping factor to be determined. The smaller-scale pneumatic actuator gives a wider range of equivalent damping factor values than the larger-scale hydraulic actuator. Comparisons between these experiments may therefore only be approximate. Damping factors at both scales are suboptimal based on prior analysis in regular waves for the flat base floats at the larger scale which implies that capture widths may be 24 % larger. However, results in regular waves do not  necessarily transfer linearly to irregular waves or to floats of different shapes. Control of mechanical damping is of course desirable to enable more precise matching between scales as well as optimisation, although this was not possible in these tests and is difficult at small scale. Nevertheless, it is notable that qualitative similarity has been observed using quite different damping actuators. The instantaneous power was measured without approximation providing accurate average power estimates. The experiments are thus valuable for exploring the power conversion potential of this system with novel features intended to give a marked improvement in capture width over existing configurations.
The size of the stern float defines the dimensions of the rest of the system and the diameter is specified as approximately 20 % of the wavelength associated with a prominent wave period, so that forcing is predominantly inertial with minimal diffraction. The natural heave period of the stern float is designed to correspond approximately with the peak spectral period by specifying ballast and hence draft. The spacing between floats is set at about half the corresponding wavelength. The mid float is intended to have a somewhat lower natural period to correspond with lower wave periods (about 20 % lower in this case) and the diameter is chosen to maximise mean power. This analysis was based simply on two heaving floats giving a diameter of three-quarters of the stern float diameter. Power is only slightly dependent on diameter for larger values. Clearly, a smaller diameter would also be advantageous for minimising fabrication material at full scale. The bow float was originally intended to provide a passive mooring point with sufficient buoyancy to support the beam to mid float. The larger of the two diameters tested was chosen to give more buoyancy without loss of overall power capture and it was not used as a mooring buoy, since it has significant motion contributing to power generation. The change in base shape on mid and stern float from flat to hemispherical and rounded with the same radius of curvature, respectively, gave a dramatic improvement in capture width emphasising the importance of drag minimisation.
The maximum capture width of over 35 % wavelength with rounded bases in irregular waves and importantly over 20 % over a broad range of peak periods for all values of spectral peakedness (γ ) and directional spreading(s) tested may be considered high values. The geometric configuration could probably be optimised further; a larger stern diameter should be considered which will influence the optimum diameter of the mid and possibly the bow float. In determining these natural heave periods, it was found that viscosity affects added mass as well as determining drag. The natural period of the stern float with the flat base in isolation was measured to be 1.12 s in a free decay test, while that calculated with added mass from diffraction analysis was 1.07 s, a difference of 4.6 % implying that 20 % of added mass may be associated with viscous effects. Note that the mid and bow floats could not be tested in isolation due to lack of stability.
In irregular waves there is a flat peak in the variation of capture width with T p with a maximum just below the heave resonant period of the stern float based on diffraction added mass. The peak is thus closer to the natural period accounting for viscous effects. That the capture width only decreases slightly for smaller T p is probably due to the lower resonant period of the mid float and it would be useful to assess optimisation of this draft further to maximise overall energy capture in irregular waves for a representative wave climate.
Reducing spectral peakedness γ from 3.3 to 1.0 reduces capture width with the rounded base floats, but has little effect with the flat base floats. Smaller H s values generally give slightly higher capture widths indicating smaller drag effects. Increasing directional spread (reducing s) reduces capture width with γ = 3.3, although there is little difference between s = 30 and s = 5. Perhaps, surprisingly directional spread results in slightly larger capture width with γ = 1.
Results have not been represented for regular waves because of uncertainties in capture width associated with reflections and occasionally cross-basin seiching, although this was not observed at larger scale and reflection was a small effect; results are reported elsewhere (Stansby et al. 2015). This is a recognised limitation associated with tank testing, particularly for multi-body systems. For small arrays of heaving floats at similar scale, the response amplitude was obtained from the initial cycles only (Weller et al. 2010).
Mathematical modelling has not been applied in this paper which is a purely experimental investigation in irregular waves. Several different modelling approaches are appropriate. Modelling based on linear diffraction coefficients for added mass, radiation damping and excitation for all relevant modes will provide an efficient approach in the frequency or time domain, but an effective drag coefficient dependent on wave height and possibly period is needed to calibrate against experiment. The effect of viscosity on added mass also needs calibrating and there may be an effect on radiation damping. This modelling is thus a valuable adjunct to experiments for fast optimisation, but is of limited value in isolation. Viscous effects are included in full Navier-Stokes solvers with a free surface, typically RANS solvers in the volume of fluid (VOF) method. Such methods with a complex fluid mesh are exceptionally demanding computationally and subtle effects of boundary layer separation with turbulence which affect drag are very difficult to predict accurately. Nevertheless, further investigation through this type of modelling would be beneficial.
While further geometric optimisation is desirable, the optimisation and design of the power take-off (PTO) is expected to be of greater importance. This is actually very difficult to achieve physically at small scale and here a calibrated mathematical model is desirable to couple with algorithms for hydraulic control (based on proven hydraulic systems) to enhance and smooth power output.
For an offshore site with a known wave climate, the characterisation of the system in terms of capture width enables annual energy yield to be readily determined. While we are concerned with capture width, the device cost per MW capacity and levelised cost of energy determine practical viability. A simple preliminary economic analysis  indicates that capital costs could be comparable to target costs for this stage of development given in Marine Energy Technology Roadmap (2010); this is competitive with offshore wind energy and without the associated environmental intrusion.

Conclusions
A three-float line absorber M4 has been proposed with float spacing of about half a wavelength to give anti-phase forcing between adjacent floats with excitation in heave, pitch and surge. The float size in terms of both diameter and draft increases from bow to stern, facilitating alignment of the line absorber with the wave direction. The heave resonance periods on mid and stern float are different to optimise power capture for a range of wave periods and the mid float diameter is sized to give maximum power output (based on an idealised heave analysis). Anti-phase surge forcing on mid and stern floats is significant, although without a hydrostatic restoring force for resonance. This forcing is exploited by the relative pitch about the hinge point well above deck level. The aim is that power from the different forms of excitation and forcing combines to give a markedly improved capture width over what has been achieved previously in wave energy conversion. This has been investigated at laboratory scales of about 1:40 and 1:8 in irregular waves with different spectral peakedness and directional spreading. Capture width is the important assessment criterion. In irregular waves, the capture width is above 20 % of the wavelength associated with the energy period, for a configuration with rounded base floats, across a broad range of wave periods and with different combinations of both directional spreading and spectral peakedness typical of an offshore site (with appropriately scaled periods); the maximum was over 35 %. This may be considered large relative to existing systems.
The tests indicate the effect of viscosity on both drag and added mass and indicate the importance of minimising drag for increasing power output. While some geometric optimisation has been undertaken, more would be beneficial. Most importantly, optimisation of the power take-off (PTO) control is required since only simple actuators are used in these reduced-scale experiments.