Stability Criteria for Input Filter Design in Converters with CPL: Applications in Sliding Mode Controlled Power Systems

: Microgrids are versatile systems for integration of renewable energy sources and non-conventional storage devices. Sliding Mode techniques grant excellent features of robustness controlling power conditioning systems, making them highly suitable for microgrid applications. However, problems may arise when a converter is set to behave as a Constant Power Load (CPL). These issues manifest in the stability of internal dynamics (or Zero Dynamics), which is determined by the input ﬁlter of the power module. In this paper, a special Lyapunov analysis is conducted to address the nonlinear internal dynamics of SM controlled power modules with CPL. It takes advantage of a Liérnad-type description, establishing stability conditions and providing a secure operation region. These conditions are translated into conductance and invariant region diagrams, turning them into tools for the design of power module ﬁlters.


Introduction
Microgrids can be succinctly understood as an interactive hybrid system that combines distributed generation modules, storage modules and loads. The number and nature of these modules, as well as the microgrid topologies, vary greatly [1][2][3]. Presently, the modules interconnection is predominantly done using AC buses, but DC microgrids are expected to increase in the coming years [4]. In particular, DC microgrids are very versatile for integration of renewable energy sources and non conventional storage devices, such as photovoltaic panels, wind turbines, lithium batteries, flow batteries, supercapacitors, fuel cells and many others (see Figure 1). The increasing development of power electronics and novel control techniques have favoured their implementation, making it possible to obtain higher efficiency and reliability. Furthermore, distribution power systems have been widespread in several applications due to their flexibility features and reduced size and cost.
In this growing scenario, where a wide range of new topologies is being developed, technological challenges are continuously arising, many of them at the level of the modules power converters and their related controllers. The nonlinear behaviour of power converters in these systems often makes conventional linear control techniques not sufficiently effective. In addition, the existence of uncertainties and system disturbances require the development of robust controllers to ensure system stability and efficiency. In this context, Sliding Mode (SM) control have been increasingly accepted for controlling power conditioning systems [5][6][7][8][9]. These control techniques proved to be particularly suitable for electronic converters, granting excellent features of robustness and finite-time convergence, which make them highly appropriate for microgrid applications [10][11][12][13][14]. However, there are certain issues to be aware of at the moment of designing and implementing those power modules. Effectively, it has been extensively reported in the literature that, when a controlled power system tightly regulates energy supplied to a DC bus, the converter behaves as a Constant Power Load (CPL) for the upstream source [15][16][17]. Such CPL acts as a negative impedance, i.e., an increase in its terminal voltage results in a decrease in its current (and vice versa), which undermines the stability margins of the system.
In the case of SM controlled power systems, these issues are manifested in the stability of the Zero Dynamics (ZD). This internal dynamics is determined by the input filter of the power module, which exists to protect the source or the upstream device against switching harmonics propagation. If ZD is not carefully analysed and its stability is not thoroughly addressed during the design of the controlled power system, then performance and robustness could be severely compromised, causing undesired oscillations or even system failures.
Diverse criteria have been proposed for the analysis of input filters dynamics with CPLs in DC microgrids. The most widely found in the literature is the small-signal criterion by eigenvalue analysis through system linearization [17,18]. Also, several methods have been presented using Nyquist approaches, for instance, Gain margin and Phase margin criterion [19], Middlebrook's criterion [20] and the opposing argument criterion [21] among the most important ones. Moreover, in [22] a passivity-based criterion is proposed for a LC filter and in [23] stability conditions are found taking into account parameters variations through bifurcation analysis.
Those methods do not guarantee large-signal stability and give no precise information about the delimitation of a stable region for secure operation of the power module [21]. So as a contribution to this field, in this paper, a special Lyapunov analysis is conducted to address the nonlinear internal dynamics of SM controlled power modules with CPL. This study is performed taking advantage of a transformation to a Liérnad-type description, establishing sufficient conditions to define the stability of the nonlinear system and providing a secure operation region. These conditions are, subsequently, condensed into conductance and invariant region diagrams, which are proposed as versatile tools for the design of power module filters.
The paper is organised as follows. In Section 2, the proposed input filter design method for power converters with CPL is developed. It is divided into two main parts. The first one, Section 2.1, deals with the search of the stability conditions via Liénard-based nonlinear Lyapunov approach. Then, the second part, Section 2.2, translates those conditions into appropriate diagrams, based on which the filter design procedure will be grounded. In Section 3, simulation results are presented and analysed, considering as application case a SM controlled Boost convert with second-order filter. Finally, in Section 4, conclusions and future lines of research are discussed.

Development of the Proposed Input Filter Design Method for Power Converters with CPL
Two steps are required to develop the proposed method for input filter design. The first one is to obtain sufficient conditions for the zero dynamics (ZD) stability, over the full nonlinear operation range of SM controlled converters with CPL. Based on those conditions, the second step consists of elaborating a filter design criteria with the help of ad hoc conductance diagram description.
In Figure 2, a schematic diagram showing the topology of the power module under study is presented. In the following subsections, both aforementioned steps are treated in detail.

Search of Stability Conditions via Liénard Based Nonlinear Lyapunov Approach
In this subsection, the ZD is thoroughly studied and sufficient conditions are obtained by mean of nonlinear Lyapunov analysis. For the sake of clarity, a path of progressive complexity is followed. Initially, a nonlinear Lyapunov approach is applied to a simple topology based on a first-order capacitive filter, and a sufficient condition for its stability is obtained. Then, the study is broadened to a widely-used second-order LC filter, describing the system in a Liénard form and analyzing its stability through a special energy-like Lyapunov function.

First-Order Capacitive Filter Topology
The electrical circuit of the capacitive filter supplying a constant power P f is shown in Figure 3. Its electrical model is defined by the first order differential equation: where C f is the filter capacitance, V oc is the open circuit source voltage, R s is the internal source resistance and v the voltage across the capacitor. The nonlinear term P f v , corresponding to the CPL current, acts as a negative impedance, when the voltage v is increased the filter output current is decreased. This negative impedance reduces the equivalent resistance of the system, leading to a reduction of the stability margin. From (1), the well-known condition for the existence of equilibrium in CPL is [24]:

Voltage source
which if strictly satisfied, results in two solutions. The first one, is the equilibrium capacitor voltage v 0 , that defines the equilibrium point of operation of the power module. The second one, has no physical meaning, and represents a voltage limit for the system stability. Additionally, by multiplying (3) by (4), it is straightforward to obtain v lim = P f R s v 0 (the operation equilibrium, given by v 0 , and the stability limit, determined by v lim , can be visualised further in Figure 4, at the end of this subsection). As previously discussed, traditional stability study of the system described by (1) has been addressed by linearization, determining the local stability near the equilibrium voltage v 0 through eigenvalue analysis. In accordance with this conventional analysis, if P f satisfies: then the eigenvalues of the linearised system around the equilibrium point have negative real part and so the trajectories of the system converge to the equilibrium voltage v 0 , for close enough initial conditions. It is well known that this linear analysis only ensures the stability for trajectories sufficiently close to v 0 , but gives no information about the attraction region. In fact, it will be shown that a power system whose eigenvalues analysis results are stable, would in practice become unstable if voltage v drops below the stability limit voltage v lim .
Effectively, through the nonlinear Lyapunov analysis presented in the sequel, it is possible to completely define the system stability, allowing establishing the actual range of v where secure operation is guaranteed.
To this end, in this simple first case study, the equilibrium point is translated to the origin through the state transformation: it leads to:ż = −a(z) with Now, the proposed energy-like Lyapunov function is locally positive-definite if condition is satisfied for some interval that includes the origin. Therefore, it follows that the first derivative of V is equal to:V which is definite-negative for a(z) = 0. The Lyapunov function, defined in (9), is shown in Figure 4, with z 0 = 0 and z lim = v lim − v 0 obtained from (3) and (4), respectively, in combination with (6). As can be appreciated, the local maximum and minimum of the function V(z) are at z lim and z 0 (corresponding to capacitor voltages v = v lim and v = v 0 , respectively). These values of z are zeros of the functionV(z), i.e., zeros of a(z), and thus zeros ofż (see (7)). Stability condition in terms of conductances. It can straightforwardly obtained by computing inequality (10) as a function of the voltage v: where g 0 is an equivalent output conductance viewed from the filter and g s is the internal conductance of the voltage source.
Stable operation region S for first-order filter. It should be remarked that, unlike condition (5) obtained via eigenvalue linear analysis, the condition given by (12) defines the stable region of the first-order filter, fully determining the admissible range for the capacitor voltage v. It establishes that if the voltage v satisfies: all system trajectories converge to the equilibrium voltage v 0 . Otherwise, the system trajectories become unstable.

Second-Order LC Filter Topology
Founded on the previous approach, in this subsection, the Zero Dynamics stability of the power module is studied for widespread second-order input filters. To this end, nonlinear Lyapunov analysis is therefore employed to obtain stability conditions in terms of equivalent conductances. Then, these conditions will help to determine an invariant region where secure operation of power modules will be completely guaranteed.
Regretfully, this is not as direct as it is for the first-order capacitive filter. To be able to perform the analysis, in this work, the ZD is firstly rewritten into a special Liénard-type form [25], by means of a linear transformation. This allows a suitable energy-like Lyapunov function capable of dealing with this nonlinear topology. Figure 5 shows a schematic model of the second-order LC filter connected to the converter acting as a constant power load.
Voltage source Figure 5. Scheme of the equivalent electrical circuit of the second-order input filter with CPL.
The system is described by: The above system presents an operation equilibrium point, defined by both the equilibrium source current i s0 and (if (2) holds) the equilibrium filter voltage v 0 : where the latter is the solution with physical meaning for the power system operation, whereas there exists a second solution which again represents a voltage limit for the system stability.
Analogously to the previous case, the equilibrium point is shifted to the origin by redefining the system states as: leading to the transformed system: To obtain a ZD representation convenient for the Lyapunov analysis, the following state linear transformation is proposed: and consequently, the transformed state system acquires the Liénard-type form [25]: Now, it can be proposed the appropriate energy-like Lyapunov function, depending on new states z 1 and z 2 : with V(0, 0) = 0. Then, similar to condition (10), if z 2 satisfies: the integral term will be positive, and in consequence, V(z 1 , z 2 ) will be positive-definite for some region around the origin.
Next, the first derivative of V is straightly computed as: which means thatV will be locally negative-definite if it satisfies a(z 2 )b(z 2 ) > 0 for some interval around the origin. Therefore, taking into account condition (28), it leads to: Hereunder, the Lyapunov conditions (28) and (31) will be rewritten in terms of equivalent conductances, preparing the ground for the completion of the proposed design method in the following subsection.
First Lyapunov condition in terms of conductances. It is obtained combining condition (28) with (25), giving: 1 where similarly to (12), g 0 is the equivalent output conductance viewed from the filter and g s is the internal conductance of the source. Second Lyapunov condition in terms of conductances. It is computed from condition (26) and (31), resulting in: where the new equivalent conductance g lc is called filter conductance.
Furthermore, both conditions can be summarised in one inequality as: or, similarly to the first-order filter, in terms of voltage v, as: However, unlike the first-order case, the above condition (alternatively (36) or (33)-(35)) is not sufficient to ensure the ZD stable behaviour, because not all trajectories will necessarily converge to (i s0 , v 0 ) for any initial conditions v > v min . In fact, with a second-order filter, the stable region also depends on the inductor current i s .
Nevertheless, (37) does guarantee the existence of an invariant region around the equilibrium point, within which all system trajectories would converge to the desired equilibrium. Moreover, in contrast to the eigenvalues approach, following the subsequent analysis this stable region can be fully determined and eventually used as a helping tool for the filter design.
Stable operation region S for second-order filter. Condition (37) is therefore used to define such invariant region in the ZD state plane, where performance and secure operation of the SM controlled power system can be ensured.
For better interpretation, Figure 6 schematically depicts the Lyapunov function V(z 1 , z 2 ), including both equilibrium solutions, these are the the desired equilibrium point of operation, p 0 = (i s0 , v 0 ) in green, and the limit point p lim = (i s0 , v lim ) in red. The rationale behind the proposed procedure is quite simple. A basin of attraction targeting to p 0 must be devised by intersecting V(z 1 , z 2 ) with an appropriate horizontal plane, Plane invar .
(b) Lyapunov function V(z 1 , z 2 ). Lateral view. It can be easily proved that the limit point p lim is an unstable equilibrium, hence it defines a plane, Plane lim (given by V(z 1 , z 2 ) = V lim , with constant V lim = V(p lim )), which naturally is an upper bound for Plane invar . Moreover, Plane lim delimited a basin-shaped, yet not necessarily attractive, Lyapunov function.
Therefore, to ensure attractiveness to the desired equilibrium, the lower (or equal) Plane invar should be obtained such that, below it, Lyapunov conditions hold (or, equivalently, condition (37) holds). Then, the proposed invariant region (in blue solid line) is delimited by the intersection of V(z 1 , z 2 ) with Plane invar , defined by V(z 1 , z 2 ) = V invar (see Figure 6a). The constant value V invar will be computed below using v min , such thatV results negative for all (z 1 , z 2 ) : V(z 1 , z 2 ) < V invar , thus the basin of attraction is attained and all ZD trajectories which start below Plane invar converge to p 0 . Therefore, the stable invariant region can be formally defined where the constant V invar is computed according to (27): with Please note that the way v min = max v lim ; P f L f R s C f v 0 is defined in (37), effectively implies that Plane lim is the upper limit plane for the Plane invar .

Proposed Filter Design Criteria through Conductance Diagram Description
The second step to conclude the method for input filter design is developed in this subsection. It is based on translating the previously obtained stability conditions into appropriate conductance diagrams. These diagrams, complemented with those of the stable operation region of the power module, are then proposed as practical tools for systematic design.

Design Criterion for the First-Order Input Filter
The case of the power module with a first-order input filter is relatively simple and can be straightforwardly designed with the assistance of condition (12) (or, equivalently, (13)). Therefore, the main contribution of this case to the paper is not the design itself, but the introduction and interpretation of the conductance diagrams, familiarising the reader with their use in the framework of a filter design procedure, aiming to the more interesting second-order topology to come.
For the case under study, the condition arrived in (12) provides a relation between the output conductance g 0 and the source conductance g s for the stability of the ZD system described in (1). This relation can be compared for different values of the output power P f in a conductance diagram, where a stable operation region can be established. Figure 7 presents the conductance diagram for a voltage source with open-circuit voltage V oc = 24 V and internal resistance R s = 0.144 Ω. The diagram is constructed using the conditions (12)-(13), obtained from the procedure described in Section 2.1.1. Each hyperbola (dashed black lines) represents an output conductance g 0 for a specific value of power P f . For the sake of illustration, the conductance g 0 for a given power P f = P f nom is highlighted in thick black line. The upper limit curve corresponds to the maximum allowed power P f max , computed through (2).
In addition, the conductance g 0 at each equilibrium point (green line) is obtained as: The intersection between the source conductance g s (horizontal red line) and each g 0 defines the limit voltage v lim , which establishes the secure operation of the power system (the locus constituted by these points for different P f corresponds to the solid red segment). Then, the region of stable operation is delimited by the boundary conformed by the latter and the maximum allowed power, depicted in solid red line. All trajectories that start inside this region (white background in the figure) will converge to the equilibrium voltage v 0 (green line). Otherwise, they will become unstable.
Please note that for the first order filter case the conductance diagram is defined by the source parameters R s and V s , thus the stability region will not be altered by the selection of the filter capacitance C f . Therefore, from the above, the design procedure for this case study is direct. Firstly, knowing the desired cut-off frequency f c for the filter and R s , capacitance C f is computed.
Then, from the stability region of the conductance diagram, for a desired power P f = P f nom specified in the design, the admissible range of v can be obtained. Alternatively, for a known admissible range of variation of the capacitor voltage v, an upper limit for the P f nom of the power module can be set. It can be appreciated in the conductance diagram that exists a trade-off between the selection of the operation power P f nom and the admissible variation range of v.

Design Criterion for Second-Order Input Filter
Building on the previous approach, in this subsection the design criterion for more prevalent second-order filter in converters with CPL is presented. As it was mentioned in Section 2.1.2, unlike the first-order case, the condition v > v min given by (37) is necessary but not sufficient for the design. Therefore, the conductance diagram remains as a helpful design tool to determine and easily visualise the viable stable operation points. However, to ensure a design that provides wholly secure operation of the power module, it must be complemented with an invariant operation region diagram. This invariant region is constructed through the proposed procedure presented in Section 2.1.2 by means of (38)-(41). The value of the constant V invar , which defines the stable region, is computed from (39), using diverse values of filter power P f and filter parameters.
To facilitate the understanding of those diagrams for this topology, a descriptive analysis is presented below, together with illustrative diagrams of conductance and of the invariant region for different values of P f , assuming fixed cut-off frequency f c and filter capacitor C f .
As with the first-order case, the conductance diagram shown in Figure 8a is constructed through the set of Equations (32)-(37). However, now, the extra condition defined by (35), which strongly depends on the filter parameters, also needs to be considered. , for different P f . However, in this case, the red solid line defined by the source conductance g s and P f max is not a stability boundary. In fact, the conductance condition in this topology, set by inequality (36), depends not only on g s , but also on the equivalent conductance g lc (dashed blue line). So, a more restrictive boundary (in thick solid blue line) arises, conformed by min(g s ; g lc ) and a critical power P f crit . The latter is an upper bound for P f , beyond which the equilibrium points of the green locus are unstable. It can be computed by Please note that the white-background region delimited by the blue boundary, does define the stable equilibria, hence it allows the visualisation of the admissible values of P f that can be selected as operation filter power.
Nevertheless, it does not provide per se a region that secures stable operation, as it was for the first-order filter. Thus, to take into account in the design procedure the admissible variation ranges of v and i s , the assistance of the invariant region diagram is required. In Figure 8b invariant regions S are plotted, parametrised for different values of filter power P f , with their boundaries computed as V(z 1 , z 2 ) = V invar in accordance with (38)-(41).
The admissible ranges of v and i s strongly depend on the selected operation power P f . Note the existing trade-off for the design. The higher the selected P f , getting closer to the P f crit , the lower the invariant region, reducing the stability tolerance to admissible variations of v and i s with respect to the nominal operation point.
In this framework, a sequential design procedure will be succinctly presented below, as an illustrative example of application.

•
Firstly, the desired cut-off frequency f c is chosen to protect the source against switching harmonics propagation. It should be carefully selected, considering that a small value of f c could lead to large components sizing, which would imply a higher cost and weight. On the other hand, there is an upper bound for f c given by the switching converter frequency, whose practical higher limit is typically selected one decade below such switching frequency. In this case study, a f c value equals to 1 kHz is taken. • Next, the operation power P f is selected, in accordance with the power module requirements. In this example, the second-order input filter is designed for supplying an electrical power up to P f = 750 W (solid black line in Figure 9a).

•
The goal of the following step is to obtain the range of selectable values for capacitor C f . These values are those which guarantee stable operation points for the f c and P f previously chosen. This step is easily fulfilled using the conductance diagram in Figure 9a, where g lc is plotted for different values of capacitor C f (blue horizontal lines).
It can be observed in Figure 9a that there exists a minimum value of capacitor C f , equals to 0.65 mF, for which the equilibrium point of P f = 750 W (black dot) would become unstable. So, designing C f greater than this so-called C f min ensures that the stable equilibria region (white background zone) will contain the desired equilibrium point.
Therefore, from the conductance diagram, the range of selectable values for the filter capacitor can be set as C f > 0.65 mF = C f min . If required, this value can be analytically computed through (35), as a function of the cut-off frequency: • Then, from those admissible values of C f , one should be chosen, according to the required variation ranges of v and i s , which are contemplated in the specifications of the power module. This design step is conducted with the help of the invariant region diagram in Figure 9b, parametrised in terms of C f . To illustrate the selection procedure, variations up to ±20% for both voltage v and current i s , are considered admissible for the normal operation of the power module (represented by the green area Figure 9b). Now, the proposed design criterion is simple. The smallest C f should be selected, such that its corresponding stable operation invariant region includes the admissible variation green area. Therefore, for this case study, C f = 850 µF is chosen (see curve in blue line).
Please note that there is a trade-off between the capacitor sizing and the allowed voltage v and current i s variations. In this particular case, all capacitances larger than 850 µF would provide secure stable operation to the power module.

•
Finally, once the required C f was selected and knowing the desired f c , the inductance L f can be easily obtained from: which for this illustrative example, results L f = 30 µH.
To conclude, it should be remarked that the previous procedure is just one of diverse possible ways to perform the filter design, which may differ depending on each particular case, available data and specifications.

Application Case: SM Controlled Boost Converter with Second-Order Filter
In this section, following the above proposed criteria, the ZD stability of a second-order input filter for a SM controlled boost converter-based power module is evaluated. An illustrative electrical diagram of the considered system is presented in Figure 10. The power module is intended to operate as a part of a microgrid, as the one presented in Figure 1, connected to a common DC Bus, in this case of 48 V. In addition, the power module is fed by the power source v s (i s ), characterized by its open circuit voltage and internal resistance, which may represent different kinds of sources (i.e., Li-Ion batteries, Flow batteries, Fuel Cells, Supercapacitors, etc.).

L b R b
Power source Figure 10. Electrical diagram of the PCS based on boost converter with second-order input filter.
In this particular case study, the power converter is assumed to be controlled by means of a Second Order Sliding Mode (SOSM) algorithm, designed to regulate a constant power to the DC bus (see details in Appendix A).
Under this CPL assumption, the filter is design following the criterion explained in Section 2.2.2. Then, in the sequel, to validate the proposed design method, the power system is evaluated under different initial conditions and filter power P f .
Recapitulating, a cut-off frequency f c = 1 kHz and an operation filter power P f = 750 W were selected. Next, the proposed design provided a filter capacitor of C f = 850 µF, considering admissible voltage and current variations up to ±20%. Lastly, L f was computed using (45), giving L f = 30 µH.
The complete set of system parameters employed for the subsequent simulations results are displayed in Table 1. It is important to remark that with this set of parameters, the resulting critical power P f crit is of about 930 W. If the system exceeds this power value, the equilibrium points become unstable (see Figure 9a). Then, to evaluate the system behaviour, several in-silico test were conducted. In the first place, the system is forced to operate under different values of constant power, until the critical power P f crit is reached. This results are shown in Figure 11. Figure 11a presents the filter power evolution for such stepped power reference. It can be seen that the delivered power results equal to the reference, until the critical power P f crit is surpassed (∼0.75 s). At this moment, the ZD becomes unstable and the filter power collapses to zero. Figure 11b depicts the corresponding filter voltage and current. Please note that as the the filter power increases, the v and i s variations during transients become bigger. Once the critical power is surpassed, the voltage drops to zero and the current saturates at its maximum (i s = V oc /R s ).
The following set of tests aims to show that region S, obtained via the proposed design, is an invariant region and that the inside trajectories converge to the desired equilibrium point. To this end, Figure 12 displays the trajectories in the i s − v plane for several initial conditions (red dots) placed near the boundary of S (750 W) (blue curve). Effectively, as was previously discussed, it can be appreciated that this region ensures ZD stability of the SM controlled module, guaranteeing its secure operation.  The final test presented in Figure 13 is intended to illustrate the unstable behaviour for an initial condition outside the invariant region S(750 W). The latter is plotted in blue line in the i s − v plane (top left figure in Figure 13). It can be observed in this example, as the ZD stability is not guaranteed, the system trajectory is not converging to equilibrium. It is important to stress that not all possible initial conditions outside this invariant region S result in unstable trajectories.
It is worth noticing that the SOSM-STA manages to reach the power reference of 750 W only for about 0.5 ms (see top right figure in Figure 13. Power reference in red line). As the ZD under this conditions results unstable, the controller is not able to remain over the sliding surface for a long time. Hence, the filter voltage v drops to zero while the filter current saturates at its maximum (i s = V oc /R s ), leading to a system failure.

Conclusions
This paper proposed design criteria for input filters operating in DC microgrid applications with CPL. The presented methods were aimed to ensure the stable operation of power modules controlled by Sliding Mode Control techniques. To this end, conductance and invariant region diagrams were proposed as practical tools for filter design. The use of these diagrams represents a handy visual assistance for the selection of the filter elements.
To develop this procedure, firstly, the Zero Dynamics of the SM controlled module were brought into a particular Liérnad-type form, through a special transformation. Then, stability conditions, in term of equivalent conductances, were obtained using a suitable energy-like Lyapunov function capable to deal with the ZD nonlinear structure. Finally, the proposed conductance diagram and invariant region were presented, allowing establishing a sequential filter design procedure.
With the purpose of validating the proposed methodology, a SM controlled Boost converter with a second-order filter was assessed by simulation. The output power of the module was assumed regulated by means of a Second Order Sliding Mode Super-Twisting Control. The module tests consisted on leading the system to different power condition until reaching critical power condition, where the system becomes unstable. Also, it was proved that, for initial conditions outside the proposed stable region, the operation of the SM controlled system was not guaranteed.
Future studies will be focused on extending the proposed methodology to nonlinear input power source. In addition, in the next stage of this research, to evaluate the design criteria against real unmodelled dynamics and perturbation, the implementation of experimental tests will be undertaken.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A. Second Order Sliding Mode Power Control
As indicated in Section 3, a SOSM controller is used in this paper, although First Order Sliding Mode (FOSM) controllers could also be considered. In particular, SOSM algorithms have been widely accepted to control power system, due to their robust features while providing chattering amelioration [26]. It is not the objective of this paper to describe the controller design in depth, but an outline is provided in this Appendix.
A well-known dynamical averaged model of the power module based on boost converter with second-order low-pass filter (see Figure 10) can be described by: The power module controller is designed to regulate the delivered electrical power P f to the DC bus.
To developed the SM controller, firstly, the control objective must be defined through the sliding variable σ as: Then, the SOSM control objective will be accomplished when system trajectories reach and remain on the second-order sliding surface σ =σ = 0.
Please note that from (A3), the sliding variable σ has relative degree one with respect to the duty cycle u. Then, a Super-Twisting Algorithm (STA) is employed, because of its direct application to relative degree one systems.
The SOSM-STA controller structure is defined as: where α and β are the STA gains. They are computed from σ(x, u) = ϕ(x, u) + γ(x) · u (A6) by bounding functions ϕ(x, u) and γ(x) (see [26,27] for more details). Then, once the SOSM controlled system reaches σ = 0, the resulting ZD will be the one described in Section 2.1.2.