1. Introduction
Marine tidal currents stand out as a promising and renewable energy resource with the advantages of high predictability and consistent availability. This reliability makes them ideal for optimised energy output [
1]. To effectively harness this energy resource, the installation of turbines in arrays proves to be the most efficient approach, to ensure maximum power extraction [
2]. In recent years, a diverse range of different turbine designs has emerged [
3]. Despite the many proposed designs of tidal stream turbines, the configuration that has garnered the most attention, especially in the context of commercial-scale array development, is the three-bladed horizontal axis tidal stream turbine [
4]. Numerous tidal stream turbine (TST) arrays have been strategically deployed in highly energetic locations worldwide, in most cases serving as generators to produce electricity for local networks [
5].
Employing a TST array stands as an effective strategy for maximising the power extraction from a given tidal energy site, while also contributing to cost reduction by sharing transmission and grid integration systems when compared to individual TST installations [
6]. However, the deployability of TSTs in an array is limited by factors such as water depth and bathymetry. Furthermore, the wake generated by TSTs can impact the downstream velocity distribution, introducing turbulence that alters the flow field, and thus affecting energy output. Consequently, investigating the spanwise and streamwise spacing between turbines is important to minimise disturbance induced by TST and maximise energy output [
7]. Several experimental and numerical studies have been conducted to improve the understanding of the TST, including aspects like design, environmental impact, and hydrodynamic performance of the TSTs, both individually and in arrays. Numerous studies [
8,
9,
10,
11,
12] have presented critical experiments suitable for the assessment of turbine performance. The data sets derived from these experiments serve as valuable benchmarks for the validation of numerical models that characterise turbine performance, providing details on local turbine flow details to be determined, which are difficult to obtain using experiments only.
Various approaches exist to numerically model TSTs to predict the response of tidal turbines in response to complex environmental inflow conditions, including their loading, performance, and wake generation. Computational fluid dynamics (CFD) is one such approach used to model TSTs. There are advantages and disadvantages to using CFD, the critical factor lies in striking a balance between detailed simulation of the physics and the computational time and resources required for accuracy [
13]. Among the CFD models, the most comprehensive and intricately detailed approach for gaining insights into the wake development of a TST is the fully resolved turbine geometry approach [
14,
15]. However, this approach requires small time steps to solve the flow field around the turbine, resulting in a high computational demand. For this reason, the fully resolved turbine geometry approach is not always a feasible option when modelling a large array of full-scale turbines to provide considerable predictions of turbine performance and wake flow, especially in environments with a wide range of highly variable flow conditions. Thus, an alternative method is needed to represent the TST in CFD simulations.
Two popular modelling approaches to represent a TST are the actuator disk (AD) and blade element momentum (BEM) theory methods. The coupling of the AD and BEM methods, known as the AD–BEM model for TST array studies, has gained traction in recent years [
4,
16,
17,
18,
19,
20]. The AD is a momentum sources model that describes the force distributions along the rotor blades to determine the overall performance of a turbine with significantly lower computational costs and fast processing time. However, the AD model is typically unable to provide realistic near-field downstream wake results within seven diameters [
21]. Consequently, BEM theory has been introduced to the AD model to better describe the momentum source term by imitating the effects of blade twist, chord, lift, and drag characteristics. The BEM–AD model requires a reasonable but inexpensive increase in computational requirements.
There were some limitations to the BEM–AD models presented in this paper. Firstly, the model is non-rotating, therefore eliminating any swirl in the flow, but studies have suggested that the effects of swirl are not significant outside the near-field downstream wake region of the TST [
22]. However, in some situations, swirl in the near wake can persist further downstream, potentially influencing the flow boundaries and distorting the wake [
23]. Secondly, it is assumed that tip vortices from a rotating turbine blade are ignored due to the actuator disk’s inability to replicate these vortices, thereby speeding up the wake velocity recovery rate downstream [
24,
25]. Thirdly, the BEM–AD model was treated as a steady-state model, keeping it in line with the RANS time-averaged nature. However, a pseudo-transient-state model can be produced to account for the various turbulence parameters [
13]. Other approaches have been used to model a TST such as moving reference frame (MRF) and sliding mesh. MRF and sliding mesh approaches can capture the rotational effects of TST such as swirl and blade tip vortices [
26] but require very fine grid resolutions around the turbine rotor to capture flow field detail; hence, this is computationally expensive [
27]. In recent years, there have been various important studies that focus on the development of vortex flow [
28,
29,
30]; these have brought MRF and sliding mesh back into the spotlight in array modelling.
Numerous studies have been carried out on TSTs using the RANS-actuator disk model to predict the velocities and turbulence intensities in the wake [
21]. Studies conducted by Harrison M.E., et al. [
31] and V.T. Nguyen [
32] have highlighted the performance of Reynolds-Averaged Navier–Stokes (RANS) actuator disk models in 3-dimensional domains. Similar studies were also carried out in an array study performed by Mycek. P [
11], Chawdhary S. [
33], and Gaurier et al. [
2] to investigate turbine–turbine wake interaction in an array to determine the optimal spacing both spanwise and streamwise. Turbine spacing is one of the most important factors influencing the performance of a tidal array, but its effects on the turbine wake interactions or turbine array efficiency are not yet well understood. Hence, an improved understanding and description of individual turbine wake development could greatly aid in optimising turbine layout or spacing that provides the maximum array efficiency [
34].
The work presented in this paper provides details on the investigation carried out on the developed BEM–AD model’s ability to predict downstream wake effects for potential future use in an array study. The paper presents an empirical formula approach to better describe the relationship between velocity profile and geometry radial effects on porosity and resistance coefficient of a porous disk, as showcased in the work as two modified variations of the BEM–AD model and a hybrid modified BEM–AD model. All models in this work are validated with experimental measurements by L.E. Myers (2013) [
8,
22]. The first variation is a modified BEM–AD model to incorporate the influence of inlet velocity profile (i.e., tidal current shear) in the disk configuration, which is known as velocity variation. The second variation, known as the radial variation, modified the BEM–AD model by dividing the actuator disk into radial annular elements; this allows the disk to have a radial variation instead of an overall averaged approach. The radial variation also included corrections such as tip loss correction, etc. to the BEM calculation. Meanwhile, the hybrid modified BEM–AD model is a combination of two of the variations discussed. In the sections below, the theory, methodology, results, and discussions of the modified BEM–AD models are detailed, and finally, some key conclusions are presented.
6. Results
In this section, the results obtained from the BEM–AD model are presented and discussed. This section contains four sub-sections detailing the capability and suitability of the two variations and hybrid models to predict downstream wake effects in the fluid flow field. The first sub-section presents the velocity variation BEM–AD model, which takes into account the velocity profile of the water column, and the results are compared to experimental measurements [
22]. The second sub-section presents the radial variation BEM–AD model, which incorporates the blade element momentum theory method into the disk, and the results are compared to experimental measurements. The third sub-section is the hybrid modified BEM–AD model, which combines both radial and velocity variations to further improve the accuracy of the results. Lastly, the final sub-section provides an overall discussion.
6.1. Modified BEM–AD Model: Velocity Variation
In this sub-section, the developed velocity variation BEM–AD model was validated against experimental measurements. The disk’s porosity and resistance coefficient take into account the velocity profile of the water column. Details of this configuration are mentioned in
Section 3.2 and
Section 4.4.
Figure 14 shows the centreline downstream velocity and turbulence intensity of the velocity variation BEM–AD model against experimental measurements. The developed model matches well with the experimental measurements, especially after a distance greater than 8D. The model tends to overpredict both velocity and turbulence intensity at downstream distances of less than 7D.
Figure 15 shows the downstream velocity and turbulence intensity profile at downstream distances of 5D, 8D, and 10D. The velocity profile at 10D matches closely with the experimental measurements, but for velocity, profile overpredictions occur at 5D and 8D. Whereas, the turbulence intensity is underpredicted at 5D, 8D, and 10D and underpredicted at the centreline region at 5D. The results show that the velocity variation model can predict downstream velocity well at distances greater than 8D. However, it is observed that for vertical depth above the centreline, the model is less accurate compared to below the centreline. The vertical depth below the centreline has a more severe change in vertical velocity profile compared to the vertical depth above the centreline, see
Section 3.2.1 for details. This results in the disk having a larger variation in porosity and resistance coefficient below the disk centreline compared to above, which causes the model to describe the downstream wake below the centreline better than above the centreline.
6.2. Modified BEM–AD Model: Radial Variation
In this sub-section, the developed radial variation BEM–AD model was validated against experimental measurements. The blade element momentum theory method is incorporated into defining the disk’s porosity and resistance coefficient. Details of this configuration are described in
Section 3.2 and
Section 4.4.
Figure 16 shows the centreline downstream velocity and turbulence intensity of the radial variation BEM–AD model with experimental measurements. The developed model matches well with the experimental measurements, and it seems to be performing better at predicting the centreline downstream wake than the velocity variation BEM–AD model, especially at a downstream distance of less than 8D.
Figure 17 shows the downstream velocity and turbulence intensity profile at downstream distances of 5D, 8D, and 10D for the radial variation BEM–AD model. The velocity profile at 10D matches well with the experimental measurements in both downstream velocity and turbulence intensity. It is noticed that downstream velocity matches closely with the experimental measurements around the depth region of 1 to −1, while the model outside this region tended to underpredict the results, and the downstream turbulence intensity was shown to be overpredicting around the depth region of 1 to −1. The high accuracy in predicting the downstream velocity at the depth region of 1 to −1 is contributed by the radial variation approach describing the porosity and resistance coefficient of the disk, while it is speculated that outside the depth region of 1 to −1, the wake prediction being less accurate might be due to the velocity profile not being taken into account while describing the disk properties.
6.3. Modified BEM–AD Model: Hybrid Model
In this sub-section, the developed velocity variation and radial variation BEM–AD models were combined to form a hybrid modified BEM–AD model in which both the blade element momentum theory radial approach and the velocity profile in describing the disk properties are taken into account; details of this configuration are described in
Section 3.2 and
Section 4.4.
Figure 18 shows the centreline downstream velocity and turbulence intensity of the hybrid modified BEM–AD model against experimental measurements. The developed model matches closely with the experimental measurements; it is shown to be performing better than both the velocity and radial variations. The combination of both variations greatly improves the ability of the hybrid model to predict downstream wake; this can be further observed in the vertical profile of downstream distance 5D, 8D, and 10D in
Figure 19.
Figure 19 shows the downstream velocity and turbulence intensity profile at a downstream distance of 5D, 8D, and 10D for the hybrid modified BEM–AD model. The velocity profile at 10D matches well with the experimental measurements in both downstream velocity and turbulence intensity. It is observed that the overall wake results greatly improved, especially for all results below the vertical depth of 1. However, the results above the vertical depth of 1 show underprediction in both downstream velocity and turbulence intensity. The reason for this might be due to the developed model not taking into account fluid surface conditions and the surface being treated as a symmetry boundary face. A relationship was observed between velocity and turbulence intensity, i.e., a higher downstream turbulence intensity will result in a higher downstream velocity. This shows that a high turbulence intensity promotes wake velocity recovery.
6.4. Overall Discussion
Overall, the modified hybrid BEM–AD model has proven to be the most accurate among all of the models proposed and investigated in this work. All three models display a similar degree of accuracy after a downstream distance greater than 8D. At the downstream distance greater than 5D, the ability of the radial variation model and modified hybrid model to capture the wake effects shows little difference. However, for a downstream distance of less than 5D, the modified hybrid BEM–AD model provides the most accurate results.
Figure 20 shows the centreline downstream velocity and turbulence intensity of the three models against experimental measurements. In
Table 5, a statistical analysis is presented for the three studied models which are compared with experimental measurements of centreline downstream velocity and turbulence intensity.
The statistical analysis presented shows that the velocity variation model has the worst correlation and the largest error among all three models, while the radial variation model and the hybrid modified model show very similar correlation and error when compared with the experimental measurements. The hybrid modified model correlates best with the experimental measurements in terms of both downstream velocity and turbulence intensity with a coefficient of determination, R2 value of 0.9917 and 0.9863, respectively. The hybrid modified model has the smallest error among all three models, with root-mean-square error values of 0.0131 and 0.0058 and mean absolute percentage error (MAPE) values of 1.31% and 2.82% for velocity and turbulence intensity, respectively.
Figure 21 and
Figure 22 show contours of velocity and turbulence intensity comparison between radial variation and hybrid modified models. Radial variation models have a more intense turbulence intensity than the hybrid modified model; this will decrease the rate of wake velocity recovery, and it is reflected in the velocity contour. Furthermore, the hybrid modified model shows a shorter contour shape than the radial variation model.
Figure 23 shows the power density of the disk for the three studied models. It is observed that the velocity variation disk has a high-power density near the edge of the disk when compared to the other two models. This is due to the velocity variation disk not taking radial blade element characteristics into account. Thus, some BEM-calculated parameters were neglected, such as tip loss correction, while the radial variation model and the hybrid modified model were able to address the tip loss issues.
Table 6 is a comparative analysis of the overall thrust and power coefficient of the disk in the three studied models compared with a BEM numerical model discussed in this work.
Table 6 shows that all three models have a lower thrust coefficient prediction than the numerical value, which were 0.7571, 0.6278, and 0.7919, respectively; and the velocity variation and hybrid modification models have a higher power coefficient prediction than the numerical value, which were 0.5229 and 0.4714, respectively; while the radial variation model has a lower power coefficient prediction than the numerical value, which was 0.4616.
The radial variation has the highest difference with the numerical value in terms of thrust coefficient, which was 21.13%, and the velocity variation has the highest difference with the numerical value in terms of power coefficient, which was 11.73%. These indicate that the velocity variation model which incorporates the inflow velocity profile in describing the disk has predicted the thrust more accurately than the radial variation model. The radial variation model, which incorporates radial blade element characteristics in describing the disk, has predicted power better than the velocity variation model. Overall, the hybrid modification model has the least difference when compared with the numerical thrust and power coefficient, which were 0.51% and 0.72%, respectively. Hence, the hybrid modification model has the advantage of both velocity and radial variation due to the incorporation of radial blade element characteristics and inflow velocity profile in describing the disk. However, the hybrid modification model requires more setup and takes longer to solve than velocity and radial variation. Hence, the velocity variation can be used in cases which require good thrust coefficient calculations, while the radial variation can be used in cases which require good power coefficient calculations.
Figure 24 presents the side-view power density contour of the hybrid modified model. The contour in
Figure 24 gives a good illustration of the available power density downstream from the disk, as observed as the wake develops further downstream from the disk, the available power from the wake also recovered.
Table 7 shows the power coefficient performance of a second disk at downstream distances of 5D, 8D, 10D, and 15D at different heights. It is observed that at a disk vertical offset distance of +1.0 D, the second disk was mostly unaffected by the downstream wake of the first disk and has a difference of less than 2.0% compared to the performance of the first disk. Whereas, it is seen from
Table 7 that a positive vertical offset of the second disk has a better performance than a negative vertical offset. At the vertical offset distance of –1.0D, the second disk was least affected by the downstream wake of the first disk but was influenced by the bottom surface and has a difference of less than 27.1% compared to the performance of the first disk.
7. Conclusions
The work presented in this paper demonstrates the ability of the modified actuator disk models to predict the wake effects of a tidal stream turbine. The hybrid modified model is the most accurate approach to predict the time-averaged velocities and turbulence intensities in the wake of a turbine; however, the pre-solver configuration of the hybrid modified model is more complex compared to the velocity and radial variation model. The results between the three models show very little difference in wake velocity at the downstream distance greater than 5D. For investigations of the far-wake region (greater than 5D), the velocity and radial variation model is a more computational efficient approach. For situations where the velocity profile is constant throughout the whole depth, a radial variation approach can be adopted, while for situations where the velocity profile is varying, a velocity variation approach can be adopted. On the other hand, the hybrid modified model is suitable for investigations for the near-wake region (less than 5D).
These studies also reveal a strong relationship between turbulence intensity and wake velocity recovery rate. A less-intense turbulence after the disk will promote the wake velocity to recover rapidly, and the intensity of the turbulence after the disk is hugely influenced by the porosity and loss coefficient of the disk domain. As a result of these findings, a much more in-depth investigation is needed to investigate the relationship of porosity and loss coefficient of the disk domain with the downstream turbulence intensity. At an offset vertical distance greater than 1D, placement of a second turbine anywhere downstream would yield little difference, while for an offset vertical distance between 0.5D and 1D, a downstream distance greater than 8D would be a better placement for a second turbine. This requires further investigation such as the effects of sitting multiple turbines in the flow field, the effect of yaw misalignment, and also predicting the array layout and power output of a real tidal current energy site. Other areas of interest for further work include an investigation to improve the BEM numerical model to address some of the underprediction issues on thrust coefficient when the tip speed ratio is greater than 8. Also, an investigation into the impact of different turbulence models in predicting the downstream wake development would be valuable.