Prediction of Seabed Scour Induced by Full-Scale Darrieus-Type Tidal Current Turbine

: Scour induced by a Darrieus-type tidal current turbine was investigated by using a joint numerical and experimental method with emphasis on the scour process of a full-scale turbine. This work proposes a new numerical method to estimate turbine scour developments, followed by model validation through experimental data in the initial stage. The small-scale numerical model was further extended to a full-scale model for the prediction of turbine scour. The numerical model consists of (1) k- 𝜔 turbulence closure, (2) a sediment transport model, and (3) a sediment slide model. The transient-state model was coupled with a morphologic model to calculate scour development. A dynamic mesh updating technique was implemented, enabling the autoupdate of data for the grid nodes of the seabed at each time step. Comparisons between the numerical results and the experimental measurements showed that the proposed model was able to capture the main features of the scour process. However, the numerical model underestimated about 15%–20% of the equilibrium scour depth than experimental data. An investigation of the temporal and spatial development of seabed scour around a full-scale Darrieus-type tidal current turbine is demonstrated. This work concludes that the proposed numerical model can effectively predict the scour process of tidal current turbines, and the rotating rotor has a significant impact on the equilibrium scour depth for full-scale turbines. and the influence of complex flow fields induced by a rotating rotor. The interaction between the turbine’s hydrodynamics and local scour around the seabed, and their results showed that axial velocity increased about 5% of initial flow, which is less than the Darrieus-type turbine in this project. This suggests that the impact of Darrieus-type turbines on scour hole development is greater than that of horizontal-axis turbines.


Introduction
Tidal current energy is a kind of clean and pollution-free marine renewable energy that can be extracted by a tidal current turbine. Scour around turbine is an important engineering issue that causes structural instability [1]. Numerous tidal current turbine projects have been carried out to harness marine energy in the past few years [2]. However, little is known about the interaction between seabed scour and turbine devices, and most previous studies have focused mainly on the hydrodynamic performance of turbines [3].
Predictions about scour depth and profiles are important to reduce the risk of turbine structural instability induced by scour. The support structure for tidal current turbines can be categorized into four main types: a gravity structure, monopile structure, jacket structure, and floating structure [4]. Of these four types, the monopile structure is a common type for tidal current turbine projects. The mechanism of turbine scour is similar to local scour in piles [5]. Seabed scour is caused by the constant change of flow around a foundation pile: (1) the flow contracts at the side of the foundation; (2) sustained downflow exists in front of the monopile foundation; (3) a horseshoe vortex is formed in front of the foundation; and (4) lee-wake vortices occur downstream. These flow changes work together to cause seabed scour. As time goes by, the downflow near the bottom of the scour hole decreases as the scour hole depth increases. Scour equilibrium is obtained when the magnitude of downflow ahead of the cylinder can no longer dislodge sediment at the surface of the scour hole [6].
Turbine scour is more profound under the turbines than around the piles for similar flow conditions due to the influence of a rotating rotor. The impact of a rotor on turbine scour has been emphasized by many researchers. Wang et al. [7] measured the slipstream characteristics of ocean stream turbines and suggested that the height of tip clearance should be at least one turbine diameter to decrease its impact on the seabed. Hill et al. [8] found that equilibrium scour depth was much deeper when the rotor was installed upstream compared to downstream installation. Zhang et al. [9] carried out a series of model tests for horizontal axis tidal current turbine scour and proposed an empirical model for scour depth prediction. Sun et al. [10] proposed semiempirical equations for temporal scour depth of Darrieus-type tidal current turbines based on experiments. Besides the aforementioned experimental works, computational fluid dynamic (CFD) simulation is also a great tool for scour prediction. Chen and Lam [11] studied slipstreams between turbines and seabeds using OpenFOAM (open source CFD software) and found that tip clearance was an important parameter in predicting horizontal axis turbine scour depth. Roulund et al. [12] coupled the flow model with a morphologic model to calculate scour at piles. Their results showed that equilibrium scour depth obtained from CFD simulation was 15%-30% smaller than experimental data due to the unsteady effects of flow not being taken into account.
The aforementioned experimental studies of turbine scour mainly used small-scale models. No full-scale experiment for turbine scour has been used due to the limitation of experimental conditions and the high cost. However, numerical simulation is a good tool to predict the scour process for largescale turbines at a relatively low cost. Some full-scale simulation works have been used to predict the hydrodynamics of tidal current turbines or the scour process. Ettema et al. [13] discussed the scale effect in pier-scour experiments and illustrated that local scour depth in a laboratory may be proportionately greater than the depths occurring in actual projects, and they suggested that CFD modeling be used to reduce the scale effects. Huang et al. [14] conducted a 3D CFD model to examine scale effects on turbulent flow and sediment scour. Dai and Lam [15] proposed a 1-MW full-scale turbine model and completed a numerical investigation for its power coefficient and structural loading. Penny et al. [16] conducted full-scale testing of the SCHOTTEL tidal turbine at Queen's University Belfast's tidal site and found that the maximum system efficiency of the turbine reached 35% in flows up to 2.5 m/s. Gonzalez-Gorbeña et al. [17] presented a study for the optimization of an instream tidal energy converter array layout and showed that large amounts of energy extraction may cause significant erosion. Burkow and Griebel [18] proposed a full three-dimensional numerical simulation model for scour at a rectangular obstacle and compared the scaled-up numerical model to a small-scale experimental model. Zhang et al. [19] developed a mathematical model in the FLOW-3D package to numerically investigate the fluid-structure interaction and its induced scour around full-scale horizontal-axis tidal current turbines. However, a limited study focused on scour induced by Darrieus-type tidal current turbines, especially full-scale turbines. To promote scour prediction in full-scale Darrieus-type tidal current turbines, a 3D CFD model and experiments were conducted to examine the temporal and spatial development of full-scale turbines.
This paper presents the progress made in the numerical simulations of turbine scour. Sediment transport and slide models using user-defined function (UDF) code associated with the fluid flow model in FLUENT (commercial CFD software) were implemented to simulate the actual scour process. After validation through experimental data, seabed scour induced by full-scale turbines was studied in detail.

Methodology
Seabed scour is a complex physical process around a Darrieus-type tidal current turbine. The scour depth depends on the size of a foundation pile and the influence of complex flow fields induced by a rotating rotor. The interaction between the turbine's hydrodynamics and local scour around the foundation was a key factor to investigate in this project. To investigate this factor, a CFD model was proposed. In this project, the computation of tidal current turbine scour was performed using the commercial package FLUENT 2018. Scale-down physical model tests and scale-up numerical simulations were carried out together for model validation. Finally, the validated model was used to investigate the seabed scour induced by a full-scale tidal current turbine. Figure 1 shows the research framework of this project.

Hydrodynamic Model
The hydrodynamic model includes geometry creation, grid generation, boundary condition setting, turbulence model selection, and result extraction, as presented in the following sections.

Governing Equations
The flow conditions around turbines can be calculated through a finite-volume numerical model that solves incompressible Reynolds-averaged Navier-Stokes equations [12], as shown in Equation (1): in which is the fluid density, is the ith component of velocity, p is the dynamic pressure, is the dynamic viscosity, is the turbulent viscosity, and is calculated by a turbulence model, as detailed in the following section Semi-implicit method for pressure-linked equation(A SIMPLE). algorithm [20] was used to calculate the flow conditions. A SIMPLE algorithm is a kind of pressurebased solver, wherein the constraint of mass conservation of the velocity field is achieved by solving a pressure equation. Pressure-velocity coupling is achieved by using face flux equations to derive an additional condition for pressure by reformatting the continuity equation [21], as given by Equation (2), where Jf is the mass flux through face. In addition, we used a second-order upwind scheme for spatial discretization: In the SIMPLE algorithm, the pressure field is calculated and the velocity field is corrected to satisfy the continuity equation [21], as given by

Turbulence Model
The k-model [22] was selected as the turbulence model to resolve the unknown Reynolds stress terms: the k-model shows great performance for shear flow spreading, as proven by the previous studies of Roulund et al. [11] and Baykal et al. [23]. Specifically, the shear stresses transport (SST) k-model, which was proposed by Menter [24], combines the positive features of both the kand k-turbulence models and was chosen for the simulation. More specifically, the kturbulence model can accurately predict the near-wall region. In contrast, the k-turbulence model can accurately predict the far field. A blending function, F, was applied to combine both models. The factor of the blending function for the k-model reduces to zero at the boundary layer, and the kmodel replaces the k-model to predict the free stream region.
Actually, three versions of the k-model are the original k-model, the k-BSL (Baseline) model, and the k-SST model. Menter [24] has compared flow past a NACA 4412 airfoil calculated by these three models. The results showed that the SST k-model gave the most accurate results for the tested adverse pressure gradient flow cases. Hence, the SST k-model was chosen to simulate the flow field with strong adverse pressure effects in this study.

Boundary Conditions
The boundaries of the computational domain are inlet, outlet, symmetry boundaries, and walls, which can be seen in Figure 2: (1) Inlet and outlet boundary conditions: For the inlet boundary, the uniform velocity inlet condition was used, and the initial pascal pressure was set as 0. The outlet boundary pressure outlet was applied for all quantities. The turbulence intensity at the inlet/outlet boundaries was 5% for medium intensity as a default value, which is suggested in Fluent Help [21]; (2) Symmetry boundaries: The symmetry boundary was applied to simulate the free surface condition. At the symmetry boundaries, Neumann conditions were applied for variables in the k-turbulence model. It was assumed that a "lid" was on the top surface. The symmetry boundary condition is an effective method to simulate the free surface, and it has been applied in many scour simulations, such as in References [12] and [23]; (3) Walls: The sides of the computational domain, surfaces of foundation, and turbine blades were set as no-slip wall boundaries. The roughness height was set as zero to simulate a smooth surface. The sediment bottom was set as a rough wall boundary. The roughness height of the sediment bottom was set as 2.5, which can simulate the true sediment bed (as suggested by Roulund) [12]. In the current study, it was assumed that the effects of link rods between the blade and monopile on the velocity and turbulence intensity around the turbines were negligible. The dimensions of the computational domain were 4.5D, 4D, and 3H in the x, y, and z directions, respectively. The flow direction was from -x to x. D is turbine rotor diameter, and H is turbine height.

Computational Mesh
The grid meshing of the numerical model was completed using an ICEM preprocessor. The overall mesh can be seen in Figure 3. A suitable mesh setting is important to ensure the accuracy of a CFD calculation. The overall mesh combined the advantages of structured grids (upper domain contains turbine rotor) and unstructured grids (lower domain between rotor and seabed). The upper structured grids can speed up the numerical calculation. The lower unstructured grids can better adapt the complicated geometry to update the bed boundary data under a scour process during the solving process. The blade surface was set as moving mesh corresponding to its adjacent cell zone. The sliding mesh modeling technique was used to simulate the rotating rotor. In this model, the computational domain consists of the rotating part and the stationary part, which is shown in Figure 4a. The transient relative motion was used to predict the true transient interaction of the flow between a stator and rotor passage. The interface position was updated at each time step, changing the relative positions of the grids on each side of the interface. The rotational motion was simulated by allowing the rotating part to rotate at a constant angular velocity.
The grid resolution close to the rotor wall needed to be further refined in order to ensure the accurate prediction of flow at the boundary layer. Hence the rotor mesh density was higher than the stator mesh density. The grid pattern of the blade boundary layer is shown in Figure 4b.

Grid Independence Analysis
The total cell numbers in the CFD mesh can influence the predictions due to the numerical approximations using RANS (Reynolds-averaged Navier-Stokes) equations. However, the increase of mesh density also requires more computational power. Hence, the grid independence test is important to obtain a suitable cell number of mesh for particular cases. A grid independent analysis was conducted to investigate the sensitivity of the grid to the change in Cp values. The turbine power coefficient Cp is given by where T is turbine torque, is turbine rotational velocity, is the density of water, A is the swept area of the rotor, and V is free stream flow velocity. Table 1 shows the cases of grid independence tests on the mesh grid numbers and Cp values to reach grid independence. The results show that Case 1 with 1,285,567 cells had a Cp value of 16.4%. The Cp value increased to 18.8% for Case 2 with an additional 394,087 cells, and the Cp value increased to 19.7% for Case 3 with an additional 367,024 cells. After Case 3, the validation of results was insignificant with further increases in the cell number. Therefore, Case 3 with 2,046,678 cells was the selected grid-independent case.

Morphological Model
The flow field can be calculated through the aforementioned hydrodynamic model, and this section gives a brief description of the morphological model that describes the sediment movement during the scour process.

Sediment Transport Model
The judgment of sediment incipience is the first stage in simulating turbine scour. The three main methods are [25] (1) judged by critical shear stress, (2) judged by critical velocity, and (3) judged by critical power. In our sediment transport model, the first method was chosen by considering the influences of turbulent flow. For no-slip wall conditions, ANSYS Fluent uses the properties of the flow adjacent to the wall/fluid boundary to predict the shear stresses on the fluid at the wall. The sediment particles are initiated when the bed shear stress is higher than a designated value. The sediment incipient shear stress , can be calculated as [26] ( ) where , is critical shear stress; s is the relative density of the sediment particle; is the medium diameter of the sediment particle; and is the critical Shields number. Equation (6) is where ( ) However, the riverbed is not always flat. The slope correction for critical shear stress is of importance. The critical shear stress on the slope can be calculated by the following equation [26]: where is the bed slope angle of repose; C is the lift coefficient; C is the drag coefficient; C /C is taken as a general empirical value of 0.85 [27]; is the angle between the horizontal axis of the slope and the flow direction; and is the angle between the slope and the horizontal plane. The slope model can be seen in Figure 5. The second stage of our sediment transport model is sediment transport rate. The sediment transport rate considering the seabed inclination angle can be calculated as [28] where is the sediment transport rate of the planar single-width volume; h is the vertical height of the seabed; is seabed shear stress; and is the shear stress component in the direction. The second section of is a correction of slope. The empirical coefficient C was suggested to be 1.5-2.3 by Brørs [28]. In this project, C = 1.5, and can be expressed as [28] ( ) ( ) Based on the equilibrium of sediment transport, the transient variations of seabed elevations can be derived in terms of the transport rates: where h is seabed elevation, and n is the void ratio of sediment. By computing the abovementioned equations, the transient elevation change of the seabed boundary can be updated according to h in real time. The sediment mass conservation is not controlled in the proposed model. The scour hole around the pile was simulated, and the sand deposition was neglected in this project.

Sand Slide
Previous research [29] has pointed out that the collapse of a seabed occurs when the local bed slope exceeds the angle of repose. As a result, shear failures occur at these locations. Roulund [12] and Liang [30] proposed two different sand slide models to solve this problem. This project extended Liang's sand slide model from 2D scour to 3D scour simulations. Figure 6 shows a diagram of the sand slide model. The 3D sand slide model works as follows: (1) Update the seabed elevation change based on the sediment transport model at each time step. After that, scan all the nodes and center points of the mesh grid on the seabed boundary; (2) If the angle between grid node A( , , ) and grid node B( , , ) is greater than the bed slope angle of repose , then the model raises or reduces the z coordinate of grid node A and B, making sure the slope angle is less than . The vertical displacements of A and B are the same, which ensures the conservation of sediment in the scour process. Equation (11) is (3) Repeat the above process 6-10 times, ensuring all of the inclination angles are under .

Updating Strategy
Dynamic meshing updating technology is used to simulate the seabed change at each time step. The procedure in the computations is as follows: (1) Calculate the flow field; (2) Calculate the sediment transport due to bed load; (3) Update the seabed boundary; (4) Check the sand slide and use the sand slide model to correct the slope angle; (5) Turn to the next time step and return to step 1. The morphologic time step in the calculation is 0.01 s initially and gradually rises to 0.1 s as the scour speed decreases. This is to improve computing efficiency. Over time, the bed shear stress becomes less as the scour hole gets deeper. Finally, the bed shear stress in the scour hole is under critical shear stress, and the seabed profile has an insignificant change over time. This means that turbine scour achieves equilibrium.

Experimental Setup
Experiments on scaled-down Darrieus-type tidal current turbines were carried out in a recirculating flume to validate the proposed turbine scour model. This flume was a purpose-made horizontal recirculating flume with water driven by a propeller, which was different from a traditional vertical recirculating flume driven by a pump [31]. The experimental installation is shown in Figure 7. The maximum scour depth of each test was measured by laser distance meter. As discussed in Section 3, the experimental data on the maximum scour depth at different tip clearances were used for analysis and validation. The main experimental parameters can be found in Table 2. More details of the experimental setup are elaborated on in the paper in Reference [10].

Numerical Model Setup
Two groups of test cases were computed to test the numerical scour model and its implementation, as illustrated in Table 3 [14]. In this project, the Reynolds number of turbine wake induced by Darrieus-type turbine models was around 25,640. Rajaratnam [32] has suggested that the effect of viscosity could be neglected as long as the Reynolds number is more than 10,000. The Reynolds number in the turbine scour tests was high enough to neglect the viscosity. As seen in Table 3, two groups of numerical test cases were simulated. The first group had two scour models (T1 and T2). Model T1 was a pile scour model without a turbine rotor. This model was aimed toward validating the numerical scour model with a comparison of Melville's experiment of scour at piers [33]. In Melville's experiment, a pile was installed at the center of a water flume under a steady current. The main parameters in the T1 model were set up the same as in Melville's experiment. The 5.08-cm diameter pile was installed inside, the medium diameter of the sand ( ) was designed at 0.385 mm, and the average flow velocity was set to be 0.25 m/s. The flow field and equilibrium scour shapes were chosen to act as validation. Model T2 was a turbine scour model, including a rotating rotor with the same pile diameter as model T1. This model was aimed at investigating the impact of a rotating rotor on seabed scour compared to single-pile scour. On the other hand, this model could also validate the scour model, including the simulation of the turbine's hydrodynamics.
The second group had four CFD models with different tip clearances from C/H = 0.25 to C/H = 1.0 (T3-T6). This group was also aimed at validating the numerical model. The scour results of these turbine models were used for comparison with the experimental model tests, as presented in Section 2.3. The scales of the numerical models in this group were 100 times bigger than the experimental setup, which was the same scaling ratio as Huang's research [14]. The scale-up models can also investigate the scale effects of turbine scour. To simulate real seabed conditions, we chose the sediment diameter as 8.0 mm, which was medium-sized, as illustrated in Zhang's research [19].

Model Validation through First Group
This section shows the flow field and equilibrium scour profile of test models in group 1. Through a comparison with Melville's experiment [33], the proposed numerical scour model was approved as effective and accurate.

Flow Field Validation
Experimental results and numerical simulations of wake of Darrieus-type wind turbines investigated by Bianchini et al. [34] were used in this study to validate the flow simulation around the turbine in model T2. The dimensionless axial velocities near the turbine efflux plane are shown in Figure 8. The computational results had good agreement with Bianchini's experimental and numerical simulation results. With turbine wake, the efflux velocity shows a little increase on the sides but decreases significantly right behind the turbine. However, the low-velocity region of tidal current turbine wake is wider than wind turbine wake. This may be due to water having a higher density than air and a stronger turbine disturbance underwater. Tidal current turbines can capture more energy and disturb more section areas of fluid than wind turbines can. The flow mixture in tidal current turbine wake is more intense. Comparisons between wakes from wind turbines and tidal current turbines show a little mismatch, but limited experimental data on Darrieus-type tidal current turbine wake are available. Despite this, these two types of turbines have analogous work mechanisms with different types of fluid. Hence, the efflux velocity can show almost the same flow trend. On the basis of the validation, the simulation results of the proposed hydrodynamic model had reasonable agreement with the experimental results. This model could be accepted for flow precalculation of turbine scour.

Scour Shape Validation
On the basis of the numerical and experimental results, the topographic contours of the local scour shapes after scour equilibrium of model T1 are shown in Figure 9. A scour hole around the foundation pile can be seen in Figure 9. The maximum scour depth occurred at the sides of the foundation. This result agrees well with pile scour experiments [33] and turbine scour experiments. In the CFD model, the horizontal length of the scour hole was 2D, 2D, and 2.5D upstream, downstream, and on the side of the pile, respectively. The horizontal extent of the scour hole was a little smaller than the experimental data, especially for downstream length.
In the simulation results, the maximum scour depth in model T1 was 5.2 cm (1.02D), and the width of the scour hole was about 1.5D. In Melville's experiment, the maximum scour depth was 6 cm (1.18D), and the width of the scour hole was about 1.8D. The numerical model underestimated the equilibrium scour depth at about 13% compared to the experimental data with the same model setup. This phenomenon showed the same trend as Roulund's CFD simulation [12]. The equilibrium scour depth obtained from that simulation was 15% smaller than the experimental scour depth. The reason was that in the iteration process in transient analysis, the seabed updates after the calculation of bed shear stress. Unsteady effects such as the fluctuating components of the horseshoe vortex and lee-wake vortex flows are not taken into account. These unsteady effects can also intensify the scour level [25]. In addition, the longitudinal length of the scour hole in the numerical results was also smaller than the experimental data due to the neglected impact of lee-wake vortex, which was the main reason for the scour hole being extended behind the foundation pile. Flow fields around rotating turbines are more complicated compared to the flow conditions of scour around a single pile. To validate the scour depth around turbines with rotating rotors, the topographic contours of the local scour shapes after scour equilibrium of model T2 are shown in Figure 10. The model T2 was set up by adding the rotating rotor on the foundation pile, as shown in Figure 10a. In model T2, the horizontal length of the scour hole was 2.5D, 4D, and 3D upstream, downstream, and on the side of the pile, respectively. The horizontal extent of the scour hole was a little larger than the scour at the single pile. The maximum scour depth was 7.8 cm, which was about 50% deeper than the scour at the single pile in model T1. This result reflects the impact of turbine rotors on seabed scour. This phenomenon of deeper scour depth induced by turbine scour than scour at piles was also mentioned in a turbine scour experiment, as discussed in Section 3.3. The experimental data also showed a slightly larger than 50% deeper scour hole when the turbine was installed 0.5H above the seabed. In summary, from the simulation results, all of the topographic bed features observed were captured: (i) A semicircular shape (in plan view) of the upstream part of the scour hole with a slope equal to the angle of repose; (ii) The formation of a gentler slope on the downstream side of the scour hole; (iii) A maximum scour depth occurring on the near sides of the foundation pile. The maximum scour depth simulated by the proposed turbine scour model had good agreement with the experimental data. The accuracy of the proposed scour simulation should be acceptable from the viewpoint of scour depth. However, note that the maximum scour depth calculated by the present numerical model was about 15% smaller than the experimental data, since unsteady effects were neglected. Hence, the maximum scour depth of full-scale turbines as discussed in Section 4 may be conservative.

Model Validation through Second Group
This section shows the equilibrium scour profile of test models in group 2. After a comparison to a scale-down turbine scour experiment, the proposed numerical scour model was approved as effective and accurate to calculate scale-up models.

Scour Depth Validation
In model group 2, the main variable was tip clearance, from C/H = 0.25 to 1.0, since tip clearance is one of the main factors influencing bed velocity, as investigated in References [11] and [35]. The topographic contours of equilibrium scour shapes of model group 2 are shown in Figure 11. It can be seen that the maximum scour depth occurred at the sides of the foundation. Generally, a lower tip clearance causes a deeper and wider scour hole. A comparison of maximum scour depths between the experiments and CFD simulations is shown in Figure 12. The scour depth St is nondimensionalized by pile diameter D. The errors between the CFD and experimental results were under 15%. Big changes in scour depth between cases of C/H = 0.5 and C/H = 0.75 occurred in the experimental and numerical results. This may have been due to the installation of a turbine rotor high enough off the seabed. The scour hole was mainly induced by accelerated flow in between the rotor and seabed. However, when the turbine was installed close to the seabed (C/H < 0.5), the flow suppression effects due to the existence of the turbine rotor and the increase in turbulence intensity around the tip regions jointly affected the increase in the scour level. Hence, the scour depth showed a big jump from C/H = 0.75 to C/H = 0.5.  Figure 13, where D is the diameter of the turbine foundation pile. In Figure 13, flow velocity U is divided by inlet velocity V on the x axis, and vertical height Z is divided by turbine height H on the y axis. It can be seen that the flow velocity decays in the region adjacent to the blade height of the rotor were due to energy extraction, and increases in the region immediately above and below the turbine were due to blockage effects. Blockage effects occur due to the incoming flow diverting the obstructed fluid to the upper and lower sides of the rotor. The diverted flow accelerates around the rotor to balance the flow mass conservation. The acceleration ratio of the flow velocity at slipstream between the turbine rotor and the seabed was about 10%-20%. In slipstream, the maximum flow velocity occurs close to the tip of the rotor. The faster flow rate approaching the seabed can be reached at a lower tip clearance. This is one of the main reasons why scour holes are deeper at a lower tip clearance.

Scour Profile Validation
The topographic contours of turbine scour shapes at C/H=0.5 were chosen to validate the simulation scour profile results. As seen in Figure 14, an outline of a scour hole is illustrated with a red line. The scour hole showed a similar outline: a semicircle upstream and a semicone downstream of the pile. However, the simulated scour width and longitudinal length were relatively lower than in the experimental data. The scour width was 3D in the simulation and 4D in the experiment, and the longitudinal length downstream was 5.5D in the simulation and 8D in the experiment. The proposed scour model could not simulate the sand dune formed behind the turbine.

Error Analysis of the Proposed Model
The quantitative comparisons between the proposed CFD model and the experimental data are shown in this section. The simulated tidal turbine scour depth was compared to the experimental measurements, as shown in Figure 15. The scour depth was normalized by the pile's diameter. A maximum variation of 18% occurred in the results from Roulund's experiment [12]. Generally, the proposed CFD model underestimated the seabed scour depth of tidal current turbines under an 18% error range. Based on the experimental results from Melville [33], Roulund et al. [12], and the current tests of turbine scour, statistical functions using the correlation coefficient (R), mean absolute error (MAE), root mean squared error (RSME), and scatter index (SI) (Equations (12)-(15)) were calculated to evaluate the accuracy of the proposed model.
where is the experimental value, is the simulated value, is the average experimental value, is the average simulated value, and M is the amount of experimental data. The specific value of error analysis of the proposed model is listed in Table 4. It can be seen that the maximum scour depth predicted by the proposed model fit well with the experimental data with a little underestimation.  In conclusion, the proposed numerical model could simulate the scour process and final scour hole profiles around the turbine with reasonable agreement. The maximum scour depth was a little lower than in the experiment, since unsteady effects were neglected with an unsteady effect. The hydrodynamic model used the RANS model to simulate the flow field and then calculate the seabed shear stresses as the initial input for the morphologic model. In the morphologic model, the seabed was updated at each time step using the dynamic meshing method. This model could capture the main factors for turbine scour and simulate the scour process with a low computational cost and an effective calculated performance. However, this model had its own shortcomings. It could not consider unsteady factors such as vortex systems. In the proposed CFD model, the calculation of sediment transport only considers steady-flow conditions. The hydrodynamic model uses turbulence intensity to express the oscillating flow components. The oscillation flow components are not brought into the calculation of seabed scour in the morphological model. The vortex system can promote turbine scour in reality. Therefore, the scour depth calculated by the CFD model reached an equilibrium depth faster than in the experiment. The scour depth was also more shallow in the CFD results than in the experiment.
On the other hand, the proposed model neglected the suspended sediments. This is because the suspended sediments showed little impact on local scour around the foundation pile. Local scour around the turbine is mainly caused by downflow and a horseshoe vortex around the pile. The suspended sediments can be taken away by tidal currents and deposited downstream.

Numerical Investigation of Full-Scale Turbine Scour
A numerical investigation of the proposed full-scale turbine, in terms of the temporal and spatial development of seabed scour, was undertaken using the developed computational methodology and models. The dimensions of the proposed full-scale turbine are shown in Table 5.

Maximum Scour Depth
Maximum scour depth is one of the main factors for turbine scour prediction and protection. The equilibrium scour hole around a full-scale turbine was simulated using the proposed model, as shown in Figure 16. Figure 16 shows the topographic contours of the local scour shapes. The maximum scour depth was 1.96D, which was about 11% less than the scale-down experiment of C/H = 0.5. To further discuss the maximum scour depth and scour hole profiles, a comparison of scourhole profiles along the upstream-downstream symmetry line of different types of scour is shown in Figure 17. The experimental and CFD data of scour at piers measured by Roulund et al. [12] and the experimental data measured by Zhang et al. [9] were selected for comparison. It can be seen that the simulated maximum scour depth (1.96D) was 27% deeper than horizontal-axis turbine scour depth (1.42D) and 34% deeper than the experimental data of scour depth around a single pile (1.3D). The results show that scour depth around tidal current turbines was deeper than scour at piles. This was due to the flow acceleration at slipstream between the rotor and seabed. Furthermore, scour induced by Darrieus turbines is also deeper than scour induced by horizontal-axis turbines with the same flow conditions. This is because Darrieus turbines have more flow-contract block degrees with bigger cross-sectional areas above the seabed. Hence, the flow between Darrieus-type turbines and the seabed has a higher speed compared to a horizontal-axis turbine.
Both CFD results of scour depth around piles and scour around tidal current turbines showed about 15% less than scour depth from the experimental data. It should be noticed that real scour depth around full-scale turbines may be more than 1.96D, as simulated by the proposed model.

Horizontal Extent of Scour Hole
The horizontal extent of a scour hole is also important for the stability of marine structures. A comparison of scour hole edges between Darrieus-type turbines and horizontal-axis turbines is shown in Figure 18. It can be seen that the scour hole width of the full-scale turbine was about 4D. The scour hole length upstream was about 6D, and the scour hole length downstream was about 6D. As a comparison, the scour hole width of the scale-down experiment was 4D. The hole lengths upstream and downstream were 4D and 8D, respectively. As discussed in Section 3.3.2, the smaller scour hole length downstream was obtained due to the unsteady effects of the wake vortex being neglected in the numerical model. However, the simulated scour hole had a bigger size upstream, interestingly. The turbine diameter was 27 m in the full-scale turbine model for this case. The influence of the rotor-swept area covered up to a 6D distance around the turbine foundation. Hence, the bigger scour hole around the full-scale turbine was greatly affected by flow suppression of the rotating rotor. Hahn et al. [36] once investigated clear-water scour due to short vertical (pressure flow) contraction in a laboratory channel. Their results showed that a scour hole was caused right below an obstacle. In our model, the operating rotor acted like an obstacle above the seabed. Hence, the turbine scour is a composition of local scour at the foundation pile and suppression scour below the rotor. This can explain why the horizontal extent of the scour hole of the simulated full-scale turbine model was more than in other tests.

Temporal Development of Scour Depth
The temporal development of scour depth of the full-scale turbine is studied in this section. As shown in Figure 19, the scour process can be classified in three stages: (1) an initial stage, where the scour depth develops quickly in the first 5 mins due to the action of downflow and a horseshoe vortex; (2) a developing stage, where after the initial stage, the scour depth develops with gradually diminishing speed until equilibrium; and (3) an equilibrium stage, where the scour depth no longer increases evidently and remains at an almost constant depth. For the proposed full-scale turbine model, the seabed scour reached equilibrium after about 30 mins, and the final depth was 1.96D. However, the equilibrium scour time and final scour depth were less than in the turbine scour experiment. This is the deficiency of the CFD model and also appeared in Roulund's CFD model [12] for scour at piers. The RANS turbulence model cannot capture the vortex at a high Reynolds number. It uses turbulence intensity to represent the flow's instability. However, this unstable factor is not used for the calculation of seabed shear stress. Therefore, our morphologic model could not consider these unsteady factors. However, the vortex system can promote turbine scour in reality. The scour depth calculated by the CFD model reached equilibrium depth faster than in the experiment. The scour depth was also more shallow in the CFD results than in the experiment.

Turbine Scour Mechanism
The turbine scour mechanism is further explained in this section based on the aforementioned discussion. Generally, turbine scour starts with an increase of flow speed around the turbine supporting pile. Large bed shear stress on both sides of the pile is caused by fluid acceleration, and sediment is washed away. The removal of sediment results in shallow holes around the piles. The energy of the horseshoe vortex is enhanced, along with growth of the scour hole, which in turns expands the scour hole. As the scour hole deepens, a weaker horseshoe vortex and lower bed shear stress exist inside the scour hole. Finally, the flow energy is insufficient to reach critical shear stress to transport the sediment.
This phenomenon is much like scour at piles. However, the existence of a rotating turbine provides a more challenging and complicated mechanism, with blockage effects and disturbance. Blockage effects lead to flow suppression passing through the narrow region in between the tip and the seabed, and turbine disturbance is due to tip vortices. This can increase the horizontal and vertical extent of a scour hole. This phenomenon was discussed in Sections 3 and 4. Our main findings of a turbine's impact on seabed scour can be concluded as follows.
First, the rotor obstructs the incoming flow. The obstructed fluid moves to the upper and lower sides of the rotor and accelerates flow velocity around the rotor based on flow mass conservation. On the basis of a forerunner's study, the maximum scour depth increases with the averaged incoming flow velocity in clear-water scour conditions [6]. To discuss the contract effect of a turbine's rotor, the axial velocity of different layers between the turbine and the seabed of a full-scale turbine model is shown in Figure 20. The axial velocity shows about a 20% increment at location -4D along the x axis. The axial velocity near the seabed is less than the initial velocity due to the boundary layer effect. However, the fluid at each layer below the turbine has a higher axial velocity than at the same level before the turbine. Chen and Lam [11] once studied the axial velocity at slipstream between a horizontal-axis tidal current turbine and the seabed, and their results showed that axial velocity increased about 5% of initial flow, which is less than the Darrieus-type turbine in this project. This suggests that the impact of Darrieus-type turbines on scour hole development is greater than that of horizontal-axis turbines. Second, the impact of a turbine rotor on scour depth is more intense with lower tip clearance. As seen in Figures 13 and 14, the maximum scour depth was 1.45D, 1.48D, 1.92D, and 1.98D for C/H = 1.0, 0.75, 0.5, and 0.25, respectively (simulated by the CFD model). The acceleration ratio of flow velocity at slipstream between the turbine rotor and the seabed was about 10%-20%. The faster flow rate approaching the seabed could be reached at a lower tip clearance. This is one of the main reasons why a scour hole is deeper at a lower tip clearance.
Third, the proposed turbine scour model underestimated the final scour depth at rates of about 15%-20%, as discussed in Section 3. This is because downflow and the vortex systems are two main reasons for local scour around foundation piles. Our model could simulate the flow velocity change around a rotating turbine and at slipstream below the turbine, but transient vortex shedding cannot be simulated by RANS models. On the basis of the comparison between the CFD model and the experimental data, we found that about 80% of turbine scour depth was caused by the change in flow velocity, and unsteady effects such as the fluctuating components of a horseshoe vortex and lee-wake vortex flows occupied about 20% of the turbine scour's impact factors.

Conclusion
A study of seabed scour around a Darrieus-type tidal current turbine with a commercial CFD code for a hydrodynamic model and a user-defined function for a morphologic model was carried out. The obtained simulation results were in good agreement with experiments. After validation, the seabed scour of full-scale tidal current turbines was simulated using the proposed CFD model. The equilibrium scour depth, scour profiles, and temporal development of scour depth were investigated. The conclusions are the following: 1). The maximum scour depth of the proposed full-scale turbine was 1.96D, and the edge of the scour hole was 4D along the y axis, -6D upstream, and 6D downstream (results of the numerical simulation). Compared to scale-down experiments, the CFD result of scour depth was about 15% smaller and the scour hole length downstream was shorter due to the neglect of unsteady factors, such as the horseshoe vortex and lee-wake vortices; 2). The scour hole upstream was much wider than in the experiment. For full-scale turbines, scour is a composition of local scour at the foundation pile and suppressed flow-induced scour below the rotor caused by the contraction effect of the rotor. This situation causes a greater horizontal and vertical extent of the scour hole; 3). The temporal development of scour depth showed that a quick scour occurred at the initial stage, and equilibrium conditions were reached after about 30 mins of simulation. This time was a little shorter than in the experiment. In addition, the scour hole was 27% deeper than the horizontal axis turbine scour depth of 1.42D and 34% deeper than in the experimental data of scour depth around piles of 1.3D after equilibrium.
This project first proposed a numerical model to simulate full-scale Darrieus-type tidal current turbines. It was proven that the proposed simulation method is a good and reliable tool to predict and discuss the spatial development of the scour process induced by tidal current turbines. It must be noted that the CFD model produces relatively conservative results of scour depth. The correction of existing empirical models of scour depth around Darrieus-type tidal current turbines by using the proposed numerical model and a flow analysis will be the next emphasis in future work.