Tidal Energy Round Robin Tests: A Comparison of Flow Measurements and Turbine Loading

: A Round Robin Tests program is being undertaken within the EC MaRINET2 initiative. This programme studies the used facility inﬂuence can have on the performance evaluation of a horizontal axis tidal turbine prototype when it is operated under wave and current conditions. In this paper, we present the design of experiments that is used throughout the work programme and the results related to the ﬂow characterisation obtained at the I FREMER wave and current circulating tank, the C NR -I NM wave towing tank and the ocean research facility F LO W AVE . These facilities have been identiﬁed to provide adequate geometric conditions to accommodate a 0.724 m diameter turbine operating at ﬂow velocities of 0.8 and 1.0 m/s. The set-up is replicated in each of the facilities with exemption of the amount of ﬂow measuring instruments. Intrinsic differences in creating wave and currents between facilities are found. Flow velocities are up to 10% higher than the nominal values and wave amplitudes higher than the target values by up to a factor of 2. These discrepancies are related to the ﬂow and wave generation methods used at each facility. When the ﬂow velocity is measured besides the rotor, the velocity presents an increase of 8% compared to the upstream measurements.


Introduction
The development of new technologies relies heavily on prototype testing activities. These tests are usually carried out in laboratories that are able to provide similar conditions to the real environment in controllable scenarios where certain parameters can be adjusted as required. In the case of ocean energy conversion research, the test environment may influence measurements related to power capture or loading, which, in consequence, influence the design of the device.
To identify the influence of the test environment on the device's performance, interlaboratory tests can be undertaken using the same testing device and subjected to the same testing conditions. These experiments are a common practice and are referred as Round Robin Tests (RRTs). For that programme, the exact same 0.724 m diameter horizontal axis tidal turbine (HATT) is tested in two tow and two flume tanks. Small discrepancies between average values of power and thrust are obtained giving confidence in the high-quality operation under current-only conditions of the four research infrastructures used at that time. Moreover, an insightful conclusion drawn from the work of [1] indicated that, even at low blockage ratio conditions of about 5%, the turbine performance within operation ranges between peak power and free-wheeling are impacted. To minimise the influence of blockage, correction algorithms are implemented based on the work from [2].
Other attempts to identify turbine model practicalities using RRTs have been reported in [3,4]. Murray et al. [3] compared the performance of a HATT in a flume and a tow tank using two types of blade materials: composites and solid aluminium blades. Small variations between power and thrust are obtained which are merely related to the turbulence levels developed in the flume. Allmark et al. [4] compared the effects of measuring the device performance using motor current measurements instead of commonly employed torque and thrust transducers. It is found that substantial efforts to characterise the motor need to be undertaken to minimise uncertainties between measurements. Both campaigns as well as that reported in [1] only considered the use of current-only conditions.
Efforts from the scientific community to drive forward the development of tidal stream technology include the estimation of performance and structural loading when these devices are subjected to diverse wave and current conditions e.g., [5][6][7][8][9][10][11][12]. As observed in Table 1, most of these studies have also been undertaken in diverse facilities under different conditions. It is also clear that there is a progression in this line of research, moving from the more controllable aspects; for instance, the use of regular waves to more complex situations including irregular waves and directionality.  [14] Flume tank. IFREMER 0.9 Collinear regular and irregular waves The research developed by the authors in Table 1, although crucial for the development of marine devices, may be difficult to relate due to the number of variables related to device geometry, facility, flow conditions, etc. Thus, the aim of this work is to evaluate the methodology implemented to assess the performance of a generic tidal stream turbine and the flow characteristics at each facility regarding the quality of the flows and waves generated. This paper complements the work done by [15], where the authors studied in detail the performance of the turbine through non-dimensional analyses using a technique to establish the input flow characteristics and blockage effects.

Design of Experiment
The purpose of this investigation is to identify the potential effects that research infrastructures may have on the evaluation of small scale tidal stream turbines performance. To meet this objective, three controlled facilities (Table 2) are selected to undertake this evaluation. Initially, the design of the experiments reported here attempts to follow the experimental process described in the EquiMar protocol [16] adapted from the ITTC Procedure 7.5-02-02-01 [17,18]. To describe the methodology that is used for the entire RRT programme, this paper focuses on the design of the experiments, the initial findings related to the quality of the flow and waves generated and the turbine performance obtained at three research facilities: the circulating flume tank of IFREMER [19], the tow tank of CNR-INM [20] and FLOWAVE TT [21]. Finally, an uncertainty analysis is presented to understand any deficiencies with the adopted methodology.

Model Specifications and Installation
The desired parameters required for the evaluation of the RRT are, mainly, flow characterisation of the facilities used and the loading and performance of the model tidal turbine. To meet the latter criteria, a 1:20 scale, three-bladed horizontal axis turbine model developed by IFREMER is used in this campaign [1]. The turbine is 0.724 m in diameter (D) and is equipped with three blades which are 307 mm long. This is representative of a 15 m diameter turbine at full scale. The blade design is based on a NACA 63-418 profile. A mounting flange to adapt the supporting stanchion or tower is located 594 mm (0.82D) from the rotor plane ( Figure 1). The turbine is rated at TSR = 4 and is designed to operate for a maximum flow velocity of 1.2 m/s. The turbine is fixed to a moving carriage on the towing tank (mounted from the top) and on a customised mounting frame in the flume tanks at IFREMER and FLOWAVE (mounted from the bottom). The turbine is kept at a constant depth of 1.0 m, to keep it in mid-water conditions according to the depth specifications of most of the tanks used in the RRTs. For a detailed description of the turbine setup, see [15].

Turbine Instrumentation
The generic turbine model is equipped with a dedicated load cell measuring the torque Q and thrust T applied on the main rotation axis of the rotor. This waterproof transducer ( Figure 2) is positioned upstream of the rotor seals to prevent parasitic friction effects. The torque and thrust transducer are custom made by the French company SIXAXES [22] in partnership with IFREMER. The measurement ranges of each of the instruments are summarised in Table 3. The shielded cables coming from the transducer are routed through a slip-ring enabling the free rotation of the cables preventing their entanglement. These low voltage signals are amplified by an electronic signal processing unit, located outside of the turbine and on the dry. The motor shaft is connected to the rotor shaft through a motor-gearbox (ratio 1:26) enabling suitable torque and rotation speed ratings.
Default calibrations provided by the manufacturer are used for the entirety of the programme. The accuracy of the instrumentation according to the service provider can be found in Table 3. These values are used to undertake an uncertainty analysis for all the test programmes.
The accuracy of the encoder is 1 RPM. However, due to the gearbox ratio (1:26), the accuracy of the rotor rotation is calculated as 1/26 = 0.04 RPM.

Flow Measurement and Locations
The flow stream is characterised in the absence of the turbine (bare-basin) and while the turbine is in operation (synchronised). An Acoustic Doppler Velocimeter (ADV) is selected to perform this measurements based on its light-weight and small dimensions. The ADV is placed in line with the rotor plane at a distance of 1.2 m (1.6D) along the longitudinal section of the tank (see Figure 3). At IFREMER and FLOWAVE, the ADV is mounted independently of the turbine. Seeding particles are deposited in the tank and continuous flow circulation permits a good dispersion throughout the flow stream.
For the experiments at CNR-INM, the ADV is mounted to the towing carriage. In addition, a seeding mast mounted on the carriage is placed 5.5D upfront of the ADV. The seeding mast consists of a cylinder beam of 15 mm in diameter supporting a tube of 8 mm.
Twenty openings along a vertical section of 800 mm releases the seeding particles. The default accuracy of the ADV is reported as ±1% of the measured value [23], i.e., ±0.01 m/s for a flow velocity of 1.0 m/s. However, an adequate operation of this device also depends on the density of particles in the water. In order to account for this, the signal correlation is preserved for the entirety of the experiments at a rate higher than 80%, hence the reason for using a seeding mast at the tow tank (CNR-INM). Furthermore, ADV measurements are corrected for the presence of spikes by using the method proposed by [24].
The use of flow measurement equipment in tow tanks and for tidal stream applications is scarce in the literature, especially when referring to ADV technology. Therefore, this research attempts to evaluate its applicability when capturing flow characteristics in both current-only and wave-current combinations as well as the measurement location.
The flow stream is characterised in each of the facilities with the ADV and the wave gauges, and with and without the turbine for all the testing conditions. For the bare-basin measurements, the position of the ADV is exactly the same as the rotor hub when operated, i.e., at x = 0 ( Figure 3). As presented by [15], depending on the facility, three or five ADV measurement points evenly spread along the vertical and horizontal axis are acquired over the turbine swept area in order to quantify any vertical or lateral velocity gradient ( Figure 4).  One identified challenge in the field of tidal energy is to adequately characterise the onset flow to calculate the turbine's non-dimensional parameters (see Section 2.5.2). An attempt to overcome this issue is undertaken by [25], where the authors recommended to carry out flow velocity measurements at least 2D upstream of the turbine when the turbulence rate is high. This procedure thus helped to better account for the velocity fluctuations perceived in synchronisation with the turbine. Another approach presented by [26] recommends to use the statistical mode of the velocity recorded instead of the mean velocity when calculating the performance coefficients. Using the mode of the flow velocity, differences in the performance coefficients of up to 12% are obtained when compared to the values obtained with the mean of the flow velocity.
In addition to the utilisation of an ADV, at IFREMER, the flow is also monitored using a Laser Doppler Velocimeter (LDV). While the ADV is able to measure three flow velocity components, the LDV available at IFREMER can capture two components: u which is related to the flow direction e x and v the horizontal component related to e y as seen in Figure 3. The LDV is installed 2D upstream from the turbine, this distance is chosen based on previous experiences from the authors in order to obtain a closer representation of the flow "perceived" by the rotor while preserving a non-intrusive measurement [25]. The focal point of the laser beam (i.e., the measurement point, 0.53 m from the probe) is set to be in-line with the turbine hub, at 1.0 m water depth. The accuracy of the LDV is ±0.1% [27]. This experiment is limited to the IFREMER facility due to difficulties arising with the transportation of the LDV to other laboratories.
A number of wave probes are also placed next to the turbine to verify the wave parameters set for each of the case scenarios. Three resistive wave probes are used at the facilities in IFREMER and FLOWAVE [15]: probes 1-3 as seen in Figure 3. A mix of resistive, ultrasound and dynamic wave probes is used in the CNR-INM tow tank: probes 3-6 as seen in Figure 3, where probes 4 and 6 are ultrasound wave gauges and probe 5 is a dynamic wave gauge. One single identical resistive wave probe is used in all facilities (probe 3). This is a crucial aspect of the test campaign, especially since the flow characteristics impacted by the addition of waves are deemed to be relevant between facilities.

Experimental Plan
The flow velocity or carriage velocity is set to 0.8 and 1.0 m/s. This flow velocity corresponds to a Reynolds numbers between 5.8 × 10 7 and 7.2 × 10 7 , based on the turbine diameter and are representative of full scale velocities of 3.5 and 4.5 m/s, respectively.
These flow speeds are chosen based on the capabilities of the laboratories which in some cases cannot operate at faster flows when waves are generated in combination with currents. According to [28][29][30][31], flow velocities around UK waters reach maximum values between 3.5 and 4.0 m/s and turbulence intensities ranging from 10 to 20%.
Four regular waves in-line with the current are chosen with wave frequencies between 0.5 and 0.7 Hz and wave amplitudes between 35 and 75 mm. The wave amplitude and frequency are selected based on information provided by the facilities where some can achieve better wave properties at certain frequencies than others, especially when combined with current. One irregular wave setting is included in the test plan with a maximum wave frequency of f max = 0.6 Hz and a significant wave height of H s = 0.1 m. At full scale, these wave parameters are associated with wave amplitudes between 0.7 and 1.5 m and wave frequencies between 0.12 and 0.18 Hz (periods between 5.5 and 7.8 s). According to [32], ordinary gravity waves with periods around 7 s have the largest relative height and are the ones a tidal turbine will be subjected to more frequently. Although the waves tested might be considered small compared to real sea states, the aim of the RRT is to repeat the same conditions in different facilities and to assess the impact on the turbine. Therefore, the selected wave conditions are theoretically replicable at all the facilities and, as it will be seen in Section 3.3, a periodic velocity component is generated. The summary of the tested flow conditions is presented in Table 4. For each of the test cases, eleven rotational speed settings are considered. These range from 0 to 148 RPM for flow velocities of 0.8 m/s and from 0 to 185 RPM for flow velocities of 1.0 m/s. Maximum values are close to free-wheeling operation. To quantify the uncertainty of the experiment, repeated tests are considered for each case for several rotational speeds. Due to time limitations and the long wait between tests at CNR-INM, the test matrix is constricted to five rotational speeds per case plus two rotational speeds for the repeated tests. Note that the repeated tests are performed at the end of the test programme to allow for a randomised data acquisition while ensuring a time-efficient test campaign.

Data Acquisition and Testing Procedure
All signals are acquired using National Instruments (USA) hardware and in-house electronics developed by IFREMER staff. The signals are mostly sampled at a frequency ( f s ) of 128 Hz, which allows the capture of the dynamic loads inherent from the turbine. This sampling frequency should suffice to describe the performance of the device adequately; however, larger sampling rates should be used for other vibration analysis.
Flow and water surface elevation measurements were synchronised with the turbine instrumentation by means of a short impulse trigger signal. The surface elevation is captured at the same sampling rate as the turbine signals. The ADV was, however, set to acquire data at a different rate depending on the facility, 128 Hz in the tow tank (CNR-INM) and 64 Hz in IFREMER and FLOWAVE.

Testing Procedure
To ensure that all the tests are performed similarly, the following test procedure is agreed between partners.
For recirculating flumes and tanks (i.e., IFREMER and FLOWAVE): 1. zero values are first recorded at stagnant conditions for 30 s; i.e., water still and parked turbine, 2.
the flow is set to the required speed, 3. the waves are imposed; i.e., initiate operation of the wave-makers (if applicable), 4.
set the rotational speed of the turbine (starting at zero RPM) 5.
For towing tanks (e.g., CNR-INM), the process is similar to the process described above. However, the wavemakers and turbine rotation are initiated before the carriage moves.
In the flumes, 128 s of measurements are acquired for the current-only conditions, 256 s for the regular wave experiments and 512 s for the irregular wave condition. These test durations are adopted in order to reach a convergence rate depending on the fluctuation of the signals that are reported in Section 3.6.
The time duration adopted in the tow tank differs due to the inherent finite length of the tank but also due to the homogeneous nature of the flow (i.e., no turbulence). For the tow-only and wave-tow conditions, 50 s and 140 s are recorded, respectively. Although shorter than the tests at the recirculating flume and tank, at a carriage speed of 1.0 m/s, about 20 wave periods are achievable with good convergence. For the irregular wave case, however, tow tank repetitions are gathered until 420 s of data are acquired. This time is equivalent to approximately 250 wave encounters at f = 0.6 Hz.

Flow Characterisation
The flow velocity records from the flow characterisation procedure are processed to define the incoming (far upstream) velocity U ∞ . As suggested by [33], the square or cubing of the velocity is performed prior to temporal and spatial averaging. The average velocity is then computed as a Disc-Integrated Velocity Average (DIVA) that serve as an input for the calculation of non-dimensional parameters of power and thrust. The DIVA method is described and employed in the associated work done in [15].
In this work, only analysis of the u flow component is explored due to the amount of data gathered during these tests. However, the analysis of the v and w flow components is presented in [15].

Turbine Performance
The turbine performance is characterised using three essential parameters: the mean hydrodynamic torque, Q in N.m, the mean thrust, T in N and the angular velocity of the turbine, usually represented by ω in rad/s, or as rotational speed using revolutions per minute RPM. In this investigation, these operational characteristics are evaluated and presented focusing on the repeatability of the experiment and the uncertainty analysis between tests.
Another common approach is to represent the performance of a tidal device with non-dimensional parameters of power coefficient C P and thrust coefficient C T , which are explored in [15].

Uncertainty Analysis
When possible, the full test cases have been repeated at each facility. This information is deemed necessary to undertake an uncertainty analysis based on a combined expanded uncertainty using precision and bias errors.
The precision uncertainty, u A is calculated based on the standard deviation obtained from the repeated tests (s) and the number of repetitions (n), as displayed in Equation (1). The bias uncertainty (u B ) is quantified based on the accuracy of the instrument as provided in Table 3.
The combined expanded uncertainty is thus calculated based on Equation (2), using a coverage factor k = 2 which gives a level of confidence of nearly 95% as reported by [34]. The combined uncertainty is presented as a percentage of mean value for better interpretation. This analysis is based on the recommendations made by [16,35]. Further calculations are required when looking at error propagation of other parameters; for example, when power and thrust coefficients are calculated, as it is described in [35].

Signal Stabilisation
To have a better understanding of how much a signal (i.e., thrust or flow velocity) deviates from the average of a full time series, the Root-Mean-Square Error (RMSE) is used [36]. The RMSE is calculated in 1 s windows. The units of the RMSE are the same as the units of the signal being analysed: where t is time in seconds, T is the total duration of the time series,ŷ [1,t] is the average of the signal from the beginning of the test to the value of t, and y T is the average of the signal for the whole time series.

Flow Velocity Characterisation (Bare-Basin)
In this section, the flow characteristics obtained at each of the facilities are studied using the data recorded in the absence of the turbine, here defined as "bare-basin". In Figure 6, time history at IFREMER shows slightly higher velocity ranges than the other two facilities. In the frequency domain, flow velocity shows higher energy levels. These are associated with the flow disruptions and turbulence created by the position of the wavemaker located upstream of the test area. The −5/3 power law is clearly present at IFREMER and FLOWAVE. However, the spectra measured at FLOWAVE contains lower energy. The lowest one is spectra measured at the CNR-INM facility. Similar to case 1, the energy spectra at CNR-INM still presents the ADV probe vibrations at f 5.5 Hz. These vibrations are actually thought to be associated with a combination of carriage speed, wave conditions, and the ADV being mounted below the carriage (see Section 2.2.2). In all facilities, a periodic component with the same energy level is clearly visible on the PSD, indicating the wave orbital velocities for case 2, with a wave frequency of 0.6 Hz. The velocities obtained using the DIVA method for the flow characterisation (barebasin) are grouped together in Table 5 as presented in [15]. Although flow characterisation for case 6 at CNR-INM was not performed, turbine performance was measured for this case. All the average velocities recorded are higher than the nominal values of 0.8 and 1.0 m/s. The highest value for the 0.8 m/s cases is for combined wave and current case 2 at IFREMER with average velocities ≈5.8% higher than the nominal velocity. For the 1.0 m/s cases, case 5 with regular waves at IFREMER has the highest value at 5.6% above the nominal flow velocity. At IFREMER, average velocities for the combined wave and current cases are up to 3% faster than the current-only cases. This could be attributed to the velocity depth profile and the contribution of the waves to the flow velocity. At each facility, average flow velocities between cases are within ∼0.5%.
At CNR-INM, however, it has been proven that the carriage speed has an accuracy of ±1 mm/s. The differences between the DIVA method (Table 5) and the carriage speed could be attributed to the accuracy of the velocity measurement technique used.

Flow Characterisation (with Operating Turbine)
ADV measurements are also obtained throughout the test campaign with the ADV placed at 1.2 m from the turbine axis (see schematic on Figure 3) when operating and in synchronisation. Figures 7 and 8 show the variation of the average and standard deviation velocity for cases 1 and 2, respectively, with the turbine rotational speed. Table 6 shows the time and case average of the streamwise velocity u for all cases at all facilities. Overall, average flow velocity at all facilities is higher than the nominal velocities, i.e., similar behaviour seen in the bare-basin measurements in Table 5. For the current-only case (case 1), FLOWAVE presents flow velocities 10% higher than the nominal velocity. With the addition of regular waves (case 2), FLOWAVE and IFREMER have flow velocities 10% higher than the target velocity. In contrast, flow velocity at CNR-INM is closer to the target value. Results obtained with the DIVA method (Table 5) present higher values than the ADV at FLOWAVE and IFREMER and lower values at CNR-INM. These differences are mainly attributed to blockage effects and rotor-induced velocities.
Notably, there is certain fluctuation in the flow velocity depending on the turbine's rotational speed. The influence from the operating turbine may explain the discrepancies of the results. It is worth noting that the LDV (Section 3) is located upstream of the turbine, where the rotor-induced velocity has a negative contribution to the upstream current, whereas the ADV is located next to the rotor, where the contribution of the rotor-induced velocity is positive due to the by-pass effect.  Table 6. Time and case average of the ADV stream-wise velocity u for all cases at all facilities. ADV mounted at x ADV = 0 and y ADV = 1200 mm as seen in Figure 3. To further investigate this discrepancy, a comparison between the ADV and the LDV data for current-only case 1 at all facilities is presented in Figure 9. At all rotational speed values, the velocity measured by the ADV is higher than the one measured by the LDV (2D upstream of the turbine) and than the nominal flow velocity in all facilities. As seen in Figure 3, the distance between probes is quite large: x LDV = −1448 mm and y LDV = 0, whereas x ADV = 0 and y ADV = 1200 mm. As the rotational speed increases, the LDV measurements decrease and the ADV measurements increase at IFREMER and CNR-INM. At FLOWAVE, however, behaviour is the opposite. As the rotational speed increases, the flow velocity decreases. At RPM = 0, the LDV has the same measurements as the nominal flow velocity. However, the ADV measurements are about 3% higher at CNR-INM, 5% higher at IFREMER and 11% higher at FLOWAVE. At the highest RPM values, the LDV measurement is 2% lower than the nominal velocity, whereas the ADV measurement is >6% higher at all facilities. Due to these differences, it is concluded that the synchronous measurement obtained with the ADV on the side of the turbine, especially with high blockage ratios, is influenced by the presence of the turbine and its rotation effect in all facilities. It is related to the bypass of the flow around the turbine increasing with the turbine rotation.

Wave Orbital Velocities and Surface Elevation
Due to the fluctuating characteristics of the waves obtained (Figures 8), similarly to [15], the instantaneous values of the velocity time series of the LDV (Figure 3) are analysed using the Hilbert transform. A Fast Fourier Transform (FFT) is also performed to quantify the amplitude of the waves. Figure 10 shows the amplitude of the Hilbert transform of the stream-wise velocity u measured by the ADV in synchronisation with the turbine parameters for case 2, presented against the rotational speed. The time-average of the horizontal wave orbital u for all the wave cases is shown in Table 7. The average and standard deviation of the Hilbert transform is higher at IFREMER than at CNR-INM by 37%. At FLOWAVE, measurements are only 8% higher. The higher average velocities at IFREMER are due to the higher wave amplitudes ( Figure 11). However, the high standard deviation values at IFREMER and FLOWAVE are due to interaction between wave and current.  The free surface elevation measured with wave probe 3 ( Figure 3) in synchronisation with the turbine parameters for case 2 using the Hilbert transform is presented in Figure 11. As seen previously in the flow velocity data in Figure 6, the wave probe data shows that, for all wave cases, the average wave amplitude is higher at IFREMER. Comparing the results with the nominal values (see Table 8), average wave amplitudes differences at IFREMER range from −20 to 16%, from −30 to −9% at CNR-INM and from −53 to −20% at FLOWAVE, versus the case number. The inability to accurately account for such broad ranges is thought to be due to the difficulty it entails to generate waves in recirculating tanks [37]. The use-or lack thereof-of wave amplitude corrections in the presence of currents or the wave generation techniques deployed in an incoming flow can explain such variations. At IFREMER, waves are generated with the addition of a wavemaker located upstream of the measurement area, at CNR-INM waves are generated in steady water and at FLOWAVE a circular array of paddles behind the current turning vanes generates the waves. In addition, uncertainties related to the resistive wave gauges may play a minor role (in the order of millimetres) in the different wave amplitudes seen between facilities. Ifremer Cnr-Inm FloWave η = 87 mm η = 66 mm η = 60 mm target: 75 mm Figure 11. Variation of the average (markers) and standard deviation (shading) of the wave amplitude η with rotational speed for case 2 at all facilities. Dashed lines are the average of the wave amplitude for case 2 at all rotational speed settings per facility. The red line indicates the target wave amplitude of case 2, data obtained with wave probe 3 (see Figure 3). Table 8. Time and case average of the wave amplitude η for all cases at all facilities. Data obtained from wave probe 3 (see Figure 3) in synchronisation with the turbine. Another method to calculate the wave heights is to use zero-up crossing techniques. In Figure 12, a visual representation of the zero-up crossing method for calculating the wave height H (i.e., η × 2) is presented. Firstly, a Loess filter using a span of 0.09% of the total number of samples is used to remove high frequency noise. Once the data are smooth, each wave period is calculated when the code detects a zero-up crossing (pink circles). Then, the wave height is calculated by measuring the distance between the maximum and minimum values within each period (grey rectangles). Figure 13 shows the average wave heights for all cases at all facilities. Wave probe number 3 in the x-axis is the wave probe used in all the facilities (Figure 3).  At CNR-INM, wave probes 3 and 4 are in the same plane, but the value of probe 4 is 15% higher than that of probe 3. This difference could be associated with the different types of wave probes. Wave probe 3 is the capacitive wave probe used at all facilities, whereas wave probe 4 is ultrasonic. Wave probes 5 (dynamic) and 6 (ultrasonic) are also in the same plane, but their measurements are the same. The difference between the capacitive wave probe and the others could be an indication of lower accuracy or calibration errors of the probe. Despite these differences, the waves generated at CNR-INM are the ones closer to the nominal values.
At IFREMER, the average wave height seems to decrease as it gets closer to the rotor. However, for the bare-basin tests (IFREMER BB), the wave height increases closer to the rotor plane. This, for all cases except for case 5, showcasing the influence of the turbine in the upstream flow.
The data retrieved from FLOWAVE seems to have the most unpredictable values. Differences between wave probes are within 10%. However, the average values can be up to 50% (case 6) short of the preset values but as close as 15% (case 3) of the predefined values. Further work needs to be undertaken to characterise in detail the flow conditions, which is limited during this test campaign. It has been shown [37] that generating high amplitude waves at high frequencies when flow velocity is higher than 0.5 m/s is difficult at FLOWAVE.
Wave heights of all cases calculated with the zero-up crossing ( Figure 13) show similar values within 10% compared to the values using the Hilbert transform approach from Figure 11.

Turbine Performance
The results presented in this section are all the physical values representing the operation of the turbine. Figures 14 and 15 show the tests for the 0.8 m/s cases.
In general, the three facilities provide similar values of power and thrust for rotational speeds below 60 RPM. It seems that, at the three facilities, the maximum power happens at around 85 RPM.
Above the rotational speed at which maximum power is reached, CNR-INM provides values of power and thrust up to 20% lower than the other facilities. This behaviour is mostly related to the flow velocity data observed in Table 5, where the average flow velocity at CNR-INM is slower, yet closer to the nominal velocity, than at the other two facilities. The different blockage ratio between the tanks is shown to be another explanation for such a difference [15].
FLOWAVE provides similar power values for most of the rotational speed settings for the current-only case than at IFREMER. For the combined wave and current cases, however, values at IFREMER are predominantly higher. Compared to CNR-INM, values at FLOWAVE and IFREMER are up to 50% higher for case 1 and up to 30% for case 2.
Thrust values at FLOWAVE and IFREMER are higher than at CNR-INM by up to 15% for all rotational speed setting available.   In most cases, the repeatability of the experiment seems to be of good quality, except for some cases where the repeat points are notably different, highlighting the requirement to repeat data-sets in turbulent flows. This discrepancies are studied further in Section 3.5.

Uncertainty Analysis
It is not a common practice to include uncertainty analyses from experimental trials within the tidal energy community yet. In fact, not many studies indicate if experiments have been replicated in order to measure any precision errors. In [13,38], errors are reported in the region of 10% (average torque and thrust) and 3% (average power and thrust) for experiments undertaken in the IFREMER flume and the tow tank at the Kelvin-hydrodynamics lab (Strathclyde University), respectively. Both studies use the same equipment, but the investigations are unrelated; i.e., Murrat et al. [38] utilises this testing campaign to validate a numerical model and use tow-only conditions, whilst Ref. [13] explores the damping effects of bend-twist composite blades when operating under effects of waves and currents.
In the present research, average values of the uncertainty analysis done for the torque, thrust and rotational speed measurements for the three facilities with the help of Equations (1) and (2) are shown in Table 9. These values comprise the mean for all the test cases including regular and irregular wave-current conditions. Overall, these parameters show low uncertainty values in most of the experiments with average uncertainties as a percentage of mean values on the order of 2%. It is evident, however, that the data gathered from FLOWAVE indicate a larger uncertainty in the experiments (between 3 and 4%). The higher uncertainty values at FLOWAVE could be associated with the higher turbulence intensity (TI) values seen at FLOWAVE as observed in [15]. The results from this section are also analysed in terms of case numbers and rotational speed (at 0.8 m/s and 1.0 m/s). Tables 10 and 11 show the average uncertainty values from all the cases for torque and thrust, respectively. It is clear that the main contribution for a large uncertainty comes from the data gathered for case 4 and from the FLOWAVE facility, which is surprising, as this is a current-only condition and, therefore, the flow is less variable [15]. Looking into more detail, the largest uncertainty for case 4 (FLOWAVE) is related to peak power (105 RPM) and free-wheel conditions (184 RPM), where combined uncertainties as a percentage of mean values are higher than 14% and 7% for the torque and thrust values, respectively. However, only three cases are repeated for this experiment and also the uncertainty for the other case performed at 1 m/s (case 6) shows low average uncertainty values for both torque and thrust measurements (<2%). This highlights the importance of performing repeated tests, which, on occasion, could arise from some malfunctioning of the device or the flow generation. One relevant aspect that is not performed in this RRT and is advisable to other experimental tests is to repeat the flow characterisation at least at the beginning of the testing practice and, at the end, again to obtain any changes arising which may be due to the degradation of the current generation devices (or similar) over time. When looking at the analysis in terms of turbine operating points; i.e., rotational speed in this work (Figure 16), the IFREMER torque data suggest that, as the turbine approaches the free-wheeling condition, the uncertainty slightly increases. The same can be said for the thrust values except for the first two points at high load conditions. A similar trend may be depicted from the two CNR-INM cases, but, due to the limited number of tests, this hypothesis can not be confirmed, especially since the FLOWAVE data do not show a clear pattern. This behaviour may be related to larger vibration from the turbine when this increases its rotational speed.  This uncertainty analysis, although on some occasions incomplete, highlights the importance of repeating experiments, especially when these are performed in highly turbulent environments. Although flow characterisation repeat tests are not performed in these tests due to time constraints, they are recommended to further investigate uncertainty values which may be related to the facility rather than the device itself.
The tests performed in environments with zero turbulence levels (CNR-INM) do not reflect the real operating environment of a tidal turbine. Nevertheless, the comparison with results from facilities where turbulent flows can be generated is very useful to separate and analyse the effects of the testing environment on turbine performance characterisations.

Turbine Stabilisation
In an attempt to understand the high uncertainty values seen at FLOWAVE (Tables 10  and 11), Figures 17 and 18 look at how the turbine loading stabilises with time when operating at TSR = 4 for cases 1 and 2 for thrust and torque, respectively. Images show repeat runs when available. These are calculated using the root-mean-square error (RMSE) over the time series (see Section 2.7). At FLOWAVE, the thrust RMSE stabilises between 6 and 8 N, whereas, at IFREMER, the thrust settles at values around 2 N except for the first few seconds of one case. At FLOWAVE, it takes the torque and thrust up to 50 s to stabilise, whereas this is obtained at around 20 s in IFREMER. From the present observations and previous observations done by [15], the turbulence intensity, the shear flow and the vertical flow component are higher at FLOWAVE than at IFREMER for case 1; thus, in the future, it would be advisable to acquire longer samples to allow the turbine to stabilise better in faster flow conditions. The other case explored is wave case 2. In contrast with the current-only cases, the time series seem to stabilise sooner at FLOWAVE than at IFREMER. Again, relating this information to the flow characteristics, the turbulence intensity and wave heights obtained at IFREMER for this specific case are greater than at FLOWAVE. The addition of waves is responsible for the larger fluctuations around the average values seen at the beginning of the tests (∼10 s).
To understand the loading stabilisation, Figure 19 looks at the convergence of the flow velocity u for case 1 at all facilities using the ADV mounted next to the rotor. At FLOWAVE and IFREMER, the turbine rotation has an influence on the flow velocity RMSE. As the rotational frequency of the turbine increases, the RMSE value decreases. In addition, RMSE values do not seem to stabilise throughout the length of the tests at FLOWAVE. At the lowest RPM, RMSE at FLOWAVE is more than three times the overall values at CNR-INM and IFREMER.

Conclusions, Recommendations and Perspectives
This paper presents the results obtained from the first three trials done for the tidal MaRINET 2 Round Robin Test at the IFREMER, CNR-INM and FLOWAVE facilities. This Round Robin Test extends the previous program, i.e., keeping constant the turbine model and testing instruments, by incorporating waves and currents in the testing protocols.
As observed in the inflow velocity analysis, the hydrodynamic conditions are close between the three facilities. However, several differences appear, especially in terms of standard-deviation, mainly due to turbulence and wave-current interaction at the flume tanks. These differences are related to the intrinsic characteristics of the wave and current generation at all the facilities.
The synchronous measurements obtained with the ADV on the turbine side are influenced by the turbine's rotation. Flow velocities recorded with the ADV can be 8% higher than the upstream LDV measurements. This phenomenon is related to the acceleration of the flow as it bypasses the turbine. Consequently, the ADV velocity measurements on the side of the turbine cannot be used for normalising the performance coefficients or as reference values.
A number of wave probes are utilised during this Round Robin tests, highlighting the quality of different wave probe types to assess surface elevation, with the ultrasound and dynamic wave gauges showing a better performance. However, this wave gauges are only used at the tow tank facility. Work related to the use of these type of probes in turbulent flumes will be assessed in the future.
Depending on the facility, a larger amount of repeatability tests and longer acquisition times are advisable in order to reduce the uncertainty of the outcomes from the turbine loading and performance and improve the convergence of the time series. This is especially important when the data are used to verify numerical models. It is advisable to repeat flow characterisation procedures which may be influenced by the flow generation systems; e.g., overheating of pumps.
Due to the finite length of the towing tank at CNR-INM, the combinations of carriage speeds and wave periods allowed for an average of 20 wave periods per run. The achievable wave periods, however, have a high repeatability between them. Thus, deeming the quality of the data comparable to the data obtained at the other facilities.
Despite the work presented here, a deeper analysis is required. In the future, processing the blade forces and moments is envisioned. In addition, further tests are planned at the Kelvin Hydrodynamic Laboratory (KHL) towing tank at the University of Strathclyde and at the Queen's University Marine Laboratory (QML) tidal site. These coming RRTs will bring new wave-current flow conditions and certainly new interesting results.