4.1. Site Layout and Characterisation
To assess whether there are any features of the seabed that could affect device performance, industry guidelines [
21] recommend surveying the bathymetry at the deployment site out to 5 equivalent rotor diameters (D) either side of the turbine, and 10 D upstream and downstream, covering an area of 10 D × 20 D. This region is shown in
Figure 6, with the area offset by 8
to align with the dominant flow directions (see
Figure 2). The turbine sits on the northernmost apex of its triangular frame with the rotor facing west when parked and non-operational. The yaw system range (±100
) was sufficient to allow the rotor to face the dominant directions of both tides, and accounted for the small offset frame installation angle from true north. The seabed ADP was housed in a small structure referred to as the Remote Acoustic Monitoring Platform (RAMP), and was cabled to the turbine in order to receive power from the shore and to enable live data transfer.
The turbine sits on a relatively flat ridge to the east of the northern portion of the trench that runs through Ramsey Sound. The depth extremities within the 10 D × 20 D area range from 31 m at just over 5 D to the north of the turbine, to 44 m near the north-west apex of this considered area. Compared with the seabed features found elsewhere in Ramsey Sound (see
Figure 1), this local depth variation is not considered significant. Both the RAMP and turbine sit at a mean depth of 35 m.
The flow characteristics of two full tidal cycles at the site are presented in
Figure 7, measured with the seabed mounted ADP using the hub height bin only (
z = 12.1 m), i.e., at a relative submergence of
= 0.35. The measurements are of 10 min time-averaged velocity (
) and turbulence intensity (
, with
denoting mean velocity and
is the root-mean-square of the velocity fluctuations during the
i-th bin). Shaded grey areas in
Figure 7 outline periods of slack water with flow velocities below a velocity threshold of 0.5 m·s
. There was an important variation of velocities between tidal phases evidencing the asymmetric flow pattern at the site. During flood tides, the flow was characterised by large velocities reaching values over 2.5 m·s
, while ebb tides attained peak velocities up to 1.45 m·s
. During slack water periods, turbulence levels greatly increased with values over
0.4, but these did not challenge the structural integrity of the turbine as the flow velocities were low and the device did not operate under such conditions.
Approach flow velocities across the measured water column during both ebb and flood tidal phases exhibited a time-dependent pattern, both in its vertical distribution and magnitude. This behaviour is observed in
Figure 8 which presents the values of
and
measured with the seabed ADP at depths between 3.1 m
30.1 m at eight different time intervals, covering both ebb and flood tides depicted in
Figure 7. At the onset of the turbine operation during ebb tide, i.e., at
, the flow overcame 1 m·s
velocity, and ebb tide peak velocities of approximately 1.5 m·s
were reached at
. The velocity distribution at the falling ebb tide, specifically at
, attained a very similar distribution to that at
, although with a slight difference found below hub height.
Figure 8b shows that during the flood phase, velocities in the rising flood profile (
) seem to follow a logarithmic distribution until hub height, but higher up in the water column the velocities do not increase much and are nearly constant at a value of 2 m·s
. At peak flood tide, i.e.,
, the maximum velocities were observed at elevations higher than the hub height with a larger region of the water column following a logarithmic distribution than at
. In the transition from flood flow to slack water, at
, the vertical distribution of velocities was very similar to
, although with a lower maximum value of 1.7 m·s
. From
Figure 8a,b it can be observed that the velocity distributions during ebb flow can be approximated to a logarithmic profile, whereas in the flood tide the velocities do not always follow such behaviour over the entire water depth, as found by others [
2,
8].
Turbulence intensity profiles, shown in
Figure 8c,d, depict again the notably asymmetric flow conditions during ebb and flood tides and variability within each tidal phase. At peak ebb flow, maximum values of
were found near the seabed and free-surface, attributed to bathymetry-induced and wave effects respectively, and minimum values were attained just above the turbine top tip, i.e.,
20–22 m. Conversely, during flood tide the lowest turbulence intensity levels were found close to the hub height whereas maximum values were again observed near the lowest and highest locations in the water column.
The variation in the spatial characteristics of turbulence suggests that the turbine’s blades face anisotropic turbulent flow conditions during its operation. The time-dependent heterogeneous distribution of
between ebb and flood tides was also found at other sites such as the Sound of Islay (UK) where the turbulence intensities can vary between 0.11 and 0.13 during flood and ebb tides respectively, or at the East river (US), which is slightly more turbulent attaining values of
equal to 0.13 (flood) and 0.18 (ebb) [
22].
Further analysis of the spatial variability of the velocities across the water column is performed in
Figure 9 with the Power Spectral Density (PSD) distribution of 10 min bin velocities during peak flood and ebb tides at hub height (
) and at bottom and top tips, i.e.,
. The
function was used to compute the spectra, dividing the signal into 20 ensembles. The spectra suggest that energy decay during the production range followed a −1 slope commonly found in boundary-layer flows [
23,
24], and the spectral energy (
) associated with the lowest frequencies was greater during flood tide, indicating that the largest flow scales approaching the turbine during this phase were more energetic than during ebb. This is expected as the flow velocities during flood were greater than in ebb (see
Figure 8).
The spectrum at
z =
during flood tide presents a lower spectral energy than those computed at higher elevations, identifying the spatial variation of the flow dynamics within the turbine’s swept area. Such flow variability at the deployment site, due to the irregular bottom bed shape (
Figure 6) among other local features, are well-identified from both the vertical profiles of time-averaged velocity statistics and power spectral energy distribution, shown in
Figure 8 and
Figure 9 respectively. Previous studies acknowledging the tidal flow asymmetry of Ramsey Sound were reported in [
1,
8,
25].
Distribution of flow velocities upstream of the turbine in the range between 1.4 m to 20.4 m measured with the turbine ADP are shown in
Figure 10, with the mean velocities and turbulence intensities calculated over 10 min bins during peak ebb and flood, i.e.,
e2 and
f2. The inflow velocities are normalised by the mean of the values obtained at an upstream range,
x, greater than one rotor diameter (12 m), which are the outermost eight ADP bins. This definition is considered the undisturbed free-stream velocity,
, since the measurements from this range are generally within 2% of each other. The velocities were subject to a considerable reduction near the rotor due to the flow moving around the turbine. This behaviour has been observed by others using a similar turbine mounted ADP arrangement [
26].
Furthermore, velocity measurements in
Figure 10 are compared with vortex sheet theory [
27] based on different thrust coefficients of the turbine, represented in terms of the induction factor (
a), which determines the upstream velocity as:
A reasonable agreement is found between field velocities and vortex sheet theory during both tidal phases at upstream distances of
6 m irrespective of the thrust coefficient attained by the turbine, which signifies Equation (7) can be used to approximate the tidal flow velocities that approach the turbine, analogous to those obtained in wind turbine research [
27,
28]. The reduction on the approach velocity was consistently found to be greater during flood tides, believed to be a consequence of the flow being disturbed by the triangular frame supporting the turbine (
Figure 3), which contributes to a larger flow deceleration at
5 m than that predicted by the vortex sheet theory. The frame only sits upstream of the rotor during flood tides (
Figure 6), with its front edge approximately 10 m upstream of the rotor. This was consistent with the location of a small trough in the approaching flood velocity and a peak in turbulence intensity.
The turbulence intensity near the rotor was also found to be consistently higher during flood, whereas further upstream (x > 15 m) the levels are not dissimilar to ebb (12–14%). These features were observed throughout the flood measurements, and provide further evidence of the turbine frame disturbing the approaching flow. Since these observations have been made on the turbine ADP which sat at an elevation of almost 8 m above the frame height, it should be expected that the flow disturbances were greater at lower elevations, i.e., in the lowermost regions of the blade rotations.
4.3. Hydrodynamic Loading Coefficients
In addition to the rotor thrust coefficient (Equation (3)), the following dimensionless parameters are used to describe the blade root bending moments,
, and axial forces,
Dimensionless rotor loading characteristics measured during ebb and flood tides are presented in
Figure 12, compared with numerically predicted curves from a Blade Element Momentum-based model described in Harrold et al. [
16]. These simulations were performed under steady conditions, with the turbine running at a fixed speed in a uniform flow field. The measurements represent the mean values from within the rotor speed stepping periods after discarding measurements from within 10 s of each step change in order to allow some settling time.
Considering firstly the results obtained during ebb conditions, the measured
values show a very good agreement with the theoretical prediction for each of the three blades, scattering closely about the numerical curve. Similarly, the measured
values on blade 1 are in good agreement with the theoretical curve and the predicted trend of decreasing
with increasing
is clear in the results from all of the blades. However, the measured
values show greater blade-to-blade variation and their magnitudes are generally lower than the theoretical prediction. This can be attributed to the derivation of
as described in
Section 3.3, in which this load was entirely dependent on the shape of the curve fitted to the measured radial bending moments. While some aforementioned measures were introduced to improve the robustness of this process, it is clear that further refinement is required to improve the derivation of
, particularly on blades 2 and 3. In contrast to this, the derivation of the
values only required a simple extrapolation and, as a result, a better agreement was found, both in terms of the blade-to-blade results and with the numerical model prediction. The difficulties encountered with the
measurements also led to lower than expected rotor thrust loading, although the shape of the
curve was captured well.
To improve the measurements of the rotor loads, ideally a secondary calibration and the relocation of the problematic sensors on the blades at radial position 1.7 m would have been performed, but this would have required the physical removal of the turbine from the site. The existing procedure could be improved by further modification of the applied weighting to measurements, such that the agreement between loads calculated on the compression and tension surfaces of the blades improves. Similarly, any modifications should attempt to improve the agreement between blade-to-blade averaged over a long enough period, i.e., a sufficient time interval to reduce any bias being introduced from the loads changing on the blades during rotation.
As for the flood tide results, fewer measurements were obtained during these conditions and only within the tip speed ratio range 1–3, limiting the data available to compare with those obtained in ebb. However, in the region that does overlap with the ebb measurements, i.e., between
= 2–3, the mean loads were consistently lower in flood conditions. For example, the
measurements mostly cluster around a value of 0.065 at
= 3 during flood, whereas in ebb they were found at values close to 0.087. It is believed that this is directly attributable to the reduced inflow and increased turbulence found during flood conditions, as observed in
Figure 10. Since the loading coefficients and
were calculated relative to the undisturbed upstream flow at hub height, the reduced inflow observed near the rotor during flood was not accounted for, nor was the slowing of velocities higher up in the water column (
Figure 8). Therefore, a reduction in mean rotor loading was experienced during flood due to the influence that the turbine frame structure and site bathymetry have on the inflow.
4.4. Instantaneous Structural Loadings
Instantaneous variation of structural loadings as a result of the fluctuating approach velocity is presented in
Figure 13 with velocities measured with the turbine ADP from bins located at distances
x = 9.4 m and 20.4 m, i.e., 0.8 D and 1.7 D, upstream of the turbine, and measurements of axial force and bending moment at the root of the three blades. These data correspond to a 2-min snapshot of the 5-min stepping period in which the turbine operated at a fixed rotational speed of
= 11.61 RPM during the ebb tide presented in
Figure 11. Velocity measurements at two different distances from the turbine are presented to observe the local flow variation within the short distance of 11 m before impinging the turbine. All the blades undergo similar magnitudes of structural loadings, although the peaks at the same rotation interval did not attain the same value, i.e., blades experienced different extreme loads during their operation as the flow is markedly uneven across the rotor swept area. Similar results of uneven load distributions among blades were found in Ouro and Stoesser [
24] who performed high-fidelity simulations of a turbine operating over irregular bathymetry.
The distribution of spectral energy of
and
across relevant frequencies is presented in
Figure 14, with the spectra computed from the loadings on one of the blades during the 5 min interval shown in
Figure 13, applying an analogous procedure as for the velocity spectra from
Figure 9. Frequencies are normalised by the turbine rotational frequency
=
(RPM)/60 s)
Hz. At frequencies in the range of
there is a peak in spectral energy in both axial load and bending moment associated with low-frequency velocity oscillations resulting from the action of large-scale flow structures. Note that this frequency coincides with the transition from the production range to the inertial sub-range obtained from the velocity spectra in
Figure 9.
An energy peak is observed at a frequency of
derived from the periodic change in the blade loadings during every rotation, which is only observed in the bending moment spectrum. The spectral energy linked at a frequency of
has a similar magnitude to peaks found at higher frequencies, despite previous data from laboratory tests showing that the maximum energy may be attained at the the turbine rotational frequency [
11]. These results can be linked to the relatively short sampling period of 5 min and due to high turbulent flow conditions in which the full-scale device operated unlike the fairly controlled conditions that are normally found during laboratory testing. Two peaks in the
and
spectra are found at frequencies of 1.35
and 1.70
owed to blade-rotation effects, and only in the latter spectrum another peak is observed at 2.0
. Finally, a large energy peak whose frequency is almost that of the data acquisition and harmonic of the rotation frequency, i.e., 16 Hz, is clearly observed from the spectra. Similar spectral energy distribution of the loadings on a prototype turbine in the frequency domain were found in Payne et al. [
11] who also discerned that peaks can vanish when the turbulence intensity in the approach flow is high.