Effects of the Current Direction on the Energy Production of a Tidal Farm: The Case of Raz Blanchard (France)

: This study aims to investigate the influence of the current direction on the energy production of a tidal turbines array. It is based on a three-dimensional (3D) numerical simulation of the flow where the turbines are represented with actuator disks. The case study consists of modelling the energy extraction of a small array of turbines (staggered and aligned layouts) placed in the Raz Blanchard (Alderney Race, France). The simulations are performed with hydrodynamic data (current magnitude and direction) representative of a mean tide, with several resistance forces and ambient turbulence intensities. The influence of the current direction on the energy production is highlighted by comparing the simulations forced with the real current direction with those in which the angle of incidence between the incoming flow and the turbine’s axis is “switched off” (bi-directional flow). When the flow is aligned with the turbines’ axis (misalignment “switched off”), the staggered layout produces more than the aligned arrangement. Comparison of the two types of simulations (misalignment switched off or not) shows that the misalignment of the flow around a predominant direction reduces the energy produced by the staggered layout and increases the production of the aligned layout. Furthermore, it suggests that the mean energy produced per machine is almost the same for both layouts. Higher turbulence intensity reduces the positive effect of the directional spreading on the aligned layout production and limits the negative effect on the staggered layout production.


Introduction
Tidal energy is seen as a promising renewable energy. In Europe, resources of tidal energy are considerable. The Raz Blanchard (Alderney Race) site, which is located between the Alderney Island and La Hague Cape (Normandy, France), capitalizes about half of the tidal current resource estimated in France [1]. This high energy potential contributes to the attractiveness of the site for industrial projects.
However, the tidal energy market is still at a nascent status and it yields a number of unknowns with regards to different aspects of industrial projects in this sector. Tidal energy research contributes to overcome a number of challenges to prove the reliability of this technology and to provide developers key information to optimize device design and to select the locations where the turbines will be deployed in order to improve power extraction. For example, [2] and [3] study, respectively, the influence of wake effect on the electrical dynamics of tidal farms, the cost assessment, and the economic viability of tidal farm development in the Raz Blanchard.

Model Presentation
The turbine representation is based on the introduction of a force exerted by a porous disk on the fluid. It is introduced as a source term S i in the RANS (Equations 1-3): ∂U i ∂x i = 0 (Einstein s notation) (2) where U i is the i th component of the mean velocity, x i is the i th component of the location vector, µ is the water dynamic viscosity, µ t is the turbulent dynamic viscosity, g i is the i th component of the gravitational acceleration, R ij is the Reynolds stress tensor, p is the pressure, ρ is the fluid density, and δij is the Kronecker symbol. The thrust force Ft representing the total action of the disk on the flow is linked to the discontinuity in pressure ∆P induced by the disk. This force can be calculated with Equation (4) according to the AD theory [19,33]. The equation involves a thrust coefficient C t , which depends on the turbine's characteristics. According to the Betz limit, a value of 8/9 for C t corresponds to the maximum power coefficient value C p =16/27 [35,36].
where S is the swept area of the blades, and U ∞ is the upstream flow velocity. However, the definition of an upstream (unperturbed) velocity U ∞ is not straightforward, especially when a turbine is placed in the wake of another turbine. We, thus, opted for the formulation proposed by Taylor (1963) [37], which estimates the thrust from a local velocity U d and a resistance coefficient K (5-7). Equation (8) describes the sink term representing the thrust force in the RANS equations.
(1 + 0.25K) 2 (7) where e is thickness of the disk. The flow is resolved with the commercial code ANSYS FLUENT 14.5. The AD formulation is added to the code using a "user-defined-function". With regards to the turbulence, we opt for the standard k-ε model (Launder and Spalding). This choice is based on an inter-comparison between several turbulence models [38], which approves that the k-ε model does not require any correction on the turbulence energy equation to give accurate velocity and turbulent intensity in the wake of a porous disk.
The inlet velocity profile is of importance as it was demonstrated by Lewis et al. (2017) [39]. Those authors characterized the effect of velocity profile on tidal stream area depending on the nature of the bottom. In our model, the inlet velocity profile should account for this effect. The profile is described in the next subsections.

Model Validation: Flow Past a Porous Disk
The model results are now validated against the experimental measurements of Harrisson et al., 2010 [40] providing measurements of the turbulence intensity and the velocity behind a non-rotating porous disk representing a tidal turbine. The experiments are carried out in a circulating current flume 21 m in length and 1.35 m in width, with a water depth of 0.3 m. The porous disk representing the turbine has a diameter D of 0.1 m and a thickness e of 0.001 m. The disk is placed at mid-depth and 20 D far from the inlet. Measurements are performed using acoustic doppler velocimetry (Myers and Bahaj [35]). Several resistance coefficient K, corresponding to several disk porosity, have been used.
The dimensions of the calculation fluid domain were defined according to the circulating water channel used by Myers and Bahaj [35] and Harrisson et al., 2010 [40] to perform their experiments. The calculation domain is 5 m long and 1.36 m wide. The domain is meshed with 16,867,840 hexahedral elements. The mesh is refined in the disk area where the cell size equals the disk thickness. Symmetry conditions are used on the lateral and top boundaries. The pressure is imposed at the outlet. At the inlet of the domain (Figure 1), the velocity U I (z) profile and the turbulent kinetic energy k(z) profile are computed with Equations (9)- (12) where U* is the friction velocity, ν is the kinematic viscosity and z 0 is the bed roughness, κ is the eddy dissipation, κ = 0.41 is the Von Karman constant, and I(z) is the turbulence intensity profile (note that, to evaluate the level of turbulence in the wake of the AD, the turbulence intensity I is calculated by using the local maximal velocity to the dimensionless turbulent kinetic energy). I H is the turbulent intensity at z = H with H the water depth, l = 0.07 H is the vortex characteristic size, and C µ = 0.09. In Equation (9), considering a smooth bottom, z 0 is assimilated to the thickness of the viscous sublayer and the three parameters velocity U*, z 0 , and χ are evaluated by Equation (13) where U in is the depth-averaged inlet velocity.
For the model validation, the depth-averaged velocity is set to 0.3 m/s and the ambient turbulent intensity is set to 5% (in accordance with the measurements). Thus, the friction velocity is U* = 7.87 10 -3 ms -1 , and z 0 = 2.5 10 -8 m is the bed roughness. Special attention was paid to the mesh size in the alignment of the disk. According to Roc et al., in 2013 [32], the results are independent of the mesh size as far the mesh size in the y and z direction are such that ∆y ≤ D/9 and ∆z ≤ D/6. In the present simulations, we opt for a much smaller mesh size, signified by ∆y = D/100 and ∆z = D/100. The residual convergence criteria for velocity and turbulence are 10 -7 . intensity is set to 5% (in accordance with the measurements). Thus, the friction velocity is U* = 7.87 10 -3 ms -1 , and z0 = 2.5 10 -8 m is the bed roughness. Special attention was paid to the mesh size in the alignment of the disk. According to Roc et al., in 2013 [32], the results are independent of the mesh size as far the mesh size in the y and z direction are such that ∆y ≤ D/9 and ∆z ≤ D/6. In the present simulations, we opt for a much smaller mesh size, signified by ∆y = D/100 and ∆z = D/100. The residual convergence criteria for velocity and turbulence are 10 -7 . Three different values of Ct are considered (Table 1). For a flow without blockage effects, the Betz limit is reached when Ct = 0.89, which corresponds to a resistance coefficient K of 2. Numerical results obtained with this value of K have been compared to the experimental values [40] measured for Ct =0.86 (Table 1). As noticed in [38], the model predictions fit the measurements either in terms of velocity or turbulence intensity. Simulations with K = 1 and K = 3 also show good velocity predictions ( Figure 2). The greatest model discrepancies are located in the near wake. It is probably due to an overproduction of the turbulent kinetic energy in the AD [38], which accelerates the mixing effect in the near wake. This overproduction seems to be more remarkable in the case of low resistance case K = 1. Noteworthy, for the Alderney Race application, is that the longitudinal spacing between two consecutive lines of turbines is 7.5D. At this distance from the disk, the satisfactory agreement between model and experimental data is observed either in terms of flow velocity or turbulence intensity ( Figure 2). Table 1. Thrust and power coefficients for three resistance coefficients K studied in this paper. The validation of these results is based on a comparison with experimental results that Harrisson [40] obtained for porous disk, with Ct close to the values usually obtained for the Betz Limit. Three different values of C t are considered (Table 1). For a flow without blockage effects, the Betz limit is reached when C t = 0.89, which corresponds to a resistance coefficient K of 2. Numerical results obtained with this value of K have been compared to the experimental values [40] measured for C t =0.86 ( Table 1). As noticed in [38], the model predictions fit the measurements either in terms of velocity or turbulence intensity. Simulations with K = 1 and K = 3 also show good velocity predictions ( Figure 2). The greatest model discrepancies are located in the near wake. It is probably due to an overproduction of the turbulent kinetic energy in the AD [38], which accelerates the mixing effect in the near wake. This overproduction seems to be more remarkable in the case of low resistance case K = 1. Noteworthy, for the Alderney Race application, is that the longitudinal spacing between two consecutive lines of turbines is 7.5D. At this distance from the disk, the satisfactory agreement between model and experimental data is observed either in terms of flow velocity or turbulence intensity ( Figure 2). Table 1. Thrust and power coefficients for three resistance coefficients K studied in this paper. The validation of these results is based on a comparison with experimental results that Harrisson [40] obtained for porous disk, with C t close to the values usually obtained for the Betz Limit.

The Present Numerical Simulations
Experimental Conditions   Table 1). Note that, to evaluate the turbulence intensity (IT) in the wake of the AD, Equation (10) was applied by using the local maximal velocity to dimensionless the turbulent kinetic energy.

Model Configuration: Wake-Field Study
The method is now applied to a pilot tidal farm located in the Raz Blanchard. Two layout configurations are considered ( More recently, Vennell et al. [41] in a review paper indicated that a minimal inter-row spacing of 5 D is recommended. Finally, Bai et al. [20] demonstrate that the appropriate lateral distance, according to the power extraction, is approximately 2.5 D, whereas the minimum inter-row spacing should be 6 D. For the present purpose, the tidal farm is constituted of three lines of turbines separated by a longitudinal distance b = 7.5 D. The lateral spacing between turbines, from center to center, is set to a = 3 D. Garrett and Cummins [42] theoretically derived the limit of power extraction from the flow by using the so-called linear momentum actuator disk theory (1D-LMADT). Nishino and Willden [43] compared 1D-LMADT to 3D RANS AD computation results obtained for a single machine in a channel and they also applied a similar approach [44,45] to a partial tidal fence. Considering the definition of the local blockage, BL (the ratio between a single device area and the local passage crosssectional area) and the global blockage, BG (the ratio between the total device area and the channel cross-sectional area), proposed in Garrett and Cummins [42] and in Nishino and Willden [44,45], the blockage values are BL = 0.1 and BG = 0.03. According to Garret and Cummins [42], a local blockage of 0.1 should give a maximum power coefficient of 0.72 by using the relationship Cpmax=16/27 (1-BL) -2 for a single device (which is beyond the limit of Betz, Cpmax =16/27).  Table 1). Note that, to evaluate the turbulence intensity (IT) in the wake of the AD, Equation (10) was applied by using the local maximal velocity to dimensionless the turbulent kinetic energy.

Model Configuration: Wake-Field Study
The method is now applied to a pilot tidal farm located in the Raz Blanchard. Two layout configurations are considered (  [41] in a review paper indicated that a minimal inter-row spacing of 5 D is recommended. Finally, Bai et al. [20] demonstrate that the appropriate lateral distance, according to the power extraction, is approximately 2.5 D, whereas the minimum inter-row spacing should be 6 D. For the present purpose, the tidal farm is constituted of three lines of turbines separated by a longitudinal distance b = 7.5 D. The lateral spacing between turbines, from center to center, is set to a = 3 D. Garrett and Cummins [42] theoretically derived the limit of power extraction from the flow by using the so-called linear momentum actuator disk theory (1D-LMADT). Nishino and Willden [43] compared 1D-LMADT to 3D RANS AD computation results obtained for a single machine in a channel and they also applied a similar approach [44,45] to a partial tidal fence. Considering the definition of the local blockage, B L (the ratio between a single device area and the local passage cross-sectional area) and the global blockage, B G (the ratio between the total device area and the channel cross-sectional area), proposed in Garrett and Cummins [42] and in Nishino and Willden [44,45] The turbine arrangement is symmetric with respect to x-and y-axes. The first line is situated 20 D behind the inlet (left). The domain is covered by a non-uniform Cartesian mesh containing 10,904,400 elements (44,391,533 nodes) with a refinement in the vicinity of the machines. At the bottom, a no-slip condition is used whereas a slip condition is applied at the top. The ratio of disk thickness over mesh size ratio is 2.5 in the x-direction and 1.25 in the y-and the z-directions.
For the present purpose, the wake-field model (Fluent model) is forced with tidal conditions representative of the Raz Blanchard. Those conditions have been extracted (at the selected point) from a regional model named Telemac2D. This model solves the shallow non-linear equations using a finite element formulation [46]. The Telemac2D results (water levels and depth-averaged current velocities) have been validated with the in-situ measurements of [47]. The Telemac2D model indicates that the annual power density, defined as 1/2ρU 3 (where U is the depth averaged velocity magnitude), is 5.22 kW/m 2 at the study site. The period retained for this study corresponds to a mean tide. Results have been extracted between 19:21 on 24 January 2005 and 07:28 on 25 January 2005. The model regional model shows that the local hydrodynamics is influenced by a large eddy that appears north of La Hague Cape, as shown in Figure 4. When the eddy develops (during flood tide), the current velocities vary in direction, which results in a large directional spreading. The current rose representation ( Figure 5) highlights a significant tidal asymmetry (in both magnitude and direction) between ebb and flood. The temporal evolution of the flow velocity and the water depth during the study period is represented in Figure 6. Those data have been used to build the vertical profiles of current velocities imposed at the boundaries of the local (fluent) model. To this end, we used Equations (9)- (12) and (14), which apply for rough seabed.
where R is the bed roughness.
During the studied period, the water depth at the extraction point varies between 52.92 and 57.15 m. Sensitivity tests indicate that the water depth has a negligible effect on the total production (<1%). We, therefore, use a unique depth value for all simulations (H = 55 m).
At the extraction point, the optimal direction (the direction maximizing the annual power extraction) was 64° with respect to the north direction. This direction was used to orientate the frame of the wake-field model (and to compute the misalignment with respect to the predominant direction). Then, periods when velocity magnitudes were smaller than 1 m/s were neglected, and the others were discretized temporally. Finally, we obtained 20 tidal conditions representative of the studied mean tide ( Table 2). In Table 2, the regimes 1 and 2 correspond to the flood and the ebb, respectively. Those 20 tidal conditions were then used to force the wake-field model. We, thus, The turbine arrangement is symmetric with respect to x-and y-axes. The first line is situated 20 D behind the inlet (left). The domain is covered by a non-uniform Cartesian mesh containing 10,904,400 elements (44,391,533 nodes) with a refinement in the vicinity of the machines. At the bottom, a no-slip condition is used whereas a slip condition is applied at the top. The ratio of disk thickness over mesh size ratio is 2.5 in the x-direction and 1.25 in the y-and the z-directions.
For the present purpose, the wake-field model (Fluent model) is forced with tidal conditions representative of the Raz Blanchard. Those conditions have been extracted (at the selected point) from a regional model named Telemac2D. This model solves the shallow non-linear equations using a finite element formulation [46]. The Telemac2D results (water levels and depth-averaged current velocities) have been validated with the in-situ measurements of [47]. The Telemac2D model indicates that the annual power density, defined as 1/2ρU 3 (where U is the depth averaged velocity magnitude), is 5.22 kW/m 2 at the study site. The period retained for this study corresponds to a mean tide. Results have been extracted between 19:21 on 24 January 2005 and 07:28 on 25 January 2005. The model regional model shows that the local hydrodynamics is influenced by a large eddy that appears north of La Hague Cape, as shown in Figure 4. When the eddy develops (during flood tide), the current velocities vary in direction, which results in a large directional spreading. The current rose representation ( Figure 5) highlights a significant tidal asymmetry (in both magnitude and direction) between ebb and flood. The temporal evolution of the flow velocity and the water depth during the study period is represented in Figure 6. Those data have been used to build the vertical profiles of current velocities imposed at the boundaries of the local (fluent) model. To this end, we used Equations (9)- (12) and (14), which apply for rough seabed.
where R is the bed roughness.
During the studied period, the water depth at the extraction point varies between 52.92 and 57.15 m. Sensitivity tests indicate that the water depth has a negligible effect on the total production (<1%). We, therefore, use a unique depth value for all simulations (H = 55 m).
At the extraction point, the optimal direction (the direction maximizing the annual power extraction) was 64 • with respect to the north direction. This direction was used to orientate the frame of the wake-field model (and to compute the misalignment with respect to the predominant direction). Then, periods when velocity magnitudes were smaller than 1 m/s were neglected, and the others were discretized temporally. Finally, we obtained 20 tidal conditions representative of the studied mean tide (Table 2). In Table 2, the regimes 1 and 2 correspond to the flood and the ebb, respectively. Those 20 tidal conditions were then used to force the wake-field model. We, thus, performed a unique non-stationary simulation with Telemac2D and 20 stationary simulations with the fluent model.        Table 2. Boundary conditions calculated with the regional model ( Figure 7). The angle is expressed in nautical convention. The incidence is given with respect to the optimal flow direction (the x axis of the local model frame). The power of each machine Pij is calculated with Equation (15) where Ud is the local velocity (along x-axis). The power Pi and the mean power Pim of a given line are given by Equation (16) Table 2. Boundary conditions calculated with the regional model ( Figure 7). The angle is expressed in nautical convention. The incidence is given with respect to the optimal flow direction (the x axis of the local model frame). The power of each machine P ij is calculated with Equation (15) where U d is the local velocity (along x-axis). The power P i and the mean power P im of a given line are given by Equation (16) and the total power of the tidal farm P farm is calculated with Equation (17). N i and N l are, respectively, the number of machines of the line i th and the total number of lines. During the flood (regime 1), the turbines of lines 2 and 3 are located downstream from line 1, whereas during the ebb (regime 2), the turbines of lines 1 and 2 are located downstream from the turbines of line 3.

Number of the Simulation
When turbines are placed in array, their wakes interact in different ways [48] and the inflow of a given turbine could be affected by one or several wakes. The most common wake interaction corresponds to a configuration where one turbine is in the wake of another. This interaction is named tandem operation. When two machines placed side by side interact and form a merged wake, the interaction is referred to as interference operation. When the inflow of a turbine is partly affected by one or more wakes, the interaction is named overlapping. To highlight different types of flow interactions occurring within the array, maps of velocity deficit (18) will be plotted.
For each layout, two scenarios are examined. The first one corresponds to the velocity and incidence angles introduced in Table 2. It will be referred as the realistic case with incidence angle (RWI). The second scenario consists of "switching off" the incidence angle (the current reverses twice per tide but remains parallel to the x-axis; it varies in magnitude). It is referred to as realistic case without incidence angle (RWOI).
The ratio between the averaged power production per line in the RWI and the RWOI cases is an indicator of the influence of the current incidence on the power. This ratio power indicator is noted RP =P im (RWI)/P im (RWOI). A value close to 1 indicates that the flow obliquity has a negligible influence on the production. A value smaller than 1 indicates that the flow obliquity reduces the production and vice versa. The energy produced by the farm during the tidal cycle is also computed as well as the energy per line. Finally, as the number of machines of the two layouts differs (9 devices in the aligned layout and 10 devices in the staggered layout), the mean energy per device has been computed, which permits an unbiased comparison.
Several values of resistance coefficient K and ambient turbulent intensity are used ( Table 3). The inlet formulation for I H =5% is related to Equation (13) whereas, for higher turbulence intensities, the relations for rough bed is used (Equation (14)).

Effect of Layout
First, the case with K = 2 and I = 5% has been considered. Figures 7 and 8 represent maps of velocity deficit (Equation (18)). The time evolution of the mean power per line (P im ) is represented in Figure 10. For both RWOI and RWI cases, three types of flow interactions (overlapping, interference, and tandem) can be observed (Figures 8 and 9).      For the A-RWOI case (aligned layout without incidence), machines of the three lines (M1j, M2j, M3j) are systematically in tandem situation (Figure 7 left). During the flood, the turbines of the two lines located downstream (lines 2 and 3) produce significantly less than the linear turbines upstream (line 1). During the ebb, the power extracted by the machines of the line 2 is lower than the power extracted by the machines of the line 3 (Figure 9c). This is due to a greater flow recovery between line 2 and line 3 than between lines 1 and 2. This higher flow recovery is imputed to a greater turbulence intensity as the turbulent intensity increases each time the flow passes through a line of turbines.
For the S-RWOI case (staggered layout without incidence), the turbines of lines 1 and 3 are placed in tandem and they are affected by the superposition of the wakes induced by the turbines of line 2  Figure 7, right). There is a significant difference in production between line 1 and line 3, which depends on the orientation of the current (considering if the line is located upstream or downstream). The turbines of the 2nd line are only slightly affected by the upstream turbines by wake overlapping. The mean power of the line 2 is similar to the one of the upstream line of machines (Figure 9b).
The staggered configuration is supposed to be advantageous compared to the aligned layout [49]. This is confirmed by our results when the flow obliquity is zero: The staggered layout leads to a greater energy production (10.57 MWh /per machine) compared to the aligned layout (8.64 MWh/per machine) ( Table 5). This is mainly due to the smaller production of the second line of machines of the aligned arrangement.
In the case of A-RWI case (Figure 8, left) (aligned layout with incidence), the overlapping configuration prevails, and the power extracted by the two downstream lines (Figure 9c) is higher than the one obtained without incidence (Figure 9d). The difference in production between the upstream line and the other lines is more significant during regime 2 than during regime 1.
With St-RWI case (Figure 8 right) (staggered layout with incidence), the wakes overlap, and the power of the second line is significantly reduced in comparison to the aligned layout case (Figure 9a). Occasionally, three machines are placed in tandem, which leads to a noticeable decrease in power (Figure 8, right bottom). By way of illustration, it occurs at t = 12.25 h when the incidence angle is close to 10 • or at t = 0 h (Figure 9f) when the incidence angle is 8.9 • . In such configurations, although the turbines are placed in a staggered layout, the flow obliquity is such that the turbines M13, M23, and M33 (or M12, M22, M31) are aligned.

Effect of Flow Obliquity
The influence of flow obliquity on the production is different when considering a staggered or an aligned layout (Figure 9). Whereas the flow obliquity increases the power of the aligned turbines, it reduces the power of the staggered turbines. It reduces the overall production of the staggered farm as well as the mean production per device ( Table 5). The energy reduction is −8.5% (from 105.72 MWh to 96.73 MWh) for the staggered layout and the energy increase is +11.3% for the aligned layout. Finally, when the current varies around the predominant direction (RWI), the gain in using a staggered layout (rather than an aligned layout) is very small (9.62 MWh/per machine for the aligned layout and 9.67 MWh/per machine for the staggered layout) ( Table 5). This indicates that the staggered arrangement is not, systematically, the most suitable layout.
The effect of the flow obliquity is complex and differs according to the type of layout, the line of machine, and the regime ( Figure 10). As an example, the production of line 1 is only slightly modified by the flow obliquity during regime 1 whereas the power of line 2 significantly decreases during the two regimes. To discriminate the effect of the flow obliquity during ebb and flood tides, the RP indicator is computed per regime for both aligned and the staggered layout (Table 4).
With regards to the line located upstream (line 1 during regime 1 and line 3 during regime 2), the value of RP is smaller than 1 (0.9 for line 1 during R1 and 0.97 for line 3 during R2). This means that the turbines produce less when the flow comes obliquely, which is due to a reduction of the normal velocity component (with respect to the devices' axis, an analogous calculation of P ij , where the x-axis component of U d is replaced by its magnitude, highlights the impact of the flow obliquity on the energy production of the aligned layout (K=2 -I=5%) (cf. Table 4). Indeed, for the RWOI case, we observe a change of 2.3% between to the two different energy calculations. Note that this gap increases until 7.8% for the RWI case. With respect to the 2nd line, the ratio RP shows opposite trends depending on the layout. For the staggered layout, the flow obliquity has a negative effect on the second line production (the ratio is 0.79 or 0.86) whereas the aligned layout creates a positive effect (the ratio is 1.17 or 1.22). For the aligned layout, the energy production is lower than the one of line 1 or 3 (Table 5A) regardless of whether the incidence is switched off or not. The effect of the flow obliquity is complex and differs according to the type of layout, the line of machine, and the regime ( Figure 10). As an example, the production of line 1 is only slightly modified by the flow obliquity during regime 1 whereas the power of line 2 significantly decreases during the two regimes. To discriminate the effect of the flow obliquity during ebb and flood tides, the RP indicator is computed per regime for both aligned and the staggered layout (Table 4).
With regards to the line located upstream (line 1 during regime 1 and line 3 during regime 2), the value of RP is smaller than 1 (0.9 for line 1 during R1 and 0.97 for line 3 during R2). This means that the turbines produce less when the flow comes obliquely, which is due to a reduction of the normal velocity component (with respect to the devices' axis, an analogous calculation of Pij, where the x-axis component of Ud is replaced by its magnitude, highlights the impact of the flow obliquity on the energy production of the aligned layout (K=2 -I=5%) (cf. Table 4). Indeed, for the RWOI case, we observe a change of 2.3% between to the two different energy calculations. Note that this gap increases until 7.8% for the RWI case. With respect to the 2nd line, the ratio RP shows opposite trends depending on the layout. For the staggered layout, the flow obliquity has a negative effect on the second line production (the ratio is 0.79 or 0.86) whereas the aligned layout creates a positive effect (the ratio is 1.17 or 1.22). For the aligned layout, the energy production is lower than the one of line 1 or 3 (Table 5A) regardless of whether the incidence is switched off or not.
Indeed, the turbines of line 2 are, most of the time, in overlapping or tandem configuration. When the current is parallel to the turbines' axis, the turbines of line 2 are in tandem, which is the worst configuration. When the flow arrives obliquely, the turbines of line 2 move from a tandem to an overlapping configuration (Figure 8, left) and, thus, produce slightly more (Figure 9c,d).
Concerning the S-RWI case, there could be a variety of wake interactions along line 2. The turbines can be located either in an unperturbed flow (the turbine located on the side (M21 or M24)) or in a funneled flow (turbines M22 and M23 benefit from the flow acceleration between the turbines of the line of turbines located upstream). They can also be in an overlapping or a tandem configuration. In the overlapping configuration (Figure 7, right), the power of line 2 is significantly reduced in   Table 5. Energy produced during the mean tide (MWh) for K = 1, 2, and 3 and t = 5% configurations. The energy is calculated either for each line or for all the turbines (farm). The column farm/N, where N is the number of machines of the farm, indicates the mean energy delivered per turbine. Energies are computed for the two scenarios (RWOI and RWI) and aligned and staggered layout. (*) The power captured by a machine is calculated by using the local velocity along x-axis (15). Replacing the local velocity along x by the local velocity magnitude puts to the fore the impact of the flow obliquity on the production of a tidal farm. Indeed, the turbines of line 2 are, most of the time, in overlapping or tandem configuration. When the current is parallel to the turbines' axis, the turbines of line 2 are in tandem, which is the worst configuration. When the flow arrives obliquely, the turbines of line 2 move from a tandem to an overlapping configuration (Figure 8, left) and, thus, produce slightly more (Figure 9c,d). Concerning the S-RWI case, there could be a variety of wake interactions along line 2. The turbines can be located either in an unperturbed flow (the turbine located on the side (M21 or M24)) or in a funneled flow (turbines M22 and M23 benefit from the flow acceleration between the turbines of the line of turbines located upstream). They can also be in an overlapping or a tandem configuration. In the overlapping configuration (Figure 7, right), the power of line 2 is significantly reduced in comparison to the aligned layout case (Figure 9a). Over the tide, the energy produced is equivalent to the one of the other lines 28 For the RWOI case, the machines placed downstream (line 3 during R1 and line 1 during R2) are in tandem configuration (aligned layout) or wake overlapping (staggered layout). Depending on the layout, the distance between two consecutive machines is b for the aligned layout and 2b for the staggered one. When the spacing between turbines is 2b, the flow recovery is nearly complete, and the turbines located downstream thus produces as much as the turbines located in the free stream. When the incidence is "switched off", the flow interaction between the wakes moves from overlapping/interference to tandem/overlapping. The flow obliquity favors the production of the aligned layout (1.22 for L3 and R1; 1.28 for L1 and R2) and reduces the production of the staggered layout (0.86 for L3 and R1; 0.97 for L1 and R2). Occasionally, three machines are placed in tandem, which leads to a significant decrease in power (Figure 8, right bottom). For instance, it occurs at t = 12.25 h when the incidence angle is close to 10 • or at t = 0 h (Figure 9f) when the incidence angle is 8.9 • . In such configurations, though the turbines are placed in a staggered layout, the flow obliquities in such that the turbines M13, M23, and M32 (or M12, M22, M31) are aligned.

Effect of Resistance Coefficient and Ambient Turbulence Level
The analysis has been extended to examine several values of the resistance coefficient and ambient turbulence level (Tables 4 and 5).
When K increases from 2 to 3, the total energy (Table 5) increase is about 40% and 48% for the A-RWOI and the S-RWOI cases, respectively. This increase is around 44% with the RWI configurations. The findings concerning the effect of the incidence obtained for the baseline case (I = 5% K = 2) are also applied to the configurations (I = 5% K = 1) and (I = 5% K = 3) (the obliquity increases the power produced by the aligned layout and reduces the power produced by the staggered layout; the mean energy produced by each machine is the same with both the aligned and the staggered layout when the current is multi-directional). However, it should be noted that the baseline trend is exacerbated with a greater value of the coefficient of resistance. The relative reduction of the energy between RWOI and RWI for the "staggered layout' increases from 8.5% when K = 2 to 9.9% when K = 3. For the aligned layout, the energy increase remains close to +11.3% for K = 1 and K = 2 and reaches +15.2% for K = 3.
The impact of the ambient turbulence intensity on the farm production is estimated by comparing simulations with different turbulence intensity values I = 5%, 10%, and 15%; here, a constant K receives the value of 2 (Table 3). It is observed that the energy produced (Table 5) decreases with turbulence intensity when it gets higher values whether the obliquity is "switched off" or not. This reduction is imputed to the modification of the vertical profile of x-velocity when the turbulent intensity increases ( Figure 11). The differences between cases RWI and RWOI are coherent with the previously described configurations. For the staggered configuration, the relative reduction of the energy between RWOI and the RWI cases is −13.2 for I = 10 and 15%. For the aligned layout, the gain in energy is +11.3 for I = 5% and +6.9% for I = 15%.

Effect of the Tidal Asymmetry
Similar to other tidal straits [7,8,10,17], the flow in the Raz Blanchard is characterized by an asymmetry of the tidal current amplitude as well as a directional spreading around a predominant direction. Depending on the current direction (ebb or flood), the turbines of lines 1 and 3 are alternately in the free stream or located downstream the array. The asymmetry of the current amplitude is 6%, which leads to a systematic difference between the averaged powers captured by L1 and L3. Moreover, the available kinetic energy fluctuates during the ebb and the flood. Our results highlight that, for both layouts, the current misalignment causes a reduction in the averaged power of about 10% for L1 (in the free stream during t during the flood tide). This decrease falls to 3% for L3 during the ebb. What is the cause of this disparity? At the studied site, the angles of incidence Similarl to the baseline case, the ratio RP is greater than 1 for the downstream lines of the aligned layout and lower than 1 for the downstream lines of the staggered layout. Nevertheless, the amplitude of the ratio for the aligned arrangement decreases with increasing values of I (Table 4). For example, considering line 2 of the aligned layout during the ebb (R2), RP varies from 1.17 in the case I = 5% to 1.02 with I = 15%. The same tendency is perceived for the downstream lines (L1 in regime R2 and L3 in regime R1). Considering the staggered arrangement, the ambient turbulence intensity does not significantly modify the ratio RP of the line 2. However, it significantly affects the downstream line of turbines (L1 in R2 and L3 in R1). For instance, the RP of line L1 in regime R2 varies from 1.01 for I = 5% to 0.87 for I = 15%). Noteworthily, the RP of the lines of turbines located in the free stream (L1 in R1 and L3 in R2) is unchanged when the turbulence intensity varies from I = 5% to 15%.

Effect of the Tidal Asymmetry
Similar to other tidal straits [7,8,10,17], the flow in the Raz Blanchard is characterized by an asymmetry of the tidal current amplitude as well as a directional spreading around a predominant direction. Depending on the current direction (ebb or flood), the turbines of lines 1 and 3 are alternately in the free stream or located downstream the array. The asymmetry of the current amplitude is 6%, which leads to a systematic difference between the averaged powers captured by L1 and L3. Moreover, the available kinetic energy fluctuates during the ebb and the flood. Our results highlight that, for both layouts, the current misalignment causes a reduction in the averaged power of about 10% for L1 (in the free stream during t during the flood tide). This decrease falls to 3% for L3 during the ebb. What is the cause of this disparity? At the studied site, the angles of incidence (with respect to the predominant direction) are relatively small (7 • on average during the studied period) during the ebb. However, during the flood, the presence of a large-scale eddy located northward of the Hague Cape deviates the current, which leads to mean and maximum angles of incidence of 12 • and 35 • (during the studied period). Indeed, as the flow obliquity, with respect to turbines' axis, increases, the axial projected area of the turbine decreases, and the kinetic energy extracted drops as well. These results are in accord with Galloway [50] and Frost's [10] considerations for a horizontal axial turbine. As a matter of fact, Frost [10] noticed a reduction of C p of 7% for an incidence angle of the incoming flow of 10 • and a value greater than 22% for an incidence angle of 20 • .
As earlier stated, the direction used to orientate the frame of the wake-field model maximizes the extraction of the hydrodynamic power. This orientation is determined from a one-year-long simulation and is, therefore, representative of a mean value. For the present study, we focus on particular tidal conditions (mean tide). Under those conditions, the mean current direction is not perfectly aligned with the yearly averaged predominant direction.

Conclusions
The influence of a time-varying current direction on the energy production of a pilot tidal farm has been analyzed using a quasi-steady RANS approach where the tidal turbines are represented as actuator disks. The interactions between wakes induced by adjacent turbines have been simulated considering two arrangements: Staggered (10 turbines) and aligned (9 turbines). The tidal farm model is forced using boundary conditions (velocity direction and magnitude) recommended by a regional model of the Alderney Race (France). Effects of the incidence between the incoming flow and the turbines' axis, the resistance coefficients, and the turbulence intensities are considered. The principal findings are listed below:

•
The flow obliquity influences, in a different way, the production of the staggered and the aligned arrays. Whereas it increases the production of the aligned layout, it reduces the production of the staggered layout. Mean energy produced per machine is almost the same for both layouts.

•
There is a significant difference between the production during the ebb and the flood tides.

•
Increasing turbulent intensity reduces the positive effect of the flow obliquity on the aligned layout production and restricts the negative effect on the staggered layout production.

•
The increase of ambient turbulence seems to slightly reduce the energy production of both layouts.
Additional simulations are therefore required to generalize the results and to improve understanding the influence of turbulence intensity influences on the energy production of a tidal farm. Our future works will consider the blockage effect, which is believed to have a strong influence on the flow dynamics in a pilot farm. Those simulations are realized with a quasi-steady model, that does not account for the Energies 2019, 12, 2478 18 of 20 rapid variation of direction and flow speed inerrant to turbulent area. Such a problem must be resolved with a finer model as presented by Bourgoin et al. [51] but the computation time is a limit to its use. The waves have an effect on the flow speed and the direction as noted by Lewis et al. [52]. Their effects on the energy production should be study further.