Next Article in Journal
Trends in CO2 Emissions from China-Oriented International Marine Transportation Activities and Policy Implications
Next Article in Special Issue
Validating a Wave-to-Wire Model for a Wave Energy Converter—Part II: The Electrical System
Previous Article in Journal
Multi-Objective Scheduling Optimization Based on a Modified Non-Dominated Sorting Genetic Algorithm-II in Voltage Source Converter−Multi-Terminal High Voltage DC Grid-Connected Offshore Wind Farms with Battery Energy Storage Systems
Previous Article in Special Issue
Optimization of the Runner for Extremely Low Head Bidirectional Tidal Bulb Turbine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Validating a Wave-to-Wire Model for a Wave Energy Converter—Part I: The Hydraulic Transmission System

1
Centre for Ocean Energy Research, Maynooth University, Maynooth, Co. Kildare, Ireland
2
Department of Mechanical Engineering, University of Bath, Claverton Down, Bath BA2 7AY, UK
*
Author to whom correspondence should be addressed.
Energies 2017, 10(7), 977; https://doi.org/10.3390/en10070977
Submission received: 24 May 2017 / Revised: 3 July 2017 / Accepted: 6 July 2017 / Published: 12 July 2017
(This article belongs to the Special Issue Marine Energy)

Abstract

:
Considering the full dynamics of the different conversion stages from ocean waves to the electricity grid is essential to evaluate the realistic power flow in the drive train and design accurate model-based control formulations. The power take-off system for wave energy converters (WECs) is one of the essential parts of wave-to-wire (W2W) models, for which hydraulic transmissions are a robust solution and offer the flexibility to design specific drive-trains for specific energy absorption requirements of different WECs. The potential hydraulic drive train topologies can be classified into two main configuration groups (constant-pressure and variable-pressure configurations), each of which uses specific components and has a particular impact on the preceding and following stages of the drive train. The present paper describes the models for both configurations, including the main nonlinear dynamics, losses and constraints. Results from the mathematical model simulations are compared against experimental results obtained from two independent test rigs, which represent the two main configurations, and high-fidelity software. Special attention is paid to the impact of friction in the hydraulic cylinder and flow and torque losses in the hydraulic motor. Results demonstrate the effectiveness of the models in reproducing experimental results, capturing friction effects and showing similar losses.

1. Introduction

Wave energy is one of the potential suppliers of clean energy in a future de-carbonised power sector along with other more developed sources, such as solar-, wind- or hydro-power. However, if wave energy is to compete with these more developed energy sources on an equal footing, it is necessary that wave energy maximises the final energy output to minimise the cost of energy.
To analyse realistic final energy generation and design accurate model-based control strategies for wave energy converters (WECs), it is crucial to examine the whole chain from waves to the electricity grid. Hence, accurate mathematical models that incorporate all the necessary components of the different conversion stages from waves to the grid, known as wave-to-wire models (W2W), are necessary. Figure 1 shows the different stages of that chain that need to be incorporated in the W2W models, including the control inputs for each stage ( α , β and γ ). The most important aspects of all conversion stages in W2W models are reviewed in [1] for different power take-off (PTO) configurations: pneumatic [2], hydraulic [3] or mechanical transmission systems [4] coupled to a rotational generator, or even direct drive systems with linear generators [5].
Although diverse PTO systems have been suggested by different wave energy researchers and developers, hydraulic PTO transmission systems appear to be the choice of the vast majority of developers, including Pelamis [6], Searev [7], Wavestar [8], Oyster [9], CETO [10] or Waveroller [11]. The present paper focuses on wave-to-wire models with hydraulic transmission systems.
The W2W model validated in the present paper is designed as a combination of inter-connected sub-models, where each sub-model can take a different time-step by implementing multirate time-integration schemes [12,13]. That way, components with faster dynamics can take shorter time-steps without adding unnecessary computations to the slower components.
The large number of components included in W2W models makes the validation of W2W models, as a whole, quite challenging, mainly due to the complexity and cost of building a physical model with all the components. Therefore, the authors decided to accomplish the validation by separately validating each conversion stage of the model, which is consistent with the way the W2W model is designed.
The hydrodynamic model of the absorption stage is validated against high-fidelity software in [14], where a heaving point absorber is examined when moving with and without control. The validation of the mathematical models for the transmission, generation and conditioning stages is divided into two parts and is presented as a series of two papers. The first paper (Part I), the present paper, shows the validation of the transmission stage, as highlighted in Figure 1, and the second paper (Part II) presents the validation for the generation and conditioning stages [15].

1.1. Hydraulic Transmission Systems in Wave Energy

Several topologies of hydraulic system have been suggested for wave energy converters in the literature. The most simple hydraulic transmission systems include a hydraulic cylinder, high- and low-pressure accumulators, a rectification bridge and a fixed-displacement hydraulic motor, as suggested in several studies, such as the W2W models presented in [7,16,17,18], and also in [19], where a W2W model for arrays is proposed.
In order to improve the performance of that ‘standard’ system, alternative configurations have been suggested in the literature. Four configurations with higher flexibility are presented in [20], where the main difference among the four configurations is the arrangement of the accumulators. Another topology with direct connection between the hydraulic cylinder and a variable-displacement motor is suggested in [21], while a configuration using a multi-chamber hydraulic cylinder connected to multiple accumulators, with different pressure settings (which, in turn, are connected to a fixed-displacement motor), is presented in [22].
All these topologies can be organised in two main groups [23]: constant- and variable-pressure configurations.

1.1.1. Constant-Pressure Configuration

The typical constant-pressure configuration is described in [24]; and generally includes a hydraulic cylinder, a set of check-valves for the rectification of flow, low- and high-pressure (LP and HP, respectively) accumulators, a hydraulic motor, a relief valve, and a boost pump powered by an electric motor, as illustrated in Figure 2. The reciprocating motion of the WEC and, as a consequence, of the cylinder piston, is rectified by means of the rectification valves, so that the hydraulic motor operates in a single quadrant, meaning that the direction of the motor rotation and torque is always the same. The HP and LP accumulators dictate the pressure difference in the cylinder and help to provide a smooth and slowly varying torque in the hydraulic motor. Relief valves only open when the pressure in the system exceeds the maximum pressure allowed by the different components of the system and the boost pump ensures that the pressure in the system never falls below a pre-defined minimum pressure level, avoiding undesirable effects such as cavitation.
Hydraulic transmission systems based on the constant-pressure configuration have, in general, quite limited controllability. The only controllable variable in such constant-pressure hydraulic systems ( α in Figure 1) is the displacement of the hydraulic motor. However, the use of large HP accumulators connected to the hydraulic motor constrains the whole system to perform in constant-pressure, meaning that variations in the motor displacement take a long time to affect the pressure in the hydraulic cylinder, and, eventually, the behaviour of the WEC. As a consequence, the overall power absorption in the wave energy converter can be rather poor, with limited possibilities for improvement. Only the use of extra (smaller) accumulators beside the hydraulic cylinder, as described in [16], can improve the power absorption, implementing control strategies such as latching or declutching, as shown in [25].
In spite of the limited controllability, the efficiency of constant-pressure hydraulic PTO systems can be reasonably high, since components can be controlled to perform close to the optimal operating point most of the time.

1.1.2. Variable-Pressure Configuration

The typical alternative to the constant-pressure configuration is shown in Figure 3, where the hydraulic cylinder is directly connected to a variable-displacement hydraulic motor. Only the LP accumulator is included in this configuration, which, along with the boost pump, avoid pressure drops in the low-pressure line, which can lead to undesirable situations, such as cavitation. A boost pump is also necessary to replenish the fluid leaks in the motor. Relief valves protect the system from over pressure, likewise in the constant-pressure configuration. In the variable-pressure PTO configuration, the reciprocating motion of the WEC and the cylinder piston is not rectified, obliging the hydraulic motor to operate, at least, in two-quadrants. The most common two-quadrant operation of the PTO system in WECs is that of a reversing torque and an irreversible rotation direction.
The absence of an HP accumulator provides a more flexible hydraulic system, increasing the control bandwidth. In the most basic form of a variable-pressure topology, cylinder ports are directly connected to the motor ports without any rectification system, i.e., the check-valves used in the constant-pressure configuration. Hence, variations in the hydraulic motor displacement rapidly affect the behaviour of the WEC, which allows for an effective control of the PTO force applied on the WEC to maximise power absorption from ocean waves. However, the motor may require four-quadrant operation capabilities to compensate for bi-directional flow.
The system shown in Figure 3 is close to the simplest version of a variable-pressure configuration. Similar configurations have been suggested in [20,21], including an energy overflow system connected to the relief valves. That way, when the cylinder produces more flow than the hydraulic motor can handle, the extra flow can be used to generate electricity, instead of leaking it to the LP accumulator, wasting energy contained in the HP fluid.
Further improvements of hydraulic transmission systems for wave energy converters combine the benefits of the constant-pressure and variable-pressure configurations, such as incorporating multiple HP accumulators with different pressure levels to achieve discrete force levels adapting the displacement of a variable-displacement motor [20]. Similarly, Ref. [22] incorporates multiple accumulators, but uses a fixed-displacement motor, where discrete force levels are achieved by means of a multi-chamber cylinder.
The mathematical models for the two configurations are significantly different, thus requiring separate validation. Therefore, experimental data, generated by using two different test-rigs, are used in the validation. The performance of the hydraulic cylinder for both hydraulic system configurations is validated against experimental results, while the well-known simulation software AMESim (v14.1, LMS/Siemens PLM Software, Leuven, Belgium) [26] is used to validate the performance of the hydraulic motor, due to the lack of experimental results for the hydraulic motor.
The paper is organised as follows: the mathematical models are described in Section 2, the experimental test-rigs are described in Section 3, results of the validation are given in Section 4; and conclusions of the validation are drawn in Section 5.

2. Mathematical Models for System Components

Mathematical equations to model the different components included in hydraulic transmission systems, i.e., hydraulic cylinders, valves, accumulators and hydraulic motors, are common to any hydraulic system configuration or topology. Therefore, the equations for the different components are given in the following subsections, regardless if such components are used in the constant- or variable-pressure configuration. The model structures used for each component are generic, in the sense that components of different characteristics, e.g., a hydraulic gear motor or a hydraulic radial piston motor, can be modelled using the same model structure. However, due to the very specific characteristics of each component, the parameters of the models may need to be fitted for each specific component.

2.1. Hydraulic Cylinder

Hydraulic cylinders, symmetric or asymmetric, consist, in general, of a piston that divides a cylinder into two chambers. Hence, pressure difference between the two chambers, imposed by the pressure in HP and LP accumulators in constant-pressure configurations and the torque-control implemented in the hydraulic motor in variable-pressure configurations, determines the cylinder force ( F P T O ) applied on the WEC. The pressure dynamics can be described by means of the continuity equation given by the following expression:
p ˙ A = β e f f ( p A ) V A + A p x p ( Q A x ˙ p A p ) ,
p ˙ B = β e f f ( p B ) V B A p x p ( x ˙ p A p Q B ) ,
where p A is the pressure in the cylinder chamber A, Q A the flow entering or exiting the cylinder chamber A when leakages are neglected, β e f f ( p A ) the effective bulk modulus of the fluid and V A the minimum volume (calculated when the piston reaches its minimum or maximum position) in chamber A. Subscripts A and B refer to chambers A and B in Figure 2 and Figure 3. In addition, A p is the piston area, and x p and x ˙ p are the piston position and velocity, respectively.
Equations (1) and (2) ignore both friction and inertia effects, which can represent between 5% and 10% of energy losses in the cylinder. As a result, the piston force only considering the pressure difference can result in an under-estimate of the force applied to the WEC. A more accurate representation of this force is given by
F P T O = A p ( p B p A ) + F f r i c + F i n ,
where F f r i c is the friction force and F i n the inertia force. F i n includes the inertial force of the piston, the rod and the oil, and the gravity force due to the mass of the piston and the rod:
F i n = x ¨ p ( M p + M r + M o i l ) + ( M p + M r ) g ,
where M p , M r and M o i l are the mass of the piston, rod and oil in the chambers, and g is the acceleration due to the gravity.

Friction Model

Friction is considered to be one of the main nonlinearities of the cylinder. Several different models have been suggested in the literature to describe friction losses in hydraulic cylinders, including static and dynamic models [27].
The vast majority of the models divide the friction model into different effects, i.e., viscous friction, Coulomb friction and static friction. A model, known as Stribeck model, combines the three effects as follows [28]:
F f r i c = σ x ˙ p + s i g n ( x ˙ p ) F c + F s t e x p - | x ˙ p | c s t ,
where σ is the viscous coefficient, F c the Coulomb friction force, F s t the static friction force and c s t the characteristic velocity of the Stribeck curve. Figure 4a illustrates the contribution of each friction effect normalized against the maximum friction force.
Once the model structure is selected, the parameters ( σ , F c , F s t and c s t ) of the model need to be defined. The most appropriate way to proceed would be fitting such parameters using experimental data from the cylinder to be implemented in real life. However, experimental results are not always available, so parameters need to be identified using other relevant information, such as manufacturers’ data. Therefore, the parameters for the present model are identified using an efficiency map as a function of velocity and pressure difference presented in [20], using the nonlinear least squares fitting technique. Figure 4b shows the efficiency map of the hydraulic cylinder, reproduced with the model parameters fitted from [20].
The friction model presented in this section assumes friction losses are symmetric, as shown in Figure 4a. However, this is not always true, especially in asymmetric cylinders. In some cases, a lag exists between the pushing (positive velocity) and retracting (negative velocity) forces. To include such asymmetry, the model structure presented in Equation (5) is still useful, but needs to be duplicated to distinguish the performance during the pushing and retracting operations, separately identifying the parameters for each operation [28].

2.2. Valves

Valves are essential components for the successful performance of hydraulic PTO systems, required to rectify or control the flow at different points of the circuit. Active and passive vales can be included in the hydraulic system, depending on the purpose of the valve, and both can be described using the orifice equation [28]:
Q v = C d A v ( Δ p v ) s i g n ( Δ p v ) 2 ρ o i l | Δ p v | ,
where Q v is the flow through the valve, C d the discharge coefficient, A v ( Δ p v ) the valve opening area as a function of the pressure difference between the outlet and inlet ports of the valve, ρ o i l the density of the hydraulic oil and Δ p v the pressure difference between the inlet and the outlet of the valve.
Only passive valves, check valves in the rectification bridge and relief valves are used in the present study, which open and close as a function of the pressure difference across the valve, as expressed in Equation (6). The valve remains closed while the pressure difference across the valve is lower than the cracking pressure ( Δ p c r a c k ). Once the the cracking pressure is reached, the valve starts to open and continues opening until the valve is fully open. The pressure at which the valve is fully open is known as the maximum pressure ( Δ p m a x ). Hence, maximum area of the valve, and cracking and maximum pressure are, in general, provided in the manufacturers’ catalogues. However, the valve opening, from fully closed to fully open, can follow different profiles. Some of the most commonly observed profiles [29] are illustrated in Figure 5. In the model validated in this paper, the step function is implemented to reproduce the same valve opening profile as in the original study, where the experimental data is taken from [30].

2.3. Accumulators

The performance of hydraulic accumulators, typically compressed gas accumulators, can be described by means of an isentropic and adiabatic process, where volume of the hydraulic fluid contained in the accumulator changes in time as follows:
V a c c ( t ) = 0 t Q a c c i n ( t ) d t ,
where Q a c c i n ( t ) is the input flow to the accumulator. Hence, following isentropic compression, the pressure in the accumulator is given by
p a c c = P p r e V t o t V g a s γ ,
where p p r e , V t o t and V g a s are the pre-charge pressure, the total volume and volume of compressed gas in the accumulator, respectively, and γ is the adiabatic index for an ideal gas.

2.4. Hydraulic Motor

As pressure increases in the cylinder chamber, oil flows from the cylinder to the motor (through the check-valves and HP accumulator, in the case of the constant-pressure configuration), where pressure and flow discharged from the cylinder are converted into mechanical torque and rotation of the shaft. Once the motor extracts the energy from the high-pressure fluid, low-pressure oil flows back to the cylinder. The torque generated in the hydraulic motor ( T M ) and the output flow ( Q M ) can be expressed, respectively, as,
T M = α D ω Δ p M - T l o s s e s ,
Q M = α D ω ω M - Q l o s s e s ,
where α is the motor displacement fraction, D ω the displacement of the hydraulic motor, ω M the rotational speed of the shaft, and Δ p M the pressure difference between the inlet and outlet ports of the hydraulic motor.

Motor Losses

Modelling mechanical ( T l o s s e s ) and volumetric losses ( Q l o s s e s ) in the motor, due to friction and leakages, is essential to accurately reproduce the performance of hydraulic motors. Different loss models have been suggested to approximate such losses, as described in [31]. The first loss model was suggested by Wilson [32], based on fixed-displacement gear pumps and motors, where the most widely used model structure separating torque and volumetric efficiencies was introduced. Wilson’s model [32] uses constant loss coefficients, and is implemented in a wave energy context in [18]. McCandish and Dorey [33] presented an analytical model, described in detail in [34], to compute losses using variable loss coefficients for fixed- and variable-displacement units. A polynomial based approach is suggested in [35], where the coefficients of the polynomial are identified using flow and torque data from experimental data.
The efficiency model implemented in the mathematical model presented in this paper is, however, the Schösser model [36,37], which includes losses due to friction and leakages in the motor. The beauty of this model is that it can be fitted using data from manufacturers, so no experimental data is required:
Q l o s s e s = Δ p M C Q 1 ,
T l o s s e s = C T 1 + C T 2 Δ p M + C T 3 ω M + C T 4 ω M 2 .
Equations (11) and (12) describe the Schlösser model. Model parameters ( C Q 1 , C T 1 , C T 2 , C T 3 and C T 4 ) are identified, in this case, using the data from manufacturers for the Sauer-Danfoss 250 cc Series 51-1 bent-axis motor (Ames, IA, USA) [38], which includes efficiency data for two specific cases, i.e., full-displacement ( α = 1 ) and partial-displacement ( α = 0 . 3 ). Figure 6a illustrates the efficiency map provided by the manufacturers and Figure 6b shows the efficiency map obtained using the fitted model, both for the full-displacement condition.
The linear least square fitting approach is used to identify the parameters of the models for each case, and once all the parameters for both cases are identified, coefficients can be expressed as a function of the displacement fraction ( α ) using linear interpolation between the two coefficients, as shown in Figure 7 for the coefficient C Q 1 . Linear interpolation is convenient to obtain the parameters of the loss model at different operating points, since manufacturers’ only provide the data for two operating points, but it should be noted that the real behaviour of the motor is not expected to be linear.
Since the manufacturers’ data is given in speed percentage, as shown in Figure 6a, the same model can be employed for different motor sizes by scaling the model as follows [21]:
Q M R S = D ω R S D ω ω r a t e d R S ω r e a t e d Q M ω r a t e d ω r e a t e d R S ω M , Δ p M ,
T M R S = D ω R S D ω T M ω r a t e d ω r a t e d R S ω M , Δ p M ,
where D ω R S is the displacement of the re-sized motor, and ω r a t e d and ω r a t e d R S are the rated rotational speed of the baseline and re-sized motors, respectively.
Several types of hydraulic motors are currently available for diverse applications. Hydraulic motors that can rapidly vary the displacement are needed in WECs, due to the variability of the resource, so high efficiency is required not only at the optimum, but also at part-load operating conditions. Different hydraulic motor topologies are suggested in different hydraulic system configurations, e.g., a bent-axis hydraulic motor in [22] or a swash-plate motor in [21], for which the model structure presented in Equations (10)–(14) is suitable. However, similarly to the case of the hydraulic cylinder, it may be necessary to identify the model parameters for each specific case.
In the mathematical model presented in this paper, pressure losses in transmission lines are neglected, which is reasonable if short transmission lines are assumed.

3. Experimental Set-Up

In order to validate constant pressure and variable pressure circuits, it is necessary to use two separate test rigs. The particulars of these test rigs are discussed in this section, but there are a number of commonalities that it is prudent to document first.
Both test rigs were set-up on the same bed plate at the University of Bath’s Centre for Power Transmission and Motion Control (PTMC). Both rigs were driven by a double ended high precision linear actuator with 250 mm of stroke and piston area of 25 cm 2 . The linear actuator was operated in position control and was powered by a power pack capable of delivering 52 L/min at pressures up to 195 bar. A Simulink Real-Time system was used for data logging and control purposes in both test rigs, using a National Instruments PCI6221 data acquisition (DAQ) board sampling (Newbury, UK) at 1 kHz and a National Instruments PCI6229 DAQ board sampling at 100 Hz on the constant- and variable-pressure rigs, respectively.
Both test rigs were also constructed from discrete components and, as a result, losses within the circuits are higher than would be expected within custom manufactured units of similar design. Maximising efficiency was not the intent of these test rigs, and so these losses were not of concern, as long as they could be quantified for the purpose of model verification.

3.1. Constant-Pressure Configuration

The hydraulic circuit used for validating the constant pressure model varies from that shown in Figure 2 due to its use of an unequal area actuator. This necessitates further changes in order to maintain similar flow rates in extension and retraction with an area ratio of approximately 2:1. Figure 8 shows the modified PTO and Table 1 details the individual components.
Hence, when the actuator is retracting, part of the fluid exiting the piston chamber (chamber A in Figure 8) is transferred, via the regenerative check valve ②, to the annulus chamber (chamber B in Figure 8). As a consequence, the rectified flow that goes to the accumulator 4 and, eventually, to the hydraulic motor 6, is equal to the annular area multiplied by the piston velocity (or half the piston area multiplied by piston velocity, due to the 2:1 area ratio). In contrast, when extending, the annulus flow is rectified to feed the motor and the LP accumulator and boost pump provide make-up flow to the piston chamber. A fixed-displacement motor was used in the PTO along with a DC generator, and a throttle valve was included in the circuit to mimic the variations of the hydraulic motor displacement.
In Figure 8, the location of the six pressure sensors and two flow sensors used in the experimental rig can be seen. The pressure sensors themselves were all Transinstrument 2000 series (Gems Sensors and Controls Ltd, Basingstoke, UK) rated to 250 bar and individually calibrated before installation. The flow meters were from the Hydac EVS3100 series (Sulzbach/Saar, Germany) capable of measuring flows between 6 L/min and 60 L/min and were delivered pre-calibrated and with signal conditioning. The physical location of these flow meters and the other major components are shown in Figure 9 .

3.2. Variable-Pressure Configuration

The hydraulic circuit used to validate the variable pressure configuration also used a fixed displacement motor, rather than the variable displacement motor shown in Figure 3. The leakage from the motor is directly returned to the circuit via an accumulator and a pair of check valves, obviating the need for a boost pump. Pressure signals are measured by means of pressure sensors from the Parker ASIC series (Cleveland, OH, USA) rated to 250 bar, which are calibrated at manufacture. Figure 10 shows the complete hydraulic circuit and Table 2 shows the details of each component.
There are two motors shown in the Figure 11, as different sized motors were tested in the PTO. However, all results discussed in the present paper are measured with the 7.8 cc/rev motor installed.

3.3. Hydraulic Motor Model in AMESim

The validation of the hydraulic system was carried out using experimental results obtained in previous projects, for which the results of the hydraulic motor were not available. Therefore, resort to a high-fidelity hydraulic model was made to validate the mathematical model of the hydraulic motor, as mentioned in Section 1. A comparison is made between the model presented in Section 2.4, referred to as the Schlösser model in the following, and two different models created with the high-fidelity AMESim software, defined as idealised and detailed models.
AMESim is a physical modelling and simulation package for mechatronic systems based on bond graph methods. It is used widely within the fluid power industry and there are numerous publications validating its results against empirical data, including [39,40,41,42].
Both models created in AMESim reproduce the behaviour of a gear motor and use the pressure at the HP and LP accumulators, and the motor’s velocity of the numerical models, as inputs, and output the torque at the motor shaft and flow at the outlet of the hydraulic motor.

3.3.1. Idealised Model

The idealised motor assumes a perfect motor with no losses and continuous flow. Therefore, motor flow in the idealised model is a linear function of the shaft velocity used as input. It is further assumed that there were no losses between the pressure sensors and the motor. This is the simplest model of a hydraulic motor and is created as a reference for future models.

3.3.2. Detailed Model

A second model is also created, which models the individual teeth of the gear motor, and so has a displacement that varies with the absolute rotational position of the motor. The pumping gears have a total of 12 driving and 12 driven teeth with a module of 2.65 mm and width of 7 mm each. These dimensions input to the gear pump model and a contact angle of 20 assumed from the available literature.
This detailed motor model is coupled to a Wilson loss model, which uses the same efficiency map as shown in Figure 6, originally used to fit the hydraulic motor model developed in Section 2.4. The motor torque measured ignores the effect of inertia within the motor or generator, and so is more variable than would actually be seen at the motor shaft.

4. Validation

The two experimental set-ups described in Section 3 allow for a complete validation of the mathematical models, analysing the performance of the different components of the hydraulic PTO systems under different operating conditions.
The inputs for the mathematical models are, for both hydraulic PTO configurations, piston position and velocity, and rotational speed of the electrical generator, as illustrated in Figure 1. Those inputs are taken directly from the experimental data and the outputs from the mathematical models are compared to the test-rig measurements and results from the high-fidelity software AMESim. Figure 12a shows the displacement and velocity of the piston and Figure 12b illustrates the rotational speed for one of the test cases simulated in the constant-pressure experimental set-up.
In the test-rig with the constant-pressure configuration, the experimental PTO was coupled to a hydrodynamic model that reproduces the behaviour of a WEC by means of a hardware-in-the-loop (HIL) system. In the experiments used for the validation in this paper, the WEC was a vertical cylinder of 2 m radius and 6 m draft, with an extended hemisphere in its lower end. More details about the WEC and the HIL system are given in [30]. Hence, the hydrodynamic model computes the position and velocity of the WEC (and, thus, postion and velocity of the piston in the hydraulic cylinder), which are used as inputs for the experimental PTO.
Three test cases, based on regular waves, were studied in the constant-pressure test rig:
  • 10 s period and 1.5 m wave amplitude,
  • 10 s period and 1 m wave amplitude, and
  • 12 s period and 1.25 m wave amplitude,
which needed to be scaled down (using a scaling factor of 10) to be implemented in the experimental PTO. For the validation of the mathematical model with the constant-pressure configuration, the hydrodynamic model of the HIL system is ignored, since the mathematical model for the hydraulic PTO is isolated, and piston position and velocity are used as inputs.
In the case of the variable-pressure configuration, no HIL system was implemented, and piston position and velocity were used as inputs. A total of 20 test cases were analysed, combining four different sinusoidal input signals and five damping coefficients for the force control, implemented through the hydraulic motor.
Hydraulic cylinders and motors are the most important components of the hydraulic systems because the performance of the whole system essentially depends on the performance of these two components. The major challenge when modelling hydraulic cylinders is to accurately reproduce losses and dynamics due to friction. In the case of the hydraulic motor, friction and leakage losses are the main challenge.
Therefore, the validation of the mathematical models focuses on the two main components, carefully examining friction effects in hydraulic cylinders, and losses in hydraulic motors.

4.1. Hydraulic Cylinder

The performance of the hydraulic cylinder strongly depends on the system topology. In the case of constant-pressure configurations, the pressure in the cylinder chambers and, as a consequence, the force applied on the WEC, are practically constant, as illustrated in Figure 13 and Figure 14. In contrast, the pressure and force of the hydraulic cylinder in variable-pressure configurations follow the profile of the velocity, as shown in Figure 15 and Figure 16.

4.1.1. Constant-Pressure Configuration

Pressure in the hydraulic cylinder is driven by the pressure in the hydraulic accumulators, which, in turn, depends on the pre-charge pressure and the oil volume, as shown in Equation (8). The pressure of the oil in the accumulator sets the threshold for the check vales, so that check valves open only when pressure in the cylinder chamber is higher than in the HP accumulator.
As one can note from Figure 13, the hydraulic cylinder in the constant-pressure configuration repeats the same cycle during the whole simulation. When the piston reaches its final position (1 in Figure 13), the pressure in the high-pressure chamber starts to reduce, falling below the accumulator pressure. At that point, check valves close and the piston starts to oscillate for a short period of time, which produces pressure fluctuations in both chambers (2 in Figure 13). As soon as the pressure in any of the chambers increases again over the pressure in the HP accumulator, check valves open again, and the piston starts to move ‘freely’, pushing the oil in the chamber, and, consequently, increasing the pressure in the chamber (3 in Figure 13). When the piston again reaches its end point, the whole process repeats in the other chamber.
The behaviour of the hydraulic cylinder in the constant-pressure configuration is similar for all the test cases analysed for the validation, where the profile of the pressure signals is practically identical and only the magnitude of the signals changes from one test case to another. Figure 13 illustrates the pressure in the chambers of the hydraulic cylinder for the 12 s period and 1.25 m amplitude test case, showing good agreement between the experimental results and the results from the mathematical model, in the sense that the profile and magnitude of the pressure signals are similar. However, the pressure signals appear to be slightly flatter at low pressures for the experimental results, while one can notice a negative overshoot for the results obtained from the mathematical model. In addition, the pressure drop across the check valves is shown to be larger for the experimental case and a phase lag exists between the experimental measurement and the numerical simulation as a result of the different pressure transitions, which are instantaneous in the mathematical model.
Similar agreement can be observed in Figure 14 for the piston force, where friction or inertia forces are also considered. The phase lag between the mathematical model and the experiments remains, as expected, in the force signals. However, it should be noted that the effect of friction is rather low in these test cases (about 1% of power loss due to friction based on the mathematical model).

4.1.2. Variable-Pressure Configuration

In the case of the variable pressure configuration, agreement between the experiments and the numerical simulations is also good for the different test cases, as illustrated in Figure 15, where the pressure signal is clearly asymmetric. This asymmetry appears due to a non-symmetrical control force signal, with a different damping coefficient for each direction of motion.
Similarly to the constant-pressure configuration, measured and simulated force signals match accurately, as illustrated in Figure 16. Apart from the total force signal, Figure 16 also shows the contributions of the different factors, such as pressure difference, friction or piston mass.
Hence, one can observe that, for the pressure force, referred to as the force considering only the contribution of the pressure difference in the cylinder ( F Δ P m o d e l ), the amplitude is lower than the experimental piston force ( F p i s d a t a ). In addition, pressure force is, as expected, asymmetric, while experimental piston force appears to be completely symmetric. However, when friction force and gravity force due to the piston mass are added to the pressure force, the total piston force obtained from the mathematical model ( F p i s m o d e l ) accurately matches the experimental piston force.
Finally, the friction force signal shows the relevance of the Coulomb and static friction over the viscous friction and, as a consequence, the impact of the friction force is especially recognisable at low velocity. When friction force suddenly changes from positive to negative, or vice versa, force signal also shows an abrupt change, which can only be reproduced by accurately modelling friction. Figure 17 illustrates the impact of the friction force at low velocity, where the force signal considering the pressure difference alone simply follows the velocity profile, while the signal that includes friction clearly shows an abrupt drop, similarly to the measured force signal.

4.2. Hydraulic Motor

Losses in the hydraulic motor also depend on the characteristics of the operation, especially pressure difference and rotational speed, as described in the Schlösser model defined in Equations (11) and (12).
The well-known simulation software, AMESim, is employed to accurately model the performance of the hydraulic motor, as described in Section 3.3. Hence, pressure results from the numerical model described in Section 2 are used as inputs for the hydraulic motor model created in AMESim to validate the Schlösser model developed in Section 2.4.
Thus, flow and torque in the hydraulic motor are analysed, comparing the results from the mathematical model described in Section 2.4.1, referred to as the Schlösser model in the following paragraphs, to the results obtained from the AMESim models.
Figure 18a,b illustrate the motor flow and torque, respectively, where the flow of the Schlösser model stays always below the ideal flow, due to the losses neglected in the ideal AMESim model, and the torque stays always above, to overcome the friction losses and supply the required torque. Compared to the detailed AMESim model, the flow and torque of the Schlösser model appear to follow the upper envelope of the flow and torque, which suggests that the efficiency is similar for the detailed AMESim and Schlösser models. However, significantly more variations can be observed in the detailed AMESim model, due to the modelling of the movement of the individual pumping elements, which cannot be captured unless a detailed model is used.

5. Conclusions

The present paper validates the mathematical models for hydraulic transmission systems used in the power take-off of wave energy converters. The different topologies suggested in the literature are divided into two main configurations, depending on the operational characteristics, as described in Section 1.1: constant-pressure and variable-pressure configurations. The mathematical equations to model the different components used in the two configurations are given in Section 2, including the main nonlinear dynamics and losses, such as friction in hydraulic cylinders or losses in the hydraulic motor.
Two experimental set-ups described in Section 3, including a set-up that reproduces the characteristics of a constant-pressure configuration and another set-up that emulates the variable-pressure configuration, are employed in the validation. That way, the validation covers the different possibilities of using hydraulic transmission systems in wave energy. In addition, a well-known software for hydraulic systems, AMESim, is used to validate the performance of the hydraulic motor. Results from the mathematical model and the experiments are compared in Section 4, which mainly focus on the performance of the hydraulic cylinder and motor.
Results of the hydraulic cylinder show very good agreement for pressure and force values in the different test cases analysed for the two configurations. In the constant-pressure configuration, the effect of the rectification valves is well reproduced by the mathematical model, including the pressure and force oscillations that arise when the piston approaches its final position. In addition, the mathematical model shows the ability to capture the asymmetry observed in the variable-pressure configuration experiments and demonstrates the effectiveness of the friction model, reproducing the effect of the static friction at low velocities.
The Schlösser model for the hydraulic motor is shown to be more realistic than the ideal model implemented in AMESim, while the efficiency of the Schlösser model appears to be similar to the detailed AMESim model, although the detailed AMESim model includes other effects, such as wave propagation that cause larger variations in the flow and torque.

Acknowledgments

This material is based upon works supported by the Science Foundation Ireland under Grant No. 13/IA/1886.

Author Contributions

Markel Penalba and John V. Ringwood participated in designing the mathematical models; Nathan P. Sell and Andrew J. Hillis provided the experimental data; Markel Penalba and Nathan P. Sell wrote the manuscript; and John V. Ringwood revised it critically for important intellectual content. Finally, all the authors gave the final approval of the version to be submitted.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

( ) A , ( ) B Subscripts that refer to cylinder chambers chambers A and B
α Hydraulic motor displacement fraction
β e f f Effective bulk modulus of the fluid
Δ p M Pressure difference between the inlet and outlet port in the hydraulic motor
Δ p v Pressure difference accross the valve
Δ p c r a c k Cracking pressure of the valve
γ Adiabatic index for an idean gas
ω M Rotational speed of the hydraulic motor shaft
ω r a t e d Rated rotational speed of the baseline hydraulic motor
ω r a t e d R S Rated rotational speed of the re-sized hydraulic motor
ρ o i l Density of the hydraulic oil
σ Viscous coefficient of the Stribeck friction model
A p Piston area
A v The valve opening area
C d Discharge coefficient
C Q 1 , C T 1 , C T 2 , C T 3 and C T 4 Parameters of the Schlösser loss model
c s t Characteristic velocity of the Stribeck curve
D ω R S Displacement of the re-sized hydraulic motor
D ω Displacement of the hydrualic motor
F c Coulomb friction force of the Stribeck friction model
F f r i c Friction force in the hydraulic cylinder
F i n Inertia force in the hydraulic cylinder
F P T O Power take-off force applied on the WEC
F s t Static friction force of the Stribeck friction model
gAcceleration due to the gravity
M p Cylinder piston mass
M r Cylinder rod mass
M o i l Mass of the oil in the cylinder chambers
pPressure in the cylinder chamber
p a c c Pressure in the hydrualic accumulator
QFlow entering or exiting the cylinder chamber, when leakages are neglected
Q M Output flow of the hydraulic motor
Q M R S Output flow of the re-sized hydraulic motor
Q v Flow through the valve
Q a c c i n Input flow to the accumulator
Q l o s s Volumetric losses in the hydraulic motor
T M Torque generated in the hydraulic motor
T M R S Torque generated in the re-seized hydrualic motor
T l o s s Mechanical losses in the hydraulic motor
VMinimum volume in the cylinder chamber
V a c c Volume of hydraulic fluid contained in the accumulator
V p r e Pre-charge pressure in the hydraulic accumulator
V t o t Total volume of the hydraulic accumulator
V g a s Volume of compressed gas in the accumulator
x, x ˙ and x ¨ piston position, velocity and acceleration

References

  1. Penalba, M.; Ringwood, J.V. A Review of Wave-to-Wire Models for Wave Energy Converters. Energies 2016, 9, 506. [Google Scholar] [CrossRef]
  2. Bailey, H.; Robertson, B.R.; Buckham, B.J. Wave-to-wire simulation of a floating oscillating water column wave energy converter. Ocean Eng. 2016, 125, 248–260. [Google Scholar] [CrossRef]
  3. Garcia-Rosa, P.B.; Vilela Soares Cunha, J.P.; Lizarralde, F.; Estefen, S.F.; Machado, I.R.; Watanabe, E.H. Wave-to-Wire Model and Energy Storage Analysis of an Ocean Wave Energy Hyperbaric Converter. IEEE J. Ocean. Eng. 2014, 39, 386–397. [Google Scholar] [CrossRef]
  4. Tedeschi, E.; Carraro, M.; Molinas, M.; Mattavelli, P. Effect of Control Strategies and Power Take-Off Efficiency on the Power Capture From Sea Waves. IEEE Trans. Energy Convers. 2011, 26, 1088–1098. [Google Scholar] [CrossRef]
  5. Amundarain, M.; Alberdi, M.; Garrido, A.J.; Garrido, I. Modeling and Simulation of Wave Energy Generation Plants: Output Power Control. IEEE Trans. Ind. Electron. 2011, 58, 105–117. [Google Scholar] [CrossRef]
  6. Henderson, R. Design, simulation, and testing of a novel hydraulic power take-off system for the Pelamis wave energy converter. Renew. Energy 2006, 31, 271–283. [Google Scholar] [CrossRef]
  7. Josset, C.; Babarit, A.; Clément, A.H. A wave-to-wire model of the SEAREV wave energy converter. Proc. Inst. Mech. Eng. Part M J. Eng. Marit. Environ. 2007, 221, 81–93. [Google Scholar] [CrossRef]
  8. Hansen, R.H.; Andersen, T.O.; Pedersen, H.C.; Hansen, A.H. Control of a 420 KN Discrete Displacement Cylinder Drive for the Wavestar Wave Energy Converter. In Proceedings of the ASME/BATH 2014 Symposium on Fluid Power and Motion Control, Bath, UK, 10–12 September 2014; p. V001T01A021. [Google Scholar]
  9. Henry, A.; Doherty, K.; Cameron, L.; Whittaker, T.; Doherty, R. Advances in the design of the Oyster wave energy converter. In Proceedings of the RINA Marine and Offshore Energy Conference, London, UK, 21–23 April 2010. [Google Scholar]
  10. Fiévez, J.; Sawyer, T. Lessons Learned from Building and Operating a Grid Connected Wave Energy Plant. In Proceedings of the 11th European Wave and Tidal Energy Conference, Nantes, France, 6–11 September 2015. [Google Scholar]
  11. Lucas, J.; Livingstone, M.; Vuorinen, M.; Cruz, J. Development of a wave energy converter (WEC) design tool—Application to the WaveRoller WEC including validation of numerical estimates. In Proceedings of the 4th International Conference on Ocean Energy, Dublin, Ireland, 17–19 October 2012. [Google Scholar]
  12. Gear, C.W.; Wells, D. Multirate linear multistep methods. BIT Numer. Math. 1984, 24, 484–502. [Google Scholar] [CrossRef]
  13. Engstler, C.; Lubich, C. Multirate extrapolation methods for differential equations with different time scales. Computing 1997, 58, 173–185. [Google Scholar] [CrossRef]
  14. Giorgi, G.; Penalba, M.; Ringwood, J.V. Nonlinear Hydrodynamic Models for Heaving Buoy Wave Energy Converters. In Proceedings of the 3rd Asian Wave and Tidal Energy Conference, Singapore, 25–27 October 2016; pp. 144–153. [Google Scholar]
  15. Penalba, M.; Cortajarena, J.; Ringwood, J.V. Validating a Wave-to-Wire Model for a Wave Energy Converter—Part II: The electrical system. Energies. submitted for publication.
  16. Ricci, P.; Lopez, J.; Santos, M.; Ruiz-Minguela, P.; Villate, J.L.; Salcedo, F.; Falcão, A.F.O. Control strategies for a wave energy converter connected to a hydraulic power take-off. IET Renew. Power Gener. 2011, 5, 234–244. [Google Scholar] [CrossRef]
  17. Cargo, C.J.; Plummer, A.R.; Hillis, A.J.; Schlotter, M. Determination of optimal parameters for a hydraulic power take-off unit of a wave energy converter in regular waves. Proc. Inst. Mech. Eng. Part A J. Power Energy 2012, 226, 98–111. [Google Scholar] [CrossRef]
  18. Cargo, C.J.; Hillis, A.J.; Plummer, A.R. Optimisation and control of a hydraulic power take-off unit for a wave energy converter in irregular waves. Proc. Inst. Mech. Eng. Part A J. Power Energy 2014, 228, 462–479. [Google Scholar] [CrossRef]
  19. Forehand, D.I.; Kiprakis, A.E.; Nambiar, A.J.; Wallace, A.R. A Fully Coupled Wave-to-Wire Model of an Array of Wave Energy Converters. IEEE Trans. Sustain. Energy 2016, 7, 118–128. [Google Scholar] [CrossRef]
  20. Kamizuru, Y. Development of Hydrostatic Drive Trains for Wave Energy Converters. Ph.D. Thesis, RWTH Aachen University, Aachen, Germany, 2014. [Google Scholar]
  21. Hansen, R.H.; Andersen, T.O.; Pedersen, H.C. Model based design of efficient power take-off systems for wave energy converters. In Proceedings of the 12th Scandinavian International Conference on Fluid Power, Tampere, Finland, 18–20 May 2011. [Google Scholar]
  22. Hansen, R.H.; Kramer, M.M.; Vidal, E. Discrete displacement hydraulic power take-off system for the wavestar wave energy converter. Energies 2013, 6, 4001–4044. [Google Scholar] [CrossRef]
  23. Plummer, A.; Cargo, C.C. Power Transmissions for Wave Energy Converters: A Review. In Proceedings of the 8th International Fluid Power Conference (IFK), Dresden, Germany, 26–28 March 2012. [Google Scholar]
  24. Antonio, F.d.O. Modelling and control of oscillating-body wave energy converters with hydraulic power take-off and gas accumulator. Ocean Eng. 2007, 34, 2021–2032. [Google Scholar]
  25. Hansen, R.H. Design and Control of the PowerTake-Off System for a Wave Energy Converter with Multiple Absorbers. Ph.D. Thesis, Aalborg Universitet, Aalborg, Denmark, 2013. [Google Scholar]
  26. AMESim. Available online: http://www.plm.automation.siemens.com/engb/products/lms/imagine-lab/amesim/platform/index.shtml (accessed on 3 April 2017).
  27. Olsson, H.; Åström, K.J.; De Wit, C.C.; Gäfvert, M.; Lischinsky, P. Friction models and friction compensation. Eur. J. Control 1998, 4, 176–195. [Google Scholar] [CrossRef]
  28. Jelali, M.; Kroll, A. Hydraulic Servo-Systems: Modelling, Identification and Control; Springer Science & Business Media: Berlin, Germany, 2012. [Google Scholar]
  29. Fisher, Emerson Process Management. Control Valve Handbook; Fisher Controls International LLC: Marshalltown, IA, USA, 2005. [Google Scholar]
  30. Cargo, C. Design and Control of Hydraulic Power Take-Off for Wave Energy Converters. Ph.D. Thesis, University of Bath, Bath, UK, 2012. [Google Scholar]
  31. Grandall, D.R. The Performance and Efficiency of Hydraulic Pumps and Motors. Master’s Thesis, University of Minnesota, Minneapolis, MN, USA, 2010. [Google Scholar]
  32. Wilson, W. Rotary-pump theory. Trans. ASME 1946, 68, 371–384. [Google Scholar]
  33. McCandlish, D.; Dorey, R. The mathematical modeling of hydrostatic pumps and motors. Proc. Inst. Mech. Eng. 1984, 198B, 165–174. [Google Scholar] [CrossRef]
  34. Dorey, R. Modeling of losses in pumps and motors. In Proceedings of the First Bath International Fluid Power Workshop, Bath, UK, 8 September 1988; pp. 71–97. [Google Scholar]
  35. Mikeska, D. A precise steady-state model of displacement machines for the application in virtual prototyping of power-split drives. In Proceedings of the 2nd International FPNI-PhD Symposium on Fluid Power, Modena, Italy, 3–6 July 2002. [Google Scholar]
  36. Schlöesser, W. Mathematical model for displacement pumps and motors. Hydraul. Power Trans. 1961, 252–257. [Google Scholar]
  37. Schlösser, W. The overall efficiency of positive-displacement pumps. In Proceedings of the 1st BHRA Fluid Power Symposium, Cranfield, UK, January 1969; pp. 34–48. [Google Scholar]
  38. Sauer-Danfoss. Series 51 and 51-1 Bent Axis Variable Displacement Motors Technical Information. Available online: https://www.bibus.hu/fileadmin/editors/countries/bihun/productdata/sauer-danfoss/documents/sauerdanfossseriess51motorscatalogueen520l0440.pdf (accessed on 24 May 2017).
  39. Borghi, M.; Milani, M.; Piraccini, M.; Flandoli, M. Analysis of a F1 Car Hydraulic Circuit Dynamic Behavio; SAE Technical Report; SAE: Warrendale, PA, USA, 2001. [Google Scholar]
  40. TIKKANEN, S. Influence of line design on pump performance. In Power Transmission and Motion Control: PTMC 2001; John Wiley & Sons: Hoboken, NJ, USA, 2001; p. 33. [Google Scholar]
  41. Manco, S.; Nervegna, N.; Rundo, M. Effects of Timing and Odd/even Number of Teeth on Noise Generation of Gerotor Lubricating Pumps for IC Engines; SAE Technical Report; SAE: Warrendale, PA, USA, 2000. [Google Scholar]
  42. Nervegna, N.; Manco, S.; Rundo, M. Variable Flow Internal Gear Pump. In Proceedings of the ASME International Mechanical Engineering Congress and Exposition, New York, NY, USA, 11–16 November 2001; American Society of Mechanical Engineering: New York, NY, USA, 2001. [Google Scholar]
Figure 1. Full conversion chain of a wave energy converter with hydraulic PTO, adapted from [1].
Figure 1. Full conversion chain of a wave energy converter with hydraulic PTO, adapted from [1].
Energies 10 00977 g001
Figure 2. The scheme of a typical constant-pressure hydraulic transmission system.
Figure 2. The scheme of a typical constant-pressure hydraulic transmission system.
Energies 10 00977 g002
Figure 3. The scheme of a typical variable-pressure hydraulic transmission system.
Figure 3. The scheme of a typical variable-pressure hydraulic transmission system.
Energies 10 00977 g003
Figure 4. Friction losses in the hydraulic cylinder: (a) the impact of different friction effects normalised against the maximum friction force; and (b) the approximated efficiency map of a cylinder, reproduced from [20].
Figure 4. Friction losses in the hydraulic cylinder: (a) the impact of different friction effects normalised against the maximum friction force; and (b) the approximated efficiency map of a cylinder, reproduced from [20].
Energies 10 00977 g004
Figure 5. Valve opening functions.
Figure 5. Valve opening functions.
Energies 10 00977 g005
Figure 6. Efficiency map of the Sauer-Danfoss 250 cc Series 51-1 bent-axis motor, (a) measures provided by [38] and (b) the reproduction of the efficiency map using the model.
Figure 6. Efficiency map of the Sauer-Danfoss 250 cc Series 51-1 bent-axis motor, (a) measures provided by [38] and (b) the reproduction of the efficiency map using the model.
Energies 10 00977 g006
Figure 7. Fitted parameter ( C Q 1 ) for the two displacement fractions available in [38] and the linear curve that fits the two parameters.
Figure 7. Fitted parameter ( C Q 1 ) for the two displacement fractions available in [38] and the linear curve that fits the two parameters.
Energies 10 00977 g007
Figure 8. Experimental constant-pressure circuit.
Figure 8. Experimental constant-pressure circuit.
Energies 10 00977 g008
Figure 9. Front and side view of constant pressure test rig.
Figure 9. Front and side view of constant pressure test rig.
Energies 10 00977 g009
Figure 10. Experimental variable-pressure circuit.
Figure 10. Experimental variable-pressure circuit.
Energies 10 00977 g010
Figure 11. Variable pressure test rig.
Figure 11. Variable pressure test rig.
Energies 10 00977 g011
Figure 12. Inputs for the mathematical model: (a) position and velocity of the piston; and (b) rotational speed of the electrical generator.
Figure 12. Inputs for the mathematical model: (a) position and velocity of the piston; and (b) rotational speed of the electrical generator.
Energies 10 00977 g012
Figure 13. Pressure in the chambers A and B of the hydraulic cylinder from the experiments and the mathematical model for the constant-pressure configuration.
Figure 13. Pressure in the chambers A and B of the hydraulic cylinder from the experiments and the mathematical model for the constant-pressure configuration.
Energies 10 00977 g013
Figure 14. Pressure in the chambers A and B of the hydraulic cylinder from the experiments and the mathematical model for the constant-pressure configuration.
Figure 14. Pressure in the chambers A and B of the hydraulic cylinder from the experiments and the mathematical model for the constant-pressure configuration.
Energies 10 00977 g014
Figure 15. Piston force from the experiments and the mathematical model for the constant-pressure configuration.
Figure 15. Piston force from the experiments and the mathematical model for the constant-pressure configuration.
Energies 10 00977 g015
Figure 16. Contribution of the different forces in the hydraulic cylinder compared to the piston force measured in the test rig with the variable-pressure configuration.
Figure 16. Contribution of the different forces in the hydraulic cylinder compared to the piston force measured in the test rig with the variable-pressure configuration.
Energies 10 00977 g016
Figure 17. The effect of friction force at low velocity, comparing measured piston force and the contribution of the different effects obtained from the numerical model for the variable-pressure configuration.
Figure 17. The effect of friction force at low velocity, comparing measured piston force and the contribution of the different effects obtained from the numerical model for the variable-pressure configuration.
Energies 10 00977 g017
Figure 18. Flow (a) and torque (b) in the hydraulic motor for the mathematical model presented in Section 2.4.1, and the detailed and ideal models implemented in AMESim.
Figure 18. Flow (a) and torque (b) in the hydraulic motor for the mathematical model presented in Section 2.4.1, and the detailed and ideal models implemented in AMESim.
Energies 10 00977 g018
Table 1. Constant pressure test rig components.
Table 1. Constant pressure test rig components.
No.DescriptionDetails
1PTO cylinder40 mm bore, 28 mm rod, 300 mm stroke
2Regenerative check valves0.35 bar preload
3Rectification check valves0.35 bar preload
4HP accumulator3.8 L, 10 bar precharge
5Relief valve100 bar cracking pressure
6Gear motor4.0 cc/rev
7DC Generator90 W at 3000 rpm
8LP accumulator1 L, 2 bar precharge
9Boost pump4 cc/rev, 7 bar relief pressure
10Throttle valve3/8" needle valve
Table 2. Variable pressure test rig components.
Table 2. Variable pressure test rig components.
No.DescriptionDetails
1PTO cylinder30 mm bore, 25 mm rod and 500 mm stroke
2Relief valves200 bar cracking pressure
3Gear motor7.8 cc/rev displacement
4Leakage check0.35 bar cracking pressure
5Accumulator5 L, 1 bar precharge

Share and Cite

MDPI and ACS Style

Penalba, M.; Sell, N.P.; Hillis, A.J.; Ringwood, J.V. Validating a Wave-to-Wire Model for a Wave Energy Converter—Part I: The Hydraulic Transmission System. Energies 2017, 10, 977. https://doi.org/10.3390/en10070977

AMA Style

Penalba M, Sell NP, Hillis AJ, Ringwood JV. Validating a Wave-to-Wire Model for a Wave Energy Converter—Part I: The Hydraulic Transmission System. Energies. 2017; 10(7):977. https://doi.org/10.3390/en10070977

Chicago/Turabian Style

Penalba, Markel, Nathan P. Sell, Andy J. Hillis, and John V. Ringwood. 2017. "Validating a Wave-to-Wire Model for a Wave Energy Converter—Part I: The Hydraulic Transmission System" Energies 10, no. 7: 977. https://doi.org/10.3390/en10070977

APA Style

Penalba, M., Sell, N. P., Hillis, A. J., & Ringwood, J. V. (2017). Validating a Wave-to-Wire Model for a Wave Energy Converter—Part I: The Hydraulic Transmission System. Energies, 10(7), 977. https://doi.org/10.3390/en10070977

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop