SOC Estimation of a Rechargeable Li-Ion Battery Used in Fuel-Cell Hybrid Electric Vehicles—Comparative Study of Accuracy and Robustness Performance Based on Statistical Criteria. Part I: Equivalent Models

Battery state of charge (SOC) accuracy plays a vital role in a hybrid electric vehicle (HEV), as it ensures battery safety in a harsh operating environment, prolongs life, lowers the cost of energy consumption, and improves driving mileage. Therefore, accurate SOC battery estimation is the central idea of the approach in this research, which is of great interest to readers and increases the value of its application. Moreover, an accurate SOC battery estimate relies on the accuracy of the battery model parameters and its capacity. Thus, the purpose of this paper is to design, implement and analyze the SOC estimation accuracy of two battery models, which capture the dynamics of a rechargeable SAFT Li-ion battery. The first is a resistor capacitor (RC) equivalent circuit model, and the second is a generic Simscape model. The model validation is based on the generation and evaluation of the SOC residual error. The SOC reference value required for the calculation of residual errors is the value estimated by an ADVISOR 3.2 simulator, one of the software tools most used in automotive applications. Both battery models are of real interest as a valuable support for SOC battery estimation by using three model based Kalman state estimators developed in Part 2. MATLAB simulations results prove the effectiveness of both models and reveal an excellent accuracy.


Literature Review
Nowadays the new technologies applied in batteries manufacturing industry "often demand more compact, higher capacity, safe and rechargeable batteries" [1]. The batteries vary by different chemistries and "generate the basic cell voltages typically in the 1.0 to 3.6 V range" [1]. The required voltages and the currents of a battery pack are obtained by adding up the number of the cells in a series

Linear Equivalent Electric Circuits and Simscape Battery Models
Finally, the last category of models is that of the linear equivalent electric circuit models (ECM), such as those discussed in the next section, Much details on these category of the models can be found in [21] for a Li-ion polymer (LiPb). The last category of models are Simscape models, described also in the next section of the present research work. The weakness and the strengths of ECM and Simscape models concerning the SOC accuracy and robustness are discussed in detail in next section.
In conclusion, this article focuses on the design and implementation of two accurate SAFT Li-ion battery models, suitable for HEV applications. For each Li-ion battery model are implemented in Part 2 three real-time SOC estimators on a MATLAB R2020a platform. The remaining sections of this paper are structured as follows. Section 2 describes the first RC ECM model attached to a SAFT Li-ion battery. Section 3 describes the second Li-ion battery model, a Simscape nonlinear model. Section 4 analyzes the SOC performance through six statistical criteria. Section 5 details the authors' contributions to this research paper.

Li-Ion RC Battery Equivalent Circuit Model-Case Study and ADVISOR Setup
The purpose of this section is to present the case study of a small urban hypothetic car (SMCAR) which is set up using the ADVISOR 3.2 version software package, one of the most used in the automotive industry. Then, in next subsections is developed and validate an accurate Li-ion battery model that describes the dynamics of a SAFT Li-ion battery with a rated capacity of 6Ah and a nominal voltage of 3.6 V. This model is a third-order RC equivalent circuit model (3RC ECM), one of the most used in HEV applications due to its simplicity, high accuracy, and fast real-time implementation [14,[17][18][19][20][21][22][23].

Li-Ion SAFT Battery and ADVISOR Small Hybrid Electric Car (SMCAR) Setup
SAFT is one of the most prestigious research companies in the US, among the most famous battery players on the commercial market in the world. It operates "under the auspices of the United States Advanced Battery Consortium (USABC) and the New Generation Vehicle Partnership (PNGV)," developing high-power lithium-ion (Li-ion) batteries over the past two decades. These batteries currently equip most HEVs and EVs [14][15][16][17][18][19]. The Li-ion battery together with other key components of a Hydrogen fuel cell electric vehicle are distributed on the car chassis as shown in Figure 1.
Batteries 2020, 6, x FOR PEER REVIEW 4 of 38 description for the behavior of a Li-Ion battery under pulsed discharge, for which the recovery effect is modeled as a decreasing exponential function of the SOC and discharge capacity" [29].

Linear Equivalent Electric Circuits and Simscape Battery Models
Finally, the last category of models is that of the linear equivalent electric circuit models (ECM), such as those discussed in the next section, Much details on these category of the models can be found in [21] for a Li-ion polymer (LiPb). The last category of models are Simscape models, described also in the next section of the present research work. The weakness and the strengths of ECM and Simscape models concerning the SOC accuracy and robustness are discussed in detail in next section.
In conclusion, this article focuses on the design and implementation of two accurate SAFT Liion battery models, suitable for HEV applications. For each Li-ion battery model are implemented in Part 2 three real-time SOC estimators on a MATLAB R2020a platform. The remaining sections of this paper are structured as follows. Section 2 describes the first RC ECM model attached to a SAFT Liion battery. Section 3 describes the second Li-ion battery model, a Simscape nonlinear model. Section 4 analyzes the SOC performance through six statistical criteria. Section 5 details the authors' contributions to this research paper.

Li-Ion RC Battery Equivalent Circuit Model-Case Study and ADVISOR Setup
The purpose of this section is to present the case study of a small urban hypothetic car (SMCAR) which is set up using the ADVISOR 3.2 version software package, one of the most used in the automotive industry. Then, in next subsections is developed and validate an accurate Li-ion battery model that describes the dynamics of a SAFT Li-ion battery with a rated capacity of 6Ah and a nominal voltage of 3.6 V. This model is a third-order RC equivalent circuit model (3RC ECM), one of the most used in HEV applications due to its simplicity, high accuracy, and fast real-time implementation [14,[17][18][19][20][21][22][23].

Li-Ion SAFT Battery and ADVISOR Small Hybrid Electric Car (SMCAR) Setup
SAFT is one of the most prestigious research companies in the US, among the most famous battery players on the commercial market in the world. It operates "under the auspices of the United States Advanced Battery Consortium (USABC) and the New Generation Vehicle Partnership (PNGV)," developing high-power lithium-ion (Li-ion) batteries over the past two decades. These batteries currently equip most HEVs and EVs [14][15][16][17][18][19]. The Li-ion battery together with other key components of a Hydrogen fuel cell electric vehicle are distributed on the car chassis as shown in Figure 1.   The key components of a Hydrogen fuel cell electric car shown in Figure 1 are described in [2] as follows: (1) The battery (auxiliary): this "powers vehicle accessories" and "provides electricity to start the car when traction battery is engaged". (2) Battery pack: this "stores energy generated from regenerative braking and provides supplemental power to the electric traction motor".
(3) Direct current-direct current (DC-DC) converter: this is an electronic device that "converts higher-voltage DC power from the traction battery pack to the lower-voltage DC power needed to run vehicle accessories and recharge the auxiliary battery". (4) Electric traction motor (FCEV): this is powered by the fuel cell and traction battery pack to drive the vehicle's wheels. It is also possible that "some of vehicles use motor generators that perform both the drive and regeneration functions". (5) Fuel-cell stack: it is "an assembly of individual membrane electrodes that use hydrogen and oxygen to produce electricity". (6) Fuel filter: this is "a nozzle from a high-pressure dispenser attaches to the receptacle on the vehicle to fill the tank". (7) Fuel hydrogen tank-it "stores the hydrogen gas on board the vehicle until it is required by the fuel cell". (8) Power electronics controller (FCEV): this is a unit that "manages the flow of electrical energy delivered by the fuel cell and the traction battery, controlling the speed of the electric traction motor and the torque it produces". (9) Thermal system (cooling) (FCEV): this "maintains a proper operating temperature range of the fuel cell, electric motor, power electronics, and other components". (10) Transmission (electric): this "transfers mechanical power from the electric traction motor to drive the wheels".
Among the Li-ion batteries of an HEV, the one with a capacity and a nominal voltage of 6 Ah and 3.6 V respectively is used for experimental validation tests, using an advanced simulator (ADVISOR) created in November 1994 by the US National Renewable Energy Laboratory (NREL). ADVISOR has so far proved to be the most suitable tools used in the design of HEV and EV systems, very well documented in [4][5][6][7]. Thanks to a wide variety of HEVs and EVs and the multitude of "real-world" driving conditions, it has gradually improved the performance until it reached version 2003-00, as well as the latest version r0116 of 24 April 2013, as mentioned in [14][15][16][17][18]. After proper installation, the ADVISOR graphical user interface (GUI) is running by typing "advisor" at the command prompt in MATLAB [14][15][16][17][18]. The ADVISOR GUI file menu has "help buttons which will either access the MATLAB help window or open a web page with appropriate context information" [15,16]. By using the ADVISOR GUI software package for design the following steps are requested: Step 1. Define a vehicle.
Step 1.1. Define the input HEV page setup shown in Figure 2, based on a large collection of HEVs types and characteristics contained by software. As a case study we consider a hypothetical SMCAR, powertrain control hybrid (hydrogen fuel cell electric vehicle) with the following characteristics [16]:  Tables 1 and 2 [14,15,[17][18][19]. Within a pack 84 cells are connected in series, with a nominal voltage capable of driving the PMDCM at approximately 300 V.
The Simulink block diagram of the transmission system and Li-ion battery storage is shown in Figure 3. As a case study we consider a hypothetical SMCAR, powertrain control hybrid (hydrogen fuel cell electric vehicle) with the following characteristics [16]:  Table 2. Li-ion mechanical and thermal characteristics SAFT battery of 6 Ah, 3.6 V (cylindrical shape) [16]. The Simulink block diagram of the transmission system and Li-ion battery storage is shown in Figure 3. As a case study we consider a hypothetical SMCAR, powertrain control hybrid (hydrogen fuel cell electric vehicle) with the following characteristics [16]: The Simulink block diagram of the transmission system and Li-ion battery storage is shown in Figure 3. Step 1.2. Drivetrain selection-selects the drivetrain configuration of the vehicle (Series, Parallel, etc.).
Step 1.5. Loading and saving vehicle configuration.
Step 2.1. Select the drive cycle-in the case study we chose the Federal Test Procedure (FTP) driving cycle used by US Environmental Protection Agency (EPA) for emissions certifications of passengers' vehicles in USA. The FTP-75 shown in Figure 4 and converted in current profile charging and discharging cycle in Figure 5 is the standard federal exhaust emissions driving cycle, which uses an Urban Dynamometer Driving Schedule (UDDS [14][15][16][17][18]). The FTP cycle has three separate phases: one cold-start phase (505 s), followed by a hot transient phase (870 s) and a hot-start phase (505 s) [14,16,18]. For a 10 min cool-down period between second phase and the third phase the engine is turned off. The first and third phase are identical. The total test time length for the FTP is 2457 s (40.95 min). The top speed is 91.25 km/h and the average speed is 25.82 km/h. The distance driven is approx. 17.7 km [14,16].
Step 2.5. Select multiple cycles to speed up the process of running many different cycles with the same initial conditions using this functionality.
Step 2.6. Choose a test procedure to select what kind of test to run.
Step 2.8. Run the simulation and wait for the results figure to popup, as shown in Figure 4 and Figure 5. Step 3. Looking for the simulation results. The first graph at the top of the Simulink configuration page, shown in Figures 4 and 5 (the first graph at the top), shows the ADVISOR FTP-75 driving cycle speed profile as the input variable. In Figure 5, the second graph at the top shows the estimated ADVISOR SOC value required to validate both models of Li-ion batteries attached to the SAFT Li-ion battery. The last chart at the bottom of Figure 5 shows the conversion of the FTP-75 driving cycle speed profile to a driving cycle current profile, required in MATLAB simulations for model validation and SOC estimators, as an input variable. The SAFT Li-ion battery electrical characteristics specifications are given in Table 1, and Table 2.  Step 2.2. Select a trip builder for a repeated cycle (if the case).
Step 2.3. Select a SOC correct options (linear or zero delta).
Step 2.4. Select interactive simulation a real-time interactive simulation interface to activate while the simulation is running.
Step 2.5. Select multiple cycles to speed up the process of running many different cycles with the same initial conditions using this functionality.
Step 2.6. Choose a test procedure to select what kind of test to run.
Step 2.8. Run the simulation and wait for the results figure to popup, as shown in Figures 4 and 5.
Step 3. Looking for the simulation results.
The first graph at the top of the Simulink configuration page, shown in Figures 4 and 5 (the first graph at the top), shows the ADVISOR FTP-75 driving cycle speed profile as the input variable. In Figure 5, the second graph at the top shows the estimated ADVISOR SOC value required to validate both models of Li-ion batteries attached to the SAFT Li-ion battery. The last chart at the bottom of Figure 5 shows the conversion of the FTP-75 driving cycle speed profile to a driving cycle current profile, required in MATLAB simulations for model validation and SOC estimators, as an input variable.
The SAFT Li-ion battery electrical characteristics specifications are given in Table 1, and Table 2.

Li-Ion-RC Equivalent Circuit Model
The 3RC ECM Li-ion battery model, as shown in Figure 6, consists of an OCV source in series with the internal Rin resistance of the battery and three parallel RC bias cells. RC cells are introduced into the circuit to capture the dynamic electrochemical behaviour of the battery and to increase the accuracy of the model. The first RC polarization cell captures the fast transient of the battery cell, and the last two RC cells capture the slow variation of the steady-state and increase the accuracy of the model. As the new technologies are largely dependent on batteries, it is important to develop accurate battery models that can be conveniently used with on-board power simulators and electronic on-board power systems, as mentioned in [14,[17][18][19][20][21][22][23].

Li-Ion-RC Equivalent Circuit Model
The 3RC ECM Li-ion battery model, as shown in Figure 6, consists of an OCV source in series with the internal Rin resistance of the battery and three parallel RC bias cells. RC cells are introduced into the circuit to capture the dynamic electrochemical behaviour of the battery and to increase the accuracy of the model. The first RC polarization cell captures the fast transient of the battery cell, and the last two RC cells capture the slow variation of the steady-state and increase the accuracy of the model. As the new technologies are largely dependent on batteries, it is important to develop accurate battery models that can be conveniently used with on-board power simulators and electronic onboard power systems, as mentioned in [14,[17][18][19][20][21][22][23].
For simulation purpose, a specific setup for the 3RC ECM Li-ion battery model parameters, shown in Table 1 or directly on the electrical scheme from Figure 1, is available to prove the effectiveness of the proposed SOC estimation strategies. This setup is achieved from a generic ECM by changing only the values of the model parameters in state-space equations.  [19]).
The Li-ion battery 3RC ECM model parameters are given in Table 3, and the OCV nonlinear model coefficients are shown in Table 4. Vs y=V(t) Figure 6. The third order 3RC electric circuit model (ECM)-Li-Ion battery represented in NI Multisim 14.1 editor (see [19]).
For simulation purpose, a specific setup for the 3RC ECM Li-ion battery model parameters, shown in Table 1 or directly on the electrical scheme from Figure 1, is available to prove the effectiveness of the proposed SOC estimation strategies. This setup is achieved from a generic ECM by changing only the values of the model parameters in state-space equations.
The Li-ion battery 3RC ECM model parameters are given in Table 3, and the OCV nonlinear model coefficients are shown in Table 4.  In a state-space representation, the continuous time 3RC ECM nonlinear model of SAFT Li-ion battery shown in Figure 6 is given by following Equations: where the components of the state vector are: x 4 = SOC is the state of charge of Li-ion battery, x 1 = V 1 is the voltage across first R 1 ||C 1 polarization cell, x 2 = V 2 denotes the voltage across the second R 2 ||C 2 polarization cell, x 3 = V 3 represents the third R 3 ||C 3 polarization cell u(t) = i(t) is the input discharging current (u(t) ≥ 0) or charging current (u(t) ≤ 0), OCV(t) represents the open-circuit voltage of Li-ion battery, and finally y(t) designates the terminal voltage of the battery. The open-circuit voltage of Li-ion battery OCV(t) given in (6) is a non-linear function of battery SOC, and contains a combination of the following three well-known generic battery models [17,[19][20][21]25]: (1) Shepherd model Batteries 2020, 6, 42 11 of 38 (2) Unnewehr universal model (3) Nernst model The performance of the generic models in terms of voltage prediction and SOC estimation is analysed in [24], and the simulations result show that the Unnewehr and Nernst models compared to the Shepherd model, criticized in literature, increase significantly the accuracy of linear ECMs, more specifically, Nernst model "showed the best performance among the three mathematical models" due to its flexibility by using two parameters (correction factors instead of one). Last, the combination of all three mathematical models in (7) and their introduction in the terminal voltage relationship (8) increases considerable the Li-ion ECM accuracy. Also, the ECM combined model proved until now that it is "amongst the most accurate formulations seen in literature from EVs/HVs field" [17,21].
Since the parameters of 3 RC ECM Li-ion model strongly depend on temperature and SOC, the combined model is beneficial due to its simplicity, accuracy, and development of BMS SOC estimators for HEVs as a "proof concept" and fast real-time implementation.
It is important to underline that the values of coefficients K 0 , K 1 , K 2 , K 3 , and K 4 , provided in Table 2, are chosen to fit the Li-ion battery model accurately according to the manufacturers' data by using a least squares curve fitting estimation method, as is suggested in [17][18][19][20][21][22][23]. The values of the resistances R 1 , R 2 and the capacitors C 1 , C 2 , as well the value of the battery nominal capacity C nom and its internal resistance are given in Table 1. The Simulink diagram of third order 3RC EMC-Li-Ion battery model that implements the Equations (3)-(8) is shown in Figure 7. combination of all three mathematical models in (7) and their introduction in the terminal voltage relationship (8) increases considerable the Li-ion ECM accuracy. Also, the ECM combined model proved until now that it is "amongst the most accurate formulations seen in literature from EVs/HVs field" [17,21]. Since the parameters of 3 RC ECM Li-ion model strongly depend on temperature and SOC, the combined model is beneficial due to its simplicity, accuracy, and development of BMS SOC estimators for HEVs as a "proof concept" and fast real-time implementation.
It is important to underline that the values of coefficients , , , , and , provided in Table 2, are chosen to fit the Li-ion battery model accurately according to the manufacturers' data by using a least squares curve fitting estimation method, as is suggested in [17][18][19][20][21][22][23]. The values of the resistances , and the capacitors , , as well the value of the battery nominal capacity and its internal resistance are given in Table 1. The Simulink diagram of third order 3RC EMC-Li-Ion battery model that implements the Equations (3)-(8) is shown in Figure 7.

RC Electric Circuit Model (ECM)-Li-Ion Battery Model-Discrete Time State Space Representation
For the design and implementation of SOC estimators it is necessary to discretize over time the continuous model of the Li-ion battery. The discrete model of 3RC ECM Li-ion model is described in a compact state space matrix representation, as:

RC Electric Circuit Model (ECM)-Li-Ion Battery Model-Discrete Time State Space Representation
For the design and implementation of SOC estimators it is necessary to discretize over time the continuous model of the Li-ion battery. The discrete model of 3RC ECM Li-ion model is described in a compact state space matrix representation, as: where is a vector of the input coefficients (input excitation terms), u(k) u(kT s ) , OCV(k) OCV(kT s ), and y(k) y(kT s ) denotes the input, respectively the OCV and the terminal output voltage at discrete time instants kT s | k∈Z + . In our MATLAB simulations, the sampling time is set to T s = 1 s, without any solver convergence problems.

Li-Ion Battery Thermal Model
Following the development in [14] the dynamics of Li-ion battery thermal model block is described by the following equations: where, the variables and the coefficients have the following significance and values: In MATLAB simulations, the battery temperature profile, and the robustness of the proposed SOC battery estimators are tested for the following approximate values, closed to a commercial battery type ICP 18,650 series [14]: An accurate simplified thermal model is provided in MATLAB R2019b library, at MATLAB/Simulink/Simscape/Battery, for a Li-ion generic battery model, implemented in Simulink Simscape as is shown in Figure 8, as is developed in [14]. The Equations (17) and (18)  An accurate simplified thermal model is provided in MATLAB R2019b library, at MATLAB/Simulink/Simscape/Battery, for a Li-ion generic battery model, implemented in Simulink Simscape as is shown in Figure 8, as is developed in [14]. The Equations (17) and (18) [14]).
Since the internal resistance of the Li-ion battery is the most sensitive to temperature developed inside the Li-ion battery, an overall Simulink model diagram block is designed that also integrates the Li-ion battery models such as is shown in Figure 9. It is essential to emphasize the fact that for performance comparison purpose, the overall Simulink model diagram shown in Figure 9, is sharing the same simplified thermal model to have an identical profile temperature and values of internal resistance ( ) and polarization ( ) . Figure 8. The detailed Simulink diagram of the Simulink Simscape thermal model block (see [14]).
Since the internal resistance of the Li-ion battery is the most sensitive to temperature developed inside the Li-ion battery, an overall Simulink model diagram block is designed that also integrates the Li-ion battery models such as is shown in Figure 9. It is essential to emphasize the fact that for performance comparison purpose, the overall Simulink model diagram shown in Figure 9, is sharing the same simplified thermal model to have an identical profile temperature and values of internal resistance R in (T) and polarization K p (T) .
inside the Li-ion battery, an overall Simulink model diagram block is designed that also integrates the Li-ion battery models such as is shown in Figure 9. It is essential to emphasize the fact that for performance comparison purpose, the overall Simulink model diagram shown in Figure 9, is sharing the same simplified thermal model to have an identical profile temperature and values of internal resistance ( ) and polarization ( ) .  It is important to remark that in Figure 9 the second block from the top of the Simulink diagram is introduced only to investigate the SOC and temperature profile evolutions delivered by the first block from the bottom side in comparison to SOC and the temperature profile delivered by the top side block. The ambient temperature profile and the output temperature of the Simulink Simscape thermal model described by Equations (15) It is important to remark that in Figure 9 the second block from the top of the Simulink diagram is introduced only to investigate the SOC and temperature profile evolutions delivered by the first block from the bottom side in comparison to SOC and the temperature profile delivered by the top side block. The ambient temperature profile and the output temperature of the Simulink Simscape thermal model described by Equations (15)- (18) are shown in Figure 10a The evolution of the battery internal resistance ( ) and polarization constant ( ) at room temperature t0 = 20 °C, is shown in Figure 11a,b.
(a) (b) Figure 10. (a) The input ambient temperature profile; (b) the output temperature as response to input ambient temperature of the thermal model block.
The evolution of the battery internal resistance R in (T) and polarization constant K p (T) at room temperature t 0 = 20 • C, is shown in Figure 11a  The output temperature profile of the Simulink Simscape thermal model for changes in ambient temperature is shown in Figure 12a, and the effects on internal battery resistance ( ) and polarization ( ) are presented in Figure 12b,c. The output temperature profile of the Simulink Simscape thermal model for changes in ambient temperature is shown in Figure 12a, and the effects on internal battery resistance R in (T) and polarization K p (T) are presented in Figure 12b,c.    Figure 13a,b are depicted the FTP current profile test (a), and the value of i-ion model SOC versus ADVISOR SOC estimate to the FTP current profile test (b) NREL ADVISOR MATLAB platform as shown in Figure 5. Using a MATLAB magnification tool on two portions of the graph, the visibility of both cu shown in Figure 13b increases considerably [26]. The Li-ion battery terminal voltage for an FTP current charging and discharging profile test, the following three, namely OCV = f(SOC) cu battery terminal voltage and its SOC, all of these three simulated for a constant discharge curren 1C-rate (6A), it can see in the Figure 14a-d.
The simulation results from last three Figure 14b-d reveal that all three battery characteris are quite close to the manufacturing specifications. It should also be noted that the OCV = f(S curve from Figure 14b is almost flat on a large portion. Therefore, the Coulomb counting metho not accurate for direct SOC measurement for Li-ion batteries. Thus, its estimation is necessary u one of the best known Kalman filtration techniques.  Using a MATLAB magnification tool on two portions of the graph, the visibility of both curves shown in Figure 13b increases considerably [26]. The Li-ion battery terminal voltage for an FTP-75 current charging and discharging profile test, the following three, namely OCV = f(SOC) curve, battery terminal voltage and its SOC, all of these three simulated for a constant discharge current at 1C-rate (6A), it can see in the Figure 14a-d.
The simulation results from last three Figure 14b-d reveal that all three battery characteristics are quite close to the manufacturing specifications. It should also be noted that the OCV = f(SOC) curve from Figure 14b is almost flat on a large portion. Therefore, the Coulomb counting method is not accurate for direct SOC measurement for Li-ion batteries. Thus, its estimation is necessary using one of the best known Kalman filtration techniques.  Using a MATLAB magnification tool on two portions of the graph, the visibility of both curves shown in Figure 13b increases considerably [26]. The Li-ion battery terminal voltage for an FTP-75 current charging and discharging profile test, the following three, namely OCV = f(SOC) curve, battery terminal voltage and its SOC, all of these three simulated for a constant discharge current at 1C-rate (6A), it can see in the Figure 14a-d.
The simulation results from last three Figure 14b-d reveal that all three battery characteristics are quite close to the manufacturing specifications. It should also be noted that the OCV = f(SOC) curve from Figure 14b is almost flat on a large portion. Therefore, the Coulomb counting method is not accurate for direct SOC measurement for Li-ion batteries. Thus, its estimation is necessary using one of the best known Kalman filtration techniques.

RC ECM Li-Ion Battery Model-MATLAB Model Validation
For validation of 3RC ECM Li-ion battery model, in the first stage is calculated the model SOC residue as a difference between the SOC values of the 3RC ECM model and the estimated values of ADVISOR SOC estimator. The SOC accuracy performance of 3RC ECM Li-ion battery model is analyzed by evaluating the SOC residual error. The residual percentage error is depictured in Figure 15.  Figure  15. The results of the MATLAB simulations shown in Figure 13b for a current FTP-75 driving cycle profile test reveal excellent SOC accuracy of the 3RC ECM compared to the estimated SOC value obtained by the ADVISOR simulator. From a quantitative point of view, this confirms the information extracted by evaluating the SOC residues generated in MATLAB and presented in Figure 15. From Figure 15 it can be seen that the SOC residue is in the range [−1.1, 0, 4], and the SOC error rate is less than 1.2%, which is an excellent result, comparable to those reported in the literature, even better. This result reveals that the 3RC ECM Li-ion battery model is very accurate in terms of SOC calculation, and the model is undoubtedly validated based on available information about its behaviour.

Li-Ion Battery Simscape Generic Model
A full representation of the generic battery model, dependent on the temperature and ageing effects, is developed by MathWorks team, as shown in the MATLAB R2019b/Simulink/Simscape/Power Systems/Extra Sources Library-Documentation.

Li-Ion Battery Simscape Generic Model-Description and Parameters' Specifications
Li-ion battery cell specifications for a Simscape model are shown in Figure 16a-c. The Li-ion battery Simscape model is more realistic and suitable to operate safely in different conditions. Also, this model is beneficial for an appropriate choice of battery chemistry and for different parameters specifications. The Simscape generic model developed by MathWorks team takes into consideration the thermal model of the battery (internal and environmental temperatures) and its ageing effects. The battery terminal voltage, current and SOC "can be visualized to monitor and control the battery SOH condition" [14].
The nominal current discharge characteristics according to a choice of a Li-ion battery which has a rated capacity of 6 Ah and a nominal voltage of 3.6 V for different x-scales (time, Ah) are shown in the Figures 17 and 18. The results of the MATLAB simulations shown in Figure 13b for a current FTP-75 driving cycle profile test reveal excellent SOC accuracy of the 3RC ECM compared to the estimated SOC value obtained by the ADVISOR simulator. From a quantitative point of view, this confirms the information extracted by evaluating the SOC residues generated in MATLAB and presented in Figure 15. From Figure 15 it can be seen that the SOC residue is in the range [−1.1, 0, 4], and the SOC error rate is less than 1.2%, which is an excellent result, comparable to those reported in the literature, even better. This result reveals that the 3RC ECM Li-ion battery model is very accurate in terms of SOC calculation, and the model is undoubtedly validated based on available information about its behaviour.

Li-Ion Battery Simscape Generic Model
A full representation of the generic battery model, dependent on the temperature and ageing effects, is developed by MathWorks team, as shown in the MATLAB R2019b/Simulink/Simscape/Power Systems/Extra Sources Library-Documentation.

Li-Ion Battery Simscape Generic Model-Description and Parameters' Specifications
Li-ion battery cell specifications for a Simscape model are shown in Figure 16a-c. The Li-ion battery Simscape model is more realistic and suitable to operate safely in different conditions. Also, this model is beneficial for an appropriate choice of battery chemistry and for different parameters specifications. The Simscape generic model developed by MathWorks team takes into consideration the thermal model of the battery (internal and environmental temperatures) and its ageing effects. The battery terminal voltage, current and SOC "can be visualized to monitor and control the battery SOH condition" [14].       The Simscape model of a generic 6 Ah and 3.6 V Li-ion battery SAFT-type without temperature and ageing effects is shown in Figure 19, the same shown in [14], p. 12. The Simscape model of a generic 6 Ah and 3.6 V Li-ion battery SAFT-type without temperature and ageing effects is shown in Figure 19, the same shown in [14], p. 12. Figure 19. The Simscape model of a generic 6 Ah and 3.2 V Li-ion battery (without temperature and aging effects (see [14], p. 12) connected to FTP-75 input current profile. A significant advantage of the Simscape model of Li-ion battery is the simplicity with which the model parameters for different chemistry and specifications are extracted as if we had access to the specifications of the battery manufacturers. The parameters of Li-ion battery choice extracted from the discharge characteristics shown in Figure 17 or Figure 18 have the following values:

Li-Ion Battery Simscape Model-Discrete Time in State Space Representation
The Simscape model parameters suggested in Figures 16-18 fit the following adopted model represented in discrete time in a unidimensional state space, like the model developed in [14], p. 21:

Li-Ion Battery Simscape Model-Discrete Time in State Space Representation
The Simscape model parameters suggested in Figures 16-18 fit the following adopted model represented in discrete time in a unidimensional state space, like the model developed in [14], p. 21: where x 1 (k) x 1 (kT s ) = SOC(kT s ) , u(k), y(k) , Q nom , η and T s have the same meaning as the variables and parameters that describe the 3RC ECM Li-ion battery model given by Equations (12)- (14). It is essential to emphasize a great advantage of the adopted Simulink Simscape model, presented in (21) and (22), consisting of a considerable model simplification and dependence only on SOC. Also, the dynamics of this model is described by the first Equation (21) which is linear and the second Equation (22) is a highly nonlinear static representation. The Simulink Simscape model of Li-ion battery that implements Equations (21) and (22) is shown in Figure 20.

Li-Ion Battery Simscape Generic Model-MATLAB Simulations Results and Model Validation
The MATLAB simulations result is shown in Figure 21a,c. In Figure 21a, the simulation result reveals an excellent SOC accuracy of the Simscape model of the Li-ion battery. This result is also supported by a small SOC residue, recorded in Figure 21b, which falls in the range [−1.4, 1]. Like the 3RC ECM Li-ion battery model, the Simscape model of the Li-ion battery based on the available information extracted from Figure 21a,b also works very well, because the SOC error percentage is less than 1.4%, compared to the typical value of 2% reported in the literature for similar applications. These results also validate this model, which is suitable to use it in the second part for real-time design and implementation on an attractive MATLAB 2020Ra environment.

Simulink Simscape Graphic Models Integrated in Fuel Cell HEV Applications-Energy Management System
This section presents some HEV applications that operate with graphic Simscape models. In this description, the Simscape "blocks language" allows much faster models of physical systems to be created within the Simulink environment, "based on physical connections that directly integrate with

Li-Ion Battery Simscape Generic Model-MATLAB Simulations Results and Model Validation
The MATLAB simulations result is shown in Figure 21a,c.

Li-Ion Battery Simscape Generic Model-MATLAB Simulations Results and Model Validation
The MATLAB simulations result is shown in Figure 21a,c. In Figure 21a, the simulation result reveals an excellent SOC accuracy of the Simscape model of the Li-ion battery. This result is also supported by a small SOC residue, recorded in Figure 21b, which falls in the range [−1.4, 1]. Like the 3RC ECM Li-ion battery model, the Simscape model of the Li-ion battery based on the available information extracted from Figure 21a,b also works very well, because the SOC error percentage is less than 1.4%, compared to the typical value of 2% reported in the literature for similar applications. These results also validate this model, which is suitable to use it in the second part for real-time design and implementation on an attractive MATLAB 2020Ra environment.

Simulink Simscape Graphic Models Integrated in Fuel Cell HEV Applications-Energy Management System
This section presents some HEV applications that operate with graphic Simscape models. In this description, the Simscape "blocks language" allows much faster models of physical systems to be created within the Simulink environment, "based on physical connections that directly integrate with In Figure 21a, the simulation result reveals an excellent SOC accuracy of the Simscape model of the Li-ion battery. This result is also supported by a small SOC residue, recorded in Figure 21b, which falls in the range [−1.4, 1]. Like the 3RC ECM Li-ion battery model, the Simscape model of the Li-ion battery based on the available information extracted from Figure 21a,b also works very well, because the SOC error percentage is less than 1.4%, compared to the typical value of 2% reported in the literature for similar applications. These results also validate this model, which is suitable to use it in the second part for real-time design and implementation on an attractive MATLAB 2020Ra environment.

Simulink Simscape Graphic Models Integrated in Fuel Cell HEV Applications-Energy Management System
This section presents some HEV applications that operate with graphic Simscape models. In this description, the Simscape "blocks language" allows much faster models of physical systems to be created within the Simulink environment, "based on physical connections that directly integrate with block diagrams and other modeling paradigms" [8]. In Simscape, the models can be parametrized using MATLAB variables and expressions and can be designed and implemented control systems for any physical system in Simulink". Users can easily integrate physical object icons into the design of Simulink diagrams or combine object models with the symbols of different physical objects. Indeed, behind each image is encapsulated the dynamic pattern of physical objects. However, it significantly eliminates the user's effort to write a lot of equations for modelling the dynamics of objects, which takes a long time, and the diagrams become much more complicated [8].

Hybrid Energy Storage of Energy Management System (EMS)-Simulink and Simscape Components Description
The hybrid energy storage (HES) of an EMS, shown in Figure 22, is a hybrid combination of three power sources, such as a fuel cell, Li-ion battery, and supercapacitor [5,7]. The control strategy of HES is implemented in Simulink Simscape to "manage the energy consumption of the hydrogen fuel, and at the same time the pulsed or transient power required (load profile) by the load should be supplied." [11]. To simplify the Simulink diagram of the EMS, are used Simscape components such as Li-ion battery, supercapacitor and FCPM that also encapsulates a hydrogen fuel stack cell Simscape model, provided by MATLAB Simulink Toolbox/Simscape. In this section, a brief presentation of this topic is given, since is only emphasized the fact that using a single power source such as Li-ion battery in driving HEV powertrains applications "has certain disadvantages such as recharging, longevity, poor power density, etc." [24]. A three-phase load profile is "emulated to consider variations in power at the different timings and simulations to see the behavior of the hybrid energy storage system (HESS) as a whole and the response of each storage system" [11]. Also, a Simscape "15 kW protecting resistor is integrated in the Simulink diagram to avoid overcharging the supercapacitor and battery systems" [7].  [5,7]).
In this diagram other Simscape components are integrated such as three DC-DC boost/buck converters blocks to interface with all three sources. The first is a 12.5 kW fuel cell DC/DC boost converter, with regulated output voltage and input current limitation, and the other are two DC/DC converters for discharging (4 kW boost converter) and respectively for charging (1.2 kW buck converter) the battery system. Normally, a "single bidirectional DC/DC converter can also be used to reduce the weight of the power system" [7]. The FC is controlled by a DC-DC boost converter, an electronic device controlled also by a signal sent by one of five control strategies conceived for this purpose inside the EMS block, as is described in [5]. Similarly, the Li-ion battery and the supercapacitor are controlled by a bidirectional DC-DC buck-boost converters, since during operation they are charging and discharging. The topology configuration and the electronic circuits are well described in [11]. The charging and the discharging cycles of the bidirectional DC-DC converters are controlled by a voltage signal provided by EMS block that adjusts the duty cycles (D) of both DC-DC converters, based on the following relationship [11]: (22) where D is the duty cycle, V out designates the output voltage of the converter, and V in denotes the input voltage. In this section is presented only briefly the most relevant MATLAB simulation results for EMS techniques to have a better insight of the behavior of all three Simscape components of hybrid power sources, i.e., FC, Li-ion battery, and SC (UC). In the Simulink Simscape diagram of fuel-cell hybrid power generation (FCHPG) shown in Figure 22, the inverter DC/AC that supplies the load is rated at 270 V DC in input, and 200 V AC, 400 Hz, 15 kVA in output. A three-phase load profile is "emulated to consider variations in power at the different timings and simulations to see the behavior of the hybrid energy storage system (HESS) as a whole and the response of each storage system" [11]. Also, a Simscape "15 kW protecting resistor is integrated in the Simulink diagram to avoid overcharging the supercapacitor and battery systems" [7].

Hybrid Energy Storage of EMS-Simulink Simscape Applications
As a practical application, the following three scenarios are implemented to reveal the behavior of HESS components: • Scenario 1: DC grid interfaces only the AC Grid and AC load, such in [7,11].
The HESS distributes the power among the energy sources according to a given energy management strategy. The MATLAB simulation results of EMS techniques are shown only for three setups, such as the state machine control strategy (SMCS), classical PI control strategy (PICS), and the equivalent consumption minimization strategy (ECMS) [5,7,25]. To obtain a sound theoretical background on the EMS design and implementation in a real-time MATLAB simulation environment, the following sources [5,7,11,24] provide valuable information. For EMS-SMCS setup shown in Figure 23, the MATLAB simulations result is presented in the Figure 24.
The HESS distributes the power among the energy sources according to a given energy management strategy. The MATLAB simulation results of EMS techniques are shown only for three setups, such as the state machine control strategy (SMCS), classical PI control strategy (PICS), and the equivalent consumption minimization strategy (ECMS) [5,7,25]. To obtain a sound theoretical background on the EMS design and implementation in a real-time MATLAB simulation environment, the following sources [5,7,11,24] provide valuable information. For EMS-SMCS setup shown in Figure 23, the MATLAB simulations result is presented in the Figure 24. In Figure 24 is depictured the powers' distribution for FC, Li-ion battery, UC, and Load profile for a hypothetical SMCAR HEV case study. The balance equation is given by: where is the power provided by FC, is the power delivered by Li-ion battery, is the power delivered by the UC to manage power peaks for vehicle acceleration and regeneration, and is the load profile (demand, total power required). From Figure 24 it is straightforward to check that Equation (23) is satisfied for each time moment. Also, it is obvious that for load power profile pecks the power delivered by UC is very sharp to cover the power demanded ( ) .  In Figure 24 is depictured the powers' distribution for FC, Li-ion battery, UC, and Load profile for a hypothetical SMCAR HEV case study. The balance equation is given by: where P FC is the power provided by FC, P Batt is the power delivered by Li-ion battery, P UC is the power delivered by the UC to manage power peaks for vehicle acceleration and regeneration, and P L is the load profile (demand, total power required). From Figure 24 it is straightforward to check that Equation (23) is satisfied for each time moment. Also, it is obvious that for load power profile pecks the power delivered by UC is very sharp to cover the power demanded (P L ). In the Figure 25a  Matlab simulation results reveal an evolution with sharp peaks for UC current and UC voltage when Li-ion battery needs to provide much more power to the DC network or during sudden acceleration and regeneration.
In the Figure 27a-c are presented the Li-ion battery current (a), battery terminal voltage (b) and battery SOC (c) respectively, according to power required. In the Figure 26a,b are depicted the UC current (a) and UC voltage variation (b) respectively, according to the load profile. Matlab simulation results reveal an evolution with sharp peaks for UC current and UC voltage when Li-ion battery needs to provide much more power to the DC network or during sudden acceleration and regeneration.
In the Figure 27a-c are presented the Li-ion battery current (a), battery terminal voltage (b) and battery SOC (c) respectively, according to power required. Matlab simulation results reveal an evolution with sharp peaks for UC current and UC voltage when Li-ion battery needs to provide much more power to the DC network or during sudden acceleration and regeneration.
In the Figure 27a-c are presented the Li-ion battery current (a), battery terminal voltage (b) and battery SOC (c) respectively, according to power required. Matlab simulation results reveal an evolution with sharp peaks for UC current and UC voltage when Li-ion battery needs to provide much more power to the DC network or during sudden acceleration and regeneration.
In the Figure 27a-c are presented the Li-ion battery current (a), battery terminal voltage (b) and battery SOC (c) respectively, according to power required.  Figure 27a,b show several peaks (positive and negative) in the evolution of the current of the Liion battery, which correspond to the charging and discharging cycles of the battery, as can be seen from the evolution of SOC in Figure 27c. The terminal voltage of the Li-ion battery also has a lot of variations in its growth, decreasing from 83 V to t = 0 V to 60 V around t = 125 [s], followed by an increase to 83V when t = 83V.MATLAB simulation results analysis for Scenario 1 In this application, it is important to analyze the power distribution shown in Figure 24. The result of the analysis provides a better perspective on how EMS works in real-time simulations. In this figure, the red colour curve represents the profile of the power load, i.e., a variable power required for the AC load in the first 350 s of real-time simulation. The blue colour curve designates the main power generated by FC source, which is the dominant source, i.e., the one that delivers the most considerable amount of power to a DC grid and is almost constant inside the 150 s window length (75,225) [s]. The brown curve refers to the second power supply source, which is a Li-ion battery that delivers power to the DC grid in a smaller and variable amount during charging and discharging cycles, compared to FC. Finally, the green colour curve refers to the third power supply source that delivers the smallest amount of power to the DC grid only during the short periods of sharp acceleration and regeneration. The power distribution balance can be easily checked for enough moments because the MATLAB Data Tips measurement tool can help to mark several points on each curve. For example, at time t = 70 s, the power delivered by each source of power supply has the following values: PFC = 6.985 kW, PBatt = 3.041 kW, PUC = 0 kW, and PL = 10.07kW The power distribution evaluated at t = 70 s verifies with enough accuracy Equation (23), since:  Figure 27a,b show several peaks (positive and negative) in the evolution of the current of the Li-ion battery, which correspond to the charging and discharging cycles of the battery, as can be seen from the evolution of SOC in Figure 27c. The terminal voltage of the Li-ion battery also has a lot of variations in its growth, decreasing from 83 V to t = 0 V to 60 V around t = 125 [s], followed by an increase to 83V when t = 83V.MATLAB simulation results analysis for Scenario 1 In this application, it is important to analyze the power distribution shown in Figure 24. The result of the analysis provides a better perspective on how EMS works in real-time simulations. In this figure, the red colour curve represents the profile of the power load, i.e., a variable power required for the AC load in the first 350 s of real-time simulation. The blue colour curve designates the main power generated by FC source, which is the dominant source, i.e., the one that delivers the most considerable amount of power to a DC grid and is almost constant inside the 150 s window length (75,225) [s]. The brown curve refers to the second power supply source, which is a Li-ion battery that delivers power to the DC grid in a smaller and variable amount during charging and discharging cycles, compared to FC. Finally, the green colour curve refers to the third power supply source that delivers the smallest amount of power to the DC grid only during the short periods of sharp acceleration and regeneration. The power distribution balance can be easily checked for enough moments because the MATLAB Data Tips measurement tool can help to mark several points on each curve. For example, at time t = 70 s, the power delivered by each source of power supply has the following values: P FC = 6.985 kW, P Batt = 3.041 kW, P UC = 0 kW, and P L = 10.07kW The power distribution evaluated at t = 70 s verifies with enough accuracy Equation (23), since: P FC + P Batt + P UC = 10.026 kW, so close to P L = 10.07 kW As Equation (23) is satisfied for each moment, it is easy to observe the behavior of all three power supply sources. The MATLAB simulation results shown in Figures 25c, 26a and 27a reveal the same trend of current evolution as that of each corresponding power supply.
As in the case of the EMS-SMCS setup, similar graphs with the same meaning are presented in Appendix A, Figures A1-A10 for second EMS-PICS setup, and in the ( [30], Figures A11-A20) for third EMS-ECMS setup. The theory and all the Simulink diagrams behind the five EMS techniques are fully documented in [5]. In our research, these EMS techniques are presented only as complementary information for interested readers, such as to give a clue, motivation and to open new research directions for future HEV developments. Nonetheless, the topic is beyond the scope of this paper, which is focused only on the modeling aspects and Li-ion battery SOC estimation techniques.

•
Scenario 2: A 100 HP, 1750 RPM asynchronous induction motor (squirrel cage) is connected to AC grid as is shown in Figure 28.
Batteries 2020, 6, x FOR PEER REVIEW 26 of 38 • Scenario 2: A 100 HP, 1750 RPM asynchronous induction motor (squirrel cage) is connected to AC grid as is shown in Figure 28. In this scenario is shown only the MATLAB simulations result related to the evolution of asynchronous induction motor speed, as can be seen in Figure 29.
From Figure 29, it is easy to observe that the speed of induction motor connected to AC grid is not controlled, but it is quite close to 1750 RPM in steady state, for a torque load of 375 Nm. In this scenario is shown only the MATLAB simulations result related to the evolution of asynchronous induction motor speed, as can be seen in Figure 29. In this scenario is shown only the MATLAB simulations result related to the evolution of asynchronous induction motor speed, as can be seen in Figure 29.
From Figure 29, it is easy to observe that the speed of induction motor connected to AC grid is not controlled, but it is quite close to 1750 RPM in steady state, for a torque load of 375 Nm.

•
Scenario 3: A 2 HP 1750 permanent magnet DC motor connected to DC side of the grid, as is shown in Figure 22 (the right topside block).
Since the block from Figure 22 encapsulates the PMDCM Simscape model, for clarity, Figure 30 shows the PMDCM Simscape model with all the details of electrical connections. From Figure 29, it is easy to observe that the speed of induction motor connected to AC grid is not controlled, but it is quite close to 1750 RPM in steady state, for a torque load of 375 Nm.

•
Scenario 3: A 2 HP 1750 permanent magnet DC motor connected to DC side of the grid, as is shown in Figure 22 (the right topside block).
Since the block from Figure 22 encapsulates the PMDCM Simscape model, for clarity, Figure 30 shows the PMDCM Simscape model with all the details of electrical connections.
Unlike the uncontrolled speed of the asynchronous induction motor (ASM), the PMDCM is connected by a negative feedback in a closed-loop to a block of proportional-integral-derivative (PID) controller, to control its speed, as is shown in Figure 31, and developed in [12]. The MATLAB PMDCM RPM speed step response is shown in Figure 32. A big advantage of the PID controller is that the PMDCM speed response converges quickly and reaches the target speed of 1500 RPM in almost 1.8 s. Figure 30. A 2HP 1750 RPM PMDCM-Simscape model (see [12]).
Unlike the uncontrolled speed of the asynchronous induction motor (ASM), the PMDCM is connected by a negative feedback in a closed-loop to a block of proportional-integral-derivative (PID) controller, to control its speed, as is shown in Figure 31, and developed in [12].
Unlike the uncontrolled speed of the asynchronous induction motor (ASM), the PMDCM is connected by a negative feedback in a closed-loop to a block of proportional-integral-derivative (PID) controller, to control its speed, as is shown in Figure 31, and developed in [12]. The MATLAB PMDCM RPM speed step response is shown in Figure 32. A big advantage of the PID controller is that the PMDCM speed response converges quickly and reaches the target speed of 1500 RPM in almost 1.8 s. The MATLAB PMDCM RPM speed step response is shown in Figure 32. A big advantage of the PID controller is that the PMDCM speed response converges quickly and reaches the target speed of 1500 RPM in almost 1.8 s. The MATLAB PMDCM RPM speed step response is shown in Figure 32. A big advantage of the PID controller is that the PMDCM speed response converges quickly and reaches the target speed of 1500 RPM in almost 1.8 s.  The behaviour of the SC connected to the DC grid to supply power to the PMDCM during The behaviour of the SC connected to the DC grid to supply power to the PMDCM during sudden changes in the load torque is described in Figure 34a for the SC current and in Figure 34b for   The PMDCM behavior during the operation is shown in Figure 35a-d.  The PMDCM behavior during the operation is shown in Figure 35a-d.  Figure 35a shows the linear evolution of the load torque with its PMDCM speed (scaling factor 0.011). Graph 35b describes the input power that changes abruptly in the step moment from 0 to 7 kW at t = 2 [s], followed by a decrease to 2.5 kW during steady state. Figure 35c shows the armature current absorbed by the PMDCM with a similar evolution trend as for the absorbed PMDCM power.  Figure 35a shows the linear evolution of the load torque with its PMDCM speed (scaling factor 0.011). Graph 35b describes the input power that changes abruptly in the step moment from 0 to 7 kW at t = 2 [s], followed by a decrease to 2.5 kW during steady state. Figure 35c shows the armature current absorbed by the PMDCM with a similar evolution trend as for the absorbed PMDCM power. Figure 35d shows the supplied armature voltage which is the same as the SC voltage, which justifies its presence to provide a large amount of energy, again protecting the Li-ion battery to make this effort. SC ensures the required voltage absorbed by PMDCM to achieve excellent speed profile tracking performance.

Statistical Criteria to Asses the Accuracy of the Models
In a general formulation, for a better understanding of how to select an accurate Li-ion battery model, as well as a high-precision SOC state estimator and an excellent prediction of the selected battery output voltage, it can use some statistical criteria performance to compute the fitting errors between a set of candidates models reported in the literature [14,17,26,27]. Selection of Li-ion battery models and Kalman filter SOC estimators can be made by using the performance criteria developed in the "recent years in statistical learning, machine learning, and big data analytics" [27]. It is essential to emphasize the fact that now there are several criteria reported in the literature for models and estimators' selection, that "receives much attention due to growing areas in machine learning, data mining and data science" [27]. Among them, the mean squared error (MSE), root mean squared error (RMSE), R 2 -squared, mean absolute squared error (MAE), standard deviation σ, the mean absolute percentage error (MAPE) [27], Adjusted R 2 , Akaike's information criterion (AIC), Bayesian information criterion (BIC), AICc are "the most common criteria that have been used to measure model performance and select the best model from a set of potential models" [27].
Both models, i.e., the 3RC ECM and Simscape, are already validated by the available information extracted for each of them from the results of MATLAB simulations shown in Figure 13b, Figure 15, Figure 21a,b, that reveal an excellent accuracy due to SOC residual percentage errors being very low. A baseline for comparison is used the estimated value of ADVISOR SOC, as mentioned in the previous sections. The accuracy of both models is better compared to SOC residual error of 2%, usually reported in the literature for similar applications. The first model records a residual error of 1.2% ( Figure 15) and the second one of 1.4% (Figure 21b).
Furthermore, to make a better delimitation between them, additional information is required. The values provided by all six statistical criteria, such as root mean squared error (RMSE), mean squared error (MSE), mean absolute error (MAE), standard deviation (std), mean absolute percentage error (MAPE) and squared, coefficient of determination (R 2 -squared), as are defined in [26,27], is valuable information that makes the difference when two models perform close in terms of accuracy. These values are presented in Table 5, for the 3RC ECM, respectively Table 6, for Simscape Li-ion battery models. All these performance criteria have lower values thus validate, without any doubt, the both models. Moreover, because for both RMSE, MSE, std, MAPE models are very close and R 2 -squared = 0.959 and 0.951, respectively, very close to 1, this is valuable information that indicates how close the values of the data set of the models are and of estimated values ADVISOR SOC. Thus, the overall performance is quite close, with a slight superiority of the 3RC ECR battery model, but the difference is still negligible..

Battery Selection Model
The lower values of all statistical performance criteria from Table 6, quite close to those from Table 5, justifies without doubts the validity of Simscape model. Based on the information collected from both Tables 5 and 6, a rigorous analysis reveals that the both models perform similarly and are suitable for building highly efficient and accurate SOC estimatiors in Part 2. However, an appropriate criterion to compare several candidate models is a hard task for any analyst since some criteria can be disadvantaged by the "model size of estimated parameters while the others could emphasis more on the sample size of a given data" [27]. As a general remark, it can be said that the Simscape Li-ion battery model excels in the following features: simplicity, friendly user-interface, and being fast to implement in real-time.
Also, the battery parameters can be extracted very easily for different chemistry and specifications, and the model is much more realistic in terms of the values of physical model parameters. However, a comparison of SOC performance, more relevant for highlighting the strengths and weaknesses of both models, is provided by the validation results of the MATLAB simulations shown in Figures 13b  and 21a,b for the same FTP-75 driving cycle profile and equal SOC value estimated by the ADVISOR simulator 3.2. Following this test, excellent SOC accuracy is observed for both models. Valuable information on SOC accuracy is extracted from residual SOC errors generated in MATLAB and analyzed in Figures 15 and 21b. The results of the MATLAB simulations show a residual SOC error in the case of the 3RC ECM Li-ion battery model in the range [−1, 1, 0, 4], and the percentage of residual SOC error is below 1.2%. In contrast, for the Simscape model, the percentage of residual SOC error is less than 1.4%, compared to the typical value of 2% reported in the literature for similar applications. The features mentioned above and in addition to the SOC percentage of excellent residual error, much lower than the typical value of 2% reported in the literature, are sometimes even better, strongly recommending the choice of the Simscape model for a wide variety of HEV and EVs applications. In conclusion, these results are encouraging for the next step to develop the most suitable SOC estimators in Part 2. The results of the MATLAB simulations will confirm that all three SOC estimators work better for a design based on choosing the Simulink model, rather than adopting 3RC ECM battery mode, even if the accuracy performance of both battery models is equally sensitive.

Conclusions
In the current research paper, the following most relevant contributions of the authors can be highlighted: • Selection for the same SAFT Li-ion battery of two models, the first model being a 3RC ECM and the second a highly non-linear MATLAB Simscape model, well known for its simplicity, excellent accuracy, practical value and suitable for real-time implementation. • Model development in continuous and discrete-time state space representation.

•
Validation of both models based on same FTP-75 driving cycle current profile test, using ADVISOR 3.2 software tool.

•
Thermal model design and Simulink implementation.

•
Opening of new research topics directions related to energy management systems, optimization techniques and HEV applications.
Based on six statistical criteria values for all three SOC estimators, as a behavior response to an FTP-75 driving cycle profile test, it was possible to decide based on SOC accuracy performance if both models are suitable to be used in Part 2 [30], for adaptive Kalman filter SOC estimators design and implementation. Furthermore, the overall performance analysis indicates that both models are accurate and suitable to be used in Part 2 [30]. In future work, our investigations will continue an improved modelling approach, by integrating the effect of degradation, temperature and SOC effects. New directions of research in energy management systems to develop power optimization techniques and for possible extensions to learning machine SOC estimation techniques will be a great challenge.