Continuous-Control-Set Model Predictive Control for Three-Level DC–DC Converter with Unbalanced Loads in Bipolar Electric Vehicle Charging Stations

: Zero-emission transportation is currently a public priority, especially in big cities. For this reason, the use of electric vehicles (EVs) is receiving much attention. To facilitate the adoption of EVs, a proper charging infrastructure together with energy management is essential. This article proposes a design guideline for a direct current (DC) charging station with bipolar properties. A bipolar system can convert a two-wire system into three wires in a microgrid system with a neutral line. The conﬁguration of the bipolar system supports different loads; therefore, the unbalanced operation is inherent to the system. The proposed bipolar DC charging station (CS) has a three-level balancing converter that reduces the step-down effort chargers. Moreover, this paper proposes the continuous-control-set model predictive control (CCS-MPC)-based balancing strategy that allows the handling of different output loads while keeping the neutral-line voltage efﬁciently regulated with improved dynamic performance compared to a traditional controller. Stability and parameter robustness analyses are also performed for the control parameter selection. To ensure the performance of the proposed method, both simulation and experimental results are presented and compared with those obtained from the traditional methods.


Introduction
Due to the rapidly increasing number of electric vehicles (EVs) and increased depletion of fossil fuels, the demand for renewable energy source (RES)-based EV charging stations (CSs) is also rising to alleviate the anxiety of EV drivers [1]. Surveys [2][3][4] depict that the number of available CSs, the charging time, and the range per charge are the biggest concerns of consumers, which also impact the sale of EVs. There is an urgent need to deploy fast CS infrastructure to increase the sale of EVs and reduce the charging time [5]. In addition, the development of RESs such as photovoltaic (PV), hydrogen fuel cell, and DC bus systems have made the distribution system more manageable than conventional distribution systems [6]. Compared to AC distribution systems, RES-based DC power systems due is its fundamental characteristics, such as high bandwidth control, faster transient response, and the ability to account in a simple way for system constraints and nonlinearities [23]. In addition, the other prominent features of MPC are the robustness of the system parameter variations, design simplicity, flexibility to include system constraints, and system non-linearities. Likewise, MPC can also help to achieve multiple control objectives by utilizing cost functions [24]. Therefore, this practice has been frequently applied to control power converters such as NPC [25][26][27] active neutral-clamped converters [28], cascade half-bridge converters [29], modular multi-level converters [30,31], and flying capacitors [32] to balance and control the current and voltage.
From the perspective of power electronic converters, MPC can attain two forms: (1) finite-set model predictive control (FCS-MPC) and (2) continuous-control-set model predictive control (CCS-MPC) [24,33]. The cost function is calculated and minimized for a finite number of switching states in the case of the FCS-MPC strategy. This strategy employs a comparator as an alternative to a modulator to determine the optimal switching frequency, which increases the computational burden and some non-desirable components [24]. However, the CCS-MPC approach can help to address the issues associated with the FCS-MPC approach more conveniently, since this control strategy predicts state variables using a discrete state model of the power converters where the state-space variables can be estimated using a conventional state observer [23]. The cost function is calculated and minimized inside the prediction horizon for obtaining vectors of future control actions. According to the receding horizon control, only the first sample of the vector is used as a space vector modulator, which leads to a fixed switching frequency; therefore, it offers less computational burden. One can always refer to a detailed comparative analysis conducted in [34,35] to compare the dynamic performance and make an appropriate choice of a power converter for a given application. Owing to the increasing use of CCS-MPC in industrial applications [35], the decoupling of switching frequency from the controller sampling time helps to operate the converter at the fixed switching frequency and enhances the control flexibility. With these benefits, it is necessary to emphasize that the FCS-MPC has been studied for a bipolar DC microgrid under a variable switching frequency, which may lead to increased computational cost [36,37]. Nevertheless, the implementation of CCS-MPC has not been very well explored for bipolar DC microgrid applications. Therefore, this paper proposes an enhanced CCS-MPC-based control strategy with a half-bridge balancing converter. This method has a high utilization of DC-DC converters to reduce the energy loss and improve the utilization of the DC power supply. Consequently, the bipolar DC CS has unique benefits in cost, power supply flexibility, and security. In a bipolar DC CS, there is coupling between the positive and the negative poles. The imbalance of the output voltages from positive and negative power supplies, loads, and line resistance can cause an unbalanced current in the neutral line, resulting in voltage deviations on the bipolar DC bus from the rated value. The proposed method has demonstrated its effectiveness in controlling power converters with a faster transient response and achieves an excellent dynamic performance without the need for any stability analysis. Moreover, the proposed CCS-MPC reduces the steady-state oscillation and increases the robustness against parameter uncertainties compared to previous methods [13][14][15][16]. In addition, the involvement of a simple cost function approach helps to avoid the tedious controller design methods such as Bode plot and Nyquist analysis approaches. CCS-MPC techniques have been widely used in different systems, but their application in bipolar DC microgrids still requires further effort and attention. The proposed control strategy opens up a new research topic for implementing the CCS-MPC in bipolar DC microgrids, as the bipolar DC bus suffers from intermittent operations due to the unbalanced power distribution between the DC buses. The three-level balancing converter combined with CCS-MPC provides a higher degree of freedom to accurately control the bipolar system during sudden load changes to achieve a better voltage balance. These asymmetrical operations can cause an unbalance in the three-phase voltage that deteriorates the power quality and leads to a premature failure The key contribution of this article is the design of a CCS-MPC controller for a bipolar EV charging station with the following properties:

•
Robustness of the output voltage over sudden load changes. • Robustness against model parameter uncertainties. • A fixed switching frequency is achievable with a reduced computational burden.

•
The implementation of a CCS-MPC control scheme maintains the bipolar DC bus with better damping performance than an optimal PI controller.

•
The steady-state oscillations are reduced compared to an optimal PI controller.
The organization of this paper is as follows: the generic topological description of the proposed DC CS is introduced in Section 2. The CCS-MPC-based control scheme with stability analysis applied to the half-bridge balancing converter is described in Section 3. Subsequently, the operating modes of the balancing converter are analyzed in Section 4. The numerical results, control parameter robustness analysis, and comparison of the proposed controller with the optimal PI controller based on the genetic algorithm (GA) are carried out in Section 5. Finally, the valuable interpretations of the proposed work are summarized in Section 6.

Topology of Bipolar DC EV Charging Station
The topology of the proposed bipolar CS architecture is exemplified in Figure 1, where a three-level balancing DC-DC converter (see Figure 1b) is employed and controlled using the CCS-MPC strategy to make this system bipolar to achieve better DC bus voltage/current regulation. The bipolar structure allows the connection of distributed power sources such as RESs (PV, wind, and fuel cell) and energy storage systems [9,38]. Conventional DC-DC converters are selected for interfacing individual power sources with the bipolar DC buses according to the capacity of the RESs. Similarly, generic bi-directional DC-DC converters [13] are employed to interface battery energy sources with bipolar DC buses, as batteries require charging and discharging. The fuel cells and PV arrays only supply power to the bipolar DC bus, so conventional boost converters are popular for interfacing them with the bipolar DC buses. Note that the detailed study of the conventional converters is exempted from this investigation, and only the three-level balancing converter is part of the investigation. CCS-MPC provides a higher degree of freedom to accurately control the bipolar system during sudden load changes to achieve a better voltage balance. These asymmetrical operations can cause an unbalance in the three-phase voltage that deteriorates the power quality and leads to a premature failure of semiconductor devices. The application of the CCS-MPC approach will allow excellent output voltage regulations during load fluctuations.
The key contribution of this article is the design of a CCS-MPC controller for a bipolar EV charging station with the following properties: • Robustness of the output voltage over sudden load changes. • Robustness against model parameter uncertainties.

•
A fixed switching frequency is achievable with a reduced computational burden.

•
The implementation of a CCS-MPC control scheme maintains the bipolar DC bus with better damping performance than an optimal PI controller.

•
The steady-state oscillations are reduced compared to an optimal PI controller.
The organization of this paper is as follows: the generic topological description of the proposed DC CS is introduced in Section 2. The CCS-MPC-based control scheme with stability analysis applied to the half-bridge balancing converter is described in Section 3. Subsequently, the operating modes of the balancing converter are analyzed in Section 4. The numerical results, control parameter robustness analysis, and comparison of the proposed controller with the optimal PI controller based on the genetic algorithm (GA) are carried out in Section 5. Finally, the valuable interpretations of the proposed work are summarized in Section 6.

Topology of Bipolar DC EV Charging Station
The topology of the proposed bipolar CS architecture is exemplified in Figure 1, where a three-level balancing DC-DC converter (see Figure 1b) is employed and controlled using the CCS-MPC strategy to make this system bipolar to achieve better DC bus voltage/current regulation. The bipolar structure allows the connection of distributed power sources such as RESs (PV, wind, and fuel cell) and energy storage systems [9,38]. Conventional DC-DC converters are selected for interfacing individual power sources with the bipolar DC buses according to the capacity of the RESs. Similarly, generic bidirectional DC-DC converters [13] are employed to interface battery energy sources with bipolar DC buses, as batteries require charging and discharging. The fuel cells and PV arrays only supply power to the bipolar DC bus, so conventional boost converters are popular for interfacing them with the bipolar DC buses. Note that the detailed study of the conventional converters is exempted from this investigation, and only the three-level balancing converter is part of the investigation.  This system can provide DC power to several loads. These loads are expected to be either conventional or fast chargers for electrified transportation with separate DC-DC converters. As the configuration of the bipolar system supports distinct loads, the unbalanced operation is, therefore, inherent to the system. The unbalanced behavior can be mitigated by supplying a part of the power to the other loads on the DC bus through an inductor. The details will be studied in the subsequent sections.

Proposed CCS-MPC Strategy for Three-Level Balancing Converter
A half-bridge balancing DC-DC converter is shown in Figure 1b, in which the inductor L is connected between the series of switches S 2 , S 3 and the neutral points of capacitors C 1 , C 2 . With the help of inductor (L), the potential shoot-through problem is avoided by transferring the unbalance power to the other load. An enhanced MPC strategy of the duty cycle is presented in Figure 2. To achieve a widespread conducive time for the freewheeling current time, an external DC bias voltage (U Bias ) is introduced to maintain the three-level states. The switches S 1 to S 4 are driven by complementary driving signals u gs1 , u gs2 , u gs3 and u gs4 with a dead time (T d ), where the freewheeling time (T f ), sampling period (T s ) and U pp is the ripple value of the carrier wave. The time T f can be calculated using the following equation: where T f can be regulated by changing the U bias voltage. The reference voltage u re f and the output voltage u c2 are equal to u in /2, as shown in Figure 2.
Mathematics 2022, 10, x FOR PEER REVIEW 5 of 24 This system can provide DC power to several loads. These loads are expected to be either conventional or fast chargers for electrified transportation with separate DC-DC converters. As the configuration of the bipolar system supports distinct loads, the unbalanced operation is, therefore, inherent to the system. The unbalanced behavior can be mitigated by supplying a part of the power to the other loads on the DC bus through an inductor. The details will be studied in the subsequent sections.

Proposed CCS-MPC Strategy for Three-Level Balancing Converter
A half-bridge balancing DC-DC converter is shown in Figure 1b, in which the inductor L is connected between the series of switches , and the neutral points of capacitors , . With the help of inductor ( ), the potential shoot-through problem is avoided by transferring the unbalance power to the other load. An enhanced MPC strategy of the duty cycle is presented in Figure 2. To achieve a widespread conducive time for the freewheeling current time, an external DC bias voltage ( ) is introduced to maintain the three-level states. The switches to are driven by complementary driving signals , , and with a dead time ( ), where the freewheeling time ( ), sampling period ( ) and is the ripple value of the carrier wave. The time can be calculated using the following equation: where can be regulated by changing the voltage. The reference voltage ( ) and the output voltage are equal to /2, as shown in Figure 2.

System Modeling
The electrical circuit diagram of the three-level balancing DC-DC converter is visualized in Figure 1b. In case and are turned on, and and are turned off, the differential equations of the large signal model [16] can be written as follows:

System Modeling
The electrical circuit diagram of the three-level balancing DC-DC converter is visualized in Figure 1b. In case S 1 and S 2 are turned on, and S 3 and S 4 are turned off, the differential equations of the large signal model [16] can be written as follows: In case S 1 and S 2 are turned off, and S 3 and S 4 are turned on, the differential equations become: Mathematics 2022, 10, 3444 6 of 23 According to the average small signal, model Equations (2) and (3) can be expressed as: The duty cycles of S 1 and S 4 can be expressed as d = D +d = t on /T s , thus, the duty cycles of S 3 and S 4 will be defined as (1 − d), where D andd are the steady-state duty ratio and the small perturbation of the duty cycle, respectively [14]. Further, it is worth noting that d is the control variable in the proposed system, andû in and U in are the small signal input voltage and the steady-state input voltage value.

Discrete Model with an Embedded Integrator
The augmented model of the converter is derived for model predictive control. At first, the model described using the following set of equations, i.e., (4) and (5), must be discretized as follows: where the sampling period is denoted by T s . The following reduced discrete state-space equations can help to obtain the linear model of the circuit: where x m = i L u c2 T is the state-space vector, u =d is the control variable, e = u in is the disturbance vector, and y m = u c2 . The matrices of disturbance vectors A m , B m , C m , E m are defined as follows: Now, an integrator is embedded in the model to eliminate the steady-state error. Applying the difference operation on both sides of (8), we get: where ∆x m , ∆e, and ∆u are the incremental state variables, disturbance variable, and control variable, respectively. Applying the same procedure for the incremental output vector yields: Note that (12) indicates the presence of an integrator. So, a new space vector containing the vectors ∆x m (k), and y(k) is expressed as follows: Finally, we obtained the augmented state-space model by combining (11) and (12), stated as follows: where the matrices of this augmented model are A m , B m , C m and E m , as follows: where I and 0 are the identity and zero square matrices, and p and q are the number of state variables and output variables, respectively.

Model Predictive Control
This approach involves the design of the embedded integrating model for the CCS-MPC controller. For the control signal, a predictive horizon of the current and output voltages of the converter is predicted inside the prediction horizon (N p ). The control block diagram of the MPC algorithm is shown in Figure 3. The inputs of this block are the input voltage (u in ), the output voltage (u c2 ), and the reference voltage (u re f ). The incremental model (12) generates the estimated vector x(k). The optimal control signal ∆u is achieved using the cost functions (28)- (31). An integrator is embedded in the model to generate the control signal u. Considering the vectors of the future incremental control signals at k j is described below: where N c is the control horizon, and the prediction of future state variables can be done with the help of given information x(k j ) for the prediction horizon, N p , with the main objective to find the control signals of the vector dimensions of N c . Assuming an incremental statespace model described by the matrices A, B, and C, and the disturbance vector ∆e(k) in (14), the future state variables can be evaluated in sequence from (16): where x k j + N p k j is the future state predicted variables at k j + N p from the control system x k j . Similarly, the predicted outputs variable can be achieved from the predicted state variables, as follows: Mathematics 2022, 10, x FOR PEER REVIEW 7 of 24 Note that (12) indicates the presence of an integrator. So, a new space vector containing the vectors ∆ ( ), and ( ) is expressed as follows: Finally, we obtained the augmented state-space model by combining (11) and (12), stated as follows: where the matrices of this augmented model are , , and , as follows: where I and 0 are the identity and zero square matrices, and p and q are the number of state variables and output variables, respectively.

Model Predictive Control
This approach involves the design of the embedded integrating model for the CCS-MPC controller. For the control signal, a predictive horizon of the current and output voltages of the converter is predicted inside the prediction horizon ( ). The control block diagram of the MPC algorithm is shown in Figure 3. The inputs of this block are the input voltage ( ), the output voltage ( ), and the reference voltage ( ). The incremental model (12) generates the estimated vector ( ) . The optimal control signal ∆ is achieved using the cost functions (28)- (31). An integrator is embedded in the model to generate the control signal u. Considering the vectors of the future incremental control signals at is described below: Cost Function Minimization (28)-(31) where is the control horizon, and the prediction of future state variables can be done with the help of given information x( ) for the prediction horizon, , with the main objective to find the control signals of the vector dimensions of . Assuming an incremental state-space model described by the matrices , , and , and the disturbance vector ∆ ( ) in (14), the future state variables can be evaluated in sequence from (16): It is clearly seen that the predicted future state variables and future control trajectory are stated in terms of x k j and ∆u, respectively. Thus, the output vector Y can be described as: Mathematics 2022, 10, 3444 8 of 23 From (18), a complex matrix form is obtained, expressed as: where F is a vector of dimension N p , and G, and H are matrices with dimension of N c × N P .
where q is the number of output variables. In our case q = 1. Now, let us define the following cost function: where R = r w I p N c pN c is a diagonal matrix and r w > 0 is the control effort used as a tuning parameter to achieve desired closed-loop performance. Moreover, Y * is the reference variable. Note that the reference is maintained practically constant inside the predictive windows; hence, Y * is expressed as: where r k j = u in 2 to achieve the voltage balance between the two DC output voltages. Then, doing the first derivative of J with respect to ∆u yields: where M = G T G + R. Then, by equalizing (25) to zero and considering the optimal vector in (26), the incremental control signal can be expressed as follows: It is noteworthy that the ∆u shows all incremental values for sampling instant k j to (k j + N c − 1). As (13)- (14) are time variants, the incremental control can be stated as a state feedback controller. In line with (15), the actual incremental can be defined as follows: (28) where: Mathematics 2022, 10, 3444 9 of 23 Note that K c = K x K r is the MPC controller gain. Furthermore, according to the principle of receding horizon, W is the N c dimensional matrix I p×q 0 p×q . . . 0 p×q . Finally, the optimal control signal can be achieved by adding the incremental value ∆u k j in the control signal of the previous instant, as follows: or equivalently: The increment of the control signal contains the voltage error, thus performing the integration in discrete form for the control signal that can achieve a zero steady-state error. Bearing this in mind, the duty cycle can be calculated as follows: where K r , K x , and K h are the gain obtained by minimizing the cost function (24). Thereby, the converter operations can be properly carried out once the switching periods are calculated, as will be discussed in Section 4.

Design Guidelines
The control design parameters of the CCS-MPC controller are selected according to the dynamic specifications of the controller. The value of the control effort parameter r w is set based on the compromise between the maximum deviation during the transient response and the settling time, according to the following specifications: (1) To achieve the stability of the prediction, horizon N p and control horizon N c are designed large enough [22]. Therefore, the stability analysis of the MPC controller is used to achieve the design specifications. With this in mind, the stability analysis of the proposed system has been developed using the control-loop poles from the small signal model described by (14) and (15). Hence, the closed-loop poles can be calculated as the function of the control effort parameter r w , the prediction and control horizons, N p and N c , and the electrical parameters of the power converter. Figure 4 shows the position of the closed-loop poles for different prediction and control horizons and sweeping the control effort, r w , from 0 to 10, where blue crosses (x) shows poles with different r w values. In Figure 4a, the values of prediction horizon N p = 5 and control horizon N c = 3 are selected, respectively. For these values of N p and N c , when the value of the control effort r w increases from 0.01 to 10, the system tends to be unstable, as shown in Figure 4a. In contrast, Figure 4b shows the poles with prediction horizon, N p = 150 and control horizon N c = 50, the poles move closer to the unitary circle. The system looks stable with the given values of control effort r w . It should be noted that the stability is ensured for the whole range of the control effort r w . Therefore, N p = 150 and N c = 50 are selected for the CCS-MPC controller's implementation. Finally, based on the given design specifications, the control effort is selected as r w = 0.1 to achieve the required transient performance. In the Figure 4b the red crosses (x) represents the closed loop poles for the designed parameters.

Operating Principle of Balancing Converter
Some simplifying assumptions are made. Firstly, the two capacitors and are equal, and the output voltage does not change during the switching period. The freewheeling diodes and switches are ideal. The turn on and off periods of the switches are not considered. As depicted in Figure 1b, when the load current of ( ) is greater than the load current of ( ), the input voltage will transfer a part of the unbalanced power to through inductor to balance the output voltages and . The inductor current is ( > 0), and during the freewheeling time, when there is no power transfer to the load , the is greater than the difference of and . In contrast, the inductor current is ( < 0).

Operation Principle under State I
The switching time diagrams of the semiconductor power switches to , the voltage stress on the power switches, and the inductor current are shown in Figure 5, together with the drain-source voltages to and the across power switches to . Altogether, there are eight modes of operation, when ( > 0) as shown in the equivalent circuits of the balancing converter in Figure 6. (2) Mode-II ( , ) As shown in Figures 5 and 6b, during the time interval, to the switch is turned on, and runs through the and . Since switches and are turned on, the voltage stress across and is equal to 0.5. Meanwhile, the voltage stress across is equal to zero.
(3) Mode-III ( , ) As shown in Figures 5 and 6c, in this mode, the switch is turned off and the will pass through the freewheeling diodes and from and .
(4) Mode-IV ( , ) As shown in Figures 5 and 6d, in this mode, the switch is turned on at the time interval , and passes through and , decreases to = − , and there is similar voltage stress on the power switches as in the third mode.

Operating Principle of Balancing Converter
Some simplifying assumptions are made. Firstly, the two capacitors C 1 and C 2 are equal, and the output voltage does not change during the switching period. The freewheeling diodes and switches are ideal. The turn on and off periods of the switches are not considered. As depicted in Figure 1b, when the load current of (I R2 ) is greater than the load current of (I R1 ), the input voltage u in will transfer a part of the unbalanced power to R 2 through inductor L to balance the output voltages u c1 and u c2 . The inductor current is (I L > 0), and during the freewheeling time, when there is no power transfer to the load R 2 , the I L is greater than the difference of I R2 and I R1 . In contrast, the inductor current is (I L < 0).

Operation Principle under State I
The switching time diagrams of the semiconductor power switches S 1 to S 4 , the voltage stress on the power switches, and the inductor current I L are shown in Figure 5, together with the drain-source voltages u s1 to u s4 and the across power switches S 1 to S 4 . Altogether, there are eight modes of operation, when (I L > 0) as shown in the equivalent circuits of the balancing converter in Figure 6.  As shown in Figures 5 and 6e, in this interval, the switch is turned on, and runs through and from .
will be increased to = . The voltage stress on and is 0.5 . From on, the next interval will start.

Operation Principle under State II
The switching time diagram is depicted in Figure 7. When is working under state II, it is seen that the first five modes during the time intervals ( , ) and ( , ) are the same as in Figure 5. Therefore, only the modes from to are discussed here and their equivalent circuit is shown in Figure 8.
(1) Mode-V ( , ) As shown in Figures 7 and 8a, in this time interval, the switches and are turned on, is decreased to zero, and = . The voltage stress on switches and will be , and the voltage stress on and is zero and , respectively.
(2) Mode-VI ( , ) As shown in Figures 7 and 8b, in this mode, is turned off and passes through and . As = 0, does not change, and the voltage stress on is half of the input voltage. The voltage stress on and is /4. Therefore, the voltage stress on freewheeling diodes and are /4 and zero, respectively.
(3) Mode-VII ( , ) (1) Mode-I (t 0 , t 1 ) As shown in Figures 5 and 6a, in this interval, the switch S 1 is turned off and I L runs through D 5 . The value of I L will not change, as the inductor voltage u L is zero, and the voltage of u s1 is equal to u c1 , and the sum of us s3 + u s4 = u c2 . Under the balance states u c1 = u c2 = 0.5 u in . Similarly, the voltage stresses across u s1 , u s3 , u s4 are 0.5u_in, 0.25u in , and 0.25u in , respectively, and the voltage stress on u D6 is u c2 − u s4 = 0.25u in .
(2) Mode-II (t 1 , t 2 ) As shown in Figures 5 and 6b, during the time interval, t 1 to t 2 the switch S 3 is turned on, and I L runs through the D 5 and S 2 . Since switches S 2 and S 3 are turned on, the voltage stress across u S1 and u s4 is equal to 0.5. Meanwhile, the voltage stress across u D6 is equal to zero.
(3) Mode-III (t 2 , t 3 ) As shown in Figures 5 and 6c, in this mode, the switch S 2 is turned off and the I L will pass through the freewheeling diodes D 3 and D 4 from S 2 and D 5 .

(4) Mode-IV (t 3 , t 4 )
As shown in Figures 5 and 6d, in this mode, the switch S 4 is turned on at the time interval t 3 , and I L passes through D 1 and D 3 , u L decreases to u L = −u C2 , and there is similar voltage stress on the power switches as in the third mode.
As shown in Figures 5 and 6c, in this mode, the switch S 4 is turned off, and the other conditions are the same as in the fourth mode.  As shown in Figures 5 and 6a, in this time interval, the switch S 3 is turned off, and I L runs through D 5 and S 2 . In addition, the value of i L does not change as u L = 0. The voltage stress on S 1 , S 3 and S 4 is 0.5u in and u D6 is zero.
As shown in Figures 5 and 6e, in this interval, the switch S 1 is turned on, and I L runs through S 1 and S 2 from D 5 . I L will be increased to u L = u C1 . The voltage stress on S 3 and S 4 is 0.5u in . From t 8 on, the next interval will start.

Operation Principle under State II
The switching time diagram is depicted in Figure 7. When I L is working under state II, it is seen that the first five modes during the time intervals (t 0 , t 4 ) and (t 9 , t 10 ) are the same as in Figure 5. Therefore, only the modes from t 4 to t 9 are discussed here and their equivalent circuit is shown in Figure 8. As shown in Figures 7 and 8e, in this mode, S1 is turned on and IL runs through D1 and D2. The value of IL decreases to negative uc1. During this mode, the voltage stress on the power switches is similar to mode-VIII. From time t10, the next working period will start, respectively.

Results and Discussion
The dynamic performance evaluation of the proposed improved topological DC bus system, simulated in the MATLAB/Simulink environment with a half-bridge three-level balancing DC-DC converter is considered. In addition, an 800 V charging station with a CCS-MPC control algorithm is developed to validate the proposed architecture. The parameters of the simulated system are listed in Table 1.  (1) Mode-V (t 4 , t 5 ) As shown in Figures 7 and 8a, in this time interval, the switches S 3 and S 4 are turned on, I L is decreased to zero, and u L = u c2 . The voltage stress on switches S 1 and S 2 will be u in 2 , and the voltage stress on u D5 and u D6 is zero and u in 2 , respectively. (2) Mode-VI (t 5 , t 6 ) As shown in Figures 7 and 8b, in this mode, S 4 is turned off and I L passes through S 3 and D 6 . As u L = 0, I L does not change, and the voltage stress on S 4 is half of the input voltage. The voltage stress on S 1 and S 2 is u in /4. Therefore, the voltage stress on freewheeling diodes u D5 and u D6 are u in /4 and zero, respectively.
(3) Mode-VII (t 6 , t 7 ) As shown in Figures 7 and 8c, in this interval, S 2 is turned on, and I L runs through S 3 and D 6 . The value of I L does not change as u L = 0. The voltage stress on S 1 and S_4 is u in /2 and on u D5 and u D6 is zero.
(4) Mode-VIII (t 7 , t 8 ) As shown in Figures 7 and 8d, in this time interval, S 3 is turned off, I L runs through D 1 and D 2 , and its value decreases positively. The voltage stress on switches S 3 and S 4 is u in /2, and the voltages stress on u D5 and u D6 is u in /2 and zero, respectively.
As shown in Figures 7 and 8e, in this mode, S 1 is turned on and I L runs through D 1 and D 2 . The value of I L decreases to negative u c1 . During this mode, the voltage stress on the power switches is similar to mode-VIII. From time t 10 , the next working period will start, respectively.

Results and Discussion
The dynamic performance evaluation of the proposed improved topological DC bus system, simulated in the MATLAB/Simulink environment with a half-bridge three-level balancing DC-DC converter is considered. In addition, an 800 V charging station with a CCS-MPC control algorithm is developed to validate the proposed architecture. The parameters of the simulated system are listed in Table 1. Table 1. System parameters.

Description Symbol Value
Capacitors   Figure 9b shows the drain-source voltages, inductor current, and output voltages when I R1 = 2.5 A and I R2 = 5 A, and u C1 and u C2 are 400.2 V and 400.1 V, respectively. In this state, I L is always greater than zero. Figure 9c,d shows the stable numerical waveforms under state II, when I L is either greater or less than zero. Figure 9c shows the switching drain-source voltages, inductor current, and output voltages when the load current I R2 is less than I R1 , where I R1 = 7.5 A, I R2 = 1.3 A, and the output voltages u c1 and u c2 are 399.82 V and 400.23 V, respectively. Likewise, Figure 9d shows similar results when I R1 = 5 A, I R2 = 3 A, and u c1 and u c2 are 400 V and 400.2 V, respectively. According to Figure 9, the simulation results show that the voltage stress on switches is half of the input voltage, which validates the theoretical analysis shown in Figures 5 and 7. The reason behind the voltage spikes of the drain-source voltages is due to the resonance when the switches are turned off. It is seen that the output voltages u c1 and u c2 are very close to 400 V when using the MPC controller in both states. Figure 9 illustrates the stable waveform under different states, where Figure 9a depicts the switches drain-source voltages, inductor current, and the output voltages. When the load current ( > ), where = 1.2 A, = 8.5 A, the output voltages and are 400.2 V and 399.8 V, respectively. Similarly, Figure 9b shows the drain-source voltages, inductor current, and output voltages when = 2.5 A and = 5 A, and and are 400.2 V and 400.1 V, respectively. In this state, is always greater than zero. Figure 9c,d shows the stable numerical waveforms under state II, when IL is either greater or less than zero. Figure 9c shows the switching drain-source voltages, inductor current, and output voltages when the load current IR2 is less than IR1, where IR1 = 7.5 A, IR2

Performance Comparison of Transient Load Response
This subsection presents a comparison between the proposed control and previous state-of-the-art controllers. For this purpose, the proposed CCS-MPC will be compared with the PI controller presented in [16]. In order to conduct a fair comparison, the PI controller was designed using a genetic algorithm to optimize the control gains, hence improving the transient response due external disturbances such as sudden load changes.
The optimal GA-PI control gains are shown in Table 1 as K p and K i for the proportional gain and integral gain, respectively.
In order to validate the dynamic performance of the proposed control scheme, one simulation test was performed. For this purpose, the controllers were tested under a sudden load change. This simulation test has the following pattern. From t = 0 s to t = 0.3 s, I R1 = 0 A and I R2 = 1 A. Then, at t = 0.3 s, I R1 increases to a value equal to 5 A, while I R2 remains constant. Finally, at t = 0.7 s, I R1 decreases to 0 A and I R2 remains constant. Figures 10 and 11 show the numerical results under transient load changes, using the proposed CCS-MPC and the GA-PI, respectively. In both cases, the u C1 and u C2 are almost equal, although they have apparent fluctuations during load transients. Moreover, the proposed CCS-MPC-based control strategy tracks the reference voltages with negligible fluctuations, while there are apparently large fluctuations in the case of the GA-PI controller due to its faster and smoother transient response with lesser control complexity than the conventional linear controller. In addition, the MPC has a lower rising time, settling time, and overshoots than the GA-PI controller. Moreover, the CCS-MPC controller reduces the dangerous oscillations in the inductor and provides a smoother operation under load transients. At the same time, the conventional PI controller has more fluctuations in the output voltages with considerable ripples before and after the load transients. The proposed CCS-MPC considerably reduces the current over the inductor, thus allowing a reduction in the filter components.
proposed CCS-MPC-based control strategy tracks the reference voltages with negligi fluctuations, while there are apparently large fluctuations in the case of the GA-PI c troller due to its faster and smoother transient response with lesser control complex than the conventional linear controller. In addition, the MPC has a lower rising time, s tling time, and overshoots than the GA-PI controller. Moreover, the CCS-MPC contro reduces the dangerous oscillations in the inductor and provides a smoother operation der load transients. At the same time, the conventional PI controller has more fluctuatio in the output voltages with considerable ripples before and after the load transients. T proposed CCS-MPC considerably reduces the current over the inductor, thus allowin reduction in the filter components.

Parameter Robustness Analysis
To confirm the parameter robustness of the proposed CCS-MPC controller against parameter uncertainty, the simulation results were compared with the designed GA-PI controller by changing the inductance and capacitance of the balancing converter.
Firstly, the proposed CCS-MPC and conventional GA-PI controller with a nominal inductance of 5 mH were tested at the 50% inductance nominal value. Figure 12 shows the output voltages and inductor current of the proposed CCS-MPC scheme with an inductance of 2.5 mH, and Figure 13 depicts the output voltages and inductor current of the GA-PI controller with a 2.5 mH inductance value. It is apparent that the conventional PI controller has an oscillatory behavior, while the CCS-MPC is more robust in case of an inductor parameter mismatch.
Secondly, to verify the robustness of the proposed CCS-MPC against the capacitance mismatch, the capacitance of the balancing converter was set to 235 μF, which is 50% of

Parameter Robustness Analysis
To confirm the parameter robustness of the proposed CCS-MPC controller against parameter uncertainty, the simulation results were compared with the designed GA-PI controller by changing the inductance and capacitance of the balancing converter.
Firstly, the proposed CCS-MPC and conventional GA-PI controller with a nominal inductance of 5 mH were tested at the 50% inductance nominal value. Figure 12 shows the output voltages and inductor current of the proposed CCS-MPC scheme with an inductance of 2.5 mH, and Figure 13 depicts the output voltages and inductor current of the GA-PI controller with a 2.5 mH inductance value. It is apparent that the conventional PI controller has an oscillatory behavior, while the CCS-MPC is more robust in case of an inductor parameter mismatch.
Secondly, to verify the robustness of the proposed CCS-MPC against the capa mismatch, the capacitance of the balancing converter was set to 235 μF, which is the nominal value. Figure 14 shows the dynamic output voltages and the inductor of the CCS-MPC control schemes, while Figure 15 shows the dynamics of the GAtroller. The proposed CCS-MPC shows robust control of the output voltages as w the inductor current. In contrast, the conventional GA-PI controller has apparent voltage deviations and inductor current ripples under the capacitance mismatch.   Secondly, to verify the robustness of the proposed CCS-MPC against the capacitance mismatch, the capacitance of the balancing converter was set to 235 µF, which is 50% of the nominal value. Figure 14 shows the dynamic output voltages and the inductor current of the CCS-MPC control schemes, while Figure 15 shows the dynamics of the GA-PI controller. The proposed CCS-MPC shows robust control of the output voltages as well as in the inductor current. In contrast, the conventional GA-PI controller has apparent output voltage deviations and inductor current ripples under the capacitance mismatch.

Experimental Results
To validate the system performance, the experimental results of the proposed control scheme were carried out using the dSPACE 1006 platform for real-time control, as depicted in Figure 16.

Experimental Results
To validate the system performance, the experimental results of the proposed control scheme were carried out using the dSPACE 1006 platform for real-time control, as depicted in Figure 16.  Figure 17 shows the transient load response of the MPC when the load is increasing, while Figure 18 shows the response of the proposed scheme when is increasing. Figures 19 and 20 show the HIL results under the transient load with the CCS-MPC and PI controller. Figure 19a shows the output voltage response with the MPC, while Fig-Figure 16. Hardware-in-loop experimental setup. Figures 17 and 18 show the hardware-in-loop (HIL) results under the transient load changes. Figure 17 shows the transient load response of the MPC when the load I R1 is increasing, while Figure 18 shows the response of the proposed scheme when I R2 is increasing.

Comparative Analysis
In Table 2, the performance comparison with regard to the proposed MPC is compared with the PI controller. Note that there is better accuracy in tracking with the MPC controller. The DC link voltages, i.e., U C1 and U C2 , are observable in the case of the MPC controller. Moreover, the proposed MPC approach also provides less peak-to-peak deviation in the voltage ripples when different test cases are performed with respect to load current variations I R1 and I R2 . As a result, less heating in the loads is evident, since the peak-topeak load ripple is lesser in magnitude. Nevertheless, the computational complexity of the proposed MPC is slightly higher than for GA-PI. However, the proper tuning of a PI controller is a tedious task to achieve for gaining high robustness. Hence, the proposed MPC is a suitable choice for meeting the objectives of better control with a three-level DC-DC converter compared to PI controllers.

Conclusions
In this paper, an enhanced CCS-MPC-based bipolar EV charging station is proposed under a MATLAB/Simulink environment, and the HIL experimental results are verified. A three-level balancing converter with an improved CCS-MPC control strategy is applied to make this system bipolar and balance the DC bus. The operation principle of the balancing converter is analyzed, and its stability analysis is developed for selecting the control parameters of the proposed control strategy. Our results indicate that this topology helps to solve the shoot-through problem by building a neutral line to balance the DC buses for different loads in EV charging stations. The proposed control strategy, with a three-level balancing converter, has lower voltage and current stresses on semi-conductor devices and can balance the DC bus voltages with extremely unbalanced load conditions. The numerical results corroborate the fact that the proposed CCS-MPC controller can better balance the bipolar DC CS under load fluctuations than the optimal GA-PI controller. In the bipolar DC CS, the proposal can complete the control with a reduced computational burden, which greatly saves the cost. However, more research is still needed in order to further optimize and validate the proposed design and thus, this paper opens a new research topic potentially leading to the implementation of the MPC in bipolar DC microgrids.
As the use of EVs is becoming more popular, the modern infrastructure of EVs is also developing toward the zero-emission goal. Therefore, the MPC approaches can be used in the design and development of bipolar-based DC EV charging stations. As an open topic for future research, the application of reinforcement learning algorithms and neuro-adaptive approaches in the design of the MPC can be studied.