Discrete Displacement Hydraulic Power Take-off System for the Wavestar Wave Energy Converter

The Wavestar Wave Energy Converter (WEC) is a multiple absorber concept, consisting of 20 hemisphere shaped floats attached to a single platform. The heart of the Wavestar WEC is the Power TakeOff (PTO) system, converting the wave induced motion of the floats into a steady power output to the grid. In the present work, a PTO based on a novel discrete displacement fluid power technology is explored for the Wavestar WEC. Absorption of power from the floats is performed by hydraulic cylinders, supplying power to a common fixed pressure system with accumulators for energy smoothing. The stored pressure energy is converted into electricity at a steady pace by hydraulic motors and generators. The storage, thereby, decouples the complicated process of wave power absorption from power generation. The core for enabling this PTO technology is implementing a near loss-free force control of the energy absorbing cylinders. This is achieved by using special multi-chambered cylinders, where the different chambers may be connected to the available system pressures using fast on/off valves. Resultantly, a Discrete Displacement Cylinder (DDC) is created, allowing near loss free discrete force control. This paper presents a complete PTO system for a 20 float Wavestar based on the DDC. The WEC and PTO is rigorously modeled from incident waves to the electric output to the grid. The resulting model of +600 states is simulated in different irregular seas, showing that power conversion efficiencies above 70% from input power to electrical power is achievable for all relevant sea conditions. 4002 Nomenclature: β eff Effective bulk modulus of a volume of fluid [Pa] ∆p Pressure difference across motor ports [Pa] η w Wave height [m] η c Efficiency of cylinder [-] η PTO PTO conversion efficiency from input to output and vice versa [-] η tot Ratio of input and output energy of the PTO [-] η DCC Average conversion efficiency of DCC [-] γ Control input to reduce power absorption [-] ν Kinematic viscosity of fluid [m 2 /s] θ arm Angular position of arm [rad] θ w Approach angle of incoming waves [rad] ρ Density of hydraulic fluid [kg/m 3 ] σ Variance [-] τ a Thermal time constant of accumulator [s] τ aH Thermal time constant of accumulator [s] τ Arch Torque on arm due to Archimedes force [Nm] τ G Torque on arm due to gravity on float and arm [Nm] τ ext Torque on arm due …


Introduction
A large group of Wave Energy Converter (WEC) concepts are based on extracting energy from ocean waves using the principle of placing buoyant bodies in the sea.As waves pass the bodies, these are forced to oscillate.This is illustrated in Figure 1, showing different version of such WECs; see [1,2] for a survey.The body movements are converted into electricity by a system called the Power Take-Off (PTO); see Figure 2a or Figure 3a.
The development of an efficient and reliable PTO system is a main challenge for WECs.The problem faced is that the movement or oscillations of bodies caused by sea waves are very slow, bidirectional and irregular.An example is given in Figure 2a, showing the resulting velocities of a 5 m diameter point absorber from the Wavestar concept.The torque, τ PTO , is the load torque applied by the PTO in order to extract energy from the absorber's pivoting motion.The graphs in Figure 2b shows the slow oscillation (≈5 s period), resulting in a requirement of a 1 MNm bi-directional load torque for extracting an average of 27 kW.In comparison, a wind turbine operating at approximately 13 RPM would produce 1.4 MW when loading the main rotor shaft with 1 MNm.Hence, the torque density requirement for wave power is immense.Furthermore, the incoming wave power is characterized by having more than a factor of 10 between mean and peak power, as seen in Figure 2c.The peak contributes heavily to the overall production [3] and may not simply be discarded.Figure 2c also shows that the energy is grouped.Resultantly, energy storage is required to store the energy peaks and use them for maintaining production between wave groups.To incorporate this, power absorption requires a PTO having a complex transmission with a very high gearing ratio in order to make ocean waves power a generator.Moreover, to extract adequate amounts of energy from waves, the PTO load force, F PTO , or torque, τ PTO , as in Figure 2b, applied by the PTO should be controlled as a function of incident wave and body movement.Finally, to fully maximize energy extraction, the PTO should be able to operate in four quadrant mode.Four quadrant mode means that all four combination of force and velocity direction should be provided, as the PTO sometimes is required to aid the float motion (reactive control strategies).An overview of required PTO characteristics are well described in [3].
The controllability and four quadrant behavior is required to compensate for the inherit off-resonance behavior of point absorbers.Frequency-wise, point absorbers are characterized by being narrow-banded with an under-damped resonance.To correct this deficiency, the force or torque, τ PTO , applied by the PTO to the absorber is controlled as a feedback of the absorber's motion, allowing the PTO to adjust the absorber's resonance frequency to match the wave frequency [4,5].Optimal energy transfer from wave to an absorber is obtained when the incident wave frequency matches the resonance frequency of the absorber.
Regarding implementation of PTO systems for wave energy, several investigations have been performed on using linear generators driven directly by the absorbers' movement.However, due to the slow linear velocities (peak linear velocity about 2 m/s), conventional permanent magnet linear generators would become very large.Typical achievable air-gap sheer-stress levels between stator and translator is about 20-25 kN/m 2 [6].Resultantly, a linear generator for a Wavestar float of 5 m diameter, requiring a load force of 400 kN, would be equivalent to a generator with 16-20 m 2 of active air-gap surface.In [7], the weight of active magnetic material alone (copper, iron laminations, magnets and back-iron) is estimated to be 1500 kg/m 2 .Neglecting the required support structure, the material requirements still is 24-30 tons, rendering the solution infeasible.Resultantly, effort is put into using a transmission combined with more conventional generators.Mechanical transmissions have been explored.However, these would be too massive, as a gearing ratio, which is a factor of 10 higher than a wind turbine transmission, is required, along with handling a bidirectional input.
Fluid power is a suitable technology for implementing the required transmission, as it is capable of producing high controllable forces at low velocities and easily "rectify" the bidirectional movement with compact actuators (cylinders).Unfortunately, fluid power systems are often characterized by poor efficiencies when operating at part load, which is crucial with the high ratio between peak and mean power in wave energy.
A conventional hydraulic transmission for wave energy is seen in Figure3b, where a cylinder operates as a pump, producing a bi-directional flow, which drives a hydraulic motor.The motor adapts to the flow and rectifies the flow into a unidirectional turning of the generator.However, an optimization of such a PTO from wave-to-grid is performed in [8] for the Wavestar converter, showing an overall PTO power conversion efficiency below 65% at the optimum point, quickly dropping to 45% in smaller waves.The system in Figure 3b also has the shortcoming of not allowing energy storage/smoothing.
PTO systems with the cylinder operating as a passive pump against a steady pressure, as in Figure 3c, have been used in [9,10], where energy smoothing may be performed using hydraulic accumulators.This allows operating the hydraulic motor and generator at a fairly constant load, yielding a PTO efficiency of up to 80%.However, the cylinder is limited to providing a Coulomb-like force load, reducing the amount of absorbed energy [11].
If the active valves are used instead of the passive check valves in Figure 3c, the absorber motion may be controlled though latching to improve energy capture [12].The latching control prolongs the natural period of the absorber motion non-linearly by locking the absorber's movement in parts of an oscillation cycle.This may be implemented by closing the valves, blocking the cylinder and absorber motion.This may yield the same or higher energy extraction as a load force controlled as a linear damper.In [13], it is suggested that the energy capture of Figure 3c may be increased by modifying with an extra accumulator with a controllable on/off valve.By, in turn, storing and releasing energy from the accumulator, the motion amplitude of the absorber may be improved.
In Figure 3c, a PTO is shown based on a hydraulic transformer [14], capable of controlling the force of the cylinder, while having a fixed system pressure with energy smoothing accumulators.However, the part load efficiency of the hydraulic transformer is poor, as it is basically two variable displacement pump/motors back-to-back.
It would be desirable to combine the positive features of the above system characteristics, i.e., having a PTO with a common fixed pressure system with accumulators for energy smoothing, combined with an efficient force control of the cylinders.This has been investigated for the Wavestar system in [15,16], where the force control of a hydraulic cylinder is based on connecting the chambers to different fixed system pressures using an arrangement of on/off vales, as shown in Figure 4a.
In Figure 4a, a cylinder with four chambers is attached to a manifold with eight on/off valves.By either having system pressure (red) or low pressure (blue) in the chambers, 16 different pressure and chamber combinations are achievable for this particular embodiment, yielding 16 different available PTO forces, as shown in Figure 4b.This corresponds to having 16 different gears available.Thus, during a wave, the PTO force is varied discretely, as shown in Figure 4c, by shifting "gear" every half second, thereby approximating a continuous force reference.The resulting PTO force, F PTO , is generated as the sum of forces produced by the different cylinder chambers: As the "gear" shifting is actually a discrete variation of the cylinder displacements using the valve system, the cylinder with manifold is referred to as a Discrete Displacement Cylinder (DDC).In [15], it was shown that despite assuming infinite fast and large switching valves, a certain amount of energy will always be lost when shifting force due to the compressibility of the fluid.To assess whether a PTO system may be feasible with this unavoidable compression loss, an estimate of the PTO efficiency was calculated for the Wavestar WEC in [15].This was performed for different wave conditions, showing that an efficiency above 90% was reachable for the DDC.In [16], it was shown that if the opening and closing time of the valves was kept less than 15 ms, the minimum loss could almost be assumed, thereby making the DDC feasible for a PTO solution.
Discrete control of a hydraulic cylinder has also been investigated for mobile hydraulics in [17], showing promising energy saving compared to, e.g., conventional load sensing systems.A PTO system with a similar approach, utilizing two asymmetric cylinders has previously been discussed in [18].The efficiency of controlling the force of the cylinders by pressure shifting was found to be between 88% and 94%, excluding the friction of the cylinder.Currently, [19] has shown a hydraulic transmission with accumulators for power smoothing operating at 70% efficiency.A similar system as in [18] is suggested in [20].The system of [20] is tested in a scaled version in [21], capable of applying a torque of 16 kNm.The PTO is tested in a test-rig, simulating regular waves.The overall efficiency is estimated to be from 69% to 80%.Estimates were given, as the test-rig was not fully operational.The PTO system presented in this paper is based on the results in [15,16] and is a complete PTO for the 20 float Wavestar C-concept shown in Figure 5b with 5 m diameter floats.The C-concept consists of semi-submerged hemisphere-shaped floats mounted on separate arms.Each arm is mounted and rotationally supported by a common platform, which is resting on the sea bed through a number of piles.One advantage of a multiple point absorber is the increased power smoothing achieved when the waves passes through the systems.The PTO components are mainly enclosed within the main structure; only the cylinder is mounted externally.The PTO is based on a hydraulic transmission, using a fluid that is biodegradable in the marine environment.For all components inside the structure, an extra level of protection is provided, such that in the rare case of leakage, the fluid is contained inside the WEC.For the cylinder outside the machine, a sump is mounted below the cylinder, connected to automated suction lines from inside the WEC.In case of registering oil in the sump, the cylinder and float is lifted into storm protection and taken out of service until maintenance.
For storm protection and extreme seas, the Wavestar concept incorporates a jacking-system, allowing the floats and platform to be lifted out of the water.Resultantly, the PTO system and structure only have to be designed for "production" waves (which, for the C5, is 3 m significant wave height).Wavestar has in 2009 installed a two-float prototype of the 600 kW C-concept at Hanstholm at the West coast of Denmark, which has been in operation ever since [22].The C-concept is a part of the strategy for integration of wave power into wind farms, as shown in Figure 5c, with the SC-concept (star concept).This reduces both establishment, installation and maintenance cost, while having an energy park with a higher power density and a more stable energy output.
A key property of the presented PTO concept is being scalable to future multi-megawatt Wavestar systems in Figure 5c.

Methods
The evaluation and presentation of the PTO is performed by first describing the layout and main components of the PTO.This description also presents the main features of PTO and explains the overall operation.
After the PTO description, the modeling of the system is attended.The modeling covers wave models, hydrodynamic-model, hydraulic system and electric system.It has not been possible to use standard component simulation software, as the system is very multi-disciplinary.Additionally, complete insight of the components models is required to avoid simulating unnecessary dynamics.Otherwise, poorly conditioned stiff systems may be obtained with slow execution times.Resultantly, modeling of all systems to the required detail level has been performed and implemented in SIMULINK R .Using this approach, the model still ends up having more than 600 states, but is able to simulate about a factor of 4 slower than real time on a reasonable work station.The reasonable execution time is important, as the model presented in the paper is used for optimizing the design.
After the modeling section, the basic control of the system is presented.This consists of two parts: The Wave Power Extraction Algorithm (WPEA), which covers how a float should be loaded to extract wave power, maximizing the PTO energy production.The presented algorithm is based on reactive control methods.The second part of the control section sketches how the PTO is controlled to implement the chosen WPEA algorithm, while ensuring high efficiency and a stable power output to the grid.
Finally, simulation results for the 20 floats system in operation are given for different sea states, evaluating the performance of the PTO system.The simulations presented are of relative short duration (approximately 100) wave periods.However, based on a number of generated short waves, a selection of short wave series have been tested to give approximately the same average power extraction and PTO performance results compared to long simulation of approximately 1000 wave periods.Resultantly, these short wave series are applicable for comparing and validating design iterations.
For complete statistical background of performance evaluation, final simulations with more subsystems and more complex wave models are required.This is, however, not a part of the scope of this paper.
For the reader not interested in the detailed modeling of the PTO, the paper is constructed such that after the PTO layout section, the reader may move directly to the control section and results.

PTO Layout
The PTO concept is illustrated in Figure 6 for a 20 float Wavestar with 5 m diameter floats.As seen in Figure 6a, each float arm is equipped with a Discrete Displacement Cylinder (DDC)-system, consisting of a multi-chambered cylinder with integrated shifting manifold, controlling the resulting force, F PTO .These DDCs provide independent force control of the floats.All DDCs supply oil flow into a common pressure line system, consisting of three pressure lines: • Low pressure line, p L : 10 bar-30 bar; • Mid-(intermediate) pressure line, p M : 70 bar-170 bar; • High pressure line, p H : 150 bar-320 bar.
These pressure lines run through the entire main tube, thereby collecting all the extracted power into common hydraulic energy storages.Hydraulic motors and generators convert the stored pressure energy to electricity.In this way, the complicated wave power absorption process and the electricity generation becomes decoupled.

Hydraulic Motors and Generators
Four main pairs of hydraulic motors and generators are placed on the high pressure line for power generation.These generators are strategically distributed across the Wavestar WEC (at float row 2, 5 and 8) to minimize pipe losses.Two parallel generators are placed at float row 5 for extra generator capacity in the center of the system during high energetic sea states.Each motor and generator pair consists of a 160 kW asynchronous generator driven by a 250 cm 3 /rev hydraulic motor, operating the generator at a torque of 1193 Nm at a pressure of 300 bar.The hydraulic motors are fixed displacement bent-axis motors, whose efficiency is about 95% in these operating conditions.
At each generator point, a storage battery of gas-loaded hydraulic accumulators is placed.These accumulators perform short-term energy storage for smoothing the extracted wave power.The size of the storage may be increased as desired to meet the required smoothing level.
As the pressure in the lines is only slowly varying, due to the accumulators, the torque produced by the motors is near constant, thereby driving the generators at a stable load.To adapt to the overall variation in power input from the waves during a sea state, the speed of the generators is controlled using a converter.Each set of generator has its own converter, allowing speed control of the generators, while feeding 690 V at 50 Hz to the grid.During less energetic seas, different combinations of generators and converters may be closed down, such that only the required number of generators are operating.

Pressure Line System and Accumulators
Power generation is only performed from the high pressure line, also having the largest storage capacity.The hydraulic motors driving the generators extract power from the high pressure line and delivers flow back to the low pressure line in a closed-circuit arrangement.The low pressure line is pressurized to avoid cavitation in the cylinders during pressure shifting in the chambers.System may include as many accumulators for power smoothening as required.
27 discrete forces available in DDC for control.
Generator speed controlled using inverters to adapt to overall available wave power Generators and motors may be disconneted during low energy sea conditions 4x250cc h All components are within machine main tube Force distribution for DDC-system Charging sys.

Flushing valve
To filter the oil in the closed-circuit system and replenish flow leakage in the motors, a charging pump is installed, supplying flow into the low pressure line.
The mid-pressure line in the system is an intermediate pressure, which is used for increasing the efficiency of the DDCs.The compression loss associated with a pressure shift in a cylinder chamber is proportional to the square of the pressure change.Resultantly, having the option of shifting to an intermediate pressure reduces this compression loss.The extra pressure line also increases the force resolution of the DDC-system from 2 3 = 16 to 3 3 = 27 values.
For energy storage and stabilizing the pressure lines, the batteries of gas-loaded hydraulic piston accumulators are placed on each line.When used for energy storage, the round-trip efficiency of accumulators are ≥94% [23].To increase utilization of the accumulators' storage capacity, the pressure lines are allowed to vary with a minimum factor of two in pressure level.At these conditions, a 50 L piston accumulator is able to store about 23 L oil.Thus, a battery of 16 may store 375 L of high pressure oil, which may keep a 250 cm 3 /rev motor running at 1500 RPM for one minute.The low pressure accumulator batteries each consists of 10 × 50 L accumulators, the mid-pressure batteries of 4 × 50 L accumulators and the high pressure of 16 × 50 L.

DDC
The core technology for enabling the PTO concept with common pressure lines is the DDC, implementing the throttle-less force control of the cylinder.The manifold of the DDC in Figure 6a is shown in Figure 6b.The cylinder has three active chambers, where each chamber may be connected to either low, mid-or high pressure.As the piston areas of the chambers varies in size, 3 3 = 27 different resulting forces, F PTO , may be generated, as illustrated in Figure 6d.
The stroke of the cylinder is 3 m.To perform the shifting of chamber pressures, a shifting manifold, as in Figure 6c, with nine fast on/off valves is used (less than 15 ms opening and closing time).The pressure drop across the valves is low, i.e., less than 3 bar, at flows corresponding to a cylinder velocity of 0.5 m/s.
The manifold is directly mounted on the cylinder to reduce volume (compression losses) and to reduce transmission phenomena in the pipe connections.The DDC-system shifts force approximately every 400 ms.To supply/consume the short flow pulses required for compression or decompression of the volumes, small accumulators are mounted directly on the DDC-manifold.
An important property of the DDC is to be able to process extreme power levels with low losses, as the cylinder chambers are directly connected to the pressure lines when in a given force configuration.Only the compression losses are constant, which, thus, become more significant at lower power levels.However, efficiency is maintained by reducing the overall pressure levels, as the high DDC forces are not required in these conditions.
Another important feature of the PTO concept with the DDC is to always be able to reduce the load force independent of the float velocity and position.In this way, the WEC may always reduce the power absorption of the floats if full load capacity is reached, similar to wind turbine pitching out of the wind.
Flow to the intermediate pressure line is only supplied and extracted by the cylinders to reduce cost.Thus, the control of the floating intermediate pressure is performed collectively by all DDCs, making sure that the net-flow into the mid-pressure system is zero over time.This means that the DDC-modules actually operate as a sort of loss-less pressure amplifier.For example, the DDC first lets the float motion supply flow into the mid-pressure line.An instant later, the DDC shifts configuration, such that the mid-pressure and float motion are now combined to produce a high pressure flow.An example of such a cycle is illustrated in Figure7a: The cylinder is supplying fluid into the mid-pressure at time instance (1); At instance (2) the DDC increases the load force, F PTO , using a configuration where the second chamber is now pumping fluid into the high pressure line, while the first chamber now consumes flow from the mid-pressure.In instance (3), the force is further increased, such that the cylinder is pumping flow into high pressure without aid from the mid-pressure.
The different configurations of pressure and chambers and the resulting cylinder forces are seen in Figure 7b.
Available cylinder force F PTO [kN] (a) (b) Considering energy extraction, the DCC approximates a continuous reference with discrete values.This may affect the amount of extracted energy.However, comparing [15] (discrete control) and [8] (continuous control), the power extraction appears unaffected.The reason for the low impact of the discrete implementation is that the 35 tons of float and arm serves as an effective low pass filter, having a break frequency around 0.28 Hz.Thus, from a load control point of view, implementation of extra low pass filters, as performed in, e.g., digital power electronics, is not required.
This concludes the overall PTO description, and the PTO modeling section follows, which also further specifies and describes the PTO.To skip the modeling section and move directly to control of the PTO, see Section 7.

Modeling
An overview of the modeled sub-systems is seen in Figure 8.These are treated in the shown order.The derived equations are implemented in Matlab/ R Simulink R and solved with a Runge-Kutta solver, running at a fixed step time of 0.5 ms.

M
Wave and Floats Absorbtion system Generation system

Wave and Float Model
An irregular wave is often characterized by two quantities, the significant wave height, H m,0 , and the peak wave period, T p .The significant wave height, H m,0 , is the average wave height of the one-third highest waves, and T p is the wave period, where most energy is concentrated.These two parameters along with an underlying assumption of the shape of the power density spectrum, S(f ), defines a sea state.In this study, the Pierson-Moskowitz (PM) spectrum is utilized [24], describing a fully developed sea state: If the PTO may control the phase of the absorber, the PM spectrum may also be the conservative choice in regards to power extraction.The PTO often tunes the absorber resonance to match the peak period; however, the power in the PM spectrum is spread across a wide band of frequencies.In comparison, an often used spectrum, the JONSWAP (Joint North Sea Wave Project) spectrum, has the energy concentrated in a more narrow frequency band (peak enhancement factor larger than one).The PM spectrum in Equation ( 2) is shown in Figure 9a.The figure also shows the three sea states, which are going to be used in the simulation.To perform time simulation, a time realization of a wave is required, complying with the spectrum.Often, the random phase method is used, where the spectrum is converted into a finite number of regular wave components, which are added with a random phase for each component.This simple method has been experienced to be adequate for estimating average power absorption [25].However, the method does not reproduce, e.g., wave grouping correctly, according to [26].Studies indicate that using the method, the mean length of wave groups is too short [26].Thus, the random phase generated waves may not accurately test the PTO's capability of, e.g., power smoothing and absorption, as this is highly dependent on the degree of wave grouping.Instead, [26] suggests that by filtering Gaussian white noise, better representation of ocean waves may be achieved, and arbitrarily long series may be generated.
The white noise method is to design a filter according to the PM-spectrum, whereby the filter shapes the power density of the white noise signal to comply with the PM-spectrum.The filers are implemented based on [11,27].
If the input process to a filter, x(t), is white noise, which is defined by having a flat spectrum, xx , and the filter, H(f ), is designed according to the PM-spectrum in Figure 9a, H(f ) = S PM (f ), then: where a white noise with variance, σ 2 xx = 1, has been used as input.Thus, by filtering white noise, x(t), with the filter, H(f ) = S PM (f ), yields a time series, y(y), having a spectrum, S yy (f ), agreeing with the PM spectrum.
The generated wave with this method is shown in Figure 9b.A spectrum analysis of the generated wave is seen in Figure 9c, showing fluctuation around the target Pierson-Moskowitz-spectrum.This is desired, as spectrum measurements of real sea waves also have a non-smooth spectrum.

Absorber Model
The Wavestar absorbers are modeled using linear wave theory, as described in, e.g., [24], yielding an adequate description in the conditions in which a WEC produces energy.The absorber is illustrated in Figure 10a, representing a single degree of freedom system expressed by its angular motion, θ arm .The angular position, θ arm , is defined to be zero when the float is horizontal, i.e., position at calm water.A derivation of the model for single absorber may be found in [28].The model here is briefly restated for a single absorber with a diameter of 5 m.

Single Absorber Model
In linear wave theory, the equation of motion of the absorber is obtained by superimposing the different effects of the wave-absorber interaction: where the torque, τ Arch , is due to Archimedes force; the torque, τ G , is caused by gravity; the torque, τ rad , is the torque experienced on the absorber from radiating waves; τ PTO is the torque applied by the PTO system; and τ ext is the wave excitation torque.The moment of inertia term, J mech , is the total (mechanical) moment inertia in the θ arm rotation.By identifying the terms in Equation Equation ( 4), the following well-known expression is obtained: where k res is the spring coefficient representing the hydro-static restoring torque; τ res , which is a linearization of the Archimedes force; τ Arch , and gravity, τ g .In Equation ( 4), the radiation term, τ rad , may be expressed as: where the function, k r (t), is the radiation-force impulse-response function and J add,∞ is the added inertia at infinite high frequency.
The wave excitation torque, τ ext , is the torque an incoming wave applies to a float held fixed, including the contribution from the float diffracting the wave.A filter for calculating the wave torque based on the wave, η w , may be obtained given as an impulse response function, h ext (t).The exciting wave torque may then be found by the convolution: The impulse response, h ext , is non-zero for t < 0, rendering the filter non-causal.This mean that the current excitation force depends on the future incident waves.This is partly due to the fact that the waves hit the float before reaching the reference point in the center of the float.Furthermore, the defined wave height is not the direct cause to the wave excitation torque, but is a quantity defined model-wise.Thus, both are caused by some unknown process, and therefore, their relation is not forced to be causal.
To solve this in the simulation, the wave excitation torque is pre-computed before simulation by discrete convolution: where the impulse response is approximated with a finite impulse response, requiring n − samples of future wave knowledge.As the excitation torque is assumed independent of the float motion, but only dependent on the incident wave height, the pre-computation of τ ext is allowed.As the entire wave signal, η w , is available beforehand, the non-causality is not a problem.

Multi-Absorber Model
In an array of closely spaced absorbers, such as the Wavestar, the absorbers interact.The absorbers interact in two ways: • By diffracting the incident wave; • By radiating waves, interacting with neighboring absorbers.
The diffraction part may be obtained by placing 20 floats in the given array and, then, identifying the force filter, h ext , for each float.The resulting diffraction effect experienced by the individual floats will be included in their force filter for calculating the exciting wave torque for each float: Considering the radiation effect, the equation of motion for each float, before given as Equation ( 5), now becomes coupled: as the matrix, K r (t), containing the radiation impulse responses is non-diagonal: e.g., the impulse response, k r 2,20 (t), describes the cross-coupling from the velocity of float number 20 to force on float number 2. Note that k r x,y (t) == k r y,x (t), rendering K r (t) symmetric.
With the focus on validating the PTO, the radiation cross terms are neglected in this paper, i.e., k r x,y = 0 for (x = y).Inclusion of the cross terms mostly affects how the absorbers should be controlled to further increase power absorption.Thus, the influence of the float dynamics on the PTO and the cylinder control may be modeled with acceptable accuracy by using a single absorber hydrodynamic model.However, it is still required to have the effect of the floats operating out of phase due to their distributed locations.This is obtained by applying individual force filters, as mentioned previously.

Model Parameters
To obtain the hydro-dynamic model parameters, WAMIT has been applied, which is a program for computing wave loads and motions of structures in waves.
The force excitation filters are shown in Figure 11 for a wave direction of θ w = 0 • and θ w = −45 • .The wave "measurement" point is shown in Figure 11c, along with a definition of float numbering and wave angle, θ w .The filters show that for θ w = 0 • , the wave reaches the floats farthest back (numbers 10 and 20), approximately 7 s later after exciting float number 1. From the plot, it is also seen that there is some amount of shadow effect, as the magnitude of the wave excitation filter reduces for each float.For the θ w = 0 • , the input for the float is pair-wise equal, whereas this is not true for the wave angle, Note that the excitation torque is independent of the float movement, i.e., the excitation torque may be calculated "off-line" by the convolution in Equation (7).
From WAMIT, all the impulse-response functions, k r x,y (t), describing the radiation-force are also obtained, including the cross terms (x = y).However, as theses are neglected, only the single absorber radiation impulse response is used; see Figure 12a.To avoid performing the convolution term, τ rad,kr = k r (t) × ω arm (t), the term is approximated as a system of Ordinary Differential Equations (ODE).This is performed using Prony's method [28].The used fifth order approximation is shown in Figure 12a.The absorber dynamics may now be written as the following transfer function: A Bode-diagram of Equation ( 12) is given in Figure 12b with the parameters in Table 1.The Bode diagram shows that the system has a resonance peak at f r = 0.285 Hz, corresponding to a period of 3.51 s.Hence, the Wavestar point absorber will have optimal absorption at a wave period of 3.51 s.At Wavestar C5 production sites, this period corresponds to the sea states with the shortest peak-period, in which a substantial part of the yearly wave energy is concentrated.

Absorption System
The PTO concept is modeled as displayed in Figure 8. First, the DDC-module is developed, which connects to the pressure line system.Secondly, the pressure line system is described, which consists of transmission pipe lines and accumulator batteries.

Discrete Displacement Cylinder-DDC
The DDC-module is shown in Figure 13a, consisting of a cylinder, a manifold with valves, accumulators and pipe-lines.

Manifold accumulators
Cylinder V aM x Min x Min x Min

Cylinder
The pressure dynamics of the cylinder chambers are described by the flow continuity equation: where V 0,1 and V 0,2 are the volumes of the connecting hoses and β eff is the effective bulk modulus of the fluid as a function of pressure.The used value of the effective bulk modulus is shown in Figure 14b.The cylinder force is calculated as: where the cylinder friction, F fric,c , is defined, such that the cylinder has a constant efficiency of η c .The function, tanh(av c ), is used instead of sign(v c ) to avoid discontinuity, where a adjusts the steepness around zero velocity.x P1 x P2 x P2 x P2 x P2 x P2 x P2 Accumulator Battery 30bar The moment arm, d a , at which the cylinder, F c , applies τ PTO is found from Figure 10c to be: where: The function is shown in Figure 10c.

Transmission Line Model
Considering the pipe connection between manifold and cylinder, the fast on/off valves that are utilized may produce steep pressure transients, exciting the transmission line dynamics.In other terms, the time required to accelerate the mass of the fluid within the hoses may not be neglected.
The applied pipe-line or transmission line model is to view the fluid within the line as a single mass column, which is accelerated by the pressure difference at the line inlet and outlet; see Figure 13b.The dynamics of the line flow velocity, v l , may then be described using the momentum equation: where p f (Q l ) is the pressure drop due to the total line resistance (fittings, orifices, hose).Pressure drop in a straight pipe/hose, p λ , may be described as: where l l is pipe length, d l is the lines inner diameter and ν is the kinematic viscosity of the fluid.The decision of flow regime is formed using Reynold number: Re: The use of the hyperbolic-tangent expression in Equation ( 20) is for creating a continues transition between linear and turbulent pressure drop for increased robustness of the simulation.Re = 2300 is used as the transition number.The pressure drops, p ξ , of fitting are described as: where ξ if a friction coefficient for a given fitting type.Thus, for a line with n fittings and k line pieces, the pressure drop, p f (Q l ), due to the total line resistance, may be written as:

Manifold with Valves and Cylinder Pipe Connection
The valve flows, Q v xx , are described by the orifice equation, here illustrated for valve vH3 (see Figure 13c): where C d = 0.65 is the discharge coefficient; ρ is the fluid density and A v (α v ) is the valve opening area.The opening and closing of the valve is described as a ramp function, which is also illustrated in Figure 15a: where A o is the valve opening area when fully opened.The areas, A o , are chosen, such that the pressure drop across a fully open valve is 5 bar at the steady state flow occurring when the cylinder displaces with 0.5 m/s.Thus, the actual value of C d is not important in this study, as it is compensated for in the chosen A o .To avoid simulating the dynamics of the pressure downstream of the valve, p v,A3 , which is a small volume with fast dynamics, an algebraic relation between p v,A3 and p MH is assumed, which is given by Equation ( 24).This yields the following combined hose and valve model to describe the flow, QvH3 :

Gas loaded accumulator
else The above "if"-condition is to avoid dividing by zero when the opening area goes to zero, i.e., the pressure resistance goes to infinity.The cylinder flows are calculated as The fitting resistance, ξ M , shown in Figure 13a, is added to take into account internal pressure drops in the manifold.The fittings, ξ c , are the cylinder inlet loss, and ξ A1 is an extra loss due to an extra bend in this particular pipe connection.

Manifold Accumulators
The operating region of charging and discharging of the gas accumulators will be a combination of adiabatic processes for the fast shifting transients and between adiabatic and isothermal processes for the slower energy smoothing cycles.Resultantly, modeling of the gas temperature and its energy exchange with the surroundings is required to get an adequate estimate on the accumulator losses.
The energy balance of the gas in the accumulator may be written as [29]: and the standard gas law: where: h : Heat exchange coefficient between gas and environment A w : Wall area m g : Gas mass T w : Wall temperature v g : Gas specific volume (V g /m g ) V g : Gas volume T : Gas temperature p a : Gas pressure c v : Gas specific heat at constant volume R : Ideal gas constant To simplify the expression, the thermal time constant, τ a , is defined [29]: yielding the equations: To improve simulation speed, the external volume seen by the accumulator is combined into the dynamics of the accumulator; see Figure 15b.The volume and accumulator dynamics are combined by letting Q V = Q a and by adding the external volume, V ext , to the accumulator oil volume.This yields the following flow continuity equation: where V a0 is the size of the accumulator; β eff is the bulk modulus of fluid in V ext and V g = V a0 when p a is below the pre-charge pressure, p a0 .Rewriting Equation (36), Equation (34) and Equation (35) yields: Thus, solved in the order, Equation (40), Equation (37) and Equation (39), these equations yield the state equations of the accumulator, including thermal losses.The initial conditions of the accumulator state are related through the standard gas law: where p a0 is the pre-charge pressure at temperature, T 0 .The pressure, p V , in the volume seen by the components connected to the accumulator may then be found by the expression for the inlet orifice with opening, A a : The inlet flow to the accumulators is given as, The external volume, V ext , is set to half the volume of the connecting transmission lines.Depending on the lines, the volume, V ext , is between 3-6 L.

Hoses from Manifold to Pressure Lines
The hoses are modeled as in Section 6.1.2.Here shown for The model parameters for the DDC-system are summarized in Table 2.The hose connection, l MM , length is 3 m for float 1-10 and 5 m for float 11-20.This is performed to spread out transmission line resonance frequencies, yielding a more steady operating pressure line system.

Pressure Line System
The pressure line system is shown in Figure 16.Accumulator batteries are placed in three locations in the system.The size of accumulator batteries is seen in Figure 14a.The three high pressure accumulator batteries each consist of 16 × 50 L piston accumulators for energy storage.Each intermediate pressure accumulator battery consist of 4 × 50 L, and the low pressure batteries consist of 10 × 50 L.To show how the pressure line system is modeled, a small section is modeled, see Figure 14a, where an example is given for a high pressure line.The same line model as in Section 6.1.2. is used: To describe the fusion of the flows from float 1 and float 11 at the T-fitting, a pressure at volume, As the accumulators are placed close together, they may model-wise be combined into one equivalent accumulator, here shown for the high pressure battery: The model parameters for the pressure line system are summarized in Table 3.

Generation System
In this section, the four identical modules of a motor, generator, inverter and grid inverter are described and modeled.

Hydraulic Motor
The hydraulic motors in the system are all 250 cc fixed displacement closed-circuit bent-axis motors.The model in this simulation is based on the efficiency plots for a Sauer-Danfoss 250 cc Series 51-1 bent-axis motor at full displacement.The efficiency plot is seen in Figure 17a.Using the Schloesser formula [30] for friction and flow loss, the following expression gives motor torque and flow inlet: The efficiency of the fitted model is seen in Figure 17b.Bent-axis motors at frame size 250 cc are allowed to operate at up to 2200 RPM continuously.To replenish the fluid leaked by the motor and to give the necessary flushing and filtering, a charge pump is installed in the system.The pump is set to be 20% of the total motor displacement; thus, D charge = 200 cc.The speed of the pump is controlled by a motor with an inverter.The power consumed for flushing is calculated as: where it is assumed that the flushing pump has an efficiency of 95%, the electric motor an efficiency of 95% and the converter an efficiency of 97%.

Generator and Inverter
The generator setup consists of an asynchronous generator and a converter for a grid connection and to enable variable speed control.The input to the generator is the hydraulic motor torque, and the output of the power system is the angular velocity of the generator, ω Gn .
The transient behavior of the electrical part of the generator is neglected.This is assumed acceptable, as the inverter-generator closed loop dynamics are fast compared to the slow load variation.An equivalent circuit for a phase of a three-phase ∆-connected induction motor is seen in Figure 18b, where γ denotes the motor slip.The slip is defined as: where n pp is the number of pole pairs and ω V is the frequency of the supply voltage.The resistor, R 2 1−γ γ , represents the mechanical power input to the generator, i.e., τ Gn ω Gn = I 2 R 2 1−γ γ ; for reference, see e.g., [31].
As the current, I 2 ≫ I M , the generator torque is given as: where V RMS is the RMS (Root-Mean-Square)-value of the line-to-line voltage.
The steady-state phase current, I P , of the generator is given as: where The electrical output power of the generator is three-times the power per phase: As the chosen hydraulic motors have their best efficiency in the range of 500-2000 RPM, a four-poled generator is used.The chosen generator in the simulation is a asynchronous high-efficiency four-pole 160 kW generator, where the parameters of the equivalent circuit have been obtained.The efficiency result is seen in Figure 18c, where the efficiency is plotted as a function of load at 50 Hz and 400 V.
Back emf Efficeincy [-] (c) The angular velocity of the generator is given by: where J Gn and J M are the moment of inertia of the generator and the hydraulic motor, respectively.The torque of the generator is controlled by an inverter.To control the torque of the generator, the expressions for finding the appropriate voltage and voltage frequency for the generator is implemented in the inverter; see Figure 19.The control of the inverter is set to optimize the generator efficiency, i.e., the voltage is reduced and the slip increased at lower speeds to reduce losses.Using the control, the generator efficiency seen in Figure 17c is obtained, widening the high efficiency region of the generator.
The inverter loses are modeled as an efficiency map η inv (ω V , P inv ).Based on the data-sheet for a 160 kW Danfoss VLT-Drive, a look-up table of the inverter efficiency in Figure19b has been constructed.Finally, the inverter to grid is assumed to have a constant efficiency of η inv,grid = 97%.The parameters are summarized in Table 4.

PTO Control
The PTO control may be divided into two parts, as shown in Figure 20.The first part is the float control, Figure 20a, which controls a DDC-system, compromising between maintaining efficiency and proper power absorption.A Wave Power Extraction Algorithm (WPEA) generates the force reference for optimizing power extraction, while taking into account the PTO efficiency.
The second part is the system control, Figure 20b, which handles the overall control of pressure lines and power generation, i.e., the numbers of active generators and their speed references.The system control may also manipulate the float control to put more or less flow into the mid-pressure line (control parameter, ϕ) or to reduce the power absorption if the pressure lines are saturated.

Wave Power Extraction Algorithm
The power extracted by the absorbers from the wave, P ext , is the product of the applied PTO torque, τ PTO , and arm velocity, θarm (or -F c • v c ).To maximize the extracted power, the float velocity should be in phase with the exiting wave torque, τ ext , i.e., the natural frequency of the float and arm should match the incoming wave.However, frequency-wise point absorbers are narrow-banded with an under-damped resonance frequency, as seen in Figure 12b.Accordingly, point absorbers are prone to operate off-resonance, as the wave period varies wave to wave.On a larger time scale, the average wave period also varies from sea state to sea state, as seen for the wave spectra in Figure 9.
To improve the behavior of point absorbers, it is well established that control of the PTO load force may be used, increasing the amount of energy extracted from waves [4,5].Basically, the torque, τ PTO , applied to the absorber is controlled as a feedback of the absorber's motion as: allowing the PTO to adjust the absorber's resonance frequency to match the wave frequency.The term, k PTO θ arm , is the reactive term and imposes the PTO to implement a virtual spring element, changing the external experienced dynamics of the absorber.Optimally, the feedback ensures that the velocity of the absorber is always in phase with the wave excitation force and the PTO damping matching the absorbers' hydrodynamic damping [4].However, this requires future knowledge of η w for some time, t future , rendering the control non-causal.The required wave prediction horizon to remove the non-causal part is approximately the same as the settling time of the system [4].
A causal and more robust control described for Wavestar in [28] is to continuously gather statistically information of the current sea state over a window of 100 or more waves.The absorbers narrow frequency response is then widened and adjusted to a fixed resonance frequency, yielding the best average abortion in the current sea state.This may be implemented by having fixed control coefficients in the feedback law Equation (60).This is a suboptimal approach, but has been applied with good performance results, as reported in [32].
Optimal adjustment of the absorbers frequency, as the above, requires the PTO to transfer energy to the absorber (assisting its movement) in parts of an oscillation cycle.Thus, the above WPEAs are termed reactive strategies and require a PTO capable of four-quadrant behavior to supply the reactive power.If reactive control is not offered by the PTO, the reactive terms of the control are removed and damping control remains, i.e., for linear damping control only the coefficient, b PTO , is non-zero in Equation (60).Linear damping is only optimal when the wave frequency and the absorber's resonance frequency match by other means.
For the Wavestar concept, it has been verified from prototype measurement [22] that increased production is obtained through reactive control; see Figure 21a.Each point is an average of extracted power over a 10 min measurement, where both linear damping and reactive control have been tested.The measured average power results are given as a function of significant wave height.Consequently, some of the scatter is due to the different wave periods.The applied control is as in Equation ( 60), where the coefficients are fixed for a given sea state.The system continuously gather statistically information of the current sea state over a window of 100 or more waves and, then, choose the coefficients yielding the best average abortion in the current sea state.

Figure 21.
(a) average extracted power measurements for a single absorber of the prototype in Figure 5a [22].Each point is a 10 min measurement; (b) comparison of power extraction of linear damping and reaction as a function of efficiency in a sea state with H m,0 = 1.75 m and T p = 4.5 s [11]; and (c), yearly production of a single absorber at the prototype site [28] as a function of PTO torque limitation and efficiency [11].
When optimizing reactive control schemes, it is important to take into account PTO power conversion efficiency, as power is lost in each conversion.Hence, having reactive power oscillating between float and PTO consumes energy.The control in Equation (60) should always optimize the power output to the grid and not the amount of extracted power.
In [8,28], it is shown how the PTO force constraints and conversion efficiency may be included into the control design.The results are a map of control parameters, b PTO and k PTO , as a function of sea state and PTO efficiency.For detailed information on the methods, see [11].The optimal parameters in [8,28] are found by numerically optimizing the output of the simulation model shown in Figure 22.The model includes the absorber dynamics, and the input to the model is an irregular sea wave, η w .The PTO is characterized by a force limitation and an efficiency.The PTO power output may then be formulated as a function of the extracted power, P ext (t): The optimal control parameters for Equation ( 60) is then found as the arguments maximizing the energy output of the PTO: The optimal values of b PTO and k PTO for two sea states are shown in Figure 22b.These results are from [11].Here, it is seen as expected, that when the PTO conversion efficiency, η PTO , is reduced, the reactive part (k PTO ) reduces and the resistive part increases.To further justify the reactive control, a comparison of reactive control and linear damping is shown in Figure 21b [11].It shows that despite having a power conversion efficiency of, e.g., 70%, 50%, more power is produced to the grid compared to linear damping.
The torque requirement of the PTO is chosen to be 1 MNm, which is a compromise made according to Figure 21c, as the PTO efficiency is expected to be about 70%.In Figure 21c, the yearly production for a single absorber has been computed in [11] by calculating the average power output for all sea states and, then, assuming a yearly wave distribution.
To summarize, given the sea state and PTO information, (H m,0 , T p , η PTO ), the force reference to the cylinder is calculated as: where d a is the moment arm and γ is a coefficient, which the system control may use to reduce the power absorption if the pressure lines are saturated.To reduce power absorption, a method could also have been to reduce the efficiency input to the WPEA algorithm, forcing the WPEA to move towards linear damping, thereby reducing power absorption and system stress.The actual PTO efficiency is not independent of the WPEA defined in Equation (63).Adjusting the reactive power may either improve or reduce component performances, depending on whether the increased power level brings the components closer to their rated/optimum power level.To this account, an initial efficiency estimate is given to find the WPEA, and through simulation of the PTO, a new efficiency is obtained, which is used to update the WPEA.A number of iterations have been performed to find the proper settings.This is also discussed in [8] for the Wavestar converter.
The force control is going to be implemented discretely, i.e., the continuous reference generated by Equation ( 63) is subjected to a quantification.To investigate the influence, comparison of a continuous control and a discrete control has been performed.
The continuous control is simulated as in Figure 22a, with the optimum parameters, a PTO efficiency of 70% and a maximum torque of ±1 MNm.The simulation of the discrete version is implemented likewise with the same WPEA parameters, but the PTO torque is quantified before being applied to the float.The available torque values besides zero are set to 9 evenly spaced values between ±1 MNm.Moreover, a limit is set on the shift rate, limiting the force to shift at maximum one time per 300 ms.Thus, after a shift has occurred, the force value is locked for 300 ms, whereafter the quantification process may choose the force closest to the continuous reference.This avoids the risk of jittering and PWM (Pulse Width Modulated)-like behavior and keeps the number of shifts at a reasonable level.The results are seen in Figure 23, showing the two power matrices of average produced power.As seen, the difference is relatively low despite the relatively rough quantization.Thus, the force approximation of the DCC having 27 force values will be adequate.The reason for the low impact of the discrete variation is that the 35 tons of the float and arm system serve as an effective low pass filter, having a break frequency around 0.28 Hz, c.f. the earlier calculated natural period of 3.5 s.Thus, from a load control point of view, implementation of extra dampers/low pass filters is not required.
Another aspect is the compliance of the mechanical structure when performing the rapid load force changes; however, these issues are not treated in this paper.One aspect to notice, though, is that the force steps are relatively small, e.g., a single step from zero to maximum is never performed.

DDC Control-Force Shifting Algorithm (FSA)
The purpose of the Force Shifting Algorithm (FSA), controlling the DDC, is to choose the appropriate force level to approximate the reference, F ref , generated by the WPEA.However, always shifting to the force closest to the continuous reference, F ref , may not energy-wise be optimal, as some force-shifts are more energy-expensive than others, due to the compression loss, E β .The minimum compression loss of a volume when shifting ∆p, independent of the pipe line dynamics [33], is given as [16]: This loss is illustrated in Figure 24a for a system consisting of four system pressures and a symmetric cylinder.The figure shows the compression loss in [kJ] from shifting between forces.The forces are sorted according to size.As illustrated, if the current force is force number 9 and an increase in force is desired, it is cheaper to shift to force number 11 instead of 10.Likewise, from number 11, it is cheaper to shift to number 14 and skip 12 and 13.Thus, to avoid doing very expensive force shifts, a more suitable strategy is to calculate the energy expense of possible force-shifts and make a compromise between tracking and energy-cost.This is the basis of the FSAstrategy developed in [15] and which is also used in this study.Note that the loss matrix changes with pressure and cylinder position.The shifting loss from shifting from force x to y is found summing the losses Equation (64) for pressure changes in the individual volumes: Note that the pressures, p Ax,old = p Ax,new , may be equal, as the pressure is not necessarily shifted in all chambers.
To allow compromising between force tracking and performing less expensive force shifts, a maximum allowed tracking error, F b , is defined, meaning that F c must stay within a band of ±F b about F ref .However, within the band, the FSA may choose the force steps with the lowest shift cost.A fixed time limit, T min , on how frequently shifting is allowed is also added to reduce tracking cost.
The FSA algorithm is illustrated in Figure 24b.The values, k − and k + , are the numbers of the forces, which are the cheapest to shift to within the band, ±F b .As seen, these are continuously updated.When the force reference comes closer to one of the forces, k − or k + , a shift is performed, and a new lock down period of T min is initiated.
Based on the optimization procedure in [15], it has been found that F b = 150 kN and T min = 350 ms for optimizing the power output of the Wavestar C5.
Thus, when a force shift is initiated, the control sends out a matrix, u c , of control values for the nine on/off valves, corresponding to the desired pressure configuration.Regarding valve timing, it was found in [16] that a small amount of overlap between opening and closing of the valves for a single volume was desirable; thus, a 3 ms overlap is used.As the valves have 12 ms opening and closing time, the signal to the opening valve is delayed 9 ms.This is illustrated in Figure 24c.
To help in controlling the pressure in the mid-pressure line, the FSA continuously identifies the force combinations that would currently supply or consume flow from the mid-pressure line.For these configurations, an "artificial" energy loss, ϕ, may be added in Equation (65) by the system control to either penalize supplying or consuming flow from the mid-pressure line.

System Control
The purpose of the system control is to: 1. Avoid the high pressure accumulator storage from depletion or saturation; 2. Keep the mid-pressure line floating between high and low pressure; 3. Ensure as steady a power production as possible, while satisfying 1 and 2; 4. Choose the proper number of generators for a given sea state; 5. Reduce power absorption when full load capacity is reached.
To maintain a stable power production, the generators are set to initially produce the expected average power, p avg,expected , in the current sea state.An initial guess is given based on the current sea state when starting production, whereafter a moving average is used based on the absorbed power over a window of 5 min.The number of active generators, k GN , is then chosen, such that generation capacity is roughly p avg,expected plus 30%.All active generators are operated at the same speed ω GN .
To avoid the high pressure accumulator storage from depleting or saturating, the power generation is increased or decreased based on the pressure in the accumulators.This is performed through the coefficient, ψ.The speed reference is then given as ω GN : where D M is the total active motor displacement in [m 3 /rad].As p avg,expected is the absorber power, the efficiency from cylinder to power out of the generator, η to-Gn , is required to calculate the average generator power.The maximum allowed speed is set according to 180 kW per generator, which is 15% overload.The lowest speed is set to 400 RPM.The high pressures are set to be between 150 bar and 300 bar; thus, the map in Figure 25 for ψ is used.Finally, if pressure is still reaching 300 bar, the value, γ, is manipulated, such that the float reduces power absorption; see Figure 25.This is similar to turbines pitching out of the wind when rated production is reached.The mid-pressure is to be floating between high and low pressure for optimizing the DDC efficiency.Hence, the "reference", p M,ref , for the mid-pressure line is: Based on this, the penalty value, ϕ, is set according to Figure 25.

Results and Discussion
The modeled PTO system is evaluated for the three sea states, 1, 2 and 3, which were defined in Figure 9.These correspond to small, medium and high production cases, respectively.To evaluate the PTO, a test wave has been created using the described white noise method.The first 100 periods of the wave have been checked to have a representative energy content complying to the PM-spectrum.Additionally, the test wave contains both high intensity peaks and calm periods between wave groups.To compare the simulations, this wave is reused for all three sea states by scaling it in time and height correspondingly.The results are shown for an incoming wave angle of θ w = 0 • .The PTO performance is nearly independent of wave direction, as the accumulators handle the wave-to-wave power fluctuation.
First, the PTO is overall evaluated for the three sea states.Afterwards, the cylinder force tracking, and control performance is inspected.

Overall System Performance
In Figure 26, approximately 6 min have been simulated of the complete Wavestar C5 for sea state 1.The power absorption is the instantaneous mechanical power input to the 20 cylinders.For sea state 1 with an average power production of about 120 kW, two generators are active, namely, generators 1 and 3, such that pipe losses are reduced.
As seen in Figure 26, the PTO system is capable of smoothing out the varying power absorption using the accumulators, whereby the power output to the grid is fairly stable.Note that at around times, t = 50 s and t = 150 s, long periods of almost no absorbed power are experienced, but the system is able to maintain output.The results show that all available "peak power" of the waves is extracted, as the high pressure, and, thereby, the accumulators' storage, never saturates [p H (t) < 320 bar].The control system stabilizes the mid-pressure between high and low pressure as desired.Note that at small waves, the high pressure level is kept a bit lower to reduce compression losses, as the highest cylinder forces are not required during smaller waves.The generators are operated at the same speed, which is seen in the bottom of Figure 26.The speed is slowly varying between 850 to 1,200 RPM; likewise, the generator input torque is also stable around 800 Nm.Thus, the operating conditions for hydraulic motors, generators and inverters are very favorable.Hence, the PTO strategy of using the hydraulic accumulators and DDCs effectively separates the generation from the complicated and highly fluctuating wave power absorption process.Looking at the power absorption of a single float, the factor between peak and average power is 12, whereas the fluctuation of output power is within ±50% of mean output.
In the simulation for sea state 2 in Figure 26, all four generators are active.Compared to the PTO behavior in sea sate 1, the storage capacity is now almost fully used, with the high pressure varying between 160 bar and 300 bar.This means that the full storage capacity of the accumulator battery is being utilized.The system completely absorbs the 10 s long 2 MW peak, while maintaining a steady production around 400 kW.
Moving to sea state 3 (Figure 27), the system is operating at full load most of the time.As seen from the pressure, the storage saturates; however, the system is kept below the maximum of 320 bar by reducing the load force on the absorber.Thus, no extra energy has to be dissipated internally in the system using, e.g., relief valves, although these are installed.Looking at the power absorption of a single float, the mean power is 45 kW and the peak power is 600 kW, yielding a factor of 12 between mean and peak power.

P DDC
To calculate efficiency accurately, it is required to take into account that the stored pressure energy in accumulators and oil volumes is not the same at the initial state and at simulation end.Resultantly, a post-simulation is performed on the pressure line system and manifold using the end states.Here, how much energy is required to fill or drain the pressure lines to get to the initial state is simulated, which is illustrated in Figure 28b.The difference in energy is then mapped into an equivalent average power input, P ∆init , to the pressure line system.
In Table 5, the total PTO efficiency, η tot , is defined as the ratio between total absorbed energy and the energy to grid, taking into account the difference in the stored energy: In the above expression, the difference in stored energy, P ∆init , is scaled with 1 η DDC , such that it seems that it has been absorbed from the waves.
From Table 5, it is seen that the total PTO efficiency, η tot , is 70%, 73%, 71% for the three sea states, respectively.Thus, despite the varying power levels and a factor of 10 between the mean and peak power, the PTO system is able to maintain efficiency.
The motors and generator all operate around 94% and 95%, as expected in the operating conditions.Most important in the results is that the DDC-system is able to maintain an overall efficiency of η DDC = 90%, which is the core technology for enabling the PTO concept.The actual power conversion efficiency is actually higher, as reactive control is being performed, making the average processed power larger than the average input power.This may be viewed from the single float power absorption, which is periodically negative.The average negative power, P ext,-, and average positive, P ext,+ , power have been calculated for float 1 in the three sea states: • Sea state 1: P ext,-= 0.57 kW and P ext,+ = 7.5 kW: Ratio of reactive and real power: c = Pext,+ Pext,-= 0.075; • Sea state 2: P ext,-= 2.5 kW and P ext,+ = 19.5 kW: Ratio of reactive and real power: c = Pext,+ Pext,-= 0.088; • Sea state 3: P ext,-= 3.1 kW and P ext,+ = 49.0 kW: Ratio of reactive and real power: c = Pext,+ Pext,-= 0.063.
Assuming equal power conversion efficiency in both direction, the relation of the actual power conversion efficiency of the DDC-system, η act , and the 'measured' average efficiency, η DDC , is given as: The actual power conversion efficiency for the manifold and cylinder is, thus, 90%, 92% and 91% for the sea states, 1, 2 and 3, respectively.

Control System Performance
In Figure 29, the operation of the cylinder and manifold of float number 1 is shown.The control tracks the force reference with the discrete steps adequately and is excellent at not performing extra unnecessary shifts, when the variation in force reference is small, for example, at time 160 s and at time 166 s.
Looking at the pressures within the three chambers, the shifting algorithm optimizes efficiency by mostly shifting one pressure at a time and using the mid-pressure as an intermediate level.The chamber, A2, is the largest chamber and, thereby, the most energy expensive to shift.This is seen in the control output, as A2 is the chamber that experiences the lowest number of pressure shifts, whereas the smallest chamber, A3, experiences the most pressure shifts.
In the zoom in on the pressure, the shifting transients in pressure may be seen.The fast oscillatory behavior in the transients are due to the fact that the flow may oscillate between the cylinder chamber and manifold accumulator, due to the pipe line between the valve and cylinder.These are most pronounced for the A1 chamber, as this chamber has the longest pipe connection, thereby increasing the effect of pipe "inductance".The slower oscillation in pressures p H and p M are likewise caused by flow oscillating between the small manifold accumulators and the large storage accumulator batteries.In Figure 30, the system control performance is shown.The mid-pressure reference, p M,ref , is set between low and high pressure, and the mid-pressure is effectively kept floating around this reference by manipulating the force shifting algorithms.This was implemented by giving a penalty for shifting to force configuration, which will take the mid-pressure further away from the reference.The penalizing "energy" input, ϕ, to the force shifting algorithms is seen in the figure.

Figure 2 .
Figure 2. (a) a Wavestar absorber; (b) velocities and PTO load in typical production waves; and (c), instantaneous power during two minutes production for a single absorber.

Figure 3 .
Figure 3. (a) the definition of Power Take-Off (PTO) system; (b-d), different fluid power based transmission in PTO systems.

Figure 4 .
Figure 4. (a) the discrete PTO system based on discrete displacement control of a multi-chambered cylinder (Discrete Displacement Cylinder (DDC)-system); (b) the forces, F PTO , the PTO may produce; and (c) how F PTO is discretely varied using the DDC-system during a wave.

Figure 6 .
Figure 6.(a) an overview of the PTO concept for the Wavestar WEC in Figure 5b; (b) the DDC-system, consisting of a multi-chambered cylinder with integrated shifting manifold for discrete throttle-less force control; (c) a sketch of the shifting manifold; and (d) the illustration of the 27 available forces.

Figure 7 .
Figure 7. (a) an example of how pressure and chamber connection may change during a wave; (b) the available cylinder forces.The numbers 1, 2 and 3 relate to which pressure is in the different cylinder areas in (a) (1 = p L , 2 = p M , 3 = p H ).

Figure 8 .
Figure 8. Modeled sub-systems of the Wavestar WEC and PTO system.

Figure 9 .
Figure 9. Wave spectra for three sea states and an example of a realization of a wave.

Figure 10 .
Figure 10.(a) definition of variables of the Wavestar absorber; (b) some dimensions for the case study; and (c) the moment arm of the cylinder.

Figure 13 .
Figure 13.(a) diagram of DDC-system, consisting of cylinder, manifold with nine on/off valves, accumulators and pipe lines.Fitting are indicated, where ξ denotes the fitting resistance coefficient; (b) illustration of hose model; and (c) illustration for combined valve and hose model.

Figure 14 .
Figure 14.(a) modeling a part of the pressure line system; and (b) the used bulk modulus model.

Figure 15 .
Figure 15.(a) the valve model; and (b) the accumulator model.

Figure 17 .
Figure 17.(a) efficiency plots (from datasheet) of a Sauer-Danfoss 250 cc Series 51-1 bent-axis motor at full displacement; (b) the efficiency of the modelled motor; and (c) the efficiency of a generator when inverter control is configured for optimal efficiency.

Figure 18 .
Figure 18.(a) Delta-connected induction motor; (b) per phase equivalent circuit; and (c) efficiency as a function of load relative to rated torque at 50 Hz and 400 V.

Figure 19 .
Figure 19.(a) torque control of the generator and inverter model; and (b) the efficiency map used for the inverter.

Figure 20 .
Figure 20.(a) float control (20 parallel system), handling the manifold control to track the generated force reference to extract wave power; and (b) system control for controlling generation and pressure lines.
power over 10 min.

Figure 22 .
Figure 22.(a) model for reactive control optimization [8]; and (b) control parameters for two sea states as a function of PTO power conversion efficiency, η PTO .

Figure 23 .
Figure 23.Power matrices of average produced power in [kW].The matrices are shown for a continuous torque control and a discrete control with ±4 force values. 0

Figure 24 .
Figure 24.(a) illustration of the shifting loses due to compression; (b) illustration of the force shifting algorithm; and (c) the opening and closing procedure of the valves during shifting.

Figure 30 .
Figure 30.Controller performance and gas temperature of accumulators.

Table 2 .
Values used for the float module, where "Par." and "Val."abbreviate parameter and value, respectively.* hose connection length is 3 m for float 1-10 and 5 m for float 11-20.

Table 3 .
Parameter values used for the pressure line system, where "Par." and "Val."abbreviate parameter and value, respectively.