# Pre-Sizing Approach of a Fuel Cell-Battery Hybrid Power System with Interleaved Converters

Victor Mercier, Adriano Ceschia, Toufik Azib, Cherif Larouci

► To cite this version:

Victor Mercier, Adriano Ceschia, Toufik Azib, Cherif Larouci. Pre-Sizing Approach of a Fuel Cell-Battery Hybrid Power System with Interleaved Converters. Energies, 2023, 16 (10), pp.4068. 10.3390/en16104068 . hal-04127116

HAL Id: hal-04127116

https://hal.science/hal-04127116

Submitted on 13 Jun 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Pre-Sizing Approach of a Fuel Cell-Battery Hybrid Power System with Interleaved Converters

Victor Mercier, Adriano Ceschia, Toufik Azib * and Cherif Larouci *

ESTACA, ESTACA'Lab, Paris-Saclay, 78066 Paris, France; victor.mercier@estaca.fr (V.M.);
adriano.ceschia@estaca.fr (A.C.)

* Correspondence: toufik.azib@estaca.fr (T.A.); cherif.larouci@estaca.fr (C.L.)

**Abstract:** This paper proposes a design methodology that is dedicated to improving the concept of a modular hybrid power chain that uses interleaved converters. The approach involves optimising the system under multi-physical constraints, where the number of cells in the interleaved converters is considered as a key modular parameter. The methodology uses analytical models to strike a balance between computation time and result accuracy. This compromise is indispensable to the construction of a smart design approach under multi-physical constraints, such as electric, efficiency, volume, and thermal constraints. The proposed approach has been applied to a hybrid fuel cell and battery power system for automotive applications; the goal is to obtain a global optimal architecture chain by optimising the number of interleaved converter cells and by determining appropriate power electronics components and the optimal sizing of sources. This constitutes the primary step for providing an effective pre-design support tool for considering architecture modularity, facilitating the use of new technologies in the early design stage. The results showed that the interleaving concept allows for better flexibility in respecting the design constraints to improve the design of hybrid power systems. The analysis also highlights the current limitations and performance of the optimisation method and suggests new areas for future work.

**Keywords:** fuel cell electric vehicle; interleaved converters; conception approach; multi-physical optimisation; modular architecture

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

check for
updates

**Citation:** Mercier, V.; Ceschia, A.; Azib, T.; Larouci, C. Pre-Sizing Approach of a Fuel Cell-Battery Hybrid Power System with Interleaved Converters. *Energies* **2023**, *16*, 4068. https://doi.org/10.3390/en16104068

Academic Editor: Antonio Cano-Ortega

Received: 3 April 2023

Revised: 8 May 2023

Accepted: 8 May 2023

Published: 12 May 2023

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

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

## 1. Introduction

Global warming and air pollution are serious issues that affect our everyday life [1]. Air pollution is a significant problem in many industrialized cities worldwide, and vehicle emissions are one of the primary contributors. The transportation sector is indeed a vital aspect of human development and economic growth [2]; however, it also has negative environmental impacts that must be addressed to ensure a sustainable future. Governments, businesses, and individuals must take action to reduce their carbon footprint and promote sustainable practices in transportation and other sectors. Fortunately, an important goal in the automotive industry continues to be the development of alternative and new solutions for mobility transition, which results in a constant evolution of the industry. This is reflected by the integration of new solutions that are, in essence, based on the electrification concept. While electric vehicles (EVs) and hybrid electric vehicles (HEVs) have certainly made progress towards reducing emissions, they still face the limitations of battery technology and charging infrastructure: autonomy, security, refueling, durability, etc., hence the importance of introducing other promising solutions for mobility transition [3,4].

The fuel cell electric vehicle (FCEV) is one of the key solutions for this transition; it is often based on the concept of hybridization, particularly the use of PEMFC technology (proton exchange membrane fuel cell) as the main source and a battery or supercapacitor as the electric assistance [5]. Compared with conventional EVs, this class of vehicle allows for a better driving range and faster refueling time, but it requires the construction of a complex
system. Improving the design process and system integration is essential for increasing the market penetration of FCEVs and for overcoming some of the challenges that are associated with this technology, such as efficiency, durability/reliability, and cost. This implies the need for design methodologies to have the ability to consider the multi-physical constraints and available technologies, topologies, and strategies (control/energy management) to achieve the system integration according to the considered performance indexes [6].

The concept of hybridization allows for the possibility of having multiple architectures. Each architecture differs from the converter(s) used; it constitutes an essential element, which is either used to increase the voltage of the energy source or lower it and have the appropriate degrees of freedom for energy management strategy. Several converters have been introduced in the literature for such an application (boost, buck, buck-boost, etc.) In [7–9], brief studies on the different power conversion topologies were presented before the designing of a multi-phase boost converter. These standard structures do not provide any degree of freedom for the sizing, energy efficiency, continuity of service, etc. Another promising concept that is emerging is named the interleaving technique; it consists of the parallel arrangement of switching cells [10–12], which allows for the power to flow through the different cells. Figure 1 illustrates the representation of an interleaved boost converter.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Diagram of the N-cell interleaved boost converter.</figcaption></figure>

This construction allows for a stress reduction on each cell (active and passive components), which increases the possibilities for technical choices and thus favors the use of more "standard" components. Another advantage of this type of architecture is the apparent frequency that is increasing inside the converter, as shown in Figure 2, which proportionally increases with the number of cells, thus reducing the need for filtering. Refs. [11,13,14] present the different benefits of an interleaving boost converter in terms of efficiency, fault tolerance, etc. The interleaving concept also allows for a current ripple reduction at the input of the converter, which preserves the energy source from high current fluctuations. This principle provides a better flexibility for respecting the design constraints, and it expands the margin between the limit and its boundary, thus making it possible for the technology to operate closer to its limits without exceeding them; this can allow for a more compact and cost-effective design.

However, interleaving implies the implementation of an appropriate control that allows for a balanced distribution between the cells and allows for the ability to perform reconfiguration to benefit from the modularity [15]. It consists of a close control that follows the same duty cycle ($\alpha$) and is regularly shifted from one cell to another ($T_d/q$).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Waveforms of a three-phase boost converter.</figcaption></figure>

It can also be noted that the complexity of the powertrain is amplified by hybridization, leading to various additional parameters, constraints, and multidisciplinary issues. This includes a strong interaction between the components, new driving conditions, and other challenges. To address this system complexity, numerous performance indexes and a range of factors must be considered, such as the most suitable components, optimal sizes, and associated energy management strategy; hence there is an industrial challenge to improve the design process and enhance the performance and efficiency. This paper deals with these challenges by developing an effective sizing approach that also addresses the modularity concept by introducing the interleaved architecture.

Compared with other pre-sizing approaches that only focus on the interleaved converter [16–18], the present design deals with a global hybrid architecture that employs a fuel cell by addressing the interaction between the sources and the associated power electronics; these interactions can have a significant impact on the performance of the hybrid powertrain and must be carefully managed to ensure optimal operation, which can be conducted by offering a more effective solution for optimising the performances and efficiency.

This paper is structured as follows: the first part will detail the existing and used topologies in the current studies of the field. Another section explains all the representations used in this work. Finally, some results and an analysis will be exposed in the last part.

## 2. State-of-the-Art

After explaining the scientific concept of this study, this section will discuss the state-of-the-art in the field considered, starting with the global architecture and then zooming in on the components, namely the sources and converters within the power chain.

### 2.1. *Architecture/Components*

In this context, many FCEV architectures are possible, which differ according to the energy sources chosen and the associated power electronics. Several topologies are possible, as explained in [19,20]. An FCEV is usually composed of a fuel cell (FC), which is considered to be an energy device with low dynamics (cannot provide power quickly), along with an energy-storage system (ESS). The ESS can comprise various equipment, such as a battery pack, ultra-capacitors, etc. For our use case, we chose to have a battery pack as a power/energy storage system to create a hybrid architecture with FC downsizing that is not supported by an ultracapacitor (power storage system). The three structures are mainly used according to [21]:
* The two-converter structure, consisting of the association of a converter for each source. Although it is the arrangement with the largest number of degrees of freedom, it remains a complex system in terms of control design.

* The direct parallel structure, which directly connects the sources to the load. This disposition, however, requires an over-sized FC to provide the demanded power to the load.

* The one-converter structure, which only places one converter at one source (usually the FC) to downsize the associated source. It presents some constraints related to the management of the battery state of charge.

Regarding the technology of the sources, many different fuel cell technologies have been developed [22]; these are typically based on the operating temperature and electrolyte type used. In the transport sector, low-temperature fuel cells are preferred, such as the PEMFC, which has a low operating temperature (30–100 °C) and a solid membrane that absorbs system vibrations [20,23–26]. In terms of modelling, there are several different approaches depending on the level of detail required and the specific application (chemical, electrochemical, efficiency, hydrogen consumption, etc.). Static models are one type of fuel cell model that are commonly used based on steady-state operating conditions. They can be relatively simple, with only a few parameters, or they can be more complex, with multiple variables and inputs [27].

On the other hand, modelling batteries is also a complex task that requires the consideration of various parameters, including the type of battery, its chemistry, and its operating conditions. The choice of an appropriate battery model depends on the application and purpose for which the battery will be used. There are several different types of battery models that can be used, each with their own advantages and disadvantages. Some common types of battery models include empirical models, equivalent circuit models, and electrochemical models, as explicated in [28].

Concerning the power converters, as in this case, both energy sources provided a DC voltage, i.e., only DC-DC converters were treated. There are different types of DC-DC converters that are well explained in the literature [29,30] (step-down, step-up, galvanically isolated or not, reversible or not, etc.). The choice of the appropriate DC-DC converter for a fuel cell/battery system depends on several parameters, including the input and output voltage, power rating, efficiency, size and weight, etc.

## 2.2. *Design Approach*

In the literature, several studies have proposed various approaches to design a hybrid system that combines different power sources; two methods have emerged in HEV design applications: empirical approaches and computational approaches, as explained in [31]. With the rise of computer science, optimisation-based approaches have become increasingly popular for the design of hybrid systems. Compared with the traditional approaches that make a solution according to the rated conditions, systemic approaches involving a sequential process of individual optimisations or multi-objective optimisation methods have also been introduced; they provide more comprehensive solutions that take into account multiple performance criteria and constraints. Their implementation depends on the compromise between the complexity of the structure, the computation time, and the exploration capabilities. On this topic, many optimisation algorithms have emerged [32]; among them, two relevant techniques were selected: the particle swarm optimisation (PSO) technique and the genetic algorithm (GA) [33]. Both techniques are rather simple to implement and fast to process and a comparative study is essential in order to make a choice between them according to the considered requirements.
## 3. Implementation and Modelling

### 3.1. *Approach Description*

Our use case was based on a hybrid fuel cell/battery architecture with two boost converters. The load was represented as a power demand based on a small city car driving through a preset circuit. The pre-sizing approach here, shown in Figure 3, involves an optimisation process with multi-physic constraints. These two algorithms (PSO and GA) have been separately tested and compared in terms of computing time, convergence speed, and accuracy; it results in a better convergence speed for the GA, although its computing time for one iteration is clearly above that of the PSO. With respect to the accuracy, the PSO provides better results and remains more accessible to manipulating and changing its internal parameters. As it remains more flexible to exploit, the PSO algorithm was used here; it realises the optimisation of the main sizing parameters, which were tested within an objective function for minimisation purposes. This optimisation was constrained by different condition limits according to the architecture's integrity and specifications. The optimisation loop will experiment a defined size for each particle according to a considered driving cycle. After reaching the end of the cycle, some physical constraints (for example, the SoC deviation throughout the simulation must be less than one) are evaluated and added to the cost function as penalties. If the constraint is not respected, it adds a cost to the objective function to exclude this solution. The load power profile is determined using the longitudinal dynamics model, which takes into account the vehicle speed and forces resulting from acceleration, drag, friction, and road slope [34].

<figure><img src="image_6.png" /><figcaption>Figure 3. Overview of the suggested design approach.</figcaption></figure>

There are several performance indexes that can be chosen as the cost function, such as the overall volume, components lifespan, efficiency, etc. For this primary step of our support tool, the objective function to minimise here was the system global volume as
it remains an industrial challenge in the development of FCEVs. Through the article, the different constraints associated with the cost function will be developed as electrical, thermal, and efficiency constraints.

To implement the optimisation problem, the first step is to model the volume of each component to evaluate the cost function. Then, the behaviour of the sources and converters are modelled. After that, the various ripples inside the converter will be described. To express the converter's efficiency, the losses in terms of the electrical losses will be developed along with the thermal behaviour, which is also considered as a constraint. The optimisation algorithm will be developed and explained in the last part. All of the following implementations will follow the described formalism contained in Table 1.

<table><caption>**Table 1.** Optimisation problem formulation.</caption><thead><tr><th>Formalism</th><th>Mathematical Definition</th></tr></thead><tbody><tr><td>Cost function</td><td>F(X) = Volume<sub>Archi</sub></td></tr><tr><td>Optimisation parameters</td><td>Converter parameters, FC parameters, and BAT parameters</td></tr><tr><td>Constraints</td><td>Ripples &lt; Ripple_max<br>Efficiency &gt; η<sub>min</sub><br>Semi-conductor temperature &lt; T_max<br>P<sub>FC_min</sub> &lt; FCMaxPower &lt; P<sub>FC_max</sub></td></tr></tbody></table>

## 3.2. Volume Cost Function

As the purpose here is to focus on the interleaving boost converter architecture, the principal criterion about this part is the overall volume, including the volume of the sources and converters. This study also gathered the constraints related to many different domains (electrical, thermal, etc.) and therefore remains as an interesting case study. This converter volume includes the passive and active components. The volume of the active components is mainly represented by the heat sink volume. Regarding the converter volume, the volume of the passive and active components volume are considered here, with the heat sink volume largely constituting the volume of the active components. Analytical models were computed to estimate the volume of inductance, capacitor, heat sink, and energy sources.

### 3.2.1. Fuel Cell Volume

The fuel cell sizing mainly depends on the application area and operating conditions (static or embedded environment). Modelling the fuel cell volume can be a complex task as it involves taking into account various factors. For PEMFC technology, it is estimated that the volume of the stack takes 30% of the overall system. Consequently, 70% is held by the auxiliaries. It is common to refer to the Ragone diagram when estimating the system volume. Using the data extracted from this diagram as well as the data from existing vehicles such as the Toyota Mirai [35], it is possible to estimate that each kW of fuel cell power required is equivalent to a volume increase of 1.5625 dm³. A fuel cell with a maximum power of 10 kW will then have a volume of about 15.6 L, i.e.,

$$Vol_{FC} = 1.5625 \times 10^{-3} P_{FCmax} \qquad (1)$$

### 3.2.2. Battery Volume

The battery pack of an electric vehicle is an assembly of serial or parallel electric cells. Therefore, based on the cell used in the battery pack of the Renault Zoé, the LG Chem E63, it is possible to estimate the volume of the battery pack as a function of the cell's dimensions, i.e.,

$$Vol_{BAT} = ((length \cdot width \cdot thickness) \cdot n_{cell_{series}}) \cdot n_{stack_{para}} \qquad (2)$$
where length = 325 mm, width = 125 mm, and thickness = 11.5 mm are the cell’s dimensions, $n_{cell,series}$ is the number of serial cells of one stack, and $n_{stack,para}$ is the number of parallel stacks.

### 3.2.3. Converter Volume

The converter volume is composed of the inductance, capacitor, and active components volume. The heat sink volume represents the volume of the active components. As several switching cells can be interconnected, each volume model presented here is given with respect to an elementary cell.

**Inductance volume:**

The volume of an inductance is proportional to the stored energy, i.e.,

$$Vol_L \propto k_L \cdot L \cdot I_{Lmoy}^2 \quad (3)$$

where $k_L$ is a constant factor that describes the relationship between the stored energy and the coil volume, $L$ is the inductance value, and $I_{Lmoy}$ is the average current flowing through the inductance. $k_L$ can be calculated from the manufacturer’s databases for different magnetic core sizes, inductance values, and current values. However, an approximation of the volume can also be written by [36]:

$$Vol_L = k_{L1} \cdot L \cdot I_{Lmoy}^2 + k_{L2} \cdot L + k_{L3} \cdot I_{Lmoy} \quad (4)$$

where the parameters $k_{L1}$, $k_{L2}$, and $k_{L3}$ depend on the shape of the coil. Here, only toroidal cores are considered because of their high volume-to-inductance ratio.

**Capacitor volume:**

The volume of a capacitor depends on the geometry and material of the component. It is also proportional to the stored energy plus a factor representing the voltage between the capacitor terminals [36]:

$$Vol_C = k_{C1} \cdot C \cdot V_C^2 + k_{C2} \quad (5)$$

where the parameter $k_{C1}$ describes the relationship of the capacitor volume and the stored energy and $k_{C2}$ is a voltage-dependent factor. These factors can be determined using the data of the manufacturers. C is the capacitance value and $V_C$ is the voltage of the capacitor.

**Active component volume:**

The heat sink is chosen to respect a maximum junction temperature of 130 °C for the semiconductors and occupy a volume that is as small as possible. Furthermore, the volume of the diode and controlled switch are negligible in front of the heat sink volume. The model used to calculate the volume of the heat sink as a function of the thermal resistance was developed from the work of [37]. This gives the following formula:

$$Vol_{Rad} = C1 \cdot e^{C2 \cdot Rth_{Rad}} \quad (6)$$

where the coefficients C1 and C2 are provided in the manufacturers’ charts [10].

### 3.3. Constraints

After illustrating the volume of the different components, the different multi-physics constraints (electrical, thermal, etc.) will be explained in this section as part of the components modelling (sources and converters).
### 3.3.1. Fuel Cell System

**Behaviour:**

There are several fuel cell modelling approaches that are based on different parameters, which depend on the level of detail and complexity required for the application (electrochemical, efficiency, hydrogen consumption, etc.), as previously explained. FC modelling allows us to identify which parameters have the greatest impact for a future optimisation. From a systemic point of view and for energy considerations, the number of cells in the FC system and the active surface area of each cell are the two most important parameters. These parameters will be used in the formulation of the optimisation problem. In this context, the FC model used here was based on the electrochemical model described in [23]. The fuel cell for one cell is modelled as:

$$V_{cell} = V_{r} - V_{act} - V_{ohm} - V_{con} \quad (7)$$

with

$$V_{r} = 1.229 - 8.5 \cdot 10^{-4}(T_{FC} - 298.15) + 4.31 \cdot 10^{-5}T_{FC} \ln(P_{H_{2}}P_{O_{2}}^{0.5}) \quad (8)$$

$$V_{act} = \frac{RT_{FC}}{z\alpha F} \ln\left(\frac{i + i_{n}}{i_{0}}\right) \quad (9)$$

$$V_{ohm} = (i + i_{n})r \quad (10)$$

$$V_{con} = -\frac{RT_{FC}}{zF} \ln\left(1 - \frac{i + i_{n}}{i_{l}}\right) \quad (11)$$

where $T_{FC}$ is the fuel cell temperature, $P_{H_2}$ and $P_{O_2}$ are the partial pressures of $H_2$ and $O_2$, $R$ is the ideal gas constant, $z$ is the number of exchanged electrons, $\alpha$ is the charge transfer coefficient, $F$ is the Faraday constant, $i$ is the current density, $i_n$ is the internal current density, $i_0$ is the current density exchanged between the electrodes, $r$ is the area-specific resistance, and $i_l$ is the limited current density. Note that the current density is a function of the active surface area.

The total voltage across the fuel cell is expressed as follows:

$$V_{FC} = N_{cell} \cdot \eta_{FC} \cdot V_{cell} \quad (12)$$

where $N_{cell}$ represents the FC stack serial cells and $\eta_{FC}$ represents the FC efficiency.

### 3.3.2. Battery Pack

**Behaviour:**

To model the battery pack, there are numerous representations that distinguish themselves from the physical constraint that they want to model. Here, temperature is not taken into account. Consequently, only electrical models are considered interesting. The battery model used is a two-time constant Thevenin model (two parallel RC networks), which provides the best accuracy, as explained in [38]. It consists of a resistor in series that is used to describe the ohmic resistance of the battery, which is associated with two parallel RC networks (polarization aspects) as shown in Figure 4 [39].
<figure><img src="image_7.png" /><figcaption>Figure 4. Second order battery modelling.</figcaption></figure>

All of the parameters of this representation are a function of the state of charge (SoC), which is calculated as:

$$ SoC(t) = SoC(0) - \frac{1}{Q_{nom}} \int i(t) dt \quad (13) $$

with SoC(0) being the initial state of charge, Qnom being the nominal capacity, and i(t) being the current. Note that the nominal capacity here is expressed in A/s.

### 3.3.3. Converter

**Behaviour:**

The average model of the converter includes two equations, one for the converter behaviour and the other reflecting the DC bus behaviour:

$$ L_i \frac{di_{L_i}}{dt} = V_{in} - V_{out}(1 - \alpha) \quad (14) $$

$$ \text{with } i = 1, \dots, q \quad (15) $$

where $L_i$ is the elementary inductance, $i_{L_i}$ is the current flowing through the inductor, $V_{in}$ is the input voltage, $V_{out}$ is the output voltage, $\alpha$ is the duty cycle, and $q$ is the number of switching cells. Regarding the DC bus, it is a function of all of the current flowing through the capacitor; hence:

$$ C \frac{dV_{out}}{dt} = (1 - \alpha) \sum_{i=1}^{q} (i_{L_i}) - i_{Load} \quad (16) $$

where C is the capacitor value and $i_{Load}$ is the current demand from the load. The capacitor is common to both sources, so it is important to keep in mind not to forget to sum all of the currents from each source.

**Electrical constraints:**

The examined electrical constraints include the input inductance current ripple, output capacitor voltage ripples, and current ripples in the elementary cell inductors. These are the main ripples that are studied in this system. Additionally, models related to electrical quantities, such as the mean, root mean square (RMS), and maximum currents, were employed to determine the losses and thermal constraints. Note that this work only focuses on the continuous conduction mode (CCM), i.e., the discontinuous conduction mode (DCM) will not be discussed here. It is also assumed that all inductance values are assumed to be the same, i.e., $L_1 = L_2 = \dots = L_i = L$. Therefore, the provided continuous conduction mode modelling is formulated as follows:

Inductance current ripple:

$$ \Delta I_L = \frac{V_{out}(1 - \alpha)\alpha}{Lf} \quad (17) $$
Output voltage ripple:

$$ \Delta V_{out} = \frac{I_{out}\alpha}{Cf} \qquad (18) $$

Input current ripple:

$$ \Delta I_{in} = \frac{\alpha \cdot (1 - q\alpha) V_{out}}{Lf} \qquad (19) $$

where $V_{out}$ is the voltage at the output of the converter, $\alpha$ is the duty cycle, $L$ is the inductance value, $f$ is the commutation frequency for one cell, $C$ is the capacitor value, $I_{out}$ is the output current of the converter, and $q$ is the number of cells. Note that because the control strategy is conducted without overlapping, (19) is only valid when $(1/q) < \alpha$. Time-domain simulations were utilised to validate these electrical models.

Efficiency:

To restrain the converter’s efficiency, loss analytical models for semiconductors and passive components (i.e., inductors and capacitors) were employed. They will be split into the semi-conductor losses at first, and then the inductance and capacitor losses will be described. The global efficiency formula will be given in a last part.

## Semi-Conductors Loss Models

Concerning the semi-conductors losses, there are three types of losses:

* Switching losses when opening the semi-conductor $P_{com_ON}$;

* Conduction losses $P_{cond}$;

* Switching losses when closing $P_{com_OFF}$;

The total losses are hence calculated as follows:

$$ P_{tot_sw} = P_{com_ON} + P_{cond} + P_{com_OFF} \qquad (20) $$

Diode commutation losses are ignored. The only losses expressed here are the conduction losses, which are represented as an equivalent circuit shown in Figure 5:

<figure><img src="image_8.png" /><figcaption>Figure 5. Equivalent circuit of a diode.</figcaption></figure>

Consequently, the diode losses are expressed as:

$$ P_{cond_d} = R_d \cdot (I_{dms})^2 + V_0 \cdot I_{dmsy} \qquad (21) $$

where $I_{dms}$ and $I_{dmsy}$ are the RMS current and average current in the diode, respectively; $R_d$ is the diode dynamic resistance and $V_0$ is the diode voltage drop. These two parameters depend on the temperature and current. They are obtained from the manufacturer’s data.

The equivalent circuit of a transistor in the conduction phase is similar to the diode’s circuit. The expression concerning the conduction losses remains almost identical except that here, the threshold voltage is negligible, thus changing Equation (21) into:

$$ P_{cond_t} = R_{ds(ON)} \cdot (I_{rms})^2 \qquad (22) $$

$$ \text{with } I_{rms} = I_L \sqrt{\alpha} \qquad (23) $$

$$ \text{thus } P_{cond_t} = R_{ds(ON)} \cdot (I_L)^2 \cdot \alpha \qquad (24) $$
where $R_{ds(ON)}$ is the switch dynamic resistance in the conduction state, $I_{trms}$ is the RMS current, and $I_L$ is the current flowing through the inductor of each cell.

The transistor commutation losses correspond to the energy lost in the switch’s activation and closing phases. By neglecting the current peaks and by assuming that the current and voltage variations are linear, we can express the switching losses as:

$$P_{com\_t} = \frac{1}{2} \cdot V_{out} \cdot I_{moy} \cdot (t_{ON} + t_{OFF}) \cdot f \quad (25)$$

where $V_{out}$ is the output voltage of the converter, $I_{moy}$ is the maximum current flowing through the transistor, $t_{ON}$ and $t_{OFF}$ are the duration of the firing and blocking phases, respectively, and $f$ is the switching frequency of the transistor.

Finally, the total semi-conductor losses are expressed by:

$$P_{tot\_sw} = P_{cond\_d} + P_{cond\_t} + P_{com\_t} \quad (26)$$

where $P_{cond\_d}$ is the diode conduction losses, $P_{cond\_t}$ is the transistor commutation losses, and $P_{com\_t}$ is the transistor conduction losses.

## Inductor Losses models

Inductance losses are composed of copper losses in the winding and core losses created by the magnetic field. Copper losses are divided into two phenomena, the skin effect losses and the proximity effect losses. The geometry of the wire (cylindrical or square) is assumed to not influence the representation of the losses. It will be modelled by:

$$P_{copper} = (K1 + K2)R_{DC}(I_{Lrms})^2 \quad (27)$$

where K1 is the resistance factor for the skin effect and K2 is for the proximity effect. $I_{Lrms}$ is the RMS current through the inductance, and $R_{DC}$ the DC resistance defined by

$$R_{DC} = \frac{\rho \cdot l}{S} \quad (28)$$

where $\rho$ is the resistivity of the material, $l$ is the length of the conductor, and $S$ is the cross-sectional area of the conductor. The K1 factor is determined using Levasseur’s formula for cylindrical conductors [40]:

$$K1 = \sqrt[6]{\left(\frac{3}{4}\right)^6 + \left(\frac{S}{Pe \cdot \delta}\right)^6 + \frac{1}{4}} \quad (29)$$

where $S$ is the cross-sectional area, $Pe$ is the perimeter of the conductor, and $\delta$ is the skin thickness, defined by

$$\delta = \sqrt{\frac{\rho}{\mu_0 \mu_r \pi f}} \quad (30)$$

where $\rho$ is the material resistivity, $\mu_0$ is the vacuum permeability, $\mu_r$ is the material permeability, and $f$ is the chopping frequency.

Core losses in an inductance are due to hysteresis and eddy current losses. These losses augment with the switching frequency and surface area of the magnetic circuit. They increase with the square of the frequency and conductivity of the material. Hysteresis losses are expressed using the formula of Steinmetz:

$$P_{core} = k V_{core} f^X B_{max}^Y \quad (31)$$

where $P_{hys}$ is the power loss, $f$ is the switching frequency, $V_{core}$ is the volume of the magnetic core, $B_{max}$ is the peak magnetic flux density, and $k$, $X$, and $Y$ are the coefficients that are specific to the material used. However, this formula is only applicable to a sinusoidal
alternating current and is therefore not suitable for our application. It then can be modified so that it is suitable in this case [41]:

$$P_{core} = k V_{core} \left[ \frac{\Delta B^m}{(2\alpha T)^n} \cdot \alpha + \frac{\Delta B^m}{(2(1-\alpha)T)^n} \cdot (1-\alpha) \right] \quad (32)$$

where $\Delta B$ is the magnetic induction change, $\alpha$ is the duty cycle, and $T$ is the switching period. The parameters $m$ and $n$ are material-specific constants, which are obviously not the same as those of Equation (31) as the core excitation is not of the same form. They are taken as $k = 1.214$, $m = 1.923$, and $n = 1.503$, as stated in [41]. When the duty cycle is 50%, this equation finds the same form as the Steinmetz Formula (31).

## Capacitor Losses

The capacitor losses are expressed using an equivalent circuit of a capacitor with an equivalent resistance as shown in Figure 6.

<figure><img src="image_9.png" /><figcaption>Figure 6. Equivalent circuit of a capacitor.</figcaption></figure>

Losses in the dielectric material are generally represented by a dissipation factor $\tan \delta_D$ constant in frequency.

$$P_C = ESR \cdot (I_{Crms})^2 \quad (33)$$

$$\text{with ESR} = \frac{\tan(\delta_D)}{C \cdot \omega} \quad (34)$$

where $\tan(\delta_D)$ is the constant frequency dissipation factor, $C$ is the capacitance of the capacitor, and $\omega = 2\pi f$ is the pulsation. The manufacturers provide charts that provide the value of these parameters, which also depend on the temperature.

## Total Losses

The efficiency of the converter is then calculated using the following model:

$$\eta = \frac{P_S}{P_S + \sum \text{losses}} = \frac{P_S}{P_S + q \cdot (P_{tot\_sw} + P_{copper} + P_{core}) + P_C} \quad (35)$$

### Thermal:

Only the major thermal stresses that are associated with semiconductors are considered here. Even if all of the components are not ideal, because of their losses, they are subject to a temperature rise. These temperature increases are negligible compared with those of semiconductors. The thermal circuit of the components is represented by a set of thermal resistances characterising each part of the component. Based on the model shown in Figure 7, when assuming a shared heat sink, it is possible to estimate the switch and diode junction temperatures ($T_{j\_sw}$ and $T_{j\_d}$) of an elementary cell. This model takes into account the switch and diode losses ($P_{sw}$ and $P_d$) as heating sources. The thermal resistances (junction case: $Rth_{jc\_sw}$ and $Rth_{jc\_d}$; case sink: $Rth_{cs\_sw}$ and $Rth_{cs\_d}$; sink air: $Rth_{sa}$) are used to simulate the thermal exchanges between the semiconductor junction and its case, between the case and the heat sink and between the heat sink and the environment.

$$T_{j\_sw} = T_a + (Rth_{jc\_sw} + Rth_{cs\_sw})P_{sw} + q \cdot Rth_{sa}(P_{sw} + P_d) \quad (36)$$
For the diode:

$$T_{j\_d} = T_a + (Rth_{jc\_d} + Rth_{cs\_d})P_d + q \cdot Rth_{sa}(P_{sw} + P_d) \quad (37)$$

<figure><img src="image_10.png" /><figcaption>Figure 7. Thermal model of a semiconductor.</figcaption></figure>

# 4. Results and Analysis

## 4.1. Use Case

Our architecture is a parallel hybrid fuel cell/battery architecture that uses a specific power converter that is connected to each source. It is well suited for optimal design and energy management, which is related to its high freedom degrees [21]. Figure 8 shows this architecture, where the fuel cell is used as the main energy source and the battery is the auxiliary energy source or ESS. The fuel cell adopted is a Ballard PEMFC. The battery pack is considered to be a lithium ion cell stack as it is the most current developed type of cell. It is treated as an auxiliary source whose main purpose is to provide peak power when demanded by the load. Indeed, the fuel cell dynamic remains quite slow and cannot provide the demanded power as fast as the battery pack would. The energy management system consists of a frequency decomposition (using a two-order Butterworth analog filter). Thus, small power fluctuations are taken by the fuel cell, and high fluctuations are absorbed by the battery.

<figure><img src="image_11.png" /><figcaption>Figure 8. Parallel hybrid power source architecture.</figcaption></figure>

The use case specifications considered are related to a city car that is similar to a Renault Zoé model, which is shown in Table 2. The load is a power demand working at a specific voltage. The chosen driving cycle is the WLTC cycle.
<table><caption>**Table 2.** Renault Zoé vehicle parameters.</caption><thead><tr><td>Parameter</td><td>Values</td></tr></thead><tbody><tr><td>M</td><td>1425 (kg)</td></tr><tr><td>L</td><td>4.804 (m)</td></tr><tr><td>l</td><td>1.730 (m)</td></tr><tr><td>H</td><td>1.568 (m)</td></tr><tr><td>S</td><td>1 × H (m²)</td></tr><tr><td>Cx</td><td>0.29</td></tr><tr><td>Cr</td><td>0.012</td></tr><tr><td>α</td><td>0</td></tr><tr><td>g</td><td>9.81 (m·s⁻²)</td></tr><tr><td>ρ</td><td>1.184 (kg·m⁻³)</td></tr></tbody></table>

The approach has been deployed and implemented in a Matlab R2020b environment, and the optimisation algorithm was self-coded to have a flexibility concerning its parameters.

The main optimisation parameters regarding the power converter are the number of cells (*q*), which are considered to be a key optimisation parameter for determining the converter architecture, switching frequency (*f*), elementary cell inductance (*L*), and output capacitance (*C*). Due to the significant impact of passive component technology on both the constraints and objective function, the technology parameters for each component are incorporated as optimisation parameters. With respect to the sources, the two parameters integrated in the optimisation process for the fuel cell are the number of cells in series and the active surface area. For the battery pack, only the number of cells that are in series are optimised here.

The optimisation aims to determine the optimal range values for minimising the objective function according to considered constraints.

## 4.2. Optimisation Results

First, to validate the interest of interleaving in a FCEV, the initial step presented here consists of analysing the benefits of the interleaving technique in a mono-source power chain. It is represented as a fuel cell and its boost converter, which is associated with the load, as shown in Figure 9:

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Mono-source power chain.</figcaption></figure>

The main specifications considered in this case are the following:

* Output network voltage: 360 V

* Output power: 20 kW

* Number of cells: 1, 2, 3, 4

Moreover, the associated problem formulation is exposed in Table 3:
<table><caption>**Table 3.** Optimisation problem formulation.</caption><thead><tr><th>Formalism</th><th>Mathematical Definition</th></tr></thead><tbody><tr><td>Cost function</td><td>F(X) = Volume<sub>Archi</sub></td></tr><tr><td>Optimisation parameters</td><td>Inductance, capacitance, frequency, magnetic density, wire diameter, winding number, heat sink thermal resistance, active surface area, number of FC cells</td></tr><tr><td>Constraints</td><td>Input current ripple &lt; 10%<br>Inductance current ripple &lt; 50%<br>Output voltage ripple &lt; 10%<br>Junction temperature MOSFET &lt; 130° C<br>Junction temperature diode &lt; 130° C<br>Efficiency &gt; 80%<br>V<sub>FC</sub> &lt; 0.9 × Output voltage<br>0.8 × P<sub>FC_max</sub> &gt; P<sub>Load</sub><br>Input converter current &lt; limited FC current<br>P<sub>FC_max</sub> &lt; 2 × P<sub>Load</sub></td></tr></tbody></table>

The results of the obtained volumes are shown in Figures 10 and 11 when the number of switching cells ($q$) is a fixed parameter that is considered as a specification.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Optimal global volume.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 11.** Converter global volume.</figcaption></figure>

According to these figures, it appears that increasing the number of switching cells results in a reduction in the overall volume of the converter. The reduction in volume is most significant when moving from a one-cell converter to a two-cell converter, which results in a reduction of almost 20% in volume compared with the previous configuration. There is a diminishing return on the volume gain as the number of switching cells is increased due to the presence of a “dead” volume in each component. This “dead” volume
cannot be made as small as possible, which limits the overall reduction in volume as more cells are added to the converter (Figure 12).

<figure><img src="image_15.png" /><figcaption>Figure 12. Optimal number of cells for different output power.</figcaption></figure>

As shown in this figure, an optimal number of cells that minimises the overall volume for each power level can be observed. As a result, for high power levels, interleaving results in a significant reduction in the converter volume. On the other hand, the reduction in volume gain is insignificant for low power levels. For the next validation, different sizing will be presented and compared. In fact, the decomposition shows a maximum power provided by the fuel cell of 40 kW with a cut-off frequency of 50 mHz, as evaluated in [42]. This will be taken as the most undersized structure ("dim1") shown in Figure 13. The other sizings are 50 kW, 60 kW, 70 kW, and 80 kW for "dim5". The parameters chosen for the converter and FC sizing of 40 kW are given in Table 4:

<table><caption>Figure 13. Comparison of different sizings with different numbers of cells data.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="5">Number of cells</th></tr><tr><th>1</th><th>2</th><th>3</th><th>4</th><th>5</th></tr><tr><th></th><th></th><th>Results</th><th>Results</th><th>Results</th><th>Results</th><th>Results</th></tr></thead><tbody><tr><td rowspan="10">Sizing</td><td rowspan="5">dim5 - 80kW</td><td>Input current ripple (%)</td><td>0.21</td><td>0.2</td><td>0.2</td><td>0.2</td><td>0.20</td></tr><tr><td>Output voltage ripple (%)</td><td>0.02</td><td>0.01</td><td>0.01</td><td>0</td><td>0.00</td></tr><tr><td>Inductance ripple (%)</td><td>0.21</td><td>0.41</td><td>0.62</td><td>0.82</td><td>1.03</td></tr><tr><td>Converter efficiency</td><td>58.91</td><td>74.02</td><td>80.93</td><td>84.9</td><td>87.47</td></tr><tr><td>Diode junction temperature (°C)</td><td>211.72</td><td>140.76</td><td>119.64</td><td>109.56</td><td>103.66</td></tr><tr><td rowspan="5">dim4 - 70kW</td><td>Input current ripple (%)</td><td>1.19</td><td>1.14</td><td>1.09</td><td>1.04</td><td>0.99</td></tr><tr><td>Output voltage ripple (%)</td><td>0.12</td><td>0.06</td><td>0.04</td><td>0.03</td><td>0.02</td></tr><tr><td>Inductance ripple (%)</td><td>1.19</td><td>2.37</td><td>3.56</td><td>4.74</td><td>5.93</td></tr><tr><td>Converter efficiency</td><td>57.32</td><td>72.77</td><td>79.95</td><td>84.10</td><td>86.80</td></tr><tr><td>Diode junction temperature (°C)</td><td>227.18</td><td>148.56</td><td>124.97</td><td>113.67</td><td>107.04</td></tr><tr><td rowspan="5">dim3 - 60kW</td><td>Input current ripple (%)</td><td>1.73</td><td>1.62</td><td>1.51</td><td>1.40</td><td>1.29</td></tr><tr><td>Output voltage ripple (%)</td><td>0.19</td><td>0.09</td><td>0.06</td><td>0.05</td><td>0.04</td></tr><tr><td>Inductance ripple (%)</td><td>1.73</td><td>3.46</td><td>5.20</td><td>6.93</td><td>8.66</td></tr><tr><td>Converter efficiency</td><td>56.24</td><td>71.85</td><td>79.25</td><td>83.51</td><td>86.30</td></tr><tr><td>Diode junction temperature (°C)</td><td>237.95</td><td>153.96</td><td>128.67</td><td>116.52</td><td>109.39</td></tr><tr><td rowspan="5">dim2 - 50kW</td><td>Input current ripple (%)</td><td>4.40</td><td>3.10</td><td>1.81</td><td>0.52</td><td>0.59</td></tr><tr><td>Output voltage ripple (%)</td><td>0.70</td><td>0.35</td><td>0.23</td><td>0.18</td><td>0.14</td></tr><tr><td>Inductance ripple (%)</td><td>4.40</td><td>8.79</td><td>13.19</td><td>17.59</td><td>21.98</td></tr><tr><td>Converter efficiency</td><td>46.45</td><td>63.39</td><td>72.08</td><td>77.41</td><td>81.00</td></tr><tr><td>Diode junction temperature (°C)</td><td>366.57</td><td>218.71</td><td>172.69</td><td>150.28</td><td>137.04</td></tr><tr><td rowspan="5">dim1 - 40kW</td><td>Input current ripple (%)</td><td>4.63</td><td>1.43</td><td>1.12</td><td>1.11</td><td>0.16</td></tr><tr><td>Output voltage ripple (%)</td><td>1.26</td><td>0.63</td><td>0.42</td><td>0.32</td><td>0.25</td></tr><tr><td>Inductance ripple (%)</td><td>4.63</td><td>9.26</td><td>13.89</td><td>18.52</td><td>23.15</td></tr><tr><td>Converter efficiency</td><td>33.62</td><td>50.25</td><td>60.17</td><td>66.75</td><td>71.41</td></tr><tr><td>Diode junction temperature (°C)</td><td>682.36</td><td>377.33</td><td>279.92</td><td>232.01</td><td>203.51</td></tr><tr><td></td><td>Switch junction temperature (°C)</td><td>1498.71</td><td>570.84</td><td>361.22</td><td>275.09</td><td>229.41</td></tr></tbody></table>
<table><caption>**Table 4.** Sizing parameters.</caption><thead><tr><th>Parameter</th><th>Values</th></tr></thead><tbody><tr><td>L</td><td>0.001 (H)</td></tr><tr><td>C</td><td>0.001 (F)</td></tr><tr><td>freq</td><td>10,000 (Hz)</td></tr><tr><td>B</td><td>1.2 (T)</td></tr><tr><td>Diameter</td><td>1 · 10<sup>-4</sup> (m)</td></tr><tr><td>Winding number</td><td>300</td></tr><tr><td>Heat sink thermal resistance</td><td>0.5 (°C/W)</td></tr><tr><td>Number of FC cells</td><td>240</td></tr><tr><td>Active surface area</td><td>140 (m<sup>2</sup>)</td></tr></tbody></table>

Note that these parameters have been chosen based on the optimisation results for an operating point of 40 kW. The following sizings are designed by increasing either the number of cells or the active surface area. Hence, the results of all constraints with respect to the sizings are presented in Figure 13.

The above table presents the constraint values of each sizing in relation to the number of switching cells. The colour red highlights the constraints that stand well below their limit. On the contrary, the colour green emphasizes the respected constraints. Finally, the colour orange designates the constraints located around their boundary that are either above or below the limit.

According to the table, the 60 kW sizing with five switching cells gives the same results as the 80 kW sizing with three cells. Both have an efficiency that is close to the boundary, but the other constraints are properly respected. Consequently, to respect all of the constraints, there are two choices: either oversize the power chain, which leads to a serious volume increase, or add several switching cells to the converter, which causes a smaller volume augmentation as seen previously.

To simplify the reading of the previous figure, Figure 14 shows a schematic representation of the data collected. It can be read as follows: the sign “-” means that at least one constraint is not respected, and “+” consequently means that almost all constraints are respected. Note that the constraints here are the same as mentioned before. The two boxes with only one “+” sign are the sizings that respect all constraints, but some of them are located just above the limit. In contrast, the two boxes with “+” have some constraints that are just under their limits. For this reason, the red boxes show that the sizing chosen does not respect the imposed constraints, whereas the green ones are the sizings that fulfill all conditions.

According to these figures, choosing a sizing of 60 kW with five switching cells for the converter gives the same results as selecting a structure of 80 kW with three or four switching cells. In comparison, a sizing of 60 kW and five switching cells gives an overall volume of 0.1265 m³ or 126.5 L, whereas a sizing of 80 kW and three cells presents a volume of 0.1625 m³ or 162.5 L. This is due to the fact that the source occupies a larger volume than the converter, thus having a greater impact. These analyses show the benefits of the interleaving principle, which permits to downsize a power chain, leading to a smaller volume of overall architecture. Thanks to this interleaving approach, a single-cell design that is considered as downsized (as it does not respect the constraints) becomes suitable from a certain level of interleaving as the constraints become respected one-by-one. The results can be improved by considering other performance indexes, such as: consumption/efficiency, reliability, etc., in keeping with real expectations. Therefore, the choice of the number of switching cells should be carefully evaluated based on a trade-off analysis between these performance indexes.
<figure><img src="image_16.png" /><figcaption>**Figure 14.** Comparison of different sizings with different numbers of cells, simplified.</figcaption></figure>

## 5. Conclusions

This paper presents a pre-sizing approach that is intended for a hybrid battery-fuel cell system with interleaved converters. The objective of the proposed approach is to assist designers in obtaining an optimal sizing of power structures in accordance with the highly restrictive specifications. By introducing the interleaving technique, the method leads to a better flexibility for respecting the design constraints and enhances the performance. The approach is based on an optimisation process using a PSO algorithm under multi-physical constraints. A use case based on a hybrid power source system was considered to analyse the benefits of the interleaving concept. Analytical models were used to represent the behaviour of the studied system, including the sources and converters, as well as the associated constraints (volume, electric, thermal, and efficiency). These models are dependent on the number of switching cells of the converters, which is a key modular design parameter.

The results of the first validation show that the overall volume of the mono-source power chain decreases with the increase in the number of switching cells, reducing the volume by up to 15%. In the second validation on a hybrid architecture with a simple energy management system, the constraints on the architecture (such as electrical, thermal, and efficiency constraints) were found to be enhanced when the number of cells increases. This means that adding switching cells to an undersized architecture operating with only one converter cell can improve its constraints without significantly increasing the overall volume. In this example, an overall volume reduction of 22% was achieved. It is also worth noting that the sources (i.e., fuel cells and batteries) take up the most significant portion of the overall volume.

The limitations mentioned in the paper are important to consider for future work. A mono-objective optimisation may not always be sufficient for capturing all aspects of system performance, and a multi-objective optimisation could be used to find trade-offs between different performance indexes such as power efficiency and component lifespan to select the most suitable design solution based on the designer's preferences. Additionally, the reactive energy management approach used in this work may not be ideal for all applications, and a more proactive approach could be explored by incorporating fault tolerance considerations to enhance the performance and the reliability of the system.
**Author Contributions:** Conceptualisation, V.M., A.C., T.A., and C.L.; methodology, V.M., A.C., and T.A.; software, V.M. and A.C.; validation, V.M. and A.C.; writing—original draft preparation, V.M.; writing—review and editing, V.M., A.C., and T.A.; supervision, T.A. and C.L. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Conflicts of Interest:** The authors declare no conflict of interest.

**Abbreviations**

The following abbreviations are used in this manuscript:

<table><tr><td>EV</td><td>Electric vehicle</td></tr><tr><td>HEV</td><td>Hybrid electric vehicle</td></tr><tr><td>FCEV</td><td>Fuel cell electric vehicle</td></tr><tr><td>PEMFC</td><td>Proton exchange membrane fuel cell</td></tr><tr><td>FC</td><td>Fuel cell</td></tr><tr><td>ESS</td><td>Energy storage system</td></tr><tr><td>SoC</td><td>State of charge</td></tr><tr><td>SoH</td><td>State of health</td></tr><tr><td>DC</td><td>Direct current</td></tr><tr><td>PSO</td><td>Particle swarm optimisation</td></tr><tr><td>GA</td><td>Genetic algorithm</td></tr><tr><td>RMS</td><td>Root mean square</td></tr><tr><td>CCM</td><td>Continuous conduction mode</td></tr><tr><td>DCM</td><td>Discontinuous conduction mode</td></tr></table>

**References**

1. Ogunbode, C.A.; Doran, R.; Böhm, G. Exposure to the IPCC special report on 1.5 °C global warming is linked to perceived threat and increased concern about climate change. *Clim. Chang.* **2020**, *158*, 361–375. [CrossRef]

2. Chakraborty, S.; Kumar, N.M.; Jayakumar, A.; Dash, S.K.; Elangovan, D. Selected Aspects of Sustainable Mobility Reveals Implementable Approaches and Conceivable Actions. *Sustainability* **2021**, *13*, 12918. [CrossRef]

3. Muneer, T.; Kolhe, M.; Doyle, A. *Electric Vehicles: Prospects and Challenges*; Elsevier: Amsterdam, The Netherlands, 2017.

4. Zhang, F.; Hu, X.; Langari, R.; Cao, D. Energy management strategies of connected HEVs and PHEVs: Recent progress and outlook. *Prog. Energy Combust. Sci.* **2019**, *73*, 235–256. [CrossRef]

5. Jayakumar, A.; Chalmers, A.; Lie, T.T. Review of prospects for adoption of fuel cell electric vehicles in New Zealand. *IET Electr. Syst. Transp.* **2017**, *7*, 259–266. [CrossRef]

6. Ceschia, A.; Azib, T.; Bethoux, O.; Alves, F. Optimal design methodology for sizing a fuel cell/battery hybrid power source. *Proc. Inst. Mech. Eng. Part A: J. Power Energy* **2021**, *235*, 3–16. [CrossRef]

7. Zivanov, M.; Sasic, B.; Lazic, M. Desing of Multiphase Boost Converter for Hybrid Fuel Cell/Battery Power Sources. In *Paths to Sustainable Energy*; Nathwani, J., Ng, A., Eds.; IntechOpen: Rijeka, Croatia, 2010; Chapter 19. [CrossRef]

8. Saadi, R.; Hammoudi, M.; Kraa, O.; Ayad, M.; Bahri, M. A robust control of a 4-leg floating interleaved boost converter for fuel cell electric vehicle application. *Math. Comput. Simul.* **2020**, *167*, 32–47. [CrossRef]

9. Sorlei, I.S.; Bizon, N.; Thounthong, P.; Varlam, M.; Carcadea, E.; Culcer, M.; Iliescu, M.; Raceanu, M. Fuel cell electric vehicles—A brief review of current topologies and energy management strategies. *Energies* **2021**, *14*, 252. [CrossRef]

10. Samosir, A.S.; Mohd Yatim, A.H. Simulation and implementation of interleaved boost DC-DC converter for fuel cell application. *Int. J. Power Electron.* *Drive Syst.* **2011**, *1*, 168–174. [CrossRef]

11. Faraj, K.S.; Hussein, J. Analysis and comparison of DC-DC boost converter and interleaved DC-DC boost converter. *Eng. Technol. J.* **2020**, *38*, 622–635. [CrossRef]

12. Kascak, S.; Prazenica, M.; Jarabicova, M.; Paskala, M. Interleaved DC/DC boost converter with coupled inductors. *Adv. Electr. Electron. Eng.* **2018**, *16*, 147–154. [CrossRef]

13. Guilbert, D.; Gaillard, A.; N'Diaye, A.; Djerdir, A. Energy efficiency and fault tolerance comparison of DC/DC converters topologies for fuel cell electric vehicles. In Proceedings of the 2013 IEEE Transportation Electrification Conference and Expo (ITEC), Metro Detroit, MI, USA, 16–19 June 2013; pp. 1–7. [CrossRef]

14. Farhani, S.; N'Diaye, A.; Djerdir, A.; Bacha, F. Design and practical study of three phase interleaved boost converter for fuel cell electric vehicle. *J. Power Sources* **2020**, *479*, 228815. [CrossRef]
15. Younsi, M.O.; Bendali, M.; Azib, T.; Larouci, C.; Marchand, C.; Coquery, G. Current-sharing control technique of interleaved buck converter for automotive application. In Proceedings of the 7th IET International Conference on Power Electronics, Machines and Drives (PEMD 2014), Manchester, UK, 8–10 April 2014; pp. 1–6. [CrossRef]

16. Bendali, M.; Larouci, C.; Azib, T.; Marchand, C.; Coquery, G. Design methodology with optimization of an interleaved buck converter for automotive application. In Proceedings of the Eurocon 2013, Zagreb, Croatia, 1–4 July 2013; pp. 1066–1072. [CrossRef]

17. Pniak, L.; Revol, B.; Quéval, L.; Magambo, J.S.N.T.; Béthoux, O. Pre-sizing of a modular high power density DC/DC converter with GaN components. *Math. Comput. Simul.* **2023**. [CrossRef]

18. Delhommais, M.; Delaforge, T.; Schanen, J.L.; Wurtz, F.; Rigaud, C. A Predesign Methodology for Power Electronics Based on Optimization and Continuous Models: Application to an Interleaved Buck Converter. *Designs* **2022**, *6*, 68. [CrossRef]

19. Xun, Q.; Liu, Y.; Holmberg, E. A comparative study of fuel cell electric vehicles hybridization with battery or supercapacitor. In Proceedings of the 2018 International Symposium on Power Electronics, Electrical Drives, Automation and Motion (SPEEDAM), Amalfi, Italy, 20–22 June 2018; pp. 389–394. [CrossRef]

20. Pan, G.; Bai, Y.; Song, H.; Qu, Y.; Wang, Y.; Wang, X. Hydrogen Fuel Cell Power System—Development Perspectives for Hybrid Topologies. *Energies* **2023**, *16*, 2680. [CrossRef]

21. Azib, T.; Bethoux, O.; Remy, G.; Marchand, C.; Berthelot, É. An innovative control strategy of a single converter for hybrid fuel cell/supercapacitor power source. *IEEE Trans. Ind. Electron.* **2010**, *57*, 4024–4031. [CrossRef]

22. Thounthong, P.; Davat, B.; Rael, S.; Sethakul, P. Fuel cell high-power applications. *IEEE Ind. Electron. Mag.* **2009**, *3*, 32–46. [CrossRef]

23. Dicks, A.; Rand, D.A.J. *Fuel Cell Systems Explained*, 3rd ed.; Wiley: Hoboken, NJ, USA, 2018.

24. Mebarki, N.; Rekioua, T.; Mokrani, Z.; Rekioua, D.; Bacha, S. PEM fuel cell/battery storage system supplying electric vehicle. *Int. J. Hydrog. Energy* **2016**, *41*, 20993–21005. [CrossRef]

25. Ajayan, S.; Selvakumar, A.I. Modeling and simulation of PEM fuel cell electric vehicle with multiple power sources. *Int. J. Recent Technol. Eng.* **2020**, *8*, 2967–2975. [CrossRef]

26. Visvanathan, V.K.; Palaniswamy, K.; Ponnaiyan, D.; Chandran, M.; Kumaresan, T.; Ramasamy, J.; Sundaram, S. Fuel Cell Products for Sustainable Transportation and Stationary Power Generation: Review on Market Perspective. *Energies* **2023**, *16*, 2748. [CrossRef]

27. Saadi, A.; Becherif, M.; Aboubou, A.; Ayad, M. Comparison of proton exchange membrane fuel cell static models. *Renew. Energ.* **2013**, *56*, 64–71. [CrossRef]

28. Saldaña, G.; San Martín, J.I.; Zamora, I.; Asensio, F.J.; Oñederra, O. Analysis of the current electric battery models for electric vehicle simulation. *Energies* **2019**, *12*, 2750. [CrossRef]

29. Kabalo, M.; Blunier, B.; Bouquain, D.; Miraoui, A. State-of-the-art of DC-DC converters for fuel cell vehicles. In Proceedings of the 2010 IEEE Vehicle Power and Propulsion Conference, Lille, France, 1–3 September 2010; pp. 1–6. [CrossRef]

30. İnci, M.; Büyük, M.; Demir, M.H.; İlbey, G. A review and research on fuel cell electric vehicles: Topologies, power electronic converters, energy management methods, technical challenges, marketing and future aspects. *Renew. Sust. Energ. Rev.* **2021**, *137*, 110648. [CrossRef]

31. Ceschia, A.; Azib, T.; Bethoux, O.; Alves, F. Multi-Criteria Optimal Design for FUEL Cell Hybrid Power Sources. *Energies* **2022**, *15*, 3364. [CrossRef]

32. Niu, M.; Wan, C.; Xu, Z. A review on applications of heuristic optimization algorithms for optimal power flow in modern power systems. *J. Mod. Power Syst. Clean Energy* **2014**, *2*, 289–297. [CrossRef]

33. Rahmat-Samii, Y. Genetic algorithm (GA) and particle swarm optimization (PSO) in engineering electromagnetics. In Proceedings of the 17th International Conference on Applied Electromagnetics and Communications, Dubrovnik, Croatia, 1–3 October 2003; pp. 1–5. [CrossRef]

34. Ceschia, A.; Azib, T.; Bethoux, O.; Alves, F. Reliability Impact in Optimal Sizing of Fuel Cell/Battery Hybrid Power Sources. *Energies* **2020**, *13*, 3510. [CrossRef]

35. Martin, E.; Papageorgopoulos, D. DOE Hydrogen and Fuel Cells Program Record; Technical Report; US Department of Energy: Washington, DC, USA, 2015.

36. Raggl, K.; Nussbaumer, T.; Kolar, J.W. Model based optimization of EMC input filters. In Proceedings of the 2008 11th Workshop on Control and Modeling for Power Electronics, Zurich, Switzerland, 17–20 August 2008; pp. 1–6. [CrossRef]

37. Ejjabraoui, K.; Larouci, C.; Lefranc, P.; Marchand, C. A pre-sizing approach of DC-DC converters, application to design a Buck converter for the automotive domain. In Proceedings of the 2009 IEEE 6th International Power Electronics and Motion Control Conference, Wuhan, China, 17–20 May 2009; pp. 517–523. [CrossRef]

38. Lin, X.; Perez, H.E.; Mohan, S.; Siegel, J.B.; Stefanopoulou, A.G.; Ding, Y.; Castanier, M.P. A lumped-parameter electro-thermal model for cylindrical batteries. *J. Power Sources* **2014**, *257*, 1–11. [CrossRef]

39. Tran, M.K.; DaCosta, A.; Mevawalla, A.; Panchal, S.; Fowler, M. Comparative study of equivalent circuit models performance in four common lithium-ion batteries: LFP, NMC, LMO, NCA. *Batteries* **2021**, *7*, 51. [CrossRef]

40. Aitakkache, M.; Enrici, P.; Matt, D.; Boubaker, N.; Piscini, L. Concept, Feasibility of Cylindrical Bar Winding for Low Voltage Permanent Magnet Synchronous Motor. *Energies* **2022**, *15*, 1507. [CrossRef]
41. Liu, J.; Wilson, T.; Wong, R.; Wunderlich, R.; Lee, F. A method for inductor core loss estimation in power factor correction applications. In Proceedings of the APEC. Seventeenth Annual IEEE Applied Power Electronics Conference and Exposition (Cat. No.02CH37335), Dallas, TX, USA, 10–14 March 2002; pp. 439–445. [CrossRef]

42. Azib, T.; Larouci, C.; Chaibet, A.; Boukhnifer, M. Online energy management strategy of a hybrid fuel cell/battery/ultracapacitor vehicular power system. *IEEJ Trans. Electr. Electron. Eng.* **2014**, *9*, 548–554. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.