Speed Tracking for IFOC Induction Motor Speed Control Using Hybrid Sensorless Speed Estimator Based on Flux Error for Electric Vehicles Application

: This paper presents hybrid sensorless speed tracking by an indirect ﬁeld-oriented control (IFOC) for an induction motor (IM). The sensorless model is based on an improved virtual estimation topology model to predict the virtual speed and ﬂux of the IM using stator current components. The hybrid sensorless model, deﬁned as a modiﬁcation of voltage with a rotor ﬂux-oriented current model, was also implemented with proportional-integral (PI) control for comparison with the conventional voltage model (CVM). The suggested adaptive mechanism for PI control in the hybrid estimator was able to compensate for the back-EMF error from the rotor ﬂux-oriented current model into the voltage model and change the air gap ﬂux of the IM. An accurate rotor ﬂux position was estimated and used to estimate the speed with low speed error. This IFOC model, with various speed change references, was tested in a simulation environment by using the MATLAB/Simulink program. The proposed hybrid estimator was tested in two different EV operations, which were reverse and forward operations. The effectiveness of the proposed estimator was analyzed for its transient and steady-state performances based on settling time, recovery time and the overshoot and speed error percentages. All the results were in good agreement in terms of the stability of the speed and current controller with minimum speed error obtained, where the average errors were 0.08% and 0.16% for high speed and lower speed, respectively.


Introduction
A high-performance induction motor (IM) requires various special characteristics, such as robustness against parameter variation, a stable system during operation, faster load disturbance rejection and also faster dynamic responses. The use of an IM in industrial applications is generally based on its main advantages, which are high reliability and a wide range of speed [1]. Due to its simple construction without magnets on its rotors, the trend of IM application will increase in the coming years because of its low cost.
The first electric vehicle (EV) was built in 1839, but the IM application in the transportation industry has been gaining acceptability in recent years [2]. EVs will be a major transportation mechanism due to the increase in concerns over global warming and due to energy efficiency plans that have led car manufacturers to deploy more IMs in EVs, compared with DC motors or permanent magnet synchronous motors (PMSMs). IMs in the EV industry are used over other types of motors because of the possibility of regenerative braking in the four-quadrant operation [3].
A key condition toward a wider application of the IM is vector control, which is mostly used to control the IM in advance with a sensorless mode. The essential feature of vector control is the ability to have decoupled control between the flux and the torque generated by the stator currents with direct and indirect field-oriented control methods, which reflect the motor's speed response [4]. Therefore, to realize the sensorless mode, an estimation technique could be used to estimate the speed of the IM, which would eliminate the use of a speed sensor and consequently reduce the maintenance cost, hence making the EV a smart system. Various estimation techniques have been developed to obtain the minimum speed error and reduce the computational burden in the estimating processes, such as the current model, voltage model, fuzzy logic and Kalman filter [5][6][7][8]. These techniques have excellent performance depending on how they are applied, but in some cases the speed estimation is still unsatisfactory due to the large speed error produced.
Vector control for motor operation is normally achieved by using a proportionalintegral (PI) controller either for the speed or current control. As discussed in [9][10][11], a PI controller is used to control the speed with a single estimator of speed and without any adaptation mechanism. However, the weakness of a PI controller is that it cannot compromise with unknown motor parameters, such as stator resistance and inductance, if they are not being considered in the input parameters, which currently are obtained from the speed input sensor. These parameters vary during operation due to temperature changes and the magnetic saturation effect, which deteriorate the IM performance and make it operate at low efficiency [12]. The issue was resolved in [13][14][15], where a tuning method (PI tuner) was used to tune the controller's response. The tuning process focused on the variable speed operation without considering the load disturbance rejection. This could harm the IM, where the cogging torque reduces the IM performance and makes it operate below the base speed [16]. Various tuning methods have been developed by introducing currently trending methods, such as fuzzy logic, Kalman filters, neural networks and neuro-fuzzy, in order to have good load disturbance rejection [17][18][19][20][21][22]. However, due to complexities, few have been applied in IM vector control. These complexities give a computational burden to the control system, which causes the overlapping loop problem between the controllers [23].
In this paper, an indirect field-oriented control (IFOC) was adopted with a sensorless model and with the PI speed and current controllers to control the speed and torque of the IM. Here, the sensorless model took the function of the speed sensor. The current and voltage coming to the IM became the inputs of the speed estimation model, which were based on the parameters of the IM. A pre-compensating method was introduced in the current controller to balance the control operation at a wide range of speed so that the ripples on the current response could be reduced. In order to resolve the load disturbance rejection, a hybrid topology based on the modified voltage model and the rotor fluxoriented current model was applied to estimate the feedback speed by the stator voltage and current components, which responded according to the load torque disturbances. The results obtained were analyzed for controller stability to seek the minimum speed error. This paper is organized as follows: In Sections 2 and 3, the dynamic modeling of the IM and the concept of field-oriented control (FOC) are described, respectively. In Section 4, the proposed hybrid estimator with the stator flux control as the adaptive mechanism is explained. In Section 5, the stability of the speed and current controllers are determined. In Section 6, the simulation results are analyzed, followed by the performance evaluation and conclusion.

Dynamic Modeling of Induction Motor
The space vector method presented in Figure 1 describes the mathematical modeling of the IM. Based on the vector method, s and r are stator and rotor circuits, respectively, and e refers to the synchronously rotating frame, while a, b and c are the axes of the phase system and d and q are direct and quadrature axes with respect to the stator components, respectively. The field angle θ e is extracted from the time integral of the angular frequency ω e . By subtracting the slip frequency ω sl from the field angle, the rotor flux position θ r can be calculated as: Machines 2022, 10, x FOR PEER REVIEW 3 of 21

Dynamic Modeling of Induction Motor
The space vector method presented in Figure 1 describes the mathematical modeling of the IM. Based on the vector method, s and r are stator and rotor circuits, respectively, and e refers to the synchronously rotating frame, while a , b and c are the axes of the phase system and d and q are direct and quadrature axes with respect to the stator components, respectively. The field angle e θ is extracted from the time integral of the angular frequency e ω . By subtracting the slip frequency sl ω from the field angle, the rotor flux position r θ can be calculated as: The state-space modeling of the IM is given in a synchronous frame, such that: where x , u and y are state vector, input vector and output vector, respectively. Meanwhile, A , B and C are the matrix system. t ds qs dr qr where i , v and ψ are current, voltage and flux, respectively, on the d q − axes with respect to their frames.  The state-space modeling of the IM is given in a synchronous frame, such that: where x, u and y are state vector, input vector and output vector, respectively. Meanwhile, A, B and C are the matrix system.
x = i ds i qs ψ dr ψ qr t , u = v ds v qs t , y = i ds i qs t where i, v and Ψ are current, voltage and flux, respectively, on the d-q axes with respect to their frames.
Meanwhile, T r is rotor time constant and σ is leakage coefficient, while R and L are resistance and inductance, respectively, with respect to their frames and L m is the magnetizing inductance. The statespace modeling of the IM creates ease in interpreting the stationary current and voltage transformation into the d-q axes. This was recreated based on the balanced condition represented by the space vector. Therefore, this model was necessary for modeling the sensorless model based on the voltage, current and motor parameters in order to create a virtual speed value.

Indirect Field-Oriented Control
Indirect field-oriented control is a vector control with high acceptability in a highperformance drive system. The principle is essentially based on the decoupling between flux and torque through controlling the stator current component, as shown in Figure 2. eling the sensorless model based on the voltage, current and motor parameters in order to create a virtual speed value.

Indirect Field-Oriented Control
Indirect field-oriented control is a vector control with high acceptability in a highperformance drive system. The principle is essentially based on the decoupling between flux and torque through controlling the stator current component, as shown in Figure 2. Based on the dynamic model of the IM, the three phase quantities are transformed into two coordinates on the d q − axes, which are the flux and torque components aligned on the q − axis and d − axis, respectively. To realize this condition, the rotor flux on the axis is set to be zero, while the d − axis reaches the nominal value of the magnetizing flux. The behavior of the drive system subjected to the rotor flux-oriented control is depicted in Figure 3. With zero flux on the q − axis, the electromagnetic torque reduces to: where p is the number of poles. The angular frequency of the rotor flux is given by: Based on the dynamic model of the IM, the three phase quantities are transformed into two coordinates on the d-q axes, which are the flux and torque components aligned on the q-axis and d-axis, respectively. To realize this condition, the rotor flux on the axis is set to be zero, while the d-axis reaches the nominal value of the magnetizing flux. The behavior of the drive system subjected to the rotor flux-oriented control is depicted in Figure 3.

Indirect Field-Oriented Control
Indirect field-oriented control is a vector control with high acceptability in a highperformance drive system. The principle is essentially based on the decoupling between flux and torque through controlling the stator current component, as shown in Figure 2. The behavior of the drive system subjected to the rotor flux-oriented control is depicted in Figure 3. With zero flux on the q − axis, the electromagnetic torque reduces to: where p is the number of poles. The angular frequency of the rotor flux is given by: With zero flux on the q-axis, the electromagnetic torque reduces to: where p is the number of poles. The angular frequency of the rotor flux is given by: where the slip frequency is obtained from the q-axis stator current, such that: The estimated rotor flux linkage can be calculated as: ψ dr e = L m i ds e 1 + T r (s) (6) Machines 2022, 10, 1089

of 21
Thus, based on Equation (1), the angular position is given by: From Figure 3, the use of a conventional voltage model's (CVM) estimator could give a large back-EMF error, since the frequency and magnitude could be different in phase. By forming the suggested hybrid model given in the next section, the phase shift could be reduced by compensating for the back-EMF error in the CVM so that a nearly zero back-EMF error can be achieved.

Proposed Hybrid Flux and Speed Estimator with Speed Controller
The proposed model of the hybrid sensorless estimator is illustrated in Figure 4, where it works together with the conventional PI speed controller based on the virtual speed generation shown in Figure 3 in order to gain the minimum rotor speed error. The combination between the modified voltage model and the rotor flux-oriented current model was to give the least back-EMF error in estimating the rotor position of the IM. An adaptive PI control mechanism was added into the flux estimator of the back EMF for error generation. In this mode, the voltage model acted as the reference model, while the rotor flux-oriented current model acted as the adjustable model, to form a new model reference in this adaptive system topology. This system formed a hybrid and interacted with the speed estimator to estimate the rotor speed feedback being fed to the speed controller.  Figure 3, the use of a conventional voltage model's (CVM) estimator co a large back-EMF error, since the frequency and magnitude could be different in By forming the suggested hybrid model given in the next section, the phase shift c reduced by compensating for the back-EMF error in the CVM so that a nearly ze EMF error can be achieved.

Proposed Hybrid Flux and Speed Estimator with Speed Controller
The proposed model of the hybrid sensorless estimator is illustrated in F where it works together with the conventional PI speed controller based on the speed generation shown in Figure 3 in order to gain the minimum rotor speed er combination between the modified voltage model and the rotor flux-oriented model was to give the least back-EMF error in estimating the rotor position of the adaptive PI control mechanism was added into the flux estimator of the back E error generation. In this mode, the voltage model acted as the reference model, w rotor flux-oriented current model acted as the adjustable model, to form a new reference in this adaptive system topology. This system formed a hybrid and in with the speed estimator to estimate the rotor speed feedback being fed to the sp troller.

Sensorless Flux Estimator
The integration of the back EMF into the flux estimator was to estimate in the rotor flux with respect to the synchronously rotating frame for rotor position tion. Two blocks, which were the voltage model and the rotor flux-oriented curren were used to synthesize the rotor flux with respect to the stationary frame for th

Sensorless Flux Estimator
The integration of the back EMF into the flux estimator was to estimate in advance the rotor flux with respect to the synchronously rotating frame for rotor position production. Two blocks, which were the voltage model and the rotor flux-oriented current model, were used to synthesize the rotor flux with respect to the stationary frame for the back-EMF error compensation. The plant was depicted as the flux estimator and can be expressed as: where M = σL s , N = L r L m and k = rotor flux at time k th . The stator flux is obtained by the stator reference current and voltage, such that: Machines 2022, 10, 1089 6 of 21 The back-EMF error in Equation (9) is given as: where G is a positive constant of proportional gain k p and integral gain k i . Thus, the adaptation mechanism law becomes: The extraction of rotor flux from the stator flux is by using: As in the rotor flux-oriented current model, the stator flux is obtained by the stator current reference, and (6) can be expressed as: where T and T −1 are Park and inverse Park transformations, respectively. The stator flux in (13) acts as the feedback to the stator flux controller, which gives the compensated back-EMF error at a reduced level.

Sensorless Speed Estimator
The speed estimator in Figure 4 is another key building block of the proposed scheme. Angular frequency was used to determine the estimated speed. By subtracting slip frequency from angular frequency, as in Equation (4), the rotor speed can be determined. In the case of a pole machine, the slip speed as seen from the rotor frame ω e − ω r_est is equal to ω e − p 2 ω r_est and the rotor speed can be expressed as: The slip speed in (14) is taken into account to reduce the rotor losses in the air gap space. Through this consideration, the IM can give a high output of mechanical power. Thus, the efficiency of the IM is increased directly proportional to the input power.

Controller's Stability on Proposed Hybrid Estimator for IFOC
In a conventional control method, a closed-loop system is where the output speed of the IM is fed back to the speed controller in order to obtain the transfer function. However, for this stability test, the closed-loop system used a current control, which was included for the transfer function derivation. Both stability analyses were tested when the new plant of the controllers was embedded with the hybrid estimator system for checking the new locations of the pole and zero.

Inner Loop Current Control
It was important to ensure that the rotor flux remained constant so that a linear structure could be achieved. The inner loop of the plant system needed to be considered to fulfill the closed-loop transfer function. Therefore, the block diagram in Figure 5 was obtained after several calculations based on the IM circuit analysis with respect to the d-axis of the current controller.
It was important to ensure that the rotor flux remained constant so that a line ture could be achieved. The inner loop of the plant system needed to be consi fulfill the closed-loop transfer function. Therefore, the block diagram in Figure 5 tained after several calculations based on the IM circuit analysis with respect to axis of the current controller. A better understanding can be obtained based on Figure 6, where the clo transfer function was interpreted into the pole-zero map plot. The closed-loop function consisted of two poles, labeled as "X", and one zero, labeled as "O", t plotted on the left-hand side of the pole-zero map. This validated that the clo transfer function of the d − axis current controller was in a stable system. The same implementation was also applied on the axis, where the slip was value that ensured that the q − axis rotor flux was zero. The closed-loop transfer for the axis current controller is shown in Figure 7.  A better understanding can be obtained based on Figure 6, where the closed-loop transfer function was interpreted into the pole-zero map plot. The closed-loop transfer function consisted of two poles, labeled as "X", and one zero, labeled as "O", that were plotted on the left-hand side of the pole-zero map. This validated that the closed-loop transfer function of the d-axis current controller was in a stable system. ture could be achieved. The inner loop of the plant system needed to be consid fulfill the closed-loop transfer function. Therefore, the block diagram in Figure 5 tained after several calculations based on the IM circuit analysis with respect to axis of the current controller. A better understanding can be obtained based on Figure 6, where the clos transfer function was interpreted into the pole-zero map plot. The closed-loop function consisted of two poles, labeled as "X", and one zero, labeled as "O", th plotted on the left-hand side of the pole-zero map. This validated that the clos transfer function of the d − axis current controller was in a stable system. The same implementation was also applied on the axis, where the slip was h value that ensured that the q − axis rotor flux was zero. The closed-loop transfer for the axis current controller is shown in Figure 7.  The same implementation was also applied on the axis, where the slip was held at a value that ensured that the q-axis rotor flux was zero. The closed-loop transfer function for the axis current controller is shown in Figure 7.
fulfill the closed-loop transfer function. Therefore, the block diagram in Figure 5 tained after several calculations based on the IM circuit analysis with respect to axis of the current controller. A better understanding can be obtained based on Figure 6, where the clo transfer function was interpreted into the pole-zero map plot. The closed-loop function consisted of two poles, labeled as "X", and one zero, labeled as "O", t plotted on the left-hand side of the pole-zero map. This validated that the clo transfer function of the d − axis current controller was in a stable system. The same implementation was also applied on the axis, where the slip was value that ensured that the q − axis rotor flux was zero. The closed-loop transfer for the axis current controller is shown in Figure 7.  The obtained closed-loop transfer function stability is shown in Figure 8. As observed, there were two poles and one zero on the left-hand side of the pole-zero map. This confirmed that the controller fulfilled the stability criterion, which meant that it was in a stable system.
The obtained closed-loop transfer function stability is shown in Figure 8 served, there were two poles and one zero on the left-hand side of the pole-zero m confirmed that the controller fulfilled the stability criterion, which meant that it stable system.

Outer Loop Speed Control
A current controller with a closed-loop time constant is ten times smaller th a rotor or mechanical time constant. Under this condition, the current controlle rapidly that it can be neglected for practical purposes. The speed controller is d based on a separately excited DC motor.
The control structure of the IM is shown in Figure 9 with the load torque, m ertia and friction factors. All these parameters could affect the PI's speed control r Initially, the load torque was assumed to be zero. Therefore, the general closed-loo fer function for the speed controller can be realized with the aid of the second-or tem, as in [24], such that: ( ) where n ω is the natural frequency and ξ is the damping ratio. The overshoot age of the system is given as: In order to have an overshoot between 2% and 5% of the speed response, t of ξ and n ω must be greater than 0.7 and 1.8 rad/s, respectively. Based on Equa ξ can be calculated as:

Outer Loop Speed Control
A current controller with a closed-loop time constant is ten times smaller than with a rotor or mechanical time constant. Under this condition, the current controller acts so rapidly that it can be neglected for practical purposes. The speed controller is designed based on a separately excited DC motor.
The control structure of the IM is shown in Figure 9 with the load torque, motor inertia and friction factors. All these parameters could affect the PI's speed control response. Initially, the load torque was assumed to be zero. Therefore, the general closed-loop transfer function for the speed controller can be realized with the aid of the second-order system, as in [24], such that: where ω n is the natural frequency and ξ is the damping ratio. The overshoot percentage of the system is given as: Machines 2022, 10, x FOR PEER REVIEW 8 of 21 The obtained closed-loop transfer function stability is shown in Figure 8. As observed, there were two poles and one zero on the left-hand side of the pole-zero map. This confirmed that the controller fulfilled the stability criterion, which meant that it was in a stable system.

Outer Loop Speed Control
A current controller with a closed-loop time constant is ten times smaller than with a rotor or mechanical time constant. Under this condition, the current controller acts so rapidly that it can be neglected for practical purposes. The speed controller is designed based on a separately excited DC motor.
The control structure of the IM is shown in Figure 9 with the load torque, motor inertia and friction factors. All these parameters could affect the PI's speed control response. Initially, the load torque was assumed to be zero. Therefore, the general closed-loop transfer function for the speed controller can be realized with the aid of the second-order system, as in [24], such that: ( ) where n ω is the natural frequency and ξ is the damping ratio. The overshoot percentage of the system is given as: In order to have an overshoot between 2% and 5% of the speed response, the value of ξ and n ω must be greater than 0.7 and 1.8 rad/s, respectively. Based on Equation (16), ξ can be calculated as: In order to have an overshoot between 2% and 5% of the speed response, the value of ξ and ω n must be greater than 0.7 and 1.8 rad/s, respectively. Based on Equation (16), ξ can be calculated as: It can be defined that the settling is almost 98% of the final response value and the duration is approximately four times the time constant of the signal. When the time constant of the second-order control system is 1/ξω rv , the steady-state settling time can be given as: Referring to the response of the closed-loop IM model under normal condition testing, the values of settling time and overshoot are 0.278 s and 0.025%, respectively. Substituting into (17), the value of ξ is 0.9352. Then, the value of ω n is obtained based on (18) at 15.3861 rad/s.
Substituting the values of ξ and ω n into (15), the general closed-loop transfer function for the speed controller can be expressed as: Based on the second-order closed-loop transfer function above, the stability of the system can be analyzed by solving the poles of the transfer function. It consists of two poles, labeled as "X", which are −14.3891 + 5.4485i and −14.3891 − 5.4485i, with a negative real part and no positive pole, as shown in Figure 10.
It can be defined that the settling is almost 98% of the final response value duration is approximately four times the time constant of the signal. When the ti stant of the second-order control system is 1/ rv ξω , the steady-state settling tim given as: Based on the second-order closed-loop transfer function above, the stabilit system can be analyzed by solving the poles of the transfer function. It consist poles, labeled as "X", which are −14.3891 + 5.4485i and −14.3891 − 5.4485i, with a real part and no positive pole, as shown in Figure 10.

Rotor Resistance Variation
The block diagram of the estimator, depicted in Figure 4, consists of the sta control obtained from the rotor position by compensating the back-EMF error. Sin position is related to speed, which is affected by rotor time constant, it depends resistance, which varies according to temperature during the operation. To test t performance of the stator flux control in the suggested hybrid estimator, rotor re was varied by 50% of the nominal value with an increment of 25%. Figures 11a the effect of rotor resistance variation on the rotor speed of two different estimator The block diagram of the estimator, depicted in Figure 4, consists of the stator flux control obtained from the rotor position by compensating the back-EMF error. Since rotor position is related to speed, which is affected by rotor time constant, it depends on rotor resistance, which varies according to temperature during the operation. To test the drive performance of the stator flux control in the suggested hybrid estimator, rotor resistance was varied by 50% of the nominal value with an increment of 25%. Figure 11a,b show the effect of rotor resistance variation on the rotor speed of two different estimators in two different modes. Both HYB and CVM estimators showed consistent performance in forward and reverse operations, with the results obtained confirming the estimator's functions. The speed characteristics of the two estimators are summarized in Table 1. The two estimators without parameter uncertainty produced similar performances. However, compared with these two estimators, the conventional voltage model, which is not discussed here, took a longer time to track the reference rotor speed (blue line). When the hybrid estimator was applied, there was a shorter tracking time, but it was not remarkable. different modes. Both HYB and CVM estimators showed consistent performance in forward and reverse operations, with the results obtained confirming the estimator's functions. The speed characteristics of the two estimators are summarized in Table 1. The two estimators without parameter uncertainty produced similar performances. However, compared with these two estimators, the conventional voltage model, which is not discussed here, took a longer time to track the reference rotor speed (blue line). When the hybrid estimator was applied, there was a shorter tracking time, but it was not remarkable.
(a) (b) Figure 11. Rotor speed response against rotor resistance variation: (a) in forward operation and (b) in reverse operation.   Based on the rotor speed obtained for the proposed hybrid sensorless estimator, the stability response of the rotor flux against rotor resistance variation is shown in forward operation in Figure 12. As can be observed, the rotor flux response consisted of a few overshoots at the beginning due to the variation of 25% R r and 50% R r . However, the value was in an acceptable range of <20% overshoot with a 0.46 damping ratio. This confirms that the proposed hybrid estimation possesses a good dynamic response against parameter variation during its operation.
Based on the rotor speed obtained for the proposed hybrid sensorless estimator, th stability response of the rotor flux against rotor resistance variation is shown in forward operation in Figure 12. As can be observed, the rotor flux response consisted of a few overshoots at the beginning due to the variation of 25% Rr and 50% Rr. However, the valu was in an acceptable range of <20% overshoot with a 0.46 damping ratio. This confirm that the proposed hybrid estimation possesses a good dynamic response against parame ter variation during its operation. For more details, Table 2 summarizes the rotor flux response. It can be seen that th response is more superior without rotor resistance variation. Hence, with the proposed stator flux control as an adaptive mechanism, the response can be stable against rotor re sistance variation with optimum overshoot and high capability of tracking rotor speed. As recorded (Table 1), it can be concluded from the values of speed error that th hybrid estimator improved motor speed performance by compensating the back-EMF er ror into the CVM via controlling the stator flux of the IM. A minimum error was achieved with respect to reference speed with no oscillations.

Stator Resistance Variation
The voltage model consists of back-EMF error compensation, as in (9), and it is indi cated that stator resistance variation should be taken into account in the study case. A mentioned above, the resistance varies due to temperature during motor operation. I could result in back-EMF error changes, consequently changing the speed's estimated value. So, stator resistance was varied as with rotor resistance in Section 6.1.1 previously As shown in Figures 13a-b, the speed estimation response operated at minimum error fo both forward and reverse operations. The proposed estimator also gave excellent dynami performance under stator resistance variation with a shorter time to track the referenc rotor speed, as tabulated in Table 3. Clearly, the speed estimation efficiently tracked th For more details, Table 2 summarizes the rotor flux response. It can be seen that the response is more superior without rotor resistance variation. Hence, with the proposed stator flux control as an adaptive mechanism, the response can be stable against rotor resistance variation with optimum overshoot and high capability of tracking rotor speed. As recorded (Table 1), it can be concluded from the values of speed error that the hybrid estimator improved motor speed performance by compensating the back-EMF error into the CVM via controlling the stator flux of the IM. A minimum error was achieved with respect to reference speed with no oscillations.

Stator Resistance Variation
The voltage model consists of back-EMF error compensation, as in (9), and it is indicated that stator resistance variation should be taken into account in the study case. As mentioned above, the resistance varies due to temperature during motor operation. It could result in back-EMF error changes, consequently changing the speed's estimated value. So, stator resistance was varied as with rotor resistance in Section 6.1.1 previously. As shown in Figure 13a,b, the speed estimation response operated at minimum error for both forward and reverse operations. The proposed estimator also gave excellent dynamic performance under stator resistance variation with a shorter time to track the reference rotor speed, as tabulated in Table 3. Clearly, the speed estimation efficiently tracked the reference with a speed error range of below 10 rpm and an acceptable range of <20% overshoot. Stator resistance variation had the least effect on IM response compared with rotor resistance variation due to the exact amount of back-EMF error compensation by the adaptation mechanism of stator flux control at the stator side. reference with a speed error range of below 10 rpm and an acceptable range of <20% overshoot. Stator resistance variation had the least effect on IM response compared with rotor resistance variation due to the exact amount of back-EMF error compensation by the adaptation mechanism of stator flux control at the stator side.  When involving temperature, the mutual inductance parameter was also taken into account in this study case. The values of stator and rotor inductance parameters in this research were the same. The variation in inductance also affect the mutual inductance of  Table 3. Rotor speed performance under stator resistance variation.

Mode
Step Info When involving temperature, the mutual inductance parameter was also taken into account in this study case. The values of stator and rotor inductance parameters in this research were the same. The variation in inductance also affect the mutual inductance of the IM. Several studies on inductance against temperature have been conducted, such as in [25], which concluded that inductance remains functional in a high temperature and there is no temperature dependence. However, [26] stated that a 1% change in resistance produced by 2.54 • C corresponds to a 0.0076% increase in inductance. In this research, mutual inductance was varied by 25% and 50% of the nominal value. This means that the temperature was very high during motor operation. As shown in Figure 14a,b, the effect of mutual inductance gave almost the same speed response in both forward and reverse operations either for the CVM or the HYB. Referring to Table 4, the HYB estimator could improve the dynamic performance of the IM compared with the CVM, as the average speed error was below 0.16%. This proves that the proposed estimator works effectively against parameter variation. the IM. Several studies on inductance against temperature have been conducted, such as in [25], which concluded that inductance remains functional in a high temperature and there is no temperature dependence. However, [26] stated that a 1% change in resistance produced by 2.54°C corresponds to a 0.0076% increase in inductance. In this research, mutual inductance was varied by 25% and 50% of the nominal value. This means that the temperature was very high during motor operation. As shown in Figures 14a-b, the effect of mutual inductance gave almost the same speed response in both forward and reverse operations either for the CVM or the HYB. Referring to Table 4, the HYB estimator could improve the dynamic performance of the IM compared with the CVM, as the average speed error was below 0.16%. This proves that the proposed estimator works effectively against parameter variation.

IM Performance under Variable Speed Trajectory
The derivative of rotor position with time was then applied to determine the rotor speed during driving, since the feedback of the speed estimator to the reference speed determines the speed error of the control system. Figure 15a,b show the speed tracking response with respect to the speed reference in order to seek the speed tracking capability. The reference speed was given in the trapezoidal signal input with a constant value during the motor's starting operation in order to see the transient behavior. The variable speed trajectory was applied in forward and reverse motoring operations with the proposed hybrid estimator of the stator flux control.

IM Performance Under Variable Speed Trajectory
The derivative of rotor position with time was then applied to determine the rotor speed during driving, since the feedback of the speed estimator to the reference speed determines the speed error of the control system. Figures 15a-b show the speed tracking response with respect to the speed reference in order to seek the speed tracking capability. The reference speed was given in the trapezoidal signal input with a constant value during the motor's starting operation in order to see the transient behavior. The variable speed trajectory was applied in forward and reverse motoring operations with the proposed hybrid estimator of the stator flux control.  Using accurate rotor position tracking, the speed error was minimized nearly to zero. A good tracking capability was obtained in both reverse and forward operations with no overshoots at the starting of the IM. As can be seen between 1.5 and 3 s, the speed errors were recorded in Table 5 from the reference values at 0.082% and 0.132%, respectively. Meanwhile, in the existing voltage model estimator, it contributed about 0.243% and 0.285% errors for both speed trajectories. With the 0.161% and 0.153% margin errors, the hybrid estimator has shown an excellent performance in estimating speeds using the stator flux control under a variable speed trajectory. This shows that the back-EMF error was completely compensated into the voltage model in order to create the correct position of the rotor.

Mode
Step Info The stability of the speed controller's response can be seen in Figure 16a-c, which present the forward motoring operation. In Figure 16a, the spiral graph shows a highly stable response during full-speed operation without any ripples. Meanwhile, as shown in Figure 16b, at a minimum low speed of 500 rpm, the current faced a few ripples with a stable response. This happened due to the integration of very small signals of angular frequency that made the rotor position shifted a few degrees. When reaching 1000 rpm, as shown in Figure 16c, the spiral response of the current had fewer ripples compared with the response at a lower speed. However, the speed estimation response was within the boundary of the speed reference's range and a high speed-tracking capability was obtained with accurate rotor position estimation.
Meanwhile, in the existing voltage model estimator, it contributed about 0.243% and 0.285% errors for both speed trajectories. With the 0.161% and 0.153% margin errors, the hybrid estimator has shown an excellent performance in estimating speeds using the stator flux control under a variable speed trajectory. This shows that the back-EMF error was completely compensated into the voltage model in order to create the correct position of the rotor.

Mode
Step Info The stability of the speed controller's response can be seen in Figures 16a-c, which present the forward motoring operation. In Figure 16a, the spiral graph shows a highly stable response during full-speed operation without any ripples. Meanwhile, as shown in Figure 16b, at a minimum low speed of 500 rpm, the current faced a few ripples with a stable response. This happened due to the integration of very small signals of angular frequency that made the rotor position shifted a few degrees. When reaching 1000 rpm, as shown in Figure 16c, the spiral response of the current had fewer ripples compared with the response at a lower speed. However, the speed estimation response was within the boundary of the speed reference's range and a high speed-tracking capability was obtained with accurate rotor position estimation.

IM Performance at Very Low Speed Under Regenerative Braking and Zero-Speed Conditions
In order to behave like in a real EV operation, the four-quadrant operation needed to be fulfilled. In the above section, forward and reverse motoring had been discussed, and is not repeated here. Forward braking and reverse braking are under the regenerative

IM Performance at Very Low Speed under Regenerative Braking and Zero-Speed Conditions
In order to behave like in a real EV operation, the four-quadrant operation needed to be fulfilled. In the above section, forward and reverse motoring had been discussed, and is not repeated here. Forward braking and reverse braking are under the regenerative mode and could be realized by slowing down the reference speed or applying a mechanical load to the IM to reduce the speed to 0 rpm. In this research, a mechanical load was applied at τ = 20 Nm to show the IM braking and operating at a very low speed. In addition, a mechanical load was also applied at τ = 40 Nm to show the IM braking until the steadystate zero-speed condition. Figure 17a,b show the response of braking in two different operations. As can be seen, the IM operated at a very low speed by following the reference speed at 100 rpm. At time τ = 1.5 s, a mechanical load was applied to show the braking system slowing down the IM until the extra-low speed of 50 rpm and zero-speed condition at 0 rpm. As shown in Table 6, the speed error obtained was less than 1 rpm for the proposed HYB estimator compared with the speed error for the CVM, which was more than 1 rpm at zero speed. This shows that the speed tracking capability successfully reduced speed estimation errors. This confirms that the proposed HYB is able to improve speed tracking capability and operate the IM under very-low-speed conditions. reference speed at 100 rpm. At time t = 1.5 s, a mechanical load was applied to show the braking system slowing down the IM until the extra-low speed of 50 rpm and zero-speed condition at 0 rpm. As shown in Table 6, the speed error obtained was less than 1 rpm for the proposed HYB estimator compared with the speed error for the CVM, which was more than 1 rpm at zero speed. This shows that the speed tracking capability successfully reduced speed estimation errors. This confirms that the proposed HYB is able to improve speed tracking capability and operate the IM under very-low-speed conditions.    This also indicates that the HYB estimator was able to operate under zero-speed conditions, where its speed estimation had the least mismatch compared with the CVM in tracking the actual zero speed. Figure 18 shows the speed estimation by the HYB, as well as the speed estimation error of under 2 rpm, approaching zero with accurate flux estimation. Furthermore, the results confirmed that the rotor flux position was at accurate reference speed, while not losing stability in very-low-speed operation. This satisfies the supremacy of the HYB's flux estimation in extracting accurate speed estimation under very-low-speed and zero-speed conditions. ditions, where its speed estimation had the least mismatch compared with the CVM in tracking the actual zero speed. Figure 18 shows the speed estimation by the HYB, as well as the speed estimation error of under 2 rpm, approaching zero with accurate flux estimation. Furthermore, the results confirmed that the rotor flux position was at accurate reference speed, while not losing stability in very-low-speed operation. This satisfies the supremacy of the HYB's flux estimation in extracting accurate speed estimation under verylow-speed and zero-speed conditions.

IM Performance Under Load Torque Disturbances
Since speed and torque are related to each other, a justification can be made under the load torque disturbances, as presented in Figures 19a and 19b for the forward and reverse operations, respectively, that torque does not give any effect to the controller. Torque disturbances were set up at times t = 0.7 s, 1.3 s and 2 s with torque values of 50 Nm, 75 Nm and 100 Nm, respectively. As depicted in Figure 19, there was a quick convergence in the speed response after the load disturbances were rejected. The quality of estimation also can be seen in the speed errors, which were relatively low as expected. The speed estimation response tracked the constant reference speed in spite of load torque disturbances both in forward and reverse operations. In addition, the dynamics of the IM was still maintained with no overshoots and oscillations during the starting of the IM.

IM Performance under Load Torque Disturbances
Since speed and torque are related to each other, a justification can be made under the load torque disturbances, as presented in Figure 19a and 19b for the forward and reverse operations, respectively, that torque does not give any effect to the controller. Torque disturbances were set up at times t = 0.7 s, 1.3 s and 2 s with torque values of 50 Nm, 75 Nm and 100 Nm, respectively. As depicted in Figure 19, there was a quick convergence in the speed response after the load disturbances were rejected. The quality of estimation also can be seen in the speed errors, which were relatively low as expected. The speed estimation response tracked the constant reference speed in spite of load torque disturbances both in forward and reverse operations. In addition, the dynamics of the IM was still maintained with no overshoots and oscillations during the starting of the IM.  Table 7 shows the summary of the torque load effect on the two estimators. Clearly, there were smaller speed drops in the hybrid (HYB) estimator compared with the CVM due to the compensated EMF error that changed the stator flux response. A faster recovery time was noticed in the proposed HYB estimator with an average of 0.078 s, which gave  Table 7 shows the summary of the torque load effect on the two estimators. Clearly, there were smaller speed drops in the hybrid (HYB) estimator compared with the CVM due to the compensated EMF error that changed the stator flux response. A faster recovery time was noticed in the proposed HYB estimator with an average of 0.078 s, which gave the IM an excellent performance. Additionally, there were no ripples on the speed response, which meant that the IM did not suffer from torque cogging. Table 7. Rotor speed performance under load torque disturbances.

Mode
Step Info

Performance Evaluation Comparison
Based on the presented simulation result, the good dynamic performance of the IM was successfully achieved by the suggested model. The HYB estimator was effectively more robust than the CVM. The results demonstrate that the system exhibits good speed estimation accuracy under different conditions. This is supported by the stability effect given in Figure 16 for speed and current controllers, which are responsible for controlling the torque and speed of the IM. An evaluation table (Table 8) is provided to compare with the findings from previous research works, as examined in [27], to clarify the improvements that have been achieved with the suggested model.  Notes: * calculated from measured and reference value. V, voltage; I, current; Wr, rotor speed; EMF, electromotive force; LS, low speed; HS, high speed; F; forward; R, reverse; PI, proportional integral; SV-PWM, space vector pulse width modulation; MED, medium; HI, high; Est, estimated; Act, actual.

Conclusions
This paper presents a hybrid estimator model, along with the accuracy analysis of speed and torque tracking of an IFOC in estimating rotor flux position and speed to control the speed of the IM. The parameters of the IM are listed in the Appendix A, Table A1. In order to enhance speed control, a new model of hybrid estimator was proposed with stator flux control as the adaptation mechanism. The proposed estimator was investigated under four cases: (1) parameter sensitivity against rotor and stator resistance and mutual inductance variation, (2) variable speed trajectory, (3) very low speed in regenerative mode and zero-speed conditions and (4) variable load torque disturbances. The analyses showed that the three controllers were confirmed in a stable operation, as evidenced by the pole-zero maps.
All results proved that parameter sensitivity did not affect speed response. All speed responses also showed good tracking capability under variable speed trajectories and had faster recovery time under torque disturbances by using the hybrid estimator as compared with the conventional model. The proposed control scheme was successfully validated in various cases. The stability analysis and reported results showed excellent performance of the IM in both forward and reverse operations at a steady state and during transients.