Optimal Tuning of Fractional Order Controllers for Dual Active Bridge-Based DC Microgrid Including Voltage Stability Assessment

In this article, three evolutionary search algorithms: particle swarm optimization (PSO), simulated annealing (SA) and genetic algorithms (GA), have been employed to determine the optimal parameter values of the fractional-order (FO)-PI controllers implemented in the dual active bridge-based (DAB) DC microgrid. The optimum strategy to obtain the parameters of these FO-PI controllers is still a major challenge for many power systems applications. The FO-PI controllers implemented in the DAB are used to control the DC link voltage to the desired value and limit the current flowing through the converter. Accordingly, the investigated control system has six parameters to be tuned simultaneously; Kp1, Ki1, λ1 for FO-PI voltage controller and Kp2, Ki2, λ2 for FO-PI current controller. Crucially, this tuning optimization process has been developed to enhance the voltage stability of a DC microgrid. By observing the frequency-domain analysis of the closed-loop and the results of the subsequent time-domain simulations, it has been demonstrated that the evolutionary algorithms have provided optimal controller gains, which ensures the voltage stability of the DC microgrid. The main contribution of the article can be considered in the successful application of evolutionary search algorithms to tune the parameters of FO-based dual loop controllers of a DC microgrid scheme whose power conditioner is a DAB topology.


Introduction
The dual active bridge (DAB) converter is considered as one of the most critical emerging topologies that have been successfully used in many electrical power applications such as DC microgrids [1,2]. The DAB offers key advantages such as a minimum size of passive filter components and high-power density when compared with other bidirectional DC-DC converters [3]. Therefore, DAB converter has gained popularity during the last decade in applications such as electric vehicles (EVs), energy management systems or DC microgrids with energy storage systems (ESS) [4][5][6][7][8].
A general average model of a DAB has been proposed in [2], whereas a simple but reliable model has been proposed in [7]. In Reference [7], the DAB converter average and small-signal models have been validated through both simulations and experimental tests for the four different control modes. The DABs have also been the object of study in [9][10][11][12][13][14][15][16][17][18][19][20]. Even though several control strategies have been explored for the DAB, the most common among practitioners is the simple phase shift modulation technique [17,18] based on conventional PI controllers [20,21]. Besides, the variable structure controllers have been employed to control the DAB's operation, as can be seen in [22].
Given the proliferation of the DC and mixed AC/DC microgrids, some research has been carried out to evaluate the capabilities of some power electronic converters for contribution of the proposed methodology. Finally, the main conclusions of the paper are summarized in Section 9.

Literature Review
Several research activities have been focused on the control techniques of the DAB during the past few years, see [1][2][3][4][5][6][7][8][9][10][11][12][13]. A comprehensive overview of different modulation strategies of DAB converters is presented in [1]. In Reference [2], an HFT DAB topology including an efficiency optimization-based modulation strategy is proposed. In Reference [3], state-feedback techniques and linear matrix inequalities technique have been employed to control DAB converter operation. The proposed approach provides a fast transient response and improves the robustness against load and voltage variations. A hybrid modulation strategy for DAB converter that improves the efficiency over a wide range of input voltage and load variations has been proposed in [4,5]. In Reference [6], a per-unit model of DAB converter has been developed, where the article also addressed an optimized control technique based on the triple-phase shift (TPS) technique. The authors of [7] proposed a simple phase-shift modulation technique based on the parabolic carrier concept that has been implemented digitally through Verilog and hardware digital circuits. Investigators in [8,9] suggested several modulation methods for the DAB converter to avoid transformer saturation.
In References [10,11], a model predictive control (MPC) approach has been employed to control the operation of DAB converter. While in [12,13] the DAB control scheme is based on direct power control concept. In References [14][15][16], the current modulation schemes enhance the transient response of the DAB and reduce the low order harmonics in the current waveform. A bi-directional power flow scheme has been proposed in [17], where the power flow between the AC main grid and DC microgrid is controlled through a DAB converter. In References [18,19], the reactive power control of a DAB converter is assessed. Reference [20] focused on reducing the current stress. While a power flow scheme has been proposed in [21] to control the power transfer between two microgrids.
Many efforts have been carried out to improve DAB converter performance for many purposes. In this vein, the authors of [37][38][39][40] have achieved considerable DAB converter modelling improvements.

Application of FO-PI Controllers
Indeed, the FO controllers have gained interest in recent years [41][42][43][44][45][46][47][48][49][50]. In Reference [41], a survey of state-of-the-art advances in FO control specifically for time-delay systems provided. In Reference [42], a detailed design and tuning of robust FO controller employed for autonomous microgrid system are presented. In Reference [43], FPGA has been utilized to implement the FO controllers, while the authors of [44] utilized DSP to implement the FO controller. Investigators of [45] presented a tuning method and a robustness analysis of a FO controller.
Similarly, the optimal design of an FO controller to tune a pneumatic servo system has been studied in [46]. In References [47,48], the FO controllers have been utilized for closed-loop control of permanent magnet synchronous motor drives. In Reference [49], a design procedure and experimental validation of a FO controller have been used to control hydraulically powered actuator operation. The FO controllers have been employed in [50] to control the back-to-back operation of some power converters.

Tuning of FO-PI Controllers
When the parameters of fractional-order PI controllers are optimally selected, the overall control system exhibits better dynamic and steady-state performance when compared with the traditional PI controllers, see [46,47]. However, two extra parameters (i.e., the integrator and differentiator orders) result in a more complicated tuning procedure. Thence, tuning of the fractional-order PI controller parameters is still of interest, see [51][52][53][54]. In Reference [51], the tuning of an FO-PI controller has been carried out using a genetic algorithm.
On the other hand, [52] proposed and derived unified fractional complex order controller expressions as a universal structure. Reference [53] presented an atom search optimization algorithm and its chaotic version to determine the optimal parameters of the FO controller. In Reference [54], hybrid PSO optimization and grey wolf optimization methods to tune FO controller parameters are used.

Conventional PID Controllers
The standard form of a conventional PID controller is expressed as Equation (1a) and represented in Figure 1, where a generalized block diagram of a single-input single-output closed-loop control system is illustrated. As is widely known, the PI controller is obtained by removing the derivative part in Equation (1a): netic algorithms [57][58][59][60], ant colony [61][62][63], particle swarm optimization [64][65][66][67][68] ogeography-based optimization [69] can be employed to search for the optimum of FO controller parameters by formulating a multi-variable optimization problem

Conventional PID Controllers
The standard form of a conventional PID controller is expressed as Equation ( represented in Figure 1, where a generalized block diagram of a single-input single closed-loop control system is illustrated. As is widely known, the PI controller is o by removing the derivative part in Equation (1a): Accordingly, the transfer function of the conventional PID controller is descri where Kp is the proportional gain, Ki is the integrator gain, Kd is the differentiator g is the system error, and u(t) is the controller output. The PID controller is still on preferable and commonly used controllers in many industrial applications [46, cently, FO controllers have been introduced to overcome the major drawbacks of t ventional PID controllers [41,45,46,49].

FO PID Controllers
Equation (3) describes the transfer function of a FO order PID controller: where λ and μ are the order and differential terms, respectively. Compared with c tional integer-order PID controller, FO controller has integral and differential or ferent than unity such that (λ and μ ≠ 1). In fact, FO controller can be consider generalized form of conventional integer-order PID controller, where the conve Accordingly, the transfer function of the conventional PID controller is described by: where K p is the proportional gain, K i is the integrator gain, K d is the differentiator gain, e(t) is the system error, and u(t) is the controller output. The PID controller is still one of the preferable and commonly used controllers in many industrial applications [46,49].
Recently, FO controllers have been introduced to overcome the major drawbacks of the conventional PID controllers [41,45,46,49].

FO PID Controllers
Equation (3) describes the transfer function of a FO order PID controller: where λ and µ are the order and differential terms, respectively. Compared with conventional integer-order PID controller, FO controller has integral and differential order different than unity such that (λ and µ = 1). In fact, FO controller can be considered as a generalized form of conventional integer-order PID controller, where the conventional PID controller is achieved by setting λ and µ to unity. Other controller functions (i.e., special cases such as P, PI or PD controllers) can be obtained by assigning λ and µ to suitable values as illustrated in Figure 2, where four distinct values are indicated through the (λ-µ) plane. In several research articles, it has been reported that FO controllers provide better response time and reduce the overshoot compared with the conventional PID controllers, see for instance [42] and [47]. However, the tuning process of an FO-PID controller is far more complex because there are many parameters to be optimized simultaneously (i.e., K p , K i , K d , λ and µ). In this paper, the derivative part is removed, and therefore the two controllers are FO-PI.

Detailed Test System
The detailed model of the DC microgrid object of study is shown in Figure 4 ing to the circuit diagram depicted in Figure 4, the whole system is compose stages. The first stage is the three-phase uncontrolled rectifier which feeds the D link. The second stage comprises a single-phase bridge AC/DC inverter who feeds the HFT transformer. In third place, the second single-phase bridge AC/D feed the DC link. As stated above, the DC bus voltage control system comprises trol loops (i.e., the outer loop for the voltage and the inner loop for the current)

Block Diagram
The block diagram of the investigated system is illustrated in Figure 3, where the DC microgrid is connected to the DAB output terminal. The voltage at the DAB output terminals is used as input for the DAB voltage regulation. The DAB converter is connected to the grid through a 3-phase six-pulse uncontrolled diode rectifier. The DC bus of the microgrid is loaded with both constant impedance and constant power load models. The control unit of the DAB converter is based on a cascaded loop with the implemented FO-PI controllers. Those parameters are tuned in an offline process by the abovementioned evolutionary search algorithms; see the procedure in Section 5.
Kp, Ki, Kd, λ and μ). In this paper, the derivative part is removed, and therefore the two controllers are FO-PI.

Block Diagram
The block diagram of the investigated system is illustrated in Figure 3, where the DC microgrid is connected to the DAB output terminal. The voltage at the DAB output terminals is used as input for the DAB voltage regulation. The DAB converter is connected to the grid through a 3-phase six-pulse uncontrolled diode rectifier. The DC bus of the microgrid is loaded with both constant impedance and constant power load models. The control unit of the DAB converter is based on a cascaded loop with the implemented FO-PI controllers. Those parameters are tuned in an offline process by the abovementioned evolutionary search algorithms; see the procedure in Section 5.

Detailed Test System
The detailed model of the DC microgrid object of study is shown in Figure 4. According to the circuit diagram depicted in Figure 4, the whole system is composed of three stages. The first stage is the three-phase uncontrolled rectifier which feeds the DC voltage link. The second stage comprises a single-phase bridge AC/DC inverter whose output feeds the HFT transformer. In third place, the second single-phase bridge AC/DC rectifier feed the DC link. As stated above, the DC bus voltage control system comprises two control loops (i.e., the outer loop for the voltage and the inner loop for the current).

Detailed Test System
The detailed model of the DC microgrid object of study is shown in Figure 4. According to the circuit diagram depicted in Figure 4, the whole system is composed of three stages. The first stage is the three-phase uncontrolled rectifier which feeds the DC voltage link. The second stage comprises a single-phase bridge AC/DC inverter whose output feeds the HFT transformer. In third place, the second single-phase bridge AC/DC rectifier feed the DC link. As stated above, the DC bus voltage control system comprises two control loops (i.e., the outer loop for the voltage and the inner loop for the current).

Test System Components
This section details the model implemented in each component of the test system. The test system illustrated in Figure 4 can be divided into several elements as follows: (i) the grid and the three-phase diode rectifier, (ii) the DAB converter and (iii) the microgrid loads.

Main Grid and Three-Phase Rectifier
The Thevenin's equivalent represents the impedance of the grid seen by the threephase diode rectifier. Therefore, the AC grid is modelled through an impedance (Rac and Lac). A coupling impedance in the DC side of the rectifier is also added (Rdc and Ldc).
The AC voltage source and the three-phase diode rectifier can be expressed as a DC source (Eg) with an equivalent impedance (RS and LS) as follows [70]: where Vg is the line-to-tine voltage of the three-phase AC system, and the last term of the equivalent series resistance in the DC side represents the effect of the commutation in the voltage drop. It is worth noting that in Figure 4, the equivalent resistance and inductance seen by the three-phase rectifier and the coupling resistance and inductance are not displayed. The data related to these parameters is provided in Appendix A.

Main Grid and Three-Phase Rectifier
In this subsection, the mathematical expressions of the DAB are recounted. The detailed model of a DAB has been proposed in [5]. Nevertheless, the average and smallsignal models are detailed in [7]. As can be seen, the DAB is formed by two single-phase

Test System Components
This section details the model implemented in each component of the test system. The test system illustrated in Figure 4 can be divided into several elements as follows: (i) the grid and the three-phase diode rectifier, (ii) the DAB converter and (iii) the microgrid loads.

Main Grid and Three-Phase Rectifier
The Thevenin's equivalent represents the impedance of the grid seen by the threephase diode rectifier. Therefore, the AC grid is modelled through an impedance (R ac and L ac ). A coupling impedance in the DC side of the rectifier is also added (R dc and L dc ).
The AC voltage source and the three-phase diode rectifier can be expressed as a DC source (E g ) with an equivalent impedance (R S and L S ) as follows [70]: where V g is the line-to-tine voltage of the three-phase AC system, and the last term of the equivalent series resistance in the DC side represents the effect of the commutation in the voltage drop. It is worth noting that in Figure 4, the equivalent resistance and inductance seen by the three-phase rectifier and the coupling resistance and inductance are not displayed. The data related to these parameters is provided in Appendix A.

Main Grid and Three-Phase Rectifier
In this subsection, the mathematical expressions of the DAB are recounted. The detailed model of a DAB has been proposed in [5]. Nevertheless, the average and smallsignal models are detailed in [7]. As can be seen, the DAB is formed by two single-phase bridges with a complementary switching cycle, where two high-frequency square-wave signals are generated. The power flow regulation across the DAB can be achieved by a shift between these two voltage signals [71][72][73][74]. These two square-wave-phase-shifted voltages create a voltage across this leakage inductance, making a current circulate through it. Thence, the current flowing through the inductor is expressed as: where V 1 and V 2 are the two square-wave voltages generated by the two bridges, L k is the leakage inductance and i lk the current that flows through this inductance. For each half period of the voltage signals, the relationship between voltages and currents is expressed as: where I 1 and I 2 are the two currents corresponding to the leakage currents when the leading and lagging bridges are switched, n is the turns ratio of the HFT (1:n), d is the phase shift angle between the two bridge voltages, and T is half of the switching period. The switching period depends on the switching frequency f s of the converter. Therefore Equation (6a,b) can also be computed as: The average DAB output current is expressed by combining Equation (6a,b) [36] i o = 1 nT where t 1 and t 2 are the sum of these two times are equal to the phase-shift for a half period. Terms t 1 , t 2 , I 1 and I 2 are given by: where V i and V o are the DAB input and output voltages. Both the input and output average currents are represented as i i and i o and are computed as follows: The average-model DAB converter is displayed in Figure 5, where the two current sources are represented by i i and i o .

+
where Ro and Co are the DC microgrid load and the capacitance of the DAB output, respectively. Note that Co is designed as C2 in Figure 4. As shown in Equation (12), the term Vo can be affected by changes in the DAB input voltage Vi, changes in the phase-shift d, and changes in the microgrid load composition. If the microgrid load is modelled as CPL, the term Ro can be rewritten as:  The small-signal model is obtained by perturbing and linearizing the averaged model of Equation (10a,b) for a particular point of operation and is defined by Equation (11a,b). The small-signal model is depicted in Figure 6.
where A, B, C and D are the four current sources of the small-signal stability model observed in Figure 6, and D is the phase shift duty ratio between the square-wave signals of both bridges. These current sources are identified as g id , g ivo , g od , g ovi and expressed by Equation (11c) through Equation (11f). See the detailed procedure to obtain these values in [36]. Given Equation (11a,b), the response of the output voltage can be predicted by a first-order model as in Equation (12), where the DAB capacitance and the load resistance are expressed as unique impedance.
where R o and C o are the DC microgrid load and the capacitance of the DAB output, respectively. Note that C o is designed as C 2 in Figure 4. As shown in Equation (12), the term V o can be affected by changes in the DAB input voltage V i , changes in the phase-shift d, and changes in the microgrid load composition. If the microgrid load is modelled as CPL, the term R o can be rewritten as:

Problem Formulation
As stated earlier, the FO controllers implemented in the DAB converter control illustrated in Figure 7 have six parameters to be tuned simultaneously. The first FO-PI controller belongs to the outer loop, which generates the current reference for the inner loop. The inner loop generates a reference for the phase-shift angle dϕ. The three evolutionary algorithms have been used to determine the optimum values of the FO-PI controller parameter values according to the objective function described in Equation (14). This objective function has to be minimized:

Particle Swarm Optimization
The PSO algorithm is an optimization technique capable of finding globally optimal solutions using the adjacent social agent interaction. In PSO, a flock of particles flies through a search space with velocity updated by inertia movement, self-cognition, and social interaction [64]. The velocity and position of each particle is updated at each iteration by Equations (15) and (16).
where xi is the position of the i-th particle in the search space, vi is the velocity of the i-th particle, w is the inertia of particles, α and β are uniformly distributed positive random vector. L is the particle best position, G is the global best position achieved by all particles, while k is the instantaneous number of each iteration [66]. In this article, the parameters used in the PSO algorithm are presented in Table 1. The flowchart is shown in Figure 8.

Problem Formulation
As stated earlier, the FO controllers implemented in the DAB converter control illustrated in Figure 7 have six parameters to be tuned simultaneously. The first FO-PI controller belongs to the outer loop, which generates the current reference for the inner loop. The inner loop generates a reference for the phase-shift angle d φ .

Problem Formulation
As stated earlier, the FO controllers implemented in the DAB converter control illustrated in Figure 7 have six parameters to be tuned simultaneously. The first FO-PI controller belongs to the outer loop, which generates the current reference for the inner loop. The inner loop generates a reference for the phase-shift angle dϕ. The three evolutionary algorithms have been used to determine the optimum values of the FO-PI controller parameter values according to the objective function described in Equation (14). This objective function has to be minimized:

Particle Swarm Optimization
The PSO algorithm is an optimization technique capable of finding globally optimal solutions using the adjacent social agent interaction. In PSO, a flock of particles flies through a search space with velocity updated by inertia movement, self-cognition, and social interaction [64]. The velocity and position of each particle is updated at each iteration by Equations (15) and (16).
where xi is the position of the i-th particle in the search space, vi is the velocity of the i-th particle, w is the inertia of particles, α and β are uniformly distributed positive random vector. L is the particle best position, G is the global best position achieved by all particles, while k is the instantaneous number of each iteration [66]. In this article, the parameters used in the PSO algorithm are presented in Table 1. The flowchart is shown in Figure 8.  The three evolutionary algorithms have been used to determine the optimum values of the FO-PI controller parameter values according to the objective function described in Equation (14). This objective function has to be minimized:

Particle Swarm Optimization
The PSO algorithm is an optimization technique capable of finding globally optimal solutions using the adjacent social agent interaction. In PSO, a flock of particles flies through a search space with velocity updated by inertia movement, self-cognition, and social interaction [64]. The velocity and position of each particle is updated at each iteration by Equations (15) and (16).
where x i is the position of the i-th particle in the search space, v i is the velocity of the i-th particle, w is the inertia of particles, α and β are uniformly distributed positive random vector. L is the particle best position, G is the global best position achieved by all particles, while k is the instantaneous number of each iteration [66]. In this article, the parameters used in the PSO algorithm are presented in Table 1. The flowchart is shown in Figure 8.

Simulated Annealing
In the SA algorithm, the process of annealing in metals [55,56] is emulated. In the annealing process, heat is applied to a solid until it melts. Then, melted solid is allowed to cool slowly in a controlled manner. The solid reaches the minimum energy state, provided that the initial temperature is sufficiently high, and the cooling time is sufficiently long [56]. The SA algorithm requires an initial temperature, final temperature and cooling rate. An objective function is customized to represent the energy of a fictitious material. Temperature T is a control parameter, which decreases gradually as the algorithm goes further from the initial step to the subsequent steps. The criterion of accepting a neighborhood solution n from the current solution c is given by the following probability function described by (17):

Simulated Annealing
In the SA algorithm, the process of annealing in metals [55,56] is emulated. In the annealing process, heat is applied to a solid until it melts. Then, melted solid is allowed to cool slowly in a controlled manner. The solid reaches the minimum energy state, provided that the initial temperature is sufficiently high, and the cooling time is sufficiently long [56]. The SA algorithm requires an initial temperature, final temperature and cooling rate. An objective function is customized to represent the energy of a fictitious material. Temperature T is a control parameter, which decreases gradually as the algorithm goes further from the initial step to the subsequent steps. The criterion of accepting a neighborhood solution n from the current solution c is given by the following probability function described by (17): The temperature decay (decreasing) scheme in its simplest form will follow (18): where T k+1 is the temperature at step k + 1, T k is the temperature at step k, and 0 ≤ α 1 ≤ 1.
The SA parameters are detailed in Table 2 and the flowchart of this algorithm in Figure 9.

Genetic Algorithm
GA is one of the stochastic optimization search techniques that rely on natural selection. The concept of GA has been borrowed from the evolution of biology [57,58]. The process of creating the next generation from the current population toward achieving a globally optimized solution is composed of three main steps: selection, crossover and mutation. The selection rules select the individuals, called parents, that contribute to the next generation's population, while crossover rules combine two parents to form children for the next generation [60]. The GA repeatedly modifies the population of individual solutions. At each step, the genetic algorithm selects individuals at random from the current population to be parents and uses them to produce the children for the next generation. Over successive generations, the population evolves toward an optimal solution based on the best fittedindividuals' fitness. The parameters of the used GA for the optimal controllers optimal tuning are summarized in Table 3, whilst the flowchart is illustrated in Figure 10.

Stability Analysis of the DAB-Based DC Microgrid
This section evaluates the stability of the test system under investigation. For that purpose, the small-signal model of the DAB has been used. In Section 6.1, the fundamental concepts used to evaluate the stability in power converters are recalled, and the stability analysis of the microgrid object of study is provided. Finally, the results of this stability assessment are detailed and discussed in Section 6.2.

Stability Criteria and Impedance Modelling
In general, stability is assured if the ratio between (Z out //Z in ) meets the Nyquist criterion [36]. Nevertheless, the phase of the resulting impedance Z out //Z in of the bode diagram must be between −90 degrees and +90 degrees at all frequencies. If the resulting parallel impedance phase is between −90 degrees and +90 degrees, the average power at the input port of the network is positive. Therefore, the system consumes energy (i.e., it is a passive system). If this impedance phase is exactly +90 degrees or −90 degrees, the average power is zero, and the system is lossless. If the phase is less than −90 degrees or greater than +90 degrees, the average power is negative, and the system injects energy (i.e., it is an active system). Therefore, the stability of a particular closed-loop system is guaranteed if the Nyquist diagram has no right half-planes (RHPs) and the contour of this diagram does not encircle the point (−1, jω) of this plane [36].
By considering the test system shown in Figure 4, the impedance seen by the DAB input terminals is designed as Z out and is represented by the grid impedance and the equivalent coupling RLC filter (4b) and (4c). The term Z out can be characterized by R S , L S and C 1 : where R S and L S are the equivalent resistance and inductances, and C 1 is the DAB input capacitor. By linearizing and perturbing the DAB averaged model, the small-signal model can be obtained [26]. The small-signal variables are denoted by uppercase letters. As can be seen in Figure 7, the phase shift angle is obtained through the closed-loop as: where the subscript * denotes a reference, G C (s) and G V (s) are the inner and outer loop transfer functions, respectively. Accordingly, I * o is the reference value of DC bus current I * DC BUS previously shown in Figure 7 and given in Equation (14). Similarly, V * o is the reference value of DC bus voltage V * DC BUS previously shown in Figure 7 and given in Equation (14) as well.
By rearranging Equation (20a,b), the following expression is obtained by: Thus, the phase-shiftd between the square-wave voltage signals of the two bridges can be expressed as a function of (21): whereV o andî o are the output voltage and current of the DAB, whereas R o is the DC load. Thereby, the voltages and currents of the DAB of the model in Figure 6, can be expressed as Equation (11a,b), where the DAB input currentî i is as follows: where V o is the rated voltage at DAB output, n is the transformer turns ratio, L k is the leakage inductance, and D is the phase shift duty ratio. The term T represents a half period, and f S is the switching frequency. If the DAB input admittance C 1 is considered part of Z out in (19), the closed-loop input impedance is obtained through: The representation of the impedance model is shown in Figure 11.
riod, and fS is the switching frequency. If the DAB input admittance C1 is considered part of Zout in (19), the closed-loop input impedance is obtained through: The representation of the impedance model is shown in Figure 11.

Discussion of the Stability Analysis
In this section, the results of the frequency-domain analysis are presented. Firstly, the small-signal system stability is evaluated considering the output impedances, which is the most unfavorable scenario for stability purposes. After that, another case study is carried out considering the grid as an ideal voltage source, and as a consequence, the DAB input voltage at DAB is considered stable. Therefore, according to these two assumptions, the first considers the line and DC coupling impedances (i.e., RDC, LDC, RAC and LAC). On the contrary, the second assumes an ideal DC source voltage where these impedances are neglected.
The stability can be demonstrated by overlapping Zin and the four Zout impedances in the same bode plot. In Figure 12, it is seen that the input impedance of the DAB is above all these output impedances, and the phase margin is below 90 degrees. The values of these four impedances are summarized in Table 4. The larger the resistor and capacitor, the higher are the gain margin and the more stable is the system.  Figure 11. Impedance model.

Discussion of the Stability Analysis
In this section, the results of the frequency-domain analysis are presented. Firstly, the small-signal system stability is evaluated considering the output impedances, which is the most unfavorable scenario for stability purposes. After that, another case study is carried out considering the grid as an ideal voltage source, and as a consequence, the DAB input voltage at DAB is considered stable. Therefore, according to these two assumptions, the first considers the line and DC coupling impedances (i.e., R DC , L DC , R AC and L AC ). On the contrary, the second assumes an ideal DC source voltage where these impedances are neglected.
The stability can be demonstrated by overlapping Z in and the four Z out impedances in the same bode plot. In Figure 12, it is seen that the input impedance of the DAB is above all these output impedances, and the phase margin is below 90 degrees. The values of these four impedances are summarized in Table 4. The larger the resistor and capacitor, the higher are the gain margin and the more stable is the system. To verify the previous explanation, the Nyquist diagram is displayed in Figure 13, where it can be seen that the four scenarios meet the stability criterion. As expected from the results in Figure 12, the most stable scenario occurs with Zout1, which is the closest to the origin in Figure 13. On the contrary, the less stable case occurs when Zout4 is considered.
Once the stability of the DAB with the three-phase rectifier and the main grid has been performed, the rest of this section aims to show the performance of the parameters obtained with the three algorithms. To that purpose, as stated earlier, the main grid is now considered an ideal source and therefore, Vin can be regarded as stable. The Nyquist diagram of the closed-loop is illustrated in Figure 14, whereas the bode plot is shown in Figure 15. From Figure 14, it can be seen that all controllers will provide a stable response. Nevertheless, the relatively most stable results are expected when the controllers are tuned with the GA. On the contrary, the parameters obtained with the PSO algorithm may result in a relatively least robust response. The bode diagram depicted in Figure 15     To verify the previous explanation, the Nyquist diagram is displayed in Figure 13, where it can be seen that the four scenarios meet the stability criterion. As expected from the results in Figure 12, the most stable scenario occurs with Z out1 , which is the closest to the origin in Figure 13. On the contrary, the less stable case occurs when Z out4 is considered.
functions according to the current and voltage controllers are displayed in Figures 16 and  17, respectively.   It is worth mentioning that the parameters in both the inner and outer loops are those obtained with the GA (i.e., K p1 = 33.55, K i1 = 0.673, λ 1 = 0.217; K p2 = 4.89; K i2 = 18.24, λ 2 = 0.32).
Once the stability of the DAB with the three-phase rectifier and the main grid has been performed, the rest of this section aims to show the performance of the parameters obtained with the three algorithms. To that purpose, as stated earlier, the main grid is now considered an ideal source and therefore, V in can be regarded as stable. The Nyquist diagram of the closed-loop is illustrated in Figure 14, whereas the bode plot is shown in Figure 15. From Figure 14, it can be seen that all controllers will provide a stable response. Nevertheless, the relatively most stable results are expected when the controllers are tuned with the GA. On the contrary, the parameters obtained with the PSO algorithm may result in a relatively least robust response. The bode diagram depicted in Figure 15

Simulation Parameters
Specifications and simulation parameters of the investigated system have been summarized in Appendix A in Table A1. One set of FO-PI controllers parameters obtained from each evolutionary search algorithms is summarized in Table 5. The simulation study is carried out based on these values. The overall DC microgrid system has been implemented in the PSIM software package, where the load models are displayed in Figure18.

Loading Conditions and Testing Scenarios
The time-domain simulations carried out with the detailed microgrid model implemented in PSIM can be divided into two major case studies. The first one tests the system's response for a step-change in the voltage reference considering a resistive load. On the other hand, the second case study analyzes a sudden load change. Within the second case study, three types of voltage-dependent load models have been considered (i.e., constant impedance, constant power, and two types of single-phase inverters feeding R and RL loads) as illustrated in Figure 18. The modelled single-phase inverters have been designed to work in two modes. The first inverter operates in square wave mode, whereas the second is a bipolar SPWM.
The responses obtained for each set of parameters obtained according to the three evolutionary algorithms are discussed and analyzed in Sections 7.3, 7.4 and 7.5, for several load models. The optimization data is analyzed in Section 7.6, whilst a quantitative analysis of the DC microgrid performance is summarized in Table 6.

Simulation Parameters
Specifications and simulation parameters of the investigated system have been summarized in Appendix A in Table A1. One set of FO-PI controllers parameters obtained from each evolutionary search algorithms is summarized in Table 5. The simulation study is carried out based on these values. The overall DC microgrid system has been implemented in the PSIM software package, where the load models are displayed in Figure 18.

System Response considering a Resistive Load
Initially, the microgrid is loaded with a constant resistive load of 20 Ω. The voltage reference is set to 100 V, and at t = 0.1 s, the reference is changed to 200 V. The response of

Loading Conditions and Testing Scenarios
The time-domain simulations carried out with the detailed microgrid model implemented in PSIM can be divided into two major case studies. The first one tests the system's response for a step-change in the voltage reference considering a resistive load. On the other hand, the second case study analyzes a sudden load change. Within the second case study, three types of voltage-dependent load models have been considered (i.e., constant impedance, constant power, and two types of single-phase inverters feeding R and RL loads) as illustrated in Figure 18. The modelled single-phase inverters have been designed to work in two modes. The first inverter operates in square wave mode, whereas the second is a bipolar SPWM.
The responses obtained for each set of parameters obtained according to the three evolutionary algorithms are discussed and analyzed in Sections 7.3-7.5, for several load models. The optimization data is analyzed in Section 7.6, whilst a quantitative analysis of the DC microgrid performance is summarized in Table 6.

System Response Considering a Resistive Load
Initially, the microgrid is loaded with a constant resistive load of 20 Ω. The voltage reference is set to 100 V, and at t = 0.1 s, the reference is changed to 200 V. The response of the DC bus voltage and the corresponding DC bus current are illustrated in Figures 19a, 20a and 21a    In addition, the voltage response obtained for the three algorithms is depicted in Figure 22, where the three evolutionary search algorithms result in a similar transient response. Moreover, the peak overshoot is less than 3% for the three algorithms. Therefore, this event cannot yet be used as a determining factor to select the best tuning strategy.
As shown in the comparison displayed in Figure 23, the controllers tuned by GA and SA cause a voltage dip of 3.75 V and 4.5 V, respectively. Meanwhile, the PSO leads to a larger voltage dip, 7.5 V. Consequently, it has been demonstrated that the relatively best dynamic response is achieved with the parameters tuned by the GA in terms of voltage stability and steady-state error. In contrast, the worst response is obtained with the parameters tuned by the PSO.  Secondly, the microgrid is subjected to a sudden load change, where the load resistance is increased from 20 Ω to 10 Ω. Meanwhile, the DC bus voltage reference is kept fixed to 200 V. The load disturbance occurs at t = 0.15 s as shown in Figures 19b, 20b and 21b for PSO, SA and GA, respectively.
The voltages displayed in these figures indicate that the FO-PI controllers tuned by GA and SA algorithms provide a smoother response than those obtained by the PSO algorithm. The percentage of steady-state errors during the microgrid load change are 0.7%, 0.9%, and 1.26% with FO-controllers tuned by GA, SA, and PSO algorithms, respectively.
In addition, the voltage response obtained for the three algorithms is depicted in Figure 22, where the three evolutionary search algorithms result in a similar transient response. Moreover, the peak overshoot is less than 3% for the three algorithms. Therefore, this event cannot yet be used as a determining factor to select the best tuning strategy. In addition, the voltage response obtained for the three algorithms is depicted in Figure 22, where the three evolutionary search algorithms result in a similar transient response. Moreover, the peak overshoot is less than 3% for the three algorithms. Therefore, this event cannot yet be used as a determining factor to select the best tuning strategy.
As shown in the comparison displayed in Figure 23, the controllers tuned by GA and SA cause a voltage dip of 3.75 V and 4.5 V, respectively. Meanwhile, the PSO leads to a larger voltage dip, 7.5 V. Consequently, it has been demonstrated that the relatively best dynamic response is achieved with the parameters tuned by the GA in terms of voltage stability and steady-state error. In contrast, the worst response is obtained with the parameters tuned by the PSO. Figure 22. Comparison between the three algorithms given a step-change in voltage reference. Figure 22. Comparison between the three algorithms given a step-change in voltage reference.
As shown in the comparison displayed in Figure 23, the controllers tuned by GA and SA cause a voltage dip of 3.75 V and 4.5 V, respectively. Meanwhile, the PSO leads to a larger voltage dip, 7.5 V. Consequently, it has been demonstrated that the relatively best dynamic response is achieved with the parameters tuned by the GA in terms of voltage stability and steady-state error. In contrast, the worst response is obtained with the parameters tuned by the PSO.

System Response considering a Constant Power Load
In this operating scenario, the DC microgrid is loaded with a constant power load with an initial power of 1 kW. At t = 0.3 s, the load power is abruptly changed to 2 kW.

System Response Considering a Constant Power Load
In this operating scenario, the DC microgrid is loaded with a constant power load with an initial power of 1 kW. At t = 0.3 s, the load power is abruptly changed to 2 kW. The corresponding system responses (i.e., the DC voltage and current) are plotted in Figure 24 for the FO-controllers tuned by PSO, SA, and GA algorithms, respectively.

System Response considering a Constant Power Load
In this operating scenario, the DC microgrid is loaded with a constant power load with an initial power of 1 kW. At t = 0.3 s, the load power is abruptly changed to 2 kW. The obtained results prove that the FO controllers provide a robust response in terms of voltage stability when a constant power load increase is considered. If the FO-PI controllers tuned by PSO are analyzed, the obtained voltage drop is 1.5 V, representing a 0.75% deviation. Meanwhile, the voltage dips are 1.3 V and 0.3 V for the SA and GA, respectively. As can be observed, the FO-PI controllers tuned by the GA provide the relatively best response in terms of voltage regulation.

Response with Single-Phase Inverters
The investigated microgrid system has been tested considering AC loads feed by single-phase inverters. In practice, these inverters can be used in many applications (e.g., including pumping systems, ventilators and air conditioning systems). In this study, two operating scenarios have been presented. In the first scenario, the inverter operates in square-wave mode, producing a fixed output voltage waveform. In the second operating scenario, the inverter is working in a bipolar-SPWM mode, where the output voltage of the inverter can be controlled by adjusting the amplitude modulation index.

Square-Wave Single-Phase Inverter
In this subsection, the response of the DC bus voltage is illustrated with the three sets of FO-PI controller parameters. Initially, the inverter is loaded with a resistance of 60 , which at = 0.25 s it is reduced to 20 . The inverter operating frequency is 60 Hz (i.e., the fundamental frequency).
The response of the microgrid has been illustrated in Figures 25-27 with the controllers' parameters tuned by PSO, SA and GA, respectively. The voltage drop during the sudden load variation proved to be 198.4 V if the controller parameters tuned by the PSO are considered, which is a deviation of 0.8% with respect to the rated DC voltage.
In Figure 25a, the inverter output current is plotted, while the resultant bus voltage is plotted in Figure 25b. The inverter output voltage is also displayed in Figure 25c. The obtained results prove that the FO controllers provide a robust response in terms of voltage stability when a constant power load increase is considered. If the FO-PI controllers tuned by PSO are analyzed, the obtained voltage drop is 1.5 V, representing a 0.75% deviation. Meanwhile, the voltage dips are 1.3 V and 0.3 V for the SA and GA, respectively. As can be observed, the FO-PI controllers tuned by the GA provide the relatively best response in terms of voltage regulation.

Response with Single-Phase Inverters
The investigated microgrid system has been tested considering AC loads feed by single-phase inverters. In practice, these inverters can be used in many applications (e.g., including pumping systems, ventilators and air conditioning systems). In this study, two operating scenarios have been presented. In the first scenario, the inverter operates in square-wave mode, producing a fixed output voltage waveform. In the second operating scenario, the inverter is working in a bipolar-SPWM mode, where the output voltage of the inverter can be controlled by adjusting the amplitude modulation index.

Square-Wave Single-Phase Inverter
In this subsection, the response of the DC bus voltage is illustrated with the three sets of FO-PI controller parameters. Initially, the inverter is loaded with a resistance of 60 Ω, which at = 0.25 s it is reduced to 20 Ω. The inverter operating frequency is 60 Hz (i.e., the fundamental frequency).
The response of the microgrid has been illustrated in Figures 25-27 with the controllers' parameters tuned by PSO, SA and GA, respectively. The voltage drop during the sudden load variation proved to be 198.4 V if the controller parameters tuned by the PSO are considered, which is a deviation of 0.8% with respect to the rated DC voltage. nics 2021, 10, x FOR PEER REVIEW 23 of 36  If the SA algorithm tune the FO-PI controllers, the voltage drop caused by the load change is 198.9 V, that is, a 0.55% decrease. Figure 26a,b show the inverter current and the DC bus voltage, respectively. Moreover, the inverter voltage is illustrated in Figure 26c.
If the obtained FO-PI controller parameters obtained by the GA are considered, the bus voltage decreased to 199.7 V during the sudden load variation, representing a voltage  If the SA algorithm tune the FO-PI controllers, the voltage drop caused by the load change is 198.9 V, that is, a 0.55% decrease. Figure 26a,b show the inverter current and the DC bus voltage, respectively. Moreover, the inverter voltage is illustrated in Figure 26c.
If the obtained FO-PI controller parameters obtained by the GA are considered, the bus voltage decreased to 199.7 V during the sudden load variation, representing a voltage  drop of 0.15% with respect to the rated value. The inverter output current is presented in Figure 27a. On the other hand, the resultant bus voltage is plotted in Figure 27b, and the inverter output voltage is displayed in Figure 27c. According to the obtained results, the FO-controller parameters tuned by the GA algorithm exhibit the minimum voltage drop during the load resistance change if the singlephase inverter operates in square-wave mode. In contrast, the relatively worst result has been achieved with PSO algorithm.

SPWM Single-Phase Inverter
The response of the DC bus voltage is illustrated in this section for the three sets of FO-PI controller parameters considering an SPWM single-phase inverter. Initially, the inverter is loaded with an RL load whose resistive and inductive values are 60  and 2 mH, respectively. Afterwards, at t = 0.25 s, the resistive term is reduced to 20 . The inverter operating frequency is 60 Hz, the switching frequency is 15 kHz and the amplitude modulation index 0.6.
The response of the microgrid has been illustrated in Figures 28-30 with controllers' parameters tuned by PSO, SA and GA, respectively.
During the transient caused by the sudden load connection, a negligible overvoltage has been identified with the parameters tuned by the PSO algorithm (200.3 V). Whereas, during the same event, the minimum recorded voltage value proved to be 198.6 V. The maximum peak-to-peak voltage ripple is 1.7 V, representing a 0.85% drop with respect to the rated voltage. In Figure 28a, the inverter output current is plotted, while the resultant Bus voltage is plotted in Figure 28b. In addition, the corresponding spectral decomposition of the DC bus voltage is plotted in Figure 28c. According to the obtained results of the spectral decomposition, the ripple of the DC bus voltage has a low-order harmonic content with negligible magnitudes. In Figure 25a, the inverter output current is plotted, while the resultant bus voltage is plotted in Figure 25b. The inverter output voltage is also displayed in Figure 25c.
If the SA algorithm tune the FO-PI controllers, the voltage drop caused by the load change is 198.9 V, that is, a 0.55% decrease. Figure 26a,b show the inverter current and the DC bus voltage, respectively. Moreover, the inverter voltage is illustrated in Figure 26c.
If the obtained FO-PI controller parameters obtained by the GA are considered, the bus voltage decreased to 199.7 V during the sudden load variation, representing a voltage drop of 0.15% with respect to the rated value. The inverter output current is presented in Figure 27a. On the other hand, the resultant bus voltage is plotted in Figure 27b, and the inverter output voltage is displayed in Figure 27c.
According to the obtained results, the FO-controller parameters tuned by the GA algorithm exhibit the minimum voltage drop during the load resistance change if the single-phase inverter operates in square-wave mode. In contrast, the relatively worst result has been achieved with PSO algorithm.

SPWM Single-Phase Inverter
The response of the DC bus voltage is illustrated in this section for the three sets of FO-PI controller parameters considering an SPWM single-phase inverter. Initially, the inverter is loaded with an RL load whose resistive and inductive values are 60 Ω and 2 mH, respectively. Afterwards, at t = 0.25 s, the resistive term is reduced to 20 Ω. The inverter operating frequency is 60 Hz, the switching frequency is 15 kHz and the amplitude modulation index 0.6.
The response of the microgrid has been illustrated in Figures 28-30 with controllers' parameters tuned by PSO, SA and GA, respectively. When the parameters tuned by the SA are considered, a negligible overvoltage has been measured during the transient caused by the load connection, see Figure 28a. Whereas, during the same event, the minimum voltage value proved to be 199.2 V. The maximum peak-to-peak voltage ripple is 1.5 V, representing a 0.75% drop with respect to the rated voltage. Figure 29a shows that the inverter output current is plotted, while the resultant DC bus voltage is plotted in Figure 29b.
In addition, the corresponding spectral decomposition of the DC bus voltage is plotted in Figure 29c. According to the obtained results of this decomposition, it is seen that the ripple of the DC bus voltage has a low-order harmonic content with negligible magnitudes. In particular, a second-order harmonic content with a magnitude of 0.109 V is observed. When the tuned parameters are obtained with the GA algorithm, the maximum overvoltage observed during the transient caused by the sudden load connection is 203 V (see Figure 30a). On the contrary, during the same event, the minimum voltage value proved to be 198.1 V. The maximum peak-to-peak voltage ripple is 4.9 V, representing a 2.45% drop with respect to the rated voltage. Accordingly, this gain is intentionally reduced to 30 and the response has been investigated again. The archived results of the DC bus voltage, current and, the corresponding spectral decomposition are plotted again in Figure 31a,b respectively.
Considering this change in the value of Kp1, the maximum detected bus voltage is 201.7 V, while the minimum detected value is 198.4 V. Thus, the peak-to-peak voltage ripple became 3.3 V instead of 4.9 V.
According to resultant spectral decomposition, the DC bus voltage has a low order harmonic content. Thus, a slight second-order harmonic component with 0.077 V magnitude is observed. During the transient caused by the sudden load connection, a negligible overvoltage has been identified with the parameters tuned by the PSO algorithm (200.3 V). Whereas, during the same event, the minimum recorded voltage value proved to be 198.6 V. The maximum peak-to-peak voltage ripple is 1.7 V, representing a 0.85% drop with respect to the rated voltage. In Figure 28a, the inverter output current is plotted, while the resultant Bus voltage is plotted in Figure 28b. In addition, the corresponding spectral decomposition of the DC bus voltage is plotted in Figure 28c. According to the obtained results of the spectral decomposition, the ripple of the DC bus voltage has a low-order harmonic content with negligible magnitudes.
When the parameters tuned by the SA are considered, a negligible overvoltage has been measured during the transient caused by the load connection, see Figure 28a. Whereas, during the same event, the minimum voltage value proved to be 199.2 V. The maximum peak-to-peak voltage ripple is 1.5 V, representing a 0.75% drop with respect to the rated voltage. Figure 29a shows that the inverter output current is plotted, while the resultant DC bus voltage is plotted in Figure 29b.
In addition, the corresponding spectral decomposition of the DC bus voltage is plotted in Figure 29c. According to the obtained results of this decomposition, it is seen that the ripple of the DC bus voltage has a low-order harmonic content with negligible magnitudes. In particular, a second-order harmonic content with a magnitude of 0.109 V is observed.
When the tuned parameters are obtained with the GA algorithm, the maximum overvoltage observed during the transient caused by the sudden load connection is 203 V (see Figure 30a). On the contrary, during the same event, the minimum voltage value proved to be 198.1 V. The maximum peak-to-peak voltage ripple is 4.9 V, representing a 2.45% drop with respect to the rated voltage. In Figure 30a, the inverter output current is plotted, while the resultant Bus voltage is plotted in Figure 30b. In addition, the corresponding spectral decomposition of the dc bus voltage is plotted in Figure 30c.
It has been observed that if an SPWM inverter is selected to fed AC loads through the DC link, the worst achieved results in terms of peak-to-peak voltage ripple are those tuned by the GA algorithm, see the plots in Figure 30.
Considering the above, the proportional gain of the voltage controller is relatively large (i.e., Kp1 33.55) compared with the corresponding values achieved by the two other algorithms (i.e., Kp 1 with PSO is 9.4 and Kp 1 with SA is 16.3).
Accordingly, this gain is intentionally reduced to 30 and the response has been investigated again. The archived results of the DC bus voltage, current and, the corresponding spectral decomposition are plotted again in Figure 31a,b respectively.

Results of the Optimization Process
The principal comparison items in this section are the following; the convergence time of each algorithm and the steady-state value of the objective function. The evolution of fitness function has also been considered as a critical factor for comparison purposes. Accordingly, the results of the convergence are provided in Figure 32, where it is seen that all algorithms converge to the same minimum final value. However, the GA algorithm converges to the final solution after 10 iterations, the PSO algorithm after 18 iterations, while the SA algorithm's convergence is achieved in 40 iterations. The time of convergence for each algorithm is summarized in Table 6. Although the GA algorithm converges to the optimal solution with fewer iterations than the other considered algorithms, the overall convergence time is the largest when compared with the PSO and SA. Considering this change in the value of Kp 1 , the maximum detected bus voltage is 201.7 V, while the minimum detected value is 198.4 V. Thus, the peak-to-peak voltage ripple became 3.3 V instead of 4.9 V.
According to resultant spectral decomposition, the DC bus voltage has a low order harmonic content. Thus, a slight second-order harmonic component with 0.077 V magnitude is observed.

Results of the Optimization Process
The principal comparison items in this section are the following; the convergence time of each algorithm and the steady-state value of the objective function. The evolution of fitness function has also been considered as a critical factor for comparison purposes. Accordingly, the results of the convergence are provided in Figure 32, where it is seen that all algorithms converge to the same minimum final value. However, the GA algorithm converges to the final solution after 10 iterations, the PSO algorithm after 18 iterations, while the SA algorithm's convergence is achieved in 40 iterations. The time of convergence for each algorithm is summarized in Table 6. Although the GA algorithm converges to the optimal solution with fewer iterations than the other considered algorithms, the overall convergence time is the largest when compared with the PSO and SA.
Accordingly, the results of the convergence are provided in Figure 32, where it is seen that all algorithms converge to the same minimum final value. However, the GA algorithm converges to the final solution after 10 iterations, the PSO algorithm after 18 iterations, while the SA algorithm's convergence is achieved in 40 iterations. The time of convergence for each algorithm is summarized in Table 6. Although the GA algorithm converges to the optimal solution with fewer iterations than the other considered algorithms, the overall convergence time is the largest when compared with the PSO and SA.  Crucially, it has been found that SA is the fastest algorithm; meanwhile, GA is the one the takes more time for each iteration. The elapsed time for each iteration in the GA algorithm is between 10 s and 12 s, between 4 s and 6 s for the PSO and between 0.5 s and 1 s for the SA.

Impact of Parameters Uncertainty on the DC Bus Voltage
The impact of parameters variation on the DC bus voltage has been studied and summarized in this subsection. Besides load variation, four uncertain parameters have been considered in this article (leakage inductance L k , inductor equivalent series resistance R k , DAB input capacitance C 1 , and DAB output capacitance C 2 ). The impact of parameters variation is studied considering a single-phase inverter operating in square-wave mode with an AC load.
The inductance (L k ) is subjected to ± 30% variation from the nominal value (170 µH). The resultant voltage dips under the extreme values have been listed in Table 7. The results indicate that the L k uncertainty has a negligible effect on the DC bus voltage. Secondly, as shown in Table 8, two additional ESR (R k ) values have been considered (i.e., 0.05 Ω and 1 Ω). The resultant voltage drops with these different ESR values are summarized in Table 8. The results indicate that the R k has a negligible effect on the DC bus voltage. Thirdly, the DAB capacitances (C 1 ) and (C 2 ) have been subjected to a ±25% variation with respect to their rated value, that is, 400 µF. The resultant voltage drops for these case studies are provided in Tables 9 and 10. The results indicate that both C 1 and C 2 have no noticeable effect on the DC bus voltage. The results summarized in Tables 6-10 indicate that the microgrid load is the most dominant uncertain parameter; meanwhile, as mentioned above, the other passive parameters of the DAB have negligible impact on the DC bus voltage of the microgrid.

Comparison with Other Existing Techniques
This article has analyzed three evolutionary algorithms to tune the FO-PI controllers of a DAB for voltage stability purposes in a DC microgrid. The following comparison highlights the contribution of our assessment:

•
In the available literature, several recent articles have used evolutionary algorithms to obtain optimal parameters in many applications, see [42][43][44][45][46][47][48]. However, any of them has been used to tune FO-PI controllers for voltage stability purposes in DC microgrids. Moreover, most of them have only used one algorithm. On the contrary, in the present paper, a comparison between three different algorithms is provided.

•
Although the recently published articles towards voltage stability in DC microgrids have considered innovative approaches to address several drawbacks, they have assumed the controller parameters a fixed value without considering the optimal tuning process. For instance, see [24][25][26][27][28][29].
Lastly, it is worth mentioning that the proposed assessment has provided a robust response for voltage stability through optimal FO-PI controller parameters, enhancing the current status within this field.

Conclusions
In this paper, three evolutionary optimization algorithms have been used to tune the FO-PI controllers implemented in a DAB-based DC microgrid. The aim is to find the optimal parameters for these controllers in terms of voltage stability.
Firstly, the DAB small-signal model has been used for the offline optimization process. Once the three sets of optimal parameters are obtained, a frequency-domain analysis model has been developed to predict the voltage stability of the DC microgrid using the passivity-based criterion.
It has been observed that the microgrid controllers with the obtained three sets of parameters for each algorithm meet the Nyquist criterion, and therefore, stable responses are expected. However, as seen in Section 6, the most stable scenario occurs with the parameters tuned with the GA and the less stable with the PSO algorithm.
To verify such analysis, the detailed model of the DC microgrid has been implemented in PSIM software, where the time-domain simulations are performed. These simulations have been carried out considering several test cases, including changes in the voltage reference and sudden load changes. Hence, as shown in Figure 23 of Section 7, the best dynamic response is achieved by the GA, whereas the least one occurs with the PSO. Thereby, the dependability of the stability analysis provided in Section 7 has been demonstrated. Moreover, it has also been found that the GA provides the lowest voltage ripple during the tested events.
In terms of computation, the GA algorithm is the one that takes more time for each iteration to converge, but the final overall time convergence is roughly 1.5 min, which is between SA and PSO.
By observing the features comparison contained in Table 6, it can be concluded that the GA is the algorithm that offers the best performance for voltage stability as well as an acceptable computation time.
Finally, the influence of the passive parameters of the DAB has also been analyzed. The results summarized in Tables 6-10 indicate that elements such as the DAB leakage inductor or DC link capacitances have no noticeable impact on the DC bus voltage. Therefore, it has been demonstrated that both the microgrid load model and its composition are the dominant parameters in the voltage stability assessment. Funding: This research did not receive any external funding.

Conflicts of Interest:
The authors declare no conflict of interest. Particle number i of the PSO algorithm G Global best position found by all particles L Local (Particle) best position L S Inductance of the AC link inductor R S Equivalent resistance of the AC link inductor v k i Velocity of the particle i at iteration k n,c Neighborhood and current solutions in SA