Next Article in Journal
Predictive Trajectory Control with Online MTPA Calculation and Minimization of the Inner Torque Ripple for Permanent-Magnet Synchronous Machines
Previous Article in Journal
Calorific Power Improvement of Wood by Heat Treatment and Its Relation to Chemical Composition
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Towards a Tidal Farm in Banks Strait, Tasmania: Influence of Tidal Array on Hydrodynamics

1
Australian Maritime College, University of Tasmania, Launceston, TAS 7250, Australia
2
School of Civil Engineering, University of Queensland, St Lucia, QLD 4072, Australia
*
Author to whom correspondence should be addressed.
Energies 2020, 13(20), 5326; https://doi.org/10.3390/en13205326
Submission received: 14 September 2020 / Revised: 1 October 2020 / Accepted: 1 October 2020 / Published: 13 October 2020
(This article belongs to the Section A3: Wind, Wave and Tidal Energy)

Abstract

:
The development of tidal energy in Australia is still a challenge with few studies performed on the characterisation of the resource, due to the difficulty to acquire data and uncertainties about the influence of this anthropogenic activity on the marine environment. Changes in flow could lead to alterations in sediment transport and have further influence on the marine habitat. A case study in a promising area, Banks Strait (Australia), was created using high resolution 2D and 3D models validated against in situ data to investigate changes to hydrodynamic conditions with two scenarios of tidal farms (100 and 300 turbines). Comparison between 2D and 3D is performed to find the best compromise between model accuracy and computational time for preliminary assessment. Changes to current speed and bed shear stress over a 35 day period were found to be localised around the tidal farms and did not extent more than 7 km from the farm (300 turbines) for both 2D and 3D. The results showed that for near field and far field, 2D models are sufficient to give a first approximation of the hydrodynamic influence of tidal farm deployment on its environment.

Graphical Abstract

1. Introduction

The development of tidal industry has made significant progress over the last decades: numerous tidal projects are under way, examples include development in United Kingdom (Pentland Firth, Shetland, Orkney), in France (Ushant, Aldernay Race) and in China (Zhoustan Archipelago) [1]. There also have been considerable advances in numerical models in the recent years, to assess the influence of tidal farms on the environment principally in 2D [2,3,4,5,6,7,8,9,10,11,12]. However, the conclusions of all these studies are site-specific and dependent on device design; thus, the effects of array size will need to be investigated in the site of interest. Generally, 2D models are used for regional modelling while 3D models are used to assess the local tidal energy resource. Nowadays researchers tend to use 3D models to assess the resource and the environmental impact, although at higher computational costs [13].
The objectives for 2030 for Australia (reducing CO 2 emissions to 26–28% below 2005 levels) and Tasmania’s plan to become 200% in renewables by 2040 are a major challenge. To meet these obligations, the Australian government is considering to speed up the deployment of ocean energy, looking towards blue economy projects to reduce fossil fuels reliance. Fourteen Australian tidal power trials have taken place over the last two decades, they are summarised in Table 1 [1,14]. To assess accurately tidal energy sites, the hydrodynamic conditions that produce effective tidal streams must be determined. Further high resolution projects are required to perform a detailed site assessments in order to evaluate the viable production, environmental influence of Australia’s tidal resource and tidal farm development.
The Australian Tidal Energy (AUSTEn) project [15] funded by the Australia Renewable Energy Agency (ARENA) aims to map Australian tidal energy resource with three components. One includes the characterization of two promising sites (the Banks Strait in Tasmania and the Clarence Strait in Northern Australia) with a focus on the energy resource potential, the influence on physical parameters and the development of fine-scale models. Banks Strait (Figure 1), a strait located in the north east of Tasmania between Swan Island and Clarke Island, has been identified as suitable for commercial tidal farm developments [16]. This 15 km wide strait connects the Tasman Sea to the south east of Bass strait: which due to its strong currents and proximity to grid-connection is a very attractive site for the development of the tidal arrays in Australia. Preliminary coarse resolution hydrodynamic modelling by the Commonwealth Scientific and Industrial Research Organisation (CSIRO) indicated strong currents of around 2.5 m/s in this area [17] and more recently [18] found an area of 9 km2 south of Clarke Island with maximum current speed higher than 2 m/s. The surveys performed by the AUSTEn project in Banks Strait in 2018 revealed important coastal features including the presence of large sand dunes near Clarke Island (Figure 1). The deployment of a tidal farm in this area may influence these natural phenomenons.
To assess a tidal energy site, the different steps are the collection of data with surveys, then the numerical modelling of the site to perform Environmental Impact Assessment (EIA). Previously in the Banks Strait region, field measurements such as bathymetry, current velocities, sub-bottom properties and suspended sediment concentration were sparse due to the energetic environment found in this area. Recently the AUSTEn campaigns filled this gap using a large set of instruments and surveys techniques. Most of studies conducted in this region [19,20,21,22,23,24,25,26] have focussed on the current flow through the Bass strait, without focussing specifically on Banks Strait region. This strait was not modelled with a resolution high enough to be reliable for accurate tidal energy site study until recently with new work found in the Banks Strait area: [18] used a 400 m mesh resolution model which was validated against an Acoustic Wave and Current Profiler (AWAC) field data with no implementation of turbines. More fine-scale modelling is therefore needed in Banks Strait in order to reduce the uncertainty of the influence of tidal energy extraction on the physical environment. To reduce uncertainty in sediment transport rates, hydrodynamic models should be calibrated to less than 10% errors in the velocity [27,28]. Recent literature and the comprehensive review of the current understanding of the hydro-environmental impacts [29] revealed the next challenges for tidal energy: more complex physics (3D), models coupling wave and current, very far field effect and one of the most important: field measurements, critical to validate model simulations.
The aim of this work was to investigate the influence of two tidal farms on physical parameters in Banks Strait using high-resolution 2D and 3D hydrodynamic ocean models. Highly calibrated unstructured mesh hydrodynamic models were validated against a large set of in situ data with implementation of two tidal farms (100 and 300 turbines). Investigation on whether 2D or 3D model dimensions are required for EIA in the near and far field are also presented. Details of the Banks Strait site and the AUSTEn surveys are presented in Section 2, with model validation results presented in Section 3 using field data from the AUSTEn campaigns. The influence in oceanographic systems of two simulated tidal farms is then described in Section 4. This study is also a solid basis for future work to determine the influence of tidal farms on sediment dynamics.

2. Modelling Methodology

This section outlines the description of the site, field data, numerical model, forcing inputs and validation process.

2.1. The Site and Data Collection

Banks Strait is located in the eastern part of Bass Strait between Clarke Island and the mainland of Tasmania. Principal harmonic constituents in the strait are the semidiurnal M2, S2 and diurnal O1, K1 with the semidiurnal M2 tide being the principal constituent [20,25]. The water mass characteristics of Banks Strait are altered by tidal stirring and local air–sea fluxes: with strong warming in summer and strong cooling in winter [26] and with saltiest water in winter [21]. Currents in the strait are primarily driven by tidal forcing, winds and shelf waves in Bass Strait. In Banks Strait, Sandery [21] described the area as well mixed because of the turbulent system due to strong tidal currents occurring. Sandery [20] also mentioned inferred flow in Banks Strait in winter with the sub-Antarctic Surface Water (SASW) intrusion.
The region of Banks Strait has been surveyed in March, July and December 2018, with the AUSTEn project [15], allowing measurements of seasonal variations. Field measurements at up to five stations using Acoustic Doppler Current Profiler (ADCP) were performed in Banks Strait in March, three in July and one in December (Figure 2). The details of ADCP deployments are described in Table 2. These instruments were configured to measure currents, to calibrate and validate the numerical model and were processed using the Integrated Marine Observing System (IMOS) Toolbox [30].
High-resolution bathymetry of the Strait was obtained using the Bathyswath v.1 (ITER system 234 kHz) between 16/03/2018 and 28/03/2018, with a resolution of approximately 5 m horizontal and 0.25 m vertical over an area of approximately 210 km2. The results of this survey showed depths varying between 10 m and 50 m with a deeper channel in the middle of Banks Strait with pockets of 70 m to 100 m. The multi-beam seafloor survey found large sand dunes that vary approximately between 5 and 10 m south-east of Clarke Island as shown in Figure 2. Field measurements of conductivity Temperature Depth (CTD) casts found vertical fluctuations of less than 5% [31], confirming the well mixed area in this energetic system and allowing the use of barotropic ocean models without any reduction in simulation accuracy. Data from the AWACS found average wave heights between 1 and 2 m with maximum values reaching over 5 m [32]. Studies of the region found currents higher than 2 m/s during ebb tides in the middle of the channel with a tidal asymmetry of 20 degrees between ebb and flood at the location of the station CW4 [18]. Maximum power density over 3 kW/m2 was found in the middle of the channel, exceeding the 2.5 KW/m2 threshold outlined for the economic viability of a site [33].

2.2. Model Domain and Forcing

High-resolution numerical models of Banks Strait were developed in 2D and 3D using the unstructured mesh MIKE21 and MIKE3 FM modelling software (DHI A/S, Hørsholm, Denmark). The software solves Reynolds-averaged Navier-Stokes (RANS) equations using the Boussinesq assumption [34]. The numerical domain extends between 38° S and 44° S and 141° E and 150° E. The mesh resolution varied in the domain from an average of 3 km at the boundary to 100 m at the approach of Banks Strait, in the strait the typical node spacing was 50 m and increased to 20 m in the area suitable for tidal farm installation and key interest features such as the sand dunes near Clarke Island (Figure 3). The model shoreline data were extracted from the Australian Bathymetry and Topography grid [35] and from Australian Hydrography Service AusEnc charts (AHS) [36] for better resolution in the Banks Strait area (north east coast of Tasmania and Furneaux Islands). Bathymetry data were also collected from Geoscience Australia (GA 2009) [35] for the larger domain, from the AHS for more precise features in the Banks Strait area (between the mainland of Tasmania and Flinders Island) and from the AUSTEn Project (multi-beam in March 2018). The collected bathymetry dataset was interpolated onto an unstructured mesh using the natural neighbour method.
Three open boundaries were forced with tidal constituents (surface elevation) from the Global tide MIKE model [37]. This tide model represents the major diurnal (S1, K1, O1, P1 and Q1), semidiurnal (M2, S2, N2 and K2) and shallow water/nonlinear (M4) constituents with a spatial resolution of 0.125° × 0.125° based on altimetry data from TOPEX/POSEIDON. The open boundaries were set far beyond Banks Strait to reduce their influence on simulation results as recommended [7,38]. The ERA5 climate reanalysis dataset [39] was used for wind components, with the collated dataset mean sea level and pressure closely matching the values from the Bureau of Meteorology measurements [40] in the Banks Strait at Swan Island ( R 2 = 0.9) when compared to the year 2018.

2.3. Model Settings

The barotropic models were run for 42 days to follow IEC recommendations [38], with 1 week spin up (to reach a steady state), from 16/03/2018 00:00 to 26/04/2018 00:00, generating a 35-day dataset. The Smagorinsky formulation was used for the horizontal eddy viscosity and the log-law formulation for vertical eddy viscosity. Changes in 3D compared to the 2D setup included the vertical layers and the bed resistance where the roughness height was used and calculated from Manning’s number (M). A Manning value of 32 was found to give the best calibration after sensitivity tests in 2D with the results shown in Section 3. This Manning value gave us a bottom drag coefficient of 0.0025 which represents a sand/gravel habitat [41]. Vertical dimension was discretised into 19 equidistant vertical layers.
Tidal farms with 100 and 300 turbines were tested to evaluate their influence on hydrodynamics. The tidal turbines were represented as a sub-grid object via an actuator disc model [34], implemented using Matlab scripts to automate the process [42]. The turbines were assumed to be aligned with the flow direction and were modelled using a fixed drag coefficient C D and thus C D = C T (thrust coefficient). The drag force is defined as:
F D = 1 2 ρ α C D A e U 0 2
where ρ is the density of water (equal to 1025 kg/m3), α a correction factor (set to 1 in this study), C D is the drag coefficient, A e is the effective area of turbine exposed to current and U 0 is the upstream current speed. A generic turbine design was used with a diameter of 20 m, with a constant C D setting of 0.85 and the hub height was set to 17 m [42] to represent a first generation of tidal stream energy converters. Each turbine was arranged in one grid cell of the numerical model for a better representation of the interactions between the turbines. The array layout was arranged with 2.5D (with D the turbine diameters) lateral spacing, recommended by the European Marine Energy Centre (EMEC) standard [43] and according to [44] inter-device spacing could be less than 3D. For the longitudinal spacing 6D was chosen as a compromise between 10D (recommended by EMEC) who allows better recovery of the flow and the characteristics (bathymetry and area of strong current) of the area.
To enable velocity profile comparisons between the 2D model and the 3D model, the power law defined as:
U z = z β h 1 / α U ¯
was used to transform the 2D results. This equation relates the velocity profile U z with bed roughness ( β = 0.4 to consider tidally energetic environment [45]), a power law coefficient α , h the water depth, U ¯ the depth average velocity and z being the height above the seabed. It should be pointed out that this equation does not take into account wind and wave-driven surface effects which could alter the velocity profile [46,47].

3. Validation

Calibration, validation techniques used for the models and results for the baseline are outlined in this section.

3.1. Calibration

3.1.1. Free Surface and Velocities

The 2D and 3D models were calibrated against five ADCP measurements (C1, CW2, CW3, CW4 and CWTb1 /cf. Table 2) from the first AUSTEn campaign in March 2018. Calibration was performed by adjusting the seabed roughness with the Manning number (M) in 2D and roughness height parameter in 3D, with seven tests to obtain the closest match of surface elevations, depth-averaged velocities (DAV) and velocity profiles between ADCP measurements and numerical model’s predictions, as shown in Figure 4. The phase shift between predicted and measured timeseries was removed following the EMEC guidelines [43]. Taylor diagrams were used to rank model skill depending on surface elevation and current speed accuracy against the field measurements at the location of the five ADCPs. These Taylor diagrams allow visual comparison of three statistical metrics (Standard deviation, Correlation coefficient and Centered Root Mean Square) and provide a summary of the performance of each simulation. For both surface elevations and DAV, the runs R4 (M = 32) and R9 (M = 30) in 2D showed good agreement with the five ADCP data (Figure 4). Three dimensional models increase computational cost so only R4 and R9 and R6 (reference for lower quality) were compared, with R4 the closest to the measurements (Figure 5a). Further comparisons between R4 and R9 were performed, examining current speed at the hub height, velocity profiles (Ebb and Flood) and harmonic analysis of the principal constituent. Ebb direction was defined for this work as the tidal current going from Bass Strait to Tasman Sea (eastwards) and flood from Tasman Sea to Bass Strait (westwards). R4 and R9 showed similar results for the velocities at the hub height layer (Figure 5b,c).
For velocity profiles, the data were decomposed in ebb and flood, with results showing good agreement with the five ADCP measurements (Figure 6). The models correlated well with all stations: both runs are pretty close to the measurements. The highest correlation between in situ data and numerical model was found for simulation run R4 at station CW2 deployed 2.2 km southwest of Clarke Island, CW4 deployed 7.5 km north of Swan Island and CWTb1 in the deep pocket at the eastwards entrance of the channel. Those results confirmed the ability of the model to represent accurately the vertical profile distribution of the velocity in Banks Strait.

3.1.2. Harmonic Analysis

Harmonic constituents were derived from the timeseries of the 35 days simulations, the dominant harmonic constituent M2 was determined using the Matlab program UTide [48] from surface elevation and current velocities for both the MIKE models and the ADCP stations. The M2 constituents found were validated for both amplitude and phase for runs R4 and R9, with similar differences (Table 3).

3.1.3. Validation

Among all comparisons, R4 had the highest correlation for both 2D and 3D depth average velocities and velocity profiles, therefore it was chosen to represent the conditions in Banks Strait. The run R4 was then run for one year in 2D and validated against measurements for three different periods with different ADCP: Autumn (22/03/2018–09/07/2018) with C1, CW2 and CWTb1; Winter (13/07/2018–06/09/2018) with CW1, CW2 bis and CW4 bis and Summer (06/12/2018–15/02/2019) with C1 bis. These large validation datasets give confidence in the ability of the model to reproduce tidal currents in Banks Strait. All periods showed good correlation for DAV, with slightly lower quality results for Winter period with a correlation coefficient of 0.86 (Figure 7). To evaluate the influence of waves on simulation accuracy in Winter, a small subset of the data was analysed: a period of smaller significant wave height in Winter (31/07/2018 to 06/08/2018) was considered from the measurements of waves for CW4 bis, leading to a better correlation coefficient: R = 0.9. The 2D models were also validated using data available for tidal range from the Australian Baseline Sea Level Monitoring Project [49] for the following stations: Portland (141.6136, −38.3439), Lorne (143.9889, −38.3439), Burnie (145.9149, −41.05) and Spring Bay (147.9308, −42.5464). Comparisons over 35 days found average correlation coefficient R = 0.95 (minimum at 0.9) indicating good agreement with the measured data.

3.2. The Hydrodynamics in Banks Strait: Baseline Results

The models were applied to develop a baseline case with no tidal turbine and hence no tidal energy extraction to evaluate the current circulation in the strait and determine suitable locations for the placement of tidal turbine farms. Similar to previous studies where the strongest currents were found in the middle of the channel south of Clarke Island [18], the 3D model results reveal mean current speeds of 1.2 m/s and maximum speeds of 2.5 m/s located at (596878/−4498311 UTM 55) for the period of 22/03/2018 to 26/04/2018 (dataset of 35 days). The average maximum for the mid water column is illustrated in Figure 8 for two spring-neap cycles, with the time series of current speed at the surface for one spring-neap cycle at the location of the maximum of velocity. The modelled currents displayed 2 eddies near the area of strong currents: with one at the vicinity of the sand dunes, the presence of this feature could be strongly linked to this eddy (Figure 9).

4. Influence of the Tidal Farms on the Hydrodynamics

To evaluate the influence of tidal farms on the hydrodynamics (current speed, bed shear stress), simulations were performed comparing the spatial influence of two tidal farms between 2D and 3D, the influence near the seabed and the hydrodynamic influence on sand dunes near Clarke Island to evaluate potential disturbance in the ecologic system. The tidal farms are located south of Clarke Island where the strongest currents were found and where the bathymetry did not exceed 50 m (operational water depth).

4.1. Comparison of Model Spatial Extent 2D/3D

Simulations were performed in 2D and 3D in Banks Strait with and without tidal farms of 100 turbines (TF100) and 300 turbines (TF300). The spatial extent of the difference in mean current speed over 35 days between the baseline (no tidal energy extraction) and the scenario TF100 was limited to less than 5km from the tidal farm location for both 2D and 3D DAV (Figure 10a,c). For the TF300 scenario, changes in mean current speed were found to exceed to more than 7 km (Figure 10b,d). The change in mean current speed was greatest at the location of the farms due to the removal of energy by tidal turbines. For the TF300 case in 2D and 3D, a slight increase of the current speed (0.04 m/s) is noticed between the farm and Clarke Island where sand dunes are present, and also south of the tidal farm, resulting from the barrier effect of the turbines on the flow. In both models (2D/3D), the scenario TF100 (Figure 10a,c) had less influence on the reduction of current speed than the scenario TF300: a maximum of 0.05 m/s against 0.12 m/s for TF300 (Figure 10b,d). The values are quite similar between 2D and 3D, the major difference is the spatial extent where 2D displays a larger influence which could be considered as a conservative measure. For more insight about the difference between the two arrays, the mean current speed over 35 days was extracted at the mid water column in the 3D scenario (Figure 11). The influence of tidal energy extraction is higher at the location of the turbines for the TF300 scenario: above 0.2 m/s reduction against around 0.1 m/s for TF100. These results are slightly lower than [5] who found a mean reduction of 0.3 m/s with a farm of 290 turbines in the Alderney Race.
A profile AA’ (location on Figure 8) parallel to the direction of the currents was extracted to examine the spatial influence of the tidal farms on depth averaged velocities for the 35 days dataset in 2D and 3D (Figure 12). The drop for the baseline case at 2 km is caused by the presence of a canyon. The 2D model baseline results for current speed were overall 5% higher than the 3D DAV results, corresponding to an increase of approximately of 0.05 m/s along the profile. A significant drop in current speed was noticed for both farms at their location, with a greater reduction for the scenario TF300 as expected due to higher levels of tidal energy extraction, equivalent to a reduction of 7% in the 3D model (Figure 12). For both dimensions, the differences were lower at approximately 13km from the tidal farm reducing to less than 1% up to 0% at 18 km towards the Tasman Sea for the two tidal farms (Table 4). The influence of the 100 turbine array was minimal in the near field with less than 1% of difference in mean current speed. For smaller scale tidal farm development, this is a positive outcome.
For 2D models, the difference in the magnitude of velocity between the baseline and scenarios with tidal farms was greater than in 3D DAV at the vicinity of the arrays, however in the far field the average difference was the same. The 2D model results match the 3D DAV results for the TF100 scenario with a maximum difference of 0.4% and are similar in the far field for the TF300 scenario (Table 4). The 2D models overestimated slightly the current speed compared to the 3D DAV results, with the difference between baseline and scenario with turbines highest at the location of the TF300. This should be considered if a preliminary resource assessment is performed in 2D. At the location of the tidal turbines, the velocities difference between 2D/3D reached an average of 0.1 m/s corresponding to 12%, which is lower than the results of [13] who found a 20% difference in an idealised channel similar to Pentland Firth. However, the difference between 2D and 3D DAV was found to be negligible in the far field of the tidal farms.
Small discrepancies were observed between 2D and 3D models in the near and far field. Given the high computational cost required by 3D models, the authors recommend that preliminary assessments be performed in 2D. These results are in accordance to previous studies in an idealised channel for the difference of dimension with a focus on the bed level change. The 2D models are a good compromise between model accuracy and time management, especially for long temporal scale assessments for a first approximation.

4.2. Influence Near the Seabed

To enable comparisons of the 2D model to the 3D model at turbine hub heights, the power law defined in Equation (1) Section 2.3 was used. Previous studies have demonstrated that the 1/7th power law ( α = 7) is generally the most suitable to capture vertical processes in coastal oceanographic research; moreover, [45] showed that α = 7 correlated well with tidal energy sites conditions in the Irish Sea. However, guidelines [43,50,51] and field data observation [52] recommend a 1/10th power law in tidal energetic sites. To examine the influence of power law factors on simulation results, two α coefficients (7 and 10) were analysed with the velocity profile for four points (inside the tidal farm then at 0.9, 4.5 and 18 km from the tidal farm towards Bass Strait) for the TF300 model as shown in Figure 13. Results indicated that the 1/10th power law correlated the best with 3D results, as the difference was 14% instead of 25% for the α = 7 model. There was a significant difference of velocity at the location of the tidal farm between 2D and 3D. The 3D model represented well the decrease of the flow in the water column due to the presence of the turbines. This velocity deficit at the hub height cannot be reproduced in the 2D model with the power law and thus lead to a misrepresentation of the current speed in 2D at the location of the tidal farm confirming the findings of [13]. It should also be noted that some tidal energy sites have more complex velocity profiles which do not fit the power law equation [53]. Future research using coupled hydrostatic and non-hydrostatic models will be able to capture flow properties at scale ranging from ocean to tidal turbines blades enabling the accurate modelling of all hydrodynamics flow properties including wake environment and velocity profiles.
Friction between the near seabed currents and the seafloor creates a layer of shear and increased turbulence above the seabed called the bottom boundary layer. This is the layer where the flow is the most influenced by the seafloor and alterations in this layer could lead to significant impact in sediment transport patterns. The profile AA’ was extracted at the bottom layer for the mean current speed (Figure 14a,b) with a coefficient α equal to 10 for the power law. An increase of current speed was noticed at the location of the tidal farm only for the 3D models: with magnitudes of 1 cm/s for TF100 and 3 cm/s for TF300 (2% and 5%, respectively) found. Away from the tidal farm, the reduction in current speeds are quite similar for 2D and 3D: for TF300 6% and 5.6% (3D, 2D, respectively) and 2% for TF100 at 2.2 km from the tidal farm. The differences fade away in the far field reaching less than 1% at 13.4 km for both dimensions. The 3D models appear to model the acceleration caused by the constriction of the flow beneath the turbines, which the 2D models cannot capture.
In the 2D model, the bed shear stress is derived from the DAV velocities, whereas the bed shear stress is determined from the bottom cell velocity in the 3D model; therefore, the profile for the bed shear stress followed the same fluctuation as in the bottom layer in 3D with an increase at the location of tidal farm (Figure 14c,d and Figure 15). Around the tidal farms the bed shear stress was reduced by 7% (0.22 Pa) for TF100 arrays in 3D, 6% for 2D, with the differences not significant in the far field. The reduction of bed shear stress for the TF300 scenario was greater reaching 20% (equivalent to 0.58 Pa) around the tidal farm in 3D and 15% (0.51 Pa) in 2D. The 2D results were similar to the 3D results for the scenario TF100 and TF300 in the far field. These results vary compared to the influence of turbines in Ramsey Sound as shown by [7] who found an average reduction of 2.3 Pa in the Ramsey sound with less devices. The peak reduction found here is 0.9 Pa, 25 times lower than [4] with 25 Pa in the Pentland Firth. These discrepancies can be explained by the faster flows occurring in these two sites and therefore a greater alteration in hydrodynamics. Influence on bed shear stress was observed with the TF100 scenario and was greater with the 300 turbines arrays with a spatial influence of 7 km instead of 4.5 km for TF100 (Figure 15), results close to the spatial extent of 12 km found by [7]. Reduction in the bed shear stress around the tidal farm suggests that accretion of sediments may appear and change the dynamics in the channel as it can be seen on maps on Figure 15. The increase of bed shear stress at the location of the farm and at the top and bottom suggest that additional scour and more erosion could occur at the location of the sand dunes between the tidal farm and Clarke Island.
A closer look was given to the region of the sand dunes near Clarke Island given the proximity to a potential tidal farm. Alteration in current speed and bed shear stress could lead to a disturbance of the benthic ecology system. Bed shear stress results in 3D and 2D were compared between the case with no tidal extraction and the two scenarios TF100 and TF300 (Figure 16). The change in the magnitude for the 100 turbines arrays was less than 0.12 Pa and less than 0.26 Pa for the 300 turbines case for both 2D/3D. The difference between 2D and 3D results were not significant: a maximum of 4.5% difference with an average difference around 0.06 Pa for all the scenarios, but the models did not display the same fluctuations. The slight change in bed shear stress at this location needs to be better understood, with more studies on the dynamics of these sand dunes. To pursue the EIA, the focus needs to be on the dynamics of sediment, especially at the location of the sand dunes which are of economic importance. They may provide sediment for the close beaches of Clarke Island, protect them from erosion during storms and may contribute to the fishing industry being a source of nutrients for many fish species [54,55].
These results reinforces previous findings [7,8] that, for preliminary assessment of oceanographic changes of a potential tidal energy sites, 2D models are suitable for examining far field parameters (current speed and bed shear stress). These 2D studies can give a first insight of the potential influence of arrays on a longer time frame given the reduced computational costs when compared to 3D models. Thus, 2D models are assets for the dynamics of sediment transport and are an effective tool to understand the influence of the removal of energy. However further work is needed to understand the capability of 2D models with more complex and challenging scenarios such as sand dunes migration. 3D models remain essential to have more insight at the location of tidal farm; however, they become less powerful in the near and far field as it gives similar results as the 2D model in terms of hydrodynamic environmental influence.

5. Conclusions

The hydrodynamics of Banks Strait was examined using two high resolution models in 2D and 3D, with resolution varying from 3 km to 20 m in the area of interest to simulate the presence of tidal farms. The developed models were successfully validated against a large dataset of ADCP field measurements, with results indicating good agreement for surface elevation, depth average velocities and velocity profiles. Validation for a full year was then performed for the 2D model with close agreement to the field data found. Those results gave confidence in the ability of these models to reproduce the hydrodynamics in Banks Strait.
The implementation of tidal turbines in the numerical models showed the footprint effect of tidal turbines on the hydrodynamics (current speed, bed shear stress). Scenario with 100 turbines showed minimal influence on the velocities in the near field, with changes in velocity only found less than 3km from the farm. For the alteration of bed shear stress, this scenario did not display any significant difference compared to the baseline, which is a relatively good outcome for the deployment of a small tidal farm. The spatial influence of the 300 turbines farm was greater reaching 7 km from the farm for the difference in the current speed. A small number of devices will not change the system in a significant way compared to the natural variability; but, as these studies have shown, large scale arrays have the potential to influence the natural process. These alterations can lead to greater changes on the dynamics of sediments.
In the far field, the difference between the two dimensions 2D/3D was negligible, confirming that 2D models are suitable for investigating the influence of tidal arrays for a preliminary environmental assessment. At the location of the tidal farm 2D models were found to have slightly higher values for the influence of tidal energy extraction and thus precaution should be taken for resource assessment. Regarding local features here in Banks Strait, more work is needed to fully understand the effect of the tidal arrays on the sand dunes compared to the natural variability.
The findings of this paper are a first step in the understanding of the influence of tidal energy extraction in Banks Strait and gave a solid basis to reduce the uncertainty in the hydrodynamics for future investigation of sediment transport in the region. Future work will focus on the creation and calibration of a sediment transport model with in situ data (multi-beam, echosounder, bottom grabs, sediment traps and sub-bottom data) collected by the AUSTEn project. The migration of sand dunes will also be of a high interest to better understand the influence of anthropogenic activities on the physical environment.

Author Contributions

C.A., conceptualization, methodology, investigation, software, formal analysis, validation, writing—original draft, visualization; P.M., supervision, software, conceptualization, writing—review and editing; J.-R.N., supervision, project administration, funding acquisition, writing—review and editing; R.C., supervision, funding acquisition, writing—review and editing; I.P., writing—review and editing, funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

The AUSTEn project is co-funded by the Australian Renewable Energy Agency through the Advancing Renewables Programme (grant G00902), the Australian Maritime College (University of Tasmania), the University of Queensland, CSIRO, our industry partners MAKO Tidal Turbines, Sabella SAS and SIMEC Atlantis Energy.

Acknowledgments

The authors wish to acknowledge DHI for their role providing MIKE21/3 software throughout this PhD research, and for their support especially Méven Huiban. The authors are also grateful to Camille Couzi for the process of the ADCP. The authors also would like to acknowledge the support of IMAREST through the award “Laurie Prandolini Research Fellowship” given to this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
ADCPAcoustic Doppler Current Profiler
AEPAnnual Energy Production
AHSAustralian Hydrography Service
ARENAAustralian Renewable Energy Agency
AUSTEnAustralian Tidal Energy
AWACAcoustic Waves And Current Profiler
CSIROCommonwealth Scientific and Industrial Research Organisation
CTDConductivity Temperature Depth
DAVDepth Average Velocities
EIAEnvironmental Impact Assessment
EMECEuropean Marine Energy Centre
IECInternational Electrotechnical Commission
IMOSIntegrated Marine Observing System
TFTidal Farm

References

  1. OES. Ocean Energy System Annual Report: An Overview of Ocean Energy Activities in 2019; OES: Paris, France, 2019. [Google Scholar]
  2. Neill, S.P.; Litt, E.J.; Couch, S.J.; Davies, A.G. The impact of tidal stream turbines on large-scale sediment dynamics. Renew. Energy 2009, 34, 2803–2812. [Google Scholar] [CrossRef]
  3. Frid, C.; Andonegi, E.; Depestele, J.; Judd, A.; Rihan, D.; Rogers, S.I.; Kenchington, E. The environmental interactions of tidal and wave energy generation devices. Environ. Impact Assess. Rev. 2012, 32, 133–139. [Google Scholar] [CrossRef]
  4. Martin-Short, R.; Hill, J.; Kramer, S.; Avdis, A.; Allison, P.; Piggott, M. Tidal resource extraction in the Pentland Firth, UK: Potential impacts on flow regime and sediment transport in the Inner Sound of Stroma. Renew. Energy 2015, 596–607. [Google Scholar] [CrossRef] [Green Version]
  5. Thiébot, J.; Bailly du Bois, P.; Guillou, S. Numerical modeling of the effect of tidal stream turbines on the hydrodynamics and the sediment transport—Application to the Alderney Race (Raz Blanchard), France. Renew. Energy 2015, 75, 356–365. [Google Scholar] [CrossRef]
  6. Fairley, I.; Masters, I.; Karunarathna, H. The cumulative impact of tidal stream turbine arrays on sediment transport in the Pentland Firth. Renew. Energy 2015, 80, 755–769. [Google Scholar] [CrossRef] [Green Version]
  7. Haverson, D.; Bacon, J.; Smith, H.C.; Venugopal, V.; Xiao, Q. Modelling the hydrodynamic and morphological impacts of a tidal stream development in Ramsey Sound. Renew. Energy 2018, 126, 876–887. [Google Scholar] [CrossRef]
  8. Robins, P.E.; Neill, S.P.; Lewis, M.J. Impact of tidal-stream arrays in relation to the natural variability of sedimentary processes. Renew. Energy 2014, 72, 311–321. [Google Scholar] [CrossRef] [Green Version]
  9. Neill, S.P.; Hashemi, M.R.; Lewis, M.J. Tidal energy leasing and tidal phasing. Renew. Energy 2016, 85, 580–587. [Google Scholar] [CrossRef] [Green Version]
  10. Ahmadian, R.; Falconer, R.A. Assessment of array shape of tidal stream turbines on hydro-environmental impacts and power output. Renew. Energy 2012, 44, 318–327. [Google Scholar] [CrossRef]
  11. Neill, S.P.; Jordan, J.R.; Couch, S.J. Impact of tidal energy converter (TEC) arrays on the dynamics of headland sand banks. Renew. Energy 2012, 37, 387–397. [Google Scholar] [CrossRef]
  12. Gallego, A.; Side, J.; Baston, S.; Waldman, S.; Bell, M.; James, M.; Davies, I.; O’Hara Murray, R.; Heath, M.; Sabatino, A.; et al. Large scale three-dimensional modelling for wave and tidal energy resource and environmental impact: Methodologies for quantifying acceptable thresholds for sustainable exploitation. Ocean Coast. Manag. 2017, 147, 67–77. [Google Scholar] [CrossRef] [Green Version]
  13. Goward Brown, A.J.; Neill, S.P.; Lewis, M.J. Tidal energy extraction in three-dimensional ocean models. Renew. Energy 2017, 114, 244–257. [Google Scholar] [CrossRef] [Green Version]
  14. Manasseh, R.; McInnes, K.L.; Hemer, M.A. Pioneering developments of marine renewable energy in Australia. Int. J. Ocean Clim. Syst. 2017, 8, 50–67. [Google Scholar] [CrossRef] [Green Version]
  15. AUSTeN. Australian Tidal Energy. 2017. Available online: http://austen.org.au/ (accessed on 30 November 2019).
  16. Penesis, I.; Hemer, M.; Cossu, R.; Hayward, J.; Nader, J.R.; Rosebrock, U.; Grinham, A.; Sayeef, S.; Osman, P.; Marsh, P.; et al. Tidal energy in Australia—Assessing resource and feasibility to Australia’s future energy mix. In Proceedings of the 4th Asian Wave and Tidal Energy Conference, Taipei, Taiwan, 9–13 September 2018; p. 507. [Google Scholar]
  17. Behrens, S.; Griffin, D.; Hayward, J.; Hemer, M.; Knight, C.; McGarry, S.; Osman, P.; Wright, J. Ocean Renewable Energy: 2015–2050: An analysis of ocean energy in Australia. North Ryde CSIRO 2012. [Google Scholar] [CrossRef]
  18. Marsh, P.; Penesis, I.; Nader, J.R.; Couzi, C.; Cossu, R. Assessment of tidal current resources in Banks. In Proceedings of the 13th European Wave and Tidal Energy Conference, Napoli, Italy, 1–6 September 2019; pp. 1–10. [Google Scholar]
  19. Wijeratne, E.M.; Pattiaratchi, C.B.; Eliot, M.; Haigh, I.D. Tidal characteristics in Bass Strait, south-east Australia. Estuar. Coast. Shelf Sci. 2012, 114, 156–165. [Google Scholar] [CrossRef]
  20. Sandery, P.A.; Kämpf, J. Winter-Spring flushing of Bass Strait, South-Eastern Australia: A numerical modelling study. Estuar. Coast. Shelf Sci. 2005, 63, 23–31. [Google Scholar] [CrossRef]
  21. Sandery, P.A.; Kämpf, J. Transport timescales for identifying seasonal variation in Bass Strait, south-eastern Australia. Estuar. Coast. Shelf Sci. 2007, 74, 684–696. [Google Scholar] [CrossRef]
  22. McIntosh, P.C.; Bennett, A.F. Open Ocean Modeling as an Inverse Problem: M 2 Tides in Bass Strait. J. Phys. Oceanogr. 1984, 14, 601–614. [Google Scholar] [CrossRef] [Green Version]
  23. McInnes, K.L.; O’Grady, J.G.; Hemer, M.; Macadam, I.; Abbs, D.J.; White, C.J.; Corney, S.P.; Grose, M.R.; Holz, G.K.; Gaynor, S.M.; et al. Climate Futures for Tasmania: Extreme Tide and Sea-Level Events Technical Report; Technical Report January; Antarctic Climate and Ecosystems Cooperative Research Centre: Hobart, Tasmania, 2011. [Google Scholar]
  24. Fandry, C.B. Development of a Numerical Model of Tidal and Wind-driven Circulation in Bass Strait. Mar. Freshw. Res. 1981, 32, 9–29. [Google Scholar] [CrossRef]
  25. Fandry, C.B.; Hubbert, G.D.; McIntosh, P.C. Comparison of predictions of a numerical model and observations of tides in bass strait. Mar. Freshw. Res. 1985, 36, 737–752. [Google Scholar] [CrossRef]
  26. Wijffels, S.E.; Beggs, H.; Griffin, C.; Middleton, J.F.; Cahill, M.; King, E.; Jones, E.; Feng, M.; Benthuysen, J.A.; Steinberg, C.R.; et al. A fine spatial-scale sea surface temperature atlas of the Australian regional seas (SSTAARS): Seasonal variability and trends around Australasia and New Zealand revisited. J. Mar. Syst. 2018, 187, 156–196. [Google Scholar] [CrossRef]
  27. Fortunato, A.B.; Bertin, X.; Oliveira, A. Space and time variability of uncertainty in morphodynamic simulations. Coast. Eng. 2009, 56, 886–894. [Google Scholar] [CrossRef]
  28. Pinto, L.; Fortunato, A.B.; Freire, P. Sensitivity analysis of non-cohesive sediment transport formulae. Cont. Shelf Res. 2006, 26, 1826–1839. [Google Scholar] [CrossRef]
  29. Nash, S.; Phoenix, A. A review of the current understanding of the hydro-environmental impacts of energy removal by tidal turbines. Renew. Sustain. Energy Rev. 2017, 80, 648–662. [Google Scholar] [CrossRef]
  30. Galibert, G. IMOS Toolbox Version 2.x.x. 2020. Available online: https://github.com/aodn/imos-toolbox (accessed on 21 March 2019).
  31. Scherelis, C.; Penesis, I.; Marsh, P.; Cossu, R.; Hemer, M.; Wright, J. Relating fish distributions to physical characteristics of a tidal energy candidate site in the Banks Strait Australia. In Proceedings of the 13th European Wave and Tidal Energy Conference, Napoli, Italy, 1–6 September 2019; pp. 1–8. [Google Scholar]
  32. Perez, L.; Cossu, R.; Couzi, C.; Penesis, I. Wave-Turbulence Decomposition Methods Applied to Tidal Energy Site Assessment. Energies 2020, 13, 1245. [Google Scholar] [CrossRef] [Green Version]
  33. Coles, D.S.; Blunden, L.S.; Bahaj, A.S. Assessment of the energy extraction potential at tidal sites around the Channel Islands. Energy 2017, 124, 171–186. [Google Scholar] [CrossRef]
  34. DHI. MIKE 21 & MIKE 3 Flow Model FM, Hydrodynamic Module Scientific Documentation; DHI: Hesholm, Denmark, 2017. [Google Scholar]
  35. Whiteway, T. Australian Bathymetry and Topography Grid. Geosci. Aust. Canberra 2009. [Google Scholar] [CrossRef]
  36. Australian Hydrographic Office Charts. Available online: http://www.hydro.gov.au/prodserv/paper/auspapercharts.htm (accessed on 1 March 2018).
  37. Yongcun, C.; Ole Baltazar, A. Improvement in global ocean tide model in shallow water regions, 2010. In Proceedings of the OSTST, Lisbon, Portugal, 18–22 October 2010. [Google Scholar]
  38. International Electrotechnical Commission. IEC/TS 62600-201: Marine Energy—Wave, Tidal and Other Water Current Converters and Characterization; International Electrotechnical Committee: Genova, Switzerland, 2015. [Google Scholar]
  39. Copernicus Climate Change Service (C3S). ERA5: Fifth generation of ECMWF Atmospheric Reanalyses of the Global Climate. Copernicus Climate Change Service Climate Data Store (CDS). 2017. Available online: https://cds.climate.copernicus.eu/cdsapp#!/home (accessed on 5 January 2019).
  40. Swan Island: Tasmania Daily Weather Observations. Available online: http://www.bom.gov.au/climate/dwo/IDCJDW7051.latest.shtml (accessed on 5 January 2019).
  41. Soulsby, R. Dynamics of Marine Sands; Thomas Telford Publishing: London, UK, 1997; Available online: https://www.icevirtuallibrary.com/doi/pdf/10.1680/doms.25844 (accessed on 5 January 2018). [CrossRef]
  42. Waldman, S.; Bastón, S.; Nemalidinne, R.; Chatzirodou, A.; Venugopal, V.; Side, J. Implementation of tidal turbines in MIKE 3 and Delft3D models of Pentland Firth & Orkney Waters. Ocean Coast. Manag. 2017, 147, 21–36. [Google Scholar] [CrossRef] [Green Version]
  43. EMEC Orkney. Environmental Impact Assessment (EIA). Guidance for Developers at the European Marine Energy Centre; Technical Report; EMEC Orkney: Stromness, UK, 2008. [Google Scholar]
  44. Thiébot, J.; Guillou, N.; Guillou, S.; Good, A.; Lewis, M. Wake field study of tidal turbines under realistic flow conditions. Renew. Energy 2020, 151, 1196–1208. [Google Scholar] [CrossRef]
  45. Lewis, M.; Neill, S.; Robins, P.; Hashemi, M.; Ward, S. Characteristics of the velocity profile at tidal-stream energy sites. Renew. Energy 2017, 114, 258–272. [Google Scholar] [CrossRef]
  46. Gunn, K.; Stock-Williams, C. On validating numerical hydrodynamic models of complex tidal flow. Int. J. Mar. Energy 2013, 3–4, e82–e97. [Google Scholar] [CrossRef]
  47. Sellar, B.; Wakelam, G. Characterisation of Tidal Flows at the European Marine Energy Centre in the Absence of Ocean Waves. Energies 2018, 11, 176. [Google Scholar] [CrossRef] [Green Version]
  48. Codiga, D. UTide Unified Tidal Analysis and Prediction Functions. MATLAB Central File Exchange. 2020. Available online: https://www.mathworks.com/matlabcentral/fileexchange/46523-utide-unified-tidal-analysis-and-prediction-functionsl (accessed on 22 January 2020).
  49. Australian Baseline Sea Level Monitoring Project Hourly Sea Level and Meteorological Data. Available online: http://www.bom.gov.au/oceanography/projects/abslmp/data/index.shtml (accessed on 1 March 2018).
  50. Hagerman, G.; Polagye, B.; Bedard, R.; Previsic, M. Methodology for Estimating Tidal Current Energy Resources and Power Production by Tidal In-Stream Energy Conversion (TISEC) Devices; EPRI North American Tidal In Stream Power Feasibility Demonstration Project, 2006. Available online: https://tethys.pnnl.gov/sites/default/files/publications/Tidal_Current_Energy_Resources_with_TISEC.pdf (accessed on 22 January 2020).
  51. Gooch, S.; Thomson, J.; Polagye, B.; Meggitt, D. Site characterization for tidal power. In Proceedings of the OCEANS 2009, Venice, Italy, 21–25 September 2009; pp. 1–10. [Google Scholar]
  52. O’Rourke, F.; Boyle, F.; Reynolds, A. Ireland’s tidal energy resource; An assessment of a site in the Bulls Mouth and the Shannon Estuary using measured data. Energy Convers. Manag. 2014, 87, 726–734. [Google Scholar] [CrossRef] [Green Version]
  53. Gunn, K.; Stock-Williams, C. Fall of Warness 3D Model Validation Report: ETI REDAPT MA1001 PM14 MD5.2; Technical Report; E.ON New Build Technology: Coventry, UK, 2012. [Google Scholar]
  54. Dyer, K.R.; Huntley, D.A. The origin, classification and modelling of sand banks and ridges. Cont. Shelf Res. 1999, 19, 1285–1330. [Google Scholar] [CrossRef]
  55. Gibbs, C.; Tomczak, M.; Longmore, A. The Nutrient Regime of Bass Strait. Mar. Freshw. Res. 1986, 37, 451–466. [Google Scholar] [CrossRef]
Figure 1. (top) Map of Australia with the location of Tasmania, and (bottom) map of Tasmania showing the location of the Banks Strait with (right) a zoom on the strait.
Figure 1. (top) Map of Australia with the location of Tasmania, and (bottom) map of Tasmania showing the location of the Banks Strait with (right) a zoom on the strait.
Energies 13 05326 g001
Figure 2. Location of Acoustic Doppler Current Profiler (ADCP) stations and high-resolution bathymetry obtained from the multi-beam survey in Banks Strait with a picture of an ADCP mounted on a frame.
Figure 2. Location of Acoustic Doppler Current Profiler (ADCP) stations and high-resolution bathymetry obtained from the multi-beam survey in Banks Strait with a picture of an ADCP mounted on a frame.
Energies 13 05326 g002
Figure 3. Domain of the numerical domain and a zoom in Banks Strait (Coordinates in UTM 55).
Figure 3. Domain of the numerical domain and a zoom in Banks Strait (Coordinates in UTM 55).
Energies 13 05326 g003
Figure 4. Taylor Diagrams for depth averaged velocities (left) and surface elevations (right) in 2D for 35 days for 5 stations (MEAS: measured values from ADCP; RMSD: Root-mean square deviation).
Figure 4. Taylor Diagrams for depth averaged velocities (left) and surface elevations (right) in 2D for 35 days for 5 stations (MEAS: measured values from ADCP; RMSD: Root-mean square deviation).
Energies 13 05326 g004
Figure 5. Taylor diagram for (a) depth averaged velocities in 3D for 35 days, (b) magnitude of velocities at the hub height in 3D for 35 days and (c) timeseries of velocities at the hub height in 3D for Acoustic Doppler Current Profiler C1 for 1 day.
Figure 5. Taylor diagram for (a) depth averaged velocities in 3D for 35 days, (b) magnitude of velocities at the hub height in 3D for 35 days and (c) timeseries of velocities at the hub height in 3D for Acoustic Doppler Current Profiler C1 for 1 day.
Energies 13 05326 g005
Figure 6. Calibration analysis for 3D current profiles for Ebb and Flood: Average for 35 days at the location of ADCP stations.
Figure 6. Calibration analysis for 3D current profiles for Ebb and Flood: Average for 35 days at the location of ADCP stations.
Energies 13 05326 g006
Figure 7. Taylor diagram for depth average velocities between ADCP stations (localisation on Figure 2) and 2D MIKE results (March 2018–February 2019).
Figure 7. Taylor diagram for depth average velocities between ADCP stations (localisation on Figure 2) and 2D MIKE results (March 2018–February 2019).
Energies 13 05326 g007
Figure 8. At the top: Statistical maximum for current speed over two spring-neap cycles at the mid water column with the location of the profile AA’ (distance from location of tidal farm in km) for the baseline case. At the bottom: Current speed at the location of a possible tidal farm (596878, −4498311 /UTM 55) over two spring-neap cycles for the 3D baseline case at the surface (no tidal energy extraction).
Figure 8. At the top: Statistical maximum for current speed over two spring-neap cycles at the mid water column with the location of the profile AA’ (distance from location of tidal farm in km) for the baseline case. At the bottom: Current speed at the location of a possible tidal farm (596878, −4498311 /UTM 55) over two spring-neap cycles for the 3D baseline case at the surface (no tidal energy extraction).
Energies 13 05326 g008
Figure 9. Spatial distribution of current direction near the area of strong currents around peak Flood. A and B highlight the presence of eddies near the area of interest.
Figure 9. Spatial distribution of current direction near the area of strong currents around peak Flood. A and B highlight the presence of eddies near the area of interest.
Energies 13 05326 g009
Figure 10. Difference in depth average velocities between Baseline and scenarios with tidal farm in 2D (a,b) and 3D (c,d) over 35 days. Negative reduction compared to the baseline is shown in blue (corresponding to increase of velocity) and positive reduction in red (corresponding to a deficit in velocity). (Legend valid for all maps/dashed lines with blue squares represent 5 km limit from the tidal farm).
Figure 10. Difference in depth average velocities between Baseline and scenarios with tidal farm in 2D (a,b) and 3D (c,d) over 35 days. Negative reduction compared to the baseline is shown in blue (corresponding to increase of velocity) and positive reduction in red (corresponding to a deficit in velocity). (Legend valid for all maps/dashed lines with blue squares represent 5 km limit from the tidal farm).
Energies 13 05326 g010
Figure 11. Difference in mean current speed over 35 days between the baseline and scenarios with tidal farm at the mid water column in 3D showing the effect of tidal energy extraction. Negative reduction compared to the baseline is shown in blue (corresponding to increase of velocity) and positive reduction in red (corresponding to a deficit in velocity).
Figure 11. Difference in mean current speed over 35 days between the baseline and scenarios with tidal farm at the mid water column in 3D showing the effect of tidal energy extraction. Negative reduction compared to the baseline is shown in blue (corresponding to increase of velocity) and positive reduction in red (corresponding to a deficit in velocity).
Energies 13 05326 g011
Figure 12. Comparison of depth averaged velocities at the profile AA’ (location Figure 8) for an average of 35 days.
Figure 12. Comparison of depth averaged velocities at the profile AA’ (location Figure 8) for an average of 35 days.
Energies 13 05326 g012
Figure 13. Velocity profile comparison due to uncertainty of alpha. The power law for four points of the profile AA’ (location on Figure 8) showing velocity deficit at the location of the tidal farm.
Figure 13. Velocity profile comparison due to uncertainty of alpha. The power law for four points of the profile AA’ (location on Figure 8) showing velocity deficit at the location of the tidal farm.
Energies 13 05326 g013
Figure 14. Change in magnitude velocity at the bottom layer (a,b) and in the bed shear stress (c,d) due to energy extraction over 35 days between baseline and TF scenarios in 2D and 3D at the profile AA.
Figure 14. Change in magnitude velocity at the bottom layer (a,b) and in the bed shear stress (c,d) due to energy extraction over 35 days between baseline and TF scenarios in 2D and 3D at the profile AA.
Energies 13 05326 g014
Figure 15. Change in the magnitude of bed shear stress (Pa) between baseline and TF scenarios in 3D (a,b) and 2D (c,d).
Figure 15. Change in the magnitude of bed shear stress (Pa) between baseline and TF scenarios in 3D (a,b) and 2D (c,d).
Energies 13 05326 g015
Figure 16. Location of the profile BB’ with bathymetry (m) of the area, Right: mean bed shear stress over 35 days for the profile BB’ in the sand dunes area for all 3D/2D scenarios.
Figure 16. Location of the profile BB’ with bathymetry (m) of the area, Right: mean bed shear stress over 35 days for the profile BB’ in the sand dunes area for all 3D/2D scenarios.
Energies 13 05326 g016
Table 1. Tidal energy trials in Australia adapted from [1,14].
Table 1. Tidal energy trials in Australia adapted from [1,14].
LocationYearDevice TypeCompanyPower RatingIndicative Depth
of Operation (m)
Darwin1996Tyson Turbine fluted cone
axial turbine (pontoon)
Northern Territory University
(NTU)
<1 kW1
Darwin1998NTU Swenson axial turbine
(pontoon)
NTU2.2 kW1–2
King Sound2000–
2014
Barrage dam
investigation
Tidal Energy Australia
and Hydro Tasmania
40 MW
Clarence river2004Aquanator reaction plates
in closed loop track
Atlantis Energy Ltd5 kW2
San Remo2006AquanatorAtlantis Resources100 kW2
San Remo2007Submersible tidal
generator
HydroGen Power
Industries
5–50 kWup to 20
Brisbane2007Submersible tidal
generator
HydroGen Power
Industries
5–50 kWup to 20
San Remo/
Stony Point
2008Floating tidal
turbine laboratory
EnGen Institute2
Corio Bay2008Solon ducted
axial turbine
Atlantis Resources160 kW10
Melbourne20093D printed multi-axis
turbine array
Cetus1 kW1
Newcastle2012Sea Urchin
axial turbine (pontoon)
Elementary Energy Tech.2 kW1–2
San Remo2012Cross-flow turbine
tidal desalination
Infra Tidal15 kW2.5
Tamar Estuary2016Ducted axial turbine
(pontoon)
Mako10 kW2
Gladstone port2018Ducted axial turbineMako
Italic font indicates information is uncertain.
Table 2. Details of ADCP.
Table 2. Details of ADCP.
Name StationType of InstrumentLongitudeLatitudeDepth (m)Date of DeploymentEnd of Data CollectedCell Size (m)Length (Days)
CW2RDI Sentinel V50 500 k Hz148.10188−40.584846.4722/03/201811/07/20180.5110
C1RDI Workshorse 300 k Hz148.23882−40.672757.9417/03/201810/07/20182.0115
CW3Nortek AWAC 1 M Hz148.07778−40.545434.9522/03/201816/06/20180.586
CW4Nortek AWAC 1 M Hz148.09241−40.666430.6715/03/201809/06/20181.085
CWTb1Nortek Signature 500 k Hz148.22626−40.667263.5722/03/201809/07/20181.0109
CW1RDI Sentinel V50 500 Hz148.05684−40.529427.1112/07/201806/09/20180.557
CW2 bisRDI Sentinel V50 500 Hz148.20132−40.70146.0812/07/201822/09/20180.572
CW4bisNortek AWAC 1 M Hz148.34497−40.729625.4213/07/201808/09/20181.054
C1 bisRDI Workshorse 300 k Hz148.12498−40.689129.0705/12/201815/02/20191.072
Table 3. Comparison between modelled and observed M2 constituents (u,v).
Table 3. Comparison between modelled and observed M2 constituents (u,v).
Phase for M2 (degree)Amplitude for M2 (m)
3D2D 3D2D
StationsOBSR4R9R4R9OBSR4R9R4R9
CW21401391391391391.251.31.251.381.31
CW31181251251251250.990.970.931.030.97
CW41611651651651651.361.31.231.381.29
C11651581571571570.850.880.850.930.89
CWTb11581541531541530.880.910.880.960.92
Table 4. Change in magnitude of current speed along the profile for tidal farm TF100 & TF300 over 35 days.
Table 4. Change in magnitude of current speed along the profile for tidal farm TF100 & TF300 over 35 days.
Difference in DAV
TF300TF100
Distance From Tidal Farm (km)3D2D3D2D
07.39%8.44%3.67%3.69%
4.54.06%5.21%1.55%1.94%
92.05%2.64%0.73%0.92%
13.40.67%0.92%0.19%0.31%
18 (towards Bass Strait)0.62%0.84%0.18%0.26%
18 (towards Tasman Sea)0.21%0.35%0%0%

Share and Cite

MDPI and ACS Style

Auguste, C.; Marsh, P.; Nader, J.-R.; Cossu, R.; Penesis, I. Towards a Tidal Farm in Banks Strait, Tasmania: Influence of Tidal Array on Hydrodynamics. Energies 2020, 13, 5326. https://doi.org/10.3390/en13205326

AMA Style

Auguste C, Marsh P, Nader J-R, Cossu R, Penesis I. Towards a Tidal Farm in Banks Strait, Tasmania: Influence of Tidal Array on Hydrodynamics. Energies. 2020; 13(20):5326. https://doi.org/10.3390/en13205326

Chicago/Turabian Style

Auguste, Christelle, Philip Marsh, Jean-Roch Nader, Remo Cossu, and Irene Penesis. 2020. "Towards a Tidal Farm in Banks Strait, Tasmania: Influence of Tidal Array on Hydrodynamics" Energies 13, no. 20: 5326. https://doi.org/10.3390/en13205326

APA Style

Auguste, C., Marsh, P., Nader, J. -R., Cossu, R., & Penesis, I. (2020). Towards a Tidal Farm in Banks Strait, Tasmania: Influence of Tidal Array on Hydrodynamics. Energies, 13(20), 5326. https://doi.org/10.3390/en13205326

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop