Next Article in Journal
Efficient One-Class False Data Detector Based on Deep SVDD for Smart Grids
Previous Article in Journal
Development of the Kinetic Equation of the Groove Corrosion Process for Predicting the Residual Life of Oil-Field Pipelines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differentiator-Based Output Feedback MPPT Controller for DFIG Wind Energy Conversion Systems with Minimal System Information

Department of Electrical and Control Engineering, Mokpo National University, Chonnam 58554, Republic of Korea
Energies 2023, 16(20), 7068; https://doi.org/10.3390/en16207068
Submission received: 2 September 2023 / Revised: 3 October 2023 / Accepted: 5 October 2023 / Published: 12 October 2023
(This article belongs to the Topic Advances in Wind Energy Technology)

Abstract

:
This paper introduces a novel differentiator-based maximum power point tracking (MPPT) controller for a wind energy conversion system (WECS) equipped with a doubly fed induction generator (DFIG). Building upon our previous algorithms, the proposed controller reduces the need for detailed system information and displays enhanced robustness against parameter variations and disturbances. The innovation lies in the elimination of the need for explicit functional forms or specific parameter values in the system’s dynamics, relying solely on relative degrees and control directions. Utilizing a higher-order switching differentiator (HOSD), this paper outlines a method for overestimating the time derivatives of system outputs, thereby simplifying both the controller design and stability analysis. Compared to existing solutions, the proposed method requires minimal information, offers simpler control law structures, and follows a systematic design approach with fewer design constants. Simulation results demonstrate the efficacy of the proposed controller in both tracking maximum power and regulating reactive power to zero, suggesting a more efficient and simplified approach to MPPT control in WECS.

1. Introduction

Due to concerns surrounding pollution, high costs, potential depletion, and other adverse impacts associated with traditional energy sources, renewable energy generation has garnered substantial interest. Among the various renewable energy technologies, the wind energy conversion system (WECS) stands out as one of the most developed and widely used, primarily because it is clean, inexhaustible, and broadly accessible. Furthermore, the doubly fed induction generator (DFIG) is deemed as a stable and efficient fixed-speed wind turbine system [1,2]. With its lower converter cost and reduced power losses, the DFIG has attracted considerable attention from numerous researchers.
While the conventional Proportional–Integral (PI) controller exhibits satisfactory performance across numerous WECS applications, efforts have been made to enhance WECS performance through the introduction of optimized PI controllers [3,4]. However, these controllers exhibit limitations, notably their inability to adapt to variations in machine parameters and uncertainties. Therefore, any deviation in operational conditions from those under which the PI controller parameters were optimized cannot assure optimal operation. This necessitates real-time PI controller parameter tuning in response to wind speed variations and potential parametric alterations, often due to technical complications such as mechanical wear or machine overheating. In a practical sense, implementing these techniques remains challenging and costly. The difficulties largely stem from determining the optimal gains of the controller needed to achieve control objectives across all operating regions and to adapt well to changes in system parameters. These shortcomings inherent in PI or optimized PI controllers have prompted the development and implementation of robust and adaptive control techniques for WECS. Proposals have been put forth for more robust control strategies to supplant PI controllers in order to enhance precision and accuracy performance [5].
Inherently, WECS exhibit strongly nonlinear dynamic equations. Furthermore, their system parameters are prone to variations due to factors such as ambient temperature changes and mechanical wear. Conventionally, the design of stabilizing controllers for nonlinear systems characterized by unstructured uncertainties including parameter variations and external disturbances have largely leveraged sliding mode control (SMC) methods [6,7] and adaptive control algorithms with universal approximators such as neural networks (NN) or fuzzy logic systems (FLS) [8,9,10,11,12,13,14,15,16]. Although these NNs and FLS approximators have found wide application in addressing system uncertainties, they require a complex structure to ensure approximation capabilities. Furthermore, they necessitate the online updating of a multitude of adaptive parameters, thereby leading to a heightened computational load and dynamic order of the controller. Recently, control strategies aiming to simplify complex control formulas without compromising their performance have been proposed. This includes prescribed performance control (PPC) techniques that guarantee a predefined tracking performance irrespective of system uncertainties, and without the need for approximation [17,18,19,20]. The PPC framework significantly simplifies the controller structure by eliminating the need for universal approximators. Yet, the steps of the backstepping design continue to be an integral part of PPC methodologies, making them vulnerable to faults or large disturbances after the transient period. More recently, the proposal of differentiator-based controllers emerged, which address system uncertainties by overestimating the time derivatives of the output tracking error, eliminating the need for universal approximations [21,22]. Despite the surge in research in this field, a large proportion of the studies target single-input single-output (SISO) nonlinear systems, leaving multi-input multi-output (MIMO) systems relatively under-explored.
Recently, SMC methods have been applied to maximum power point tracking (MPPT) control of WECS [23,24]. However, these control strategies are constructed under the assumption that both the formulas and parameters of system dynamics are fully known, in order to cancel out the system’s nonlinearities, making this approach somewhat restrictive. To alleviate these restrictions, adaptive fuzzy or neuro-controllers are introduced in [25,26,27,28,29,30]. While these control techniques, which are based on universal approximators, offer the benefit of not necessitating prior knowledge about the nonlinear functions in the system’s dynamic equation, they come with the complexity of intricate control laws and formulas for updating adjustable parameters. This complexity poses significant challenges for the actual implementation of the control algorithm. In [31,32], the PPC schemes are utilized for MPPT control of WECS. Nonetheless, the PPC algorithm is susceptible to faults or abrupt disturbances that manifest in a steady state. In this paper, building upon the control algorithm presented in [21,22], we propose a novel differentiator-based MPPT controller for WECS equipped with DFIG. The advantage of the proposed controller is twofold: it obviates the need for information about the nonlinear functions present in the system’s dynamic equation, and it can adeptly handle both parameter variations and disturbances. When designing the controller, we assume that the time derivatives of outputs are aggregated unknown functions of the inputs and time. Consequently, apart from information regarding relative degrees and control direction, specific functional expressions and parameter values of the system dynamics are not considered in the controller’s formulation. Employing HOSD [33], we over-estimate the time derivatives of the system outputs, facilitating both a more streamlined controller structure and a more straightforward stability analysis. The advantages of the controller presented in this paper, relative to existing research, are as follows:
  • The proposed controller demands significantly less information about the system’s dynamic equation. The control formulation relies solely on the relative degrees between inputs and outputs, the directions of control inputs, and the measured output values.
  • Owing to the absence of universal approximators, the structures of the control laws are comparatively straightforward, while still ensuring the asymptotic stability of the outputs.
  • The proposed output feedback control algorithm offers a cohesive and systematic approach for crafting control laws.
  • The quantity of design constants is minimal in comparison to that of other methods.
Simulation results are provided to showcase the efficacy of the proposed controller and the consistency inherent in its design.
The organization of this paper is as follows. Section 2 provides the dynamic equations for DFIG-WECS and derives the state equations. Section 3 introduces a differentiator-based output feedback MPPT controller. Section 4 presents simulation results, including a comparison of the performance between the proposed controller and a conventional PI controller. Finally, Section 5 offers conclusions. The symbols used in describing the DFIG-WECS and their meanings in this paper are summarized in Table 1.

2. Dynamic Model of the WECS with DFIG

2.1. Model of WECS

The mechanical power generated by the wind turbine is described by equation
P t = 1 2 ρ A C p ( λ , β ) v w 3 ,
where ρ is the air density, A is the area swept by the blades, C p ( λ , β ) is the power coefficient, λ is the tip speed ratio, β is the blade pitch angle, and v w is the wind speed. The tip speed ratio is defined as
λ = v t v w = ω R v w ,
where v t is the tip speed of the blade, ω is the rotational velocity of the wind turbine’s rotor, and R is the radius of the blade. The conceptual diagram illustrating these parameters is depicted in Figure 1.
The power coefficient serves as a metric for quantifying the efficiency of the wind turbine in converting the kinetic energy of the wind into mechanical energy. It is defined as the ratio of the power harnessed by the wind turbine to the total power available in the wind. In this study, we employ the following equation for the power coefficient, as adapted from the literature [34]:
C p ( λ , β ) = c 1 c 2 λ β c 3 β c 4 e c 5 λ β + c 6 λ
and
1 λ β = 1 λ + 0.008 β 0.035 β 3 + 1 .
Here, coefficients c i ( i = 1 , , 6 ) are obtained from [34] and are detailed in the simulation section. The aerodynamic torque of turbine T a is given by equation
T a = P t ω = 1 2 ω ρ π R 2 C p ( λ ) v w 3 .
The first-order dynamic model of the wind turbine can be expressed as
J ω ˙ = T a D ω n g T e m ,
where J represents the total inertia, which is the sum of the inertia from the turbine and the generator:
J = J t + n g 2 J g .
J t is the inertia of the turbine, J g is the inertia of the generator, n g is the transmission ratio of the gearbox, D is the damping constant of the wind turbine, and T e m is the electromagnetic torque, which is defined subsequently. Overall schematic of the proposed control system is illustrated in Figure 2 and typical power coefficient curves are depicted in Figure 3 for various blade pitch angles.
At low wind speeds, λ should be maintained at its optimal value λ to maximize conversion efficiency. Therefore, the desired rotor speed ω d can be calculated using the current wind speed as follows:
ω d = λ R v w .
The primary control objective for the WECS under consideration is to ensure that the angular speed of the blade ω tracks the desired value given by Equation (8). This is commonly referred to as MPPT in the subsequent discussion.

2.2. Mathematical Model of DFIG

In the equations that follow, v s d , v s q , v r d , v r q , i s d , i s q , i r d , i r q , ϕ s d , ϕ s q , ϕ r d , and ϕ r q represent the d and q components of the stator and rotor voltages, currents, and fluxes, respectively. The commonly used electrical equations for the DFIG in the Park reference frame are given as follows:
v s d = R s i s d + d ϕ s d d t ω s ϕ s q ,
v s q = R s i s q + d ϕ s q d t + ω s ϕ s d ,
v r d = R r i r d + d ϕ r d d t ω r ϕ r q ,
v r q = R r i r q + d ϕ r q d t + ω r ϕ r d .
Here, R s and R r denote the resistances of the stator and rotor, respectively. Symbols ω s and ω r represent the stator and rotor electrical angular speeds in the synchronous reference frame. Since the rotor is directly connected to the blade through a gearbox, equation ω r = n g ω holds true. The dynamics of the stator and rotor fluxes are described by the following equations:
ϕ s d = L s i s d + L m i r d ,
ϕ s q = L s i s q + L m i r q ,
ϕ r d = L r i r d + L m i s d ,
ϕ r q = L r i r q + L m i s q ,
where L s and L r denote the inductances of the stator and rotor, respectively, and L m represents the mutual inductance. The electromagnetic torque is described by
T e m = p r ( ϕ s d i s q ϕ s q i s d ) ,
where p r is the number of pole pairs. Furthermore, the stator’s active and reactive powers are given by
P s = v s d i s d + v s q i s q ,
Q s = v s d i s q v s q i s d .
The state equations are derived based on the assumption that both stator and rotor variables are referred to the stator reference Park frame [23,34]. Given this orientation, the following relationships hold:
ϕ s d = ϕ s ,
ϕ s q = 0 ,
and
v s d = 0 ,
v s q = v s ,
where ϕ s represents the total flux and v s signifies the total voltage of the stator. From these, i s d and i s q can be derived using Equations (13) and (14):
i s d = 1 L s ( ϕ s L m i r d ) ,
i s q = L m L s i r q .
Utilizing the following equations, which arise from Equations (16) and (25),
d ϕ r q d t = L r d i r q d t + L m d i s q d t ,
d i s q d t = L m L s d i r q d t ,
the dynamic equation for i r q can be derived from Equation (12) as follows:
v r q = R r i r q + d ϕ r q d t + ω r ϕ r d = R r i r q + L r d i r q d t + L m d i s q d t + ω r L r i r d + L m i s d = R r i r q + L r L m 2 L s d i r q d t + ω r L r L m 2 L s i r d + L m L s ϕ s = R r i r q + σ L r d i r q d t + ω r σ L r i r d + L m L s ϕ s ,
where σ = 1 L m 2 L r L s , leading to
d i r q d t = R r σ L r i r q ω r i r d + L m σ L r L s ϕ s + 1 σ L r v r q .
Dynamic equation for ϕ s is induced from (9) using (24) as
d ϕ s d t = R s i s d = R s L s ϕ s + R s L m L s i r d .
Utilizing the following equations, derived from (15) and (24),
d ϕ r d d t = L r d i r d d t + L m d i s d d t ,
d i s d d t = 1 L s d ϕ s d t L m L s d i r d d t ,
the dynamic equation for i r d can be derived from (11) as follows:
v r d = R r i r d + d ϕ r d d t ω r ϕ r q = R r i r d + L r d i r d d t + L m d i s d d t ω r L r i r q + L m i s q = R r i r d + L r d i r d d t + L m 1 L s d ϕ s d t L m L s d i r d d t ω r L r i r q L m 2 L s i r q = R r i r d + L r L m 2 L s d i r d d t + L m L s d ϕ s d t + ω r L r L m 2 L s i r q = R r i r d + σ L r d i r d d t + L m L s R s L s ϕ s + R s L m L s i r d + ω r σ L r i r q = R r + R s L m 2 L s 2 i r d + σ L r d i r d d t L m R s L s 2 ϕ s + ω r σ L r i r q .
This leads to
d i r d d t = 1 σ L r R r + R s L m 2 L s 2 i r d + ω r i r q + L m R s σ L r L s 2 ϕ s + 1 σ L r v r d .
The resulting four dynamic equations that constitute the state-space representation of the WECS with DFIG are collected as follows:
d ω d t = 1 J ( T a D ω n g T e m ) ,
d i r q d t = R r σ L r i r q ω r i r d + L m σ L r L s ϕ s + 1 σ L r v r q ,
d i r d d t = 1 σ L r R r + R s L m 2 L s 2 i r d + ω r i r q + L m R s σ L r L s 2 ϕ s + 1 σ L r v r d ,
d ϕ s d t = R s L s ϕ s + R s L m L s i r d .
Here, T e m = p r ϕ s L m L s i r q , as indicated by (17). The active and reactive powers, represented by (18) and (19), are reformulated as
P s = v s L m L s i r q ,
Q s = v s L s ( ϕ s L m i r d ) .
In this case, v s = ω s ϕ s R s L m L s i r q . Note that if R s is assumed to be approximately zero, then d ϕ s d t 0 , and the dynamics become identical to those presented in [23]. In this paper, however, R s is non-zero, and ϕ s is a time-varying signal governed by (38).

3. Design of Output Feedback Controllers

The state-space equations include two control inputs, v r q and v r d . The primary control objective is to drive ω to ω d for optimal power extraction from the wind turbine. The secondary objective is to regulate reactive power Q s to zero. From Equations (35)–(38), it is evident that the second-order derivative of ω is a function of v r q , and the first-order time derivative of Q s involves v r d . Based on these observations, the following dynamic equations are formulated:
ω ¨ = f 1 ( v r q , t ) ,
Q ˙ s = f 2 ( v r d , t ) ,
where functions f 1 ( · ) and f 2 ( · ) are considered to be unknown. These functions are explicit functions of time t, and may involve state variables ( ω ( t ) , i r d ( t ) , i r q ( t ) , ϕ s ( t ) ) as well as time-varying parameters, external disturbances, and measurement errors that are difficult to model precisely. It is also inferred that the inputs v r q and v r d can be employed to control ω ( t ) and Q s ( t ) independently. That is, these inputs and outputs are decoupled, enabling the independent design of two controllers. Subsequent subsections present the design of control laws to achieve these objectives.
In what follows, the two-norm of vector x is denoted by | x | , and the absolute value of the scalar v is also indicated by | v | . Notation a ( t ) 0 is employed to signify lim t a ( t ) = 0 , indicating that a ( t ) converges to zero as t approaches infinity. Similarly, a ( t ) b ( t ) denotes that a ( t ) asymptotically approaches b ( t ) as t approaches infinity, or lim t a ( t ) = b ( t ) .

3.1. MPPT Control

In this paper, the differentiator-based output feedback control scheme, as discussed in the introduction and detailed in [21,22,35], is employed for MPPT control of the WECS as represented by Equations (35)–(38). It is evident that the relative degree between the outputs ω and v r d is two, as v r d first appears in the second-order time derivative of ω . To construct the feeding signal a 1 ( t ) for the higher-order switching differentiator (HOSD), a second-order linear filter is required, as given by
w ˙ 11 = w 11 + w 12 , w ˙ 12 = w 12 + v r q .
Here, w 11 and w 12 are the state variables of the Linear Time-Invariant (LTI) filter (43). Wedefine the tracking error e 1 as
e 1 = ( ω ω d ) g 1 ,
where g 1 > 0 is a gain constant to be determined. Subsequently, signal a 1 ( t ) is given by
a 1 ( t ) = e 1 w 11 .
To formulate the control law, the following HOSD is introduced, as per [21,35]:
Lemma 1. 
Let a 1 ( t ) be a signal whose time derivatives are piecewise-bounded, such that | a ˙ 1 | L 11 and | a ¨ 1 | L 12 , for some positive constants L 11 , L 12 . Consider dynamics
α ˙ 11 = 10 L 1 ϵ 11 + σ 11 σ ˙ 11 = L 1 sgn ( ϵ 11 ) } ,
α ˙ 12 = 7 L 1 ϵ 12 + σ 12 σ ˙ 12 = L 1 sgn ( ϵ 12 ) } ,
where ϵ 11 = a 1 ( t ) α 11 and ϵ 12 = σ 11 α 12 . If L 1 is chosen sufficiently large such that L 1 > max { L 11 , L 12 } , then σ 11 a ˙ 1 and σ 12 a ¨ 1 .
The comprehensive proof is presented in [33]. From Lemma 1, the following equations hold:
σ 11 = a ˙ 1 + δ 11 ( t ) = e ˙ 1 w ˙ 11 + δ 11 ( t ) = e ˙ 1 + w 11 w 12 + δ 11 ( t ) ,
σ 12 = a ¨ 1 + δ 12 ( t ) = e ¨ + w ˙ 11 w ˙ 12 + δ 12 ( t ) = e ¨ w 11 + 2 w 12 v r q + δ 12 ( t ) ,
where δ 11 ( t ) and δ 12 ( t ) are asymptotically vanishing estimation errors; that is, δ 11 ( t ) 0 and δ 12 ( t ) 0 . If the error vector is defined as e 1 = [ e 1 , e ˙ 1 ] T , its estimate can be determined from Equation (48) as follows:
e ^ 1 = e σ 11 w 11 + w 12 .
Subsequently, the output feedback control input v r q is proposed as
v r q = σ 12 ( w 11 2 w 12 ) k 1 T e ^ 1 ,
where k 1 = [ k 11 , k 12 ] T is determined so that the polynomial s 2 + k 12 s + k 11 is Hurwitz.
Theorem 1. 
The control input v r q given by (51), in conjunction with the HOSD as described in Lemma 1, and the LTI filter (43) ensures that ω asymptotically tracks ω d .
Proof. 
From Equations (51) and (49), the following equality can be straightforwardly derived:
v r q = σ 12 ( w 11 2 w 12 ) k 1 T e ^ 1 = ( e ¨ 1 w 11 + 2 w 12 v r q + δ 12 ( t ) ) ( w 11 2 w 12 ) k 1 T e 1 k 12 δ 11 ( t ) = e ¨ 1 + v r q k 1 T e 1 δ 12 ( t ) k 12 δ 11 ( t ) .
Defining d 1 ( t ) = δ 12 ( t ) k 12 δ 11 ( t ) , we can induce
e ¨ 1 = k 1 T e + d 1 ( t ) .
This can be rewritten in vector form as
e ˙ 1 = A e 1 + b d 1 ( t ) ,
where
A = 0 1 k 11 k 12 , b = 0 1 .
Positive definite matrices P and Q exist such that A T P + P A + Q = 0 . The time derivative of the Lyapunov function V 1 = e 1 T P e 1 is given by
V ˙ 1 = e 1 T Q e 1 + 2 e 1 T P b d 1 ( t ) λ m i n ( Q ) | e 1 | 2 + 2 | e 1 | λ m a x ( P ) | d 1 ( t ) | ,
where λ min ( · ) and λ max ( · ) denote the minimum and maximum eigenvalues of a matrix, respectively. From this inequality, it follows that if | e 1 | > μ | d 1 ( t ) | , where μ = 2 λ max ( P ) λ min ( Q ) , then V ˙ 1 < 0 . Given that d 1 ( t ) converges to zero asymptotically, it can be concluded that | e 1 | is also asymptotically stable. □

3.2. Regulating Q s to Zero

The second control objective outlined in [23] involves regulating the reactive power Q s to zero using v r d . The relative degree between the output Q s and the input v r d is one, simplifying the control law compared to the MPPT control law discussed in the preceding subsection. To generate the feeding signal a 2 ( t ) for the second HOSD, we adopt the following first-order linear filter:
w ˙ 21 = w 21 + v r d ,
where w 21 is a state variable of this filter. We let the tracking error be defined as
e 2 = ( Q s 0 ) g 2 .
Here, g 2 > 0 is a tunable amplification parameter. Signal a 2 ( t ) is then given by
a 2 ( t ) = e 2 w 21 .
We introduce the following HOSD that is capable of observing Q ˙ s .
Lemma 2 
([36]). Let a 2 ( t ) be a signal whose time derivatives are bounded in the piecewise sense such that | a ˙ 2 | L 2 for some positive constants L 2 . Consider the following dynamics:
α ˙ 21 = 10 L 2 ϵ 21 + σ 21 σ ˙ 21 = L 2 sgn ( ϵ 21 ) , }
where ϵ 21 = a 2 ( t ) α 21 . If L 2 is chosen sufficiently large such that L 1 > L 2 , then σ 21 a ˙ 2 .
The detailed proof is presented in [33]. From Lemma 2, the following equation holds:
σ 21 = a ˙ 2 + δ 21 ( t ) = e ˙ 2 w ˙ 21 + δ 21 ( t ) = e ˙ 2 + w 21 v r d + δ 21 ( t ) ,
where δ 21 ( t ) is asymptotically vanishing estimaion error, i.e., δ 21 ( t ) 0 . Then, the control input v r d is proposed as
v r d = σ 21 + w 21 k 2 e 2 ,
where k 2 > 0 is a design constant.
Theorem 2. 
The control input v r d as described in (62), in conjunction with HOSD (2) and the LTI filter (57), ensures that the reactive power Q s ( t ) asymptotically converges to zero.
Proof. 
The following equality can be straightforwardly derived from Equations (62) and (61):
v r d = σ 21 + w 21 k 2 e 2 = ( e ˙ 2 + w 21 v r d + δ 21 ( t ) ) + w 21 k 2 e 2 = e ˙ 2 + v r d k 2 e 2 δ 21 ( t ) .
From the last equality, the following relationship is deduced:
e ˙ 2 = k 2 e 2 δ 21 ( t ) .
We let V 2 = e 2 2 2 k 2 serve as the second Lyapunov function. Its time derivative is
V ˙ 2 = e 2 e ˙ 2 k 2 = e 2 2 e 2 δ 21 ( t ) k 2 | e 2 | 2 + | e 2 | | δ 21 ( t ) | k 2 .
This implies that if | e 2 | > | δ 21 ( t ) | k 2 , then V ˙ 2 < 0 . Given that δ 21 ( t ) converges to zero asymptotically, | e 2 | also asymptotically converges to zero. □
Remark 1. 
It is worth noting that the proposed control scheme does not require information about system parameters or nonlinear functions in the WECS dynamics. Only measured values ( v w , ω, and Q s ) and the information about λ to calculate ω d are needed. Other details such as the exact structure of the C p ( · ) function, precise system parameters, and the configuration of the dynamic equations are not necessary for the formulation of control laws. Moreover, the control strategy is robust to unknown or varying system parameters and unmodeled dynamic structures within the WECS.
Remark 2. 
The proposed controller, based on differentiator-based control techniques, shares similar advantages with data-driven control [37] in that it does not require a precise dynamic model of the system under control. However, while data-driven control involves an identification process based on observed system behavior data, differentiator-based control eliminates the need for such an identification phase. It enables the immediate design of output feedback controllers using minimal information, such as relative degree and control direction.

4. Simulations

In this section, the effectiveness of the proposed control strategy is validated through simulations using a 1.5 MW DFIG-WECS model as described in [23]. All simulations were conducted using Python libraries, specifically NumPy, SciPy, and Matplotlib [38]. The parameters relevant to this WECS model are provided in Table 2, while the design parameters of the proposed controllers are detailed in Table 3. In addition to the proposed controller, simulation experiments were also conducted using a PI controller to compare performance. The control equations for the PI controller are as follows:
v r q = 1000 e 1 + 200 0 t e 1 d t ,
v r d = 1000 e 2 + 200 0 t e 2 d t .
Here, the error signals e 1 and e 2 are defined as in Equations (44) and (58), respectively, with the values of g 1 and g 2 as given in Table 3. The gains for the PI controller were selected empirically through multiple simulation runs, prioritizing output tracking performance over transient response.
As demonstrated in Figure 4, the wind speed model incorporates both variability and rapid fluctuations, providing a realistic test scenario for evaluating the proposed control scheme. As the machinery operates over an extended period, factors such as wear and tear can lead to changes in inertia coefficients or damping constants. Additionally, variations in ambient temperature can alter the resistance and inductance of the generator. To showcase the controller’s performance, a simulation was executed where the constants J, D, R s , R r , L s , L r , and L m were changed by an extreme 40% at the 100 s mark, and the results were presented to illustrate its robustness under such conditions.
Figure 5 depicts the ω ( t ) ( = y 1 ) along with its desired value ω d ( t ) , illustrating rapid and accurate tracking performance of the proposed controller. Moreover, the transient response of the proposed controller is significantly better compared to that of the PI controller. The output tracking performance of the proposed controller closely aligns with ω d ( t ) , whereas the PI controller struggles to accurately track ω d ( t ) at inflection points.
Figure 6 presents the regulating performance for reactive power Q s ( = y 2 ) , confirming the controller’s ability to regulate Q s to zero effectively. In these results, the proposed controller stands out for its oscillation-free transient response and more pronounced regulation performance. Additionally, as can be seen in Figure 5, both controllers exhibit insensitivity to parameter variations in ω at t = 100 s, while Q s experiences a spike before immediate return to proper regulation. Figure 7 shows the actual power coefficient C p ( λ , 0 ) , which almost perfectly tracks the optimal command C p ( λ , 0 ) 0.48 by the proposed controller. In contrast, the PI controller exhibits significant oscillations in its transient response and shows inferior steady-state performance. In Figure 8, the trajectories of the active power P s are depicted, and Figure 9 illustrates the control voltage inputs v r d and v r q . In Figure 9a, it can be observed that the control inputs v r d for both the PI and the proposed controllers are sensitive to the parameter value change occurring at 100 s, leading to spikes in Q s at that moment.

5. Conclusions

This paper presents a novel differentiator-based MPPT controller for WECS equipped with DFIG. The work builds upon existing control algorithms [21,22] and offers several advantages, such as reduced reliance on detailed system information and enhanced robustness to parameter variations and disturbances. Specifically, our approach eliminates the necessity for explicit functional expressions or parameter values in the dynamic equations of the system, focusing instead solely on relative degrees and control directions. Employing HOSD to overestimate the time derivatives of system outputs, our methodology facilitates both a simplified controller architecture and a more straightforward stability analysis. In comparison with existing research, the proposed controller is distinguished by its minimal information requirements, simplified control law structures, systematic design approach, and a reduced number of design constants. Simulation results substantiate the effectiveness of the proposed control algorithm in both tracking maximum power and asymptotically regulating the reactive power to zero. Overall, this contribution offers a more efficient and less complicated approach to MPPT control in WECS.
The differentiator-based controller proposed in this paper is considered to be directly applicable to the MPPT control of permanent magnet synchronous generator WECS, and this is reserved as a topic for future research. One limitation of the proposed controller is that it requires measurements of both the system output and wind speed. Therefore, subsequent research should focus on developing controller strategies that eliminate the need for wind speed sensors.

Funding

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (Ministry of Science and ICT, MSIT) (No. 2021R1A2C1094914).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Cardenas, R.; Pena, R.; Alepuz, S. Overview of control systems for the operation of DFIGs in wind energy applications. IEEE Trans. Ind. Electron. 2013, 60, 2776–2798. [Google Scholar] [CrossRef]
  2. Jadhav, H.T.; Roy, R. A comprehensive review on the grid integration of doubly fed induction generator. Int. J. Electr. Power Energy Syst. 2013, 49, 8–18. [Google Scholar] [CrossRef]
  3. El-Naggar, M.F.; Mosaad, M.I.; Hasanien, H.M.; AbdulFattah, T.A.; Bendary, A.F. Elephant herding algorithm-based optimal PI controller for LVRT enhancement of wind energy conversion systems. Ain Shams Eng. J. 2021, 12, 599–608. [Google Scholar] [CrossRef]
  4. Mossad, M.I.; Sabiha, N.A.; Abu-Siada, A.; Taha, I.B.M. Application of Superconductors to Suppress Ferroresonance Overvoltage in DFIG-WECS. IEEE Trans. Energy Convers. 2021, 37, 766–777. [Google Scholar] [CrossRef]
  5. Hamane, B.; Doumbia, M.L.; Bouhamida, M.; Draou, A.; Chaoui, H.; Benghanem, M. Comparative study of PI, RST, sliding mode and fuzzy supervisory controllers for DFIG based wind energy conversion system. Int. J. Renew. Energy Res. 2015, 5, 1174–1185. [Google Scholar]
  6. Levant, A. Higher-order sliding modes, differentiation and output-feedback control. Int. J. Control 2003, 76, 924–941. [Google Scholar] [CrossRef]
  7. Levant, A.; Livne, M. Weighted Homogeneity and Robustness of Sliding Mode Control. Automatica 2016, 72, 186–193. [Google Scholar] [CrossRef]
  8. Tong, S.; Liu, C.; Li, Y. Fuzzy-adaptive decentralized output-feedback control for large-scale nonlinear systems with dynamical uncertainties. IEEE Trans. Fuzzy Syst. 2010, 18, 845–861. [Google Scholar] [CrossRef]
  9. Tong, S.C.; Li, Y.M.; Zhuang, H.-G. Adaptive neural network decentralized backstepping output-feedback control of nonlinear large-scale systems with time-delays. IEEE Trans. Neural Netw. 2011, 22, 1073–1086. [Google Scholar] [CrossRef]
  10. Li, J.; Chen, W.; Li, J.-M. Adaptive NN output-feedback decentralized stabilization for a class of large-scale stochastic nonlinear strict-feedback systems. Int. J. Robust Nonlinear Control 2011, 21, 452–472. [Google Scholar] [CrossRef]
  11. Tong, S.; Li, Y.; Liu, Y. Adaptive fuzzy output feedback decentralized control of pure-feedback nonlinear large-scale systems. Int. J. Robust Nonlinear Control 2014, 24, 930–954. [Google Scholar] [CrossRef]
  12. Yang, Y.; Yue, D.; Xue, Y. Decentralized adaptive neural output feedback control of a class of large-scale time-delay systems with input saturation. J. Frankl. Inst. 2015, 352, 2129–2151. [Google Scholar] [CrossRef]
  13. Jiang, X.; Mu, X.; Hu, Z. Decentralized Adaptive Fuzzy Tracking Control for a Class of Nonlinear Uncertain Interconnected Systems with Multiple Faults and DoS Attack. IEEE Trans. Fuzzy Syst. 2020, 29, 3130–3141. [Google Scholar] [CrossRef]
  14. Wang, H.; Xiaoping, P.; Bao, J.; Xie, X.-J.; Li, S. Adaptive Neural Output-Feedback Decentralized Control for Large-Scale Nonlinear Systems with Stochastic Disturbances. IEEE Trans. Neural Netw. Learn. Syst. 2020, 31, 972–983. [Google Scholar] [CrossRef] [PubMed]
  15. Ma, Z.; Ma, H. Decentralized Adaptive NN Output-Feedback Fault Compensation Control of Nonlinear Switched Large-Scale Systems with Actuator Dead-Zones. IEEE Trans. Syst. Man Cybern. Syst. 2020, 50, 3435–3447. [Google Scholar] [CrossRef]
  16. Wang, H.; Chen, B.; Lin, C.; Sun, Y. Neural-network-based decentralized output-feedback control for nonlinear large-scale delayed systems with unknown dead-zones and virtual control coefficients. Neurocomputing 2020, 424, 255–267. [Google Scholar] [CrossRef]
  17. Bechlioulis, C.P.; Rovithakis, G.A. A low-complexity global approximation-free control scheme with prescribed performance for unknown pure feedback systems. Automatica 2014, 50, 1217–1226. [Google Scholar] [CrossRef]
  18. Zhang, J.-X.; Yang, G.-H. Low-Complexity Tracking Control of Strict-Feedback Systems with Unknown Control Directions. IEEE Trans. Autom. Control 2019, 64, 5175–5182. [Google Scholar] [CrossRef]
  19. Zhang, J.-X.; Yang, G.-H.; Ding, W. Global Output-Feedback Prescribed Performance Control of Nonlinear Systems with Unknown Virtual Control Coefficients. IEEE Trans. Autom. Control 2022, 67, 6904–6911. [Google Scholar] [CrossRef]
  20. Guo, Z.; Oliveira, T.R.; Guo, J.; Wang, Z. Performance-Guaranteed Adaptive Asymptotic Tracking for Nonlinear Systems with Unknown Sign-Switching Control Direction. IEEE Trans. Autom. Control 2023, 68, 1077–1084. [Google Scholar] [CrossRef]
  21. Park, J.-H.; Kim, S.-H.; Park, T.-S. Approximation-Free Output-Feedback Non-Backstepping Controller for Uncertain SISO Nonautonomous Nonlinear Pure-Feedback Systems. Mathematics 2019, 7, 456. [Google Scholar] [CrossRef]
  22. Park, J.-H.; Kim, S.-H.; Lee, D.-H. Decentralized Output-Feedback Controller for Uncertain Large-Scale Nonlinear Systems Using Higher-Order Switching Differentiator. IEEE Access 2021, 9, 21227–21235. [Google Scholar] [CrossRef]
  23. Wang, J.; Bo, D.; Miao, Q.; Li, Z.; Wu, X.; Lv, D. Maximum power point tracking control for a doubly fed induction generator wind energy conversion system based on multivariable adaptive super-twisting approach. Int. J. Electr. Power Energy Syst. 2021, 124, 106347. [Google Scholar] [CrossRef]
  24. Xiong, L.; Li, P.; Wu, F.; Ma, M.; Khan, M.W.; Wang, J. A coordinated high-order sliding mode control of DFIG wind turbine for power optimization and grid synchronization. Int. J. Electr. Power Energy Syst. 2019, 105, 679–689. [Google Scholar] [CrossRef]
  25. Shihabudheen, K.V.; Raju, S.K.; Pillai, G.N. Control for grid-connected DFIGbased wind energy system using adaptive neuro-fuzzy technique. Int. Trans. Electr. Energy Syst. 2018, 28, e2526. [Google Scholar] [CrossRef]
  26. Kumar, R.; Agrawal, H.P.; Shah, A.; Bansal, H.O. Maximum power point tracking in wind energy conversion system using radial basis function based neural network control strategy. Sustain. Energy Technol. Assess. 2019, 36, 100533. [Google Scholar] [CrossRef]
  27. Aissaoui, H.E.; Ougli, A.E.; Tidhaf, B. Neural Networks and Fuzzy Logic Based Maximum Power Point Tracking Control for Wind Energy Conversion System. Adv. Sci. Technol. Eng. Syst. J. 2021, 6, 586–592. [Google Scholar] [CrossRef]
  28. Sami, I.; Ullah, S.; Amin, S.U.; Al-Durra, A.; Ullah, N.; Ro, J.S. Convergence Enhancement of Super-Twisting Sliding Mode Control Using Artificial Neural Network for DFIG-Based Wind Energy Conversion Systems. IEEE Access 2022, 10, 97625–97641. [Google Scholar] [CrossRef]
  29. Yan, S.; Gu, Z.; Park, J.H.; Xie, X. Adaptive Memory-Event-Triggered Static Output Control of T–S Fuzzy Wind Turbine Systems. IEEE Trans. Fuzzy Syst. 2022, 30, 3894–3904. [Google Scholar] [CrossRef]
  30. Yan, S.; Gu, Z.; Park, J.H.; Xie, X. Sampled Memory-Event-Triggered Fuzzy Load Frequency Control for Wind Power Systems Subject to Outliers and Transmission Delays. IEEE Trans. Cybern. 2023, 53, 4043–4053. [Google Scholar] [CrossRef]
  31. Jing, F.; Liu, W.; Sun, H.; Yu, F.; Wang, K. Prescribed performance finite-time control of wind energy conversion systems with input constraint and system uncertainty. Int. Trans. Electr. Energy Syst. 2021, 32, e13215. [Google Scholar] [CrossRef]
  32. Lv, G.; Cong, S.; Liu, Y.; Peng, Z.; Wang, D. Maximum Power Tracking Control of Wind Energy Conversion Systems Based on Prescribed Performance Function and Extended State Observer. In Proceedings of the 2018 3rd International Conference on Advanced Robotics and Mechatronics (ICARM), Singapore, 18–20 July 2018; pp. 829–834. [Google Scholar] [CrossRef]
  33. Park, J.-H.; Kim, S.-H.; Park, T.-S. Asymptotically Convergent Higher-Order Switching Differentiator. Mathematics 2020, 8, 185. [Google Scholar] [CrossRef]
  34. Boukhezzar, B.; Siguerdidjane, H. Nonlinear control with wind estimation of a DFIG variable speed wind turbine for power capture optimization. Energy Convers. Manag. 2009, 50, 885–892. [Google Scholar] [CrossRef]
  35. Park, J.-H.; Kim, S.-H.; Park, T.-S. Differentiator-Based Output-Feedback Controller for Uncertain Nonautonomous Nonlinear Systems with Unknown Relative Degree. IEEE Access 2020, 8, 172593–172600. [Google Scholar] [CrossRef]
  36. Park, J.-H.; Kim, S.-H.; Park, T.-S. Asymptotically convergent switching differentiator. Int. J. Adapt. Control Signal Process. 2019, 33, 557–566. [Google Scholar] [CrossRef]
  37. Hou, Z.; Chi, R.; Gao, H. An Overview of Dynamic-Linearization-Based Data-Driven Control and Applications. IEEE Trans. Ind. Electron. 2017, 64, 4076–4090. [Google Scholar] [CrossRef]
  38. Hunter, J.D. Matplotlib: A 2D graphics environment. Comput. Sci. Eng. 2007, 9, 90–95. [Google Scholar] [CrossRef]
Figure 1. Concept diagram of the paramters.
Figure 1. Concept diagram of the paramters.
Energies 16 07068 g001
Figure 2. Overall schematic of the DFIG-based generator side control system.
Figure 2. Overall schematic of the DFIG-based generator side control system.
Energies 16 07068 g002
Figure 3. Typical C p curves with respect to pitch angle.
Figure 3. Typical C p curves with respect to pitch angle.
Energies 16 07068 g003
Figure 4. Waveform of the quickly varying wind speed v w ( t ) used in the simulation.
Figure 4. Waveform of the quickly varying wind speed v w ( t ) used in the simulation.
Energies 16 07068 g004
Figure 5. Trajectories of ω ( t ) and ω d ( t ) .
Figure 5. Trajectories of ω ( t ) and ω d ( t ) .
Energies 16 07068 g005
Figure 6. Trajectories of Q s ( t ) .
Figure 6. Trajectories of Q s ( t ) .
Energies 16 07068 g006
Figure 7. Trajectory of C p ( λ , 0 ) .
Figure 7. Trajectory of C p ( λ , 0 ) .
Energies 16 07068 g007
Figure 8. Trajectories of P s ( t ) .
Figure 8. Trajectories of P s ( t ) .
Energies 16 07068 g008
Figure 9. Trajectories of (a) v r d and (b) v r q ( t ) .
Figure 9. Trajectories of (a) v r d and (b) v r q ( t ) .
Energies 16 07068 g009
Table 1. Symbols of DFIG-WECS.
Table 1. Symbols of DFIG-WECS.
NotationDescription
v w ( t ) wind speed
J t inertia of the turbine
J g inertia of the generator
n g gear ratio
Jtotal inertia (= J t + n g 2 J g )
Ddamping constant of the wind turbine
ρ air density
Rradius of the blade
Athe area swept by the blades
ω ( t ) rotational velocity of the wind turbine’s rotor
λ ( t ) tip-speed ratio (= ω R / v w )
λ optimal value of tip-speed ratio λ
β blade pitch angle
C p ( λ , β ) power coefficient function defined as (3)
c 1 , , c 6 system constants in the C p ( λ , β ) function
ω s stator electrical angular speed
R s stator resistance
R r rotor resistance
L s stator inductance
L r rotor inductance
L m mutual inductance
p r number of pole pairs
i r d ( t ) , i r q ( t ) d- and q-axis currents of the generator’s rotor
v r d ( t ) , v r q ( t ) d- and q-axis voltages of the generator’s rotor
ϕ s ( t ) d-axis flux of the generator’s stator
P s ( t ) , Q s ( t ) active and reactive powers
Table 2. Parameters of WECS.
Table 2. Parameters of WECS.
NotationValueDescription
c 1 0.5176constant in (3)
c 2 116constant in (3)
c 3 0.4constant in (3)
c 4 5constant in (3)
c 5 21constant in (3)
c 6 0.0068constant in (3)
ω s 100 π stator electrical angular speed
R s 0.005stator resistance
R r 0.228rotor resistance
L s 0.407stator inductance
L r 0.299rotor inductance
L m 0.0016mutual inductance
p r 4number of pole pairs
J 4.4532 × 10 5 total inertia
D400damping constant
ρ 1.08air density
R35radius of the blade
n g 43.165gear ratio
λ 8.1072optimal value of λ
Table 3. Design parameters of the controllers.
Table 3. Design parameters of the controllers.
NotationValueDescription
L 1 1000constant in (46) and (47)
g 1 20000constant in (44)
k 1 [ 10 6 , 2000 ] T vector in (51)
L 2 1000constant in (60)
g 2 300constant in (58)
k 2 1000constant in (62)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Park, J.-H. Differentiator-Based Output Feedback MPPT Controller for DFIG Wind Energy Conversion Systems with Minimal System Information. Energies 2023, 16, 7068. https://doi.org/10.3390/en16207068

AMA Style

Park J-H. Differentiator-Based Output Feedback MPPT Controller for DFIG Wind Energy Conversion Systems with Minimal System Information. Energies. 2023; 16(20):7068. https://doi.org/10.3390/en16207068

Chicago/Turabian Style

Park, Jang-Hyun. 2023. "Differentiator-Based Output Feedback MPPT Controller for DFIG Wind Energy Conversion Systems with Minimal System Information" Energies 16, no. 20: 7068. https://doi.org/10.3390/en16207068

APA Style

Park, J. -H. (2023). Differentiator-Based Output Feedback MPPT Controller for DFIG Wind Energy Conversion Systems with Minimal System Information. Energies, 16(20), 7068. https://doi.org/10.3390/en16207068

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop