Research article 12 Sep 2018
Research article  12 Sep 2018
Control design, implementation, and evaluation for an infield 500 kW wind turbine with a fixeddisplacement hydraulic drivetrain
 ^{1}Delft Center for Systems and Control, Faculty of Mechanical Engineering, Delft University of Technology, Mekelweg 2, 2628 CD Delft, the Netherlands
 ^{2}DOT B.V., Raam 180, 2611 WP Delft, the Netherlands
 ^{1}Delft Center for Systems and Control, Faculty of Mechanical Engineering, Delft University of Technology, Mekelweg 2, 2628 CD Delft, the Netherlands
 ^{2}DOT B.V., Raam 180, 2611 WP Delft, the Netherlands
Correspondence: Sebastiaan Paul Mulders (s.p.mulders@tudelft.nl)
Hide author detailsCorrespondence: Sebastiaan Paul Mulders (s.p.mulders@tudelft.nl)
The business case for compact hydraulic wind turbine drivetrains is becoming ever stronger, as offshore wind turbines are getting larger in terms of size and power output. Hydraulic transmissions are generally employed in highload systems and form an opportunity for application in multimegawatt turbines. The Delft Offshore Turbine (DOT) is a hydraulic wind turbine concept replacing conventional drivetrain components with a single seawater pump. Pressurized seawater is directed to a combined Pelton turbine connected to an electrical generator on a central multimegawatt electricity generation platform. This paper presents the control design, implementation, and evaluation for an intermediate version of the ideal DOT concept: an infield 500 kW hydraulic wind turbine. It is shown that the overall drivetrain efficiency and controllability are increased by operating the rotor at maximum rotor torque in the belowrated region using a passive torque control strategy. An active valve control scheme is employed and evaluated in nearrated conditions.
The drivetrain of horizontalaxis wind turbines (HAWTs) generally consists of a rotor–gearbox–generator configuration in the nacelle, which enables each wind turbine to produce and deliver electrical energy independent of other wind turbines. While the HAWT is a proven concept, the turbine rotation speed decreases asymptotically and torque increases exponentially with increasing blade length and power ratings (Burton et al., 2011). As offshore wind turbines are getting ever larger, this results in a lower rotation speed and higher torque at the rotor axis. This inevitably leads to design challenges when scaling up conventional turbine drivetrains for the highload subsystems (Kotzalas and Doll, 2010). The increased loads primarily affect highweight components in the turbine, such as the generator, bearings, and gearbox, and makes repair and replacement a costly and challenging task (Ragheb and Ragheb, 2010; Spinato et al., 2009). Furthermore, due to the contribution of all nacelle components to the total nacelle mass, the complete wind turbine support structure and foundation are designed to carry this weight for the entire expected lifetime, which in turn leads to extra material, weight, and thus total cost of the wind turbine (EWEA, 2009; Fingersh et al., 2006).
In an effort to reduce turbine weight, maintenance requirements, complexity, and thus the levelized cost of energy (LCOE) for offshore wind, hydraulic drivetrain concepts have been considered in the past (Piña Rodriguez, 2012). Integration of hydraulic transmission systems in offshore wind turbines seems to be an interesting opportunity, as they are generally employed in highload applications and have the advantage of a high powertoweight ratio (Merritt, 1967). It is concluded in Silva et al. (2014) that hydrostatic transmissions could lower drivetrain costs, improve system reliability, and reduce the nacelle mass.
Various hydraulic turbine concepts have been considered in the past. The first 3 MW wind turbine with a hydrostatic power transmission was the SWT3, developed and built from 1976 to 1980 by Bendix (Rybak, 1981). After deployment of the turbine with 14 fixeddisplacement oil pumps in the nacelle and 18 variabledisplacement motors at the tower base, it proved to be overly complex, inefficient and unreliable. With the aim of eliminating the power electronics entirely by application of a synchronous generator, Voith introduced the WinDrive in 2003 that provides a hydraulic gearbox with a variable transmission ratio (Muller et al., 2006). In 2004, the ChapDrive hydraulic drivetrain solution was developed with the aim of driving a synchronous generator by a fixeddisplacement oil pump and variabledisplacement oil motor (Chapple et al., 2011; Thomsen et al., 2012). Although the company acquired funding from Statoil for a 5 MW concept, the company ceased operations. In cooperation with Hägglunds, Statoil modeled a drivetrain with hydrostatic transmission, consisting of a single oil pump connected to the rotor and six motors at ground level (Skaare et al., 2011, 2013). Each motor can be enabled or disabled to obtain a discrete variable transmission ratio to drive a synchronous generator. In 2005, Artemis Ltd. developed a digital displacement pump, meaning that it can continuously adjust its volume displacement in a digital way by enabling and disabling individual cylinders (Artemis Intelligent Power, 2018; Rampen, 2006). In 2010, Mitsubishi acquired Artemis Intelligent Power and started testing of its SeaAngel 7 MW turbine with hydraulic power drive technology in 2014 (Sasaki et al., 2014).
The abovedescribed full hydraulic concepts aim to eliminate power electronics from the turbine for the use of a synchronous generator and therefore use a mechanism to vary the hydraulic gear ratio. However, to date, none of the concepts have made their way to a commercial product. All concepts use oil as the hydraulic medium because of the favorable fluid properties and component availability, but therefore also need to operate in a closed circuit. Closedcircuit operation for an offshore wind application using oil is required to minimize the risk of environmental pollution, but also to abandon the need for a continuous fresh oil supply to the circuit. Furthermore, often an additional cooling circuit is needed when losses in hydraulic components are significant and natural heat convection to the surroundings is insufficient.
A novel and patented hydraulic concept with an opencircuit drivetrain using seawater as the hydraulic medium is the Delft Offshore Turbine (DOT) (Van der Tempel, 2009), and it is shown in Fig. 1. The open circuit is enabled by the use of preconditioned seawater and alleviates the need for a cooling circuit by the continuous fresh supply. The DOT concept only requires a single seawater pump directly connected to the turbine rotor. The pump replaces components with high maintenance requirements in the nacelle, which reduces the weight, support structure requirements, and turbine maintenance frequency. This effect is clearly visualized in Fig. 2. All maintenancecritical components are located at sea level, and the centralized generator is coupled to a Pelton turbine. Turbines collectively drive the Pelton turbine to harvest the hydraulic into electrical energy. A feasibility study and modeling of a hydraulic wind turbine based on the DOT concept is performed in Diepeveen (2013). Hydraulic wind turbine networks employing variabledisplacement components are modeled and simulated in Jarquin Laguna (2017). Using the DOT concept, a wind turbine drivetrain to generate electricity and simultaneously extract thermal energy has also been proposed in Buhagiar and Sant (2014).
This paper presents the first steps in realizing the integrated hydraulic wind turbine concept by fullscale prototype tests with a retrofitted Vestas V44 600 kW wind turbine, the conventional drivetrain of which is replaced by a 500 kW hydraulic configuration. A spear valve is used to control the nozzle outlet area, which in effect influences the fluid pressure in the hydraulic discharge line of the water pump and forms an alternative way of controlling the reaction torque to the rotor. Preliminary results of the infield tests were presented earlier in Mulders et al. (2018a).
The main contribution of this paper is to elaborate on mathematical modeling and controller design for a hydraulic drivetrain with fixeddisplacement components, subject to efficiency and controllability maximization of the system. The control design is based on steadystate and dynamic turbine models, which are subsequently evaluated on the actual infield retrofitted 500 kW wind turbine. Furthermore, a framework for the modeling of fluid dynamics is provided, a parameter study on how different design variables influence the controllability is given, and future improvements to the system and control design are proposed.
The organization of this paper is as follows. In Sect. 2, the DOT configuration used during the infield tests is explained, and drivetrain components are specified. Section 3, which involves drivetrain modeling, is divided into two parts: a steadystate drivetrain model is derived in Sect. 3.1, and a drivetrain model including fluid dynamics is presented in Sect. 3.2. Controller design is presented in Sect. 4, and the steadystate model is used in Sect. 4.1 for the design of a passive control strategy for belowrated operation. In Sect. 4.2 the dynamic model is used to derive an active control strategy for the region between below and aboverated operation (nearrated region). Because the infield tests are performed prior to theoretical model derivation and control design, preliminary conclusions from these tests are incorporated. In Sect. 5, infield test results are presented for the evaluation of the overall control design. Finally, in Sect. 6, conclusions are drawn and an outlook for the DOT concept is given.
In this section, the intermediate prototype DOT turbine on which infield tests are performed is described in Sect. 2.1. Subsequently, the drivetrain components used for the intermediate concept are discussed in Sect. 2.2.
2.1 The intermediate DOT500 prototype
At the time of writing, a lowspeed hightorque seawater pump required for the ideal DOT concept is not commercially available. This pump is being developed by DOT, enabling the ideal concept in later stages of the project (Nijssen et al., 2018). An intermediate setup using offtheshelf components is proposed to speed up development and test the practical feasibility. A visualization of the DOT500 setup is given in Fig. 3.
A Vestas V44 600 kW turbine is used and its drivetrain is retrofitted into a 500 kW hydraulic configuration. The original Vestas V44 turbine is equipped with a conventional drivetrain consisting of the main bearing, a gearbox, and a 600 kW threephase asynchronous generator. The blades are pitched collectively by means of a hydraulic cylinder driven by a HPU (hydraulic power unit) with a safety pressure accumulator. The drivetrain of the Vestas turbine is replaced by a hydraulic drivetrain. This means that the gearbox and generator are removed from the nacelle and replaced by a single oil pump coupled to the rotor via the main bearing.
In the retrofitted hydraulic configuration, a lowspeed oil pump is coupled to the rotor and its flow hydraulically drives a highspeed oilmotor–waterpump combination at the bottom of the turbine tower. The oil circuit acts as a hydraulic gearbox between the rotor and the water pump. The water circuit, including the Pelton generator combination, is depicted in Fig. 1 as an ideal setup. It is known and taken into account that the additional components and energy conversions result in a reduced overall efficiency: the infield tests show a belowrated efficiency in the range of 30 %–45 %, whereas in the aboverated region a drivetrain efficiency of 45 % is attained. However, the DOT500 allows for prototyping and provides a proof of concept for faster development towards the ideal DOT concept. From this point onwards, all discussions and calculations will refer to the intermediate DOT500 setup, including the described oil circuit (Diepeveen et al., 2018). A prototype was erected in June 2016 at Rotterdam Maasvlakte II, the Netherlands.
2.2 Drivetrain component specification
The drivetrain component specifications are summarized in Table 1 and a schematic overview and photograph of the DOT500 setup is presented in Fig. 3. An exhaustive description of the component matching process for the DOT500 is described in Diepeveen (2013). The components are described according to the enumerated labels in the figure.

Oil pump. The rotor drives a Hägglunds CB840, which is a fixeddisplacement radial piston motor. The motor is used here as pump and is referred to as the oil pump in the remainder of this paper. The pump is supplied with sufficient flow and constant feed pressure of 21 bar to keep the piston bearings in continuous contact with the cam ring (Hägglunds, 2015). Load pins are integrated into the suspension of the pump in the nacelle to measure the pump torque.

Oil motor. The highpressure hydraulic discharge line of the oil pump drives a BoschRexroth A6VLM variabledisplacement axial piston oil motor (BoschRexroth, 2012), configured here with a (maximum) fixed displacement.

Water pump. The oil motor is mechanically coupled to a Kamat K801205G fixeddisplacement water plunger pump (KAMAT, 2017). An external centrifugal pump supplies the water pump with sufficient flow and feed pressure of around 2.6 bar.

Spear valve. The pressure in the water pump discharge line is controlled by variablearea orifices in the form of two spear valves. The valves are used to control the system torque and rotor speed in below and nearrated operating conditions and form highspeed water jets towards the Pelton turbine. A schematic visualization of this system is shown in Fig. 4. The spear positions are measured and individually actuated by DC motors. The spear valves are actuated in such a way that only full on–off spear position actuation is possible. This means that either the spear moves forwards–backwards at full speed or the spear valve position remains at its current position. A deadband logic controller is implemented to enable position control of the valve within a predefined band around the set point.

Pelton turbine. When the water flow exits the spear valves, the hydrostatic energy in the highpressure line is converted to a hydrodynamic highvelocity water jet. The momentum of the jet exerts a force on the Pelton turbine buckets (Zhang, 2007). Pelton turbines are highly specialized pieces of equipment and are designed to meet specific conditional requirements (Brekke, 2001). The Sy Sima 315 MW turbine in Norway, for 88.5 bar of head pressure, is to date the largest known (Cabrera et al., 2015). The custommanufactured Pelton turbine employed in the DOT500 setup is designed for the nominal pressure and the optimal speed conditions of the connected electrical generator.
The pitch circle diameter (PCD) of the custommanufactured wheel is 0.85 m and the nominal speed of the turbine is in the range of 1230–1420 RPM. According to the manufacturer, the optimal operational pressure is in the range of 60–80 bar for a nominal flow of approximately 58 L min^{−1} using two spear valves. The optimal ratio between the tangential Pelton and water jet speed is approximately 1∕2 (Thake, 2000) and is maintained by speed control of the asynchronous generator using a filtered pressure measurement. Further considerations on optimal Pelton operation under varying conditions are given below.

Generator. The mechanical rotational energy of the Pelton turbine is harvested by a mechanically coupled generator. As a grid connection was unavailable at the test location, the electrical energy excess is dissipated by a brake resistor.
In addition to the components described in the abovegiven enumeration, auxiliary hardware is present and required for operation of the turbine. A more detailed description is presented in Diepeveen et al. (2018); however, a short summary of the relevant components and remarks is given.

Cooling system. Because of the significant mechanical losses, the working medium needs to be cooled to enable longterm operation of the turbine. The pump and motor in the oil circuit are equipped with flushing lines, which lead to the oil reservoir. A parallel forcedconvection cooling circuit regulates the oil temperature.

Water circuit. After the highvelocity water jet has hit the Pelton buckets, it falls back into the first water reservoir. A second water reservoir is connected to the first by a highvolume line to prevent a disturbed flow to the water pump. From the second reservoir the water is fed to the water pump by a centrifugal pump. The highpressure circuit includes a pressure relief valve.
Furthermore, a remark has to be made on the operational strategy of the Pelton turbine. It is concluded from earlier tests and manufacturer data sheets that Pelton efficiency characteristics are mainly a function of flow and to a lesser extend of line pressure. The flow to the Pelton turbine is proportional to the wind turbine rotor speed, which results in suboptimal operation of the Pelton. The Pelton speed control strategy described above aims to extract the maximum amount of energy given the conditions it is subjected to. For now, this is a design choice, and further research needs to be conducted to elaborate on Pelton design and efficiency maximization given the varying operational conditions.
The theory for model derivation of the DOT500 hydraulic drivetrain is presented in this section. As control actions influence the turbine operating behavior to the point at which the hydrodynamic water jet exits the spear valve, modeling of the turbine drivetrain is performed up to that point. After the water flow exits the spear valve, the aim of operating the Pelton turbine–generator combination at maximum efficiency is a decoupled control objective from the rest of the drivetrain. Considerations on this aspect are given in Sect. 2.2.
A simplified hydraulic diagram of the setup is given in Fig. 5. The components considered in the model derivation are shown in black, whereas auxiliary systems are presented in gray. The symbols in this figure are specified throughout the different parts of this section. In this section, pressures are generally given as a pressure difference Δp over a component, but when the pressure with respect to the atmospheric pressure p_{0} is intended, the Δ indication is omitted.
The organization is divided into two parts. First, a steadystate drivetrain model is derived in Sect. 3.1. This model is used later for the derivation of a passive torque control strategy. Secondly, in Sect. 3.2, a drivetrain model including oil line dynamics is derived and is used for the design of an active spear valve control strategy.
3.1 Steadystate drivetrain modeling
A steadystate model of the drivetrain is derived for hydraulic torque control design in belowrated operating conditions. Mathematical models of hydraulic wind turbines have been established, but mostly incorporate a drivetrain with variabledisplacement components (Buhagiar et al., 2016; Jarquin Laguna, 2017). In Skaare et al. (2011, 2013) a more simple, robust, and efficient drivetrain with a discrete hydraulic gear ratio is proposed by enabling and disabling hydraulic motors. Recently, a feedback control strategy for wind turbines with digital fluid power transmission was described in Pedersen et al. (2017). However, only fixeddisplacement components are considered and modeling such a DOT drivetrain is described in Diepeveen (2013). The model derivation in this section incorporates the components employed in the actual DOT500 setup. A simplified wind turbine model is introduced in Sect. 3.1.1. This model is complemented in Sect. 3.1.2 by analytic models of drivetrain components.
3.1.1 Simplified wind turbine model
The Newton law for rotational motion is employed as a basis for modeling the wind turbine rotor speed dynamics:
where J_{r} is the rotor inertia, ω_{r} the rotor rotational speed, τ_{r} the mechanical torque supplied by the rotor to the lowspeed shaft, and τ_{sys} the system torque supplied by the hydraulic oil pump to the shaft. The rotor inertia J_{r} of the rotor is not publicly available. However, an estimation of the rotor inertia is obtained using an empiric relation on blade length given in Rodriguez et al. (2007), resulting in a value of 6.6×10^{5} kg m^{2}. Moreover, experiments were performed on the actual turbine and confirm this theoretical result (Jager, 2017). The torque supplied by the rotor (Bianchi et al., 2006) is given by
where the density of air ρ_{air} is taken as a constant value of 1.225 kg m^{−3}, U is the velocity of the upstream wind, and R is the blade length of 22 m. The power coefficient C_{p} represents the fraction between the captured rotor power P_{r} and the available wind power P_{wind} and is a function of the blade pitch angle β and the dimensionless tipspeed ratio λ given by
The power coefficient C_{p} is related to the torque coefficient by ${C}_{\mathit{\tau}}(\mathit{\lambda},\mathit{\beta})={C}_{\mathrm{p}}(\mathit{\lambda},\mathit{\beta})/\mathit{\lambda}$ such that Eq. (2) can be rewritten as
The rotor power and torque extraction capabilities from the wind are characterized in respective power and torque coefficient curves. These curves of the actual DOT500 rotor are generated by mapping the actual blade airfoils and applying blade element momentum (BEM) theory (Burton et al., 2011); they are given in Fig. 6 at the blade finepitch angle. The finepitch angle β_{0} indicates the blade angle resulting in maximum rotor power extraction in the belowrated operating region (Bossanyi, 2000). The theoretical maximum rotor power and torque coefficients equal ${C}_{\mathrm{p},\mathrm{max}}=\mathrm{0.48}$ and ${C}_{\mathit{\tau},max}=\mathrm{7.2}\times {\mathrm{10}}^{\mathrm{2}}$ at tipspeed ratios of 7.8 and 5.9, respectively. The complete power, torque, and thrust coefficient data set is available as an external supplement under Mulders et al. (2018b).
The system torque τ_{sys} is supplied by the hydraulic drivetrain to the rotor lowspeed shaft. This torque is influenced by the components in the drivetrain, which all have their own energy conversion characteristics expressed in efficiency curves. All components are off the shelf and their combined efficiency characteristics influence the operating behavior of the turbine.
Hydraulic components are known for their high torquetoinertia ratio and have high acceleration capabilities as a result (Merritt, 1967). In typical applications of a hydraulic transmission, the fairly low rotational inertia of pumps and motors is still relevant. However, the considered wind turbine drivetrain is driven by a rotor with a large inertia J_{r} compared to the drivetrain components. Referring to the specification sheet of the oil motor (BoschRexroth, 2012), it is stated that the unit has a moment of inertia of J_{b}=0.55 kg m^{2}. The resulting reflected inertia to the rotor of ${J}_{b\to \mathrm{r}}=\mathrm{0.55}/{G}^{\mathrm{2}}=\mathrm{1533.312}$ kg m^{2} is still negligible, where G^{−1} represents the hydraulic gear ratio of 52.8. Furthermore, a particular study on this aspect has been carried out in Kempenaar (2012), in which it is concluded that the inclusion of component dynamics does not result in significantly improved model accuracy. For the reasons mentioned, the pumps and motor included in the drivetrain are assumed to have negligible dynamics, and the power conversion (flow speed, torque pressure) is given by static relations.
3.1.2 Analytic drivetrain components description
A flow diagram of the modeling strategy is presented in Fig. 7. To obtain an expression for the system torque τ_{sys}, the complete hydraulic flow path with its volumetric losses is modeled first. When the flow path reaches the spear valve at the water discharge to the Pelton turbine, the simulation path is reversed to calculate the effect of all component characteristics to the line pressures. The spear valve allows for control of the water discharge pressure, the effect of which propagates back to the system torque τ_{sys}. The highpressure oil flow by the oil pump is proportional to the rotor speed,
where V_{p,h} is the pump volumetric displacement, and η_{v,h} the volumetric efficiency. The volumetric efficiency indicates the volume loss as a fraction of the total displaced flow as a function of the component operating conditions. However, as volumetric efficiency data are unavailable for most of the components and the aim is to provide a simplified model of the hydraulic drivetrain, volumetric losses are considered as a constant factor of the displaced flow. The resulting oil flow Q_{o} drives the oil motor, which results in a rotational speed of the motor shaft subject to volumetric losses:
where V_{p,b} is the oil motor volumetric displacement, and η_{v,b} the volumetric efficiency of the oil motor. As the water pump is mechanically coupled to the oil motor axis, its rotational speed is equal to ω_{b}. The expression relating the rotational speed to the water pump discharge water flow Q_{w} is given by
where V_{p,k} is the volumetric displacement of the water pump, and η_{v,k} the volumetric efficiency of the water pump. The pressure in the water discharge line is controlled by a spear valve of which a visualization is given in Fig. 8 with the spear position coordinate system. The effective nozzle area is variable according to the relation
where $\left\{s\subset \mathbb{R}\phantom{\rule{0.25em}{0ex}}\phantom{\rule{0.25em}{0ex}}\mathrm{0}\le s\le {s}_{\mathrm{max}}\right\}$ represents the position of the spear in the circular nozzle cross section, D_{nz} is the nominal nozzle diameter, α the spear coning angle, and N_{s} indicates the amount of spear valves on the same line. Modeling multiple spear valves by N_{s} assumes equal effective nozzle areas for all valves. The maximum spear position (fully open) is given by
and a mapping for spear position to effective nozzle area for different nozzle diameters D_{nz} is given in Fig. 9. The spear valve is closed for all cases at position s=0 mm.
The spear valve converts the hydrostatic water flow into a highspeed hydrodynamic water jet that exerts a thrust force on the buckets of the Pelton turbine (Zhang, 2007). Using the Bernoulli equation for incompressible flows (White, 2011), an expression for the discharge water pressure is obtained:
with the flow and effective nozzle area $\left\{{Q}_{\mathrm{w}},{A}_{\mathrm{nz}}\right\}\subset {\mathbb{R}}^{+}$. As observed in the above given relation, the pressure can be controlled by varying the feed flow and spear position, as the latter influences the effective nozzle area A_{nz}. The discharge coefficient C_{d} is introduced to account for pressure losses due to the geometry and flow regime at the nozzle exit (Al'tshul' and Margolin, 1968). The discharge coefficient of an orifice is defined as the ratio between the vena contracta area and the orifice area (Bragg, 1960). The vena contracta is the point at which the streamlines become parallel, which usually occurs downstream of the orifice at which the streamlines are still converging.
The pressure in the water discharge line propagates back into the system and is used as a substitute for conventional wind turbine torque control. A relation for the mechanical torque at the axis between the water pump and oil motor is given by
where η_{m,k} is the mechanical efficiency of the water pump as a function of the rotational speed and pressure difference over the pump:
where p_{k,f} is a known and constant feed pressure. The torque τ_{b} is used to calculate the pressure difference over the oil motor and pump by
where η_{m,b} is the mechanical efficiency of the oil motor. It is assumed that the pressure at the discharge outlet of the oil motor Δp_{b} is constant as the feed pressure to the oil pump is regulated. Finally, the system torque supplied to the rotor lowspeed shaft is given by
Using the relations derived, a passive strategy for belowrated turbine control is presented in Sect. 4.1.
3.2 Dynamic drivetrain modeling
In contrast to the steadystate model presented previously, this section elaborates on the derivation of a drivetrain model including fluid dynamics for validation of the control design in the nearrated region. First, preliminary knowledge on fluid dynamics is given in Sect. 3.2.1, whereafter a dynamic DOT500 drivetrain model is presented in Sect. 3.2.2.
3.2.1 Analysis of fluid dynamics in a hydraulic line
The dynamics of a volume in a hydraulic line are modeled in this section. For this, analogies between mechanical and hydraulic systems are employed for modeling convenience. A full derivation is given in Appendix B, and only the results are given in this section. The system considered, representing the highpressure discharge oil line, is a cylindrical control volume ${V}_{\mathrm{H}}=A{L}_{\mathrm{l}}=\mathit{\pi}{r}_{\mathrm{l}}^{\mathrm{2}}{L}_{\mathrm{l}}$ with radius r_{l} and length L_{l} excited by a pressure $\mathrm{\Delta}p={p}_{\mathrm{in}}{p}_{\mathrm{out}}$. The net flow into the control volume is defined as $Q={Q}_{\mathrm{in}}{Q}_{\mathrm{out}}$. For a hydraulic expression with pressure Δp as the external excitation input and the flow Q as output, one obtains
where L_{H}, R_{H}, and C_{H} are the hydraulic induction, resistance, and capacitance (Esposito, 1969), respectively, and are defined in Appendix A. The inverse result of Eq. (15) is obtained (Murrenhoff, 2012) with flow Q as the external excitation and Δp as output:
Finally, the differential equation defined by Eq. (15) is expressed as a transfer function,
and the same is done for Eq. (16):
The transfer functions defined in Eqs. (17) and (18) show the characteristics of an inverted notch with +1 and −1 slopes on the left and right side of the natural frequency, respectively. This physically means that exciting the system pressure results in a volume velocity change predominantly at the system natural frequency for the former case. An illustrative Bode plot is given in Sect. 4.2.1. Exciting the flow results in amplification and transmission to pressure in a wider frequency region. This effect is a result of the inverse proportionality between the damping coefficients ζ_{Q} and ζ_{p} (see Appendix B).
3.2.2 Drivetrain model derivation
A dynamic model of the DOT500 drivetrain is derived by application of the theory presented in the previous section. The drivetrain is defined from the rotor up to the spear valve, and the following assumptions are made.

Because of the high torquetoinertia ratio of hydraulic components (Merritt, 1967), the dynamics of oil pumps and motors are disregarded and taken as analytic relations.

Because of the longer line length and higher compressibility of oil compared to the shorter water column, the highpressure oil line is more critical for control design, and a dynamic model is implemented for this column only.

The fluids have a constant temperature.
The dynamic system is governed by the following differential equations:
where K_{e} is the equivalent bulk modulus, including the fluid and line compressibility defined in Eq. (A5), and 𝒱 is the net volume inflow to the considered oil line between the oil pump discharge and oil motor feed port. For convenience, mechanical model quantities are expressed hydraulically in terms of fluid flows and pressure differences over the components. Therefore, the rotor inertia J_{r} is expressed in terms of fluid induction and is combined with the hydraulic induction term into ${L}_{\mathrm{R}}^{*}$.
Both the spear position and pitch angle are modeled by a firstorder actuator model:
where t_{s} and t_{β} are the time constant for the spear valve and pitch actuators, respectively, and the phase loss at the actuator bandwidth is assumed to account for actuation delay effects.
The abovegiven dynamic equations are written in a statespace representation as
It is seen that the pressure difference over the oil pump and motor appear as inputs, but these quantities cannot be controlled directly. For this reason, linear expressions of the rotor torque and spear valves are defined next. The rotor torque is linearized with respect to the tipspeed ratio, pitch angle, and wind speed:
where $\widehat{(\cdot )}$ indicates a value deviation from the operating point, and $\stackrel{\mathrm{\u203e}}{(\cdot )}$ is the value at the operating point (Bianchi et al., 2006). Furthermore,
where the quantities ${k}_{{\mathit{\omega}}_{\mathrm{r}}}$, k_{β}, and k_{U} represent the intrinsic speed feedback gain, the linear pitch gain, and the linear wind speed gain, respectively. The intrinsic speed feedback gain can also be expressed as a function of the tipspeed ratio by
For aerodynamic rotor stability, the value of k_{λ} needs to be negative. In Fig. 10 the intrinsic speed feedback gain ${k}_{\mathit{\lambda}}(\mathit{\lambda},\stackrel{\mathrm{\u203e}}{\mathit{\beta}},\stackrel{\mathrm{\u203e}}{U})$ is evaluated as a function of the tipspeed ratio at the finepitch angle β_{0}. For incorporation of the linearized rotor torque in the drivetrain model, Eq. (25) is expressed in the pressure difference over the oil pump:
where the conversions of the required quantities are given by
Similarly, the water line pressure as defined in Eq. (10) is linearized with respect to the spear position and flow through the valve:
where
The pressure difference over the oil motor is defined in terms of the water line pressure, which is a function of the spear position:
where the mechanical and volumetric conversion factors from oil to water pressure and flow are defined as
The system defined in Eq. (24) is now presented as a linear statespace system of the form
With substitution of the rotor torque and water pressure approximations defined by Eqs. (32) and (37) in Eq. (24), the state A, input B, wind disturbance B_{U}, and output C matrices are given by
with the state, input, and output matrices
and
The dynamic model derived in this section is used in Sect. 4.2 to create an active spear valve torque control strategy in the nearrated region.
In this section designs are presented for control in the below and nearrated operating region. A schematic diagram of the overall control system is given in Fig. 11. It is seen that the turbine is controlled by two distinct proportional–integral (PI) controllers, a spear valve torque, and blade pitch controller acting on individual rotor speed set point errors e_{s} and e_{β}, respectively. As both controllers have a common control objective of regulating the rotor speed and are implemented in a decentralized way, it is ensured that they are not active simultaneously. The gainscheduled pitch controller is designed and implemented in a similar way as in conventional wind turbines (Jonkman et al., 2009) and is therefore not further elaborated in this paper. The spear valve torque controller, however, is nonconventional and its control design is outlined in this section.
For the belowrated control design a passive torque control strategy is employed, a description of which is given in Sect. 4.1. Subsequently, in Sect. 4.2, the infield active spear valve control implementation is evaluated using the dynamic drivetrain model.
4.1 Passive belowrated torque control
The passive control strategy for belowrated operation is described in this section. Conventionally, in belowrated operating conditions, the power coefficient is maximized by regulating the tipspeed ratio at λ_{0} using generator torque control. Generally, the maximum power coefficient tracking objective is attained by implementing the feedforward torque control law
where K_{r} is the optimal mode gain in Nm (rad s^{−1})^{−2}.
As the DOT500 drivetrain lacks the option to directly influence the system torque, hydraulic torque control is employed using spear valves. An expression for the system torque for the hydraulic drivetrain is derived by substitution of Eqs. (11) and (13) in Eq. (14),
and by substituting Eqs. (10) and (12) an expression as a function of the spear position is obtained:
Now substituting Eqs. (5) and (6) in Eq. (7) results in an expression relating the water flow to the rotor speed:
Combining Eqs. (47) and (48) and disregarding the water pump feed pressure p_{k,f} gives
Rotor speed variations cause a varying flow through the spear valve, which results in a varying system pressure and thus system torque, regulating the tipspeed ratio of the rotor. The aboveobtained result shows that when K_{s} is constant, the tipspeed ratio can be regulated in the belowrated region by a fixed nozzle area A_{nz}. Under ideal circumstances, it is shown in Diepeveen and JarquinLaguna (2014) that a constant nozzle area lets the rotor follow the optimal power coefficient trajectory and is called passive torque control. This means that no active control is needed up to the nearrated operating region. For this purpose, the optimal mode gain K_{s} of the system side needs to equal that of the rotor K_{r} in the belowrated region.
However, for the passive strategy to work, the combined drivetrain efficiency needs to be consistent in the belowrated operating region. As seen in Eq. (49), the combined efficiency term is a product of the consequent volumetric divided by the mechanical efficiencies of all components as a function of their current operating point. To assess the consistency of the overall drivetrain efficiency, different operating strategies are examined in Sect. 4.1.1. Subsequently, a component efficiency analysis is given in Sect. 4.1.2.
4.1.1 Operational strategies
Because hydraulic components are in general more efficient in highload operating conditions (Trostmann, 1995), it might be advantageous for a hydraulic wind turbine drivetrain to operate the rotor at a lower tipspeed ratio. Operating at a lower tipspeed ratio results in a lower rotational rotor speed and a higher torque for equal wind speeds, but at the same time decreases the rotor power coefficient C_{p}. By sacrificing rotor efficiency, the resulting higher operational pressures might lead to maximization of the total drivetrain efficiency. For these reasons, an analysis of this tradeoff is divided into two cases.

Case 1: operating the rotor at its maximum power coefficient C_{p,max}.

Case 2: operation at the maximum torque coefficient C_{τ,max}.
Referring back to the rotor power–torque curve in Fig. 6 and substituting the values for operation at C_{p,max} and C_{τ,max} in Eq. (45), optimal mode gain values of ${K}_{\mathrm{r},p}=\mathrm{1.00}\times {\mathrm{10}}^{\mathrm{4}}$ Nm (rad s^{−1})^{−2} and ${K}_{\mathrm{r},\mathit{\tau}}=\mathrm{2.05}\times {\mathrm{10}}^{\mathrm{4}}$ Nm (rad s^{−1})^{−2} are found for cases 1 and 2, respectively. The result of evaluating the rotor torque path in the belowrated region for the two cases is presented in Fig. 12. Due to the lower tipspeed ratio in case 2 , the rotor speed is lower for equal wind speeds or a higher wind speed is required for operation at the same rotor speed, resulting in a higher torque. An efficiency evaluation of the proposed operational cases using actual component efficiency data is given in the next section.
4.1.2 Drivetrain efficiency and stability analysis
This section presents the available component efficiency data and evaluates steadystate drivetrain operation characteristics for the two previously introduced operating cases. The component efficiency characteristics primarily influence the steadystate response of the wind turbine, as shown in Eq. (49). To perform a fair comparison between the two operating cases, the rotor efficiency is normalized with respect to case 1, resulting in a constant efficiency factor of 0.85 for operating case 2. Detailed efficiency data are available for the oil pump and motor; however, as no data for the efficiency characteristics of the water pump are available, a constant mechanical and volumetric efficiency of ${\mathit{\eta}}_{\mathrm{m},k}=\mathrm{0.83}$ and ${\mathit{\eta}}_{\mathrm{m},k}=\mathrm{0.93}$ are assumed, respectively.
The oil pump is supplied with total efficiency data η_{t,h} as a function of the (rotor) speed ω_{r} and the supplied torque τ_{sys} (Hägglunds, 2015). An expression relating the mechanical, volumetric, and total efficiency is given by
where η_{v,h} is taken as 0.98, and the result is presented in Fig. 13 (left). The plotted data points (dots) are interpolated on a mesh grid using a regular grid linear interpolation method from the Python SciPy interpolation toolbox (SciPy.org, 2017). Operating cases 1 and 2 are indicated by the solid lines. The same procedure is performed for the data supplied with the oil motor, the result of which is presented in Fig. 13 (right), in which η_{v,b} is taken as 0.98. As concluded from the efficiency curves, hydraulic components are generally more efficient in the lowspeed, hightorque, and/or highpressure region. It is immediately clear that for both the oil pump and the motor, operating the rotor at a lower tipspeed ratio (case 2) is beneficial from a component efficiency perspective.
The drivetrain efficiency analysis for both operating cases is given in Fig. 14. The lack of efficiency data at lower rotor speeds in the left plot of Fig. 14 (case 1) is due to the unavailability of data at lower pressures. From the resulting plot it is concluded that the overall drivetrain efficiency for case 2 is higher and more consistent compared to case 1. The consistency of the total drivetrain efficiency is advantageous for control, as this enables passive torque control to maintain a constant tipspeed ratio. As a result of this observation, the focus is henceforth shifted to the implementation of a torque control strategy tracking the maximum torque coefficient (case 2).
It should be stressed that this operational strategy is beneficial for the considered drivetrain, but can by no means be generalized for other wind turbines with hydraulic drivetrains. As the overall efficiency of hydraulic components is a product of mechanical and volumetric efficiency, a more rigorous approach would be to optimize the ideal belowrated operational trajectory subject to all component characteristics. However, to perform a more concise analysis, only the two given trajectories are evaluated.
A stability concern for operation at the maximum torque coefficient needs to be highlighted. For stable operation, the value of k_{λ} needs to be negative. As shown in Fig. 10, the stability boundary is located at a tipspeed ratio of 5.9. Operation at a lower tipspeed ratio results in unstable turbine operation and deceleration of the rotor speed to standstill. However, as concluded in Schmitz et al. (2013), hydraulic drivetrains can compensate for this theoretical instability, allowing for operation at lower tipspeed ratios. Therefore, the case 2 torque control strategy is designed for the theoretical calculated minimum tipspeed ratio of 5.9, and infield test results need to confirm the practical feasibility of the implementation.
4.2 Active nearrated torque control
A feedback hydraulic torque control is derived for nearrated operation in this section. To this end, active spear position control is employed to regulate the rotor speed. The effect on fluid resonances is analyzed, as these are possibly excited by an increased rotor speed control bandwidth. The infield tests with corresponding control implementations are performed prior to the theoretical dynamic analysis of the drivetrain. For this reason, the control design and tunings used infield are evaluated and possible improvements are highlighted. Section 4.2.1 and 4.2.2 define the modeling parameters of the oil column and spear valve actuator, which is used in Sect. 4.2.3 for spear valve torque controller design.
4.2.1 Defining the hydraulic model parameters
The highpressure oil line in the DOT500 drivetrain is considered to contain SAE30 oil, with a density of ρ_{o}=900 kg m^{−3} and an effective bulk modulus of ${K}_{\mathrm{f},\mathrm{o}}=\mathrm{1.5}$ GPa. The hydraulic line is cylindrical with a length of L_{l}=50 m, a radius of r_{l}=43.3 mm, and a bulk modulus of K_{l}=0.80 GPa (Hružík et al., 2013). According to Eq. (A5), the equivalent bulk modulus becomes K_{e}=0.52 GPa. The dynamic viscosity of SAE30 oil is taken at a fixed temperature of 20 ^{∘}C, at which it reads a value of μ=240 mPa s. With these data the hydraulic inductance, resistance, and capacitance have calculated values of ${L}_{\mathrm{H}}=\mathrm{7.64}\times {\mathrm{10}}^{\mathrm{6}}$ kg m^{−4}, ${R}_{\mathrm{H}}=\mathrm{8.69}\times {\mathrm{10}}^{\mathrm{6}}$kg m^{−4} s^{−1}, and $\mathrm{1}/{C}_{\mathrm{H}}=\mathrm{1.77}\times {\mathrm{10}}^{\mathrm{9}}$ kg m^{−4} s^{−2}, respectively. Using Eq. (B4), the flow is calculated to be laminar as Re=1244 with oil flow at a nominal rotor speed of 1478 L min^{−1}, and thus a correction factor ${f}_{\mathrm{c}}=\mathrm{4}/\mathrm{3}$ is applied to the hydraulic inductance.
By substitution of the calculated values in Eq. (17), a visualization of the transfer function frequency response is given in Fig. 15 at a range of hydraulic line lengths. In this Bode plot, it is shown that the line length has a great influence on the location of the natural frequency and damping coefficient. A longer line shifts the frequency to lower values and increases the damping coefficient.
4.2.2 Modeling spear valve characteristics
For determining the time constant t_{s} of the spear valve actuator model defined in Eq. (22), a generalized (pseudorandom) binary noise (GBN) identification signal (Godfrey, 1993) is supplied to one of the spear valve actuators. From this test it is seen that the actuator has a fixed and ratelimited positioning speed of ${\dot{s}}_{\mathrm{max}}=\mathrm{0.44}$ mm s^{−1} and shows no observable transient response.
Because of the nonlinear ratelimited response, an actuator model is parameterized for the worstcase scenario. This is done by evaluating the response at maximum actuation amplitude and determining the corresponding bandwidth such that closedloop reference position tracking is ensured. As concluded from infield experiments, the spear position control range in the nearrated operating region is 1.5 mm, and the corresponding time constant for reference tracking is given by t_{s}=1.69 s.
The spear position relates in a nonlinear fashion to the applied system torque as a consequence of the spear valve geometry presented in Figs. 8 and 9. Therefore, in Fig. 16, an evaluation of the spear valve pressure gradient with respect to the spear position k_{s,s} is given. This is done for distinct nozzle head diameters at a range of effective nozzle areas. It is shown that the spear pressure gradient with respect to the position is higher for larger nozzle head diameters at equal effective areas.
4.2.3 Torque control design and evaluation
The active spear valve torque control strategy employed during the infield tests is now evaluated. A fixedgain controller was implemented for rotor speed control in the nearrated region. As the goal is to make a fair comparison and evaluation of the infield control design, the same PI controller is used in this analysis.
The dynamic drivetrain model presented in Sect. 3.2 is further analyzed. The linear statespace system in Eq. (40) is evaluated at different operating points in the nearrated region. The operating point is chosen at a rotor speed of ${\mathit{\omega}}_{\mathrm{r},s}=\mathrm{27}$ RPM, which results in a water pump discharge flow of ${\stackrel{\mathrm{\u203e}}{Q}}_{\mathrm{w}}=\mathrm{2965}$ L min^{−1}, taking into account the volumetric losses. For the entire nearrated region, a range of wind speeds and corresponding model parameters are computed and listed in Table 2. An analysis is performed on the singleinput singleoutput (SISO) openloop transfer system with the speed error e_{s} as input and rotor speed as output, including the spear valve PI torque controller used during field tests. The gains of the PI controller were ${K}_{\mathrm{p}}=\mathrm{3.8}\times {\mathrm{10}}^{\mathrm{3}}$ m (rad s^{−1})^{−1} and ${K}_{\mathrm{i}}=\mathrm{6.6845}\times {\mathrm{10}}^{\mathrm{4}}$ m rad^{−1}.
By inspection of the state A matrix, various preliminary remarks can be made regarding system stability and drivetrain damping. First it is seen that for the (2, 2) element, the hydraulic resistance R_{H} influences the intrinsic speed feedback gain ${k}_{{Q}_{\mathrm{in}}}^{*}$. It was concluded in Sect. 4.1.2 that the rotor operation is stable for negative values of ${k}_{{\mathit{\omega}}_{\mathrm{r}}}$ and thus ${k}_{{Q}_{\mathrm{in}}}^{*}$. Thus, the higher the hydraulic resistance, the longer the (2, 2) element stays negative for decreasing tipspeed ratios, resulting in increased operational stability. This effect has been observed earlier in Schmitz et al. (2012, 2013), in which it was shown that turbines with a hydraulic drivetrain are able to attain lower tipspeed ratios. This is in accordance with Eqs. (B7) and (B8), in which it is shown that the resistance term only influences the damping. Furthermore, the spear valve pressure feedback gain ${k}_{\mathrm{s},{Q}_{\mathrm{w}}}$ provides additional system torque when the rotor has a speed increase or overshoot, resulting in additional damping to the (3, 3) element in the state matrix.
During analysis, an important result is noticed and is shown by discarding and including the spear valve pressure feedback gain ${k}_{\mathrm{s},{Q}_{\mathrm{w}}}$ in the linearized statespace system. The openloop Bode plot of the system excluding the term is presented in Fig. 17 and including the term in Fig. 18. It is noted that the damping term completely damps the hydraulic resonance peak in the oil column as a result of the intrinsic flow pressure feedback effect. At the same time, the attainable bandwidth of the hydraulic torque control implementation is limited. This bandwidthlimiting effect becomes more severe by applying longer line lengths and thus larger volumes in the discharge line, as it increases the hydraulic induction. Fortunately, various solutions are possible to mitigate this effect and are discussed next.
The effect of increasing different model parameters is depicted in Fig. 18. In the DOT500 setup an intermediate oil circuit is used, which is omitted in the ideal DOT concept. Seawater has a higher effective bulk modulus compared to oil and this has a positive effect on the maximum attainable torque control bandwidth for equal line lengths. Moreover, a faster spear valve actuator has a positive influence on the attainable control bandwidth.
From both figures it is also observed that the system dynamics change according to the operating point. At higher wind speeds, the magnitude of response results is higher. The effect was seen earlier in Fig. 16, in which the pressure gradient with respect to the spear valve position is higher for smaller effective nozzle areas. It is concluded that the sizing of the nozzle diameter is a tradeoff between the minimal achievable pressure and its controllability with respect to the resolution of the spear positioning and actuation speed.
The loopshaped frequency responses attain a minimum and maximum bandwidth of 0.35 and 0.43 rad s^{−1} with phase margins (PM) of 64 and 58^{∘}, respectively. For later control designs a more consistent control bandwidth can be attained by gainscheduling the controller gains on a measurement of the water pressure or spear position. Closedloop step responses throughout the nearrated region (for different wind speeds) are shown in Fig. 19. It is seen that the controller stabilizes the system, and the control bandwidth increases for higher wind speeds.
This section covers the implementation and evaluation of the derived control strategies on the realworld infield DOT500 turbine. In accordance with the previous section, a distinction is made between passive and active regulation for below and nearrated wind turbine control, respectively. To illustrate the overall control strategy from infield gathered test data, operational visualizations are given in Sect. 5.1. Evaluation of the effectiveness of the passive belowrated and active nearrated torque control strategy is presented in Sect. 5.2.
5.1 Turbine performance characteristics and control strategy
To illustrate the overall control strategy, threedimensional rotor torque and rotor speed plots are shown as a function of spear position and wind speed in Fig. 20. By fixing the valve position at a range of positions, making sure that sufficient data are collected throughout all wind speeds and binning the data accordingly, a steadystate drivetrain performance mapping is derived. Both figures show the data binned in predefined spear valve positions and wind speeds. This is done on a normalized scale, on which 0 % is the maximum spear position (larger nozzle area) and 100 % the minimum spear position (smaller nozzle area). The absolute difference between the minimum and maximum spear position is 1.5 mm. The spear position is the only control input in the belowrated region and is independent from other system variables. During data collection, the pitch system regulates the rotor speed up to its nominal set point value when an overspeed occurs.
In both figures the control strategy is indicated by colored lines. For belowrated operation (red), the spear valve position is kept constant: flow fluctuations influence the water discharge pressure and thus the system torque. In nearrated conditions (blue), the spear position is actively controlled by a PI controller. Under conditions of a constant regulated water flow corresponding to ${\mathit{\omega}}_{\mathrm{r},s}=\mathrm{27}$ RPM, the controller continuously adjusts the spear position and thus water discharge pressure to regulate the rotor speed. Once the turbine reaches its nominal power output, the rotor limits wind energy power capture using gainscheduled PI pitch control (green), maintaining the rotor speed at ${\mathit{\omega}}_{\mathrm{r},\mathit{\beta}}=\mathrm{28}$ RPM.
5.2 Evaluation of the control strategy
Previously, in Sect. 4.1, drivetrain characteristics are deduced from prior component information throughout the wind turbine operating region. Characteristic data from the rotor, the oil pump, and the oil motor are evaluated to generate a hydraulic torque control strategy. Due to the predicted consistent mechanical efficiency of the hydraulic drivetrain, the nozzle area can be fixed and no active torque control is needed in the belowrated region. This results from the rotor speed being proportional to water flow and relates to system torque according to Eq. (49). From the analysis it is also concluded that operating at a lower tipspeed ratio results in a higher and more consistent overall efficiency for this particular drivetrain.
As concluded in Sect. 4.1.2, stable turbine operation is attained when the rotor operates at a tipspeed ratio such that k_{λ} is negative, and Fig. 10 shows that the predicted stability boundary is located at a tipspeed ratio of λ=5.9. Using the data obtained from infield tests, a mapping of the attained tipspeed ratios as a function of the spear valve position and rotor speed is given in Fig. 21. An anemometer on the nacelle and behind the rotor measures the wind speed. As turbine wind speed measurements are generally considered less reliable (Østergaard et al., 2007) and the effect of induction is not included in this analysis, the obtained results serve as an indication of the turbine behavior. The dashed line indicates the fixed spear position of 70 % and is chosen as the position for passive torque control in the belowrated region. The attained tipspeed ratio averages are presented in Fig. 21a, and Fig. 21b shows a twodimensional visualization of the data indicated by the dashed line, including 1 standard deviation. Closing the spear valve further for operation at an even lower tipspeed ratio and higher water pressures resulted in a slowly decreasing rotor speed and thus unstable operation.
In Fig. 21, it is shown that the calculated tipspeed ratio is regulated around a mean of 5.5 for belowrated conditions. Although the attained value is lower than the theoretical calculated minimum tipspeed ratio of 5.9, stable turbine operation is attained during infield tests. A plausible explanation is that the damping characteristics of hydraulic components compensate for instability as predicted in Sect. 4.2.3.
Figure 22 shows reaction torque measurements by the load pins in the suspension of the oil pump to estimate the attained rotor torque during belowrated operation. From the tipspeed ratio heat map and the rotor torque measurements it is concluded that the case 2 (maximum rotor torque coefficient) strategy works on the actual turbine, and the passive strategy regulates the torque close to the desired predefined path. However, as seen in both figures for lower rotor speeds, the tipspeed ratio attains lower values and the rotor torque increases. An explanation for this effect is the decreased mechanical water pump efficiency, the efficiency characteristics of which, as stated earlier, are unknown. An analysis of the water pump efficiency is performed using measurement data and shows nonconstant mechanical efficiency characteristics: the efficiency drops rapidly when the rotor speeds is below 15 RPM.
Finally, the active spear valve torque control strategy is evaluated. The aim is to regulate the rotor speed to a constant reference speed in the nearrated operating region. Infield test results are given in Fig. 23, in which the environmental conditions were such that the turbine operated around the nearrated region. All values are normalized for convenient presentation. It is shown that active spear valve control combined with pitch control regulates the wind turbine for (near)rated conditions in a decentralized way. Small excursions to the belowrated region are observed when the spear valve position saturates at is minimal normalized value. The spear position tracks the control signal reference and shows that the strategy has sufficient bandwidth to act as a substitute for conventional torque control.
This paper presents the control design for the intermediate DOT500 hydraulic wind turbine. This turbine with a 500 kW hydraulic drivetrain is deployed infield and served as a proof of concept. The drivetrain included a hydraulic transmission in the form of an oil circuit, as at the time of writing a lowspeed hightorque seawater pump was not commercially available and is being developed by DOT.
First it is concluded that for the employed drivetrain, operating at maximum rotor torque, instead of maximum rotor power, is beneficial for drivetrain efficiency maximization. This results not only in an increased overall efficiency, but comes with an additional advantage of the efficiency characteristics being consistent for the considered drivetrain. From a control perspective, a consistent overall drivetrain efficiency is required for successful application of the belowrated passive torque control strategy. Another benefit of the hydraulic drivetrain is the added damping, enabling operation at lower tipspeed ratios. It is shown using infield measurement data that the passive strategy succeeds in tracking the torque path corresponding to maximum rotor torque for a large envelope in the belowrated region. For a smaller portion, the combined drivetrain mechanical efficiency drops, which results in deviation from the desired trajectory.
Secondly, a drivetrain model including the oil dynamics is derived for spear valve control design in the nearrated region. It is shown that by including a spear valve as a control input, the hydraulic resonance is damped by the flowinduced spear valve pressure feedback. This intrinsic pressure feedback effect also limits the attainable torque control bandwidth. However, this limiting effect can be coped with by using a stiffer fluid, a decreased hydraulic line volume, or a faster spear valve actuator. The sizing of the nozzle head diameter influences the pressure sensitivity with respect to the spear position and affects the attainable control bandwidth by spear valve actuation speed constraints and positioning accuracy. Infield test results show the practical feasibility of the strategy including spear valve and pitch control inputs to actively regulate the wind turbine in the (near)rated operating region. Future control designs will be improved by including a control implementation taking into account the varying spear valve pressure gradient. This will result in a higher and more consistent system response.
The ideal DOT concept discards the oil circuit and only uses water hydraulics with an internally developed seawater pump. As a result, the control design process is simplified and the overall drivetrain efficiency should be greatly improved. Future research will focus on the design of a centralized control implementation for DOT wind turbines acting in a hydraulic network.
The data set and code used in this paper are available under https://doi.org/10.5281/zenodo.1405387 (Mulders et al., 2018b).
Hydraulic induction. The hydraulic induction L_{H} resembles the ease of acceleration of a fluid volume and is related to the fluid inertia I_{f} by
with the assumption that the flow speed profile is radially uniform (Akers et al., 2006). For this reason, a distinction should be made between laminar and turbulent flows in circular lines: the induction of a laminar flow is generally corrected by a factor ${f}_{\mathrm{c}}=\mathrm{4}/\mathrm{3}$, whereas a turbulent flow does not need correction with respect to the fluid inertia I_{f} (Bansal, 1989).
Hydraulic resistance. The hydraulic resistance dissipates energy from a flow in the form of a pressure decrease over a hydraulic element. In most cases, hydraulic resistances are taken as an advantage by means of control valves. For example, by adjusting a valve set point, one adjusts the resistance to a desired value. Mathematically, the hydraulic resistance relates the flow rate to the corresponding pressure drop,
analogous to an electrical circuit in which the voltage over a resistive element equals the current times the resistance. The hydraulic resistance for a hydraulic line with a circular cross section and a laminar flow is
which is a constant term independent of the flow rate. For a turbulent fluid flow, the computation of the resistance is more involved and results in a quantity that is dependent on the flow rate and effective pipe roughness. For simulation purposes this would require reevaluation of the resistance in each time step or for each operating point during linear analysis. Such a nonlinear timevariant (NLTV) system is employed in Buhagiar et al. (2016), updating the resistive terms in each iteration for a hydraulic variabledisplacement drivetrain with seawater under turbulent conditions.
Hydraulic capacitance. Due to fluid compressibility and line elasticity, the amount of fluid can change as a result of pressure changes in a control volume. The effective bulk modulus K_{f} of a fluid is defined by the pressure increase to the relative decrease in the volume:
The equivalent bulk modulus (Merritt, 1967) of a compressible fluid without vapor or entrapped air in a flexible line with bulk modulus K_{l} is defined as
Subsequently, the pressure change with respect to time is
and thus the hydraulic capacitance C_{H} is directly proportional to the volume amount and gives the pressure change according to a net flow variation into a control volume.
For modeling the dynamics of a volume in a hydraulic line, analogies between mechanical and hydraulic systems are employed for convenience. First, the differential equation for a standard mass–damper–spring system driven by an external force F is given by
For conversion to a hydraulic equivalent expression, the driving mechanical force is substituted by F=ΔpA, the control volume mass is taken as $m=\mathit{\rho}{V}_{\mathrm{H}}=\mathit{\rho}A{L}_{\mathrm{l}}$, and the fluid inflow velocity defined as $\dot{x}=Q/A$. By rearranging terms, one obtains
which is further simplified into
where L_{H}, R_{H}, and C_{H} are the hydraulic induction, resistance, and capacitance (Esposito, 1969), respectively, and are defined in Appendix A. The former two of these three quantities depends on the flow Reynolds number, which shows whether the inertial or viscosity terms are dominant in the Navier–Stokes equations (Merritt, 1967). The Reynolds number is defined as
where D_{l}=2r_{l} is the line diameter, and μ the fluid dynamic viscosity. For Reynolds numbers larger than 4000 the flow is considered as turbulent and the inertial terms are dominant, whereas for values smaller than 2300 the viscosity terms are deemed dominant.
For evaluation of the natural frequency ω_{n} and damping coefficient ζ for the considered system, the characteristic equation by neglecting the external excitation force (Δp=0) is defined as
Evaluating the quantities ω_{n} and ζ results in
The inverse result of Eq. (B3) is obtained (Murrenhoff, 2012), with flow Q as the external excitation and Δp as output:
Now by evaluating the characteristic equation
and using Eq. (B6), it is seen that the natural frequency remains unchanged with the result obtained in Eq. (B7), but the definition of the damping coefficient changes:
Finally, the differential equation defined by Eq. (B3) is expressed as a transfer function in
and the same is done for Eq. (B9):
The authors declare that they have no conflict of interest.
The research presented in this paper was part of the DOT500 ONT project,
which was conducted by DOT in collaboration with Delft University of
Technology and executed with funding received from the Ministerie van
Economische zaken via TKI Wind op Zee, Topsector Energie.
Edited by: Luciano Castillo
Reviewed by: five
anonymous referees
Akers, A., Gassman, M., and Smith, R.: Hydraulic power system analysis, CRC Press, Boca Raton, Florida, 2006. a
Al'tshul', A. D. and Margolin, M. S.: Effect of vortices on the discharge coefficient for flow of a liquid through an orifice, Hydrotechnical Construction, 2, 507–510, 1968. a
Artemis Intelligent Power: available at: http://www.artemisip.com/sectors/wind/ (last access: 6 September 2018), 2018. a
Bansal, R.: A Text Book of Fluid Mechanics and Hydraulic Machines, in: MKS and SI Units, Laxmi publications, New Delhi, India, 1989. a
Bianchi, F. D., De Battista, H., and Mantz, R. J.: Wind turbine control systems: principles, modelling and gain scheduling design, Springer, Science & Business Media, 2006. a, b
BoschRexroth: Axial Piston Variable Motor A6VM – Sales information/Data sheet, Tech. rep., BoschRexroth, 2012. a, b, c
Bossanyi, E. A.: The design of closed loop controllers for wind turbines, Wind Energy, 3, 149–163, 2000. a
Bragg, S.: Effect of compressibility on the discharge coefficient of orifices and convergent nozzles, J. Mech. Eng. Sci., 2, 35–44, 1960. a
Brekke, H.: Hydraulic turbines: design, erection and operation, Norwegian University of Science and Technology (NTNU), Trondheim, Norway, 2001. a
Buhagiar, D. and Sant, T.: Steadystate analysis of a conceptual offshore wind turbine driven electricity and thermocline energy extraction plant, Renew. Energ., 68, 853–867, 2014. a
Buhagiar, D., Sant, T., and Bugeja, M.: A comparison of two pressure control concepts for hydraulic offshore wind turbines, J. Dyn. Syst.T. ASME, 138, 081007, https://doi.org/10.1115/1.4033104, 2016. a, b
Burton, T., Jenkins, N., Sharpe, D., and Bossanyi, E.: Wind energy handbook, John Wiley & Sons, Chichester, UK, 2011. a, b
Cabrera, E., Espert, V., and Martínez, F.: Hydraulic Machinery and Cavitation: Proceedings of the XVIII IAHR Symposium on Hydraulic Machinery and Cavitation, Springer, Berlin, Germany, 2015. a
Chapple, P., Dahlhaug, O., and Haarberg, P.: Turbine driven electric power production system and a method for control thereof, Patent US007863767B2, 2011. a
Diepeveen, N.: On the application of fluid power transmission in offshore wind turbines, PhD thesis, Delft University of Technology, Delft, the Netherlands, 2013. a, b, c
Diepeveen, N. and JarquinLaguna, A.: Wind tunnel experiments to prove a hydraulic passive torque control concept for variable speed wind turbines, IOP Publishing, J. Phys. Conf. Ser., 555, https://doi.org/10.1088/17426596/555/1/012028, 2014. a
Diepeveen, N., Mulders, S., and van der Tempel, J.: Field tests of the DOT500 prototype hydraulic wind turbine, in: 11th International Fluid Power Conference, 530–537, Aachen, 2018. a, b
Esposito, A.: A simplified method for analyzing hydraulic circuits by analogy, Mach. Des., 41, 173–177, 1969. a, b
EWEA: The economics of wind energy, EWEA, Brussels, Belgium, 2009. a
Fingersh, L., Hand, M., and Laxson, A.: Wind turbine design cost and scaling model, Tech. rep., National Renewable Energy Lab (NREL), Golden, Colorado, 2006. a
Godfrey, K.: Perturbation signals for system identification, Prentice Hall International (UK) Ltd., 1993. a
Hägglunds: Compact CB – Product Manual, Lohr am Main, Germany, 2015. a, b, c
Hružík, L., Vašina, M., and Bureček, A.: Evaluation of bulk modulus of oil system with hydraulic line, in: EPJ Web of Conferences, Vol. 45, 01041, EDP Sciences, Les Ulis, France, 2013. a
Jager, S.: Control Design and DataDriven Parameter Optimization for the DOT500 Hydraulic Wind Turbine, Master thesis, Delft University of Technology, 2017. a
Jarquin Laguna, A.: Centralized electricity generation in offshore wind farms using hydraulic networks, PhD thesis, Delft University of Technology, 2017. a, b
Jonkman, J., Butterfield, S., Musial, W., and Scott, G.: Definition of a 5MW reference wind turbine for offshore system development, Tech. rep., National Renewable Energy Lab (NREL), Golden, Colorado, 2009. a
KAMAT: Quintuplex Plunger Pump K800005G, Website, available at: https://www.kamat.de/en/plungerpumps/k800005g.html (last access: 6 September 2018), 2017. a, b
Kempenaar, A.: Small Scale Fluid Power Transmission for the Delft Offshore Turbine, Master's thesis, TU Delft, 2012. a
Kotzalas, M. N. and Doll, G. L.: Tribological advancements for reliable wind turbine performance, Philos. T. Roy. Soc. Lond. A, 368, 4829–4850, 2010. a
Merritt, H. E.: Hydraulic control systems, John Wiley & Sons, Hoboken, New Jersey, 1967. a, b, c, d, e
Mulders, S., Diepeveen, N., and van Wingerden, J.W.: Control design and validation for the hydraulic DOT500 wind turbine, in: 11th International Fluid Power Conference, 200–209, RWTH Aachen University, Aachen, 2018a. a
Mulders, S. P., Diepeveen, N. F. B., and van Wingerden, J.W.: Data set: Control design, implementation and evaluation for an infield 500 kW wind turbine with a fixeddisplacement hydraulic drivetrain, https://doi.org/10.5281/zenodo.1405387, 2018b. a, b
Muller, H., Poller, M., Basteck, A., Tilscher, M., and Pfister, J.: Grid compatibility of variable speed wind turbines with directly coupled synchronous generator and hydrodynamically controlled gearbox, in: Sixth International Workshop on LargeScale Integration of Wind Power and Transmission Networks for Offshore Wind Farms, Delft University of Technology, Delft, the Netherlands, 307–315, 2006. a
Murrenhoff, H.: Grundlagen der Fluidtechnik: Teil 1: Hydraulik, Shaker, Herzogenrath, Germany, 2012. a, b
Nijssen, J., Diepeveen, N., and Kempenaar, A.: Development of an interface between a plunger and an eccentric running track for a lowspeed seawater pump, International Fluid Conference 2018, RWTH Aachen University, Aachen, Germany, 2018. a
Østergaard, K. Z., Brath, P., and Stoustrup, J.: Estimation of effective wind speed, in: Journal of Physics: Conference Series, Vol. 75, IOP Publishing, Bristol, UK, 2007. a
Pedersen, N. H., Johansen, P., and Andersen, T. O.: Optimal control of a wind turbine with digital fluid power transmission, Nonlinear Dynam., 91, 1–17, 2017. a
Piña Rodriguez, I.: Hydraulic drivetrains for wind turbines: Radial piston digital machines, Master's thesis, Delft University of Technology, 2012. a
Ragheb, A. and Ragheb, M.: Wind turbine gearbox technologies, in: Nuclear & Renewable Energy Conference (INREC), 2010 1st International, 1–8, IEEE, Piscataway, New Jersey, 2010. a
Rampen, W.: Gearless transmissions of large wind turbinesthe history and future of hydraulic drives, Artemis IP Ltd., Midlothian, UK, 2006. a
Rodriguez, A. G., Rodríguez, A. G., and Payán, M. B.: Estimating wind turbines mechanical constants, in: Proc. Int. Conf. Renewable Energies and Power Quality (ICREPQ'07), 27–30, Vigo, Spain, 2007. a
Rybak, S.: Description of the 3 MW SWT3 wind turbine at San Gorgonio Pass California, in: Fifth Bien. Wind Energy Conference and Workshop (WW5), Vol. 1, 193–206, Bendix Corporation Energy, Environment and Technology Office, San Gorgonio Pass, California, 1981. a
Sasaki, M., Yuge, A., Hayashi, T., Nishino, H., Uchida, M., and Noguchi, T.: Large capacity hydrostatic transmission with variable displacement, in: The 9th International Fluid Power Conference, Vol. 9, RWTH Aachen University, Aachen, Germany, 2014. a
Schmitz, J., Diepeveen, N., Vatheuer, N., and Murrenhoff, H.: Dynamic transmission response of a hydrostatic transmission measured on a test bench, in: EWEA 2012, EWEA, Copenhagen, Denmark, 2012. a
Schmitz, J., Vukovic, M., and Murrenhoff, H.: Hydrostatic transmission for wind turbines: An old concept, new dynamics, in: ASME/BATH 2013 Symposium on Fluid Power and Motion Control, American Society of Mechanical Engineers, Sarasota, Florida, 2013. a, b
SciPy.org: Interpolation (scipy.interpolate.RegularGridInterpolator), available at: https://docs.scipy.org/doc/scipy0.19.1/reference/interpolate.html (last access: 1 July 2017), 2017. a
Silva, P., Giuffrida, A., Fergnani, N., Macchi, E., Cantù, M., Suffredini, R., Schiavetti, M., and Gigliucci, G.: Performance prediction of a multiMW wind turbine adopting an advanced hydrostatic transmission, Energy, 64, 450–461, 2014. a
Skaare, B., Hörnsten, B., and Nielsen, F. G.: Energy considerations for wind turbines with hydraulic transmission systems, EWEA OFFSHORE, Amsterdam, the Netherlands, 2011. a, b
Skaare, B., Hörnsten, B., and Nielsen, F. G.: Modeling, simulation and control of a wind turbine with a hydraulic transmission system, Wind Energy, 16, 1259–1276, 2013. a, b
Spinato, F., Tavner, P., Van Bussel, G., and Koutoulakos, E.: Reliability of wind turbine subassemblies, IET Renew. Power Gen., 3, 387–401, 2009. a
Thake, J.: The microhydro Pelton turbine manual, Practical Action Publishing, https://doi.org/10.3362/9781780445519, 2000. a
Thomsen, K. E., Dahlhaug, O., Niss, M., and Haugset, S.: Technological advances in hydraulic drive trains for wind turbines, Enrgy. Proced., 24, 76–82, 2012. a
Trostmann, E.: Water hydraulics control technology, CRC Press, Boca Raton, Florida, 1995. a
Van der Tempel, J.: Energy extraction system, has water pump attached to rotor, windmill for pumping water from sea, water system connected to water pump, for passing water pumped from sea, and generator connected to water system, 2009. a
White, F. M.: Fluid Mechanics, McGrawHill, New York, 7th Edn., 2011. a
Zhang, Z.: Flow interactions in Pelton turbines and the hydraulic efficiency of the turbine system, P. I. Mech. Eng. AJ. Pow., 221, 343–355, 2007. a, b
 Abstract
 Introduction
 The DOT500 – prototype turbine with offtheshelf components
 Theory and model derivation of the hydraulic drivetrain
 Controller design
 Implementation of control strategy and infield results
 Conclusions
 Data availability
 Appendix A: Definition of hydraulic induction, resistance, and capacitance for Sect. 3.2.1
 Appendix B: Model derivation of a hydraulic control volume for Sect. 3.2.1
 Competing interests
 Acknowledgements
 References
 Abstract
 Introduction
 The DOT500 – prototype turbine with offtheshelf components
 Theory and model derivation of the hydraulic drivetrain
 Controller design
 Implementation of control strategy and infield results
 Conclusions
 Data availability
 Appendix A: Definition of hydraulic induction, resistance, and capacitance for Sect. 3.2.1
 Appendix B: Model derivation of a hydraulic control volume for Sect. 3.2.1
 Competing interests
 Acknowledgements
 References