Rajae Gaamouche ·Abdelbari Redouane·Imad El harraki ·Bouchra Belhorma ·Abdennebi El Hasnaoui

Abstract This paper presents a contribution related to the control of nonlinear variable-speed marine current turbine (MCT)without pitch operating below the rated marine current speed. Given that the operation of the MCT can be divided into several operating zones on the basis of the marine current speed, the system control objectives are different for each zone. To deal with this issue, we develop a new control approach based on a linear quadratic regulator with variable generator torque. Our proposed approach enables the optimization of the rotational speed of the turbine, which maximizes the power extracted by the MCT and minimizes the transient loads on the drivetrain. The novelty of our study is the use of a real profile of marine current speed from the northern coasts of Morocco. The simulation results obtained using MATLAB Simulink indicate the effectiveness and robustness of the proposed control approach on the electrical and mechanical parameters with the variations of marine current speed.

Keywords Marinecurrentturbine .Two-massmodel .Tipspeedratio .Linearization .Optimalcontrol .Linearquadraticregulator(LQR)

1 Introduction

The majority of renewable energy sources,such as wind and solar energies, have a fluctuating character, which alters the quality of the power injected into the electrical grid (Kumar et al. 2016; Anvari et al. 2016). In fact, the integration of these renewable sources into the electrical grid can lead to the major need for the adjustment of voltage and frequency to overcome eventual electrical disturbances (e.g., voltage droops and short circuits) (Antonova et al. 2012). By contrast, marine current energy is more stable than wind and solar energies and represents a promising solution to deal with the aforementioned issues (Melikoglu 2018; Chen et al.2018).Indeed,marine current energy mainly originates from the tides, which are a direct result of the combination of the Earth’s rotation and the gravitational forces exerted on a water body by the Moon and Sun (Hodur 1997). Thus,predicting the particular location of the main component of these currents for long time periods is possible(Hodur 1997;Rourke et al. 2010).

In addition, according to Thiébaut and Sentchev(2015), marine current resources are abundant in such a scale that their use for energy production can satisfy the world’s energy demand. Marine current turbines (MCTs)are evolving rapidly because more reliable oceanographic data are available and all of the advanced techniques developed, tested, and applied for wind turbines can be adapted to MCTs even if the dynamics of tides are totally different from those of wind (Fox et al. 2018).

In recent years, numerous studies of the hydrodynamic aspects (Mycek et al. 2014; Frost et al. 2015; Blackmore et al. 2016) and the design of specific generators for the marine environment (Benelghali et al. 2012; Seck et al.2018;Chen et al.2019)have been conducted.Meanwhile,other studies focus on the electronic part of power (Pham et al. 2017; Einrí et al. 2019; Omkar et al. 2019; Qian et al. 2019).

To improve the energy conversion performance of MCTs and ensure that tidal energy resources are costeffective, rotor speed regulation, which can be classified into two types, i.e., fixed and variable speeds, is investigated in this study. Jena and Rajendran (2015)and Kettache (Khettache 2019) conducted a comparative study between fixed-speed and variable-speed wind generators. The variable-speed wind turbine is more advantageous than the fixed-speed wind turbine because it increases the energy efficiency, improves the quality of the kinetic energy produced, and stabilizes the fluctuations of the voltage and power of the electrical grid.

In the context of wind turbine technologies, many command algorithms have been designed and developed in the past decades to optimize the operation of turbines to maximize the energy efficiency. These controllers can be classified into four categories, namely, optimal control, nonlinear control, predictive control, and fuzzy logic.

Liu et al. (2016) proposed a fuzzy logic controller for the wind turbine system to achieve the objective of maximum power extraction based on a two-mass model.Bassi et al. (Bassi and Mobarak 2017) examined the variablespeed wind turbine in full and partial load segments through predictive control. This approach considers the maximum power point, with the goal of maximizing the power generated by the wind turbine.

In the litrature, many controllers based on nonlinear control have been used to model wind turbines.Boukhezzar and Siguerdidjane (2010, 2011) used a nonlinear controller to deal with the wind power capture optimization problem while restricting transient loads on the drivetrain components. Prasad et al. (2019) developed a new control approach for the generations and loads in the integrated wind power system. This control approach is based on nonlinear control.

Finally, optimal control has been widely used to model and design the variable-speed wind turbine. Optimal control can be applied in two quadratic linear forms, namely,linear quadratic regulator (LQR) and linear quadratic Gaussian (LQG).

Barrera-Cardenas and Molinas (2012) applied the LQG control to model wind turbines. The proposed LQG control has shown useful properties, good performance, and robustness in controller design which has been applied to wind energy converter systems.Kumar and Stol (2010) proposed and tested the use of Simulink to simulate the LQR controller for wind turbines to achieve better rotor speed regulation. To maximize the energy generated from wind, Fakharzadeh et al.2013) presented a linear control law using the LQR approach. Bayat and Bahmani (2017) addressed both the problems of power regulation and wind turbine control using the LQR approach feedback. Mahmoud and Oyedeji (2016) proposed a new voltage control scheme based on the LQR control design for a grid-connected wind farm. The proposed solution can be conveniently utilized for multi-input multi-output systems.

Although the control design is widely used in wind turbine technology, to the best of our knowledge, only a few studies considered the control of MCTs. Zhou et al.(2013) proposed two control strategies (i.e., speed control and torque control) for the power limitation of the MCT when the speed of the marine current exceeds the nominal value. The torque control strategy limits the generator power to a certain value during the dynamic process. However, the proposed control strategy limits only the power for the high-speed marine current. It does not deal with the entire range of variations of the marine current speed. Toumi et al. (2017) reported on speed control using only the classical proportional–integral correctors. The design of this type of control is robust in the case where the MCT is coupled with the permanent magnet synchronous generator. Nevertheless,this type of control cannot maximize the energy efficiency of the MCT. To cope with these drawbacks, in this study, we propose an efficient LQR design for maximizing the power efficiency and energy capture of the MCT.Our decision to use this controller is motivated by the results obtained in the context of wind turbine technology.

The main contributions of this study are summarized as follows:

· We have used a profile of the daily variation of marine current speed in the northern coasts of Morocco;

· We have developed a hydrodynamic model based on the experimental results of Mason-Jones et al.;

· We have modeled flexible transmission using the twomass model;

· We have applied a linearization technique to regulate the system around a specific operating point;

· We have developed a new optimal control based on the LQR approach that achieves the optimization of the operation of the turbine;

· We have obtained the simulation results of the assembled system using MATLAB Simulink.

The remainder of this paper is organized as follows: in Section 2, we describe the MCT system and the aerodynamic and two-mass models. In Section 3, we present the linearization model for the determination of a state space.In Section 4, we develop a new optimal control based on the LQR approach to find a compromise between the optimization of the energy generated by the turbine and the reduction of the load at the level of two-mass mechanical shaft. In Section 5, we provide an overview of the system control objective implementation of the MCT system in the low-speed and high-speed areas of the marine current. In Section 6, we validate our optimal control performance using MATLAB Simulink. Finally, in Section 7, we conclude this paper and indicate the direction of future work.

2 Marine Current Turbine System Modeling

The marine current energy conversion chain,which is a combination of three subsystems,is modeled in Figure 1:

1) The aerodynamic subsystem(Figure 1(a))consists of several blades(mostly three blades)and a hub.The blades of the MCT extract the kinetic energy from the marine current and convert it into mechanical energy.

2) The mechanical subsystem(Figure 1(b))is composed of the gearbox (also called the drivetrain). This subsystem transforms the rotation rate of the shaft from low rotational speed at the rotor side into high rotational speed at the generator side.

3) The electromechanical subsystem(Figure 1(c))is composed of the generator and a power electronics module, which converts the mechanical energy at the turbine into electrical energy.This subsystem will not be modeled in this work.In fact, the dynamics of electrical machines and power electronics systems are faster than that of the other parts of the MCT. Consequently, the system will be considered a mechanical structure.Therefore,the dynamics of the generator can be disregarded. The Temvalue of the generator is considered equal to the reference value Tem-ref.Peis equal to the product of Temand Ωg:

Figure1 MCT model. (a) Aerodynamic system. (b) Drive train. (c)Generator

We describe the model of each subsystem in the following subsections:

2.1 Model of the Marine Current and Aerodynamic Subsystem

2.1.1 Model of the Resource

Before installing the MCT, it is necessary to model the resources of the considered site. For this, the speed of the marine current can be predicted using two methods(Gaamouche et al. 2018):

1) Direct measurement: This method can be performed using conventional devices (e.g., helical current meters and acoustic doppler current profilers),a model proposed by the SHOM (French Naval Hydrographic and Oceanographic Service),or marine ships.

2) Modelization:Several techniques,such as harmonic analysis method, double cosine method, TideSim, and Tide 2D,can be used to model marine currents.

The lack of data of marine current speed in Morocco leads to difficulties in utilizing this resource.Figure 2 presents the variation of marine current speed for a certain day in April 2018 in the northern coasts of the Moroccan kingdom near Spanish borders.The data measured by the European Marine Observation and Data Network were utilized in this study.This figure shows that we can launch such project in this region.

2.1.2 Aerodynamic Modeling

The aerodynamic power recovered by the rotor of the MCT can be expressed as follows:

In this study,we considered that the MCT is not pitched.Moreover,Cpcan be approximated by an equation depending only on λ,which is defined as the ratio between Ωtand Vm:

Mason-Jones et al.(2012)provided a relevant set of experimental data to validate the theoretical and numerical methods for MCT.In this study,we focused on the Cpcurve according to λ,which is obtained from these experimental data when the turbine rotates at a uniform speed of 3.08 m/s over a diameter of 20 m. We notice that the data presented in Table 1 are extracted from this curve.On the basis of this table,we proposed the interpolation function of Cpto model the turbine.This function can be expressed as follows:

Cpcan be considered a strong nonlinear function of λ.Afterward,we implemented our interpolation function for modeling the turbine using the MATLAB Simulink environment. To validate our proposed model, in this work, we selected the turbine “Guinard Energies P66 3.5 kW.” The obtained curve of Cpis shown in Figure 3. We notice that the maximum Cpvalue is 0.39,which corresponds to a tip speed ratio of 3.67 to a marine current speed of 3 m/s.This value is considered the optimal tip speed ratio(λopt)to achieve maximum power point tracking under rated marine current speeds.

2.2 Model of the Mechanical Subsystem

Figure2 Marine current speed profile

Table1 Experimental data of Cp and λ

The mechanical subsystem (Figure 1(b)) is composed of the gearbox (drivetrain). The gearbox adjusts the speed of the turbine to that of the generator through two shafts, namely,the slow shaft on the turbine side and the fast shaft on the generator side. In the literature, the two types of mechanical transmission models are rigid transmission and flexible transmission.The flexible transmission model attracted our attention because the mechanical coupling between the turbine and the electric machine is modeled using the two-mass model,as shown in Figure 1(b). The two masses are connected to a flexible shaft characterized by the elasticity coefficient of the drive shaft of the blades k and the coefficient of friction of the shaft relative to the gearbox d. For wind turbines, this twomass model is sufficient to properly represent the dynamic behavior of the turbine (Fakharzadeh et al. 2013). For this purpose, the MCT drive system will be modeled using the two-mass model.Then,we write the following equations for the low-speed drivetrain:

Figure3 Cp curve of the MCT

3 Linearized and State Representation of the MCT System

Given that determining the optimal control for nonlinear systems is difficult, deriving a solution for an equivalent linear control system is necessary. For the MCT system,nonlinearity is detected in the aerodynamic torque Taer.For this, we need to adopt a linearization approach for Taer(Eq. (8)) according to Ωt.

The operating point corresponding to the marine current speed Vmis variable.Consequently,linearizing the MCT system around different operating points corresponding to the marine current speed Vmis possible.Thus,we can obtain the derivative of the system from an operating point given the marine current speed Vm:In Figure 4, we can distinguish four operating areas of a variable-speed MCT,in which the index a1of the matrix Ai,jdepends on the variation of the marine current:

1) Zones 1 and 4: The MCT does not provide any power because the marine current speed is lower than the start speed and higher than the rated speed. Thus, the MCT will be stopped.

2) Zone 2:For the low marine current speed(Vm≤VΩt?r):

4 LQR Optimal Control

The objective of optimal control has two main orientations.The first orientation aims to minimize energy, whereas the second orientation seeks to reduce the convergence time of the system.The general objective is to find the optimal control that minimizes the criterion that varies according to the orientation adopted.A large variety of optimal control techniques have been applied to the MCT in a permanent attempt to improve its function and benefit as much as possible from the energy that it can produce. Optimal control is the state feedback of a nonlinear time-invariant system with LQR that has evolved significantly in recent years.The principle of the LQR command is shown in Figure 5 (Khargonekar et al.1990; Anderson et al. 2007; Athans and Falb 2013; Levine 2018).

Figure4 Control regions for a turbine controller

Figure5 Principle of the LQR command

In general,the system model can be expressed as the following state space equation:

where x(t)€?ndenotes the state vector,u(t)€?ndenotes the control vector, y(t) €?qdenotes the output vector, A is the matrix of evolution or state,B is the control or input matrix,and C is the output matrix or measured value.Moreover,the pair (A, B) is assumed to denote that the system is controllable. The process of the feedback regulator needs to follow the state space equation. The simplified problem of the LQR is to find the matrix of the corrector K,which minimizes the function of the cost (or the criterion of performance),as follows:

5 Control Objectives

In this part,we developed the control laws that can be applied to zones 2 and 3 (Figure 4) on the basis of the marine current speed:

· When Vm≤VΩt?r, the main control objectives are as follows:The MCT starts to generate energy at a certain marine current speed.As result,the operation of the MCT is related to the low marine current speed. In this area, the control objective is to operate the turbine at maximum efficiency. To ensure that the power coefficient is maintained at the optimal value Cpmax=Cp(λ0), λ must reach its optimum value λ0(Figure 3). This means that the control in this area acts on the electromagnetic torque of the generator, which in turn acts on the rotational speed and reduces its variations from the reference value Ωt?ref,as expressed in Eq.(30).

The controller must minimize the fluctuations of the mechanical torque Tmec,which affects the quality of the electrical power generated(although this effect is less important in this area than in other areas). The reference equation for maximizing energy conversion is expressed as follows:

· When VΩt?r≤Vm≤Vm?r,the main control objectives are as follows:When the marine current reaches the value of VΩt?r, the MCT corresponds to the intermediate marine current speeds. In this area, the rotational speed of the turbine reaches its nominal value.The main control objective is to reduce the variations of Ωtfrom the nominal value Ωt?rwhile acting on the electromagnetic torque.The reference equation for this area is expressed as follows:

Figure6 Block diagram of the system controlled during operation

Figure 6 presents the block diagram of the system controlled during operation.

6 Simulation Results

To validate our proposed control strategy based on the LQR approach, we conducted simulations using the MATLAB Simulink software. The simulation results of the MCT are depicted in Figure 7.

The considered system is a variable-speed, fixedpitch MCT with a nominal power rating of 3.5 kW.Table 2 lists the values of the parameters used in the simulation.

Table2 Simulation parameters

The aim of the proposed control law is to minimize the quadratic criterion previously expressed in Eq. (35). The weight coefficients used in the quadratic criterion J are derived as follows:

The resolution of Riccatti’s equation(i.e.,Eq.(34))for the investigated system defines the control gain K,which is previously expressed in Eq.(29).The result is the optimal state feedback,which is expressed in Eq.(37)or(38)in the following form:

Figure7 Simulation of the MCT system with LQR regulator using MATLAB Simulink

Fig.8 MCT speed variation and state and output variables in the partial load regime(zone 2)

Figure9 MCT speed variation and state and output variables in the partial load regime(zone 3)

To illustrate the effectiveness of our proposed control system for a variable-speed MCT, we conducted two simulations of zones 2 and 3. The results are presented in Figure 4.

6.1 Partial Load Operation(Zone 2)

In this simulation, we considered the low-speed MCT that corresponds to Vm≤VΩt?r.Figure 8 shows the different parameters that we have optimized to maximize the energy collected by the turbine.Figure 8(a)presents the variation of the marine current velocity profile used to determine in which region control will be applied. Moreover, Figure 8 (d) shows that the proposed LQR control follows the reference values nearly perfectly,with a static error of 0.Moreover,the gain provided by this control is equal to K=[0.2079 ?02896 ?0.0353]when Vm=1 m/s.

Figure 8 (b), (c), and (e) present the variations of three parameters,namely,rotor speed,generator speed,and mechanical torque, respectively. The LQR controller ensures that the system maintains the power coefficient at the optimal value Cpmax.As a result,the parametres can be optimized in zone 2.Figure 8(f)shows the variation of the electrical power according to generator speed and electromagnetic torque.Figure 8(f)clearly shows that LQR maximizes the electrical power.

6.2 Partial Load Operation(Zone 3)

In this simulation,the rotational speed of the MCT reaches its rated value that corresponds to VΩt?r≤Vm≤Vm?r. Figure 9 shows the different parameters that we have optimized to maintain the turbine at its optimum speed.In addition,these parameters enable the generator to turn at a speed that does not fluctuate rapidly(minimize the transient load).Figure 9(a)presents the variation of the marine current velocity profile used to determine in which region control will be applied.In this case,a marine current profile with the average value of Vm=3 m/s was imposed.

Figure 9(d)shows that our proposed LQR control can be applied to zone 3,with a good performance in terms of electromagnetic torque. Furthermore, from this figure, we notice that the static error is always equal to 0 and the gain provided by this control is equal to K=[0.1851?0.2899?0.0337]at Vm=3 m/s.

Figure 9 (b), (c), and (e) show the variations of three parameters,namely,rotor speed,generator speed,and mechanical torque, respectively. We notice that the LQR controller can ensure that the system maintains the rotational speed at the nominal value. Consequently, the simulation results indicate that these three parameters have been optimized with the use of LQR control.In addition,Figure 9(f)presents the variation of the output power of the MCT.We notice that the value of this parameter decreases from the nominal value to 3.5 kW because we used the LQR control.

On the basis of the simulation results obtained for zones 2 and 3,we confirm that our proposed LQR control can be set and used to act only on the electromagnetic torque to automatically regulate the different parameters of the MCT system,such as rotor speed,generator speed,and mechanical torque.

7 Conclusions

In this study,we proposed the LQR optimal controller approach coupled with the reference model of the outputs for the variable-speed MCT. Our solution maximizes power capture where the tidal speed is greater or less than the nominal speed.

In the first step,we introduced the differential equations for modeling the state space representation, which connects the outputs and controls to ensure better precision,and the reference model of the outputs whose main role is to impose the desired dynamic on the response of the system’s outputs,including the power produced,whose response should be rapid.

In the second step, the LQR command, which involves minimizing the quadratic criterion that takes into consideration all of the variables of the system, was presented.Furthermore,the LQR command can ensure the best compromise between the desired performance and the feedback of the command.This principle of criterion minimization leads us to determine the optimal control structure of the state feedback for the linearized MCT system around different operating points.

The obtained simulation results enabled us to conclude that our linear control can provide a satisfactory performance in terms of following the reference values for any marine current speed nearly perfectly. In future work, we will focus on the application of more sophisticated control strategies to deal with the influence of any marine current speed by considering the entire system.

NomenclatureVm,marine current speed,m/s;ρ,water density,kg/m3;Rt,rotor radius,m;P,aerodynamic power,W;Taer,aerodynamic torque,N·m;Pe,electrical power,W;λ,tip speed ratio;Cp,power coefficient;Ωt,rotor speed,rad/s;Ωg,generator speed,rad/s;Ωtt,low-speed shaft,rad/s;Jt, rotor inertia, kg·m2; Jg, generator inertia, kg·m2; d, damping coefficient,N·m/(rad·s);k,stiffness coefficient,N·m/rad;Ωt?r,rotor rated speed; G, gearbox gain; Tem, generator torque, N·m; Tmec, mechanical torque, N·m; MCT, marine current turbine; LQR, linear quadratic regulator

Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing,adaptation,distribution and reproduction in any medium or format,as long as you give appropriate credit to the original author(s)and the source,provide a link to the Creative Commons licence,and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence,unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use,you will need to obtain permission directly from the copyright holder. To view a copy of this licence,visit http://creativecommons.org/licenses/by/4.0/.