On a near-optimal control approach for a wave energy converter in irregular waves

Real-time smooth reactive control and optimal damping of wave energy converters in irregular waves is difficult in part because the radiation impulse response function is real and causal, which constrains the frequency-dependent added mass and radiation damping according to the Kramers–Kronig relations. Optimal control for maximum energy conversion requires independent synthesis of the impulse response functions corresponding to these two quantities. Since both are non-causal (one being odd and other even), full cancellation of reactive forces and matching of radiation damping requires knowledge or estimation of device velocity into the future. To address this difficulty and the non-causality of the exciting force impulse response function, this paper investigates the use of propagating-wave surface elevation up-wave of the device to synthesize the necessary forces. Long-crested waves are assumed, and the approach is based on the formulations of Naito and Nakamura [2] and Falnes [22]. A predominantly heaving submerged device comprised of three vertically stacked discs driving a linear power take-off is studied. The overall formulation leads to smooth control that is near-optimal, given the approximations involved in the time-shifting of the non-causal impulse response functions and the consequent up-wave distances at which wave surface elevation is required. Absorbed power performance with the near-optimal approach is compared with two other cases, (i) when single-frequency tuning is used based on non-real time adjustment of the reactive and resistive loads to maximize conversion at the spectral peak frequency, and (ii) when no control is applied with damping set to a constant value. Simulation results for wave spectra over a range of energy periods and significant wave heights are compared for the three situations studied. While practical implementation presents engineering challenges, in terms of time-averaged absorbed power, unconstrained near-optimal control is found to perform significantly better than single-frequency tuning in the spectra with longer energy periods (>10 s for the present device), and somewhat better in the spectra with shorter energy periods (here ≤10 s).     

Third order wave force on axisymmetric bodies

The third order triple-frequency wave load on fixed axisymmetric bodies by monochromatic waves is considered within the frame of potential theory. Waves are assumed to be weak non-linearity and a perturbation method is used to expand velocity potentials and wave loadings into series according to a wave steepness of kA. Integral equation method is used to compute velocity potentials up to second order in wave steepness. The third order triple-frequency wave loads are computed by an indirect method and an efficient method is applied to form the third order forcing term on the free surface quickly. The method can be used to compute third order triple-frequency surge force, heave force and pitch moment on any revolution bodies with vertical axes. The comparison with Malenica and Molin's results is made on surge force on a uniform cylinder, and comparison with experimental results is made on third order surge force, heave force and pitch moment on a truncated cylinder. More numerical computations are carried out for third order forces and moments on a uniform cylinder, truncated cylinders and a hemisphere.     

Constrained near-optimal control of a wave energy converter in three oscillation modes

This paper investigates the performance of a small axisymmetric buoy under wave-by-wave near optimal control in surge, heave, and pitch modes in long-crested irregular waves. Wave prediction is obtained using a deterministic propagation model. The paper describes the overall formulation leading up to the derivation of the feedforward control forces in surge and heave, and the control moment in pitch. The radiation coupling between surge and pitch modes is accounted for in the model. Actuation is relative to deeply submerged reaction masses. Heave oscillations are constrained by the swept-volume limit. Oscillation constraints are also applied on the surge and pitch oscillations. The paper discusses time-domain simulations for an irregular wave input with and without the present control. Also discussed are results obtained over a range of irregular wave conditions derived for energy periods from 7 s to 17 s, and a significant wave height of 1 m. It is found that, while the gains in power capture enabled by the present control are significant, the actuation forces are also very large, given the small size of the buoy. Further, due to the small size, heave is found to be the dominant contributor to power capture, with relatively modest contributions from surge and pitch.     

On non-causal impulse response functions related to propagating water waves

Particular attention is drawn to the non-causality of two impulse response functions, one of which relates the excitation force on an immersed body to the incident wave elevation at the body's reference position, while the other relates the incident wave elevations at two different positions along the line of wave propagation. An explanation is proposed for the non-causality of the impulse response functions, in spite of the fact wave propagation is a causal process. An indication is given of how far ‘upstream’ the incident wave elevation should be measured in order to be able to know, with reasonable accuracy, the current heave excitation force on a floating truncated cylinder with vertical axis, given current and past wave-elevation measurements. This provides a method for wave prediction, which is required for optimum control of the oscillation of the immersed body.     

Deep ocean wave energy conversion using a cycloidal turbine

A lift based wave energy converter, namely, a cycloidal turbine, is investigated. 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. In the two-dimensional limit, i.e. for large spans of the hydrofoil (or an array of these), the geometry of the converter is suitable for wave termination of straight crested Airy waves. Results from two-dimensional potential flow simulations, with thin hydrofoils modeled as either a point vortex or discrete vortex panel, are presented. The operation of the cycloidal turbine both as a wave generator as well as a wave-to-shaft energy converter interacting with a linear Airy wave is demonstrated. The impact on the performance of the converter for design parameters such as device size, submergence depth, and number of hydrofoils is shown. For optimal parameter choices, simulation results demonstrate inviscid energy conversion efficiencies of more than 99% of the incoming wave energy to shaft energy. This is achieved using feedback control to synchronize the rotational rate, blade pitch angle, and phase of the cycloidal wave energy converter to the incoming wave. While complete termination of the incoming wave is shown, the remainder of the energy is lost to harmonic waves traveling in the up-wave and down-wave directions.     

Hydrodynamic Analysis and Power Conversion for Point Absorber WEC with Two Degrees of Freedom Using CFD

Point absorber wave energy device with multiple degrees of freedom (DOF) is assumed to have a better absorption ability of mechanical energy from ocean waves. In this paper, a coaxial symmetric articulated point absorber wave energy converter with two degrees of freedom is presented. The mechanical equations of the oscillation buoy with power take-off mechanism (PTO) in regular waves are established. The three-dimensional numerical wave tank is built in consideration of the buoy motion based upon the CFD method. The appropriate simulation elements are selected for the buoy and wave parameters. The feasibility of the CFD method is verified through the contrast between the numerical simulation results of typical wave conditions and test results. In such case, the buoy with single DOF of heave, pitch and their coupling motion considering free (no PTO damping) and damped oscillations in regular waves are simulated by using the verified CFD method respectively. The hydrodynamic and wave energy conversion characteristics with typical wave conditions are analyzed. The numerical results show that the heave and pitch can affect each other in the buoy coupling motion, hydrodynamic loads, wave energy absorption and flow field. The total capture width ratio with two coupled DOF motion is higher than that with a single DOF motion. The wave energy conversion of a certain DOF motion may be higher than that of the single certain DOF motion even though the wave is at the resonance period. When the wave periods are high enough, the interaction between the coupled DOF motions can be neglected.     

Influence of wave approach angle on TLP's response

The current study focuses on the response analysis of triangular tension leg platform (TLP) for different wave approach angles varying from 0° through 90° and its influence on the coupled dynamic response of triangular TLPs. Hydrodynamic loading is modeled using Stokes fifth-order nonlinear wave theory along with various other nonlinearities arising caused by change in tether tension and change in buoyancy caused by set down effect. Low frequency surge oscillations and high frequency tension oscillations of tethers are ignored in the analysis. Results show that wave approach angle influences the coupled dynamic response of triangular TLP in all degrees of freedom except heave. Response in roll and sway degrees of freedom are activated which otherwise are not present in TLP's response to unidirectional waves. Pitch and roll responses are highly stochastic in nature indicating high degree of randomness. Variation in surge, sway and heave responses are nonlinear and are not proportional to change in wave height for the same period.     

Estimating the potential of ocean wave power resources

The realistic assessment of an ocean wave energy resource that can be converted to an electrical power at various offshore sites depends upon many factors, and these include estimating the resource recognizing the random nature of the site-specific wave field, and optimizing the power conversion from particular wave energy conversion devices. In order to better account for the uncertainty in wave power resource estimates, conditional probability distribution functions of wave power in a given sea-state are derived. Theoretical expressions for the deep and shallow water limits are derived and the role of spectral width is studied. The theoretical model estimates were compared with the statistics obtained from the wave-by-wave analysis of JONSWAP based ocean wave time-series. It was shown that the narrow-band approximation is appropriate when the variability due to wave period is negligible. The application of the short-term models in evaluating the long-term wave power resource at a site was illustrated using wave data measured off the California coast. The final example illustrates the procedure for incorporating the local wave data and the sea-state model together with a wave energy device to obtain an estimate of the potential wave energy that could be converted into a usable energy resource.     

Implementation and evaluation of alternative wave breaking formulas in a coastal spectral wave model

This paper describes methods and results of research for incorporating four different parameterized wave breaking and dissipation formulas in a coastal wave prediction model. Two formulations assume the breaking energy dissipation to be limited by the Rayleigh distribution, whereas the other two represent the breaking wave energy by a bore model. These four formulations have been implemented in WABED, a directional spectral wave model based on the wave action balance equation with diffraction, reflection, and wave–current interaction capabilities. Four parameterized wave breaking formulations are evaluated in the present study using two high-quality laboratory data sets. The first data set is from a wave transformation experiment at an idealized inlet entrance, representing four incident irregular waves in a slack tide and two steady-state ebb current conditions. The second data set is from a laboratory study of wave propagation over a complex bathymetry with strong wave-induced currents. Numerical simulation results show that with a proper breaking formulation the wave model can reproduce laboratory data for waves propagating over idealized or complicated bathymetries with ambient currents. The extended Goda wave breaking formulation with a truncated Rayleigh distribution, and the Battjes and Janssen formulation with a bore model produced the best agreement between model and data.     

Processing irregular wave measurements to enhance point absorber power capture performance

The paper aims at introducing practical methods for power capture performance enhancement of a heaving wave energy converter in irregular seas. The optimum control solution requires tuning to wave frequency based on wave force information. However, identification of the wave frequency in irregular seas is considered to be a complex and difficult task. This is partly due to technical difficulties in determination of the wave force. Besides, there are no clear guidelines for identification of wave frequency from an irregular sea state based wave force information. In a typical application, one of the available sources of information about the wave properties is the wave elevation record. The proposed approach presents a method for estimation of the wave frequency information from the wave elevation data by using signal processing and filtering techniques. The proposed method uses filters to generate an estimation of wave force information, which is used to identify the local wave frequency by method of a time-series analysis of the data. This wave frequency information is then used in tuning the device. The details of the proposed techniques, the model of the wave energy converter, the simulated sea states and the related simulation results are also presented.     

Laboratory study of wave and turbulence characteristics in narrow-band irregular breaking waves

Wave elevations and water particle velocities were measured in a laboratory surf zone created by the breaking of a narrow-band irregular wave train on a 1/35 plane slope. The incident waves form wave groups that are strongly modulated. It is found that the waves that break close to the shoreline generally have larger wave-height-to-water-depth ratios before breaking than the waves that break farther offshore. After breaking, the wave-height-to-water-depth ratio for the individual waves approaches a constant value in the inner surf zone, while the standard deviation of the wave period increases as the still water depth decreases. In the outer surf zone, the distribution of the period-averaged turbulent kinetic energy is closely correlated to the initial wave heights, and has a wider variation for narrow-band waves than for broad-band waves. In the inner surf zone, the distribution of the period-averaged turbulent kinetic energy is similar for narrow-band waves and broad-band waves. It is found that the wave elevation and turbulent kinetic energy time histories for the individual waves in a wave group are qualitatively similar to those found in a spilling regular wave. The time-averaged transport of turbulent kinetic energy by the ensemble-averaged velocity and turbulence velocity under the irregular breaking waves are also consistent with the measurements obtained in regular breaking waves. The experimental results indicate that the shape of the incident wave spectrum has a significant effect on the temporal and spatial variability of wave breaking and the distribution of turbulent kinetic energy in the outer surf zone. In the inner surf zone, however, the distribution of turbulent kinetic energy is relatively insensitive to the shape of the incident wave spectrum, and the important parameters are the significant wave height and period of the incident waves, and the beach slope.     

Stochastic control of wave energy converters with constrained displacements for optimal power absorption

An semi-analytical solution is derived for the optimal control of the power take-off of a single-degree of freedom heave point absorber with constraints on the displacement. At first the control force is derived during states, where the displacement constraint is active. This results in an open-loop control law dependent on the external wave load on the absorber. Next, the analytical solution for the optimal control in the unconstrained state is indicated, which turns out to be of the closed loop type with feedback from the present displacement and acceleration and from future velocities. The derived control law contains an undetermined constant, which is calibrated at the interface to the previous constrained state. The approach requires the estimation of the wave load during the constrained states, and the prediction of the future velocity response during unconstrained states. An algorithm has been devised in the paper for handling these problems. The theory has been validated against numerical solutions obtained by nonlinear programming.     

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.     

Wave energy conversion under constrained wave-by-wave impedance matching with amplitude and phase-match limits

This paper investigates an approach to limit the fullness of ‘tuning’ provided by wave-by-wave impedance matching control of wave energy devices in irregular waves. A single analytical formulation based on the Lagrange multiplier approach of Evans [1] is used to limit the velocity amplitude while also limiting the closeness of the phase match between velocity and exciting force. The paper studies the effect of the present technique in concurrently limiting the device velocity and the required control/actuation force. Time domain application requires wave-profile prediction, which here is based on a deterministic propagation model. Also examined in the time domain is the effect of possible violation of the displacement constraint, which for many designs implies impacts at hard stops within the power take-off mechanism. Time domain simulations are carried out for a 2-body axisymmetric converter (with physical end-stops) in sea states reported for a site off the US east coast. It is found that the approach leads to effective power conversion in the less energetic sea states, while as desired, considerable muting of the optimal response is found in the larger sea states. Under the assumptions of this work, the end-stop collisions are found to have a minor effect on the power conversion. The present approach could be used to guide the design of power take-off systems so that their displacement stroke, maximum force, and resistive and reactive power limits are well-matched to the achievable performance of a given controlled primary energy converter.     

Interaction of oblique waves with infinite number of perforated caissons

An analytic solution based on the division of the fluid domain is developed for the interaction of obliquely incident waves with infinite number of perforated caissons. The whole fluid domain is firstly divided into infinite sub-domains according to the division of structures, and subsequently eigenfunction expansion is employed to represent the velocity potential in each domain. A phase relation is utilized for the analysis of wave oscillation in each caisson, and the character of structure geometry is considered in setting up the mathematical model of reflection waves. The reflection waves from the present analysis include many propagation waves traveling in different directions when the incident wave frequency is high. Benchmark examinations show that the continuous condition of water particle velocity is satisfied at the front walls of caissons, and the reflection coefficients keep agreement with the energy conservation relation very well when porous effect parameter is infinite. Numerical results show that the reflection coefficients of obliquely incident waves are smaller when the length of caissons is shorter at low frequency. The wave reflection coefficients and the wave forces normal to caissons decrease and the wave forces along caissons increase with the increase of the wave incident angle.     

Accurate determination of internal kinematics from numerical wave model results

The internal kinematics for surface waves propagating over a locally constant depth are expressed as convolution integrals. Given the wave kinematics at the still water level (SWL), this provides explicit and exact potential flow expressions for the internal kinematics as convolutions in space with appropriate impulse response functions. These functions are derived in closed form and they are shown to decay exponentially. This effectively reduces the limits of the convolution integral to a horizontal distance of approximately three water depths from the water column of interest. The SWL kinematics must be provided within this region. The source of SWL kinematics may, e.g. be one of the recently developed highly accurate Boussinesq-type formulations. The method is valid for multidirectional, irregular waves of arbitrary nonlinearity at any constant water depth.     

Linear irregular wave generation in a numerical wave tank

In the design of any floating or fixed marine structure, it is vital to test models in order to understand the fluid/structure interaction involved. A relatively inexpensive method, compared to physical model testing, of achieving this is to numerically model the structure and the wave conditions in a numerical wave tank. In this paper, a methodology for accurately replicating measured ocean waves in a numerical model at full scale is detailed. A Fourier analysis of the measured record allows the wave to be defined as a summation of linear waves and, therefore, Airy's linear wave theory may be used to input the wave elevation and associated water particle velocities. Furthermore, a structure is introduced into the model to display the ability of the model to accurately predict wave–structure interaction. A case study of three individual measured waves, which are recorded at the Atlantic marine energy test site, off the west coast of Ireland, is also presented. The accuracy of the model to replicate the measured waves and perform wave–structure interaction is found to be very high. Additionally, the absolute water particle velocity profile below the wave from the numerical model is compared to a filtered analytical approximation of the measured wave at a number of time-steps and is in very good agreement.     

船载波浪发电设备设计及其性能研究

为解决小型无人船由于船体空间紧凑,携带燃料较少导致的续航力差、航程短等问题,本文以三体船为平台,将摆式波浪能发电技术与多体船相结合,提出了一种新型船载波浪发电设备设计方案。该发电设备可由主侧船体之间的波浪运动响应差异驱动,在海上航行时捕获波浪能并转化为电能,以此来补充船上能源。通过数值软件对三体船进行性能预报,结果表明,在四级海况下迎浪航行时,三体船平均波浪能捕获功率可达3.57 kW,波浪能一级转换效率为7.32%,能够有效的补充船上能源。同时,安装船载波浪发电设备具备减摇减荡的效果,提升了三体船的航行稳定性,进一步体现了该装置的可行性和实用性。     

Wave prediction using wave rider position measurements and NARX network in wave energy conversion

Several control methods of wave energy converters (WECs) need prediction in the future of wave surface elevation. Prediction of wave surface elevation can be performed using measurements of surface elevation at a location ahead of the controlled WEC in the upcoming wave. Artificial neural network (ANN) is a robust data-learning tool, and is proposed in this study to predict the surface elevation at the WEC location using measurements of wave elevation at ahead located sensor (a wave rider buoy). The nonlinear autoregressive with exogenous input network (NARX NN) is utilized in this study as the prediction method. Simulations show promising results for predicting the wave surface elevation. Challenges of using real measurements data are also discussed in this paper.     

Triangular Configuration Tension Leg Platform behaviour under random sea wave loads

This study investigates the dynamic response of a Triangular Configuration Tension Leg Platform (TLP) under random sea wave loads. The random wave has been generated synthetically using the Monte-Carlo simulation with the Peirson–Moskowitz (P–M) spectrum. Diffraction effects and second-order wave forces have not been considered. The evaluation of hydrodynamic forces is carried out using the modified Morison equation with water particle kinematics evaluated using Airy's linear wave theory. Wave forces are taken to be acting in the surge degree-of-freedom. The effect of coupling of various structural degrees-of-freedom (surge, sway, heave, roll, pitch and yaw) on the dynamic response of the TLP under random wave loads is studied. Parametric studies for random waves with different H_s and T_z under the presence of current have also been carried out. For the orientation of the TLP, surge, heave and pitch degrees-of-freedom responses are influenced significantly. The surge power spectral density function (PSDF) indicates that the mean square response is affected by the amplification at the natural frequency of the surge degree-of-freedom and also at the peak frequency of the wave loading. The PSDF of the heave response shows higher peak values near the surge frequency and near the peak frequency of the wave loading. Surge response, therefore, influences heave response to the maximum. Variable submergence seems to be a major source of nonlinearity and significantly enhances the responses in surge, heave and pitch degrees-of-freedom. In the presence of current, the response behaviour of the TLP is altered significantly introducing a non-zero mean response in all degrees-of-freedom.