Adaptive Takagi–Sugeno Fuzzy Model Predictive Control for Permanent Magnet Synchronous Generator-Based Hydrokinetic Turbine Systems

: This paper presents a sensorless model predictive torque control strategy based on an adaptive Takagi–Sugeno (T–S) fuzzy model for the design of a six–phase permanent magnet synchronous generator (PMSG)–based hydrokinetic turbine systems (PMSG-HTs), which not only provides clean electric energy and stable energy-conversion e ﬃ ciency, but also improves the reliability and robustness of the electricity supply. An adaptive T–S fuzzy model is ﬁrst formed to characterize the nonlinear system of the PMSG before a model predictive torque controller based on the T–S fuzzy model for the PMSG system is employed to indirectly control the stator current and the stator ﬂux magnitude, which improves the performance in terms of anti–disturbance, and achieves maximum hydropower tracking. Finally, we consider two types of tidal current, namely the mixed semidiurnal tidal current and the northwest European shelf tidal current. The simulation results demonstrate that the proposed control strategy can signiﬁcantly improve the voltage–support capacity, while ensuring the stable operation of the PMSG in hydrokinetic turbine systems, especially under uneven tidal current speed conditions.


Introduction
With the international oil crisis and the continuous deterioration of the global environment, many countries are stepping up their development of renewable energy sources, which include tidal, wind, solar, geothermal, and marine energy, as well as biomass, and biofuels. According to Statista [1], electricity from renewable sources is expected to experience the largest growth, of up to 21.66 trillion kilowatt hours in 2050, from almost 7 trillion kilowatt hours in 2018. Indeed, the exploitation of wind and marine energy has dramatically increased in recent years. Here, the former has the disadvantages of the length of operating time and the dependence on the seasonal wind speed variations, while, given that the oceans contain rich tidal current energy and that the tidal current is highly predictable since the high-and low-tide cycles are well known, the latter can provide a stable output and sufficient power to the grid. In addition, the International Renewable Energy Agency (IRENA) statistics report [2] shows the global renewable generation capacity amounted to 2537 GW at the end of 2019, where the hydropower accounted for the largest share of the global total, with a capacity of 1190 GW. Hydroelectric power is a reliable renewable energy source that contributes to the stability and stable electricity delivery by reason of its flexible generation potential. Recently, new technology has applied ocean power from ocean tides, waves, currents, salinity, and ocean temperature difference to increase the amount of electricity generated by water. Due to ocean energy occupying no landspace, advantages of the circular economy and pollution free, it is a high-potential renewable energy source [3].
After the Fukushima nuclear accidents in 2011, the Taiwan government decided to halt construction of the new nuclear power plant, and phase out nuclear power by 2025. Since the Taiwan government's announcement of the "Renewable Energy Development Act" in 2009 [4], the Taiwan Ministry of Economic Affairs (MOEA) has been promoting the development of various renewable energies, especial on wind and hydro power. Moreover, Taiwan is surrounded by the ocean and boasts a coastline of more than 1500 kilometers. Even the Kuroshio flows more than 3000 kilometers through the Taiwan east coast; hence, the ocean energy has great potential for development in Taiwan [5,6]. Therefore, a tidal energy conversion system (TECS) in an offshore floating platform is applicable for the development of an independent renewable energy system. However, Taiwan is still in the development stage for ocean-based renewable energy because of both lack of experience and great technical difficulty [7]. Additionally, the uncertainty surrounding the tidal flow causes random fluctuations in hydropower and may further endanger the operational reliability and efficiency of the grid. The cost of using these technologies is site specific and very dependent on turbine technology. In view of this, the main purpose of this paper is to develop the control strategy of a sensorless PMSG-based hydrokinetic turbine systems (PMSG-HTs), and then to explore the effect of the turbine on the tidal velocity and preliminary results on the power output. Not only can that establish technology in developing hydrokinetic power generation for Taiwan, but also provide (to primarily evaluate) the cost-effectiveness of ocean power generation for the Taiwan government, as well as potential environmental impacts and future perspectives.
The hydrokinetic turbine systems is an electromechanical device that converts kinetic energy from flowing water into electricity. The purpose of the system's control is to maximize the energy yield in view of the variable water velocity, such that the maximum power point for the hydrokinetic turbine occurs at different turbine rotor speeds. At present, a permanent magnet synchronous generator (PMSG) is widely employed as the hydrokinetic turbine driven device in the hydropower systems, largely due to its capacity for reaching the maximum power tracking operation, the flexible additional control via converters, and the high thrust at low speed, as well as its high reliability and the lower maintenance costs. From the control point of view, various control strategies have been proposed for the PMSG-based renewable energy conversion systems, including fuzzy logic control [8], robust control [9], field-orientated control (FOC) [10,11], discrete torque control (DTC) [12], proportional-integral (PI) control [13][14][15], proportional-integral-derivative (PID) control [16], sliding mode control [17], model reference adaptive control (MRAC) [18,19], self-adaptive global harmony search (SGHS) [20], and instantaneous-maximum efficiency tracking (i-MET) [21]. However, obtaining the maximum hydropower using the control schemes based on traditional feedback controllers often depends on the accuracy of the mathematical model of the hydrokinetic turbine systems, implying that it tends to be difficult to establish an accurate dynamic model of the hydropower generation system that incorporates a PMSG because of the hydropower system, often a complex nonlinear system, and may be prone to instability and oscillatory behavior from uncertainty tidal flow. Therefore, an adaptive Takagi-Sugeno (T-S) fuzzy model is presented to express the larger nonlinear tidal-based hydropower system via linear combination. In addition, because this paper considers the sensorless PMSG-based hydro turbine, i.e., the angular speed in PMSG is unknown, and persistent disturbances from the tidal current, the T-S fuzzy process is suitably employed to approximate the nonlinear behavior using input/output data of sensorless PMSG-based hydro turbine.
The model predictive control (MPC) scheme, which is also known as model predictive direct torque control (MP-DTC) [12,22], has been widely applied in the control of hydropower generation systems in recent years, since, unlike the conventional feedback control scheme, the MPC remains independent of the system's mathematical model, while it also has the benefit of rolling optimization and online feedback. Using the stator flux along with torque estimation and prediction, the hydrokinetic turbine systems can achieve good transit and static performance in a wide range of regions, including Energies 2020, 13, 5296 3 of 18 maximum torque per ampere control (MTPA) in low-speed regions, and field-weakening control in high-speed regions. However, the internal and external disturbances related to hydrokinetic turbine systems, which include the modeling error and the external uncertainties in ocean current measurements, not only result in inaccurate feedback correction in the predictive MPC model, but also degrade the control performance and the stability of the hydrokinetic turbine systems. In addition, the drawbacks of the traditional MP-DTC scheme may lead to the problem of a large amount of required calculations due to the weighing factor tuning work. Furthermore, to track the maximum power from the tidal current speed below the rated speed, so as to maximize the electricity generated from the tidal current, capturing the rotor speed, and/or position of the generator, is vital. While this information can be easily obtained by utilizing the speed/position sensor, according to [8], more than 14% of the system failures are directly related to sensor failure, while more than 40% are indirectly related to sensor failure. Moreover, to achieve the real-time speed control of the PMSG, the speed of the generator must be derived as a feedback signal. However, only a small number of specialized speed sensors are appropriate for an unfavorable hydropower operating environment, while these are generally extremely expensive. In view of this, numerous researchers have shifted their focus to sensorless generator control strategies for the PMSG [23][24][25][26][27][28][29]. Here, Abdelrahem et al. presented a sensorless FOC scheme [27], also known as vector control, to control the generator speed and to derive the maximum electromagnetic torque with the minimum current, specifically in variable-speed wind turbine systems. The experimental results demonstrate that the power generation efficiency is outstanding under varied operating conditions and PMSG parameter variations. However, the FOC method often involves high parameter dependence. In fact, the accurate dynamic model of a highly complex PMSG-based power generator system is generally difficult to establish due to the nonlinearity, the strong coupling, the multi-varieties, and the uncertainties surrounding the tidal flow or wind. To estimate the nonlinear characteristics of the generator system, various state observer-based methods, such as robust extended Kalman filter (REKF) [28] and unscented Kalman filter (UKF) [29], have been proposed for the sensorless control in PMSGs. However, while these methods can help to reduce the effect of overall error on the estimation accuracy, the process of iterative calculation is relatively time consuming and difficult to compute.
In this paper, an adaptive Takagi-Sugeno (T-S) fuzzy model predictive control (ATSFMPC) scheme is proposed for the design of a six-phase PMSG-HTs. The main objective is to maintain the maximum hydrokinetic energy extraction, i.e., tracking a specified rotational speed of the PMSG when random fluctuations due to the uncertainty of tidal flow are incorporated. To derive the reference torque and stator flux of the PMSG at the maximum power point, an ocean current meter is employed to measure the ocean flow speed before the maximum power point tracking (MPPT) scheme [30] is used to estimate the reference angular speed of the generator. Next, using the MTPA approach [31], the reference torque and stator flux of the PMSG could be extracted as the tracking target for the model predictive torque control (MPTC) strategy. In addition, to represent the PMSG-based nonlinear hydrokinetic turbine systems, the adaptive T-S fuzzy model is presented to establish the approximate model. The designed scheme is capable of concurrent estimations of the stator current and voltage along the d-axis and the q-axis. Finally, a model predictive torque controller for the PMSG system is employed to indirectly control the stator current and the stator flux magnitude, which improves the control performance and achieves maximum hydropower tracking. The simulation results demonstrate the effectiveness of the proposed control scheme.
The main contributions of the proposed control scheme are as follows: The adaptive Takagi-Sugeno (T-S) fuzzy model is used to express the larger nonlinear sensorless PMSG-based hydro turbine system via linear combination. Based on the T-S fuzzy model, the proposed adaptive T-S Fuzzy model predictive torque control (ATSFMPC) scheme is then utilized to indirectly control the electromagnetic torque and stator flux by selecting an adequate voltage vector selection strategy, which not only provides fast and dynamic response, but also achieves the maximum hydro power tracking for sensorless PMSG-based hydro turbine system. To extract the optimal reference torque and stator flux of PMSG driver at maximum available power of hydropower system, the maximum power point tracking (MPPT) and maximum torque per ampere (MTPA) approaches are used to derive the optimal reference torque and stator flux as tracking target of the proposed ATSFMPC strategy.
To demonstrate the performance of the proposed ATSFMPC approach for PMSG-HTs, two real-world tidal current profiles cases are considered. The simulation results show that the proposed ATSFMPC scheme can efficiently provide the stable output power of hydrokinetic turbine and ensure maximum power output without the effect of the disturbance from the regular tidal currents.

Problem Formulation
The basic constituents of the tidal energy conversion system (TECS) [8,32] include the blades, a turbine generator, an inverter, and a grid, as shown in Figure 1. The turbine generator largely functions to convert tidal energy into electricity. First, the tidal current generates a fluid dynamic torque on the turbine, which is then passed through the drivetrain to the generator rotor shaft. The resulting differential pressure over the blade area produces a force, meaning a variable moment of inertia around the axis of rotation can be generated. Following this, the PMSG is connected to a six-phase inverter that rectifies the current from the PMSG to charge a DC-link capacitor. A second inverter is then used to appropriately boost the DC voltage level before it is converted to AC and the yielded power is exported to the grid. ampere (MTPA) approaches are used to derive the optimal reference torque and stator flux as tracking target of the proposed ATSFMPC strategy.
To demonstrate the performance of the proposed ATSFMPC approach for PMSG-HTs, two realworld tidal current profiles cases are considered. The simulation results show that the proposed ATSFMPC scheme can efficiently provide the stable output power of hydrokinetic turbine and ensure maximum power output without the effect of the disturbance from the regular tidal currents.

Problem Formulation
The basic constituents of the tidal energy conversion system (TECS) [8,32] include the blades, a turbine generator, an inverter, and a grid, as shown in Figure 1. The turbine generator largely functions to convert tidal energy into electricity. First, the tidal current generates a fluid dynamic torque on the turbine, which is then passed through the drivetrain to the generator rotor shaft. The resulting differential pressure over the blade area produces a force, meaning a variable moment of inertia around the axis of rotation can be generated. Following this, the PMSG is connected to a sixphase inverter that rectifies the current from the PMSG to charge a DC-link capacitor. A second inverter is then used to appropriately boost the DC voltage level before it is converted to AC and the yielded power is exported to the grid. In the following sections, we provide a brief introduction to tidal energy conversion systems, including the tidal turbine model and the PMSG model.

Tidal Turbine Model
A hydrokinetic turbine operates under the same principles as a wind turbine and is based on similar design features. First, kinetic energy from the moving tidal current is directly converted into mechanical energy by the rotor blades, which can be expressed as where w P is the turbine output power in the horizontal axis (W), ρ is the fluid density  hydrodynamic losses depend on the rotor construction, and it can be defined as follows [32][33][34]: where 1 6 C C are the coefficients, 3 (1 ( 0.08 )) (0.035 (1  )) is the tip speed ratio, and m ω is the mechanical rotation speed (rad/s) of the rotor shaft. To obtain the maximum power from a range of tidal currents to achieve MPPT, i.e., to identify the maximum power operating point and to produce maximum power by regulating the hydrokinetic turbine systems, the power coefficient must be maintained at its maximum max ( ) , which is achieved by ensuring the In the following sections, we provide a brief introduction to tidal energy conversion systems, including the tidal turbine model and the PMSG model.

Tidal Turbine Model
A hydrokinetic turbine operates under the same principles as a wind turbine and is based on similar design features. First, kinetic energy from the moving tidal current is directly converted into mechanical energy by the rotor blades, which can be expressed as where P w is the turbine output power in the horizontal axis (W), ρ is the fluid density (kg/m 3 ), A T = πr 2 T is the area swept by the rotor blade (m 2 ), r T is the blade radius or tip radius (m) , and v w is the tidal current speed (m/s). C p (λ, β) is the power coefficient and a nonlinear function of the tip speed ratio (λ) and the pitch angle of the turbine (β), which can reflect to various hydrodynamic losses depend on the rotor construction, and it can be defined as follows [32][33][34]: where C 1 ∼C 6 are the coefficients, α = (1/(λ + 0.08β)) − (0.035/(1 + β 3 )), λ = ω m r T /v w is the tip speed ratio, and ω m is the mechanical rotation speed (rad/s) of the rotor shaft. To obtain the maximum power from a range of tidal currents to achieve MPPT, i.e., to identify the maximum power operating Energies 2020, 13, 5296 5 of 18 point and to produce maximum power by regulating the hydrokinetic turbine systems, the power coefficient must be maintained at its maximum C max p = C p (λ opt ), which is achieved by ensuring the blade tip speed ratio remains equal to the optimal value λ opt at a fixed value. In other words, the maximum output power of the generator can be derived when the hydrokinetic turbine runs at the optimal blade tip speed ratio λ opt . Hence, the mechanical reference rotation speed ω re f m can be represented as ω where v w is the tidal current speed.
In this paper, we consider two real-world tidal current speed profiles as input for the hydrokinetic turbine systems, namely the northwest European shelf [24] and the Pentland Firth [25]. Remark 1. to achieve sensorless MPPT for hydrokinetic power generation systems, the rotational speed of the hydrokinetic turbines should be regulated in real time according to different tidal current speeds. Hence, the maximum power P max w can be obtained from Equation (1) when the power coefficient C p approaches a maximum value at a measured tidal current speed. In this paper, the relationship between λ and C p in the fixed β value from Equation (2) is as shown in Figure 2a, where the coefficients in Equation (2) are defined as C 1 = 0.5167, C 2 = 116, C 3 = 0.4, C 4 = 5, C 5 = 21, and C 6 = 0.0068, respectively. Meanwhile, Figure 2b shows the hydrokinetic turbine power in terms of the different maximum amplitudes of tidal current speed.
where w v is the tidal current speed.
In this paper, we consider two real-world tidal current speed profiles as input for the hydrokinetic turbine systems, namely the northwest European shelf [24] and the Pentland Firth [25].
Remark 1: to achieve sensorless MPPT for hydrokinetic power generation systems, the rotational speed of the hydrokinetic turbines should be regulated in real time according to different tidal current speeds. Hence, the maximum power max w P can be obtained from Equation (1) when the power coefficient p C approaches a maximum value at a measured tidal current speed. In this paper, the relationship between λ and p C in the fixed β value from Equation (2) is as shown in Figure 2a, where the coefficients in Equation (2) are defined as 1 0.5167 Figure 2b shows the hydrokinetic turbine power in terms of the different maximum amplitudes of tidal current speed.  Notation 1: in this paper, we only consider the MPPT performance of the generator side; hence, the uncertainties of the grid side are assumed to be equal to those of the generator side. As Figure 2a shows, the power coefficient of the hydrokinetic turbine is larger when the pitch angle is smaller. For example, when the pitch angle β = 0 • , the maximum power coefficient is C max p = 0.48 and the optimal speed ratio is λ op T = 8.1. Therefore, the optimal blade tip speed ratio can be obtained by selecting the mechanical rotation speed ω m and tidal current speed v w , such that the maximum power coefficient C max p and the maximum hydrokinetic turbine power can also be derived. According to the optimum speed ratio λ opt and the maximum power coefficient C max p , the hydrokinetic turbine output power in Equation (1) in terms of the various maximums of tidal current speed v max w can be illustrated as in Figure 2b.

Permanent Magnet Synchronous Generator Model
Notation 1: in this paper, we only consider the MPPT performance of the generator side; hence, the uncertainties of the grid side are assumed to be equal to those of the generator side.

Permanent Magnet Synchronous Generator Model
A surface-seated PMSG is widely used in direct-drive hydrokinetic turbine systems [8,17,33,34]. The mathematical model of a PMSG in a direct quadrature (d-q) reference frame, without considering the armature reaction effect and the saliency, is expressed by and is the generator side converter output voltage vector, with u d and u q being the d-axis and q-axis voltages, respectively, T is the generator current vector, with i d and i q being the d-axis and q-axis currents, respectively, R s is the stator resistance, ω e is the electrical generator rotation speed, and ψ dq (t) = ψ d (t), ψ q (t) T is the flux linkage in the stator of the generator with a d-axis and a q-axis. According to [26], it is well known that the magnet flux linkage with the d-axis and q-axis can be defined as ψ d (t) = L d i d (t) + ψ r and ψ q (t) = L q i q (t), respectively. Hence, Equation (4) can be rewritten in the synchronous d-q reference frame as where L d and L q are the stator d-axis and q-axis inductances, respectively, and ψ r is the permanent-magnet flux linkage. Meanwhile, the dynamic of the mechanical hydrokinetic turbine systems is given by where J r is the overall rotor inertia, B m is the viscous friction coefficient, T m is the mechanical torque of the hydrokinetic turbine, i.e., T m = P w /ω m ; T e is the electromagnetic torque, which is given by where P is the number of pole pairs on the rotor.
Remark 2. in this paper, to reduce the complicating features of the mechanism of PMSG, the air-gap of the PMSG is assumed to be smooth with no slotting effect, i.e., L d = L q . Hence, the electromagnetic torque from Equation (7) can be rewritten as T e (t) = ξi q (t), where ξ = 1.5Pψ r .
Remark 3. the model reference adaptive system (MRAS) [35] strategy is often used to deal with the sensorless speed control of the generator. By comparing the stator currents between the reference model and the adaptive model, with the error given to the adaptive mechanism based on a basic PI controller until the error between the two models is reduced to zero, the rotation speed of the generator can be estimated. Therefore, it is a fairly straightforward task to estimate the electrical generator rotation speed ω e (t) and to then derive ω e (t) = Pω m (t).

Adaptive Takagi-Sugeno Fuzzy Model Predictive Controller Design
In this paper, a novel control architecture for a six-phase PMSG-based hydrokinetic turbine systems (PMSG-HTs) is proposed with the aim of maintaining the maximum hydrokinetic energy extraction under the uncertainty of tidal flow, further improving the control performance and achieving MPPT. Figure 3 presents the block diagram of the proposed adaptive T-S fuzzy model predictive control In this paper, a novel control architecture for a six-phase PMSG-based hydrokinetic turbine systems (PMSG-HTs) is proposed with the aim of maintaining the maximum hydrokinetic energy extraction under the uncertainty of tidal flow, further improving the control performance and achieving MPPT. Figure 3   Notation 2: to implement the following MPTC strategy, the discrete-time domain is used.

Stator Flux Reference Value Estimation
To achieve the follow-up MPTC under MTPA operation, the appropriate stator flux reference value at the maximum power point of the hydrokinetic power generation system must be computed from the relevant torque reference. In the MTPA strategy, a proportional-integral (PI) controller is first used to calculate the torque reference ref e T [14,35], which is is the rotor speed error. It should be noted that the PI controller parameters were chosen based on a trial and error method to ensure good control performance. Next, the d-q axis stator reference currents were determined using the MTPA operation, which can be defined as follows [31,36]: Based on the two-axis stator reference currents, the stator flux reference can be determined as Notation 2: to implement the following MPTC strategy, the discrete-time domain is used.

Stator Flux Reference Value Estimation
To achieve the follow-up MPTC under MTPA operation, the appropriate stator flux reference value at the maximum power point of the hydrokinetic power generation system must be computed from the relevant torque reference. In the MTPA strategy, a proportional-integral (PI) controller is first used to calculate the torque reference T re f e [14,35], which is T re f e (k) = (K P + K I/s)∆ω m (k) (8) where ∆ω m (k) = ω re f m (k) − ω m (k) is the rotor speed error. It should be noted that the PI controller parameters were chosen based on a trial and error method to ensure good control performance. Next, the d-q axis stator reference currents were determined using the MTPA operation, which can be defined as follows [31,36]: Based on the two-axis stator reference currents, the stator flux reference can be determined as Energies 2020, 13, 5296 8 of 18

Adaptive Takagi-Sugeno Fuzzy Model
In this study, an adaptive T-S fuzzy model is utilized to represent the nonlinear PMSG-HTs. Besides the model approximation effect, the T-S fuzzy model can also linearize the nonlinear PMSM system. Therefore, the future states of the system can be predicted, i.e., the stator flux and the electromagnetic torque. Here, the T-S fuzzy rules are defined as where R i (i = 1, . . . , l) represents the i th fuzzy rule, l refers to the number of rules, s 1 (k), . . . , s n (k) are the input variables of the T-S fuzzy model, M i j ( j = 1, . . . , n) are the input fuzzy sets, F i (·) is an arbitrary function, andx i (k) is the output in the THEN statement of R i . According to Equation (5), Equation (6), and Remark 3, the input of the T-S fuzzy model can be defined according to the , ω e (k − 1) denote the d-axis and q-axis currents and the electrical generator rotation speed, respectively. Thus, the aggregated output of the T-S fuzzy model isx where ε i (k) = To obtain the predictive states of the PMSG, it was important to compute the discrete-time model. Therefore,x(k) could be rewritten asx . Since the T-S fuzzy model is an adaptive fuzzy system, to establish the adaptive laws via tuning weights, the formulation of the cost function was defined as where e T−S (k) = x(k) −x(k) is the estimated error and P is the weight matrix. Consequently, by using the gradient descent algorithm, the weights could be adjusted, such that the cost function in Equation (14) was minimized. Thus, the tuning weights formula can be expressed as where γ is the learning rate and Θ(k) is the weighting matrix of the T-S fuzzy model. Thus, the future states of the PMSG could be obtained via the T-S fuzzy model.

Model Predictive Torque Control
This section focuses on developing the main control task in the MPTC scheme for sensorless PMSG-HTs. First, the torque reference and stator flux reference values based on the MPPT and MTPA Energies 2020, 13, 5296 9 of 18 schemes were estimated from Equation (8) and Equation (10), respectively. In addition, the d-axis and q-axis currents and the electrical generator rotation speed of the sensorless PMSG were used to obtain the predicted state sequence of its predecessor based on our T-S fuzzy model. Thus, the main characteristics of the MPTC scheme include using an adequate voltage vector selection strategy to indirectly control the stator flux and the electromagnetic torque simultaneously.
Based on the discrete-time T-S fuzzy model of the PMSG as defined in Equation (13), the predictive states can be described asX (k + 1) =Fx(k) +ĜU(k) and U(k) = [u(k), . . . , u(k + N c − 1)] T denote the predicted state vector and the control vector, respectively, while N p , N c are the output and control horizon at each time step k. Next, on the basis of the PMSM model, defined in Equation (4), the stator flux in the d-q reference frame within the k th sampling period could be obtained. Here, the prediction of the stator flux at the (k + 1) th sampling period via the root mean square (RMS) was derived as (17) while the prediction of the electromagnetic torque at the (k + 1) th sampling period could be expressed as whereî q andî d are the generator current vector in the d-q axis of our ATSFM scheme. (13) was used to obtain a prediction of the stator current at the next sampling instant from Equation (16). After obtainingî dq (k + 1), both the torque and the flux at the (k + 1) th instant could then be obtained according to Equations (4), (17), and (18).

Remark 4. the T-S fuzzy model of the PMSG defined in Equation
In this paper, pulse-width-modulation (PWM) voltage source converters with two levels (2L-VSCs) were applied to control each set of the three-phase windings of the generator. Each converter comprised of three legs, with the two switches in each leg operating in a complementary mode. Each of the two power switches were controlled via inverse logical signals, with "1" and "0" expressing the "on" and "off" states of the switch, respectively. In general, the cost function for PMSG drives is selected such that both the torque and the flux at the end of the cycle are as close as possible to the required value. Hence, the proposed finite time optimal control problem can be expressed as where λ ψ represents the weighting factor (here, the value was set to be λ T = T e /ψ s , where T e and ψ s are the rated value of torque and stator flux, respectively), U S (k) represents the switching state of phase S at the present cycle (where S is the switching vector, i.e., S = S a S b S c ), and V i is the inverter voltage of the PMSG system. The system also involved six different voltage vectors. All of the possible vectors from V 1 to V 6 for U S (k) at the current sampling instant could be easily evaluated such that Equation (19) was minimized. Finally, using the grid-side inverter to control the DC-link voltage v dc ensures that the maximum power will be produced. For clarity, the flowchart of the proposed ATSFMPC control scheme is shown in Figure 4. ), and i V is the inverter voltage of the PMSG system. The system also involved six different voltage vectors. All of the possible vectors from at the current sampling instant could be easily evaluated such that Equation (19) was minimized. Finally, using the grid-side inverter to control the DC-link voltage dc v ensures that the maximum power will be produced. For clarity, the flowchart of the proposed ATSFMPC control scheme is shown in Figure 4.
Torque and Stator Flux Prediction, i.e., Eqs. (17),  of the PMSG can further be extracted. Thus, the above reference torque and stator flux of PMSG driver will be as the tracking target for the proposed ATSFMPC strategy. Next, an adaptive Takagi-Sugeno (T-S) fuzzy model is used to express the highly nonlinear behavior of the larger, nonlinear sensorless PMSG-based hydro turbine system using input/output data of sensorless PMSG-based hydro turbine. Based on the T-S fuzzy model, the predicted state sequence of PMSG can then be obtained. Finally, the cost function for PMSG drives is selected, such that both the torque and the flux at the end of the cycle are as close as possible to the reference value. By minimizing the cost, the adequate voltage vector in the voltage-source inverter can then be determined in order to ensure the maximum power will be produced. First, to obtain the optimal reference torque and stator flux of PMSG driver at maximum available power of hydropower system, a tidal current meter is used to measure the tidal current speed, and then the maximum power point tracking (MPPT) is employed to estimate the mechanical reference rotation

Simulation Results and Discussion
speed ω re f m when the hydrokinetic turbine runs at the optimal blade tip speed ratio, i.e., perform on the maximum power operating point. According to Remark 3 and PI controller, the reference torque T re f e can be derived. Using the MTPA approach, the reference stator flux ψ re f s of the PMSG can further be extracted. Thus, the above reference torque and stator flux of PMSG driver will be as the tracking target for the proposed ATSFMPC strategy. Next, an adaptive Takagi-Sugeno (T-S) fuzzy model is used to express the highly nonlinear behavior of the larger, nonlinear sensorless PMSG-based hydro turbine system using input/output data of sensorless PMSG-based hydro turbine. Based on the T-S fuzzy model, the predicted state sequence of PMSG can then be obtained. Finally, the cost function for PMSG drives is selected, such that both the torque and the flux at the end of the cycle are as close as possible to the reference value. By minimizing the cost, the adequate voltage vector in the voltage-source inverter can then be determined in order to ensure the maximum power will be produced.

Simulation Results and Discussion
To verify the performance of the proposed control scheme for PMSG-based hydrokinetic turbine systems (PMSG-HTs), two real-world tidal current profile cases were considered in terms of the northwest European shelf [37] and the Pentland Firth [38]. In addition, considering the future practicability of the hydropower system, the classical PI current control scheme [13][14][15] is, therefore, utilized to compare with the proposed ATSFMPC approach for two realistic tidal current profiles in this paper. Taking a PMSG-HTs as an example, the detailed parameters are provided in Table 1. First, a realistic tidal current profile over the northwest European shelf area [37] was considered to demonstrate the effectiveness of the proposed control strategy, wherein the tidal speed profile can be expressed as follows: where t 1r and t 2r represent the initial time and the terminal time, respectively. The northwest European shelf tidal current speed profile is illustrated in Figure 5a. In addition, to maintain the maximum output power of the PMSG-HTs, the MPPT strategy was utilized to derive the mechanical reference rotation speed ω

Case 1. Northwest European Shelf Tidal Current Speed Profile
First, a realistic tidal current profile over the northwest European shelf area [37] was considered to demonstrate the effectiveness of the proposed control strategy, wherein the tidal speed profile can be expressed as follows: 12 12 6sin , and 30 where 1r t and 2r t represent the initial time and the terminal time, respectively. The northwest European shelf tidal current speed profile is illustrated in Figure 5a. In addition, to maintain the maximum output power of the PMSG-HTs, the MPPT strategy was utilized to derive the mechanical reference rotation speed ref m  , as shown in Figure 5b.  Figure 6 shows the simulation results of the d-axis and q-axis rotor currents for the proposed controller and the PI controller. Here, it is clear that the ATSFMPC strategy retained a more stable rotor current without the effect of the disturbance from the regular tidal currents.  Figure 6 shows the simulation results of the d-axis and q-axis rotor currents for the proposed controller and the PI controller. Here, it is clear that the ATSFMPC strategy retained a more stable rotor current without the effect of the disturbance from the regular tidal currents.
(a) (b) Figure 6. Response of d-axis and q-axis rotor currents for proportional-integral (PI) and ATSFMPC for Northwest European shelf.
Meanwhile, Figure 7 illustrates the performance comparison of the proposed ATSFMPC strategy and the conventional PI controller. To extract the maximum hydropower, the power coefficient should be maintained at its maximum value. Figure 7a shows that the maximum power coefficient max p C was largely achieved by the proposed controller despite the regular tidal current. Figure 7b shows that simulation compared the DC-link output voltage for the conventional PI controller and the proposed ATSFMPC scheme at the 200 V operation mode. As can be observed, voltage fluctuation of the proposed ATSFMPC scheme was smaller than PI controller when there were random fluctuations due to the uncertainty of tidal current speed, meaning the proposed ATSFMPC scheme can contribute to maintaining a more stable DC-link output voltage of the grid to increase the electric quality. Figure 7c,d shows the output voltages in terms of the d-q axis and the turbine power, respectively. Here, it is clear that the proposed ATSFMPC scheme provided a more stable stator voltage than the PI control method. Meanwhile, Figure 8d shows the hydrokinetic turbine power generated online from the variation current and voltage. Furthermore, obviously, as seen in Figures  5 and 6, the proposed ATSFMPC scheme has faster response and good performance when tidal current speed is changed instantly (sixth cycle) because its prediction optimal output derives future possible control input via evaluation of future resulting states. Conversely, the results indicate that the proposed ATSFMPC strategy demonstrated a superior performance and provided a more stable output to the grid. Moreover, Table 2 shows that the proposed ATSFMPC scheme can achieve closer to the 200 V operation mode for Northwest European shelf, and it has a smaller rms relative error of DC-link voltage values, and, further, it will effectively enable more efficient and stable hydrokinetic energy production. Meanwhile, Figure 7 illustrates the performance comparison of the proposed ATSFMPC strategy and the conventional PI controller. To extract the maximum hydropower, the power coefficient should be maintained at its maximum value. Figure 7a shows that the maximum power coefficient C max p was largely achieved by the proposed controller despite the regular tidal current. Figure 7b shows that simulation compared the DC-link output voltage for the conventional PI controller and the proposed ATSFMPC scheme at the 200 V operation mode. As can be observed, voltage fluctuation of the proposed ATSFMPC scheme was smaller than PI controller when there were random fluctuations due to the uncertainty of tidal current speed, meaning the proposed ATSFMPC scheme can contribute to maintaining a more stable DC-link output voltage of the grid to increase the electric quality. Figure 7c,d shows the output voltages in terms of the d-q axis and the turbine power, respectively. Here, it is clear that the proposed ATSFMPC scheme provided a more stable stator voltage than the PI control method. Meanwhile, Figure 8d shows the hydrokinetic turbine power generated online from the variation current and voltage. Furthermore, obviously, as seen in Figures 5 and 6, the proposed ATSFMPC scheme has faster response and good performance when tidal current speed is changed instantly (sixth cycle) because its prediction optimal output derives future possible control input via evaluation of future resulting states. Conversely, the results indicate that the proposed ATSFMPC strategy demonstrated a superior performance and provided a more stable output to the grid. Moreover, Table 2 shows that the proposed ATSFMPC scheme can achieve closer to the 200 V operation mode for Northwest European shelf, and it has a smaller rms relative error of DC-link voltage values, and, further, it will effectively enable more efficient and stable hydrokinetic energy production.

Case 2. Pentland Firth Tidal Current Speed Profile
In this case, the following realistic Pentland Firth tidal current speed profile [38] is defined as: where t 1r , t 2r , t 3r , and t 4r represent the initial, interval, and terminal times, respectively. The Pentland Firth tidal current speed profile is illustrated in Figure 8a. The mechanical reference rotation speed ω re f m was obtained via the MPPT strategy, as shown in Figure 8b. t represent the initial, interval, and terminal times, respectively. The Pentland Firth tidal current speed profile is illustrated in Figure 8a Figure 9 shows the response comparison between the proposed ATSFMPC strategy and the conventional PI controller. As the figure shows, the generated d-axis and q-axis currents for the hydro turbine continuously followed the tidal current speed pattern. In addition, the proposed ATSFMPC strategy largely provided a more stable output than the PI control method, while the maximum torque/ampere was achieved via the MTPA operation and the q-axis directly tracked the reference electromagnetic torque. Meanwhile, Figure 10 shows the simulation comparison of the proposed ATSFMPC strategy and the conventional PI controller. As shown in Figure 10a,b, the power coefficient and the DC-link voltage in the proposed control scheme exhibited less oscillations than in the conventional method. This means that the proposed control scheme provided a more efficient power generation, which was also maintained at maximum power.  Figure 9 shows the response comparison between the proposed ATSFMPC strategy and the conventional PI controller. As the figure shows, the generated d-axis and q-axis currents for the hydro turbine continuously followed the tidal current speed pattern. In addition, the proposed ATSFMPC strategy largely provided a more stable output than the PI control method, while the maximum torque/ampere was achieved via the MTPA operation and the q-axis directly tracked the reference electromagnetic torque. Meanwhile, Figure 10 shows the simulation comparison of the proposed ATSFMPC strategy and the conventional PI controller. As shown in Figure 10a,b, the power coefficient and the DC-link voltage in the proposed control scheme exhibited less oscillations than in the conventional method. This means that the proposed control scheme provided a more efficient power generation, which was also maintained at maximum power.  Moreover, Figure 10c,d shows the d-q axis output voltages and the turbine power, respectively. The results indicate that the proposed strategy demonstrated a superior performance and provided a more stable output to the grid. Finally, Table 3 shows that the proposed ATSFMPC scheme can Moreover, Figure 10c,d shows the d-q axis output voltages and the turbine power, respectively. The results indicate that the proposed strategy demonstrated a superior performance and provided a more stable output to the grid. Finally, Table 3 shows that the proposed ATSFMPC scheme can achieve closer to the 200 V operation mode for two Pentland Firth tidal currents, and it has a smaller relative RMS-error of DC-link voltage values, and, further, it will effectively enable a more efficient and stable hydrokinetic energy production. In other words, the proposed ATSFMPC approach can be applied to effectively suppress random fluctuations of the tidal current for supplying the stable electrical power with the hydro turbine at all times.

Conclusions
In this paper, an adaptive T-S fuzzy model predictive control (ATSFMPC) strategy has been proposed for a sensorless PMSG-based hydrokinetic turbine systems (PMSG-HTs) when there are random fluctuations due to the uncertainty of tidal current speed. To estimate the optimal reference torque and stator flux of the PMSG driver, an ocean current meter is first employed to measure the tidal current speed, and then the maximum power point tracking (MPPT) scheme is applied to predict and track the maximum power point in order to maximize the output power. Meanwhile, the maximum torque per ampere (MTPA) approach is used to derive the optimal reference torque and stator flux as a tracking target for the following proposed ATSFMPC strategy. Moreover, an adaptive T-S fuzzy model is implemented to deal with the nonlinearity of a PMSG-HTs. Thus, the nonlinear system can be represented by several linear sub-systems. Based on the T-S fuzzy model, the model predictive torque control (MPTC) strategy is then utilized to indirectly control the electromagnetic torque and stator flux by selecting an adequate voltage vector selection strategy, which not only improves the control performance, but also achieves maximum hydro power tracking. Eventually, two real-world tidal current speed profiles are utilized to verify the performance of the proposed control scheme and compared with the classical proportional-integral (PI) control strategy. According to the numerical simulation results, the proposed ATSFMPC scheme represents an improved performance, compared with the PI controller. The ATSFMPC strategy fully meets the desired objectives of ensuring stable tidal current conversion efficiency and output power of hydrokinetic turbine. We believe that our approach not only can establish technology in developing hydrokinetic power generation for Taiwan, but also holds great promise for solving a variety of practical engineering problems and challenges.