Next Article in Journal
3D–3D Computations on Submerged-Driftwood Motions in Water Flows with Large Wood Density around Driftwood Capture Facility
Next Article in Special Issue
Risk Assessment of Water Resources and Energy Security Based on the Cloud Model: A Case Study of China in 2020
Previous Article in Journal
Seismic Fragility Evaluation of Simply Supported Aqueduct Accounting for Water Stop’s Leakage Risk
Previous Article in Special Issue
The Effects of Triclosan on Physiological and Photosynthetic Characteristics of Chlorella vulgaris
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Study on the Influence of Suspended Matter by the Foundation Construction of Different Offshore Wind Turbines in the Taiwan Sea Area

1
Department of Harbor and River Engineering, College of Engineering, National Taiwan Ocean University, Keelung 20224, Taiwan
2
Department of Shipping Technology, College of Maritime, National Kaohsiung University of Science and Technology, Kaohsiung 805, Taiwan
3
Center for General Education, College of Education & Communication, Tzu Chi University, Hualien 97004, Taiwan
*
Author to whom correspondence should be addressed.
Water 2021, 13(10), 1405; https://doi.org/10.3390/w13101405
Submission received: 31 March 2021 / Revised: 10 May 2021 / Accepted: 13 May 2021 / Published: 18 May 2021
(This article belongs to the Special Issue Water Environmental System Analysis)

Abstract

:
The objective of this study was to adequately examine potential wave fields, flow fields, and suspended load changes in different wind turbine foundations. Accordingly, this study applied the hydrodynamic model to simulate waves, currents, and suspended load in the study area. The simulation results are based on the assumption that dredging and rubble bed trimming were performed for 8 h and that the per foundation setting operation was completed in 2 h. The influence on the tripile and jacket was larger than that on the monopile, and the influence time was longer. However, due to the influence of tidal currents on the sea, the suspended load also became more acceptable than the initial concentration. From a macroscopic perspective, the different foundations did not sufficiently affect the study area. From a microscopic perspective, changes in the suspended load were only limited to areas surrounding the piles after the installation of the wind turbines.

1. Introduction

Due to the increasing environmental awareness in recent years, people have begun to consider the problems of air pollution and nuclear waste treatment caused by fossil fuels and nuclear power generation, respectively. Therefore, determining economically beneficial and environmentally friendly power generation methods has become an imperative task for Taiwan. Over the past 10 years, Taiwan has been actively developing terrestrial wind farms. Currently, a total of 24 wind farms have been established on the island. Nevertheless, favorable locations for developing wind farms have been exhausted, and wind turbine generators because considerable noise engendered by wind shear; therefore, establishing offshore wind farms will become a future trend for Taiwan.
Taiwan is a nation surrounded by the ocean and has ideal conditions for the development of marine applications. In recent years, the needs for economic development, the improvement in marine clean energy technology, the rise in marine recreational activities, the development of new types of coasts and oceans, and the rise in marine environmental awareness highlight the problems of space use. Based on the issues of global warming and clean energy demand in recent years and under the pressure of abolishing nuclear energy and reducing carbon emissions, offshore wind power has become one of the options for Taiwan’s energy transformation. The offshore wind field in the Taiwan Strait has excellent resources. For example, 4C Offshore: Marine Consultants uses the observation data of the average wind speed to show that 13 of the best observation sites in the world are located on the western coast of Taiwan [1]. According to the selection results of offshore wind power planning site selection [2], the Changhua sea area in western Taiwan has a capacity of as high as 62.6% for offshore wind power installations, ranking first in Taiwan’s offshore wind power generation capacity distribution area.
The layout of offshore structures is closely related to sedimentation because the layout of engineering structures changes the sediment bed dynamic conditions of the original flow field and wave field. In order to avoid the influence of the wind turbine’s wake, the layout of the generator is generally set to use the distance of up to 4 times the parallel wind direction and the distance of 10 times the vertical wind direction as the interval. This arrangement is mainly used to avoid the kinetic energy loss caused by downstream wind turbines [3,4,5]. In addition, referring to the relevant environmental protection specifications of the Environmental Impact Assessment Inquiry System for the development of offshore wind farms in the sea area around Taiwan, with the future development of large-scale wind turbines, the offshore wind turbines planned by Taiwan’s offshore wind power developers are separated by about 1 km. The development utilized in the comprehensive planning of ocean space was carried out in the long term, prioritizing public interest regarding the use of the sea area for offshore wind power [6,7].
In order to take into account, the comprehensive evaluation of the use of space in the sea area, this study uses the hydrodynamic model under the influence of wave–current interaction to simulate the wave–current pattern around the structure. In addition, studies on the impact of offshore wind farms on marine ecology have mostly focused on the impact of fish, birds, and cetaceans, and there are relatively few studies on the impact on the marine environment and marine basic productivity. Among them, the three most important environmental factors affecting the primary productivity of the ocean are the abundance of light in the marine environment, the stability of water bodies, and the abundance of seawater nutrients. Due to the construction of offshore wind turbines which may increase suspended loads in the seawater and obscure the light in the sea, this effect may weaken the photosynthesis ability of phytoplankton in the seawater, leading to a decline in primary productivity. Accordingly, this study simulated the suspended load changes occurring after the installation of different wind turbine foundations. Numerical modeling was performed to examine the influence of the offshore wind farm set on its neighboring ocean environment. The results can serve as a reference for subsequent engineering construction activities and long-term geomorphological change analysis.

2. Methodology

2.1. Hydrodynamic Model

In this study, the tidal effect is caused by the tide, which is added to the hydrodynamic model. The governing equation of the two-dimensional hydrodynamic model can be derived through the following process [8,9,10,11,12]: first, the Navier–Stokes equation is simplified by assuming an incompressible fluid and hydrostatic pressure distribution. The conservation of mass and momentum equations are then acquired on the basis of a depth integral equation and appropriate boundary conditions:
η t + x U h + η + y V h + η = 0
U t + U U x + V U y = f v g η x + 1 ρ τ x x x + τ y x y + 1 ρ h + η τ s x τ b x 1 ρ h + η S x x x + S y x y
V t + U V x + V V y = f v g η x + 1 ρ τ x y x + τ y y y + 1 ρ h + η τ s y τ b y 1 ρ h + η S x y x + S y y y
where η is the water surface elevation (m) [13]; h represents the distance from the static water level to the bed (m); g is the gravitational constant (m/s2); U and V represent the mean current velocity (m/s) of the water depth in the fixed coordinates of the x and y axes:
U = 1 h + η h η u d z , V = 1 h + η h η v d z
Shear stresses τ x x , τ x y , τ y x , and τ y y include viscous stress caused by fluid viscosity and Reynold’s stress caused by turbulent effects. As the value of viscous stress compared with Reynold’s stress is very small, viscous stress is ignored generally, and only Reynolds stress is considered to represent the momentum exchange between fluids:
τ x x = ρ E v U x , τ x y = ρ E v U y , τ y x = ρ E v V x , τ y y = ρ E v V y
The vortex viscosity coefficient E v is obtained from the semi-empirical formula of the Prandtl mixing length theory [14]:
E v = k v g d + h U 2 + V 2 6 C c
The sea surface wind shear components τ s x and τ s y are the components of the sea surface wind shear in the x and y directions [15]:
τ s x = ρ k w W 2 cos a ;   τ s y = ρ k w W 2 sin a
k w = 1.2 × 10 5 , W W c 1.2 × 10 6 + 2.25 × 10 6 1 W c W 2 , W > W c
The bottom friction stresses τ b x and τ b y are the components in the x and y directions [16]:
τ b x = ρ E r U U 2 + V 2 ; τ b y = ρ E r V U 2 + V 2
where the coefficient of the bottom friction is F r = g / C c 2 .
S x x , S x y , S y x , and S y y represent the components of the radiation stress caused by waves (kg/ms2), which can be estimated using the method proposed by [8]:
S x x S x y S y x S y y = ρ g H a 2 g n 1 + cos 2 θ 1 2 n 2 sin 2 θ n 2 sin 2 θ n 1 + sin 2 θ 1 2
where E ¯ is the total wave energy per unit time and area of section. Under airy wave theory,
E ¯ = ρ g H 2 8
The boundary conditions of the hydrodynamic model are shown in Figure 1. The water level change includes the level rise and fall caused by waves and tides.
The tide is a sine function changing from right to left, T L is the time difference when the tide reaches the left and right boundaries, and T s t a r t is the phase difference:
Left boundary:
η L = A L t sin 2 π T t t + T t + T s t a r t , T t = L y g h M a x
Right boundary:
η R = A R t sin 2 π T t t + T s t a r t
Offshore boundary:
η 0 = A R t + A L t + A R t N y j N y 1 sin 2 π T t t + T t N y j N y 1 + T s t a r t
The water level ξ ¯ caused by waves is based on [13], ignoring the reflection effect. Equation (15) is the water level descent outside the surf zone, and Equation (16) is the water level uplift in the surf zone:
ξ d ¯ = H 2 8 k sinh 2 k h cos θ 2 / 3
d ξ u ¯ d x = K d h d x ,   K = 1 1 + 8 / 3 γ 2
The left boundary of the velocity is shown in Equation (17), the right boundary is shown in Equation (18), the offshore boundary is shown in Equation (19), and the longshore boundary is shown in Equation (20).
U j = 1 = U j = 2 ,   V y j = 1 = 0
U j = N Y = U j = N Y 1 ,   V y j = N Y = 0
U j = N X = U j = N X 1 ,   U y j = N X = 0
U = 0 ,   V = 0
The stability of the hydrodynamics calculation must satisfy Δ t 2 Δ s / g h M a x , where Δ s is the grid size. The maximum value of two adjacent time steps is less than the allowable error with Equation (21), and then the calculation of the next time step can be performed.
M a x η i j k + 1 η i j k ε η η i j k , ε η = 0.0001 M a x U i j k + 1 U i j k ε U U i j k , ε U = 0.0001 M a x V i j k + 1 V i j k ε V V i j k , ε V = 0.0001

2.2. Wave Model

The data required to calculate radiation stress are the wave height, period, and direction, which can be obtained from various types of wave models addressing dissimilar physical problems. When both waves and tides are present, there is substantial diversity in both the spatial and temporal scales; thus, if waves and tides are considered to be on the same space–time scale, obtaining practical engineering applications for large sea areas is difficult. Accordingly, the calculation efficiency can be substantially improved if the force exerted by the wave field is regarded as steady over a certain period of time, thus ignoring instantaneous changes in wave motion. Considering refraction, diffraction, and wave breaking caused by the transmission of deep-water waves to shallow-water areas, this study’s calculation of wave patterns was carried out by using the mild slope equation of the current effect [17]:
D 2 φ D t 2 + U D φ D t C C g φ + σ 2 k 2 C C g φ = 0
where U is the ambient current, is the horizontal gradient operator, φ is the two-dimensional velocity potential, k is the wave number, C and C g are the phase and group speed of the waves, and σ is the dispersion relation given by σ 2 = g k tan k h . Under the assumption of an irrotational field, with single-frequency linear surface waves, the potential energy of a wave can be expressed as follows:
φ x , y , z , t = f z , h φ x , y , t
where f z , h = cosh k h + z cosh k h .
In a single periodic harmonic motion, Equation (22) can be rewritten as follows:
φ x , y , t = Re a e i s e i ω t
The following expression can be obtained by substituting Equation (24) into Equation (22) for the real part and imaginary part:
1 a C C g U a U + U 1 a 2 a + 1 C C g C C g a k 2 + s 2 = 0
a 2 σ U + C g = 0
Equations (25) and (26) are the equations of motion for wave interactions before breaking waves. When the current velocity U is known, it solves the system of linear equations in two unknown parabolic simultaneous equations and obtains the amplitude a x , y and the wave numbers s . When U = 0 , Equations (25) and (26) become:
1 a 2 a x 2 + 2 a y 2 + 1 c c g a c c g + k 2 s 2 = 0
a 2 C C g s = 0
where a represents the wave amplitude (m) and S represents the phase function. RCPWAVE solves the mild slope equation of the parabolic type in a stable and fast manner with minimum computational calculation; hence, the model is reasonably applicable for making engineering application calculations for large sea areas.
In addition, energy is dissipated in the surf zone, and the energy expression of Equation (26) must be modified. Based on energy flux theory, ignore the effect of bottom friction [18]:
d E C g d x = ε ,   ε = 1 2 ρ V e k H B 2 ,   V e = V e B H B 2 c h B γ h B ,   V e B = 5 S B g 8 k B ρ 1 1 C 0 ,   S B = tan β 1 + 3 r 2 2
where c is the ratio of the radiation to the water depth of the recovery zone. According to [18], c = 0.17 when the wave recovery zone is not obvious in a gentle slope.
In the area of wave–current interaction, the energy dissipated by the nearshore current within the surf zone is small and negligible; therefore, the energy amplitude expression according to Equation (29) can be expressed as follows:
E σ U + C g = 5 16 ρ g 2 k B σ 2 tan β 1 + 3 r 2 2 1 1 C 0 H B 2 c h B r h B H B 2
With Equation (29), the energy in the surf zone is expressed, and Equation (30) is modified as follows in Equation (31).
a 2 σ U + C g = 2 g ρ E σ U + C g = 5 8 g 2 k B σ tan β 1 + 3 r 2 2 1 1 c r H B 2 c h B r h B H B 2
In Equations (29)–(31), the subscript B indicates the value at the surf zone. As the phase function of ϕ is x x , t = s x ω t , the wave number obtained from the deformed mild slope equation can be expressed as follows:
k = x = s
To obtain s from Equations (25), (26) or (31), it is necessary to know the direction of the wave. There are only two equations to solve a , s , θ . The linearity of the wave phase function gradient is assumed to be irrotational by Equation (15), and the convergence conditions of the wave model are given by Equation (16).
× s = 0 s = s cos θ i + s sin θ j x s sin θ y s cos θ = 0
H n o w H o l d ε H H n o w , ε H = 0.001 H 1 n o w H 1 o l d ε k H 1 n o w , ε k = 0.001 H 2 n o w H 2 o l d ε k H 2 n o w , ε k = 0.001

2.3. Particle Tracking Model

In this study, the particle tracking model is used to simulate the diffusion transmission. The model can be used to investigate the relationship between the influences of environmental factors on the suspended load and to understand the impact of changes in environmental factors caused by the construction of coastal structures on the transmission. The motion behavior of the simulated independent particle in the dynamic environment is an important tool for understanding the ocean pollution diffusion process, while the particle tracking model replaces the suspended load with the water quality point and treats the diffusion process of the independent particle as random motion [19]. This study refers to the research of [20,21] from the power point of view:
Q x = Q c u + U r , Q y = Q c v + V r
Q c = A 1 f c u + U r 2 + v + V r 2 + A 2 U H 2 U H C 2 g
U H = τ m ρ = f 2 U M a x , U M a x = π H T sinh k h
where u and v are the velocities of the current in the x and y directions caused by wave motion; U r and V r represent the speed of the source in the x and y directions, respectively; U H is the maximum shear velocity caused by water particles on the seabed under the action of waves; U M a x is the maximum velocity of water particles on the seabed under the action of waves; T is the wave period; H is the wave height; h is the water depth; f and f c are the friction coefficients of wave motion and average flow, respectively.
f = 0.00251 × exp 5.21 × exp A k s 0.19 , A k s > 1.57 0.3 , A k s 1.57
f c = g C c 2
where A b is the half width of the wave orbit; K s is the bottom roughness; U H C is the critical shear velocity of water particles on the bottom. U H C refers to [22]
U H C = 8.41 × d 50 11 / 32
where d 50 is the average particle size on the bottom, which is between 0.565 mm > d 50 > 0.065 mm ; A 1 is the sediment transport coefficient caused by the current; A 2 is the sediment transport coefficient caused by the wave motion.
In this particle tracking model, the vertical integral equations of the continuous equation and the equation of motion in the depth direction can be expressed as follows:
ζ t + H U x + H V y = 0
U t + U U x + V U y f V + g ζ x + g U U 2 + V 2 0.5 C 2 H 1 ρ H τ x s E x x 2 U x 2 + E x y 2 U y 2 = 0
V t + U V x + V V y + f U + g ζ y + g V U 2 + V 2 0.5 C 2 H 1 ρ H τ y s E y y 2 V y 2 + E y x 2 V x 2 = 0
In the two-dimensional mode, the position of each time step can be expressed as
X n = X n 1 + U Δ t ; Y n = Y n 1 + V Δ t
where n and n 1 each represent the old and new time steps; u and v are horizontal and vertical flow; and Δ t is unit time.

3. Model Set-Up

3.1. Modeling Procedure and Validation

Figure 2 displays a satellite image of the studied area, which was located at the offshore Changhua Coastal Industrial Park and had a water depth of −17 to −50 m. Figure 3 displays the calculation range of the model and the applied topographical data, which entailed a rectangular region with a length of 33 km and a width of 30 km, and the calculation ranged from the south of Changhua Coastal Industrial Park in the north to Zhuoshui River in the south. The topography data employed by the model comprised the data measured in 2014 and 200 m water depth data from the National Science Council. Geographic data were collected through the following procedures: after the distribution of the plane wave field was calculated, the mechanisms of radiation stress and tide-level variation were applied to derive data regarding mixed wave fields, including inshore currents. Table 1 lists the configuration settings of the hydrodynamic model. The model was verified as described below.
The currents simulated by the model generally flowed parallel to the coastline (Figure 4 and Figure 5). Specifically, the current flowed from the southwest to the northeast during flood tides, and from the northeast to the southwest during ebb tides. Figure 6 compares the model simulation results with the actual current measurement. The calculated inshore flow velocity, offshore flow velocity, flow direction, and changes in water level were all highly consistent with the measured data. The results, in accordance with the actual situation, indicate that the model successfully reproduced the characteristics of the current in the sea area adjacent to the study area. After the boundary conditions and related parameters adopted by the hydrodynamic model were determined, the model was employed to calculate the hydrodynamic and suspended load transport characteristics in different offshore wind turbines.

3.2. Suspended Load Characteristics at the Different Foundations

The case study of the offshore wind turbines in this area is shown in Figure 7, and it used monopile, tripile, and jacket foundations. In the monopile, the pile diameter is 6 m; the diameter is 3 m in the tripile foundation, and the interval between piles is 30 m; the diameter is 3 m in the jacket, and the interval between piles is 20 m. The calculated water depth near the wind turbine foundation is approximately 20 m, and the simulated conditions were simulated in summer marine conditions (wave height below 0.5 m). We assumed that dredging and rubble bed trimming were performed for 8 h and that the per foundation setting operation was completed in 2 h. We referred to the conditions that have a large impact on pollution in the foundation construction, such as trailing suction hopper dredgers (1500 m3/h), and carried out continuous construction for analysis under the abovementioned working hours. The simulation results are shown in Figure 8. Take the wind turbine foundation as the center, and take the simulation results in four directions with a radius of 50 m. As this study is based on the assumption that dredging and rubble bed trimming were performed for 8 h and that the per foundation setting operation was completed in 2 h, the degree of influence on the tripile and jacket was larger than that on the monopile, and the influence time was longer. In Figure 8, the maximum value of the monopile is 0.012 kg/s, for the tripile, it is 0.018 kg/s, and for the jacket, it is 0.017 kg/s. The period of time to complete the construction of the maritime engineering was the time period with the highest concentration of suspended load. After that, because there was no source, the suspended load gradually spread, and the influence was gradually reduced. In addition, due to the influence of tidal currents on the sea, the suspended load also became more acceptable than the initial concentration.

4. Conclusions

With different wind turbine foundations as a research subject, this study simulated the suspended load of summer conditions. The simulation indicated that the degree of influence on the tripile and jacket was larger than that on the monopile, and the influence time was longer. When the construction of the maritime engineering was completed, this period had the highest concentration of suspended load. After that, because there was no source, the suspended load gradually spread, and the influence was gradually reduced. In addition, due to the influence of tidal currents on the sea, the suspended load also became more acceptable than the initial concentration.
The validated model can be used to simulate the effect of suspended loads on offshore wind farms. In addition to attempting to understand the major suspended load mechanism and its impact within the studied area, subsequent studies can employ the current velocity and suspended sediment concentration during different return periods for scenario simulation. These should be implemented according to the analysis of trends in the targeted area under various circumstances.

Author Contributions

Conceptualization, H.-Y.W. and Y.-C.C.; methodology, Y.-C.C. and S.-S.H.; software, Y.-C.C. and H.-Y.W.; validation, S.-S.H., H.-Y.W. and Y.-C.C.; formal analysis, H.-Y.W.; investigation, H.-Y.W.; resources, Y.-C.C. and S.-S.H.; data curation, H.-Y.W. and Y.-C.C.; writing—original draft preparation, H.-Y.W. and Y.-C.C.; writing—review and editing, S.-S.H., H.-Y.W. and Y.-C.C.; visualization, Y.-C.C. and S.-S.H.; supervision, Y.-C.C. and S.-S.H.; project administration, Y.-C.C.; funding acquisition, Y.-C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 4C Offshore. Global Wind Speed Rankings. 2014. Available online: http://www.4coffshore.com/windfarms/windspeeds.aspx (accessed on 11 December 2014).
  2. Ministry of Science and Technology. Ocean Data Bank. 2018. Available online: http://www.odb.ntu.edu.tw/ (accessed on 20 September 2018).
  3. Elliott, D.L.; Barnard, J.C. Observations of wind turbine wakes and surface roughness effects on wind flow variability. Sol. Energy 1990, 45, 265–283. [Google Scholar] [CrossRef]
  4. Burton, T.; Jenkins, N.; Sharpe, D.; Bossanyi, E. Wind Energy Handbook, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2001; pp. 15–18. [Google Scholar]
  5. Jin, W.J. Numerical Simulation of Wind Turbine Wakes Based on Actuator Line Method in NEK5000. Master’s Thesis, Royal Institute of Technology, Stockholm, Sweden, 2013. [Google Scholar]
  6. Schomerus, T.; Runge, K.; Nehls, G. Strategic environmental assessment for the utilisation of offshore wind energy. In Proceedings of the 2nd Scientific Conference on the Use of Offshore Wind Energy by the Federal Ministry for the Environment, Berlin, Germany, 20–21 February 2007; pp. 57–73. [Google Scholar]
  7. Sivas, D.; Caldwell, M. A new vision for Californian ocean governance: Comprehensive ecosystem-based marine zoning. Stanf. Environ. Law J. 2008, 27, 209–270. [Google Scholar]
  8. Longuet-Higgins, M.S.; Stewart, R.W. Radiation in water waves: A physical discussion with applications. Deep-Sea Res. 1964, 11, 529–563. [Google Scholar] [CrossRef]
  9. Chiang, Y.C.; Hsiao, S.S.; Lin, M.C. A Quasi-3D Sediment Transport Modeling for Coastal Morphodynamics. In Proceedings of the 21st International Offshore and Polar Engineering Conference, Maui, HI, USA, 19 June 2011; pp. 1053–1058. [Google Scholar]
  10. Chiang, Y.C.; Hsiao, S.S.; Fang, H.M.; Wang, H.Y. A Study of The Sandbar Evolution by Tidal Current Effect. J. Coast. Ocean Eng. 2013, 13, 275–294. [Google Scholar]
  11. Chiang, Y.C.; Hsiao, S.S.; Fang, H.M.; Wang, H.Y. A Study of the Asymmetric Wave Parameterize. Int. J. Eng. Technol. 2017, 9, 239–243. [Google Scholar]
  12. Wang, H.Y.; Fang, H.M.; Hsiao, S.S. Morphological Characteristics of Tidal Inlets Subject to a Short Term Typhoon Event: A Case Study in Lanyan River Estuary. J. Mar. Sci. Technol. 2018, 26, 552–561. [Google Scholar]
  13. Matushevskiy, T.V. Method of separation of wave induced and turbulent motions in the near-surface layer of the ocean. In On Problems of Ocean Dynamics; Gidrometeoizdat: Leningrad, Russia, 1975; pp. 124–151. (In Russian) [Google Scholar]
  14. Falconer, R.A. Numerical modeling of tidal circulation in harbors. J. Waterw. Port Coastal Ocean Div. 1980, 106, 31–48. [Google Scholar] [CrossRef]
  15. Dean, R.G.; Dalrymple, R.A. Water Wave Mechanics for Engineers and Scientists; Prentice-Hall, Inc.: Hoboken, NJ, USA, 1984. [Google Scholar]
  16. Ponce, V.M.; Yabusaki, S.B. Modeling circulation in depth-averaged flow. J. Hydraul. Div. 1981, 107, 1501–1518. [Google Scholar] [CrossRef]
  17. Kirby, J.T. A parabolic equation for the combined refraction-diffraction of Stokes waves by mildly varying topography. J. Fluid Mech. 1983, 136, 453–466. [Google Scholar] [CrossRef] [Green Version]
  18. Mizuguchi, M. An heuristic model of wave height distribution in surf zone. In Proceedings of the 17th Coastal Engineering Conference, Sydney, Australia, 23–28 March 1980; pp. 278–289. [Google Scholar]
  19. Csanady, G.T. Turbulent Diffusion in the Environment; D. Reidel Publishing Company: Amsterdam, The Netherlands, 1973. [Google Scholar]
  20. Komar, P.D.; Inman, D.L. Longshore sand transport on beaches. J. Geophys. Res. 1970, 75, 5914–5927. [Google Scholar] [CrossRef]
  21. Watanabe, A.; Maruyama, K. Numerical modelling of nearshore wave field under combined refraction, diffraction and breaking. Coast. Eng. Jpn. 1986, 29, 19–39. [Google Scholar] [CrossRef]
  22. Sawaragi, T.; Deguchi, I. Distribution of sand transport rate across the surf zone. In Proceedings of the 16th Coastal Engineering Conference, Hamburg, Germany, 27 August–3 September 1978; Volume 2, pp. 1596–1613. [Google Scholar]
  23. Vries, W.E.D.; Vemula, N.K.; Passon, P.; Fischer, T.; Kaufer, D.; Matha, D.; Schmidt, B.; Vorpahl, F. Final Report WP 4.2: Support Structure Concepts for Deep Water Sites: Deliverable D4.2.8 (WP4: Offshore Foundations and Support Structures). 2011. Available online: https://www.researchgate.net/profile/Bjoern-Schmidt/publication/241886865_Final_report_WP_42_Support_Structure_Concepts_for_Deep_Water_Sites_Deliverable_D428_WP4_offshore_foundations_and_support_structures/links/5693cbb508ae3ad8e33b301a/Final-report-WP-42-Support-Structure-Concepts-for-Deep-Water-Sites-Deliverable-D428-WP4-offshore-foundations-and-support-structures.pdf (accessed on 15 April 2021).
Figure 1. Hydrodynamic model boundary condition diagram.
Figure 1. Hydrodynamic model boundary condition diagram.
Water 13 01405 g001
Figure 2. Satellite image showing the studied sea area.
Figure 2. Satellite image showing the studied sea area.
Water 13 01405 g002
Figure 3. Topographical settings of the model.
Figure 3. Topographical settings of the model.
Water 13 01405 g003
Figure 4. Current velocity and direction during flood tides.
Figure 4. Current velocity and direction during flood tides.
Water 13 01405 g004
Figure 5. Current velocity and direction during ebb tides.
Figure 5. Current velocity and direction during ebb tides.
Water 13 01405 g005
Figure 6. Comparison between the data obtained using the numerical model and the measured ocean currents: (a) tidal; (b) current velocity; (c) current direction.
Figure 6. Comparison between the data obtained using the numerical model and the measured ocean currents: (a) tidal; (b) current velocity; (c) current direction.
Water 13 01405 g006
Figure 7. Different foundation types: (a) monopile; (b) tripile; (c) jacket. The offshore wind foundation type images are from [23].
Figure 7. Different foundation types: (a) monopile; (b) tripile; (c) jacket. The offshore wind foundation type images are from [23].
Water 13 01405 g007
Figure 8. Simulation results of different foundation types: (a) monopile; (b) tripile; (c) jacket.
Figure 8. Simulation results of different foundation types: (a) monopile; (b) tripile; (c) jacket.
Water 13 01405 g008aWater 13 01405 g008b
Table 1. Calculation conditions for model calibration.
Table 1. Calculation conditions for model calibration.
AreaGrid SizeNumber of Grid Points
33 × 30 km50 × 50 m660 × 600
Coordinate of the origin (TWD67)Angle of deviation (counterclockwise from the north)Time step size
147892, 265658325°1.5 s
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hsiao, S.-S.; Wang, H.-Y.; Chiang, Y.-C. A Study on the Influence of Suspended Matter by the Foundation Construction of Different Offshore Wind Turbines in the Taiwan Sea Area. Water 2021, 13, 1405. https://doi.org/10.3390/w13101405

AMA Style

Hsiao S-S, Wang H-Y, Chiang Y-C. A Study on the Influence of Suspended Matter by the Foundation Construction of Different Offshore Wind Turbines in the Taiwan Sea Area. Water. 2021; 13(10):1405. https://doi.org/10.3390/w13101405

Chicago/Turabian Style

Hsiao, Sung-Shan, Hsing-Yu Wang, and Yun-Chih Chiang. 2021. "A Study on the Influence of Suspended Matter by the Foundation Construction of Different Offshore Wind Turbines in the Taiwan Sea Area" Water 13, no. 10: 1405. https://doi.org/10.3390/w13101405

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