An improved two-step soil-structure interaction modeling method for dynamical analyses of offshore wind turbines

The detailed modeling of soil-structure interaction is often neglected in simulation codes for offshore wind energy converters. This has several causes: On the one hand, soil models are in general sophisticated and have many degrees of freedom. On the other hand, for very stiff foundations the effect of soil-structure interaction could often be discounted. Therefore, very simple approaches are utilized or the whole structure is assumed to be clamped at the seabed. To improve the consideration of soil-structure interaction, a six-directional, coupled, linear approach is proposed, which contains an implementation of soil-structure interaction matrices in the system matrices of the whole substructure. The aero-hydro-servo-elastic simulation code FAST has been modified for this purpose. Subsequently, a 5 MW offshore wind energy converter with pile foundation is regarded in two examples.     

Analytical investigation of hydrodynamic performance of a dual pontoon WEC-type breakwater

Based on the linear potential flow theory and matching eigen-function expansion technique, an analytical model is developed to investigate the hydrodynamics of two-dimensional dual-pontoon floating breakwaters that also work as oscillating buoy wave energy converters (referred to as the integrated system hereafter). The pontoons are constrained to heave motion independently and the linear power take-off damping is used to calculate the absorbed power. The proposed model is verified by using the energy conservation principle. The effects of the geometrical parameters on the hydrodynamic properties of the integrated system, including the reflection and transmission coefficients and CWR (capture width ratio, which is defined as the ratio of absorbed wave power to the incident wave power in the device width). It is found that the natural frequency of the heave motion and the spacing of the two pontoons are the critical factors affecting the performance of the integrated system. The comparison between the results of the dual-pontoon breakwater and those of the single-pontoon breakwater shows that the effective frequency range (for condition of transmission coefficient K_T < 0.5 and the total capture width ratio η_total > 20%) of the dual-pontoon system is broader than that of the single-pontoon system with the same total volume. For the two-pontoon system, the effective frequency range can be broadened by decreasing the draft of the front pontoon within certain range.     

Experimental investigation on the flow induced vibration of an equilateral triangle prism in water

A series of flow induced vibration (FIV) experiments for an equilateral triangle prism elastically mounted in a water channel are performed with different system stiffness at constant damping and mass. An amplitude variation coefficient is proposed to describe FIV stationarity in the present study. The FIV of the prism can be divided into three primary regions based on the amplitude and frequency responses, which are the vortex induced vibration (VIV) branch, the transition branch from VIV to galloping, and the galloping branch. The transition branch occurs at the reduced velocity in the range of 7.8 < U_r = U/(f_n,air·D) < 10.4, accompanied with a relatively rapid increase in amplitude and a precipitous drop in frequency and vibration stationarity. In addition, the reduced velocity where the transition region is initiated is independent of the system stiffness. The maximum amplitude reaches 3.17 D in the galloping branch. The ratio of the response frequency to the natural frequency of the prism in air remains locked to approximately 0.65 throughout the fully developed galloping branch. Large amplitude responses in an infinite range of flow velocities, excellent vibration stationarity and steady vibration frequencies, which are characteristics of the galloping of the prism, have a positive impact on improving energy conversion.     

Wave climate scatter performance of a cycloidal wave energy converter

A lift based cycloidal wave energy converter (WEC) was investigated using potential flow numerical simulations in combination with viscous loss estimates based on published hydrofoil data. This type of wave energy converter consists of a shaft with one or more hydrofoils attached eccentrically at a radius. The main shaft is aligned parallel to the wave crests and submerged at a fixed depth. The operation of the WEC as a wave-to-shaft energy converter interacting with straight crested waves was estimated for an actual ocean wave climate. The climate chosen was the climate recorded by a buoy off the north-east shore of Oahu/Hawaii, which was a typical moderate wave climate featuring an average annual wave power P_W = 17 kWh/m of wave crest. The impact of the design variables radius, chord, span and maximum generator power on the average annual shaft energy yield, capacity factor and power production time fraction were explored. In the selected wave climate, a radius R = 5 m, chord C = 5 m and span of S = 60 m along with a maximum generator power of P_G = 1.25 MW were found to be optimal in terms of annual shaft energy yield. At the design point, the CycWEC achieved a wave-to-shaft power efficiency of 70%. In the annual average, 40% of the incoming wave energy was converted to shaft energy, and a capacity factor of 42% was achieved. These numbers exceeded the typical performance of competing renewables like wind power, and demonstrated that the WEC was able to convert wave energy to shaft energy efficiently for a range of wave periods and wave heights as encountered in a typical wave climate.     

The Mediterranean response to different space–time resolution atmospheric forcings using perpetual mode sensitivity simulations

The Mediterranean basin features a semi-enclosed sea, where interactions and feedbacks between the atmosphere and the Sea at various temporal and spatial scales play a predominant role in the regional climate. This study analyzes the Mediterranean Sea response in sensitivity experiments conducted by driving the NEMO-MED12 oceanic model in perpetual mode with various atmospheric forcings, all produced by the WRF non-hydrostatic mesoscale atmospheric model, but differing by their resolutions: two horizontal resolutions (20 km at basin scale and 6.7 km in the North-Western [NWE] area) and two temporal resolutions (daily and three-hourly). The atmospheric fields available from August 1998 to July 1999 are in good agreement with estimates derived from satellite data. The heat budget of the Mediterranean Sea represents an heat loss of 5 W/m² and the annual freshwater budget is ?1.04 m, in agreement with climatologies. An increase in the spatial resolution in the NWE area modifies the modeled circulation from ?10% to +15% for the SST, from ?30% to +50% for the SSS, from ?10% to +30% for the MLD and from ?10% to +30% for the EKE in surface. The increase in the wind speed with a better chanelling by the land orography enhances in particular the oceanic convection process in the NWE area. On the other hand, the increase in the temporal resolution reduces the convection process, because of the diurnal restratification of the oceanic upper layer. It also reduces the surface parameters high-frequency variability, whereas it increases the EKE values in surface, due to the rapid response to the wind.     

Large capacity multi-float configurations for the wave energy converter M4 using a time-domain linear diffraction model

The moored three-float line absorber WEC M4 has been developed to optimise power capture through experiments and linear diffraction modelling. With the progression down wave from small to medium to large floats, the device heads naturally into the wave direction. The bow and mid floats are rigidly connected by a beam and a beam from the stern float is connected to the hinge point above the mid float for power take off (PTO). Increasing the bow to mid float spacing to be more than 50% greater than the mid to stern float spacing has been found to improve power capture. To increase power capture further and potentially reduce electricity generation cost the number of mid floats and stern floats is increased while maintaining a single bow float for mooring connection. The bow and mid floats still form a rigid body while the stern floats may respond independently. A time domain linear diffraction model based on Cummins method has been applied to configurations of 121, 123, 132, 133, and 134 floats where the numbers indicate the number of floats: bow, mid, stern. This shows how power capture is increased while response remains similar. We only consider uni-directional (long-crested) waves with narrow band width typical of swell. By considering scatter diagrams for various offshore sites capacities may range from 3.7 MW to 17.3 MW for the eight float system with a capacity factor of 1/3 while the cost of electricity assuming capital cost to be a fixed multiple of steel cost is reduced from that for the three-float system.     

Numerical investigation into the application of response conditioned waves for long-term nonlinear fatigue analyses of rigid hulls

Previous studies of response conditioned wave methods have been focused on their applicability to the prediction of extreme nonlinear wave-induced load effects. The results showed that theses methods can be used to accurately and efficiently predict the nonlinear short-term probability distributions for rigid hull responses. This has led us to investigate how response conditioned wave methods can be used for long-term nonlinear fatigue analyses, and with which accuracy this can be done. In this paper we present the results from our investigation. The studies were performed using a container vessel with a length between perpendiculars of 281 m. Calculations were done with a nonlinear strip theory method in which the hull of the vessel was assumed to be rigid. The most likely response wave (MLRW) method was used to condition the waves. Only head seas were considered. We found that the MLRW method under-predicted the long-term fatigue damage by 3%. The method, however, required a simulation time that was approximately three orders of magnitude less than that required for a conventional long-term nonlinear analysis based on random irregular waves. A preliminary investigation showed that due to lacking springing and whipping contributions the MLRW method under-predicted the fatigue damage for a flexible hull by approximately 50%. Several comments about a more accurate extension of the proposed method to flexible hulls are included.     

湛江某海上风电桩基局部冲刷特征研究

海上风电场桩基局部冲刷是工程设计与运行阶段的重要参数之一。基于湛江某海上风电场桩基3次现场局部冲刷实测数据,进行冲刷坑最大深度、冲刷坑半径和冲淤变化特征的分析与研究;根据桩基局部冲刷经验公式,采用工程海域实测海洋水文动力学数据进行最大冲刷深度与冲刷半径的计算,并进行公式计算值的对比与分 析。结果表明:桩基础在防冲刷设施的保护下,3次实测最大冲刷深度基本稳定为4.0 m,最大冲刷深度与桩径之比为0.57。而经验公式的最大冲刷深度与桩径之比均超过了1.1,说明桩基防冲刷设施取得了一定的效果,冲刷坑半径的计算值与现场实测值吻合较好。建议海上风电场在运行阶段进一步加强桩基冲刷坑监测与防护。     

The influence of sea state on formation speed of alongshore variability in surf zone sand bars

The formation time of alongshore morphological variability in surf zone sand bars has long been known to differ from one beach to the other and from one post-storm period to another. Here we investigate whether the type of sea state, i.e. distant swell waves or locally generated short period wind sea, affects the formation time of the emerging alongshore topographic variability.A numerical modeling approach is used to examine the emergence of alongshore variability under different shore-normal wave forcing. A research version of Delft3D, operating on the time-scale of wave groups, is applied to a schematised bathymetry with a single bar. The model is then used to investigate several wave scenarios, examining the impact of peak period, frequency spread and directional spread on the formation time of alongshore variability.Results show that an increase in wave period has a large effect, changing the formation time up to O (250%) in case the wave period is changed from a representative value for the Dutch coast (T_p ~ 5–6 s) to an Australian South East coast value (T_p ~ 10–12 s). In contrast, modifications in the directional and frequency spread of the wave field result only in a minor change in the formation time.Examination of hydrodynamics and potential sediment transport shows that the variations in formation time are primarily related to changes in the magnitude of the time-averaged flow conditions. Variations in the magnitude of very low frequency (f < 0.004 Hz) or infragravity (0.004 < f < 0.04 Hz) surf zone flow velocities do not affect the mean sediment transport capacity. Consequently the formation speed of patterns is primarily governed by positive feedback between mean flow and morphology, and low frequency flow fluctuations are of minor importance.These findings indicate that the development of alongshore topographic variability may be faster at swell dominated open coasts, primarily due to the occurrence of longer period swell. Also, at a given site, the arrival of a long wave period swell after a storm can accelerate the emergence of variability.     

A numerical study of the effect of hurricane wind asymmetry on storm surge and inundation

The influence of the asymmetric structure of hurricane wind field on storm surge is studied with five types of numerical experiments using a three-dimensional storm surge model. The results from the case of Hurricane Floyd (1999) show that Floyd-induced peak surge would have been much higher had the region of maximum wind rotated 30–90° counterclockwise. The idealized cases (the hypothetical hurricanes) with a wind speed asymmetry of 20 m s^?1 show that the peak (negative) surge varied from 4.7 to 6.0 m (?5 to ?5.7 m) or equivalent to ?8.8% and 16.3% (2.8% and ?10.4%) differences as compared to the control experiment. The area of flooding varied from 3552 to 3660 km². The results from two other idealized cases of varying degree of wind speed asymmetry further show that with decreasing (increasing) asymmetry of wind fields, the variations of peak surge and peak negative surge caused by the rotation of wind fields decrease (increase) accordingly. The results suggest that in storm surge simulations forced by winds derived from balanced models, considerable uncertainty in storm surge and inundation can result from wind asymmetries. This is true even if all other storm parameters, including maximum wind speed, the radius of maximum winds (storm size), minimum central pressure, storm translation speed, drag coefficient, and model settings (domain size and resolution) are identical. Thus, when constructing ensemble and probabilistic storm surge forecasts, uncertainty in wind asymmetry should be considered in conjunction with variations in storm track, storm intensity and size.     

Consumption and production on scales of a few days of inorganic carbon,nitrate and oxygen by the planktonic community: results of continuous measurements at the Dyfamed Station in the northwestern Mediterranean Sea (May 1995)

Continuous measurements between 0 and 200 m depth were performed every 2 h over two separate periods of four days at a station in the open northwestern Mediterranean Sea (Dyfamed Station) during the Dynaproc cruise in May 1995. Estimates of the daily variations in profiles of temperature, partial pressure of CO₂, oxygen, chlorophyll a and nutrients were obtained. The distributions of the various physical and chemical properties were clearly different during the two time series, which were separated by a period of 11 days during which a wind event occurred. The mean daily utilization or production due to biological processes of dissolved inorganic carbon (DIC), nitrate+nitrite and oxygen were calculated along isopycnals using a vertical diffusion model. Between the surface and about 20 m depth, DIC was consumed and O₂ released during the two time series while the nitrate+nitrite concentrations as well as supplies were zero. After the wind event, the O₂ : C : N ratios of consumption (or production) were, on average, near the Redfield ratios, but during the first time series, the C : N utilization ratio between 20 and 35 m was two to three times that of Redfield stoichiometry and the oxygen release was low. The integrated net community production (NCP) in terms of carbon was equivalent during the two time series, whereas the chlorophyll a biomass was twice as high, on average, during the first time series but did decrease. These results imply that the production systems were different during the two periods. The first time series corresponds to a period at the end of production, due to the nutrient depletion in the euphotic layer. The formation of degradation products of the living material in dissolved organic form is probably important as indicated by the high C : N utilization ratios. The second time series corresponds to a reactivation of the primary production due to the upward shift of nutrients after the wind event.     

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.     

An octo-generator for energy harvesting based on the piezoelectric effect

A renewable energy harvester using the piezoelectric effect is developed for the ocean tidal and wind flow. The harvester is made of connected driving blades to an octo-generator, which has a rotator with n blades and a stator attached by eight mass-spring-piston-cylinder-piezoelectricity devices. The resonance and force magnification are utilized to increase the power output of the harvester. A corresponding mathematical model is developed to calculate the root mean square of the generated electric power. The simulation results indicate that the generated power is largely enhanced when the near-resonant condition is established. The power increases with increases in the magnetic flux density, the large-to-small diameter ratio of the cylinder, the size of magnetic bar face, and decreases in the gap between two magnetic faces and the size of the piezoelectric bar face. A generated power of 5 kW is realized by the harvester working under an ocean tidal speed, V = 1.75 m/s, and its geometric and material properties of driving length L = 7.5 m, spring constant k_v = 65000 N/m, gap between the two magnets s = 0.0015 m, large to small diameter ratio of the cylinder z = 6, and magnetic flux density Br = 1.45 T.     

A simplified calculation method for axial cyclic degradation of offshore wind turbine foundations in clay

Abstract

Pile foundation is the most popular option for the foundation of offshore wind turbines. The degradation of stiffness and bearing capacity of pile foundation induced by cyclic loading will be harmful for structure safety. In this article, a modified undrained elastic–plastic model considering the cyclic degradation of clay soil is proposed, and a simplified calculation method (SCM) based on shear displacement method is presented to calculate the axial degradated capacity of a single pile foundation for offshore wind turbines resisting cyclic loadings. The conception of plastic zone thickness R_p is introduced to obtain the function between accumulated plastic strain and displacement of soil around pile side. The axial ultimate capacity of single piles under axial cyclic loading calculated by this simplified analysis have a good consistency with the results from the finite element analysis, which verifies the accuracy and reliability of this method. As an instance, the behavior of pile foundation of an offshore wind farm under cyclic load is studied using the proposed numerical method and SCM. This simplified method may provide valuable reference for engineering design.     

A re-appraisal of the total biomass and annual production of Antarctic krill

Despite much research on Euphausia superba, estimates of their total biomass and production are still very uncertain. Recently, circumpolar krill databases, combined with growth models and revisions in acoustics have made it possible to refine previous estimates. Net-based databases of density and length frequency (KRILLBASE) yield a summer distributional range of ～19×10⁶ km² and a mean total abundance of 8×10¹⁴ post-larvae with biomass of 379 million tonnes (Mt). These values are based on a standardised net sampling method but they average over the period 1926–2004, during which krill abundance has fluctuated. To estimate krill biomass at the end of last century we combined the KRILLBASE map of relative krill density around Antarctica with an acoustics-derived biomass estimate of 37.3 Mt derived for the Scotia Sea area in 2000 by the Commission for the Conservation of Antarctic Marine Living Resources (CCAMLR). Thus the CCAMLR 2000 survey area contains 28% of the total stock, with total biomass of ～133 Mt in January–February 2000. Gross postlarval production is estimated conservatively at 342–536 Mt yr^?1, based on three independent methods. These are high values, within the upper range of recent estimates, but consistent with the concept of high energy throughput for a species of this size. The similarity between the three production estimates reflects a broad agreement between the three growth models used, plus the fact that, for a given population size, production is relatively insensitive to the size distribution of krill at the start of the growth season. These production values lie within the envelope of what can be supported from the Southern Ocean primary production system and what is required to support an estimated predator consumption of 128–470 Mt yr^?1. Given the range of recent acoustics estimates, plus the need for precautionary management of the developing krill fishery, our net-based data provide an alternative estimate of total krill biomass.     

Marine dissolved organic matter: can its C : N ratio explain carbon overconsumption?

Carbon overconsumption, i.e. the consumption of inorganic carbon relative to inorganic nitrogen in excess of the Redfield ratio at the sea surface, was examined in relation to the dynamics of dissolved organic carbon and nitrogen (DOC and DON) in the northeast Atlantic. We observed the presence of N-poor dissolved organic matter (DOM) in surface water during summer, requiring the consumption of inorganic carbon and nitrogen in a ratio exceeding the Redfield ratio. The C : N ratio of bulk DOM is not only different from the Redfield ratio but also variable, i.e. no fixed conversion factor of C and N exists where DOM is important in C and N transformations. The existence of N-poor DOM is recognized as a feature typical of oligotrophic systems. At the same time, the C : N ratios of particles conform to Redfield stoichiometry as does deep-ocean chemistry. The implications of this finding are discussed, the conclusion being that, while DOM buildup contributes to CO₂ drawdown seasonally, its impact on long-term carbon and nitrogen balance of the ocean is small.     

Research on truncation method of FPSO and offloading system in model test

The model test method of the FPSO and offloading system is investigated by using the development mode of “FPSO + CALM + TANKER” working in a 1700-m depth of offshore West Africa. An equivalent design based on static and dynamic similarity criteria for oil offloading line (OOL) is discussed, and a type of creative method for the equivalent design of OOL in a model test is proposed. Based on the static similarity criterion, the truncated design of the FPSO mooring system in water depth and horizontal directions is carried out. After that, a relevant static optimization is conducted. Meanwhile, to avoid interference between the FPSO mooring system and CALM mooring system, a horizontal equivalent design for the CALM mooring system is provided. On this basis, the model test scheme is conducted. Time domain coupled analyses for the whole system before and after truncation are later performed. After comparison, it is observed that the calculated results of the truncated system are basically consistent with those of the prototype system, and the design of the model test scheme is demonstrated to be robust and reliable.     

Experiments with flexible shrouds to reduce the vortex-induced vibration of a cylinder with low mass and damping

Experiments employing a low-mass-damping cylinder have been conducted to determine the vortex-induced vibration (VIV) response of four suppressors of the flexible-shroud family. The VIV suppressors were inspired in the concept of the Ventilated Trousers (VT), a flexible shroud composed of a flexible net fitted with three-dimensional bobbins. Reynolds number varied between 5 × 10³ and 25 × 10³, while reduced velocity varied from 2 to 26. The VIV dynamic response showed that the VT suppressed the peak amplitude of vibration down to 40% of that of a bare cylinder. Other flexible shrouds also achieved suppression, but not as efficiently. Drag was reduced during the VIV synchronization range, but remained above the value for a bare static cylinder thereafter. Spectral analysis of displacement and lift revealed that, depending on the geometry and distribution of the bobbins, the flexible shroud can develop an unstable behavior, capturing energy from the wake and sustaining vibrations for higher reduced velocities. PIV measurements of the wake revealed that the entrainment flow through the mesh is necessary to extend the vortex-formation length of the wake; this mechanism only occurs for the VT mesh.     

Flow induced vibrations of a sharp edge square cylinder in the wake of a circular cylinder

The response of an oscillating circular cylinder at the wake of an upstream fixed circular cylinder was classified by different researchers as galloping, wake induced galloping or wake induced vibration. Furthermore it is already known that a sharp edge square cylinder would undergo galloping if it is subjected to uniform flow. In this study the influence of the wake of a fixed circular cylinder on the response of a downstream square cylinder at different spacing ratios (S/D = 4, 8, 11) is experimentally investigated. The subject appears not to have received previous attention. The lateral displacements, lift forces and the pressure data from gauges mounted in the wake of the oscillating cylinder are recorded and analyzed. The single degree of freedom vibrating system has a low mass-damping parameter and the Reynolds number ranges from 7.7 × 10² to 3.7 × 10⁴.In contrast to that for two circular cylinders in tandem arrangement, the freely mounted downstream square cylinder displays a VIV type of response at all spacing ratios tested. There is no sign of galloping or wake induced galloping with the square cylinder. With increase at the spacing ratio the cross-flow oscillations decrease. It is shown that the vortices arriving from the upstream fixed circular cylinder play a major role on the shedding mechanism behind the downstream square cylinder and cause the square cylinder to shed vortices with frequencies above Strouhal frequency of the fixed square cylinder (St = 0.13). The VIV type of oscillations in the downstream square cylinder is most probably caused by the vortices newly generated behind the square cylinder.     

The effect of yaw angle on VIV suppression for an inclined flexible cylinder fitted with helical strakes

Experimental studies were carried out to investigate the response features of an inclined flexible bare cylinder as well as a straked cylinder in a towing tank, with the main purpose of further improving the understanding of the effect of yaw angle on vortex-induced vibration (VIV) suppression. Four yaw angles (a = 0°, 15°, 30°, 45°), which is defined as the angle between the cylinder axis and the plane orthogonal to the oncoming fluid flow, were tested. The cylinder model was towed along the tank to generate a uniform fluid flow. The towing velocity was in the range of 0.05–1.0 m/s with an interval of 0.05 m/s. The corresponding Reynolds number ranged from 800 to 16000. The strakes selected for the experiments had a pitch of 17.5D and a height of 0.25D, which is generally considered as the most effective configuration for VIV suppression of a flexible cylinder in water. The experimental results indicate that VIV suppression effectiveness of the inclined flexible straked cylinder is closely related to the yaw angle. The displacement amplitudes are significantly suppressed in both cross-flow (CF) and in-line (IL) directions at a = 0°. However, with increasing yaw angle, the suppression efficiencies of the CF and IL displacement amplitudes gradually decrease. In addition, the CF dominant frequencies of the straked cylinder obviously deviate from those of the bare cylinder at a = 0° and 15°. This deviation is substantially alleviated with increasing yaw angle. The IL dominant frequencies show less dependency on the yaw angle. Similar trends are also observed on the dominant modes of vibration and the mean drag coefficients.