Numerical analysis of blockage correction methods for tidal turbines

The current work presents the analysis of six different blockage correction methods for axial-flow turbines from the literature and compares them to a method proposed by the authors using blade-resolved CFD simulations. The analysis of the different methods is performed by comparing the corrected blocked-flow thrust and power results of a tidal turbine to a quasi free-flow condition, obtaining a quantitative measure of the error for every method at different blockages and tip-speed ratios. This study shows that not every available method for blockage correction remains valid for the range of turbine thrusts that may occur in blocked flows. Of the seven methods analysed, only three were found to provide a reliable estimate of the power and thrust across the range of blockage and tip-speed ratios analysed in this paper, including the one proposed by the authors. Of the remaining correction methods, those that could be used as a first approximation for low-blockage conditions are identified.


List of symbols
Velocity at the turbine plane a Axial induction factor β Blockage ratio C Channel cross-sectional area P +/− 1 Static pressure at both sides of the turbine plane according to the actuator disc theory U T Tank or channel undisturbed flow velocity U F Equivalent free-stream velocity A R Cross-sectional area of the wake in a reference plane A W Cross-sectional area of the wake where the static pressure equalizes with the bypass δ F Empirical factor for the proposed blockage correction method U 2 Velocity in the wake where static pressure equalizes with the bypass U 3 Velocity in the bypass of the turbine plane U 4 Velocity in the bypass where static pressure equalizes with the wake E C X Normalised error for thrust and power MAE X Average magnitude of the error values γ Blade twist angle

Introduction
Blockage is the interaction effect between the flow past a body and the surrounding boundaries in constrained flow conditions. Its effects have been discussed since the early application of laboratories to study aero-and hydrodynamic phenomena (Glauert 1933). Wind tunnels, towing tanks and other fluid-dynamics laboratories are prone to induce interaction effects not seen in full-scale open flows due to the close proximity of the walls to the body. Blockage effects are also present for full-scale hydrokinetic turbines where the swept area of the turbines can represent a substantial fraction of the tidal channel cross-sectional area. Blockage effects modify the thrust and power performance of the turbines. Blockage interaction effects are typically observed, in the particular case of turbines, as an increase in both power and thrust forces when compared to the open-flow condition. This is caused by different effects such as an increase in flow speed around the body, a change of pressure in the wake, differences in flow-development for lifting surfaces such as foil sections or even longitudinal pressure gradients associated with the tank boundary layer and fluid losses (Glauert 1933;Pope and Harper 1966).
Early attempts to explain the behaviour of a turbine were made by means of the one-dimensional actuator-disc theory. This concept was employed to derive the energy extraction limit for an open-flow condition, set by a maximum power coefficient C Pmax = 16/27 known as the Lanchester-Betz limit (Lanchester 1915;Betz 1920). These efforts were expanded by Garrett and Cummins (2007) showing that the maximum power coefficient for an actuator disk in a constrained flow is C Pmax = 16/27(1 − β) −2 , where β is the blockage expressed as a ratio between the turbine swept area and the channel cross-sectional area, showing that an increase in power extraction is dependent on the blockage ratio.
A more complex approach to the problem of modelling a turbine has been done by considering conservation of linear and angular momentum in a series of concentric annuli over the rotor swept area. The axial and tangential forces of the turbine on the flow are determined by coupling momentum theory with blade element theory by applying bi-dimensional lift and drag curves as functions of angle of attack for the turbine aerofoils. Overall turbine torque and thrust are determined by integrating these forces across the blade span. This approach is known as the Blade Element Momentum method (BEM) (Glauert 1947), and has been widely employed for turbine and propeller design and load assessment (e.g. Hansen et al. 2006;Wang et al. 2016). The traditional derivation does not consider blockage and is limited to steady two-dimensional analysis of the flow over the aerofoils and thus extensions have been proposed to consider effects such as tip-losses (see, e.g. Burton et al. 2001;Wimshurst and Willden 2017), turbulent wake effects (e.g. Buhl 2005), and even transient phenomena (e.g. Hansen et al. 2004). The effects of blockage have recently been incorporated into BEM theory by Vogel et al. (2018) by modifying the linear momentum equations. The modified BEM theory shows reasonable agreement to within 5% of comparable blade resolved simulations.
Turbine modelling can also be carried out by means of computational fluid dynamics (CFD) simulations. These methods consist in obtaining numerical solutions to the Navier-Stokes equations, with different models of the physics and resolution, both in spatial and temporal dimensions, providing insights about the devices and the surrounding fluid as well. CFD models are capable of accounting for blockage effects by directly modelling the flow itself, but at a high cost in terms of time and computational effort.
Different approaches for modelling actual turbines have been implemented within CFD codes such as actuator disks, lines or CFD-BEM models. These are capable of providing different levels of detail depending on the implementation and requirements, ranging from simple wake simulations to actual turbine design for particular cases (see, e.g. Shen et al. 2012;Hunter et al. 2015;Schluntz and Willden 2015;Wimshurst and Willden 2017). The most complex approach, however, consists of fully blade-resolved simulations, where the blade surfaces are modelled as solid boundaries. This approach provides insights about the development of the flow across both the chord and blade-span dimensions, accounting inherently for effects such as spanwise flow, flow separation, etc. (see, e.g. Tampier et al. 2017;Wimshurst and Willden 2016).
From an engineering perspective, blockage is relevant as it changes significantly the conditions in which a turbine operates, and corrections are likely to be required. Blockage corrections may be required to translate laboratory-scale experimental results to full-scale conditions, or to calculate the expected power and loads on tidal turbines deployed in blocked-flow conditions. For these reasons, the availability of quick blockage correction methods that can be employed for results obtained from different sources (experiments, CFD or BEM) is valuable when assessing performance and behaviour under different operational and deployment conditions. The present work analyses a series of steady-state, bladeresolved, Reynolds-averaged Navier-Stokes (RANS) CFD simulations of a horizontal axis tidal turbine under different blockage ratios representing conditions ranging from highly blocked flow to quasi-open-water over a range of different tip-speed ratios. The CFD results are used to compare the blockage correction methods described by Glauert (1933), Maskell (1963), Pope and Harper (1966), Mikkelsen and Sørensen (2002), Bahaj et al. (2007) and Werle (2010) and a correction method proposed by the authors by applying them to the different blocked CFD results and comparing them with those of the quasi open-flow condition.
(a) Turbine frontal view and hydrofoil section cuts.

Analysis cases
A tidal turbine designed by Schluntz and Willden (2015) and modified by Wimshurst and Willden (2016) in the root region was simulated under different operational conditions. This rotor, referred to as Sch15B, was designed and optimised to operate under a blockage ratio of β = 0.196, with β defined as the ratio of the turbine swept area to the channel crosssectional area.
The Sch15B rotor is a three-bladed 20-m diameter turbine developed with a single RISØ A1-24 aerofoil and a thickened trailing edge from 25% to the span up to the tip, as described in Wimshurst and Willden (2016). A render of the full turbine and its hydrofoil cross sections are shown for illustration purposes in Fig. 1a. The hydrofoil section is centred at the quarter-chord location and the blade has solidity σ and twist γ distributions that are shown in Fig. 1b, c, respectively, with the solidity defined as from which N is the number of blades, c the local chord and r the local radius. Five sets of blockage simulations were prepared, with β ∈ [0.01, 0.05, 0.10, 0.20, 0.40]. These cases cover a wide range, from the quasi-open-flow condition (β = 0.01) (Gar-  Table 1. For each blockage condition, the turbine was simulated at four different tip-speed ratios (defined in Eq. 2) λ ∈ [4.0, 5.0, 6.0, 7.0]. This covers a range of likely operating conditions for tidal turbines limited by flow separation at lower tip-speed ratios and cavitation inception at the higher end (Wimshurst et al. 2018).
Several quantities are extracted from the CFD simulations to enable comparison with and between the different blockage correction methods. Thrust and torque are obtained by integrating the static pressures over the turbine surfaces, and then converted to non-dimensional power and thrust coefficients according to Eqs. 3 and 4, where ρ is the fluid density, T is turbine thrust and the power P = Qω is product of the rotor torque Q and rotational speed ω. U is the undisturbed flow speed far upstream of the turbine and A is the rotor swept area. (3) In addition to these results, the axial induction factors a are also obtained from the simulations, with a defined as where U 1 is the average flow velocity at the turbine plane.
Rearranging the previous equation, a can be easily calculated as The axial induction factor was extracted from the simulations by obtaining the velocity U 1 averaging the axial velocities in annular sections located at 1.5 local chords upstream and downstream from the turbine plane. Finally, qualitative information is also obtained from the CFD in the form of velocity profiles, providing insights about the wake development and interaction effects between the rotor and the boundaries.

Blockage correction methods
The blockage β is defined as the ratio of the cross-sectional area of a tank or channel C where a turbine is operating and the turbine swept area A: Blockage correction methods typically seek to obtain a relationship between the upstream tank velocity U T (see Fig. 2) and the equivalent upstream velocity U F in free-flow conditions where the turbine would be operating at the same conditions such as angle of attack, mass-flux through the rotor, turbine RPM and turbine thrust (Bahaj et al. 2007). Most blockage correction methods are based on traditional actuator disc theory adapted for the particular case of a blocked flow and attempt to obtain an equivalent free-flow velocity U F as function of different variables such as blockage β, thrust coefficient C T , wake expansion or induction factors. In this kind of approach, the turbine is represented as an actuator disc across which a static pressure discontinuity occurs, and thrust can be defined as in the following equation, where P + 1 and P − 1 are the average static pressures on both sides of the actuator disc: Figure 2 represents the flow past an idealised turbine in a blocked environment. The undisturbed flow velocity in the blocked environment U T at the inlet is considered uniform across the cross-sectional area of the tunnel C. The flow speed reduces to U 1 through the rotor at the turbine plane, and by conservation of mass is accelerated to U 3 in the bypass region. At some distance downstream of the turbine plane the area of the wake before mixing is defined as A W , with a reduced flow speed U 2 , a recovered pressure P 2 , and a flow speed in the bypass region that is further accelerated to U 4 .
For each blockage correction method, a relation between the tank velocity U T and a corresponding free-stream velocity U F is given, leading to the following corrections for λ, C T and C P determined for blocked conditions: where the subscript c is denotes corrected or equivalent freeflow values. In the following sections, each of the blockage correction methods is presented and briefly discussed.

Glauert (1933)
The blockage correction of Glauert (1933), intended initially for aircraft propellers, considers an equivalent free-stream velocity at equal thrust in the tank or wind tunnel. The resulting correction: can be used for propellers and turbines, but is limited to thrust coefficients C T < 1. This can be problematic for highly loaded turbines such as tidal turbines where C T can be significantly greater than unity, and thus limits the applicability of the method for tidal applications.

Maskell (1965)
In 1965, Maskell (1963) proposed a blockage correction method for bluff bodies, which has also been used for turbines and aerodynamic profiles. The blockage correction: with (1 − k 2 ) −1 = −C pb , the base-pressure coefficient as defined in Maskell (1963). Maskell recommended an empirical value of 2.5 for the (1 − k 2 ) −1 term due to its relatively small variation over the range of blockage ratios covered in that study. Due to the large range of blockage ratios in the present study, a linear regression was made from the data presented in the original reference giving (1 − k 2 ) = 0.3551 − 5.1050β. Other modifications to the (1 − k 2 ) term have been proposed for other geometries or different aspect ratios, e.g. Mercker (1986) Hackett and Cooper (2001) and Leuschen and Mebarki (2012).

Pope and Harper (1966)
Pope and Harper (1966) presented a blockage correction which has been widely used for low-speed wind tunnel tests for blockages β from 1 to 10%. The correction utilises an empirical factor t to account for blockage effects. The authors suggest t = 1/4β, although this does not necessarily apply to all geometries or for complex wake flows.
The correction where u = 1−a. This correction requires the induction factor a, which can be obtained directly from CFD simulations such as those described in Sect. 2 or indirectly from the thrust coefficient. The latter can be calculated from momentum theory for low thrust coefficients, but has to be corrected empirically for larger thrust coefficients, as described for example in Burton et al. (2001). Bahaj et al. (2007) presented a blockage correction method based on an actuator disc model. The correction

Bahaj et al. (2007)
is solved iteratively, with the term U 1 obtained from continuity arguments as described in Bahaj et al. (2007). The method only requires blockage and and thrust information (β and C T ), and has been used for hydrokinetic turbines by several authors, e.g.

Werle (2010)
Werle (2010) presented a blockage correction, aimed as an easily applicable correction for measurements near the maximum extractable power. The correction only requires the blockage β, but should not be used to correct C T and C P as for the other methods as described in Eqs. 10 and 11. Alternatively, Werle proposes the following corrections for C T and C P : It should be noted that Werle recommended this approach only as a reasonable first approximation near the maximum power coefficient.

Continuity correction
A correction based on continuity arguments and CFD simulations is proposed by the authors to provide a comparison with the correction methods listed above and to expand the capabilities of the methods available from the literature by introducing an empirical factor. We are looking for a free-stream velocity U F which would give us an equivalent thrust to that of the blocked-flow condition in an unconstrained flow. The effect of blockage is to increase turbine thrust for a given through turbine flow speed, thus U F > U T . We also know that the velocity in the bypass between the wake and tunnel, downstream of the turbine, is greater than U T , so we assume that there is an arbitrary plane where the velocity in the bypass would be equal to U F .
Then we can define the cross-sectional area of the wake in this reference plane as where δ F is an empirical factor. Values for δ F vary such that Defining the wake velocity in this reference plane as U R and then applying the continuity equation: Using Eq. 5, we can then define U R as Following continuity in the bypass flow, and recalling that in this reference plane the velocity in the bypass is equal to U F , then which, using Eqs. 7 and 21, can be rearranged as Then through a continuity balance between the turbine plane and the channel inflow, and using the definition of Eq. 25: Finally, Eq. 27 can be rearranged to define the correction factor as As for the Mikkelsen and Sørensen (2002) correction, this method requires the induction factor a, which can be obtained directly from CFD simulations or estimated empirically. The influence of the empirical factor δ F will be analysed in more detail in Sect. 6.

Error definition
To assess the performance of the different methods, the blockage-corrected results were analysed and compared to a quasi free-flow condition, with error metrics as defined below. The β = 0.01 simulations represented the case of quasi-free-flow condition as blockage effects are described as negligible once β becomes smaller than 2% (Garrett and Cummins 2007).
For the analysis, the normalised error is calculated as  where the subscript X corresponds to thrust or power (T or P), subscript c corresponds to the corrected results, and subscript 0 to the unblocked results (β = 0.01). The relative error is obtained for each method along the corrected range of λ c for each blockage ratio. The reader should note that the tip-speed ratio is corrected by the different methods in addition to the power and thrust coefficients. To obtain the corresponding C X 0 values for the quantification of error, a third-order spline interpolation was employed to estimate the unblocked values at λ c . A global overview of error for each method was obtained from the average magnitude of the error values for E C T and E C P :

Reynolds-averaged Navier-Stokes equations
The dataset for this study was obtained by numerically solving the steady-state Reynolds-averaged Navier-Stokes equations (RANS) under constant density and constant temperature assumptions (Ferziger and Peric 2001) using the commercial software ANSYS Fluent 19.0. The multiple reference frame (MRF) model introduced by Luo et al. (1994) was employed to model the turbine rotation under steadystate conditions. This approach allows the flow behaviour to be predicted using two rotating frames of reference, one that rotates at the turbine speed and the second one being stationary. The velocities, under this model, are defined in terms of the absolute velocity v, the relative velocity v r and the whirl velocity u r = ω × r, with ω the rotational speed and r the position vector relative to the centre of the local frame of reference. These quantities are related as The governing equations are a modified version of the Navier-Stokes equations that include centrifugal and Coriolis acceleration terms. In a relative velocity formulation (Fluent 2018), these can be written as where t is the time variable, ρ the flow density, p the static pressure, T the stress tensor and F the external body forces. Equations 32 and 33 are solved in each region. If the region is stationary (ω = 0), the equations are reduced to their stationary form. At the interfaces between regions with different reference frames, an implicit coupling is performed using Eq. 31, transforming the velocities and velocity gradients to a common reference frame, while scalar quantities remain invariant.
The k−ω SST turbulence model described by Menter (1994), with the model coefficients updated as described in Menter et al. (2003), was employed to close the RANS equations. The k−ω SST turbulence model has been found to compare well to higher order turbulence closures in numerical studies and experiments of turbines (Abolghasemi et al. 2016;Shives and Crawford 2016). The non-dimensional wall distance y + was maintained in the wall-modelling range (y + ∈ [30, 300]) near fluid boundaries by appropriately scaling the size of wall-adjacent cells.
To solve the discretised problem, we employed a pressurebased coupled solver, as described in Fluent (2018). Convergence was achieved after 15,000-20,000 iterations, with reduction of the continuity and velocity residual by at least six orders of magnitude, and five orders of magnitude for the turbulence scalars.

Computational domain and mesh study
A solution domain that exploits the azimuthal symmetry of the problem was defined using periodic boundary conditions requiring only one third of the problem to be modelled. The 120 • wedge consists in two structured mesh subdomains referred to as the rotatory and stationary subdomains. The meshes were built using ICEM 19.0.
The rotatory domain was based on the work presented by Wimshurst and Willden (2016). The turbine blade was modelled with 128 cells over the chord, and 110 in the span-wise direction, remaining within the mesh-independent region.   The domain follows a Y-Mesh topology, resulting in approximately 4.2 million cells for the 120 • wedge. The mesh on the blade surface is shown in Fig. 3.
A mesh sensitivity study also analysed the influence of different parameters inside the O-grid layer surrounding the blade (see Fig. 5) was presented by Zilic de Arcos et al. In particular, the number of layers of the O-grid around the profile n l , the thickness of the first layer Y and the growth rate G r were studied. Following Zilic de Arcos et al. (2019), a growth rate of 1.05, 25 layers and a first-layer thickness of 2.5 × 10 −4 m were selected for the analysis due to the low mesh sensitivity as well as avoiding the buffer region between the wall-modelling and wall-solving approaches (5 < y + < 30). The rest of the domain was adjusted accordingly to provide smooth transitions with the O-grid region.
The mesh in the stationary domain follows an azimuthal O-grid strategy. Considering that the interest of this study is placed on the forces over the blade, the resolution of the stationary domain is of lower importance. Nonetheless, a mesh independence study was performed using the model described in this study by simulating three meshes with different numbers of cells n c . The mesh convergence assessment method presented by Celik and Karatekin (1997) was employed to determine that the meshes are within the asymptotic convergence region. Specifically, evaluating the parameter R (Eq. 34) for the thrust and power coefficients indicates that thrust displays an oscillatory convergence (R < 0) whereas power shows a monotonic convergence (0 < R < 1).
with F i being the thrust or power coefficient values and the subindex i = 1, 2, 3 indicating the coarse, medium and fine meshes, respectively.
Following this analysis, the mesh convergence method uses Richardson's extrapolation to approximate the analysed quantity F ext in the limit where the mesh typical size h i is zero, based on the apparent order of the method n and the refinement ratios a i such that F ext = a n 2 F 1 − F 2 / a n 2 − 1 Finally, the discretization error is described as Notice that, due to the oscillatory convergence of thrust, a sign change is used in Eq. 35 according to Celik and Karatekin (1997).
A detailed discussion of this grid-convergence assessment methodology, as well as the relationship between e ext and the grid-convergence index (Roache 1994) can be found in Celik and Karatekin (1997). The results of the mesh convergence for C T and C P , as well as the relative error E R and the extrapolated error e ext according to the described methodology are provided in Table 2.
Considering the mesh sensitivity analysis results, the medium-resolution configuration was deemed sufficient and was employed for the rest of the analysed cases. A longitudinal cut of this mesh can be seen in Fig. 4.
The boundary conditions were configured following (Wimshurst and Willden 2016): turbine surfaces were defined as no-slip smooth walls; the inlet defined with a steady uniform inflow velocity of 4.5 m/s, a turbulent length scale of 0.70 times the turbine diameter and a turbulent intensity of 10%; the outer boundary as a symmetry condition to model a frictionless wall; the outlet defined by a constant pressure equal to the undisturbed flow; and axial symmetry was exploited using non-conformal periodic boundary conditions on the walls of the 120 • wedge (Fig. 5).
The modelling approach based on the k − ω SST turbulence model, the periodic symmetry condition, the MRF steady-state approach, the boundary conditions and other model characteristics have been previously validated through comparison of the simulated results of the MEXICO rotor to the experiments published by Schepers et al. (2012) and Boorsma and Schepers (2014) in terms of spanwise force distribution and sectional pressure distributions, showing a good agreement for both quantities. The full details of this study are presented in Wimshurst and Willden (2017).

Uncorrected CFD results for blockage variation
As discussed in Sect. 2 and detailed in Table 1, CFD simulations were carried out for a blockage range from 0.01 ≤ β ≤ 0.40. Figure 6 shows the thrust and power coefficients for each blockage condition as functions of the uncorrected tip-speed ratio λ. The lowest blockage ratio (β = 0.01) is considered as a quasi-free-flow condition and is used as a reference to benchmark the correction methods.
Transverse and streamwise contours of the axial velocity for the different blockage ratios are presented in Figs. 7 and 8 at a tip-speed ratio of λ = 7.0. It can be observed from these figures that the flow in the bypass region is accelerated as blockage increases. The mass flux through the rotor also increases with blockage ratio when λ is held constant as there is greater resistance to acceleration of the bypass flow in constrained flows. The greater shear stress developed between the bypass and through-turbine flows means that there is greater recovery of the flow speed in the wake of the turbine as blockage increases at comparable streamwise positions.

Proposed method
As mentioned in Sect. 3.7, the proposed continuity method requires an empirical correction factor δ F . For this study, two different cases were evaluated: one with δ F defined as a function of β and another one where δ F = 0.1 is considered constant for all the different cases.
The curve shown in Fig. 9 shows the function δ F = f (β) that minimises the mean average error defined in Sect. 4 at each blockage ratio across the full range of tip-speed ratios.

Results with blockage correction
The presented results for β = 0.01 to 0.40 were corrected using each one of the presented methods, as shown in Fig. 10. This set of graphs gives a first overview of the effectiveness of each one of the methods. Notice that, for the Glauert correction, results are shown only when C T < 1.

Error analysis
The results for the error are shown in Table 3 for each correction method and in Fig. 11 with the power and thrust errors plotted as functions of both λ c and β.
From Table 3, it can be observed that the proposed method displays the smallest error both if used with the optimum δ F (β) or the constant δ F = 0.1 value, followed closely by Bahaj's and Mikkelsen and Sørensen methods. Pope and Harper and Maskell corrections also display an acceptable performance if a low blockage is considered. For the Glauert method, even if only results where C T < 1 are considered, a higher total error can be observed. The Werle method, as proposed in Werle (2010), also demonstrates a high level of error, and is, therefore, not recommended as a correction method for any blockage and tip-speed ratio range.
To study the dependency of the mean absolute error over blockage β and reference tip-speed ratio λ, results were plotted over these variables, respectively. 1 In Fig. 11, MAE T and MAE P are shown, each as a function of β and λ, for each correction method.
The results were only plotted up to 15% error, considering that methods providing larger errors should not be taken into account for blockage correction. Therefore, most of Werle's correction points exceed the plotting range. Also, the Glauert correction, due to the limited number of data points available, was omitted on this comparison.
From the plots, a blockage dependency is observed for both MAE T and MAE P on most methods (Fig. 11), which can be expected considering that larger corrections are necessary for higher blockage ratios.
The Maskell correction shows a more complex dependency of blockage (Fig. 11). This becomes evident from Fig. 10, where an overestimation of C T c and C Pc is observed for β = 0.4, contrary to the small underestimation which can be observed for the remaining blockage ratios.
Regarding the Pope and Harper correction, a blockage and tip-speed ratio dependency is observed. This method, despite its simplicity, could be useful as an estimation for relatively small blockage values.  The Mikkelsen and Sørensen correction shows, in addition to a small value of MAE T and MAE P , a very small βand λ-dependency along with the smallest dispersion of all methods (Table 3), and is, therefore, recommended as a suitable blockage correction method for a wide range of blockage and tip-speed ratios, especially if results of induction factors are available and must not be estimated.
For the Bahaj correction method, a relatively small β and λ dependency can be observed, with a relatively small overall error for both thrust and power coefficient corrections. Considering that this correction only requires β and C T from experiments (or computations), it can be especially suitable for practical applications as no induction factors are required.
The method proposed in our work is shown for two cases of the empirical factor, with δ F = 0.1 and δ F = f (β).
These have a similar behaviour, with the β-dependent δ F case showing a relatively small improvement in the accuracy of the prediction over the constant δ F case. Both cases display higher levels of error at the highest blockage condition, although the magnitude of the error is similar to the Bahaj and Mikkelsen and Sørensen methods.
The inclusion of the empirical factor is also interesting from a different perspective. Systematic series of δ F could be developed for, not only different blockage conditions, but also different channel shapes and aspect ratios, something that is not considered by any of the presented methods. Additional evaluation is required to further validate the proposed correction method.

Conclusions
The use of a RANS-CFD tool to model the performance of a hydrokinetic turbine over a wide range of blockage conditions can provide useful data to analyse existing blockage correction methods, propose improvements to these or to develop new methods for specific applications.
From the presented results, the authors recommend the use of the Mikkelsen and Sørensen method for the correction of the blockage for horizontal axis turbines subject to experimental test or numerical simulations under conditions similar to the presented here, considering its low error, low dispersion and low dependency on β and λ. If the induction factors are not available, the Bahaj correction is recommended instead, as it shows an even smaller error than the previous method but a larger spread on the corrected results. Despite showing very encouraging results, the proposed method should be validated for other application cases before it can be recommended for its general use.
It is also recommended, considering the presented results, to avoid using the Werle and Glauert methods, providing that other simple and more reliable methods such as those by Maskell or Pope and Harper are available, when quick preliminary corrections at low blockages are required.
Caution must be taken as this study only addressed the problem of homogeneous blockage. Future work should analyse series of different blockage ratios for different turbines, tank sections (e.g. rectangular sections) or other arrangements such as non-cylindrical channels. The proposed model, with its empirical component, could be a viable option to address the problem of non-homogeneous blockage.