Offshore wind farm impacts on surface waves and circulation in Eastern Lake Ontario

A coupled wave and hydrodynamic model was applied to the Kingston Basin of eastern Lake Ontario, a region with bathymetric variability due to channels and shoals, to assess the potential impacts on surface waves and wind-driven circulation of an offshore wind farm. The model was used to simulate a series of storm events with time-varying wind forcing and validated against wave, current and water level observations. The wind farm was simulated by adding semi-permeable structures in the surface wave model to represent the turbine monopiles, and by adding an energy loss term to the fluid momentum equations in the hydrodynamic model to represent the added drag of the monopiles on the flow. The results suggest that the wind farm would have a small influence on waves and circulation throughout the wind farm area, with spatial variability due to focussing of wave energy and re-direction of the flow. Overall, the results indicate that the wave height in coastal areas will be minimally affected with changes in significant wave height predicted to be < 3%. Larger changes to the strength of circulation occur inside the wind farm region with localized changes in current magnitude of up to 8 cm s^− 1. The results of this study may help to understand the impacts of future offshore wind farms and other offshore structures in the Great Lakes.     

Impact of parameterized lee wave drag on the energy budget of an eddying global ocean model

The impact of parameterized topographic internal lee wave drag on the input and output terms in the total mechanical energy budget of a hybrid coordinate high-resolution global ocean general circulation model forced by winds and air-sea buoyancy fluxes is examined here. Wave drag, which parameterizes the generation of internal lee waves arising from geostrophic flow impinging upon rough topography, is included in the prognostic model, ensuring that abyssal currents and stratification in the model are affected by the wave drag.An inline mechanical (kinetic plus gravitational potential) energy budget including four dissipative terms (parameterized topographic internal lee wave drag, quadratic bottom boundary layer drag, vertical eddy viscosity, and horizontal eddy viscosity) demonstrates that wave drag dissipates less energy in the model than a diagnostic (offline) estimate would suggest, due to reductions in both the abyssal currents and stratification. The equator experiences the largest reduction in energy dissipation associated with wave drag in inline versus offline estimates. Quadratic bottom drag is the energy sink most affected globally by the presence of wave drag in the model; other energy sinks are substantially affected locally, but not in their global integrals. It is suggested that wave drag cannot be mimicked by artificially increasing the quadratic bottom drag because the energy dissipation rates associated with bottom drag are not spatially correlated with those associated with wave drag where the latter are small. Additionally, in contrast to bottom drag, wave drag is a non-local energy sink.All four aforementioned dissipative terms contribute substantially to the total energy dissipation rate of about one terawatt. The partial time derivative of potential energy (non-zero since the isopycnal depths have a long adjustment time), the surface advective fluxes of potential energy, the rate of change of potential energy due to diffusive mass fluxes, and the conversion between internal energy and potential energy also play a non-negligible role in the total mechanical energy budget. Reasons for the <10% total mechanical energy budget imbalance are discussed.     

多方向不规则波传播变形数值模拟

在推广的缓坡方程数学模型基础上建立了多方向不规则波数学模型,综合考虑了波浪折射、绕射、反射、底摩擦和风能输入等因素。基于线性波浪理论,将波浪方向谱在频率和方向上按等能量分割法离散后,分别计算各组成波的传播变形,再计算合成波要素。缓坡方程数学模型采用改进的ADI法求解,计算效率高,稳定性好。采用椭圆形浅滩不规则波模型试验结果和单突堤不规则波绕射理论解对数学模型进行了验证,数值模拟结果和试验值及理论解符合良好。利用该模型进行了某港港内波浪折射、绕射和反射的联合数值模拟,给出了合理的港内波高分布。     

An Extended Mild-Slope Equation

On the assumption that the vortex and the vertical velocity component of the current aresmall,a mild-slope equation for wave propagation on non-uniform flows is deduced from the basichydrodynamic equations,with the terms of (V_hh)~2 and (V_h~2)h included in the equation.The terms of bot-tom friction,wind energy input and wave nonlinearity are also introduced into the equation.The wind en-ergy input functions for wind waves and swells are separately considered by adopting Wen′s(1989)empiri-cal formula for wind waves and Snyder′s observation results for swells.Thus,an extended mild-slope equa-tion is obtained,in which the effects of refraction,diffraction,reflection,current,bottom friction,wind en-ergy input and wave nonlinearity are considered synthetically.     

海上风电对海洋生态环境与海洋生物资源的综合影响研究进展

海上风电具有就近消纳方便、发电效率高和不消耗化石能源等特点,在低碳经济发展背景下,加快海上风电开发已成为全球各国促进能源结构转型与可持续发展的普遍共识。但海上风电在建设及运营过程中所产生的噪音和磁场对海洋环境和生物的影响尚不明确。本文系统梳理了全球海上风电发展现状,分析了海上风电开发对海洋生态环境与生物资源的综合影响,从生理、行为和分子三个层面重点分析了海上风电所产生的噪音和磁场对海洋生物的潜在影响,以期为科学利用海上风电提供参考。     

Experimental simulation of shoreface nourishments under storm events: A morphological,hydrodynamic, and sediment grain size analysis

This study focuses on barred beach shoreface nourishments physically simulated in a wave flume. The attack of a schematic storm on three different nourishments is analysed. The apex and waning storm phases lead respectively to offshore and onshore sediment transports. Nourishments in the trough and on the outer bar feed the bar and increase wave dissipation offshore. The bar acts as a wave filter and reduces shore erosion (lee effect). In contrast, nourishment on the beach face leads mostly to shore feeding and reconstruction (feeder effect). With successive nourishments, the beach face clearly becomes steeper and onshore sediment transport is reduced during moderate wave climates. The surface grain size analysis reveals marked variations. Coarser sediments are sorted on the bar and the upper beach face. These locations correspond to large wave dissipation zones during the storm apex.     

Depth- and current-induced effects on wave propagation into the Altamaha River Estuary, Georgia

A study of sea surface wave propagation and its energy deformation was carried out using field observations and numerical experiments over a region spanning the midshelf of the South Atlantic Bight (SAB) to the Altamaha River Estuary, GA. Wave heights on the shelf region correlate with the wind observations and directional observations show that most of the wave energy is incident from the easterly direction. Comparing midshelf and inner shelf wave heights during a time when there was no wind and hence no wave development led to an estimation of wave energy dissipation due to bottom friction with corresponding wave dissipation factor of 0.07 for the gently sloping continental shelf of the SAB. After interacting with the shoaling region of the Altamaha River, the wave energy within the estuary becomes periodic in time showing wave energy during flood to high water phase of the tide and very little wave energy during ebb to low water. This periodic modulation inside the estuary is a direct result of enhanced depth and current-induced wave breaking that occurs at the ebb shoaling region surrounding the Altamaha River mouth at longitude 81.23°W. Modelling results with STWAVE showed that depth-induced wave breaking is more important during the low water phase of the tide than current-induced wave breaking during the ebb phase of the tide. During the flood to high water phase of the tide, wave energy propagates into the estuary. Measurements of the significant wave height within the estuary showed a maximum wave height difference of 0.4 m between the slack high water (SHW) and slack low water (SLW). In this shallow environment these wave–current interactions lead to an apparent bottom roughness that is increased from typical hydraulic roughness values, leading to an enhanced bottom friction coefficient.     

基于数值模式评估细尺度参数化方案在南大洋背风波生成源地的适用性

验证基于GM(Garret-Munk)大洋内波普适谱的细尺度参数化方案在不同海域的适用性,对于湍流混合研究来说非常重要。包含背景GM内波场的高分辨率数值模式被用于评估细尺度参数化方案在背风波生成源地处的适用性。细尺度参数化方案主要包括基于剪切的G89 (Gregg 1989)方案、基于应变的W93 (Wijesekera 1993)方案、基于剪切和应变的GHP (Gregg-Henyey-Polzin)方案以及对GHP方案中的频率矫正项作出变换的IH (Ijichi-Hibiya)方案。计算结果显示,背风波的生成伴随着海底上方近惯性内波的增强,使得内波场的动能与势能的比值相较于GM内波场偏大。在这种情况下,基于剪切的G89方案会因为高估内波场的总能量而高估耗散率。反之,基于应变的W93方案会因为低估内波场的总能量而低估耗散率。计算结果还显示,已经考虑了内波谱变形的GHP方案仍然会高估耗散率,但IH方案能比较准确地估算耗散率。     

An approximate solution for the wave energy shadow in the lee of an array of overtopping type wave energy converters

In this study we investigate how the wave energy deficit in the lee of an array of overtopping type wave energy converting devices (WECs), redistributes with distance from the array due to the natural variability of the wave climate and wave structure interactions. Wave directional spreading has previously been identified as the dominant mechanism that disperses the wave energy deficit, reducing the maximum wave height reduction with increasing distance from the array. In addition to this when waves pass by objects such as an overtopping type WEC device, diffracted waves re-distribute the incident wave energy and create a complex interference pattern. The effect of wave energy redistribution from diffraction on the wave energy shadow in the near and far field is less obvious. In this study, we present an approximate analytical solution that describes the diffracted and transmitted wave field about a single row array of overtopping type WECs, under random wave conditions. This is achieved with multiple superpositions of the analytical solutions for monochromatic unidirectional waves about a semi-infinite breakwater, extended to account for partial reflection and transmission. The solution is used to investigate the sensitivity of the far field wave energy shadow to the array configuration, level of energy extraction, incident wave climate, and diffraction. Our results suggest that diffraction spreads part of the wave energy passing through the array, away from the direct shadow region of the array. This, in part, counteracts the dispersion of the wave energy deficit from directional spreading.     

An extended time-dependent numerical model of the mild-slope equation with weakly nonlinear amplitude dispersion

In the present paper, by introducing the effective wave elevation, we transform the extended ellip- tic mild-slope equation with bottom friction, wave breaking and steep or rapidly varying bottom topography to the simplest time-dependent hyperbolic equation. Based on this equation and the empirical nonlinear amplitude dispersion relation proposed by Li et al. (2003), the numerical scheme is established. Error analysis by Taylor expansion method shows that the numerical stability of the present model succeeds the merits in Song et al. (2007)’s model because of the introduced dissipation terms. For the purpose of verifying its performance on wave nonlinearity, rapidly vary- ing topography and wave breaking, the present model is applied to study: (1) wave refraction and diffraction over a submerged elliptic shoal on a slope (Berkhoff et al., 1982); (2) Bragg reflection of monochromatic waves from the sinusoidal ripples (Davies and Heathershaw, 1985); (3) wave transformation near a shore attached breakwater (Watanabe and Maruyama, 1986). Comparisons of the numerical solutions with the experimental or theoretical ones or with those of other models (REF/DIF model and FUNWAVE model) show good results, which indicate that the present model is capable of giving favorably predictions of wave refraction, diffraction, reflection, shoaling, bottom friction, breaking energy dissipation and weak nonlinearity in the near shore zone.     

SWAN模型中不同风拖曳力系数对风浪模拟的影响

本文以荷兰哈灵水道海域为实验区域,通过敏感性实验,研究了在14 m/s、31.5 m/s和50 m/s（分别代表一般大风、强热带风暴和强台风的极端条件）定常风速下SWAN模型中不同风拖曳力系数对风浪模拟的影响程度。结果表明,对于近岸浅水区域（水深小于20 m）,风拖曳力系数计算方案的选择对有效波高影响较小,而且当风速增加到一定程度后,波浪破碎成为影响波高值的主要因素;对于深水区域（水深大于30 m）,一般大风条件下风拖曳力系数计算方案的选择对有效波高影响仍然较小,随着风速的继续增大,风拖曳力系数计算方案的选择对有效波高的影响逐渐显著。对于平均周期,风拖曳力系数计算方案的选择和风速的改变对其影响均较小,而由水深变浅导致的波浪破碎对其影响较为显著。根据敏感性实验结果,本文对SWAN模型中风拖曳力系数计算方案的选择做出如下建议:计算近岸浅水区域风浪场或深水区域一般大风条件风浪场时,其风拖曳力系数可以直接采用模型默认选项;而对于深水区域更大风速条件,可首先采用模型默认选项试算,然后结合当地海域实测波浪资料进行修正。     

Spectral distribution of wave energy dissipation by salt marsh vegetation

Spectral energy dissipation of random waves due to salt marsh vegetation (Spartina alterniflora) was analyzed using field data collected during a tropical storm. Wave data (significant wave heights up to 0.4 m in 0.8 m depth) were measured over a two-day period along a 28 m transect using 3 pressure transducers. The storm produced largely bimodal spectra on the wetland, consisting of low-frequency swell (7–10 s) and high-frequency (2–4.5 s) wind-sea. The energy dissipation varied across the frequency scales with the largest magnitude observed near the spectral peaks, above which the dissipation gradually decreased. The wind-sea energy dissipated largely in the leading section of the instrument array in the wetland, but the low-frequency swell propagated to the subsequent section with limited energy loss. Across a spectrum, dissipation did not linearly follow incident energy, and the degree of non-linearity varied with the dominant wave frequency. A rigid-type vegetation model was used to estimate the frequency-dependent bulk drag coefficient. For a given spectrum, this drag coefficient increased gradually up to the peak frequency and remained generally at a stable value at the higher frequencies. This spectral variation was parameterized by employing a frequency-dependent velocity attenuation parameter inside the canopy. This parameter had much less variability among incident wave conditions, compared to the variability of the bulk drag coefficient, allowing its standardization into a single, frequency-dependent curve for velocity attenuation inside a canopy. It is demonstrated that the spectral drag coefficient predicts the frequency-dependent energy dissipation with more accuracy than the integral coefficient.     

海上风电升压站平台高程确定标准的合理性研究

为保证海上风电升压电站建设的经济合理与安全可靠,合理确定海上风电升压电站平台高程十分必要。文中从波浪与潮位的遭遇组合、最大波高取值与现行相关标准的比较、最大波峰高度计算的合理性等方面,全面分析了确定海上风电升压站平台高程各组成项取值标准的合理性,研究认为现行标准明显偏高。建议海上升压站平台底部高程按"100年一遇极端高水位+重现期50年波列累积频率1%的最大波峰高度+安全超高"确定。结合工程实例计算分析,按本文建议可使海上升压站平台高程明显降低,从而节省工程造价,还可减轻升压站工程对周边风机的遮蔽影响,以达到多发电量的效果。     

Examination of wind-wave interaction source term in WAVEWATCH III with tropical cyclone wind forcing

Results of drag coefficient(CD) from field observations and laboratory wave tank experiments indicate that the operational wave model can overestimate wind energy input under high wind conditions. The wind-wave interaction source term in WAVEWATCH Ⅲ has been modified to examine its behavior with tropical cyclone wind forcing. Using high resolution wind input,numerical experiments under idealized wind field and tropical cyclone Bonnie(1998) were designed to evaluate performance of the modified models. Both experiments indicate that the modified models with reduced CD significantly decrease wind energy input into the wave model and then simulate lower significant wave height(SWH) than the original model. However,the effects on spatial distribution of SWH,mean wavelength,mean wave direction,and directional wave spectra are insignificant. Due to the reduced wind energy input,the idealized experiment shows that the modified models simulate lower SWH than the original model in all four quadrants. The decrease in the front quadrants is significantly larger than that in the rear quadrants;it is larger under higher winds than lower winds. The realistic experiment on tropical cyclone Bonnie shows that the modified model with the various downward trends of CD in high winds creates a simulation that agrees best with scanning radar altimeter observations.     

On the interaction of a solitary wave and a submerged dike

The unsteady, two-dimensional Navier–Stokes equations and the exact free surface boundary conditions were solved to study the interaction of a solitary wave and a submerged dike. A piston-type wavemaker was set up in the computational domain to produce the incident solitary waves. The incident wave and the associated boundary layer flow in a wave tank with a flat bed were compared with the analytical solutions to verify the accuracy of this numerical scheme. Effects of the incident wave height and the size of the dike on the wave transformation, the flow fields, and the drag forces on the dike were discussed. Our numerical results showed that even though the induced local shear stress on the top surface of the dike is large at some particular locations, the resultant pressure drag is much larger than the friction drag. The primary vortex generated at the lee side of the dike and the secondary vortex at the right toe of the dike may scour the bottom and cause a severe problem for the dike.     

中咀湾避风港工程波浪分析及其应用

中咀湾是一个天然的避风良港,一般情况下外海波浪影响很小,主要受局部风场产生的局部风浪影响。本文采用曹宏生在Massel的扩展缓坡方程基础上推导出来的考虑陡变地形和能量耗散效应的缓坡方程为控制方程,结合固边界的反射边界条件,构成波浪传播变形的联合折射、绕射和反射的数学模型。文中将此方法运用在中国台州市大陈岛中咀湾避风港中,用波浪数学模型计算极端高水位和设计高水位时3种波况分别在3组重现期时工程海域的波浪要素,提供防波堤的堤前波高,并分析比较此处实心式和透空式防波堤的防浪性能。     

A Model for Nonlinear Irregular Waves in Coastal Area

Based on the refraction-diffraction theory of irregular waves in the waters of slowly-varying cur-rents and depths,and the generation dissipation theory of wind wave,a model for nonlinear irregular wavesin coastal area is developed.In light of the specific conditions of coastal wave character and engineering ap-plication,a practical mathematical model for the nonlinear irregular waves is presented.with directional spec-trum in coastal area.Coast effect,refraction,whitecapping.bottom friction.current,wind and nonlinearaction are considered in this model.The numerical methods and schemes for wave refraction ray,energy con-servation of propagation,energy balance of the generation and dissipation of wind waves have been studied.Finally,the model is used for the directional wave spectrum computation in the Daya Bay.Compared withthe measured data with 956 wave bouys in the Daya Bay,the model results are in good agreement with themeasured results.     

Numerical simulation and preliminary analysis of typhoon waves during three typhoons in the Yellow Sea and East China Sea

In this study,typhoon waves generated during three typhoons(Damrey(1210),Fung-wong(1416),and Chan-hom(1509))in the Yellow Sea and East China Sea were simulated in a simulating waves nearshore(SWAN)model,and the wind forcing was constructed by combining reanalyzed wind data with a Holland typhoon wind model.Various parameters,such as the Holland fitting parameter(B)and the maximum wind radius?,were investigated in sensitivity experiments in the Holland model that affect the wind field construction.Six different formulations were considered and the parameters determined by comparing the simulated wind results with in-situ wind measurements.The key factors affecting wave growth and dissipation processes from deep to shallow waters were studied,including wind input,whitecapping,and bottom friction.Comparison with in-situ wave measurements suggested that the KOMEN scheme(wind input exponential growth and whitecapping energy dissipation)and the JONSWAP scheme(dissipation of bottom friction)resulted in good reproduction of the significant wave height of typhoon waves.A preliminary analysis of the wave characteristics in terms of wind-sea and swell wave revealed that swell waves dominated with the distance of R to the eye of the typhoon,while wind-sea prevailed in the outer region up to six to eight times the R values despite a clear misalignment between wind and waves.The results support the hypothesis that nonlinear wave-wave interactions may play a key role in the formation of wave characteristics.     

Propagation of solitary waves over a submerged permeable breakwater

We study the interactions between a non-breaking solitary wave and a submerged permeable breakwater experimentally and numerically. The particle image velocimetry (PIV) technique is employed to measure instantaneous free surface displacements and velocity fields in the vicinity of a porous dike. The porous medium, consisting of uniform glass spheres, is mounted on the seafloor. Due to the limited size of each field of view (FOV) for high spatial resolution purposes, four FOVs are set in order to form a continuous flow field around the structure. Quantitative mean properties are obtained by ensemble averaging 30 repeated instantaneous measurements. The Reynolds decomposition method is then adopted to separate the velocity fluctuations for each trial to estimate the turbulent kinetic energy. In addition, a highly accurate two-dimensional model with the volume of fluid interface tracking technique is used to simulate an idealized volume-averaged porous medium. The model is based on the Volume-Averaged Reynolds Averaged Navier–Stokes equations coupled with the non-linear k–ε turbulence closure solver. Comparisons are performed between measurements and numerical results for the time histories of the free surface elevation recorded by wave gauges and the spatial distributions of free surface displacement with the corresponding velocity and turbulent kinetic energy around the permeable object imaged by the PIV system. Fairly good agreements are obtained. It is found that the measured and modeled turbulent intensities on the weather side are much larger than those on the lee side of the object, and that the magnitude of the turbulent intensity increases with increasing wave height of a solitary wave at a constant water depth. The verified numerical model is then used to estimate the energy reflection, transmission and dissipation using the energy integral method by varying the aspect ratio and the grain size of the permeable obstacle.     

海滩修复工程影响下的低能海岸波浪能量时空分布特征研究

本研究基于第三代海浪模式SWAN(Simulating Wave Nearshore),对茅尾海及其邻近海域波浪场进行了为期la的数值模拟,利用实测资料验证了该模型的可靠性.根据模型计算结果分析了茅尾海海域波浪要素的时空分布特征,在此基础上进一步探讨了波浪能量的输入耗散过程以及海滩修复对波浪能量空间分布的影响.研究发现...