• For Contributors +
• Journal Search +
Journal Search Engine
ISSN : 1229-3431(Print)
ISSN : 2287-3341(Online)
Journal of the Korean Society of Marine Environment and Safety Vol.25 No.7 pp.982-990
DOI : https://doi.org/10.7837/kosomes.2019.25.7.982

# A Study on the Tidal Energy Yield Capability according to the Yaw Angle in Jangjuk Strait

Bao Ngoc Tran*,Min Seon Choi**,Changjo Yang**†
* Doctoral Student, Graduate School of Mokpo National Maritime University, Mokpo 58628, Korea
** Professor, Division of Marine System Engineering, Mokpo National Maritime University, Mokpo 58628, Korea
* First Author : ngoctranbao.hn@gmail.com, 061-240-7472
Corresponding Author : cjyang@mmu.ac.kr, 061-240-7228
December 6, 2019 December 26, 2019 December 27, 2019

## Abstract

The interest of researchers and governments in exploiting tidal energy resources is increasing. Jangjuk strait is a place with high tidal energy density potential and is therefore appropriate for the constructing of a tidal turbine farm. In this study, a numerical approach is presented to evaluate the current flow and power potential in Jangjuk strait with an ADCIRC model. Then, the tidal field characteristics are utilized as input parameters for tidal resource calculation with an in-house program. The 1 MW scale tidal energy converter devices are employed and arranged in 4 layouts to investigate the annual energy yield as well as flow deficit due to the wake effect at the surveyed area. The best-performed array generates an annual energy yield up to 12.96 GWh/year (without considering the wake effect); this value is reduced by 0.16 GWh/year when accounting for the energy loss caused by the flow deficit. Moreover, by altering the turbine yaw angle during the flood and ebb tides, the impacts of this factor on the energy extraction are analyzed. This indicates that the turbine array attains the maximum tidal power when the turbine yaw angle is at 346o and 164o (clockwise, to the North) for the spring and neap tide in turns.

# 장죽수도에서의 요각변화에 따른 조류에너지 생산량에 관한 연구

뜨란 바우억*,최민선**,양창조**†
* 목포해양대학교 대학원 박사과정 , ** 목포해양대학교 기관시스템공학부 교수

## 초록

최근 점점 더 많은 연구자와 정부에서 해양에너지 자원 개발에 대한 관심이 고조되고 있다. 장죽수도는 조류에너지 밀도가 높 아 조류 발전소를 건설하기에 적합한 잠재적 후보지 중 하나이다. 따라서 본 연구에서는 ADCIRC 모델을 이용하여 장죽수도의 조류자원 의 잠재량을 평가하기 위한 수치적 접근방식을 제시하고, 내부 코드를 이용하여 조석 특성을 입력 매개변수로 활용하여 1 MW급 규모의 조류에너지 변환장치를 대상으로 4개의 레이아웃으로 배열하고 후류 효과로 인한 연간 에너지 생산량에 관한 수치해석을 수행하였다. 그 결과 효율이 가장 좋은 배치는 연간 최대 12.96 GWh/year의 에너지를 생성할 수 있으며, 이 값은 후류 효과로 인한 에너지 손실을 고려하 면 연간 0.16 GWh씩 감소될 수 있음을 보였다. 또한, 창낙조 때 터빈 요 각도를 변경함으로써 이 요소가 에너지 추출에 미치는 영향을 분 석하였으며, 터빈 배열은 터빈 요 각도가 346°와 164°(북쪽에서 시계 방향으로)일 때 대조기와 소조기에서 차례로 최대 조류 에너지를 얻 을 수 있었다.

Ministry of Oceans and Fisheries
20170333

## 1. Introduction

In our modern world, the demand for energy consumption has increasingly grown year by year, making a big challenge for energy sources since fossil fuel is becoming exhausted. A series of global issues such as climate change, global warming, and lacking energy impulse the attempt to reduce the use of fossil fuel and to invest more effort in new sources of renewable energy. Ocean-related sources including tidal energy, wave energy, and thermal energy provided a massive and perennial amount of power can be the feasible solutions for the thorny energy problem. Although there are still some weaknesses related to sediment accumulation and ecological change, it is undeniable that the advantages of the tidal energy are overwhelming the drawbacks. It is not only highly predictable relative to wind energy with higher rate of power absorbability (Esteban and Leary, 2012), more environmentally friendly due to fewer construction requirements for large-scale tidal turbines farmer but also less sensitive to the natural habitat of marine creatures.

The energy consumption of Korea is ranked the 8th largest amount in the world with 282 million tons of oil equivalent in 2016 (Korea Energy Agency, 2018). These statistics combined with the fact that Korea is a peninsula surround by sea by three sides and possess huge resources of tidal energy make them spend more progress on harvesting this potentially renewable source. To utilize the tidal resource, two important missions needed to be executed are to investigate the tidal current characteristics as well as the potentiality of tidal power at interested regions and to select the optimal configuration for tidal energy converters array which can maximize the power output of the energy projects. In recent years, several research applied different methods have been carried out to accomplish these above essential tasks. By using the experimental method, Byun (Byun et al., 2013) observed measured data at 264 stations to estimate the prospective tidal current energy resource of the south and west coast of Korea and categorized locations into six groups according to their annual power density. A study on the characteristics of tidal elevation and current in the Chusu Bay based on numerical modeling experiment was conducted (Jung et al., 2013). They found that in Chunsu Bay the size of the tidal residual currents was 1-7 cm/s in the northern part and 10-30 cm/s in the south with the highly complex distribution which largely depends on the geographic factors such as bathymetry and coastline. The numerical approach using tidal-related simulation software has also proved to be an alternative solution for tidal energy estimation. Ribal (Ribal et al., 2017) employed ADCIRC (Advanced Circulation) model and Green’s function to perform an assessment of tidal energy around Buton Island, Indonesia. Meanwhile, Martin investigated the flow characteristic of the tidal energy converters farm with the three-dimensional OpenFOAM model (Martin and Tao, 2018). Nguyen optimized the tidal farm around Jangjuk-sudo by testing 3 layouts of 28 tidal turbines with different longitudinal and lateral spacings. The energy yield, as well as wake pattern for each layout, was well predicted by the numerical modeling method (Nguyen et al., 2016).

Jangjuk strait was famous for its relatively high depth (25-45 m), medium current speed and was estimated to possess an annual energy density from 6.9 to 8.8 MWh/m2 (Byun et al., 2013;Ko et al., 2018). The interested area in our research is located between Jangjuk island and Jindo island which is remarked within the dark dot rectangular in Fig. 1. The tidal simulations are carried out in a region with the longitude between 126°04'00''E and 126°10'30''E and the latitude between 34°19'10''N to 34°22'30''N. In this paper, we present a numerical approach to analyze the tidal current flow at Jangjuk strait and then investigate the effect of tidal turbine arrangements and the turbine yaw angles on the energy yield capability. The former target is performed with a hydrodynamics model - ADCIRC which has been tested and intensively employed for tidal flow simulations in decade (Bacopoulos, 2006;Blain et al., 2002). While the influence of array layout and yaw angles are analyzed with the aid of a program relying on indispensable inputs including flow field data (tidal current’s direction and magnitude, ambient turbulence intensity), turbines’ characteristic, devices installation constraints, and model specifications. The annual energy yield of a tidal turbine array, as well as visualization of the wake effect due to devices mutual interactions, will be predicted. The flow speed during flood and ebb-tides obtained by the computational method will be compared with the gauged data at the real site provided by KHOA (Korea Hydrographic and Oceanographic Administration, 2010). After that, the energy yield and power loss are considered as key factors to determine the effects of turbines’ arrangement and yaw angles on the extracted power of the tidal turbine farm.

## 2. Methodology of numerical approach

### 2.1 Tidal flow field simulation with ADCIRC

In this study, the tidal elevation and the currents speed are simulated with a two-dimensional, depth-integrated, barotropic time-dependent longwave, hydrodynamic circulation model - ADCIRC. This model can be applied to computational domains encompassing the ocean, continental shelves, coastal seas, and small-scale estuarine systems. The model’s formulation is in the form of the Generalized Wave Continuity Equation based on the shallow water equations, hence solves the continuity and the momentum equations (with advection, Coriolis, eddy viscosity, and surface stress terms included). ADCIRC utilizes the finite element method in either 2D or 3D depth-integrated form on a triangular unstructured grid with Cartesian or spherical coordinates.

#### 2.1.2 Model set up

To prepare for simulations in ADCIRC, the data of coastlines, islands, and bathymetry must be imported to the model. When the geometrical data are available, the finite element grid can be generated with specific setting factors that are appropriate to domain size and computational accuracy. Fig. 2 displays the unstructured grid of the simulated area with depth integration. The mesh formation consists of 19,657 elements and 10,303 nodes in a calculation domain area of 85.8 km2. The size of the minimum and the maximum grids are 20 m (for islands) and 500 m (for the open ocean), respectively. Meanwhile, the coastline nodes’ spacing is about 50 m. In terms of boundary conditions for the ocean, Table 1 contains set up parameters employed in the model. Eight major tidal constituents are applied to Jangjuk strait based on the suggestions that they are fulfill the demand for calculating the tidal elevation and current speeds. The eight major tidal constituents are divided into two groups, the first one includes K1, O1, P1, Q1 for diurnal tides which are one cycle per day, the second group are K2, M2, N2, S2 for semidiurnal tides which are two cycles per day. Eight principal constituents and their periods are formed by the gravitational attraction between the earth and the moon and sun. The detail of the tidal constituents can be found in the study done by Militello (Militello and Zundel, 1999). The model runs in a period of 30 days from the 1st to 31st July 2019 with 1 day ram time and a time step of 5 seconds and the calculated results are obtained every 30 minutes.

### 2.2 Tidal potential energy assessment

#### 2.2.1 Calculation process

An in-house software is employed to provide the comprehensive and detailed assessment of the potential energy captured by tidal arrays. This software allows us to design a tidal array to achieve the maximum energy production within the geometric and environmental constraints of the site. It uses mathematical models to predict the flow throughout and around a tidal turbine array. Fig. 3 describes the required inputs to calculate the energy yield for a tidal turbine array, classified into 3 main components, namely: flow field data, tidal energy converter (TEC) features, and selected specifications in model setup. The data of the flow field including tidal current’s direction and intensity over the surveyed location can be attained with the hydrodynamics output of the ADCIRC model. Whereas the design parameters, as well as performance characteristics of the tidal turbines, are the results of the specific researches involved in tidal turbine’s blade design and behavior testing (Hoang and Yang, 2014).

Along with the provided information on the flow and tidal turbines, the program utilizes a proper wake model to determine the downstream wake deficit which considerably affects the energy absorbability of the downstream devices. Energy yield value, which is calculated by Equation 1, is a key output of the model, enabling the effect of wake losses and array efficiency to be evaluated for different array layouts.

$E A r r a y = ∑ k = 1 N D i r ∑ j = 1 N υ ∑ i = 1 N t P j k i ⋅ f j k$
(1)

where EArray is the gross energy output for the array, concerning the wake effect on the incident flow, Pjk and fjk is the power output and the frequency of occurrence, Nt is the total number of devices in the array and i is the device index, Nv is the number of flow speeds (bin) for each flow direction and j is the flow speed index related to the speed bin, NDir is the total number of flow speed directions and k is the index for flow direction.

#### 2.2.2 Initial inputs and set up parameters

In order to specify the flow resource in the hydrodynamic data and investigate the local potential energy yield, a series of initial conditions should be inputted, including the coordinate and depth of simulated area, long-term resource at a certain point, power law, and turbulence intensity profiles. The long-term resource is typically based on a harmonic model of the flow currents which is displayed in Fig. 4. It presents the principal directions and magnitude of the flood and ebb-tides at the long-term reference point (a point at the feasible turbines installation place). At the researched site, the dominant directions during the flood and ebb tides are 346° and 164° to the North in turn while the bin speeds range from 0.2 m/s to 3.2 m/s.

Besides, the parameters of turbines’ characteristics and constraints for device installation also play essential rules in the energy yield evaluation. In this study, a model of 1 MW power scale tidal stream turbine is employed. The tidal turbines with a rotor diameter of 22 m, rated velocity of 2.33 m/s are deployed near Jangjuk strait, obeyed the regulations in Table 2. The constraints should be defined with the purpose to make sure that the deployment of turbines at any location is feasible. Otherwise, the violated turbine will be switched off, unable to be considered in calculating the total energy yield of an array. Fig. 5 shows the depth information of the expected installing area and 4 layouts of turbines’ arrangement, namely layout A, B, C, and D. With the installation depth of around 30 m, four tidal turbines are arranged in 4 array configurations to investigate their extracted output power. Hence, it is able to analyze the impact of the wake deficit and turbines’ position to find out the best layout should be applied to the surveyed region.

## 3. Results and Discussion

### 3.1 Demonstration of the tidal flow field

Fig. 6 illustrates the tidal current distribution around Jangjuk island during the flood and ebb-tides. As can be seen, due to a narrow channel between two islands the flow is accelerated in both tides, accounting for a maximum depth-averaged velocity of 3.2 m/s. Besides, the geographic features at this site make the flow considerably change its direction compared to the surrounding area as mentioned in Fig. 4. The high current speed gives rise to a high potential energy density and consequently a higher extracted power. It induces this place to become a feasible location for the construction of a tidal turbine farm.

The ADCIRC model’s calculation is validated by a comparison between simulation results and measured data provided by KHOA at an observed point. The tidal gauge station is located at the coordinate 34°21'18''N, 126°4'44''E while the tidal simulations have been conducted for 1 month from 1st to 30th June 2019 as can be seen in Fig. 7. For each source of data, the results are recorded every 30 minutes. The maximum flow velocity of the observation data from KHOA is slightly higher than that of computed data with the ADCIRC model. Nevertheless, two lines show a good agreement in both speed magnitude as well as the tidal phase. This reveals the accuracy of the numerical model and ADCIRC outputs can be used as the initial conditions for energy yield calculation.

In the next step, the flow field information obtained from the ADCIRC model is used for resource calculation to investigate the available power across the site. The root-mean-cube flow speed and the energy density, which are two essential factors indicating the potential power, are depicted in Fig. 8. The root mean cube speed (Vrmc) is defined by Equation 2 while the energy density over a certain period time is estimated by Equation 3:

$V r m s = V 3 ¯ 3$
(2)

$E d = 1 2 ∫ 0 T ρ V 3 d t$
(3)

where Vrmc is the root mean cube speed (m/s); V is the flow velocity (m/s); T is the time of surveyed period (hour); Ed is the energy density (Wh/m2).

A close relationship between Vrmc and Energy density is obvious in Fig. 8 since the higher RMC speed is, the higher energy density is observed. The maximum energy density at the narrow channel is around 15 kWh/m2, being relatively good for tidal energy exploitation. In general, the numerical visualization of the tidal flow speed and tidal energyresource provides a quick effective method to evaluate the potential energy at the interested region. These useful functions facilitate the process of looking for the best location with high energy density for the construction of tidal turbine array, thus shorten the required time and progress compared to real site measurement.

Moreover, by conducting the energy yield calculation for different configurations of the turbine array, it is able to appreciate the impacts of wake effect on the captured energy each turbine. Fig. 9 illustrates the wake flow visualization of 4-turbine arrays arranged in 4 layouts in two directions of flood and ebb-tides. Theoretically, placing turbines at high flow velocity places results in more output power. Nevertheless, in reality, the turbines must be positioned in a manner that satisfies the installation constraints including maximum slope, minimum depth, minimum tip clearance, and lateral and longitudinal spacing. Fig. 10 shows the energy yield with and without the wake effect. The most significant feature is that layout C and layout D experience better performance than two others, both for gross energy and net energy. It may be caused by the fact that there are more turbines installed at higher power density region in layout C and layout D. As for the gross energy, layout C extracts the maximum value of 12.96 GWh/year while the minimum number belongs to layout B with 1.1 GWh/year less.

Nonetheless, due to a great loss brought about the wake effect of 3.4 % (presented in Table 3), the net energy captured by layout A is the smallest, comprising only 11.69 GWh/year. It reveals that turbines in layout A are the most impacted by the wake, especially device No. 2 with a remarkable loss of 0.199 GWh/year. In contrast, layout B is almost uninfluenced by the wake as the rate of loss is only 0.01 %. Concerning the power loss among turbines in an array, device No. 4 witnesses the least amount of energy deficit while device No. 2 suffers from the highest loss which means it is most vulnerable to the flow attenuation caused by other turbines in the array.

### 3.3 Impact of the yaw angle on the tidal energy absorbability

The angle between the turbine plane and the flow direction plays an important role in determining the volume of extracted power by dint of its effect on the incident angle to the blade. Therefore, the influence of the turbine yaw angle (defined as the angle between the North and turbine hub, clockwise) will be analyzed by evaluating the net energy yield at various yaw angles. Based on the results in section 3.2, layout C is selected to perform the investigation. For each case, the angle of turbine is changed either in flood-tide or ebb-tide. Fig. 11a and Fig. 11b express the energy yield each device when altering the yaw angle during flood and ebb-tides, respectively. A feature observed from the chart is that each turbine experiences individual trends and intensity of variation in output power. Device No. 4 seems most effected by yaw angle change whereas the impact on device No. 1 is negligible. In Fig. 12, the yaw angle of 346°-164° (flood-ebb) is the position at which the tidal array generates the highest power. At these angles, the turbines’ plane is perpendicular to the current direction. It is obvious that the net yield energy is in inverse proportion to the energy loss. For instance, turbine arrays with yaw angle 320° and 360° (spring-tide) possess the biggest losses of 1.67 % and 1.6 % in turn, hence they produce the lowest power production of 12.66 GWh/year. A similar point occurs to the array when varying the yaw angle during the neap-tide. The obtained energy peaks at the angle of 164° and gradually decreases even if the angle goes up or down.

## 4. Conclusion

In this paper, the tidal current conditions at Jangjuk strait was simulated to analyze its flow field behavior as well as the potential tidal energy resource. The efficiency of the tidal turbine arrangement and the influence of the yaw angle on energy yield are investigated to obtain some following points:

The area of interest possesses a number of favorable conditions including averaged-depth current speed of 3.2 m/s, sea depth around 30 m, showing high potential to build a small scale tidal turbine farm for harvesting tidal current energy.

The flow deficit and energy loss caused by the wake effect are visualized and presented. Among 4 layouts of turbine array, the stagger layout C performs the best with the net annual energy yield of 12.8 GWh/year while layout A generates the lowest output power.

The turbine yaw angle is one of the factors impacting on the amount of extracted energy. The energy production reaches the highest value at the yaw angle of 346° and 164° during flood-tide and ebb-tide where the turbine plane is perpendicular to the current.

## Acknowledgement

This study is financially supported by Mokpo National Maritime University and is a part of the project titled "Establishment of sea test-bed for tidal current energy converters", funded by the Ministry of Oceans and Fisheries, Korea (Project No: 20170333).

## Figure

Region of interest in Jangjuk strait.

Mesh formation for the region of interest.

Annual energy calculation for the tidal array.

Tidal dominant directions at the reference point.

Sea depth and layouts for turbine installation.

Tidal currents distribution.

Comparison of the current speed between numerical and measured data at the observation point.

Tidal resource at investigated area.

Wake pattern illustration during flood and ebb-tides in 4 cases of turbines’ arrangements.

Gross & Net energy yield for different layouts.

Energy yield absorbed by individual TEC.

Array annual energy yield at various yaw angles.

## Table

Tidal turbine's set up parameters

Energy loss due to wake effect

## Reference

1. Bacopoulos, P. (2006), Analysis, Modeling, and Simulation of the Tides in the Loxahatchee River Estuary (Southeastern Florida), University of Central Florida.
2. Blain, C. , R. Preller, and A. Rivera (2002), Tidal prediction using the advanced circulation model (ADCIRC) and a relocatable PC-based system. Oceanography, Vol. 15, No. 1, pp. 77-87.
3. Byun, D. S. , E. H. Deirdre, and W. J. Jeong (2013), Tidal Current Energy Resources off the South and West Coasts of Korea: Preliminary Observation-Derived Estimates, Energies, Vol. 6, pp. 566-578.
4. Esteban, M. and D. Leary (2012), Current developments and future prospects of offshore wind and ocean energy, Appl. Energy, Vol. 90, pp. 128-136.
5. Hoang, A. D. and C. J. Yang (2014), Design and performance evaluation of a 10kW scale counter-rotating wind turbine rotor, J. Korean Society of Marine environment & Safety, Vol. 20, No. 1, pp. 104-112.
6. Jung, K. Y. , Y. J. Ro, and B. J. Kim (2013), Characteristics of Tidal Current and Tidal Residual Current in the Chunsu Bay, Yellow Sea, Korea based on Numerical Modeling Experiments, J. Korean Society of Coastal and Ocean Engineers, Vol. 25, pp. 207-218.
7. Ko, D. H. , J. S. Park, and K. S. Lee (2018), Assessment of Tidal current energy potential at Uldolmok in the Southwestern coast of Korea, J. Coast. Res., Vol. 85, pp. 1301-1305.
8. Korea Energy Agency (2018), Status and forecast of national energy, In 2018 New & Renewable Energy, pp. 58-59.
9. Korea Hydrographic and Oceanographic Agency (2010), A Study on Development of Tidal Current Energy Map(I): Based on Measured Tidal Current Speed.
10. Martin, N. and L. Tao (2018), Three dimensional tidal turbine array simulations using OpenFOAM with dynamic mesh, Ocean Eng., Vol. 147, pp. 629-646.
11. Militello, A. and A. K. Zundel (1999), Surface-Water modeling system tidal constituents toolbox for ADCIRC, Coastal Engineering Technical Note, Vol 4, No. 21, pp. 1-8.
12. Nguyen, M. H. , B. G. Kim, and C. J. Yang (2016), Tidal Farming Optimization around Jangjuk-sudo by Numerical Modelling, J. Fluid Machinery, Vol. 19, No. 4, pp. 54-62.
13. Ribal, A. , K. Amir, S. Toaha, J. Kusuma, and K. Ddin (2017), Tidal Current Energy Resource Assessment Around Buton Island, Southeast Sulawesi, Indonesia, Int. J. Renewable Energy Research, Vol. 7, No. 2, pp. 857-865.