Modelling of a wave energy converter array with a nonlinear power take-off system in the frequency domain

This paper presents a nonlinear frequency domain model and uses this to assess the performance of a wave energy converter (WEC) array with a nonlinear power take-off (PTO). In this model, the nonlinear PTO forces are approximated by a truncated Fourier series, while the dynamics of the WEC array are described by a set of linear motion equations in the frequency domain, and the hydrodynamic coefficients are obtained with the boundary element method. A single heave absorber is firstly investigated to establish the accuracy of the new model in capturing the nonlinear behaviour of the pumping system. Subsequently, simulations of a 2D array with 18 WECs and a pillar in the centre (representing the tower of a wind turbine) are carried out to understand wave interference effects. Several optimisation strategies are proposed to improve the overall performance of the WEC array. These results demonstrate a computationally effective method for accounting for nonlinear effects in large WEC arrays. The proposed approach may potentially be applied for developing control algorithms for the adaptability of a 2D array to incoming wave excitation.     

Added mass and damping of a sphere section in heave

The Norwegian wave-power buoy¹ consists of a half-immersed floating sphere which is open to the sea at the bottom end. It is a two degree-of-freedom device involving the independent motion of the outer rigid sphere and the pressure across the internal free surface. A simpler model of the device is to represent the oscillatory flow through the bottom opening as another rigid body motion being that of the curved surface which would complete the sphere.The wave-induced forces on this surface and also the outer spherical surface due to independent oscillations of either are determined semi-analytically using a simple extension of the method of Havelock,² recently simplified and generalised by Hulme,³ for the complete half-immersed sphere. In particular the 2 × 2 added mass and damping matrices are determined as a function of frequency and relative size of bottom opening to sphere radius. These quantities are essential in any theoretical analysis of a multi-degree of freedom wave-energy device.     

Optimal configurations of wave energy device arrays

In this paper, the influence of the spatial configuration of a wave energy device array upon total power output is investigated. Hydrodynamic interactions are computed using a method capable of producing the linear wave theory solution to arbitrary accuracy. The overall performance of devices with two different power take-off arrangements is maximised at one incident wave frequency and direction by altering the formation of the array. Minimisation of the power is also carried out in a third case in order to demonstrate potential array-related losses. The optimisation is applied using two different approaches in each case: the Parabolic Intersection (PI) method and a Genetic Algorithm (GA). The former is a heuristic technique that has been devised for this study to enable rapid array construction using only simple calculations. The latter is an existing method, applied here with a novel crossover operator. Although considerably more computational effort is required, superior results may be obtained using the GA compared to the PI method. All of the arrays are subsequently analysed under incident waves of different frequency and direction, the resulting behaviour explained in terms of certain geometrical features of the arrangements.     

Design optimisation of Propeller Boss Cap Fins for enhanced propeller performance

Economic pressures and regulatory requirements have brought about a great interest in improving ship propulsion efficiency. This can be exercised by installing Energy Saving Devices (ESD) such as Propeller Boss Cap Fins (PBCF). This paper demonstrates an approach for optimising PBCF by using Computational Fluid Dynamics (CFD) analysis. The conducted Reynolds-averaged Navier-Stokes (RANS) CFD open water model tests were validated by comparison with experimental data until the simulation was deemed satisfactory within the capabilities and limitations of the model. A design and optimisation procedure was defined to analyse the impact of ESDs on propeller efficiency and then used to evaluate the influence of alternative geometric parameters and locations of the PBCF on the hub. This analysis was done at full scale using high fidelity CFD-based RANS methods. Outcomes of the study include a design and optimisation process that can be used for the analysis of other ESDs on the market. The influences of various PBCF geometry were examined with optimal solutions presented for the analysis case. Results indicated a net energy efficiency improvement of 1.3% contributing to a substantial minimisation of cost and energy consumption. A reduction in the hub vortex was also clearly identified and presented.     

点吸收式波浪能发电装置姿态特性分析

点吸收式波浪能发电装置是一种最简单的振荡体式波浪能发电装置, 但其安装成本高、生存能力较差。本文针对点吸收式波浪能发电装置的姿态稳定性问题, 开展了其在波浪作用下的运动姿态和发电功率之间的关系研究。首先介绍了点吸收波浪能发电装置的工作原理; 然后,根据我国南海海域的自然资源条件, 划定波况范围, 利用相似理论在实验室中模拟波浪参数,选定工况, 建立模型, 设计测量系统, 开展物理模型试验; 最后, 利用试验结果分析了发电装置的最佳发电周期、波高对装置发电功率的影响, 装置姿态对发电功率的影响等。本文为点吸收式波浪能发电装置设计及测试提供了参考。     

Optimisation of an idealised ocean model,stochastic parameterisation of sub-grid eddies

An optimisation scheme is developed to accurately represent the sub-grid scale forcing of a high dimensional chaotic ocean system. Using a simple parameterisation scheme, the velocity components of a 30 km resolution shallow water ocean model are optimised to have the same climatological mean and variance as that of a less viscous 7.5 km resolution model. The 5 day lag-covariance is also optimised, leading to a more accurate estimate of the high resolution response to forcing using the low resolution model.The system considered is an idealised barotropic double gyre that is chaotic at both resolutions. Using the optimisation scheme, we find and apply the constant in time, but spatially varying, forcing term that is equal to the time integrated forcing of the sub-grid scale eddies. A linear stochastic term, independent of the large-scale flow, with no spatial correlation but a spatially varying amplitude and time scale is used to represent the transient eddies. The climatological mean, variance and 5 day lag-covariance of the velocity from a single high resolution integration is used to provide an optimisation target. No other high resolution statistics are required. Additional programming effort, for example to build a tangent linear or adjoint model, is not required either.The focus of this paper is on the optimisation scheme and the accuracy of the optimised flow. However the forcing can provide insights in the design of deterministic and stochastic parameterisations. In the present study, we found that the stochastic parameterisation correcting the model variance is associated with the spatial pattern of eddy-decorrelation timescales rather than the spatial pattern of the amplitude of the variance. The method can be applied in future investigations into the physical processes that govern barotropic turbulence and it can perhaps be applied to help understand and correct biases in the mean and variance of a more realistic coarse or eddy-permitting ocean model. The method is complementary to current parameterisations and can be applied at the same time without modification.     

Limitations on the overlap-correlator method imposed by noise andsignal characteristics

Limitations on the performance of the overlap-correlator method of forming a passive synthetic aperture are derived. The technique uses the overlap of the array in sequential positions to estimate a series of phase correction factors that compensate for the motion of the array over time. It is of primary interest to optimize this overlap with respect to the effects of random noise. By minimizing the variance of the estimates of the set of phase correction factors, it is found that the optimal overlap is one-half the length of the physical array. Using this optimal overlap, the bounds on the usable spatial response are then determined as a function of signal-to-noise ratio and the number of hydrophones in the physical array. The ability of the overlap-correlator algorithm to synthesize a coherent aperture is investigated for the case of multiple sources in the absence of noise     

基于CSA-PLS算法的养殖水体水质快速高光谱预测反演模型研究

养殖水体水质的优劣直接影响养殖对象的成长,准确、快速、全面地掌控养殖水环境的水质参数变化情况具有重要意义。传统的水质指标监测方法都通过人工采样的方式,不仅耗费时间长,且只能体现局部水体情况。针对这些问题,提出了一种乌鸦搜索算法(CSA)结合偏最小二乘回归(PLSR)的高光谱特征波段筛选方法,快速构建回归模型,实现光谱数据的精准预测反演。以连片的养殖小区为研究对象,采集养殖水体样本并拍摄同时期的高光谱影像数据。首先对提取的采样点光谱数据利用多种数据变换方法分别预处理;其次利用这些数据,对水质指标总氮(TN)、氨氮(NH₄⁺-N)、总磷(TP)和化学需氧量(COD)分别构建全波段的SVR和AdaBoost回归模型,同时与提出的CSA-PLS自动筛选波段方法和传统的连续投影算法(SPA)筛选波段后构建的模型进行比较分析;最后根据决定系数(R²)和均方根误差(REMS)选出适合各水质指标的最优模型。从实验结果可以看出,所提波段筛选方法的AdaBoost模型预测结果优于SVR和传统SPA方法提取特征波段后构建的模型,与全波段最优模...     

一种波浪能实验装置水动力学分析与优化设计

将造波水槽内二维浮体牵引弹簧回复液压缸的受力系统简化为弹簧—质量—阻尼器系统,建立数学模型,并根据牛顿第二定律得到运动方程式。采用基于简单格林函数的边界元方法对所研究浮体的水动力学系数和波浪力进行计算,对于施加给液压系统的不同外部阻尼值,由运动方程可得到相应的浮体垂荡运动位移。为求浮体对液压系统做功的最大值,在给定条件下着重对外部阻尼系数进行了优化。     

Effects of sensor placement on acoustic vector-sensor arrayperformance

We consider the role played by the sensor locations in the optimal performance of an array of acoustic vector sensors, First we derive an expression for the Cramer-Rao bound on the azimuth and elevation of a single far-field source for an arbitrary acoustic vector-sensor array in a homogeneous wholespace and show that it has a block diagonal structure, i.e., the source location parameters are uncoupled from the signal and noise strength parameters. We then derive a set of necessary and sufficient geometrical constraints for the two direction parameters, azimuth and elevation, to be uncoupled from each other. Ensuring that these parameters are uncoupled minimizes the bound and means they are the natural or “canonical” location parameters for the model. We argue that it provides a compelling array design criterion. We also consider a bound on the mean-square angular error and its asymptotic normalization, which are useful measures in three-dimensional bearing estimation problems. We derive an expression for this bound and discuss it in terms of the sensors' locations. We then show that our previously derived geometrical conditions are also sufficient to ensure that this bound is independent of azimuth. Finally, we extend those conditions to obtain a set of geometrical constraints that ensure the optimal performance is isotropic     

Optimal Active Control of Wave-Induced Vibration for Offshore Platforms

An obvious motivation of this paper is to examine the effectiveness of the lateral vibration control of a jacket type offshore platform with an AMD control device, in conjunction with H2 control algorithm, which is an optimal frequency domain control method based on minimization of H2 norm of the system transfer function. In this study, the offshore platform is modeled numerically by use of the finite element method, instead of a lumped mass model. This structural model is later simplified to be single-degree-of-freedom (SDOF) system by extracting the first vibration mode of the structure. The corresponding "generalized" wave force is determined based on an analytical approximation of the first mode shape function, the physical wave loading being calculated from the linearized Morison equation. This approach facilitates the filter design for the generalized force. Furthermore, the present paper also intends to make numerical comparison between H2 active control and the corresponding passive control using a T     

Near-optimal control of a wave energy device in irregular waves with deterministic-model driven incident wave prediction

This paper investigates wave-by-wave control of a wave energy converter using incident wave prediction based on up-wave surface elevation measurement. The goal of control is to approach the hydrodynamically optimum velocity leading to optimum power absorption. This work aims to study the gains in energy conversion from a deterministic wave propagation model that accounts for a range of group velocities in deriving the prediction. The up-wave measurement distance is assumed to be small enough to allow a deterministic propagation model, and further, both wave propagation and device response are assumed to be linear. For deep water conditions and long-crested waves, the propagation process is also described using an impulse response function (e.g. [1]). Approximate low and high frequency limits for realistic band-limited spectra are used to compute the corresponding group velocity limits. The prediction time into the future is based on the device impulse response function needed for the evaluation of the control force. The up-wave distance and the duration of measurement are then determined using the group velocity limits above.A 2-body axisymmetric heaving device is considered, for which power capture is through the relative heave oscillation between the two co-axial bodies. The power take-off is assumed to be linear and ideal as well as capable of applying the necessary resistive and reactive load components on the relative heave oscillation. The predicted wave profile is used along with device impulse response functions to compute the actuator force components at each instant. Calculations are carried out in irregular waves generated using a number of uni-modal wave spectra over a range of energy periods and significant wave heights. Results are compared with previous studies based on the use of instantaneous up-wave wave-profile measurements, both without and with oscillation constraints imposed. Considerable improvements in power capture are observed with the present approach over the range of wave conditions studied.     

Towed array shape estimation using Kalman filters-experimentalinvestigations

The problem of estimating the shape of a towed array instrumented with either depth sensors, compasses, or both in a discrete-time state-space formulation is treated in a companion paper by D. A. Gray et al. (to appear), in which the state-space representation is derived from a dynamical model of the propagation of tow-point-induced motion down the array. A Kalman filter is derived to recursively estimate the shape of this towed array, and solutions to the Riccati equation are used to predict the mean square error of the Kalman filter array shape estimates. The present study investigates the performance of this Kalman filter approach as an array shape estimator using both simulated examples and sea trial data. Fundamental to the Kalman filter approach is the model that describes the dynamical behavior of the towed array. The results of an experimental program that was undertaken to validate this model are also presented     

Evaluation of manoeuvring coefficients of a self-propelled ship using a blade element momentum propeller model coupled to a Reynolds averaged Navier Stokes flow solver

The use of an unsteady computational fluid dynamic analysis of the manoeuvring performance of a self-propelled ship requires a large computational resource that restricts its use as part of a ship design process. A method is presented that significantly reduces computational cost by coupling a blade element momentum theory (BEMT) propeller model with the solution of the Reynolds averaged Navier Stokes (RANS) equations. The approach allows the determination of manoeuvring coefficients for a self-propelled ship travelling straight ahead, at a drift angle and for differing rudder angles. The swept volume of the propeller is divided into discrete annuli for which the axial and tangential momentum changes of the fluid passing through the propeller are balanced with the blade element performance of each propeller section. Such an approach allows the interaction effects between hull, propeller and rudder to be captured. Results are presented for the fully appended model scale self-propelled KRISO very large crude carrier 2 (KVLCC2) hull form undergoing static rudder and static drift tests at a Reynolds number of 4.6×10⁶ acting at the ship self-propulsion point. All computations were carried out on a typical workstation using a hybrid finite volume mesh size of 2.1×10⁶ elements. The computational uncertainty is typically 2–3% for side force and yaw moment.     

锑原位掺杂SiO2@PDA-Sb复合材料修饰电极用于海水中Zn离子检测

本研究提出了一种用于检测海水中锌离子(Zn²⁺)含量的电化学检测方法。该方法首先应用水热法合成SiO₂@PDA-Sb复合材料,然后将该复合材料修饰到玻碳电极上制备SiO₂@PDA-Sb-Nafion/GCE,利用差分脉冲阳极溶出伏安法(DPASV)对珠江口海水水样中Zn²⁺进行测定。研究结果表明,Zn²⁺在SiO₂@PDA-Sb-Nafion/GCE上具有较强的溶出峰。在最优条件下,SiO₂@PDA-Sb-Nafion/GCE对Zn²⁺浓度在1～1 000 nmol/L范围内可实现灵敏、准确的检测,Zn²⁺的检测出限为0.71 nmol/L。加标回收率实验显示Zn²⁺加标回收率为93.19%～100.12%,表明该方法可应用于现场海水样品Zn²⁺测定。本方法具有电极制作简单、稳定性和抗干扰性良好,能够提高现有测定方法的检测限和精确度,在现场海水检测...     

The effect of temporal wave averaging on the performance of an empirical shoreline evolution model

The effect of using time-averaged wave statistics in a simple empirical model for shoreline change is investigated. The model was first calibrated with a six-year time series of hourly wave conditions and weekly shoreline position at the Gold Coast, Australia. The model was then recalibrated with the hourly waves averaged over intervals up to 1 year. With wave averaging up to 2 days, model performance was approximately constant (squared correlation r² ~ 0.61–0.62), with only small changes in the values of empirical model parameters (e.g. the beach response coefficient c varied by less than 4%). With between 2 and 40 day averaging, individual storms are not resolved; model skill decreased only modestly (r² ~ 0.55), but c varied erratically by up to 40% of the original value. That is, optimal model coefficients depend on wave averaging, an undesirable result. With increased averaging (> 40 days) seasonal variability in the wave field is not resolved well and model skill declined markedly. Thus, temporal averaging of wave conditions increases numerical efficiency, but over-averaging degrades model performance and distorts best-fit values of model free parameters.     

基于水声环境空间中多模态深度融合模型的目标识别方法研究

随着对水下目标特性研究的深入和声学探测技术的发展,基于单模态的阵列式信息融合或基于空间信息的分布式信息融合的水下目标识别方法研究已有一定成果,但针对复杂海况导致单一物理场或单一融合层次的系统识别性能提高有限等方面影响的水下目标识别方法研究还有所不足,因此,开展基于多模态深度融合模型的水下目标识别方法研究可利用模态互补,共享信息而提升识别率。文中在国内外研究基础上,深入研究了基于到达时差法和多模态方法组合的检测方法,初步形成了基于水声环境空间中多模态深度融合模型的识别框架,开展了海洋中典型自然与人为事件的信号分析与特征提取,并在此基础上,设计新型基于海底基站的被动识别系统。该系统同步记录和由位置等组成的时间序列标记声、磁和压数据,可实现高精度、高分辨率的识别。本研究可满足未来海洋观测对高性能水下目标探测、定位和跟踪系统的迫切需要,为海洋安全监管、海洋突发事件应急响应等领域提供新的技术手段和科学参考。     

Seasonal cycle of circulation in the Antarctic Peninsula and the off-shelf transport of shelf waters into southern Drake Passage and Scotia Sea

The seasonal cycle of circulation and transport in the Antarctic Peninsula shelf region is investigated using a high-resolution (∼2 km) regional model based on the Regional Oceanic Modeling System (ROMS). The model also includes a naturally occurring tracer with a strong source over the shelf (radium isotope ²²⁸Ra, t_1/2=5.8 years) to investigate the sediment Fe input and its transport. The model is spun-up for three years using climatological boundary and surface forcing and then run for the 2004–2006 period using realistic forcing. Model results suggest a persistent and coherent circulation system throughout the year consisting of several major components that converge water masses from various sources toward Elephant Island. These currents are largely in geostrophic balance, driven by surface winds, topographic steering, and large-scale forcing. Strong off-shelf transport of the Fe-rich shelf waters takes place over the northeastern shelf/slope of Elephant Island, driven by a combination of topographic steering, extension of shelf currents, and strong horizontal mixing between the ACC and shelf waters. These results are generally consistent with recent and historical observational studies. Both the shelf circulation and off-shelf transport show a significant seasonality, mainly due to the seasonal changes of surface winds and large-scale circulation. Modeled and observed distributions of ²²⁸Ra suggest that a majority of Fe-rich upper layer waters exported off-shelf around Elephant Island are carried by the shelfbreak current and the Bransfield Strait Current from the shallow sills between Gerlache Strait and Livingston Island, and northern shelf of the South Shetland Islands, where strong winter mixing supplies much of the sediment derived nutrients (including Fe) input to the surface layer.     

Numerical modeling of gas hydrate emplacements in oceanic sediments

We have implemented a 2-dimensional numerical model for simulating gas hydrate and free gas accumulation in marine sediments. The starting equations are those of the conservation of the transport of momentum, energy, and mass, as well as those of the thermodynamics of methane hydrate stability and methane solubility in the pore-fluid. These constitutive equations are then integrated into a finite element in space, finite-difference in time scheme. We are then able to examine the formation and distribution of methane hydrate and free gas in a simple geologic framework, with respect to the geothermal heat flow, fluid flow, the methane in-situ production and basal flux. Three simulations are performed, leading to the build up of hydrate emplacements largely linear through time. Models act primarily as free gas accumulators and are relatively inefficient with respect to hydrate emplacements: 26–33% of formed methane are converted to hydrate. Seepage of methane across the sea-floor is negligible for fluid flow below 2. 10⁻¹¹ kg/m²/s. At 5.625 10⁻¹¹ kg/m²/s however, 9.7% of the formed methane seeps out of the model. Moreover, along strike variation arising in the 2-dimensional model are outlined. In the absence of focused flow, the thermodynamics of hydrate accumulation are primarily one-dimensional. However, changes in free methane compressibility (density) and methane solubility (the intrinsic dissolved methane flux) subtlety impact on the formation of a free gas zone and the distribution of the hydrate emplacements in our 2-dimensional simulations.     

Numerical modeling of gas hydrate emplacements in oceanic sediments

We have implemented a 2-dimensional numerical model for simulating gas hydrate and free gas accumulation in marine sediments. The starting equations are those of the conservation of the transport of momentum, energy, and mass, as well as those of the thermodynamics of methane hydrate stability and methane solubility in the pore-fluid. These constitutive equations are then integrated into a finite element in space, finite-difference in time scheme. We are then able to examine the formation and distribution of methane hydrate and free gas in a simple geologic framework, with respect to the geothermal heat flow, fluid flow, the methane in-situ production and basal flux. Three simulations are performed, leading to the build up of hydrate emplacements largely linear through time. Models act primarily as free gas accumulators and are relatively inefficient with respect to hydrate emplacements: 26–33% of formed methane are converted to hydrate. Seepage of methane across the sea-floor is negligible for fluid flow below 2. 10⁻¹¹ kg/m²/s. At 5.625 10⁻¹¹ kg/m²/s however, 9.7% of the formed methane seeps out of the model. Moreover, along strike variation arising in the 2-dimensional model are outlined. In the absence of focused flow, the thermodynamics of hydrate accumulation are primarily one-dimensional. However, changes in free methane compressibility (density) and methane solubility (the intrinsic dissolved methane flux) subtlety impact on the formation of a free gas zone and the distribution of the hydrate emplacements in our 2-dimensional simulations.