# Hydrodynamics and Tidal Turbine Generator Stability Analysis in Several Wave Variations

Authors: Muhammad Ikhwan; Syamsul Rizal; Marwan Ramli; Zainal A. Muchlisin; Said Munzir
Journal: International Journal of Differential Equations
 (2021)

---

## Abstract

The development of tidal turbines continues to be carried out by many researchers, including the incorporation of a control system for optimization purposes. This paper attempts to assess the stability of two mechanical systems in a tidal turbine: a propeller harvesting kinetic energy and a d-q diagram system on a permanent-magnet synchronous generator (PMSG). The method employed is the representation of a phase plane profile with a stable eigenvalue. The critical value of the turbine’s rotations per minute provides some points of equilibrium. The effect of the angular velocity singular on the modified system was also investigated. There is no cutoff control for the generator rotational speed in the case of weak currents, according to the results. The combination of the three tidal turbine components results in a shift in the equilibrium point. Although PMSG has an infinite equilibrium point along the lineId = 0, the effect of the rotor angular velocity prevents all of these points from being in equilibrium. Finally, in this study, the rotor angular velocity caused by the speed and type of ocean currents are only the upper and lower limits. The stability of the various wave variations is within this range.


---

## Body

## 1. Introduction

This research was driven by a hydrodynamic laboratory’s need to produce various types of waves in the wave container, including extreme waves. These waves will be used to test an oceanic object before it is operated, either in the middle of the ocean or near the land. Wind pressure, moving objects in water, and large-scale objects such as celestial bodies’ gravity all contribute to the formation of waves. The ocean’s kinetic energy is the most promising renewable energy as a substitution for the use of fossil energy. One of the energies that come from the sea is tidal energy, which has received many researchers’ attention because of its regularity and certainty [1, 2]. In both short and long term, tidal energy can be predicted with high accuracy. As a result, the power output of a specific tidal generator at a specific location can also be accurately estimated [3]. The advantages of tidal currents compared to wind resources are that seawater is denser than air and the ocean currents produce a more predictable resource than wind resource [4, 5].The tidal turbine system uses a horizontal axis tidal current turbine system. This type of turbine has a better working efficiency than a turbine that rotates vertically at a performance coefficient ofλ = 5 [6]. A horizontal tidal turbine captures hydrodynamic energy and converts it to electricity using a generator. The generator commonly used is PMSG which has a weakness in the vulnerability to damage caused by the absence of a protection system between the turbine and generator [7]. Another study states that PMSG protection can be adopted from competing generators by adding a gearbox system [8]. In the development of PMSG, the control system at the propeller pitch angle is the best step to control the tidal turbine [9, 10]. Pitch angle can assist to decrease a variety of errors, including turbine laying errors that do not match the tidal current direction [11].Sangiuliano and Mastrantonis [12] investigated the benefits of tidal energy in Scotland. Offshore renewable energy has helped the transition from fossil energy to environmental-friendly energy. These drivers have prompted Scotland in developing and implementing the first ocean plan for consuming energy form tidal currents. Tidal turbine research has been carried out on the southwest coast of Yemen using a tidal current turbine (TCT) that uses ocean currents. Towers are used to harvest wind potential energy at maximum resource simultaneously. The combination structured between the wind turbine and tidal turbine is a cost-effective way that harvested energy [13]. Eme et al. [14] showed that the power product rate is immediately proportional to a hydrokinetic turbine’s current speed. Therefore, the installation cost was drastically reduced by about $ 5,400 per kW. It is easy to assess the lack of mechanical and engine technology familiar to most turbine activists. Research examining technical efficiency has also been conducted by Wang and Wang [15]. The study compared tidal turbines, their advantages and disadvantages under various conditions, and several future modifications.The installation of the existing tidal-stream turbines is carried out in energetic locations with flow rates of more than 2 m/s. Turbines operating in Mexico and the Philippines have current speeds below 1.5 m/s [16]. Another study stated that a constant average speed of 0.81 m/s could be used to install tidal turbine [17]. Tidal turbines were simulated along the Moroccan coast (550 kM), and Nachtane et al. [18] recommend the three most strategic locations for tidal energy harvesting. The technical problems associated with the hydrodynamic design of the rotor blades due to marine turbine design are some of the techniques used for wind turbines. Arribas [4] describes a rotor design’s projected hydrodynamics to extract maximum power from the tides. The theory of blade momentum is used for rotor modelling to overcome various aspects and limitations when applied to marine rotor designs. The use of tidal turbines can erode the seabed where the power plant is installed. Sun et al. [19] simulated the scouring effect of turbine-induced currents, followed by validation of the model through experimental data at an early stage. Simulation scale models were expanded to large-scale models for turbine scour estimation. However, the simulation misses the estimated about 15–20% based on field-based data [20].Tian et al. [21] reviewed and tested the prototype in a tower chamber. The field experimental data expressed the TCT is limited by a maximum speed of 0.25 m/s and had an efficiency of about 0.33. Fluid dynamics is based on the Reynolds Averaged Navier Stokes equation, and we carried out the computational simulation of the generator performances in both unstable and stable conditions. These models were mostly utilized in the distribution of polluted matter and flood calculations, but the model could be modified for waves and current simulation. The article studied the available one-to-three-dimensional computer application to analyze their validity to examine the TCT’s energy that is free in ocean [22].Li et al. [8] use a full-scale tidal current turbine to verify the design method with high reliability, and a TCT of 600 kW was used in that study. The turbine consists of a two-propeller rotor, a low rate gearbox, and a generator. Within the logical analysis and test results, satisfactory verification is achieved. Small-power turbines are frequently used on a laboratory scale with a generation method similar to tidal energy generation. Tides are often viewed from the gravity of celestial bodies against the Earth and wind [23–25]. The increase in wave amplitude is greatly influenced by the composition of the wave and its propagation as it moves from the middle of the sea to the land [26–28]. As a consequence, the observation of the turbine location takes into account the current strength, as indicated by the water ripples on the surface, as well as the internal current strength [29, 30].


## 2. Materials and Methods

There are many possible tidal current turbine arrangements. The planner can decide between straight impulse systems or gear-equipped systems. The two models of TCT are synchronous generator systems and indirect generator systems. The network’s connection can be via a full-power converter, partial rated, or directly connected converter. The use of PMSG in a direct blade and generator is an engaging and frequently used vehicle by researchers. When used together with abundant power, the dynamo is separated from the network/grid, and a fully unsteady velocity process can be performed. Figure1 gives the configuration regarded in the paper.Figure 1
PMSG-based TCT [8].The system consists of a generator and maximum power converter, which consists of a back-to-back voltage source converter. This topology allows a two-way flow of current, and the others have completely regulated above the generator twist. An alternative is to practice a diode rectifier over PMSG, simplifying command, and reducing costs. The difficulties of using rectifier are often associated with power. Diode rectifiers have been reported to offer large harmonic distortion to the dynamo that affect efficiency and can produce harmonic torque.The hydrodynamic torque developed by the turbine rotor in Figure1 is [31]:(1)τh=0.5ρAv3Cpβ,λωr.For modelling purposes,Cp is power extraction efficiency. Tip speed ratio (TSR) λ is determined based on current and rotor velocity. The first conditions are adjusted so that the rotor velocity is not zero. The estimated λ and β are later drafted for relevant values using Cp listed in the table. The lookup table is for a generic 1 MW TST rotor. After obtaining the values, the hydrodynamic torque is calculated. The model only turns the incident currents into a hydrodynamic twist and does not include dynamic flow effects, as seabed and tower influences.The hydrodynamic twist caused by the rotor revolution is assigned for the dynamo via the impulse train. A direct-drive system is used as a configuration indicating that the rotor is straightly related to a single bar that ends to a dynamo. For TCT modelling, a doubled mass system is governed in which both the bar and rotor are deemed whole combined mass, and the dynamo is the examined mass held together. The couple masses are attached within an adjustable bar with an inflexibilityKdt and an attenuation Ddt. One point of the bar is pushed by a turbine rotor, providing torque τh. The other point of the bar is saddled by the dynamo, which generates a twist τg. The following system of equations describes the drive train:(2)τh=Jrd2θrdt2+Kdtθr−θg+Ddtdθrdt−dθgdt,τg=Jgd2θgdt2+Kdtθg−θr+Ddtdθgdt−dθrdt,dθrdt=ωr,dθgdt=ωg,ωg=Ppωr,where θr and θg are the angular positions of the rotor and generator bar, ωr, and ωg are the speed of rotor and generator bar, and Jr and Jg are the moment of inertia for the rotor and generator, respectively. Pp is the number of pairing poles that made ωg > ωr. The electrical equalization used to the generator model is pictured in a rotating source framework. The d-axis is adapted simultaneously with the position of the rotor vector, and the q-axis affects the d-axis by a quarter circle degree. In this way, representing the electrical equation reduces the twist composing and flux magnetizing elements of the diesel. Therefore, it supports the improvement of a regulated plan that accommodates the self-governing regulator of twist.Also, it is an area where surface vertical TCTs are hit with maximum average flow during flood and ebb cycles [32]. Turbulence is mostly insignificant when the TCT responds in regular harmonic waves. Martinez et al. [33] have shown that the usual flow thrust for the two turbine models and the same flow conditions is about twice as high as the thrust induced by turbulence. The turbulence effect is neglected for analysis purposes. With this approach, the turbine’s power P in waves can be explained by the density-area-speed equation [17]:(3)P=0.5ρAv3Cpλ,β,where ρ express the fluid density of sea water, area A = πr2 follows the circle cross section with the radius of the turbine r, v is the current speed, Cp is the turbine efficiency, β represents the blade pitch angle, and λ is the velocity ratio at the defined propeller tip shown as follows [34]:(4)λ=ωrRsv,where ωr is the rotor angular velocity.The TCT system uses a motor rotation system, gearbox, and generator. This system is structured with a control system structure with the assumption that combining three elements in one unit, the equation can be formulated as in [35]. The gearbox has a significant impact on determining the angular velocity value of each shaft. Equation (2) applies linearly in the absence of the gearbox [36]. However, if a gearbox is added, equation (2) is modified by viscous friction to the following equation:(5)Jeqω˙r=τg−τhg−Vfωr,where Jeq is the inertia equivalency of Jr and Jg, where Jeq = Jg + Jr/ng2 and ng is the gear ratio. From this structure, Vf is the friction of the viscous coefficient. Torque τhg is obtained from τh/ng, and τg is derived from the d-q diagram of the generator as follows [37]:(6)LdI˙d=ud−RsId+LqωgIq,(7)LqI˙q=uq−LdωgId−RsIq−ϕmωg,(8)τg=32Ppϕm+Ld−LqIdIq,where Ld and Lq are the direct inductance and its quadrature stator and ud and uq are the related voltage with the inductances. Id and Iq represent the stator currents, and φm is the permanent magnet flux. The hydrodynamic and generator parameters can be seen in Table 1.Table 1
Hydrodynamic and generator parameters [35].
Parameter symbolValueParameter symbolValueΡ1025 psuφm6.8973 VsA355 m2Vf0.3v2.4 m/sJeq100 kg m2Cp0.39Rs4.61 mΩPp60 pairsL886.48μH


## 3. Results and Discussion

### 3.1. Model Linearization and Parameter Simplification

By paying attention to the relationship betweenωr and ωg in equation (2), equations (6) and (7) are modified to ωr. The aim is to make it easier to get the τg value. The generator rotation is adjusted optimally by reducing its angular velocity. The implication is that a hydrodynamic torque regulator is required based on the speed of the propeller. Equation (5) can be obtained optimally by solving the following series of differential equations:(9)dωrdt=1Jeqτg−τhg−Vfωr,dIddt=udL−RsLId+PpωrIq,dIqdt=uqL−PpωrId−RsLIq−PpLϕmωr,where there is a simplification of the parameter Ld = Lq = L.Furthermore, model (9) is linearized with the Jacobian matrix J as follows:(10)J=ρAv3Cp2Jeqngωr2−VfJeq0−3Ppϕm2JeqPpIq−RsLPpωr−PpId−PpLϕm−Ppωr−Pp.MatrixJ is used as the linear form of equation (9) and is useful for obtaining eigenvalues.


### 3.2. Equilibrium and Eigenvalues

The equilibrium point is obtained from the condition of the first derivative of each state variable with a value of zero. It means that the condition is on a straight gradient and is considered stable.(11)1Jeqτg−τhg−Vfωr=0,udL−RsLId+PpωrIq=0,uqL−PpωrId−RsLIq−Pp=0.The equilibrium point is proposed at the limit conditionsωr of 13 rpm or 81.68 rad/s. Equation (8) shows that the rotor flux is constant and the d-axis current component is always equal to zero. As a result, the electric torque is proportional to the current along the q-axis. Solution Id and Iq on equation (11) is at the equilibrium point (Id, Iq) = (0, 0). Some points, the proposed theory and equilibrium points, have stability that is checked by the eigenvalue. The eigenvalues formed can be seen in Table 2.Table 2
The eigenvalues at some equilibrium points.
Test point (ωr, Id, Iq)Eigen real numberEigen imaginary numberStability note(−81.68, 0, 0)0.01620Not stable−5.2003 × 10−3−4.9009 × 103Not equilibrium−5.2003 × 10−34.9009 × 103Extreme point(0, 0, 0)——Not defined by division by zero(81.68, 0, 0)0.01620Not stable−5.2003 × 10−3−4.9009 × 103Not equilibrium−5.2003 × 10−34.9009 × 103Extreme point(1.9777 × 10−10 − 8.6672 × 10−5i, 2.3766 × 105 − 1.9116 × 10−13i, 5.4230 × 10−1 + 2.3766 × 105i)−0.00520Asymptotically stable−0.01040Equilibrium point−1.7021 × 10107.7680 × 104(1.9777 × 10−−10 + 8.6672 × 10−5i, 2.3766 × 105 + 1.9116 × 10−13i, 5.4230 × 10−1 − 2.3766 × 105i)−0.00520Asymptotically stable−0.01040Equilibrium point−1.7021 × 1010−7.7680 × 104(−1.9777 × 10−10 + 2.0645 × 102i, −0.0078, −1.91161 × 10−13 + 9.4334 × 10−8i)1.2387 × 1041.8712 × 10−8Not stable−0.00600Equilibrium point−1.2387 × 104−1.8712 × 10−8(−1.9777 × 10−10 − 2.0645 × 102i, −0.0078, −1.91161 × 10−13 − 9.4334 × 10−8i)1.2387 × 104−1.8712 × 10−8Not stable−0.00600Equilibrium point−1.2387 × 1041.8712 × 10−8The eigenvalues correspond toωr, Id, and Iq. The three eigenvalues at the negative pole mean that this model is stable around the equilibrium point. When the state variable is reduced to Id and Iq, all points near Id = 0 have negative eigenvalues. The system in equation (11) is rarely used in research involving the combination of rotors, gearboxes, and generators. This is because obtaining a stable point is difficult. Negative eigenvalues are asymptotically stable for all pairs of Iq with Id = 0. However, because of the system’s influence of ωr, only a very large position Id can produce a stable equilibrium point. This condition is outside the range of Id movements, so this equilibrium point has never been reached in practice.The rotation speed of the generator determines the amount of energy produced in the generator system. Several studies have shown that the success of the TCT in harvesting tidal energy has been widely proven in theory, simulation, and in the field. The energy under review is only limited to the propeller until it reaches the generator side converter. Meanwhile, the grid side converter has other problems that are not reviewed in this article. Based on the review results, the turbines are mostly placed in narrow straits or canals, and in some cases, they are placed in shallow waters. The desired power efficiency is in the range of 0.33–0.39, with a speed limit of less than 2 m/s.The phase plane is away from the lineId = 0 and runs well in its surroundings. The relationship between Id and ωr cannot be seen because the phase plane is not formed, while Iq and ωr result in a separate pattern by limit ωr = [−81.68, 81.68]. The value of ωr depends on Iq due to the simplification of the parameter L which causes the Id value to stick to equation (8) (see Figure 2).Figure 2
The stability of the generator based on the stator current d-q; the black quiver indicates the direction of change inId–Iq, and the color shade indicates ωr.Colour change in Figure2 describes the ωr range from the lower to the upper limit. The PMSG used is limited to ωr = [−81.68, 81.68], but the figure shows possible out-of-range values. The green area has a low angular velocity and is located on the line Iq = 0. Changes in stator current continue to grow away from the point Id = 0, indicating that this is the optimal area for producing constant energy. Outside of the green zone, we refer to it as the cutoff zone or extreme.The cases of high and extreme waves with extraordinary strength cannot be explained on the stability curve. However, the power more significant than the 13 rpm rotation of the turbine can be described if the system cutoff is not applied in real conditions. Extreme waves are waves with great strength and are a combination of waves that have the same phase and reinforce each other [27]. If the system has no cutoff, this wave can damage the system, when ωr < −81.68 and ωr > 81.68. The twist on the curve becomes very large from the steady-state. Figure 3 shows the curve movement is no longer the same if ωr is very large. The electric current strength of the generator in full negative and positive conditions is the opposite of each other.Figure 3
Hydrodynamic simulation at low-velocity current (2.4 m/s): (a) hydrodynamic torqueτh, (b) generator torque τg, (c) d-axis stator current id, (d) q-axis stator current iq, and (e) rotor angular speed ωr.
(a)(b)(c)(d)(e)The stator current changes that occur in the cutoff zone are more dynamic. For extracting energy, only a limited range is provided. Moving slightly to theId > 0 area, the change immediately enlarges and tends to flatten out in very large ωr. That is, the effect of a change in Id is greater than the effect of a change in Iq, resulting in contradictory events when compared to equation (8), which is only influenced by Iq [36]. According to these findings (Figure 2 and Table 2), the Id boundary should be updated, so that the equilibrium region is near the line Id = 0 or in the vicinity of the singular value ωr.The simulation in the green area is shown using the parameters in Table1. The simulation is run without the use of any input interference, such as rotor and generator speed control. The blade-side torque change appears normal in Figures 3(a) and 3(b), but the torque generator remains small. Since friction has no effect on the torque on the blade, the value remains very high. The value of viscosity friction occurs when τh is converted to τg, and even though the gearbox has increased the value of τg, the value of τg will remain small. As a result, the current can expand to greater than 2.4 m/s, which is known as the minimum current velocity when hitting the blade in other studies. Figures 3(c) and 3(d) show that Id and Iq remain stable. As previously explained in Figure 2, Iq can be increased around Iq = 0. In the same case, Iq is increased and controlled by the blade pitch angle, so that Iq approaches −2, the rotor angular speed approaches 13 rpm, and the power reaches 1 MW [35].


## 3.1. Model Linearization and Parameter Simplification

By paying attention to the relationship betweenωr and ωg in equation (2), equations (6) and (7) are modified to ωr. The aim is to make it easier to get the τg value. The generator rotation is adjusted optimally by reducing its angular velocity. The implication is that a hydrodynamic torque regulator is required based on the speed of the propeller. Equation (5) can be obtained optimally by solving the following series of differential equations:(9)dωrdt=1Jeqτg−τhg−Vfωr,dIddt=udL−RsLId+PpωrIq,dIqdt=uqL−PpωrId−RsLIq−PpLϕmωr,where there is a simplification of the parameter Ld = Lq = L.Furthermore, model (9) is linearized with the Jacobian matrix J as follows:(10)J=ρAv3Cp2Jeqngωr2−VfJeq0−3Ppϕm2JeqPpIq−RsLPpωr−PpId−PpLϕm−Ppωr−Pp.MatrixJ is used as the linear form of equation (9) and is useful for obtaining eigenvalues.


## 3.2. Equilibrium and Eigenvalues

The equilibrium point is obtained from the condition of the first derivative of each state variable with a value of zero. It means that the condition is on a straight gradient and is considered stable.(11)1Jeqτg−τhg−Vfωr=0,udL−RsLId+PpωrIq=0,uqL−PpωrId−RsLIq−Pp=0.The equilibrium point is proposed at the limit conditionsωr of 13 rpm or 81.68 rad/s. Equation (8) shows that the rotor flux is constant and the d-axis current component is always equal to zero. As a result, the electric torque is proportional to the current along the q-axis. Solution Id and Iq on equation (11) is at the equilibrium point (Id, Iq) = (0, 0). Some points, the proposed theory and equilibrium points, have stability that is checked by the eigenvalue. The eigenvalues formed can be seen in Table 2.Table 2
The eigenvalues at some equilibrium points.
Test point (ωr, Id, Iq)Eigen real numberEigen imaginary numberStability note(−81.68, 0, 0)0.01620Not stable−5.2003 × 10−3−4.9009 × 103Not equilibrium−5.2003 × 10−34.9009 × 103Extreme point(0, 0, 0)——Not defined by division by zero(81.68, 0, 0)0.01620Not stable−5.2003 × 10−3−4.9009 × 103Not equilibrium−5.2003 × 10−34.9009 × 103Extreme point(1.9777 × 10−10 − 8.6672 × 10−5i, 2.3766 × 105 − 1.9116 × 10−13i, 5.4230 × 10−1 + 2.3766 × 105i)−0.00520Asymptotically stable−0.01040Equilibrium point−1.7021 × 10107.7680 × 104(1.9777 × 10−−10 + 8.6672 × 10−5i, 2.3766 × 105 + 1.9116 × 10−13i, 5.4230 × 10−1 − 2.3766 × 105i)−0.00520Asymptotically stable−0.01040Equilibrium point−1.7021 × 1010−7.7680 × 104(−1.9777 × 10−10 + 2.0645 × 102i, −0.0078, −1.91161 × 10−13 + 9.4334 × 10−8i)1.2387 × 1041.8712 × 10−8Not stable−0.00600Equilibrium point−1.2387 × 104−1.8712 × 10−8(−1.9777 × 10−10 − 2.0645 × 102i, −0.0078, −1.91161 × 10−13 − 9.4334 × 10−8i)1.2387 × 104−1.8712 × 10−8Not stable−0.00600Equilibrium point−1.2387 × 1041.8712 × 10−8The eigenvalues correspond toωr, Id, and Iq. The three eigenvalues at the negative pole mean that this model is stable around the equilibrium point. When the state variable is reduced to Id and Iq, all points near Id = 0 have negative eigenvalues. The system in equation (11) is rarely used in research involving the combination of rotors, gearboxes, and generators. This is because obtaining a stable point is difficult. Negative eigenvalues are asymptotically stable for all pairs of Iq with Id = 0. However, because of the system’s influence of ωr, only a very large position Id can produce a stable equilibrium point. This condition is outside the range of Id movements, so this equilibrium point has never been reached in practice.The rotation speed of the generator determines the amount of energy produced in the generator system. Several studies have shown that the success of the TCT in harvesting tidal energy has been widely proven in theory, simulation, and in the field. The energy under review is only limited to the propeller until it reaches the generator side converter. Meanwhile, the grid side converter has other problems that are not reviewed in this article. Based on the review results, the turbines are mostly placed in narrow straits or canals, and in some cases, they are placed in shallow waters. The desired power efficiency is in the range of 0.33–0.39, with a speed limit of less than 2 m/s.The phase plane is away from the lineId = 0 and runs well in its surroundings. The relationship between Id and ωr cannot be seen because the phase plane is not formed, while Iq and ωr result in a separate pattern by limit ωr = [−81.68, 81.68]. The value of ωr depends on Iq due to the simplification of the parameter L which causes the Id value to stick to equation (8) (see Figure 2).Figure 2
The stability of the generator based on the stator current d-q; the black quiver indicates the direction of change inId–Iq, and the color shade indicates ωr.Colour change in Figure2 describes the ωr range from the lower to the upper limit. The PMSG used is limited to ωr = [−81.68, 81.68], but the figure shows possible out-of-range values. The green area has a low angular velocity and is located on the line Iq = 0. Changes in stator current continue to grow away from the point Id = 0, indicating that this is the optimal area for producing constant energy. Outside of the green zone, we refer to it as the cutoff zone or extreme.The cases of high and extreme waves with extraordinary strength cannot be explained on the stability curve. However, the power more significant than the 13 rpm rotation of the turbine can be described if the system cutoff is not applied in real conditions. Extreme waves are waves with great strength and are a combination of waves that have the same phase and reinforce each other [27]. If the system has no cutoff, this wave can damage the system, when ωr < −81.68 and ωr > 81.68. The twist on the curve becomes very large from the steady-state. Figure 3 shows the curve movement is no longer the same if ωr is very large. The electric current strength of the generator in full negative and positive conditions is the opposite of each other.Figure 3
Hydrodynamic simulation at low-velocity current (2.4 m/s): (a) hydrodynamic torqueτh, (b) generator torque τg, (c) d-axis stator current id, (d) q-axis stator current iq, and (e) rotor angular speed ωr.
(a)(b)(c)(d)(e)The stator current changes that occur in the cutoff zone are more dynamic. For extracting energy, only a limited range is provided. Moving slightly to theId > 0 area, the change immediately enlarges and tends to flatten out in very large ωr. That is, the effect of a change in Id is greater than the effect of a change in Iq, resulting in contradictory events when compared to equation (8), which is only influenced by Iq [36]. According to these findings (Figure 2 and Table 2), the Id boundary should be updated, so that the equilibrium region is near the line Id = 0 or in the vicinity of the singular value ωr.The simulation in the green area is shown using the parameters in Table1. The simulation is run without the use of any input interference, such as rotor and generator speed control. The blade-side torque change appears normal in Figures 3(a) and 3(b), but the torque generator remains small. Since friction has no effect on the torque on the blade, the value remains very high. The value of viscosity friction occurs when τh is converted to τg, and even though the gearbox has increased the value of τg, the value of τg will remain small. As a result, the current can expand to greater than 2.4 m/s, which is known as the minimum current velocity when hitting the blade in other studies. Figures 3(c) and 3(d) show that Id and Iq remain stable. As previously explained in Figure 2, Iq can be increased around Iq = 0. In the same case, Iq is increased and controlled by the blade pitch angle, so that Iq approaches −2, the rotor angular speed approaches 13 rpm, and the power reaches 1 MW [35].


## 4. Conclusions

This study discovered a positive relationship between angular velocity and q-axis stator current in low-speed PMSG generators, whereas d-axis stator current had no effect on the other two state variables. The d-axis stator current provided a stable condition around its zero value. This condition changes when the rotor’s angular speed reaches the limit of −13 to 13 rpm, a condition in which the PMSG tries to remain stationary while rapidly moving away from the point of equilibrium. Although the eigenvalues at each pole are negative, the real value formed is very small and close to 0. When the rotor’s angular velocity is greater than −13 to 13 rpm, the d-axis tolerance of the stator current becomes smaller. Current changes are becoming larger and more volatile. The result may be detrimental to the turbine structure. The gearbox system also performs well; at low-current speed simulations, the torque on both sides of the gearbox can vary significantly, necessitating a minimum value of ocean current speed.



---

*Source: 6682407-2021-06-16.xml*