# Energy management of electric vehicle using a new strategy based on slap swarm optimization and differential flatness control

Houssam Eddine Ghadbane, Said Barkat, Ali Djerioui, Azeddine Houari, Mihai Oproescu, Nicu Bizon

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Houssam Eddine Ghadbane, Said Barkat, Ali Djerioui, Azeddine Houari, Mihai Oproescu, et al. Energy management of electric vehicle using a new strategy based on slap swarm optimization and differential flatness control. Scientific Reports, 2024, 14 (1), pp.3629. 10.1038/s41598-024-53396-3. hal-04531516

HAL Id: hal-04531516

https://hal.science/hal-04531516

Submitted on 3 Jun 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Energy management of electric vehicle using a new strategy based on slap swarm optimization and differential flatness control

Houssam Eddine Ghadbane, Said Barkat, Ali Djerioui, Azeddine Houari, Mihai Oproescu, Nicu Bizon

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Houssam Eddine Ghadbane, Said Barkat, Ali Djerioui, Azeddine Houari, Mihai Oproescu, et al. Energy management of electric vehicle using a new strategy based on slap swarm optimization and differential flatness control. Scientific Reports, 2024, 14 (1), pp.3629. 10.1038/s41598-024-53396-3. hal-04531516

HAL Id: hal-04531516

https://hal.science/hal-04531516v1

Submitted on 3 Jun 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

OPEN

# Energy management of electric vehicle using a new strategy based on slap swarm optimization and differential flatness control

Houssam Eddine Ghadbane<fnref n="1" />,*, Said Barkat²,*, Ali Djerioui², Azeddine Houari³, Mihai Oproescu⁴ & Nicu Bizon⁴,⁵,⁶

Optimal energy management of electric vehicles using slap swarm optimization and differential flatness control has been proposed. A battery-supercapacitor power system is adopted. Each source is connected in parallel to the DC-bus using DC-DC bidirectional converters and supplies a synchronous reluctance motor (SynRM) based drive. The proposed EMS fundamental forces lie in using a combination of complementary proprieties of two approaches, a Slap Swarm optimization Algorithm and Differential Flatness (DF). With a fast optimization mechanism, the Slap Swarm optimization algorithm allows adapting in real-time conditions the DF gains to optimize the system performances. On its side, DF uses predefined trajectories respecting the physical proprieties of the system, which is a powerful tool to guarantee the dynamic constraints of the sources when ensuring desired robust control proprieties. To check the feasibility and performance of the suggested EMS, comprehensive processor-in-the-loop co-simulations of the electric vehicle were carried out using the C2000 launchxlf28379d DSP board. The main goal of the proposed EMS is to guarantee the DC-bus stabilization, reducing the DC-bus voltage ripples (Δv = 5 V) and the voltage overshoots 15 V (3.2%), respect the source dynamics, and satisfy the SynRM motor power demand. Furthermore, the algorithm minimizes induced harmonics by the drive (10.49%), reducing the battery current ripple by 17.15A, thereby enhancing the battery lifecycle.

Nowadays, the influence of the transportation system's growth on global warming and climate change has become more and more tangible and has impacted our everyday lives<fnref n="1" />. According to², integrating the road transport sector into the EU ETS (European Emissions Trading System) is expected to be a cost-effective CO₂ reduction method. The transport sector accounted for roughly 27% of total EU greenhouse gas (GHG) emissions in 2017, with light commercial vehicles (LCVs) contributing to approximately 9% of EU transport GHG emissions³. According to the European Environment Agency (EEA), average CO₂ emissions from new light-duty vehicles (LDVs) rose in 2017 for the first time since 2010 and continued to rise in 2018⁴. The European Commission published a proposal to outlaw the sale of petrol and diesel automobiles in the EU from 2035⁵. In this context, transportation decarbonization is primarily recognized as one of the priorities for considerably decreasing harmful emissions.

The facts presented here confirm that the automobile industry has adopted the electrified mobility vision. Several automakers have openly stated their commitment to an electric vehicle future, expecting the EV market to be poised for rapid expansion. For example, Volvo Company stated that cars with internal combustion engines have no long-term future. Volvo is firmly dedicated to becoming an all-electric vehicle manufacturer, with the transition expected to occur by 2030⁶. By 2035, General Motors intends to eliminate tailpipe emissions from new

<footnote n="1">¹Electrical Engineering Laboratory of Guelma (LGEG), Electrotechnical and Automatic Engineering Department, Université, 8 Mai 1945, 24000 Guelma, Algeria. ²Electrical Engineering Laboratory, Electrical Engineering Departement, University of M'sila, M'sila, Algeria. ³IREENA Laboratory, Industrial Engineering Department, University of Nantes, Nantes, France. ⁴Piteşti University Centre, National University of Science and Technology POLITEHNICA Bucharest, 110040 Pitesti, Romania. ⁵Doctoral School, Polytechnic University of Bucharest, 313 SplaiuIIndependententei, 060042 Bucharest, Romania. ⁶İCSI Energy Department, National Research and Development Institute for Cryogenic and Isotopic Technologies, 240050 RamnicuValcea, Romania. email: ghadbane.houssameddine@univ-guelma.dz; said.barkat@univ-msila.dz</footnote>
<figure><img src="image_3.png" /></figure>

OPEN

# Energy management of electric vehicle using a new strategy based on slap swarm optimization and differential flatness control

Houssam Eddine Ghadbane<fnref n="1" />,*, Said Barkat²,*, Ali Djerioui², Azeddine Houari³, Mihai Oproescu⁴ & Nicu Bizon⁴,⁵,⁶

Optimal energy management of electric vehicles using slap swarm optimization and differential flatness control has been proposed. A battery-supercapacitor power system is adopted. Each source is connected in parallel to the DC-bus using DC-DC bidirectional converters and supplies a synchronous reluctance motor (SynRM) based drive. The proposed EMS fundamental forces lie in using a combination of complementary proprieties of two approaches, a Slap Swarm optimization Algorithm and Differential Flatness (DF). With a fast optimization mechanism, the Slap Swarm optimization algorithm allows adapting in real-time conditions the DF gains to optimize the system performances. On its side, DF uses predefined trajectories respecting the physical proprieties of the system, which is a powerful tool to guarantee the dynamic constraints of the sources when ensuring desired robust control proprieties. To check the feasibility and performance of the suggested EMS, comprehensive processor-in-the-loop co-simulations of the electric vehicle were carried out using the C2000 launchxlf28379d DSP board. The main goal of the proposed EMS is to guarantee the DC-bus stabilization, reducing the DC-bus voltage ripples (Δv = 5 V) and the voltage overshoots 15 V (3.2%), respect the source dynamics, and satisfy the SynRM motor power demand. Furthermore, the algorithm minimizes induced harmonics by the drive (10.49%), reducing the battery current ripple by 17.15A, thereby enhancing the battery lifecycle.

Nowadays, the influence of the transportation system's growth on global warming and climate change has become more and more tangible and has impacted our everyday lives<fnref n="1" />. According to², integrating the road transport sector into the EU ETS (European Emissions Trading System) is expected to be a cost-effective CO₂ reduction method. The transport sector accounted for roughly 27% of total EU greenhouse gas (GHG) emissions in 2017, with light commercial vehicles (LCVs) contributing to approximately 9% of EU transport GHG emissions³. According to the European Environment Agency (EEA), average CO₂ emissions from new light-duty vehicles (LDVs) rose in 2017 for the first time since 2010 and continued to rise in 2018⁴. The European Commission published a proposal to outlaw the sale of petrol and diesel automobiles in the EU from 2035⁵. In this context, transportation decarbonization is primarily recognized as one of the priorities for considerably decreasing harmful emissions.

The facts presented here confirm that the automobile industry has adopted the electrified mobility vision. Several automakers have openly stated their commitment to an electric vehicle future, expecting the EV market to be poised for rapid expansion. For example, Volvo Company stated that cars with internal combustion engines have no long-term future. Volvo is firmly dedicated to becoming an all-electric vehicle manufacturer, with the transition expected to occur by 2030⁶. By 2035, General Motors intends to eliminate tailpipe emissions from new

<footnote n="1">¹Electrical Engineering Laboratory of Guelma (LGEG), Electrotechnical and Automatic Engineering Department, Université, 8 Mai 1945, 24000 Guelma, Algeria. ²Electrical Engineering Laboratory, Electrical Engineering Departement, University of M'sila, M'sila, Algeria. ³IREENA Laboratory, Industrial Engineering Department, University of Nantes, Nantes, France. ⁴Piteşti University Centre, National University of Science and Technology POLITEHNICA Bucharest, 110040 Pitesti, Romania. ⁵Doctoral School, Polytechnic University of Bucharest, 313 SplaiulIndependentei, 060042 Bucharest, Romania. ⁶İCSI Energy Department, National Research and Development Institute for Cryogenic and Isotopic Technologies, 240050 RamnicuValcea, Romania. ⁷email: ghadbane.houssameddine@univ-guelma.dz; said.barkat@univ-msila.dz</footnote>
light-duty cars⁷. Following a year in which Volkswagen electric vehicle sales quadrupled, the company noticed that 2020 represented a turning point in customer mood⁸. EV conventional power system mainly involves batteries and power electronic converters⁹. It is well established that these power converters' currents contain harmonics. The detrimental effect of current harmonics on the aging of Li-ion cells was studied based on the results of the tests conducted by the authors in¹⁰⁻¹⁵.

To provide better performance in terms of power quality, durability, and reliability, hybrid power systems (HPSs) based EVs are suggested in the literature¹⁶⁻¹⁸. HPS comprises a battery as the primary source and a supercapacitor (SC) as an auxiliary. The battery provides most of the motor power due to its high energy density. At the same time, the SC enhances the power quality by supplying the transient periods due to its high power density. The electric machine has played a significant role throughout the evolution of EVs, regardless of the precise types of EVs provided. Different electric machines have been used in the EV design; each has advantages and disadvantages. Synchronous Reluctance Motor (SynRM) has recently gained more attention in EV applications due to its increased efficiency, reduced weight, and moderate cost¹⁹⁻²³. In addition, SynRM can provide torque enhancement by 34% more than the classical machines, which is suitable for the EV industry¹⁹. Despite the outstanding advantages described above, SynRM has several relevant disadvantages for traction applications, including torque ripple and SynRM current harmonics that affect battery lifetime (battery aging)²¹. The primary challenge is to provide an appropriate energy management strategy (EMS) that enhances the performance of the HPS and extends the battery lifecycle while considering different constraints.

However, the main focus lies in evaluating the efficacy of the chosen EMS. To attain this objective, multiple EMSs have been documented in academic literature. Broadly, these EMSs can be categorized into optimization, rule-based, and learning-based techniques, as mentioned in²⁴,²⁵. Optimization-driven approaches employ optimization theory tools to solve the issue, aiming for an optimal power load distribution across the battery-supercapacitor storage system to maximize component lifespans. The methods focused on optimization can be classified as offline (global scale) or online (real-time scale). Offline optimization EMS involves determining the best control solution for a predefined condition, like a speed profile. There are two methods for this: Direct methods encompass disciplined optimum control; Indirect methods involve approaches such as the calculus of variations²⁴, Pontryagin's Maximum Principle (PMP)²⁶, Pontryagin's Minimum Principle²⁷, stochastic dynamic programming (SDP)²⁸, and dynamic programming (DP)²⁹. The real-time application of online optimization includes determining the most efficient energy distribution in a hybrid system using specified data. The cost function factors in the system's current status, operational costs, and emissions. Model predictive control (MPC)³⁰, equivalent consumption minimization strategy (ECMS)³¹, and external energy maximization strategy (EEMS)³² fall under this category. In artificial intelligence, learning-based techniques, particularly reinforcement learning (RL)³³ and deep learning (DL)³⁴, leverage machine learning advancements. These methods have exhibited effectiveness across diverse domains, notably in image categorization, leading to widespread implementation in energy management³⁵,³⁶. However, the unavailability of databases required for training models poses a challenge, given the lack of extensive research on this emerging subject. These databases offer no assurance of compatibility with data beyond the provided training set. Rule-based strategies are management strategies based on a sequence of IF-THEN situations. This kind of management strategy can be categorized into two types: deterministic strategies, such as state machine control (SMC) strategy²⁷,²⁸, and intelligent strategies, such as fuzzy logic-based EMS²⁹,³⁰. The fundamental disadvantage of these strategies is that developing their roles requires the designer's experience, which is not always available. They also suffer from problems related to the abrupt transitions between various operating modes that present a real challenge for conventional controllers to keep the desired power quality and system stability. A linear proportional-integral (PI) based management strategy was reported in²⁴⁻³⁵. Authors in³⁷ proposed a PI-based voltage control strategy for the common DC-bus of hybrid EVs. However, PI controller design based on linearization around an operating point suffers from robustness problems. The flatness-based management strategy (FLAT) has been widely used to manage the power flow in a multi-sources power system³⁸,³⁹. This strategy can provide excellent performance, but defining the trajectory parameters is challenging; the definition of flat outputs and the dependence on the reference model decreases these performances in case of essential parameter changes. The traditional approaches may offer restricted effectiveness. This research aims to evaluate the optimization of trajectory generation parameters utilizing a metaheuristic optimization technique known as the Salp Swarm Algorithm (SSA). The main benefits of SSA over other optimization strategies are faster solution convergence and fewer parameters with faster-resolving capability⁴⁰. Based on the advantages of the above-mentioned optimization algorithms, this paper reports a new optimized flatness-based management strategy to manage the power flow in an EV power system while reducing battery aging. The considered hybrid power system is based on a Li-ion battery and supercapacitor as energy storage elements. The main goals of the suggested EMS are to provide DC-bus stabilization, respect the dynamics of the sources, and satisfy the SynRM motor power demand. Additionally, the algorithm allows the minimization of the effect of harmonics generated by the motor on the battery current. Thus, the main contributions can be listed as follows:

* Stabilizing the DC-bus voltage;

* Ensuring the proper use of the power flows between the battery and SC while respecting the SoC limits;

* Preventing the SynRM motor from the harmonics that impact the energy storage systems' lifetimes;

* Reducing the battery current harmonics induced by the electric drive;

* Minimizing the battery current ripples generated by both the SynRM and PWM inverter;

* Combining different strategies to obtain the best performances.

The control algorithm and EMS for the proposed hybrid power system are validated by processor-in-the-loop (PIL) technology in Matlab/Simulink using embedded Matlab functions. The final generated control algorithm
light-duty cars⁷. Following a year in which Volkswagen electric vehicle sales quadrupled, the company noticed that 2020 represented a turning point in customer mood⁸. EV conventional power system mainly involves batteries and power electronic converters⁹. It is well established that these power converters' currents contain harmonics. The detrimental effect of current harmonics on the aging of Li-ion cells was studied based on the results of the tests conducted by the authors in¹⁰⁻¹⁵.

To provide better performance in terms of power quality, durability, and reliability, hybrid power systems (HPSs) based EVs are suggested in the literature¹⁶⁻¹⁸. HPS comprises a battery as the primary source and a supercapacitor (SC) as an auxiliary. The battery provides most of the motor power due to its high energy density. At the same time, the SC enhances the power quality by supplying the transient periods due to its high power density. The electric machine has played a significant role throughout the evolution of EVs, regardless of the precise types of EVs provided. Different electric machines have been used in the EV design; each has advantages and disadvantages. Synchronous Reluctance Motor (SynRM) has recently gained more attention in EV applications due to its increased efficiency, reduced weight, and moderate cost¹⁹⁻²³. In addition, SynRM can provide torque enhancement by 34% more than the classical machines, which is suitable for the EV industry¹⁹. Despite the outstanding advantages described above, SynRM has several relevant disadvantages for traction applications, including torque ripple and SynRM current harmonics that affect battery lifetime (battery aging)²¹. The primary challenge is to provide an appropriate energy management strategy (EMS) that enhances the performance of the HPS and extends the battery lifecycle while considering different constraints.

However, the main focus lies in evaluating the efficacy of the chosen EMS. To attain this objective, multiple EMSs have been documented in academic literature. Broadly, these EMSs can be categorized into optimization, rule-based, and learning-based techniques, as mentioned in²⁴,²⁵. Optimization-driven approaches employ optimization theory tools to solve the issue, aiming for an optimal power load distribution across the battery-supercapacitor storage system to maximize component lifespans. The methods focused on optimization can be classified as offline (global scale) or online (real-time scale). Offline optimization EMS involves determining the best control solution for a predefined condition, like a speed profile. There are two methods for this: Direct methods encompass disciplined optimum control; Indirect methods involve approaches such as the calculus of variations²⁴, Pontryagin's Maximum Principle (PMP)²⁶, Pontryagin's Minimum Principle²⁷, stochastic dynamic programming (SDP)²⁸, and dynamic programming (DP)²⁹. The real-time application of online optimization includes determining the most efficient energy distribution in a hybrid system using specified data. The cost function factors in the system's current status, operational costs, and emissions. Model predictive control (MPC)³⁰, equivalent consumption minimization strategy (ECMS)³¹, and external energy maximization strategy (EEMS)³² fall under this category. In artificial intelligence, learning-based techniques, particularly reinforcement learning (RL)³³ and deep learning (DL)³⁴, leverage machine learning advancements. These methods have exhibited effectiveness across diverse domains, notably in image categorization, leading to widespread implementation in energy management³⁵,³⁶. However, the unavailability of databases required for training models poses a challenge, given the lack of extensive research on this emerging subject. These databases offer no assurance of compatibility with data beyond the provided training set. Rule-based strategies are management strategies based on a sequence of IF-THEN situations. This kind of management strategy can be categorized into two types: deterministic strategies, such as state machine control (SMC) strategy²⁷,²⁸, and intelligent strategies, such as fuzzy logic-based EMS²⁹,³⁰. The fundamental disadvantage of these strategies is that developing their roles requires the designer's experience, which is not always available. They also suffer from problems related to the abrupt transitions between various operating modes that present a real challenge for conventional controllers to keep the desired power quality and system stability. A linear proportional-integral (PI) based management strategy was reported in²⁴⁻³⁵. Authors in³⁷ proposed a PI-based voltage control strategy for the common DC-bus of hybrid EVs. However, PI controller design based on linearization around an operating point suffers from robustness problems. The flatness-based management strategy (FLAT) has been widely used to manage the power flow in a multi-sources power system³⁸,³⁹. This strategy can provide excellent performance, but defining the trajectory parameters is challenging; the definition of flat outputs and the dependence on the reference model decreases these performances in case of essential parameter changes. The traditional approaches may offer restricted effectiveness. This research aims to evaluate the optimization of trajectory generation parameters utilizing a metaheuristic optimization technique known as the Salp Swarm Algorithm (SSA). The main benefits of SSA over other optimization strategies are faster solution convergence and fewer parameters with faster-resolving capability⁴⁰. Based on the advantages of the above-mentioned optimization algorithms, this paper reports a new optimized flatness-based management strategy to manage the power flow in an EV power system while reducing battery aging. The considered hybrid power system is based on a Li-ion battery and supercapacitor as energy storage elements. The main goals of the suggested EMS are to provide DC-bus stabilization, respect the dynamics of the sources, and satisfy the SynRM motor power demand. Additionally, the algorithm allows the minimization of the effect of harmonics generated by the motor on the battery current. Thus, the main contributions can be listed as follows:

* Stabilizing the DC-bus voltage;

* Ensuring the proper use of the power flows between the battery and SC while respecting the SoC limits;

* Preventing the SynRM motor from the harmonics that impact the energy storage systems' lifetimes;

* Reducing the battery current harmonics induced by the electric drive;

* Minimizing the battery current ripples generated by both the SynRM and PWM inverter;

* Combining different strategies to obtain the best performances.

The control algorithm and EMS for the proposed hybrid power system are validated by processor-in-the-loop (PIL) technology in Matlab/Simulink using embedded Matlab functions. The final generated control algorithm
codes are embedded in a real target DSP (C2000 launchxl-f28379d). In this, the integrated DSP control algo-
rithm manages all components of the hybrid power system (Bat-SC-SynRM), modeled in the host PC using
SimPowerSystems (SPS) Toolbox.

The remainder of this work is arranged as follows: Sect. 2 presents modeling the electric traction chain’s parts. Section 3 discusses the suggested EMS based on the optimal differential flatness. The results and interpretations of the co-simulations are presented in Sect. 4. The study’s conclusion is presented in Sect. 5.

System description and modeling

Detailing the study and the proposed diagram

Finding the specific numerical values for trajectory generation parameters presents a challenge because there
is not an exact model replicating the physical system. To address this, an approach involves improving these
parameters through metaheuristic optimization algorithms. The main concept involves creating random poten-
tial solutions within a confined search area. These solutions will be directed to the HPS, where their behavior
will be evaluated by calculating the integral square error (ISE) between the reference and actual DC bus and
supercapacitor voltage and battery and the supercapacitor current. The optimizers will adjust and update the
potential solutions based on this ISE, which represents the fitness value. The optimized adaptive Flatness-based
EMS is described in Figure 1.

HPS topology

Three hybrid power system topologies for electric vehicles, as illustrated in Fig. 2, have been commonly reported in the literature. These topologies with distinctive qualities can be categorized as passive, semi-active, and active<sup>41</sup>.

The architecture of the studied HPS is based on the fully active topology, as shown in Fig. 3. Since each ESS is controlled independently, this architecture provides the maximum level of controllability. It makes it easier for the energy management plan to use the complementing traits of HPS and HES. Additionally, it supports

<figure><img src="image_4.png" /><figcaption>Figure 1. SSA-FLAT-based EMS.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The HPS topologies: (a) passive topology; (b,c) semi-active topological structure; (d) fully active topology.</figcaption></figure>
codes are embedded in a real target DSP (C2000 launchxl-f28379d). In this, the integrated DSP control algo-
rithm manages all components of the hybrid power system (Bat-SC-SynRM), modeled in the host PC using
SimPowerSystems (SPS) Toolbox.

The remainder of this work is arranged as follows: Sect. 2 presents modeling the electric traction chain’s parts. Section 3 discusses the suggested EMS based on the optimal differential flatness. The results and interpretations of the co-simulations are presented in Sect. 4. The study’s conclusion is presented in Sect. 5.

System description and modeling
Detailing the study and the proposed diagram

Finding the specific numerical values for trajectory generation parameters presents a challenge because there
is not an exact model replicating the physical system. To address this, an approach involves improving these
parameters through metaheuristic optimization algorithms. The main concept involves creating random poten-
tial solutions within a confined search area. These solutions will be directed to the HPS, where their behavior
will be evaluated by calculating the integral square error (ISE) between the reference and actual DC bus and
supercapacitor voltage and battery and the supercapacitor current. The optimizers will adjust and update the
potential solutions based on this ISE, which represents the fitness value. The optimized adaptive Flatness-based
EMS is described in Figure 1.

HPS topology

Three hybrid power system topologies for electric vehicles, as illustrated in Fig. 2, have been commonly reported in the literature. These topologies with distinctive qualities can be categorized as passive, semi-active, and active<sup>41</sup>.

The architecture of the studied HPS is based on the fully active topology, as shown in Fig. 3. Since each ESS is controlled independently, this architecture provides the maximum level of controllability. It makes it easier for the energy management plan to use the complementing traits of HPS and HES. Additionally, it supports

<figure><img src="image_4.png" /><figcaption>Figure 1. SSA-FLAT-based EMS.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** The HPS topologies: (a) passive topology; (b,c) semi-active topological structure; (d) fully active topology.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>Figure 3. Schematic structure of the proposed EV.</figcaption></figure>

implementing a wide range of control mechanisms$^{42}$. The traction chain comprises two energy storage systems (battery and SC), their converters (bidirectional DC-DC), and a synchronous reluctance motor drive.

## Power system modeling

The switching and average models are the two DC-DC converters models extensively utilized in the literature. The switching model is primarily used to investigate different forms of pulse width-modulated systems in terms of power losses and harmonics. It requires a short sampling period, making simulation time-consuming; hence, the following average model was adopted in this study.

$$ \left\{ \begin{aligned} L_{bat} \frac{di_{bat}}{dt} &= V_{bat} - V_{bus}d_{batt} - r_{Bat}i_{batt} \\ L_{sc} \frac{di_{sc}}{dt} &= V_{sc} - V_{bus}d_{sc} - r_{sc}i_{sc} \\ c_{bus} \frac{dV_{bus}}{dt} &= d_{batt}i_{bat} + d_{sc}i_{sc} - i_{Load} \end{aligned} \right. \quad (1) $$

where $V_{sc}$ and $V_{bus}$ are the supercapacitor and DC-bus voltage, respectively; $r_{sc}$ and $r_{batt}$ denote the internal converter resistors; $L_{sc}$ and $L_{batt}$ are the converter inductors; $d_{sc}$ and $d_{batt}$ represent the duty cycle ratios; $C_{bus}$ represents the DC-bus capacitance; $i_{Load}$, $i_{sc}$, and $i_{batt}$ are the load, the supercapacitor, and the battery currents, respectively.

The electromagnetic energy stored in the DC bus capacitance $E_{bus}$ and the capacitive energy of the SC $E_{sc}$ are given by:

$$ E_{Bus} = 0.5 c_{bus} V_{bus}^2 \quad (2) $$

$$ E_{sc} = 0.5 c_{sc} V_{sc}^2 \quad (3) $$

The total stored electromagnetic energy $E_T$ is expressed as follows:

$$ E_T = 0.5 c_{bus} V_{bus}^2 + 0.5 c_{sc} V_{sc}^2 \quad (4) $$

According to Fig. 1, the differential equation for power balance is$^{43}$:

$$ \dot{E}_{Bus} = P_{Bato} + P_{sc} - P_{Load} \quad (5) $$

The supplied powers from the sources, including converter losses, are expressed by $P_{Bato}$ and $P_{sc}$ as follows:

$$ P_{Bato} = P_{Bat} - r_{Bat} \left( \frac{P_{Bat}}{V_{Bat}} \right)^2 \quad (6) $$

$$ P_{sc} = P_{sc} - r_{sc} \left( \frac{P_{sc}}{v_{sc}} \right)^2 \quad (7) $$

where

$$ P_{sc} = V_{sc} i_{sc} = i_{sc} \left( 2E_{sc} / c_{sc} \right)^{1/2} \quad (8) $$
<figure><img src="image_6.png" /><figcaption>Figure 3. Schematic structure of the proposed EV.</figcaption></figure>

implementing a wide range of control mechanisms$^{42}$. The traction chain comprises two energy storage systems (battery and SC), their converters (bidirectional DC-DC), and a synchronous reluctance motor drive.

## Power system modeling

The switching and average models are the two DC-DC converters models extensively utilized in the literature. The switching model is primarily used to investigate different forms of pulse width-modulated systems in terms of power losses and harmonics. It requires a short sampling period, making simulation time-consuming; hence, the following average model was adopted in this study.

$$ egin{cases} L_{bat} \frac{di_{bat}}{dt} = V_{bat} - V_{bus}d_{batt} - r_{bat}i_{batt} \\ L_{sc} \frac{di_{sc}}{dt} = V_{sc} - V_{bus}d_{sc} - r_{sc}i_{sc} \\ c_{bus} \frac{dV_{bus}}{dt} = d_{batt}i_{bat} + d_{sc}i_{sc} - i_{Load} \end{cases} \quad (1) $$

where $V_{sc}$ and $V_{bus}$ are the supercapacitor and DC-bus voltage, respectively; $r_{sc}$ and $r_{batt}$ denote the internal converter resistors; $L_{sc}$ and $L_{batt}$ are the converter inductors; $d_{sc}$ and $d_{batt}$ represent the duty cycle ratios; $C_{bus}$ represents the DC-bus capacitance; $i_{Load}$, $i_{sc}$, and $i_{batt}$ are the load, the supercapacitor, and the battery currents, respectively.

The electromagnetic energy stored in the DC bus capacitance $E_{bus}$ and the capacitive energy of the SC $E_{sc}$ are given by:

$$ E_{Bus} = 0.5 c_{bus} V_{bus}^2 \quad (2) $$

$$ E_{sc} = 0.5 c_{sc} V_{sc}^2 \quad (3) $$

The total stored electromagnetic energy $E_T$ is expressed as follows:

$$ E_T = 0.5 c_{bus} V_{bus}^2 + 0.5 c_{sc} V_{sc}^2 \quad (4) $$

According to Fig. 1, the differential equation for power balance is$^{43}$:

$$ \dot{E}_{Bus} = P_{Bato} + P_{sc0} - P_{Load} \quad (5) $$

The supplied powers from the sources, including converter losses, are expressed by $P_{Bato}$ and $P_{sc0}$ as follows:

$$ P_{Bato} = P_{Bat} - r_{Bat} \left( \frac{P_{Bat}}{V_{Bat}} \right)^2 \quad (6) $$

$$ P_{sc0} = P_{sc} - r_{sc} \left( \frac{P_{sc}}{v_{sc}} \right)^2 \quad (7) $$

where

$$ P_{sc} = V_{sc}i_{sc} = i_{sc} \left( 2E_{sc} / c_{sc} \right)^{1/2} \quad (8) $$
$$P_{Load} = V_{bus}i_{Load} = i_{Load}(2E_{Bus}/c_{bus})^{1/2} \quad (9)$$

## Synchronous reluctance motor and vehicle dynamics modeling

There are currently several main rotor structures of SynRM: solid rotor, flux barrier rotor, axially laminated rotor, and magnet-assisted rotor. The structures suitable for traction applications have flux barriers since they meet performance, robustness, cost, and manufacturing requirements²⁰. The SynRM electromechanical model is presented in the (d-q) frame as follows:

$$\frac{d}{dt} \begin{bmatrix} i_{sd} \\ i_{sq} \\ \Omega \\ \theta \end{bmatrix} = \begin{bmatrix} \frac{-R_s}{L_d} i_{sd} + \frac{p L_d}{L_d} i_{sq} \Omega \\ \frac{-R_s}{L_d} i_{sq} - \frac{p L_d}{L_d} i_{sd} \Omega \\ \frac{3p(L_d - L_q)i_{sd}i_{sq}}{2J} - \frac{f}{J} \Omega - T_L \\ \frac{p \Omega}{J} \end{bmatrix} + \begin{bmatrix} \frac{1}{L_d} & 0 \\ 0 & \frac{1}{L_q} \\ 0 & 0 \\ 0 & 0 \end{bmatrix} \begin{bmatrix} v_{sd} \\ v_{sq} \end{bmatrix} \quad (10)$$

where $v_{sq}$ and $v_{sd}$ are the stator's quadrature and direct axis voltages, respectively. $i_{sq}$ and $i_{sd}$ represent their corresponding currents. The $L_d$ and $L_q$ are the d-q magnetizing inductances. $\Omega$ represents the mechanical rotation speed, $R_s$ represents the stator resistance, $p$ denotes the number of pair poles, $T_L$, is the rotor's torque load, $f$ and $J$ represent the viscous friction coefficient, and inertia moment.

The equation representing the dynamics of an electric motor is formulated as follows¹⁸:

$$T_e - T_{Load} = J_{ve} \frac{d\Omega}{dt} \quad (11)$$

The torque of the load is provided by:

$$T_{Load} = \frac{T_{Lwheel}}{i} = \frac{r}{i} \left[ (m_v g \mu \operatorname{sign}(v)) + (0.5 \gamma C_d A_f (v \pm v_w)^2) \pm (m_v g \sin \alpha) \pm (k_m m_v \frac{dv}{dt}) \right] \quad (12)$$

where $T_{Lwheel}$, $i$, and $r$ represent the load torque on the wheels, the transmission ratio, and the wheel radius, $v$ and $v_w$ are the vehicle and wind speed, respectively, $\alpha$ represents the angle of the slope, $m_v$ is vehicle mass, $\mu, k_m$, and $C_d$ are the tire rolling resistance, the rotational inertia, and aerodynamic drag coefficient, respectively, $\gamma$ and $A_f$ represent the air density and the frontal vehicle area. Respectively, $g$ represents the earth's gravity.

## Classification and comparative analysis of the energy management approach

EMSs can be partitioned into three main methods: optimization-based methods, rule-based methods, or more recently, learning-based methods, as depicted in Fig. 4.

<figure><img src="image_7.png" /><figcaption>Figure 4. Classification of EMS of HPS.</figcaption></figure>
$$P_{Load} = V_{bus}i_{Load} = i_{Load}(2E_{Bus}/c_{bus})^{1/2} \quad (9)$$

## Synchronous reluctance motor and vehicle dynamics modeling

There are currently several main rotor structures of SynRM: solid rotor, flux barrier rotor, axially laminated rotor, and magnet-assisted rotor. The structures suitable for traction applications have flux barriers since they meet performance, robustness, cost, and manufacturing requirements²⁰. The SynRM electromechanical model is presented in the (d-q) frame as follows:

$$\frac{d}{dt} \begin{bmatrix} i_{sd} \\ i_{sq} \\ \Omega \\ \theta \end{bmatrix} = \begin{bmatrix} \frac{-R_s}{L_d} i_{sd} + \frac{p L_d}{L_d} i_{sq} \Omega \\ \frac{-R_s}{L_d} i_{sq} - \frac{p L_d}{L_d} i_{sd} \Omega \\ \frac{3p(L_d - L_q)i_{sd}i_{sq}}{2J} - \frac{f}{J} \Omega - T_L \\ \frac{p \Omega}{J} \end{bmatrix} + \begin{bmatrix} \frac{1}{L_d} & 0 \\ 0 & \frac{1}{L_q} \\ 0 & 0 \\ 0 & 0 \end{bmatrix} \begin{bmatrix} v_{sd} \\ v_{sq} \end{bmatrix} \quad (10)$$

where $v_{sq}$ and $v_{sd}$ are the stator's quadrature and direct axis voltages, respectively. $i_{sq}$ and $i_{sd}$ represent their corresponding currents. The $L_d$ and $L_q$ are the d-q magnetizing inductances. $\Omega$ represents the mechanical rotation speed, $R_s$ represents the stator resistance, $p$ denotes the number of pair poles, $T_L$, is the rotor's torque load, $f$ and $J$ represent the viscous friction coefficient, and inertia moment.

The equation representing the dynamics of an electric motor is formulated as follows¹⁸:

$$T_e - T_{Load} = J_{ve} \frac{d\Omega}{dt} \quad (11)$$

The torque of the load is provided by:

$$T_{Load} = \frac{T_{Lwheel}}{i} = \frac{r}{i} \left[ (m_v g \mu \operatorname{sign}(v)) + (0.5 \gamma C_d A_f (v \pm v_w)^2) \pm (m_v g \sin \alpha) \pm (k_m m_v \frac{dv}{dt}) \right] \quad (12)$$

where $T_{Lwheel}$, $i$, and $r$ represent the load torque on the wheels, the transmission ratio, and the wheel radius, $v$ and $v_w$ are the vehicle and wind speed, respectively, $\alpha$ represents the angle of the slope, $m_v$ is vehicle mass, $\mu, k_m$, and $C_d$ are the tire rolling resistance, the rotational inertia, and aerodynamic drag coefficient, respectively, $\gamma$ and $A_f$ represent the air density and the frontal vehicle area. Respectively, $g$ represents the earth's gravity.

## Classification and comparative analysis of the energy management approach

EMSs can be partitioned into three main methods: optimization-based methods, rule-based methods, or more recently, learning-based methods, as depicted in Fig. 4.

<figure><img src="image_7.png" /><figcaption>Figure 4. Classification of EMS of HPS.</figcaption></figure>
<table><caption>Table 1. Comparative analysis of EMS⁴⁷.</caption><thead><tr><td>EM strategy</td><td>Advantages</td><td>Limitations</td><td>Refs.</td></tr></thead><tbody><tr><td>Deterministic rule-based</td><td>Simple, reliable, and robust, low computation complexity, easy to implement</td><td>Not adaptive and poor parametric calibration</td><td>48</td></tr><tr><td>Frequency separation</td><td>Simple and easy implementation</td><td>Not robust, filter design is difficult</td><td>49</td></tr><tr><td>Fuzzy logic Control</td><td>Robust and good with model uncertainties and state variations, real-time implementable</td><td>Dependent on membership functions, optimal control is not guaranteed</td><td>50</td></tr><tr><td>Dynamic programming</td><td>Global optimal solution is found, optimal control, easy to solve nonlinear optimization problems</td><td>Not real-time implementable requires heavy computational burden</td><td>51</td></tr><tr><td>Pontryagin's principle</td><td>Easy to adapt and simple to implement, no additional controllers are required</td><td>The computation burden is heavy</td><td>52,53</td></tr><tr><td>Instantaneous optimization</td><td>The optimal value is found at each instant</td><td>Not guaranteed to be optimal,</td><td>54,55</td></tr><tr><td>Model predictive control</td><td>Has potential for real-time implementations; Easy to handle constraints directly in the design procedure</td><td>Model accuracy may be compromised by using a linearized model May need large memory for heavy computations</td><td>30</td></tr><tr><td>Reinforcement Learning (RL)</td><td>RL can adapt to unknown environments; DL handles complex nonlinear relationships</td><td>RL requires significant exploration time. DL might lack explainability and require extensive data</td><td>56</td></tr><tr><td>Neural network</td><td>It is more robust to new information, real time implementable</td><td>Large amount of training data is required; stability is not guaranteed</td><td>57</td></tr></tbody></table>

## The suggested EMS

The suggested EMS aims to control the motor speed effectively and stabilize the voltage of the DC bus. Indirect Field Oriented Control (IFOC) is used to control the motor speed. This strategy provides the d-q voltage references for the State Space Vector PWM (SVPWM) system. Hybridization between the differential flatness theory and metaheuristic optimization is used to control the bus voltage optimally.

## SynRM drive IFOC design

The SynRM's Indirect Field Oriented Control (IFOC), as shown in Fig. 5, comprises three control loops: an external speed control loop, an internal current control loop for $i_{sd}$ and another loop for $i_{sq}$. The first loop output generates the quadrature reference current, which is compared to the current value resulting from the current measurements. Their error is applied to the input of the $i_{sq}$ current controller. Similarly, there is a current regulation loop of direct current $i_{sd}$. The outputs of the two internal loops $i_{sq}$ and $i_{sd}$ are applied to a decoupling block generating the reference voltages $v_{sq}^{ref}$ and $v_{sd}^{ref}$. By passing from the reference (d-q) to the reference $(\alpha, \beta)$, the two voltages reference $v_{sa}^{ref}$, $v_{sb}^{ref}$ required by the Space Vector Modulation (SVM) block.

## The proposed optimal differential flatness-based EMS

The proposed EMS is based on optimal differential flatness. It is divided into two parts: a lower-level controller and a higher-level controller. The higher-level controller generates the power reference for each source. The power

<figure><img src="image_8.png" /><figcaption>Figure 5. Diagram of SynRM Field-oriented control.</figcaption></figure>
<table><caption>**Table 1.** Comparative analysis of EMS⁴⁷.</caption><thead><tr><td>EM strategy</td><td>Advantages</td><td>Limitations</td><td>Refs.</td></tr></thead><tbody><tr><td>Deterministic rule-based</td><td>Simple, reliable, and robust, low computation complexity, easy to implement</td><td>Not adaptive and poor parametric calibration</td><td>48</td></tr><tr><td>Frequency separation</td><td>Simple and easy implementation</td><td>Not robust, filter design is difficult</td><td>49</td></tr><tr><td>Fuzzy logic Control</td><td>Robust and good with model uncertainties and state variations, real-time implementable</td><td>Dependent on membership functions, optimal control is not guaranteed</td><td>50</td></tr><tr><td>Dynamic programming</td><td>Global optimal solution is found, optimal control, easy to solve nonlinear optimization problems</td><td>Not real-time implementable requires heavy computational burden</td><td>51</td></tr><tr><td>Pontryagin's principle</td><td>Easy to adapt and simple to implement, no additional controllers are required</td><td>The computation burden is heavy</td><td>52,53</td></tr><tr><td>Instantaneous optimization</td><td>The optimal value is found at each instant</td><td>Not guaranteed to be optimal,</td><td>54,55</td></tr><tr><td>Model predictive control</td><td>Has potential for real-time implementations; Easy to handle constraints directly in the design procedure</td><td>Model accuracy may be compromised by using a linearized model May need large memory for heavy computations</td><td>30</td></tr><tr><td>Reinforcement Learning (RL)</td><td>RL can adapt to unknown environments; DL handles complex nonlinear relationships</td><td>RL requires significant exploration time. DL might lack explainability and require extensive data</td><td>56</td></tr><tr><td>Neural network</td><td>It is more robust to new information, real time implementable</td><td>Large amount of training data is required; stability is not guaranteed</td><td>57</td></tr></tbody></table>

## The suggested EMS

The suggested EMS aims to control the motor speed effectively and stabilize the voltage of the DC bus. Indirect Field Oriented Control (IFOC) is used to control the motor speed. This strategy provides the d-q voltage references for the State Space Vector PWM (SVPWM) system. Hybridization between the differential flatness theory and metaheuristic optimization is used to control the bus voltage optimally.

## SynRM drive IFOC design

The SynRM's Indirect Field Oriented Control (IFOC), as shown in Fig. 5, comprises three control loops: an external speed control loop, an internal current control loop for $i_{sd}$ and another loop for $i_{sq}$. The first loop output generates the quadrature reference current, which is compared to the current value resulting from the current measurements. Their error is applied to the input of the $i_{sq}$ current controller. Similarly, there is a current regulation loop of direct current $i_{sd}$. The outputs of the two internal loops $i_{sq}$ and $i_{sd}$ are applied to a decoupling block generating the reference voltages $v_{sq}^{ref}$ and $v_{sd}^{ref}$. By passing from the reference (d-q) to the reference $(\alpha, \beta)$, the two voltages reference $v_{sa}^{ref}$, $v_{sb}^{ref}$ required by the Space Vector Modulation (SVM) block.

## The proposed optimal differential flatness-based EMS

The proposed EMS is based on optimal differential flatness. It is divided into two parts: a lower-level controller and a higher-level controller. The higher-level controller generates the power reference for each source. The power

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Diagram of SynRM Field-oriented control.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Figure 6. Proposed optimal differential flatness-based EMS.</figcaption></figure>

references created by the inverse dynamic in Eq. (16) and Eq. (18) are divided by the measured supercapacitor and battery voltages to generate the reference currents for the supercapacitor and battery converters. The lower-level controller adjusts for DC-bus current variations and compensates for the DC-bus current harmonics using a supercapacitor, which results in improved energy quality and battery lifecycle enhancement. Figure 6, illustrates the operational control structure.

### Energy regulation

In the considered HPS, the flat model is expressed by its flat output $y = [y_1 y_2]^T$, control variable $u = [u_1 u_2]^T$, and state variable $x = [x_1 x_2]^T$,

Where:

$$y = \begin{bmatrix} E_{Bus} \\ E_{T} \end{bmatrix}, u = \begin{bmatrix} P_{SC}^{ref} \\ P_{Bat}^{ref} \end{bmatrix}, x = \begin{bmatrix} V_{Bus} \\ V_{SC} \end{bmatrix} \quad (13)$$

From Eq. (2), the state variable $x_1$ representing $V_{bus}$ can be expressed as follows:

$$x_1 = (2y_1/c_{bus})^{1/2} = \lambda_1(y_1) \quad (14)$$

From Eq. (4), the state variable $x_2$ representing $V_{SC}$ can be expressed as follows:

$$x_2 = (2(y_2 - y_1)/c_{SC})^{1/2} = \lambda_2(y_1, y_2) \quad (15)$$

The SC power reference $P_{SCref}$, considered the first input control variable, $u_1$, is derived from Eqs. 5, 7, and 9.
<figure><img src="image_9.png" /><figcaption>Figure 6. Proposed optimal differential flatness-based EMS.</figcaption></figure>

references created by the inverse dynamic in Eq. (16) and Eq. (18) are divided by the measured supercapacitor and battery voltages to generate the reference currents for the supercapacitor and battery converters. The lower-level controller adjusts for DC-bus current variations and compensates for the DC-bus current harmonics using a supercapacitor, which results in improved energy quality and battery lifecycle enhancement. Figure 6, illustrates the operational control structure.

### Energy regulation

In the considered HPS, the flat model is expressed by its flat output $y = [y_1 y_2]^T$, control variable $u = [u_1 u_2]^T$, and state variable $x = [x_1 x_2]^T$,

Where:

$$y = \begin{bmatrix} E_{Bus} \\ E_{T} \end{bmatrix}, u = \begin{bmatrix} P_{SC}^{ref} \\ P_{Bat}^{ref} \end{bmatrix}, x = \begin{bmatrix} V_{Bus} \\ V_{SC} \end{bmatrix} \quad (13)$$

From Eq. (2), the state variable $x_1$ representing $V_{bus}$ can be expressed as follows:

$$x_1 = (2y_1/c_{bus})^{1/2} = \lambda_1(y_1) \quad (14)$$

From Eq. (4), the state variable $x_2$ representing $V_{SC}$ can be expressed as follows:

$$x_2 = (2(y_2 - y_1)/c_{SC})^{1/2} = \lambda_2(y_1, y_2) \quad (15)$$

The SC power reference $P_{SCref}$, considered the first input control variable, $u_1$, is derived from Eqs. 5, 7, and 9.
$$u_1 = 2P_{\text{SC max}} \left[ 1 - \left( 1 - \frac{\dot{y}_1 + \sqrt{\frac{2y_1}{c_{\text{bus}}}} i_{\text{Load}} - P_{\text{Bato}}}{P_{\text{SC max}}} \right)^{1/2} \right] = \vartheta_1(y_1, \dot{y}_1) = P_{\text{SCref}} \quad (16)$$

with $P_{\text{SCmax}}$ is the SC converter's maximum power limit. It can be defined as

$$P_{\text{SC max}} = V_{\text{SC}}^2 / 4r_{\text{SC}} \quad (17)$$

Concerning the second input, the power reference of the battery $P_{\text{Battref}}$ is calculated from Eqs. 5, 6, and 9 as follows:

$$u_2 = 2P_{\text{Bat max}} \left[ 1 - \left( 1 - \frac{\dot{y}_2 + \sqrt{\frac{2y_2}{c_{\text{bus}}}} i_{\text{Load}}}{P_{\text{Bat max}}} \right)^{1/2} \right] = \vartheta_2(y_1, \dot{y}_2) = P_{\text{Bat ref}} \quad (18)$$

where $P_{\text{Battmax}}$ is defined as the maximum limiting power from the battery converter.

$$P_{\text{Bat max}} = V_{\text{bat}}^2 / 4r_{\text{Bat}} \quad (19)$$

Above all, the flat output ($y_1 = E_{\text{bus}}$) is the most crucial variable to control. For this, a classical PI controller is used to ensure that this flat variable is under control. Assuming the SC control loop is significantly faster than the battery control loop$^{44}$, the DC-bus power stated in Eq. (7) can be approximated as:

$$\dot{E}_{\text{Bus}} = P_{\text{Sco}} \quad (20)$$

The transfer function is represented as a straightforward integrator. A PI regulator is used to control this part$^{44}$. Recognizing that $y_1 = E_{\text{bus}}$, as shown below:

$$\dot{y}_1 = \frac{1}{s} \left( k_p^{V_{\text{bus}}} + \frac{k_i^{V_{\text{bus}}}}{s} \right) (y_1 - y_{1-\text{ref}}) \quad (21)$$

where $k_p^{V_{\text{bus}}}, k_i^{V_{\text{bus}}}$ are the integral and the proportional gains chosen so that the closed loop characteristic polynomial is expressed as follows:

$$p(s) = s^2 + \lambda_1 s + \lambda_0 \quad (22)$$

Clearly, the error $e_1 = y_1 - y_{1-\text{ref}}$ meets:

$$\ddot{e}_1 + k_p^{V_{\text{bus}}} \dot{e}_1 + k_i^{V_{\text{bus}}} e_1 = 0 \quad (23)$$

By matching the characteristic polynomial $p(s)$ to a desired one $p_{\text{des}}(s)$, given by Eq. (22), with pre-specified root positions, an adequate choice of controller parameters can be calculated by Eq. (25) and Eq. (26).

$$p_{\text{des}}(s) = s^2 + 2\xi\omega_n s + \omega_n^2 \quad (24)$$

$$k_p^{V_{\text{bus}}} = 2\xi\omega_n \quad (25)$$

$$k_i^{V_{\text{bus}}} = \omega_n^2 \quad (26)$$

where $\omega_n$ represents the natural frequency, and $\xi$ represents the dumping ratio.

The SC energy control loop depends on total energy management. A linearizing feedback control rule is used to achieve an exponential asymptotic tracking of the trajectory$^{44}$ as follows:

$$(\dot{y}_2 - \dot{y}_{2-\text{ref}}) + k_p^{V_{\text{SC}}} (y_2 - y_{2-\text{ref}}) = 0 \quad (27)$$

$$k_p^{V_{\text{SC}}} = \omega_n^2 \quad (28)$$

where $k_p^{V_{\text{SC}}}$ is the proportional gain of the SC voltage controller.

### Current control

The flat output $y = [y_3 y_4]^T$, control variable $u = [u_3 u_4]^T$, and state variable $x = [x_3 x_4]^T$ may be written as follows:

$$y = \begin{bmatrix} i_{\text{bat}} \\ i_{\text{SC}} \end{bmatrix}, u = \begin{bmatrix} d_{\text{bat}} \\ d_{\text{sc}} \end{bmatrix}, x = \begin{bmatrix} i_{\text{bat}} \\ i_{\text{SC}} \end{bmatrix} \quad (29)$$

The control vector variables $u_3, u_4$ are evaluated from Eq. (1) and Eq. (29) as follows:
$$u_1 = 2P_{\text{SC max}} \left[ 1 - \left( 1 - \frac{\dot{y}_1 + \sqrt{\frac{2y_1}{c_{\text{bus}}}} i_{\text{Load}} - P_{\text{Bato}}}{P_{\text{SC max}}} \right)^{1/2} \right] = \vartheta_1(y_1, \dot{y}_1) = P_{\text{SCref}} \quad (16)$$

with $P_{\text{SCmax}}$ is the SC converter's maximum power limit. It can be defined as

$$P_{\text{SC max}} = V_{\text{SC}}^2 / 4r_{\text{SC}} \quad (17)$$

Concerning the second input, the power reference of the battery $P_{\text{Battref}}$ is calculated from Eqs. 5, 6, and 9 as follows:

$$u_2 = 2P_{\text{Bat max}} \left[ 1 - \left( 1 - \frac{\dot{y}_2 + \sqrt{\frac{2y_2}{c_{\text{bus}}}} i_{\text{Load}}}{P_{\text{Bat max}}} \right)^{1/2} \right] = \vartheta_2(y_1, \dot{y}_2) = P_{\text{Bat ref}} \quad (18)$$

where $P_{\text{Battmax}}$ is defined as the maximum limiting power from the battery converter.

$$P_{\text{Bat max}} = V_{\text{bat}}^2 / 4r_{\text{Bat}} \quad (19)$$

Above all, the flat output ($y_1 = E_{\text{bus}}$) is the most crucial variable to control. For this, a classical PI controller is used to ensure that this flat variable is under control. Assuming the SC control loop is significantly faster than the battery control loop$^{44}$, the DC-bus power stated in Eq. (7) can be approximated as:

$$\dot{E}_{\text{Bus}} = P_{\text{Sco}} \quad (20)$$

The transfer function is represented as a straightforward integrator. A PI regulator is used to control this part$^{44}$. Recognizing that $y_1 = E_{\text{bus}}$, as shown below:

$$\dot{y}_1 = \frac{1}{s} \left( k_p^{V_{\text{bus}}} + \frac{k_i^{V_{\text{bus}}}}{s} \right) (y_1 - y_{1-\text{ref}}) \quad (21)$$

where $k_p^{V_{\text{bus}}}, k_i^{V_{\text{bus}}}$ are the integral and the proportional gains chosen so that the closed loop characteristic polynomial is expressed as follows:

$$p(s) = s^2 + \lambda_1 s + \lambda_0 \quad (22)$$

Clearly, the error $e_1 = y_1 - y_{1-\text{ref}}$ meets:

$$\ddot{e}_1 + k_p^{V_{\text{bus}}} \dot{e}_1 + k_i^{V_{\text{bus}}} e_1 = 0 \quad (23)$$

By matching the characteristic polynomial $p(s)$ to a desired one $p_{\text{des}}(s)$, given by Eq. (22), with pre-specified root positions, an adequate choice of controller parameters can be calculated by Eq. (25) and Eq. (26).

$$p_{\text{des}}(s) = s^2 + 2\xi\omega_n s + \omega_n^2 \quad (24)$$

$$k_p^{V_{\text{bus}}} = 2\xi\omega_n \quad (25)$$

$$k_i^{V_{\text{bus}}} = \omega_n^2 \quad (26)$$

where $\omega_n$ represents the natural frequency, and $\xi$ represents the dumping ratio.

The SC energy control loop depends on total energy management. A linearizing feedback control rule is used to achieve an exponential asymptotic tracking of the trajectory$^{44}$ as follows:

$$(\dot{y}_2 - \dot{y}_{2-\text{ref}}) + k_p^{V_{\text{SC}}} (y_2 - y_{2-\text{ref}}) = 0 \quad (27)$$

$$k_p^{V_{\text{SC}}} = \omega_n^2 \quad (28)$$

where $k_p^{V_{\text{SC}}}$ is the proportional gain of the SC voltage controller.

### Current control

The flat output $y = [y_3 y_4]^T$, control variable $u = [u_3 u_4]^T$, and state variable $x = [x_3 x_4]^T$ may be written as follows:

$$y = \begin{bmatrix} i_{\text{bat}} \\ i_{\text{SC}} \end{bmatrix}, u = \begin{bmatrix} d_{\text{bat}} \\ d_{\text{sc}} \end{bmatrix}, x = \begin{bmatrix} i_{\text{bat}} \\ i_{\text{SC}} \end{bmatrix} \quad (29)$$

The control vector variables $u_3, u_4$ are evaluated from Eq. (1) and Eq. (29) as follows:
$$u_3 = \frac{1}{V_{bus}}(L_{bat}\dot{y}_3 - V_{bat} + r_{batt}y_3) = v_1(y_3, \dot{y}_3) = d_{bat} \quad (30)$$

$$u_4 = \frac{1}{V_{bus}}(L_{sc}\dot{y}_4 - V_{sc} + r_{sc}y_4) = v_2(y_4, \dot{y}_4) = d_{sc} \quad (31)$$

The first current control law, $y_{3-ref}$, defines the set-point for the battery current. The closed-loop control law is written as follows$^{37,45}$:

$$(\dot{y}_3 - \dot{y}_{3-ref}) + k_p^{bat}(y_3 - y_{3-ref}) + k_i^{bat} \int (y_3 - y_{3-ref}) dt = 0 \quad (32)$$

where $k_p^{bat}$ and $k_i^{bat}$ represent the controller's parameters. Let us consider the following desired dynamic polynomial$^{46}$:

$$p_1(s) = s^2 + 2\xi_{i1}\omega_{ni1}s + \omega_{ni1}^2 \quad (33)$$

By matching the derivative of Eq. (32) with the desired dynamic polynomial $p_1(s)$ with predetermined root locations, the appropriate controller parameters are expressed by:

$$k_p^{bat} = \xi_{i1}\omega_{ni1} \quad (34)$$

$$k_i^{bat} = \omega_{ni1}^2 \quad (35)$$

The second current control law based on feedback regulation is written as follows:

$$(\dot{y}_4 - \dot{y}_{4-ref}) + k_p^{SC}(y_4 - y_{4-ref}) + k_i^{SC} \int (y_4 - y_{4-ref}) dt = 0 \quad (36)$$

By matching the derivative of Eq. (36) with the following desired characteristic polynomial $p_2(s)$ with pre-specified root positions, the controller parameters can be obtained as in Eq. (38) and Eq. (39).

$$p_2(s) = s^2 + 2\xi_{i2}\omega_{ni2}s + \omega_{ni2}^2 \quad (37)$$

$$k_p^{SC} = 2\xi_{i2}\omega_{ni2} \quad (38)$$

$$k_i^{SC} = \omega_{ni2}^2 \quad (39)$$

where $y_{3-ref}$ and $y_{4-ref}$ are the required inductor current references; $k_i^{sc}$, $k_i^{bat}$ are the integral gains of the supercapacitor and battery current controllers, respectively. $k_p^{sc}$ and $k_p^{bat}$ represent the proportional gains of the supercapacitor and the battery current controllers, respectively. $\omega_{ni1}$, $\omega_{ni2}$ and $\xi_{i1}$, $\xi_{i2}$ represent the natural frequencies and the dumping factors, respectively.

**Salp Swarm Algorithm:**

Slap swarm algorithm (SSA) was proposed by Mirjalili in$^{32}$ and inspired by the behavior of slaps in the ocean. It is distinguished primarily by its high precision and quick convergence capabilities. There are two sorts of agents in the agent set: leaders and followers, and their movements can be described as:

$$LP(n) = \begin{cases} FP(n) + c_1((ub - lb)c_2 + lb) & \text{if } c_3 \ge \frac{1}{2} \\ FP(n) - c_1((ub - lb)c_2 + lb) & \text{if } c_3 < \frac{1}{2} \end{cases} \quad (40)$$

where $LP(n)$ and $FP(n)$ are the leader and food position at iteration $(n)$, respectively, and $c_2$ and $c_3$ are arbitrary variables $[0,1]$. $u_b$ and $l_b$ are the higher and lower search space boundaries, respectively. The coefficient $c_1$ is the most significant parameter that affects the algorithm performance since it balances exploration and exploitation. It is described as:

$$c_1 = 2e^{-(4n/N_{max})^2} \quad (41)$$

where $N_{max}$ is the maximum number of iterations.

The movement of the follower can be expressed as:

$$FP_i(n) = \frac{1}{2}(FP_i(n-1) + FP_{i-1}(n)) \quad (42)$$

where $FP_i(n)$ is the position of the $i$-th follower. This last one updates its position based on its position and the position of the previous salp. The objectives are attenuating the battery current harmonics and DC voltage ripples, reducing the voltage overshoots, and ensuring the stable operation of the HPS. Thus, the cost function can be formulated as the integral square error (ISE) as:
$$u_3 = \frac{1}{V_{bus}}(L_{bat}\dot{y}_3 - V_{bat} + r_{batt}y_3) = v_1(y_3, \dot{y}_3) = d_{bat} \quad (30)$$

$$u_4 = \frac{1}{V_{bus}}(L_{sc}\dot{y}_4 - V_{sc} + r_{sc}y_4) = v_2(y_4, \dot{y}_4) = d_{sc} \quad (31)$$

The first current control law, $y_{3-ref}$, defines the set-point for the battery current. The closed-loop control law is written as follows$^{37,45}$:

$$(\dot{y}_3 - \dot{y}_{3-ref}) + k_p^{bat}(y_3 - y_{3-ref}) + k_i^{bat} \int (y_3 - y_{3-ref}) dt = 0 \quad (32)$$

where $k_p^{bat}$ and $k_i^{bat}$ represent the controller's parameters. Let us consider the following desired dynamic polynomial$^{46}$:

$$p_1(s) = s^2 + 2\xi_{i1}\omega_{ni1}s + \omega_{ni1}^2 \quad (33)$$

By matching the derivative of Eq. (32) with the desired dynamic polynomial $p_1(s)$ with predetermined root locations, the appropriate controller parameters are expressed by:

$$k_p^{bat} = \xi_{i1}\omega_{ni1} \quad (34)$$

$$k_i^{bat} = \omega_{ni1}^2 \quad (35)$$

The second current control law based on feedback regulation is written as follows:

$$(\dot{y}_4 - \dot{y}_{4-ref}) + k_p^{SC}(y_4 - y_{4-ref}) + k_i^{SC} \int (y_4 - y_{4-ref}) dt = 0 \quad (36)$$

By matching the derivative of Eq. (36) with the following desired characteristic polynomial $p_2(s)$ with pre-specified root positions, the controller parameters can be obtained as in Eq. (38) and Eq. (39).

$$p_2(s) = s^2 + 2\xi_{i2}\omega_{ni2}s + \omega_{ni2}^2 \quad (37)$$

$$k_p^{SC} = 2\xi_{i2}\omega_{ni2} \quad (38)$$

$$k_i^{SC} = \omega_{ni2}^2 \quad (39)$$

where $y_{3-ref}$ and $y_{4-ref}$ are the required inductor current references; $k_i^{sc}$, $k_i^{bat}$ are the integral gains of the supercapacitor and battery current controllers, respectively. $k_p^{sc}$ and $k_p^{bat}$ represent the proportional gains of the supercapacitor and the battery current controllers, respectively. $\omega_{ni1}$, $\omega_{ni2}$ and $\xi_{i1}$, $\xi_{i2}$ represent the natural frequencies and the dumping factors, respectively.

**Salp Swarm Algorithm:**

Slap swarm algorithm (SSA) was proposed by Mirjalili in$^{32}$ and inspired by the behavior of slaps in the ocean. It is distinguished primarily by its high precision and quick convergence capabilities. There are two sorts of agents in the agent set: leaders and followers, and their movements can be described as:

$$LP(n) = \begin{cases} FP(n) + c_1((ub - lb)c_2 + lb) & \text{if } c_3 \ge \frac{1}{2} \\ FP(n) - c_1((ub - lb)c_2 + lb) & \text{if } c_3 < \frac{1}{2} \end{cases} \quad (40)$$

where $LP(n)$ and $FP(n)$ are the leader and food position at iteration $(n)$, respectively, and $c_2$ and $c_3$ are arbitrary variables $[0,1]$. $u_b$ and $l_b$ are the higher and lower search space boundaries, respectively. The coefficient $c_1$ is the most significant parameter that affects the algorithm performance since it balances exploration and exploitation. It is described as:

$$c_1 = 2e^{-(4n/N_{max})^2} \quad (41)$$

where $N_{max}$ is the maximum number of iterations.

The movement of the follower can be expressed as:

$$FP_i(n) = \frac{1}{2}(FP_i(n-1) + FP_{i-1}(n)) \quad (42)$$

where $FP_i(n)$ is the position of the $i$-th follower. This last one updates its position based on its position and the position of the previous salp. The objectives are attenuating the battery current harmonics and DC voltage ripples, reducing the voltage overshoots, and ensuring the stable operation of the HPS. Thus, the cost function can be formulated as the integral square error (ISE) as:
<figure><img src="image_10.png" /><figcaption>Figure 7. PIL co-simulation strategy: (a) PIL scheme, (b) PIL platform.</figcaption></figure>

$$f = \min \left( \int_{0}^{t} \sqrt{\varepsilon dt} \right) \qquad (43)$$

For the hybrid power systems (HPSs), there are four errors written as follows:

$$\left\{ \begin{array}{l} \varepsilon_{V_{bus}} = V_{bus}^{ref} - V_{bus} \\ \varepsilon_{V_{SC}} = V_{SC}^{ref} - V_{SC} \\ \varepsilon_{ibatt} = i_{bat}^{ref} - i_{batt} \\ \varepsilon_{i_{SC}} = (i_{SC}^{ref} - i_{h}) - i_{SC} \end{array} \right. \qquad (44)$$

where $\varepsilon_{V_{bus}}$, $\varepsilon_{SC}$ are the DC-bus and the supercapacitor voltage errors, respectively, $\varepsilon_{ibatt}$ and $\varepsilon_{i_{SC}}$ represent the battery and the supercapacitor current errors, $i_{SC}^{ref}$ and $i_{SC}$ represent the SC's reference and measured currents; $i_{bat}^{ref}$ and $i_{batt}$ are the battery's reference and measured currents, respectively; $V_{bus}^{ref}$ and $V_{bus}$ denotes the DC bus's reference and measured voltages, respectively; $V_{SC}^{ref}$ and $V_{SC}$ are the SC's reference and measured voltages; and $i_{h}$ is the harmonic current.

Relying on the adopted cost function representing the ISE, the SSA calculates the necessary controllers' parameters $k_{p}^{V_{bus}}, k_{i}^{V_{bus}}, k_{p}^{V_{SC}}, k_{i}^{V_{SC}}, k_{p}^{SC}, k_{i}^{SC}, k_{p}^{bat}, k_{p}^{bat}$.

## PIL implementation technique description and performing steps

The PIL co-simulation technique allows the verification and validation of the proposed control algorithms by generating code onto the embedded processor core and running these algorithms in a real environment based on the C2000 launchxl-f28379d DSP board. During PIL co-simulation, the implemented control algorithm is linked to a computer on which the physical system model is carried out. Subsequently, it is possible to evaluate the performance of the system in order to assess and improve some essential factors such as storage capacity,
<figure><img src="image_10.png" /><figcaption>Figure 7. PIL co-simulation strategy: (a) PIL scheme, (b) PIL platform.</figcaption></figure>

$$f = \min \left( \int_{0}^{t} \sqrt{\varepsilon dt} \right) \qquad (43)$$

For the hybrid power systems (HPSs), there are four errors written as follows:

$$\left\{ \begin{aligned} \varepsilon_{V_{bus}} &= V_{bus}^{ref} - V_{bus} \\ \varepsilon_{V_{SC}} &= V_{SC}^{ref} - V_{SC} \\ \varepsilon_{ibatt} &= i_{bat}^{ref} - i_{batt} \\ \varepsilon_{i_{SC}} &= (i_{SC}^{ref} - i_{h}) - i_{SC} \end{aligned} \right. \qquad (44)$$

where $\varepsilon_{V_{bus}}$, $\varepsilon_{SC}$ are the DC-bus and the supercapacitor voltage errors, respectively, $\varepsilon_{ibatt}$ and $\varepsilon_{iSC}$ represent the battery and the supercapacitor current errors, $i_{SC}^{ref}$ and $i_{SC}$ represent the SC's reference and measured currents; $i_{bat}^{ref}$ and $i_{batt}$ are the battery's reference and measured currents, respectively; $V_{bus}^{ref}$ and $V_{bus}$ denotes the DC bus's reference and measured voltages, respectively; $V_{SC}^{ref}$ and $V_{SC}$ are the SC's reference and measured voltages; and $i_{h}$ is the harmonic current.

Relying on the adopted cost function representing the ISE, the SSA calculates the necessary controllers' parameters $k_{p}^{V_{bus}}, k_{i}^{V_{bus}}, k_{p}^{V_{SC}}, k_{i}^{V_{SC}}, k_{p}^{SC}, k_{i}^{SC}, k_{p}^{bat}, k_{p}^{bat}$.

## PIL implementation technique description and performing steps

The PIL co-simulation technique allows the verification and validation of the proposed control algorithms by generating code onto the embedded processor core and running these algorithms in a real environment based on the C2000 launchxl-f28379d DSP board. During PIL co-simulation, the implemented control algorithm is linked to a computer on which the physical system model is carried out. Subsequently, it is possible to evaluate the performance of the system in order to assess and improve some essential factors such as storage capacity,
code size, and execution of the algorithm according to the required time. As indicated in Fig. 7, during the prototyping of the PIL, based on a fixed simulation time, the power part of the power system is simulated in the Matlab/Simulink platform. At each step, the C2000 launchxl-f28379d DSP board receives the signals from the computer, implements control algorithms, and sends the control commands back to the computer to control the power system. At this point, a PIL co-simulation cycle is performed. The data exchange between the computer and the DSP board is synchronized using the serial communication of the DSP board. To perform the PIL, the following steps are needed to be carried out:

* Connecting the C2000 launchxl-f28379d DSP board to the computer,

* Tuning several settings and configuration parameters from the select hardware implementation tab found on Matlab/Simulink,

* Selecting the appropriate hardware port using the defined DSP board under the Hardware implementation,

* Selecting target hardware resources and selecting the correct device name of the adopted DSP board again,

* Choosing the external mode to set up the serial communication.

After that, a configuration of the PIL procedure should be performed based on multiple line codes in the command window of the Matlab software as follows:

* Calling of the system model,

* Setting the com port number,

* Defining the board rates, which represent how fast the computer and the DSP board will communicate,

* Enabling the serial communication for the PIL co-simulation,

* Generating the PIL model that will be used for the PIL procedure.

A trial implementation in a real-time interface (RTI) utilizing a DSpace card can be conducted to enhance the proposed energy management in the electric vehicle (EV) system. This experimental evaluation aims to assess the performance of the enhanced energy management system.

The experimental implementation in a real-time interface (RTI) using a DSpace card for the proposed energy management optimization in an electric car system is a promising avenue for advancing research in sustainable transportation. By leveraging the computational capabilities of a DSpace card within the vehicle's infrastructure, the energy management system can dynamically optimize power distribution among various components, such as the battery, motor, and auxiliary systems, in real-time.

## Results and discussion

### Simulation part results

In this section, MATLAB-Simulink is utilized to construct the HPS model and validate the suggested EMS. The HPS is simulated employing the urban drive cycle ECE-15, where the HPS and SSA parameters are summarized in Table 1,2, respectively.

Notwithstanding the change in the driving cycle, the EV speed response is shown in Fig. 8.a, showing a good follow-up. The torque curve in Fig. 8.b, shows that the motor creates max torque when the vehicle's velocity attains the reference path. When it reaches a steady state, the motor provides lower torque, only enough to compensate for total load torque.

The suggested energy management technique can stabilize the DC-bus voltage immediately, as depicted in Fig. 9, despite the significant load power variations. Compared to the classical differential flatness and the PSO differential flatness, the proposed EMS based on the SSA differential flatness reduces DC-bus voltage ripples and overshoots. Indeed, for a maximum load step of 21 kW, the DC-bus voltage overshoot is minimized by 15V (3.2%) and 3V (0.6%) compared to the classical differential flatness and PSO differential flatness.

<table><caption>**Table 2.** HPS parameters.</caption><thead><tr><td>Parameters</td><td>Value</td></tr></thead><tbody><tr><td>r<sub>batt</sub> Ω</td><td>0.1</td></tr><tr><td>r<sub>sc</sub> Ω</td><td>0.1</td></tr><tr><td>L<sub>bat</sub> mH</td><td>2</td></tr><tr><td>L<sub>sc</sub> mH</td><td>2</td></tr><tr><td>V<sup>ref</sup><sub>bus</sub> V</td><td>500</td></tr><tr><td>V<sup>ref</sup><sub>SC</sub> V</td><td>200</td></tr><tr><td>V<sup>ref</sup><sub>bat</sub> V</td><td>100</td></tr><tr><td>c<sub>SC</sub> F</td><td>80</td></tr><tr><td>c<sub>bat</sub> Ah</td><td>340</td></tr><tr><td>c<sub>bus</sub> μF</td><td>2000</td></tr></tbody></table>
code size, and execution of the algorithm according to the required time. As indicated in Fig. 7, during the prototyping of the PIL, based on a fixed simulation time, the power part of the power system is simulated in the Matlab/Simulink platform. At each step, the C2000 launchxl-f28379d DSP board receives the signals from the computer, implements control algorithms, and sends the control commands back to the computer to control the power system. At this point, a PIL co-simulation cycle is performed. The data exchange between the computer and the DSP board is synchronized using the serial communication of the DSP board. To perform the PIL, the following steps are needed to be carried out:

* Connecting the C2000 launchxl-f28379d DSP board to the computer,

* Tuning several settings and configuration parameters from the select hardware implementation tab found on Matlab/Simulink,

* Selecting the appropriate hardware port using the defined DSP board under the Hardware implementation,

* Selecting target hardware resources and selecting the correct device name of the adopted DSP board again,

* Choosing the external mode to set up the serial communication.

After that, a configuration of the PIL procedure should be performed based on multiple line codes in the command window of the Matlab software as follows:

* Calling of the system model,

* Setting the com port number,

* Defining the board rates, which represent how fast the computer and the DSP board will communicate,

* Enabling the serial communication for the PIL co-simulation,

* Generating the PIL model that will be used for the PIL procedure.

A trial implementation in a real-time interface (RTI) utilizing a DSpace card can be conducted to enhance the proposed energy management in the electric vehicle (EV) system. This experimental evaluation aims to assess the performance of the enhanced energy management system.

The experimental implementation in a real-time interface (RTI) using a DSpace card for the proposed energy management optimization in an electric car system is a promising avenue for advancing research in sustainable transportation. By leveraging the computational capabilities of a DSpace card within the vehicle's infrastructure, the energy management system can dynamically optimize power distribution among various components, such as the battery, motor, and auxiliary systems, in real-time.

## Results and discussion
### Simulation part results

In this section, MATLAB-Simulink is utilized to construct the HPS model and validate the suggested EMS. The HPS is simulated employing the urban drive cycle ECE-15, where the HPS and SSA parameters are summarized in Table 1, respectively.

Notwithstanding the change in the driving cycle, the EV speed response is shown in Fig. 8.a, showing a good follow-up. The torque curve in Fig. 8.b, shows that the motor creates max torque when the vehicle's velocity attains the reference path. When it reaches a steady state, the motor provides lower torque, only enough to compensate for total load torque.

The suggested energy management technique can stabilize the DC-bus voltage immediately, as depicted in Fig. 9, despite the significant load power variations. Compared to the classical differential flatness and the PSO differential flatness, the proposed EMS based on the SSA differential flatness reduces DC-bus voltage ripples and overshoots. Indeed, for a maximum load step of 21 kW, the DC-bus voltage overshoot is minimized by 15V (3.2%) and 3V (0.6%) compared to the classical differential flatness and PSO differential flatness.

<table><caption>**Table 2.** HPS parameters.</caption><thead><tr><td>Parameters</td><td>Value</td></tr></thead><tbody><tr><td>r<sub>batt</sub> Ω</td><td>0.1</td></tr><tr><td>r<sub>sc</sub> Ω</td><td>0.1</td></tr><tr><td>L<sub>bat</sub> mH</td><td>2</td></tr><tr><td>L<sub>sc</sub> mH</td><td>2</td></tr><tr><td>V<sup>ref</sup><sub>bus</sub> V</td><td>500</td></tr><tr><td>V<sup>ref</sup><sub>SC</sub> V</td><td>200</td></tr><tr><td>V<sup>ref</sup><sub>bat</sub> V</td><td>100</td></tr><tr><td>c<sub>SC</sub> F</td><td>80</td></tr><tr><td>c<sub>bat</sub> Ah</td><td>340</td></tr><tr><td>c<sub>bus</sub> μF</td><td>2000</td></tr></tbody></table>
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 8.** EV traction side simulation results: (a) EV linear speed, (b) EV load torque (TL) and SynRM electromagnetic torque (Te).</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 9.** DC-bus voltage.</figcaption></figure>
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 8.** EV traction side simulation results: (a) EV linear speed, (b) EV load torque (TL) and SynRM electromagnetic torque (Te).</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 9.** DC-bus voltage.</figcaption></figure>
The battery provides most of the load power to the motor during the acceleration phases, which causes the battery SOC to decrease. It receives energy during deceleration phases when the motor torque is negative, thus increasing the SOC battery, as shown in Fig. 10.

From Fig. 11, which depicts the battery current, the proposed EMS can improve the power quality by reducing the current ripple, which is the main cause of battery aging, thus enhancing the battery life cycle.

Table 3 provides the comparison study results between the proposed EMS based on optimal adaptive and classical differential flatness strategies. In this table, the battery current ripple differences $ΔI_1 = ΔI_{C\_FLAT} - ΔI_{SSA\_FLAT}$ and $ΔI_2 = ΔI_{PSO\_FLAT} - ΔI_{SSA\_FLAT}$ are evaluated, where $ΔI_{C\_FLAT}$ and $ΔI_{SSA\_FLAT}$, $ΔI_{PSO\_FLAT}$ are the battery current ripple at the periods ($P_n$) for classical FLAT and SSA\_FLAT, PSO\_FLAT, respectively. The battery current ripple is minimized by 15.5A for a maximum load step of 21kW.

Figure 12, shows that the SSA differential flatness minimizes the battery current THD to 10.49% instead of 77.39% for the classical differential flatness strategy and 34.52 for the PSO differential flatness strategy. This finding suits the purpose of the proposed EMS, which is the reduction of battery current harmonics leading to battery lifecycle enhancement.

As shown in Fig. 13, the battery supplies the motor power and absorbs it during braking periods, while the SC works to assist the battery during transient periods (acceleration and deceleration phases), which accords with the management strategy.

Figure 14, shows that the suggested EMS algorithm exhibits quicker supercapacitor voltage dynamics with fewer ripples than the classical differential flatness technique. The quick response of the proposed management system leads to enhancing the stability of the SC voltage.

<figure><img src="image_14.png" /><figcaption>Figure 10. Battery SoC (%).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure 11. Battery Current.</figcaption></figure>

<table><caption>Table 3. Optimizer parameters.</caption><thead><tr><th>gains</th><th>k<sub>p</sub><sup>Vbus</sup></th><th>k<sub>i</sub><sup>Vbus</sup></th><th>k<sub>p</sub><sup>VSC</sup></th><th>k<sub>i</sub><sup>VSC</sup></th><th>k<sub>p</sub><sup>SC</sup></th><th>k<sub>i</sub><sup>SC</sup></th><th>k<sub>p</sub><sup>bat</sup></th><th>k<sub>i</sub><sup>bat</sup></th></tr></thead><tbody><tr><td>Lower bound (× 10)</td><td>100<sup>^2</sup></td><td>418.879</td><td>0.01</td><td>8.7730 × 10<sup>4</sup></td><td>8.3576</td><td>9.7478 × 10<sup>3</sup></td><td>27.8253</td><td></td></tr><tr><td>Upper bound (× 0.5)</td><td>100<sup>^2</sup></td><td>418.879</td><td>0.01</td><td>8.7730 × 10<sup>4</sup></td><td>8.3576</td><td>9.7478 × 10<sup>3</sup></td><td>27.8253</td><td></td></tr></tbody><tfoot><tr><td colspan="9">Search Agents number = 30 Maximum number of iteration = 150</td></tr></tfoot></table>
The battery provides most of the load power to the motor during the acceleration phases, which causes the battery SOC to decrease. It receives energy during deceleration phases when the motor torque is negative, thus increasing the SOC battery, as shown in Fig. 10.

From Fig. 11, which depicts the battery current, the proposed EMS can improve the power quality by reducing the current ripple, which is the main cause of battery aging, thus enhancing the battery life cycle.

Table 3 provides the comparison study results between the proposed EMS based on optimal adaptive and classical differential flatness strategies. In this table, the battery current ripple differences $ΔI_1 = ΔI_{C\_FLAT} - ΔI_{SSA\_FLAT}$ and $ΔI_2 = ΔI_{PSO\_FLAT} - ΔI_{SSA\_FLAT}$ are evaluated, where $ΔI_{C\_FLAT}$ and $ΔI_{SSA\_FLAT}$, $ΔI_{PSO\_FLAT}$ are the battery current ripple at the periods ($P_n$) for classical FLAT and SSA\_FLAT, PSO\_FLAT, respectively. The battery current ripple is minimized by 15.5A for a maximum load step of 21kW.

Figure 12, shows that the SSA differential flatness minimizes the battery current THD to 10.49% instead of 77.39% for the classical differential flatness strategy and 34.52 for the PSO differential flatness strategy. This finding suits the purpose of the proposed EMS, which is the reduction of battery current harmonics leading to battery lifecycle enhancement.

As shown in Fig. 13, the battery supplies the motor power and absorbs it during braking periods, while the SC works to assist the battery during transient periods (acceleration and deceleration phases), which accords with the management strategy.

Figure 14, shows that the suggested EMS algorithm exhibits quicker supercapacitor voltage dynamics with fewer ripples than the classical differential flatness technique. The quick response of the proposed management system leads to enhancing the stability of the SC voltage.

<figure><img src="image_14.png" /><figcaption>Figure 10. Battery SoC (%).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure 11. Battery Current.</figcaption></figure>

<table><caption>Table 3. Optimizer parameters.</caption><thead><tr><th>gains</th><th>k<sub>p</sub><sup>Vbus</sup></th><th>k<sub>i</sub><sup>Vbus</sup></th><th>k<sub>p</sub><sup>VSC</sup></th><th>k<sub>i</sub><sup>VSC</sup></th><th>k<sub>p</sub><sup>sc</sup></th><th>k<sub>i</sub><sup>sc</sup></th><th>k<sub>p</sub><sup>bat</sup></th><th>k<sub>i</sub><sup>bat</sup></th></tr></thead><tbody><tr><td>Lower bound (× 10)</td><td>100<sup>^2</sup></td><td>418.879</td><td>0.01</td><td>8.7730 × 10<sup>4</sup></td><td>8.3576</td><td>9.7478 × 10<sup>3</sup></td><td>27.8253</td><td></td></tr><tr><td>Upper bound (× 0.5)</td><td>100<sup>^2</sup></td><td>418.879</td><td>0.01</td><td>8.7730 × 10<sup>4</sup></td><td>8.3576</td><td>9.7478 × 10<sup>3</sup></td><td>27.8253</td><td></td></tr></tbody><tfoot><tr><td colspan="9">Search Agents number = 30 Maximum number of iteration = 150</td></tr></tfoot></table>
<figure><img src="image_16.png" /><figcaption>**Figure 12.** Harmonics spectrum of battery current: (a) for SSA differential flatness strategy; (b) for PSO differential flatness strategy; (c) for classical differential flatness strategy.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 13.** Load, battery, and SC power curves</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 12.** Harmonics spectrum of battery current: (a) for SSA differential flatness strategy; (b) for PSO differential flatness strategy; (c) for classical differential flatness strategy.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 13.** Load, battery, and SC power curves</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>Figure 14. Supercapacitor voltage.</figcaption></figure>

## Co-simulation results

To approve and evaluate the performance of the suggested EMS, the system was modeled using embedded Matlab functions and co-simulated using the C2000 launchxl-f28379d DSP board through the processor-in-the-loop (PIL). The co-simulation is performed utilizing a reduced-time version of the ECE-15 urban drive cycle depicted in Fig. 15, based on the parameters listed in Table 4, Tables A1 and A4.

The proposed EMS stabilizes the bus voltage shown in Fig. 16, against load variations. Compared to classical differential FLAT, using the SSA differential flatness, the maximum overshoot voltage is 4 V instead of 10 for the PSO differential flatness and 15 V for the classical FLAT strategy at the maximum power demanded by the electric vehicle (about 21 kW). For the classical FLAT strategy, the SSA differential flatness minimizes the voltage ripple at ($\Delta V = 0.5$ V) instead of ($\Delta V = 4$ V).

From Fig. 17, representing the current of the battery, the ripple is minimized by the SSA differential flatness to ($\Delta I = 3$ A) instead of ($\Delta I = 12$ A) for PSO differential flatness and ($\Delta I = 20$ A) for the classical differential flatness strategy at the maximum power demanded by the electric vehicle (about 21 kW).

As illustrated in Fig. 18, the battery provides the average power demanded by the traction system and receives energy during the braking phases. The supercapacitor assists the battery during transient periods (acceleration and deceleration phases), which agrees with the adopted energy management strategy. These results confirm the efficiency of the suggested EMS in managing both energy storage systems.

<figure><img src="image_19.png" /><figcaption>Figure 15. Reduced time version of the ECE-15 urban drive cycle</figcaption></figure>

<table><caption>Table 4. Comparison based on battery current ripple between classical and proposed differential flatness.</caption><thead><tr><th>Time (P<sub>n</sub>) (s)</th><th>P<sub>1</sub> = [0,11]</th><th>P<sub>2</sub> = [28,49]</th><th>P<sub>3</sub> = [61,85]</th><th>P<sub>4</sub> = [143,155]</th></tr></thead><tbody><tr><td>ΔI<sub>C_FLAT</sub> (A)</td><td>20</td><td>15</td><td>18</td><td>20</td></tr><tr><td>ΔI<sub>PSO_FLAT</sub> (A)</td><td>10</td><td>6</td><td>8.7</td><td>15</td></tr><tr><td>ΔI<sub>SSA_FLAT</sub> (A)</td><td>3</td><td>2.8</td><td>2.5</td><td>2.85</td></tr><tr><td>ΔI<sub>1</sub></td><td>17</td><td>12.2</td><td>15.5</td><td>17.15</td></tr><tr><td>ΔI<sub>2</sub></td><td>10</td><td>9</td><td>9.3</td><td>5</td></tr></tbody></table>
<figure><img src="image_18.png" /><figcaption>Figure 14. Supercapacitor voltage.</figcaption></figure>

## Co-simulation results

To approve and evaluate the performance of the suggested EMS, the system was modeled using embedded Matlab functions and co-simulated using the C2000 launchxl-f28379d DSP board through the processor-in-the-loop (PIL). The co-simulation is performed utilizing a reduced-time version of the ECE-15 urban drive cycle depicted in Fig. 15, based on the parameters listed in Table 4, Tables A1 and A4.

The proposed EMS stabilizes the bus voltage shown in Fig. 16, against load variations. Compared to classical differential FLAT, using the SSA differential flatness, the maximum overshoot voltage is 4 V instead of 10 for the PSO differential flatness and 15 V for the classical FLAT strategy at the maximum power demanded by the electric vehicle (about 21 kW). For the classical FLAT strategy, the SSA differential flatness minimizes the voltage ripple at (ΔV = 0.5 V) instead of (ΔV = 4 V).

From Fig. 17, representing the current of the battery, the ripple is minimized by the SSA differential flatness to (ΔI = 3 A) instead of (ΔI = 12 A) for PSO differential flatness and (ΔI = 20 A) for the classical differential flatness strategy at the maximum power demanded by the electric vehicle (about 21 kW).

As illustrated in Fig. 18, the battery provides the average power demanded by the traction system and receives energy during the braking phases. The supercapacitor assists the battery during transient periods (acceleration and deceleration phases), which agrees with the adopted energy management strategy. These results confirm the efficiency of the suggested EMS in managing both energy storage systems.

<figure><img src="image_19.png" /><figcaption>Figure 15. Reduced time version of the ECE-15 urban drive cycle</figcaption></figure>

<table><caption>Table 4. Comparison based on battery current ripple between classical and proposed differential flatness.</caption><thead><tr><th>Time (P<sub>n</sub>) (s)</th><th>P<sub>1</sub> = [0,11]</th><th>P<sub>2</sub> = [28,49]</th><th>P<sub>3</sub> = [61,85]</th><th>P<sub>4</sub> = [143,155]</th></tr></thead><tbody><tr><td>ΔI<sub>C_PLAT</sub> (A)</td><td>20</td><td>15</td><td>18</td><td>20</td></tr><tr><td>ΔI<sub>PSO_PLAT</sub> (A)</td><td>10</td><td>6</td><td>8.7</td><td>15</td></tr><tr><td>ΔI<sub>SSA_PLAT</sub> (A)</td><td>3</td><td>2.8</td><td>2.5</td><td>2.85</td></tr><tr><td>ΔI<sub>1</sub></td><td>17</td><td>12.2</td><td>15.5</td><td>17.15</td></tr><tr><td>ΔI<sub>2</sub></td><td>10</td><td>9</td><td>9.3</td><td>5</td></tr></tbody></table>
<figure><img src="image_20.png" /><figcaption>**Figure 16.** DC-bus voltage.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 17.** Battery Current</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure 18.** Load, battery, and supercapacitor power waveforms</figcaption></figure>

Compared to the PSO and classical differential flatness strategy, the proposed EMS algorithm shows faster supercapacitor voltage dynamics with fewer ripples, as depicted in Fig. 19. During the acceleration and deceleration phases, the SC provides the electric vehicle power until the power provided by the battery reaches its reference, which improves the DC-bus voltage stability.

The suggested technique adjusts its control parameters in accordance with the measured cost function, resulting in enhanced system performance for battery current harmonics, DC-bus voltage ripples, and voltage overshoots.

## Conclusion and future works

This work proposed a new optimized energy management strategy for a battery/supercapacitor-based hybrid power system dedicated to an electric vehicle. Based on the optimal difference flatness, the proposed management technique seeks to manage both sources' power well depending on the load demand. The main goal of this energy
<figure><img src="image_20.png" /><figcaption>**Figure 16.** DC-bus voltage.</figcaption></figure>

<figure><img src="image_21.png" /><figcaption>**Figure 17.** Battery Current</figcaption></figure>

<figure><img src="image_22.png" /><figcaption>**Figure 18.** Load, battery, and supercapacitor power waveforms</figcaption></figure>

Compared to the PSO and classical differential flatness strategy, the proposed EMS algorithm shows faster supercapacitor voltage dynamics with fewer ripples, as depicted in Fig. 19. During the acceleration and deceleration phases, the SC provides the electric vehicle power until the power provided by the battery reaches its reference, which improves the DC-bus voltage stability.

The suggested technique adjusts its control parameters in accordance with the measured cost function, resulting in enhanced system performance for battery current harmonics, DC-bus voltage ripples, and voltage overshoots.

## Conclusion and future works

This work proposed a new optimized energy management strategy for a battery/supercapacitor-based hybrid power system dedicated to an electric vehicle. Based on the optimal difference flatness, the proposed management technique seeks to manage both sources' power well depending on the load demand. The main goal of this energy
<figure><img src="image_23.png" /><figcaption>Figure 19. Supercapacitor voltage.</figcaption></figure>

management is to optimize the power quality by reducing the current harmonics while satisfying the SynRM motor power demand, which positively impacts the battery lifecycle. Compared to the classical differential flatness strategy, the proposed optimal adaptive differential flatness strategy can protect the battery against peak current during the acceleration and deceleration phases and significantly reduce the battery current harmonics and DC-bus voltage ripples (Δv=5 V), as well as the voltage overshoots 15 V (3.2%) for a load power of 21 kW. In addition, the online updating approach also enhances the power system behavior under unknown load changes, improving its robustness and efficiency. The obtained co-simulation results using the C2000 launchxl-f28379d DSP board confirm the effectiveness of the suggested energy management strategy.

The proposed EMS holds the potential for broader applicability across various power systems, including those reliant on fuel-cell technology within hybrid electric vehicles. Our ongoing research involves exploring alternative objective functions within the proposed EMS framework and investigating the efficacy of diverse optimization algorithms. These aspects will be further elucidated in our forthcoming publications.

Received: 10 November 2023; Accepted: 31 January 2024

Published online: 13 February 2024

## References

1. Yao, Z., Wang, Y., Liu, B., Zhao, B. & Jiang, Y. Fuel consumption and transportation emissions evaluation of mixed traffic flow with connected automated vehicles and human-driven vehicles on expressway. *Energy* **230**, 120766 (2021).

2. Heinrichs, H., Jochem, P. & Fichtner, W. Including road transport in the EU ETS (European Emissions Trading System): A model-based analysis of the German electricity and transport sector. *Energy* **69**, 708–720 (2014).

3. Greenhouse gas emissions from transport in Europe—European Environment Agency. https://www.eea.europa.eu/data-and-maps/indicators/transport-emissions-of-greenhouse-gases/transport-emissions-of-greenhouse-gases-12.

4. Average CO₂ emissions from new cars and vans registered in Europe increased in 2018, requiring significant emission reductions to meet the 2020 targets. https://climate.ec.europa.eu/news-your-voice/news/average-co2-emissions-new-cars-and-vans-registered-europe-increased-2018-requiring-significant-2020-06-03_en.

5. EU proposes effective ban for new fossil-fuel cars from 2035 - Times of India. https://timesofindia.indiatimes.com/auto/news/eu-proposes-effective-ban-for-new-fossil-fuel-cars-from-2035/articleshow/84433196.cms.

6. Volvo says it will be 'fully electric' by 2030, move car sales online. https://www.cnbc.com/2021/03/02/volvo-says-it-will-be-fully-electric-by-2030-move-car-sales-online.html.

7. CES 2021: GM says auto sector at 'inflection point' toward zero-emission future | S&P Global Market Intelligence. https://www.spglobal.com/marketintelligence/en/news-insights/latest-news-headlines/ces-2021-gm-says-auto-sector-at-inflection-point-toward-zero-emission-future-62081278.

8. How Volkswagen's Sins Fueled Its Redemption - The New York Times. https://www.nytimes.com/2021/03/19/business/volkswagen-electric-cars-share-price.html.

9. Araria, R., Negadi, K., Boudiaf, M. & Marignetti, F. Nonlinear control of DC-DC converters for battery power management in electric vehicle application 82–88. https://doi.org/10.15199/48.2020.03.20. (2020).

10. Soares, R., Bessman, A., Wallmark, O., Lindbergh, G. & Svens, P. Measurements and analysis of battery harmonic currents in a commercial hybrid vehicle. In *2017 IEEE Transportation Electrification Conference Expo, ITEC 2017* 45–50. https://doi.org/10.1109/ITEC.2017.7993245. (2017).

11. Brand, M. J., Hofmann, M. H., Schuster, S. S., Keil, P. & Jossen, A. The influence of current ripples on the lifetime of lithium-ion batteries. *IEEE Trans. Veh. Technol.* **67**, 10438–10445 (2018).

12. Bessman, A., Soares, R., Wallmark, O., Svens, P. & Lindbergh, G. Aging effects of AC harmonics on lithium-ion cells. *J. Energy Storage* **21**, 741–749 (2019).

13. Uddin, K., Moore, A. D., Barai, A. & Marco, J. The effects of high frequency current ripple on electric vehicle battery performance. *Appl. Energy* **178**, 142–154 (2016).

14. Juang, L. W. et al. Investigation of the influence of superimposed AC current on lithium-ion battery aging using statistical design of experiments. *J. Energy Storage* **11**, 93–103 (2017).

15. Uno, M. & Tanaka, K. Influence of high-frequency charge-discharge cycling induced by cell voltage equalizers on the life performance of lithium-ion cells. *IEEE Trans. Veh. Technol.* **60**, 1505–1515 (2011).

16. Lü, X. et al. Energy management of hybrid electric vehicles: A review of energy optimization of fuel cell hybrid power system based on genetic algorithm. *Energy Convers. Manag.* **205**, 112474 (2020).

17. Baba, M. A., Labbadi, M., Cherkaoui, M. & Maaroufi, M. Fuel cell electric vehicles: A review of current power electronic converters topologies and technical challenges. In *IOP Conference Series Earth Environmental Science* **785**, 1–29 (2021).
<figure><img src="image_23.png" /><figcaption>**Figure 19.** Supercapacitor voltage.</figcaption></figure>

management is to optimize the power quality by reducing the current harmonics while satisfying the SynRM motor power demand, which positively impacts the battery lifecycle. Compared to the classical differential flatness strategy, the proposed optimal adaptive differential flatness strategy can protect the battery against peak current during the acceleration and deceleration phases and significantly reduce the battery current harmonics and DC-bus voltage ripples (Δv = 5 V), as well as the voltage overshoots 15 V (3.2%) for a load power of 21 kW. In addition, the online updating approach also enhances the power system behavior under unknown load changes, improving its robustness and efficiency. The obtained co-simulation results using the C2000 launchxl-f28379d DSP board confirm the effectiveness of the suggested energy management strategy.

The proposed EMS holds the potential for broader applicability across various power systems, including those reliant on fuel-cell technology within hybrid electric vehicles. Our ongoing research involves exploring alternative objective functions within the proposed EMS framework and investigating the efficacy of diverse optimization algorithms. These aspects will be further elucidated in our forthcoming publications.

Received: 10 November 2023; Accepted: 31 January 2024

Published online: 13 February 2024

## References

1. Yao, Z., Wang, Y., Liu, B., Zhao, B. & Jiang, Y. Fuel consumption and transportation emissions evaluation of mixed traffic flow with connected automated vehicles and human-driven vehicles on expressway. *Energy* **230**, 120766 (2021).

2. Heinrichs, H., Jochem, P. & Fichtner, W. Including road transport in the EU ETS (European Emissions Trading System): A model-based analysis of the German electricity and transport sector. *Energy* **69**, 708–720 (2014).

3. Greenhouse gas emissions from transport in Europe—European Environment Agency. https://www.eea.europa.eu/data-and-maps/indicators/transport-emissions-of-greenhouse-gases/transport-emissions-of-greenhouse-gases-12.

4. Average CO₂ emissions from new cars and vans registered in Europe increased in 2018, requiring significant emission reductions to meet the 2020 targets. https://climate.ec.europa.eu/news-your-voice/news/average-co2-emissions-new-cars-and-vans-registered-europe-increased-2018-requiring-significant-2020-06-03_en.

5. EU proposes effective ban for new fossil-fuel cars from 2035 - Times of India. https://timesofindia.indiatimes.com/auto/news/eu-proposes-effective-ban-for-new-fossil-fuel-cars-from-2035/articleshow/84433196.cms.

6. Volvo says it will be 'fully electric' by 2030, move car sales online. https://www.cnbc.com/2021/03/02/volvo-says-it-will-be-fully-electric-by-2030-move-car-sales-online.html.

7. CES 2021: GM says auto sector at 'inflection point' toward zero-emission future | S&P Global Market Intelligence. https://www.spglobal.com/marketintelligence/en/news-insights/latest-news-headlines/ces-2021-gm-says-auto-sector-at-inflection-point-toward-zero-emission-future-62081278.

8. How Volkswagen's Sins Fueled Its Redemption - The New York Times. https://www.nytimes.com/2021/03/19/business/volkswagen-electric-cars-share-price.html.

9. Araria, R., Negadi, K., Boudiaf, M. & Marignetti, F. Nonlinear control of DC-DC converters for battery power management in electric vehicle application 82–88. https://doi.org/10.15199/48.2020.03.20. (2020).

10. Soares, R., Bessman, A., Wallmark, O., Lindbergh, G. & Svens, P. Measurements and analysis of battery harmonic currents in a commercial hybrid vehicle. In *2017 IEEE Transportation Electrification Conference Expo, ITEC 2017* 45–50. https://doi.org/10.1109/ITEC.2017.7993245. (2017).

11. Brand, M. J., Hofmann, M. H., Schuster, S. S., Keil, P. & Jossen, A. The influence of current ripples on the lifetime of lithium-ion batteries. *IEEE Trans. Veh. Technol.* **67**, 10438–10445 (2018).

12. Bessman, A., Soares, R., Wallmark, O., Svens, P. & Lindbergh, G. Aging effects of AC harmonics on lithium-ion cells. *J. Energy Storage* **21**, 741–749 (2019).

13. Uddin, K., Moore, A. D., Barai, A. & Marco, J. The effects of high frequency current ripple on electric vehicle battery performance. *Appl. Energy* **178**, 142–154 (2016).

14. Juang, L. W. et al. Investigation of the influence of superimposed AC current on lithium-ion battery aging using statistical design of experiments. *J. Energy Storage* **11**, 93–103 (2017).

15. Uno, M. & Tanaka, K. Influence of high-frequency charge-discharge cycling induced by cell voltage equalizers on the life performance of lithium-ion cells. *IEEE Trans. Veh. Technol.* **60**, 1505–1515 (2011).

16. Lü, X. et al. Energy management of hybrid electric vehicles: A review of energy optimization of fuel cell hybrid power system based on genetic algorithm. *Energy Convers. Manag.* **205**, 112474 (2020).

17. Baba, M. A., Labbadi, M., Cherkaoui, M. & Maaroufi, M. Fuel cell electric vehicles: A review of current power electronic converters topologies and technical challenges. In *IOP Conference Series Earth Environmental Science* **785**, 1–29 (2021).
18. Eddine, G. H., Said, B., Houari, A., Dieroui, A. & Mesbahi, T. Integral sliding mode control of synchronous reluctance machine based electric vehicle powered by battery/supercapacitor hybrid source. In 2022 19th International Multi-Conference Systems Signals Devices 2133–2138. https://doi.org/10.1109/SSD54932.2022.9955985. (2022).

19. Ansari, M. M. A., Cronje, W. A. & Meyer, A. Evaluation of a reluctance synchronous motor: For use in an electric mine shuttle vehicle (EMSV). In 2012 IEEE International Electric Vehicle Conference IEVC 2012. https://doi.org/10.1109/IEVC.2012.6183195. (2012).

20. Taghavi, S. & Pillay, P. A sizing methodology of the synchronous reluctance motor for traction applications. *IEEE J. Emerg. Sel. Top. Power Electron.* **2**, 329–340 (2014).

21. Kumar, G. V., Chuang, C. H., Lu, M. Z. & Liaw, C. M. Development of an electric vehicle synchronous reluctance motor drive. *IEEE Trans. Veh. Technol.* **69**, 5012–5024 (2020).

22. Herrera, D. B., Galvan, E. & Carrasco, J. M. Synchronous reluctance motor design based EV powertrain with inverter integrated with redundant topology. In *IECON 2015—41st Annual Conference IEEE Industrial Electronics Society* 3851–3856. https://doi.org/10.1109/IECON.2015.7392700. (2015).

23. Kim, Y. H. & Lee, J. H. Optimum design of ALA-SynRM for direct drive electric valve actuator. *IEEE Trans. Magn.* **53**, 1–4 (2017).

24. Peng, H., Li, J., Löwenstein, L. & Hameyer, K. A scalable, causal, adaptive energy management strategy based on optimal control theory for a fuel cell hybrid railway vehicle. *Appl. Energy* **267**, 114980 (2020).

25. Rimpas, D. et al. Materials today: Proceedings energy management and storage systems on electric vehicles: A comprehensive review. *Mater. Today Proc.* **61**, 813–819 (2022).

26. Bauer, S., Suchaneck, A. & Puente León, F. Thermal and energy battery management optimization in electric vehicles using Pon-tryagin's maximum principle. *J. Power Sour.* **246**, 808–818 (2014).

27. Song, K., Wang, X., Li, F., Sorrentino, M. & Zheng, B. Pontryagin's minimum principle-based real-time energy management strategy for fuel cell hybrid electric vehicle considering both fuel economy and power source durability. *Energy* **205**, 118064 (2020).

28. Tona, T. L. F. V. P. Stochastic dynamic programming based energy management of HEV's: An experimental validation. In *IFAC Proceedings Volumes* Vol. 47 (IFAC, 2014).

29. Zhang, S. & Xiong, R. Adaptive energy management of a plug-in hybrid electric vehicle based on driving pattern recognition and dynamic programming. *Appl. Energy* **155**, 68–78 (2015).

30. Lü, X. et al. Hybrid electric vehicles: A review of energy management strategies based on model predictive control. *J. Energy Storage* **56**, 106112 (2022).

31. Lei, Z. et al. An adaptive equivalent consumption minimization strategy for plug-in hybrid electric vehicles based on traffic information. *Energy* **190**, 116409 (2020).

32. Dhifli, M., Lashab, A., Guerrero, J. M. & Al-turki, Y. A. An efficient external energy maximization-based energy management strategy for a battery/supercapacitor of a micro grid system enhanced intelligent energy management system for a renewable energy-based AC microgrid. *Energies* https://doi.org/10.3390/en13123268 (2020).

33. Ye, Y., Zhang, J., Pilla, S., Rao, A. M. & Xu, B. Application of a new type of lithium-sulfur battery and reinforcement learning in plug-in hybrid electric vehicle energy management. *J. Energy Storage* **59**, 106546 (2023).

34. Li, W. et al. Deep reinforcement learning-based energy management of hybrid battery systems in electric vehicles. *J. Energy Storage* **36**, 102355 (2021).

35. Wu, C. et al. The application of machine learning based energy management strategy in multi-mode plug-in hybrid electric vehicle, part I: Twin delayed deep deterministic policy gradient algorithm design for hybrid mode. *Energy* **262**, 125084 (2023).

36. Qi, C., Song, C., Xiao, F. & Song, S. Generalization ability of hybrid electric vehicle energy management strategy based on reinforcement learning method. *Energy* **250**, 123826 (2022).

37. Roche, M., Shabbir, W. & Evangelou, S. A. Voltage control for enhanced power electronic efficiency in series hybrid electric vehicles. *IEEE Trans. Veh. Technol.* **66**, 3645–3658 (2017).

38. Ferahtia, S., Djerioui, A., Zeghlache, S. & Houari, A. A hybrid power system based on fuel cell, photovoltaic source and supercapacitor. *SN Appl. Sci.* **2**, 1–11 (2020).

39. Yodwong, B., Thounthong, P., Guilbert, D. & Bizon, N. Differential flatness-based cascade energy/current control of battery/supercapacitor hybrid source for modern e-vehicle applications. *Mathematics* **8**, 704 (2020).

40. Mirjalili, S. et al. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. *Adv. Eng. Softw.* **114**, 163–191 (2017).

41. Ferahtia, S. et al. Optimal adaptive gain LQR-based energy management strategy for battery-supercapacitor hybrid power system. *Energies* **14**, 1660 (2021).

42. Babu, T. S. et al. A comprehensive review of hybrid energy storage systems: Converter topologies, control strategies and future prospects. *IEEE Access* **8**, 148702–148721 (2020).

43. Thounthong, P. et al. Performance investigation of high-energy high-power densities storage devices by li-ion battery and supercapacitor for fuel cell/photovoltaic hybrid power plant for autonomous system applications. In *IEEE Industry Applications Society 51st Annual Meeting IAS 2015, Conference Record*. https://doi.org/10.1109/IAS.2015.7356844 (2015).

44. Thounthong, P., Tricoli, P. & Davat, B. Performance investigation of linear and nonlinear controls for a fuel cell/supercapacitor hybrid power plant. *Int. J. Electr. Power Energy Syst.* **54**, 454–464 (2014).

45. Mukherjee, N. & Strickland, D. Control of cascaded DC-DC converter-based hybrid battery energy storage systems-part I: Stability issue. *IEEE Trans. Ind. Electron.* **63**, 2340–2349 (2016).

46. Thounthong, P. et al. Nonlinear differential flatness-based speed/torque control with state-observers of permanent magnet synchronous motor drives. *IEEE Trans. Ind. Appl.* **54**, 2874–2884 (2018).

47. Sankarkumar, R. S. & Natarajan, R. Energy management techniques and topologies suitable for hybrid energy storage system powered electric vehicles: an overview. *Int. Trans. Electr. Energy Syst.* **31**, 1–30 (2021).

48. Zhang, Q., Deng, W., Zhang, S. & Wu, J. A rule based energy management system of experimental battery/supercapacitor hybrid energy storage system for electric vehicles. *J. Control Sci. Eng.* **2016**, 18–20 (2016).

49. Xun, Q., Roda, V., Liu, Y., Huang, X. & Costa-Castelló, R. An adaptive power split strategy with a load disturbance compensator for fuel cell/supercapacitor powertrains. *J. Energy Storage* **44**, 103341 (2021).

50. Feng, N., Ma, T. & Chen, C. Fuzzy energy management strategy for hybrid electric vehicles on battery state-of-charge estimation by particle filter. *SN Appl. Sci.* **4**, 256 (2022).

51. Song, Z., Hofmann, H., Li, J., Han, X. & Ouyang, M. Optimization for a hybrid energy storage system in electric vehicles using dynamic programming approach. *Appl. Energy* **139**, 151–162 (2015).

52. Yi, F. et al. Energy management strategy for hybrid energy storage electric vehicles based on Pontryagin's minimum principle considering battery degradation. *Sustain* **14**, 1214 (2022).

53. Fu, J., Fu, Z. & Song, S. Energy management strategy integrating self-adaptive adjustment and Pontryagin's minimum principle-based optimization for fuel-cell hybrid electric vehicle. *Optim. Control Appl. Methods* **43**, 1059–1075 (2022).

54. Wieczorek, M. & Lewandowski, M. A mathematical representation of an energy management strategy for hybrid energy storage system in electric vehicle and real time optimization using a genetic algorithm. *Appl. Energy* **192**, 222–233 (2017).

55. Hussaian Basha, C. H. et al. Design of GWO based fuzzy MPPT controller for fuel cell fed EV application with high voltage gain DC-DC converter. *Mater. Today Proc.* **92**, 66–72 (2023).
18. Eddine, G. H., Said, B., Houari, A., Dieroui, A. & Mesbahi, T. Integral sliding mode control of synchronous reluctance machine based electric vehicle powered by battery/supercapacitor hybrid source. In 2022 19th International Multi-Conference Systems Signals Devices 2133–2138. https://doi.org/10.1109/SSD54932.2022.9955985. (2022).

19. Ansari, M. M. A., Cronje, W. A. & Meyer, A. Evaluation of a reluctance synchronous motor: For use in an electric mine shuttle vehicle (EMSV). In 2012 IEEE International Electric Vehicle Conference IEVC 2012. https://doi.org/10.1109/IEVC.2012.6183195. (2012).

20. Taghavi, S. & Pillay, P. A sizing methodology of the synchronous reluctance motor for traction applications. *IEEE J. Emerg. Sel. Top. Power Electron.* **2**, 329–340 (2014).

21. Kumar, G. V., Chuang, C. H., Lu, M. Z. & Liaw, C. M. Development of an electric vehicle synchronous reluctance motor drive. *IEEE Trans. Veh. Technol.* **69**, 5012–5024 (2020).

22. Herrera, D. B., Galvan, E. & Carrasco, J. M. Synchronous reluctance motor design based EV powertrain with inverter integrated with redundant topology. In *IECON 2015—41st Annual Conference IEEE Industrial Electronics Society* 3851–3856. https://doi.org/10.1109/IECON.2015.7392700. (2015).

23. Kim, Y. H. & Lee, J. H. Optimum design of ALA-SynRM for direct drive electric valve actuator. *IEEE Trans. Magn.* **53**, 1–4 (2017).

24. Peng, H., Li, J., Löwenstein, L. & Hameyer, K. A scalable, causal, adaptive energy management strategy based on optimal control theory for a fuel cell hybrid railway vehicle. *Appl. Energy* **267**, 114980 (2020).

25. Rimpas, D. et al. Materials today: Proceedings energy management and storage systems on electric vehicles: A comprehensive review. *Mater. Today Proc.* **61**, 813–819 (2022).

26. Bauer, S., Suchaneck, A. & Puente León, F. Thermal and energy battery management optimization in electric vehicles using Pon-tryagin's maximum principle. *J. Power Sour.* **246**, 808–818 (2014).

27. Song, K., Wang, X., Li, F., Sorrentino, M. & Zheng, B. Pontryagin's minimum principle-based real-time energy management strategy for fuel cell hybrid electric vehicle considering both fuel economy and power source durability. *Energy* **205**, 118064 (2020).

28. Tona, T. L. F. V. P. Stochastic dynamic programming based energy management of HEV's: An experimental validation. In *IFAC Proceedings Volumes* Vol. 47 (IFAC, 2014).

29. Zhang, S. & Xiong, R. Adaptive energy management of a plug-in hybrid electric vehicle based on driving pattern recognition and dynamic programming. *Appl. Energy* **155**, 68–78 (2015).

30. Lü, X. et al. Hybrid electric vehicles: A review of energy management strategies based on model predictive control. *J. Energy Storage* **56**, 106112 (2022).

31. Lei, Z. et al. An adaptive equivalent consumption minimization strategy for plug-in hybrid electric vehicles based on traffic information. *Energy* **190**, 116409 (2020).

32. Dhifli, M., Lashab, A., Guerrero, J. M. & Al-turki, Y. A. An efficient external energy maximization-based energy management strategy for a battery/supercapacitor of a micro grid system enhanced intelligent energy management system for a renewable energy-based AC microgrid. *Energies* https://doi.org/10.3390/en13123268 (2020).

33. Ye, Y., Zhang, J., Pilla, S., Rao, A. M. & Xu, B. Application of a new type of lithium-sulfur battery and reinforcement learning in plug-in hybrid electric vehicle energy management. *J. Energy Storage* **59**, 106546 (2023).

34. Li, W. et al. Deep reinforcement learning-based energy management of hybrid battery systems in electric vehicles. *J. Energy Storage* **36**, 102355 (2021).

35. Wu, C. et al. The application of machine learning based energy management strategy in multi-mode plug-in hybrid electric vehicle, part I: Twin delayed deep deterministic policy gradient algorithm design for hybrid mode. *Energy* **262**, 125084 (2023).

36. Qi, C., Song, C., Xiao, F. & Song, S. Generalization ability of hybrid electric vehicle energy management strategy based on reinforcement learning method. *Energy* **250**, 123826 (2022).

37. Roche, M., Shabbir, W. & Evangelou, S. A. Voltage control for enhanced power electronic efficiency in series hybrid electric vehicles. *IEEE Trans. Veh. Technol.* **66**, 3645–3658 (2017).

38. Ferahtia, S., Djerioui, A., Zeghlache, S. & Houari, A. A hybrid power system based on fuel cell, photovoltaic source and supercapacitor. *SN Appl. Sci.* **2**, 1–11 (2020).

39. Yodwong, B., Thounthong, P., Guilbert, D. & Bizon, N. Differential flatness-based cascade energy/current control of battery/supercapacitor hybrid source for modern e-vehicle applications. *Mathematics* **8**, 704 (2020).

40. Mirjalili, S. et al. Salp Swarm Algorithm: A bio-inspired optimizer for engineering design problems. *Adv. Eng. Softw.* **114**, 163–191 (2017).

41. Ferahtia, S. et al. Optimal adaptive gain LQR-based energy management strategy for battery-supercapacitor hybrid power system. *Energies* **14**, 1660 (2021).

42. Babu, T. S. et al. A comprehensive review of hybrid energy storage systems: Converter topologies, control strategies and future prospects. *IEEE Access* **8**, 148702–148721 (2020).

43. Thounthong, P. et al. Performance investigation of high-energy high-power densities storage devices by li-ion battery and supercapacitor for fuel cell/photovoltaic hybrid power plant for autonomous system applications. In *IEEE Industry Applications Society 51st Annual Meeting IAS 2015, Conference Record*. https://doi.org/10.1109/IAS.2015.7356844 (2015).

44. Thounthong, P., Tricoli, P. & Davat, B. Performance investigation of linear and nonlinear controls for a fuel cell/supercapacitor hybrid power plant. *Int. J. Electr. Power Energy Syst.* **54**, 454–464 (2014).

45. Mukherjee, N. & Strickland, D. Control of cascaded DC-DC converter-based hybrid battery energy storage systems-part I: Stability issue. *IEEE Trans. Ind. Electron.* **63**, 2340–2349 (2016).

46. Thounthong, P. et al. Nonlinear differential flatness-based speed/torque control with state-observers of permanent magnet synchronous motor drives. *IEEE Trans. Ind. Appl.* **54**, 2874–2884 (2018).

47. Sankarkumar, R. S. & Natarajan, R. Energy management techniques and topologies suitable for hybrid energy storage system powered electric vehicles: an overview. *Int. Trans. Electr. Energy Syst.* **31**, 1–30 (2021).

48. Zhang, Q., Deng, W., Zhang, S. & Wu, J. A rule based energy management system of experimental battery/supercapacitor hybrid energy storage system for electric vehicles. *J. Control Sci. Eng.* **2016**, 18–20 (2016).

49. Xun, Q., Roda, V., Liu, Y., Huang, X. & Costa-Castelló, R. An adaptive power split strategy with a load disturbance compensator for fuel cell/supercapacitor powertrains. *J. Energy Storage* **44**, 103341 (2021).

50. Feng, N., Ma, T. & Chen, C. Fuzzy energy management strategy for hybrid electric vehicles on battery state-of-charge estimation by particle filter. *SN Appl. Sci.* **4**, 256 (2022).

51. Song, Z., Hofmann, H., Li, J., Han, X. & Ouyang, M. Optimization for a hybrid energy storage system in electric vehicles using dynamic programing approach. *Appl. Energy* **139**, 151–162 (2015).

52. Yi, F. et al. Energy management strategy for hybrid energy storage electric vehicles based on Pontryagin's minimum principle considering battery degradation. *Sustain* **14**, 1214 (2022).

53. Fu, J., Fu, Z. & Song, S. Energy management strategy integrating self-adaptive adjustment and Pontryagin's minimum principle-based optimization for fuel-cell hybrid electric vehicle. *Optim. Control Appl. Methods* **43**, 1059–1075 (2022).

54. Wieczorek, M. & Lewandowski, M. A mathematical representation of an energy management strategy for hybrid energy storage system in electric vehicle and real time optimization using a genetic algorithm. *Appl. Energy* **192**, 222–233 (2017).

55. Hussaian Basha, C. H. et al. Design of GWO based fuzzy MPPT controller for fuel cell fed EV application with high voltage gain DC-DC converter. *Mater. Today Proc.* **92**, 66–72 (2023).
56. Xiong, R., Cao, J. & Yu, Q. Reinforcement learning-based real-time power management for hybrid energy storage system in the plug-in hybrid electric vehicle. *Appl. Energy* **211**, 538–548 (2018).

57. Automatisés, S. et al. Artificial neural networks based power management for a battery/supercapacitor and integrated photovoltaic hybrid storage system for electric vehicles. *J. Eur. Syst. Autom.* **56**, 139–151 (2023).

**Author contributions**

Conceptualization, H.E.G., S.B., and A.H.; methodology, A.D., and A.H.; software, H.E.G., S.B., and A.D.; validation, S.B., and A.H.; formal analysis, M.O., and N.B.; investigation, H.E.G., S.B., A.H., and A.D.; resources, M.O., and N.B.; data curation, M.O., and N.B.; writing—original draft preparation, H.E.G., and S.B., writing—review and editing, A.H., A.D., and N.B.; visualization, M.O., and N.B.; supervision, A.H.; project administration, M.O., and N.B.; funding acquisition, M.O., and N.B. All authors have read and agreed to the published version of the manuscript.

**Funding**

The research was carried out with partial support from Subprogramme 1.1., Institutional performance-Projects to finance excellence in RDI, Contract No. 19PFE/30.12.2021.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-53396-3.

Correspondence and requests for materials should be addressed to H.E.G. or S.B.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_24.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024
56. Xiong, R., Cao, J. & Yu, Q. Reinforcement learning-based real-time power management for hybrid energy storage system in the plug-in hybrid electric vehicle. *Appl. Energy* **211**, 538–548 (2018).

57. Automatisés, S. et al. Artificial neural networks based power management for a battery/supercapacitor and integrated photovoltaic hybrid storage system for electric vehicles. *J. Eur. Syst. Autom.* **56**, 139–151 (2023).

## Author contributions

Conceptualization, H.E.G., S.B., and A.H.; methodology, A.D., and A.H.; software, H.E.G., S.B., and A.D.; validation, S.B., and A.H.; formal analysis, M.O., and N.B.; investigation, H.E.G., S.B., A.H., and A.D.; resources, M.O., and N.B.; data curation, M.O., and N.B.; writing—original draft preparation, H.E.G., and S.B., writing—review and editing, A.H., A.D., and N.B.; visualization, M.O., and N.B.; supervision, A.H.; project administration, M.O., and N.B.; funding acquisition, M.O., and N.B. All authors have read and agreed to the published version of the manuscript.

## Funding

The research was carried out with partial support from Subprogramme 1.1., Institutional performance-Projects to finance excellence in RDI, Contract No. 19PFE/30.12.2021.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-53396-3.

Correspondence and requests for materials should be addressed to H.E.G. or S.B.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_24.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024