Generalized Characterization Methodology for Performance Modelling of Lithium-Ion Batteries

: Lithium-ion (Li-ion) batteries are complex energy storage devices with their performance behavior highly dependent on the operating conditions (i.e., temperature, load current, and state-of-charge ( SOC )). Thus, in order to evaluate their techno-economic viability for a certain application, detailed information about Li-ion battery performance behavior becomes necessary. This paper proposes a comprehensive seven-step methodology for laboratory characterization of Li-ion batteries, in which the battery’s performance parameters (i.e., capacity, open-circuit voltage ( OCV ), and impedance) are determined and their dependence on the operating conditions are obtained. Furthermore, this paper proposes a novel hybrid procedure for parameterizing the batteries’ equivalent electrical circuit (EEC), which is used to emulate the batteries’ dynamic behavior. Based on this novel parameterization procedure, the performance model of the studied Li-ion battery is developed and its accuracy is successfully veriﬁed (maximum error lower than 5% and a mean error below 8.5 mV) for various load proﬁles (including a real application proﬁle), thus validating the proposed seven-step characterization methodology.


Introduction
Recent developments of lithium-ion (Li-ion) batteries based on new and improved chemistries have resulted in batteries with high performance, long lifetime and increased safety [1,2]. Thus, Li-ion batteries have become the key energy storage technology for e-mobility applications [3,4]. Furthermore, energy storage systems based on Li-ion batteries are evaluated in different projects where they are used in grid ancillary service applications [5][6][7] or for renewables' grid integration [6,8,9]. Consequently, Li-ion batteries are expected to become the major player in utility-scale applications as stated in different surveys [10,11]. Nevertheless, Li-ion batteries are characterized by a high investment cost in comparison to other energy storage technologies (e.g., lead acid and nickel metal hydride (NiMH) batteries) and their performance behavior is highly influenced by the operating conditions: temperature, load current and state-of-charge (SOC). Thus, in order to gain from the aforementioned advantages, the Li-ion batteries have to be operated in an efficient and cost-effective manner.
This objective can be achieved by relying on accurate performance models, which are able to estimate accurately the dynamic behavior of Li-ion batteries independently of the operating conditions; thus, precise energy management strategies can be defined and accurate sizing of the Li-ion batteries based systems can be realized. Moreover, by relying on accurate performance models, Li-ion batteries can be tested by simulations considering different scenarios, thereby reducing the laboratory testing efforts that are usually cost demanding and time consuming.
Li-ion battery. Section 3 presents the proposed seven-step methodology for characterizing and parameterizing the performance model of the studied Li-ion battery. Based on the obtained results from the characterization test, the performance model of the Li-ion battery is developed in Section 4. In Section 5, extensive results obtained from the validation of the proposed performance model are presented, while concluding remarks are given in Section 6.

Experimental Set-Up
At present, a wide variety of Li-ion battery chemistries are available on the market that are, depending on their characteristics, suitable for a broad range of applications from electro-mobility to back-up power and renewables' grid integration [32]. Nevertheless, none of these chemistries is superior to the others from all perspectives (e.g., energy density, specific energy, cost per cycle, etc.) [33]. Therefore, the characterization methodology presented in this work and the proposed performance model can be applied independently of the type of the selected Li-ion battery chemistry. For exemplifying the proposed characterization methodology, a commercial 2.5 Ah battery based on the lithium iron phosphate/graphite (further referred to as LFP/C), is used in this work.
During the whole characterization process, the battery was placed inside a temperaturecontrolled chamber, as illustrated in Figure 1, and its temperature was continuously monitored using a PT100 sensor. Furthermore, the temperature of the chamber was adjusted in order to obtain the desired temperature on the surface of the LFP/C battery; thus, the temperature values mentioned in this work are the ones of the battery and not the ambient temperature values of the chamber.

Lithium-Ion Battery Characterization Methodology
The main objective of the proposed characterization methodology was to measure the performance parameters of the studied battery and to determine their dependence on the operating conditions; the operating conditions cover the sum of the conditions given by the load current, temperature, and SOC. Based on the obtained results, which are presented throughout this section, the performance model of the LFP/C battery was developed and parameterized (see Section 4).
In order to determine the performance behavior of the tested LFP/C battery, the seven-step methodology, which is presented in Figure 2, was proposed. This methodology is composed of various tests, which were carried out following the sequence illustrated in Figure 2; the first six tests are presented throughout the next subsections, while the last step (i.e., verification test) is extensively discussed in Section 5.

Lithium-Ion Battery Characterization Methodology
The main objective of the proposed characterization methodology was to measure the performance parameters of the studied battery and to determine their dependence on the operating conditions; the operating conditions cover the sum of the conditions given by the load current, temperature, and SOC. Based on the obtained results, which are presented throughout this section, the performance model of the LFP/C battery was developed and parameterized (see Section 4).
In order to determine the performance behavior of the tested LFP/C battery, the seven-step methodology, which is presented in Figure 2, was proposed. This methodology is composed of various tests, which were carried out following the sequence illustrated in Figure 2; the first six tests are presented throughout the next subsections, while the last step (i.e., verification test) is extensively discussed in Section 5.

Preconditioning Test
The objective of the preconditioning test was to remove any possible passivation to which the battery was subjected to, between the manufacturing time and the initial tests. Moreover, an additional goal of this test was to stabilize the battery capacity since the solid electrolyte interface of the batteries is not shaped after the manufacturing process and its porosity/shape changes meaningfully after the first few cycles. Thus, the preconditioning test was composed of five successive charge-discharge cycles performed with 1C-rate (i.e., 2.5 A) at 25 °C [21]. The results obtained for the case of the discharging capacity of the tested battery are presented in Figure 3.  The measured capacity of the LFP/C battery was stable, with a tendency of slight monotonic increase of approximately 0.6% during the five performed cycles; nevertheless, the battery was considered preconditioned since its capacity did not change more than 3% during two consecutive discharges [34].

Relaxation Test
Usually, the behavior of the Li-ion batteries is altered by the concentration gradients of the ionic charge carriers immediately after the load current is switched off [35]. Consequently, a relaxation period, which will allow the battery to reach thermodynamic stability, has to be applied between the switch off of the load current and the desired measurement. In order to determine the optimal length of this relaxation period, the behavior of the OCV of the LFP/C battery was investigated; the battery was fully charged and then discharged at three different SOCs (i.e., 80%, 50% and 20%), where it was kept at OCV conditions for 24 h, respectively [21]. The voltage values measured, with one-second resolution, during the relaxation period were related to the OCV value measured after 24 h (when the battery was considered quasi-stabilized); thus, the OCV error was computed according to Equation (1).

Preconditioning Test
The objective of the preconditioning test was to remove any possible passivation to which the battery was subjected to, between the manufacturing time and the initial tests. Moreover, an additional goal of this test was to stabilize the battery capacity since the solid electrolyte interface of the batteries is not shaped after the manufacturing process and its porosity/shape changes meaningfully after the first few cycles. Thus, the preconditioning test was composed of five successive charge-discharge cycles performed with 1C-rate (i.e., 2.5 A) at 25 • C [21]. The results obtained for the case of the discharging capacity of the tested battery are presented in Figure 3.

Preconditioning Test
The objective of the preconditioning test was to remove any possible passivation to which the battery was subjected to, between the manufacturing time and the initial tests. Moreover, an additional goal of this test was to stabilize the battery capacity since the solid electrolyte interface of the batteries is not shaped after the manufacturing process and its porosity/shape changes meaningfully after the first few cycles. Thus, the preconditioning test was composed of five successive charge-discharge cycles performed with 1C-rate (i.e., 2.5 A) at 25 °C [21]. The results obtained for the case of the discharging capacity of the tested battery are presented in Figure 3.  The measured capacity of the LFP/C battery was stable, with a tendency of slight monotonic increase of approximately 0.6% during the five performed cycles; nevertheless, the battery was considered preconditioned since its capacity did not change more than 3% during two consecutive discharges [34].

Relaxation Test
Usually, the behavior of the Li-ion batteries is altered by the concentration gradients of the ionic charge carriers immediately after the load current is switched off [35]. Consequently, a relaxation period, which will allow the battery to reach thermodynamic stability, has to be applied between the switch off of the load current and the desired measurement. In order to determine the optimal length of this relaxation period, the behavior of the OCV of the LFP/C battery was investigated; the battery was fully charged and then discharged at three different SOCs (i.e., 80%, 50% and 20%), where it was kept at OCV conditions for 24 h, respectively [21]. The voltage values measured, with one-second resolution, during the relaxation period were related to the OCV value measured after 24 h (when the battery was considered quasi-stabilized); thus, the OCV error was computed according to Equation (1). The measured capacity of the LFP/C battery was stable, with a tendency of slight monotonic increase of approximately 0.6% during the five performed cycles; nevertheless, the battery was considered preconditioned since its capacity did not change more than 3% during two consecutive discharges [34].

Relaxation Test
Usually, the behavior of the Li-ion batteries is altered by the concentration gradients of the ionic charge carriers immediately after the load current is switched off [35]. Consequently, a relaxation period, which will allow the battery to reach thermodynamic stability, has to be applied between the switch off of the load current and the desired measurement. In order to determine the optimal length of this relaxation period, the behavior of the OCV of the LFP/C battery was investigated; the battery was fully charged and then discharged at three different SOCs (i.e., 80%, 50% and 20%), where it was kept at OCV conditions for 24 h, respectively [21]. The voltage values measured, with one-second resolution, during the relaxation period were related to the OCV value measured after 24 h (when the battery was considered quasi-stabilized); thus, the OCV error was computed according to Equation (1).
where ε OCV represents the OCV error and V i represents the voltage measured at every second during the 24 h relaxation. Based on the obtained results, which are summarized in Table 1, it was concluded that the OCV value is only slightly influenced by a relaxation time, which varies between 15 min and 24 h. Thus, a relaxation time of one hour was considered enough for the LFP/C battery to reach quasi-thermodynamic stability and to allow for accurate measurement of the desired performance parameters.

Capacity Test
The capacity of Li-ion batteries is dependent on the operating temperature and on the applied load current [15,17]. Moreover, as presented by various researchers, the capacity of Li-ion batteries degrades over time [36][37][38]; however, this analysis is out of the scope of the present paper. The objective of the capacity test was to determine the dependence of the charged and discharged capacity of the LFP/C battery on the two aforementioned parameters. Thus, the capacity of the battery was measured for six different C-rates (i.e., C/4, C/2, 1C, 2C, 3C and 4C) and at four different temperatures (i.e., 15 • C, 25 • C, 35 • C and 45 • C), following the procedure which is presented in [21].
The dependence of the LFP/C battery's discharging capacity on the C-rate and on the temperature is illustrated in Figure 4. The discharged capacity of the LFP/C battery measured at 25 • C decreases as the applied current increases; however, the decrease of the capacity for the considered C-rate interval (C/4-4C) is lower than 1%. This behavior suggests a Peukert number very close to 1, which is in good agreement with the results reported for a similar battery chemistry by Omar et al. [39]. Furthermore, as presented in Figure 4, for the considered temperature interval, a slight increase of the LFP/C battery capacity with increasing temperature was measured, which might have been caused by the increased electronic and ionic conductivity in the electrode and electrolyte at high temperatures; similar results were reported in [17,39,40].
where εOCV represents the OCV error and Vi represents the voltage measured at every second during the 24 h relaxation. Based on the obtained results, which are summarized in Table 1, it was concluded that the OCV value is only slightly influenced by a relaxation time, which varies between 15 min and 24 h. Thus, a relaxation time of one hour was considered enough for the LFP/C battery to reach quasithermodynamic stability and to allow for accurate measurement of the desired performance parameters.

Capacity Test
The capacity of Li-ion batteries is dependent on the operating temperature and on the applied load current [15,17]. Moreover, as presented by various researchers, the capacity of Li-ion batteries degrades over time [36][37][38]; however, this analysis is out of the scope of the present paper. The objective of the capacity test was to determine the dependence of the charged and discharged capacity of the LFP/C battery on the two aforementioned parameters. Thus, the capacity of the battery was measured for six different C-rates (i.e., C/4, C/2, 1C, 2C, 3C and 4C) and at four different temperatures (i.e., 15 °C, 25 °C, 35 °C and 45 °C), following the procedure which is presented in [21].
The dependence of the LFP/C battery's discharging capacity on the C-rate and on the temperature is illustrated in Figure 4. The discharged capacity of the LFP/C battery measured at 25 °C decreases as the applied current increases; however, the decrease of the capacity for the considered C-rate interval (C/4-4C) is lower than 1%. This behavior suggests a Peukert number very close to 1, which is in good agreement with the results reported for a similar battery chemistry by Omar et al. [39]. Furthermore, as presented in Figure 4, for the considered temperature interval, a slight increase of the LFP/C battery capacity with increasing temperature was measured, which might have been caused by the increased electronic and ionic conductivity in the electrode and electrolyte at high temperatures; similar results were reported in [17,39,40].

Open-Circuit Voltage versus State-of-Charge Test
Similar to most of the Li-ion battery parameters, the OCV is strongly dependent on the operating conditions (i.e., temperature, SOC, etc.). Therefore, the goal of this test was to measure the OCV-SOC characteristic of the LFP/C battery at different temperatures for both charging and discharging conditions; the OCV was measured in steps of 5% SOC following the methodology described in [21].
For instance, the OCV-SOC characteristic of the tested battery, measured at 25 • C, is presented in Figure 5. For the measured characteristic, a hysteresis effect is observed, which is more pronounced for the 10%-40% SOC interval. This behavior of the OCV is intrinsic for Li-ion batteries based on active materials, which perform a two-phase transition lithium insertion and extraction process, as is the case of the LFP material [41]. The hysteresis effect is mainly influenced by the previous history (e.g., charging/discharging C-rate and relaxation time) of the LFP/C battery and can be modelled following various approaches reported in literatures [42,43].

Open-Circuit Voltage versus State-of-Charge Test
Similar to most of the Li-ion battery parameters, the OCV is strongly dependent on the operating conditions (i.e., temperature, SOC, etc.). Therefore, the goal of this test was to measure the OCV-SOC characteristic of the LFP/C battery at different temperatures for both charging and discharging conditions; the OCV was measured in steps of 5% SOC following the methodology described in [21].
For instance, the OCV-SOC characteristic of the tested battery, measured at 25 °C, is presented in Figure 5. For the measured characteristic, a hysteresis effect is observed, which is more pronounced for the 10%-40% SOC interval. This behavior of the OCV is intrinsic for Li-ion batteries based on active materials, which perform a two-phase transition lithium insertion and extraction process, as is the case of the LFP material [41]. The hysteresis effect is mainly influenced by the previous history (e.g., charging/discharging C-rate and relaxation time) of the LFP/C battery and can be modelled following various approaches reported in literatures [42,43].

Hybrid Pulse Power Characterization Test
The impedance represents a very important performance parameter of the Li-ion batteries since it defines the power capability of the battery. Moreover, it is the parameter which describes the dynamic behavior of Li-ion batteries. As most of the Li-ion battery parameters, the impedance shows a non-linear behavior depending on the operating conditions such as SOC, temperature, C-rate [44,45]. Furthermore, the battery impedance is highly dependent on the battery's state-of-health (SOH) [44]; however, this is out of the scope of the present work. In order to account for all the aforementioned dependences, the impedance of the LFP/C battery was measured following a modified version of the hybrid pulse power characterization (HPPC) test. The measurement profile applied to the LFP/C battery is presented in Figure 6 and consists of consecutive charging and discharging DC current pulses of different C-rates (i.e., 0.1C, 0.5C, 1C, 2C, 3C and 4C); the pulse duration was set to 18 s [34] and the relaxation time between each two consecutive pulses was set to 15 min.

Hybrid Pulse Power Characterization Test
The impedance represents a very important performance parameter of the Li-ion batteries since it defines the power capability of the battery. Moreover, it is the parameter which describes the dynamic behavior of Li-ion batteries. As most of the Li-ion battery parameters, the impedance shows a non-linear behavior depending on the operating conditions such as SOC, temperature, C-rate [44,45]. Furthermore, the battery impedance is highly dependent on the battery's state-of-health (SOH) [44]; however, this is out of the scope of the present work. In order to account for all the aforementioned dependences, the impedance of the LFP/C battery was measured following a modified version of the hybrid pulse power characterization (HPPC) test. The measurement profile applied to the LFP/C battery is presented in Figure 6 and consists of consecutive charging and discharging DC current pulses of different C-rates (i.e., 0.1C, 0.5C, 1C, 2C, 3C and 4C); the pulse duration was set to 18 s [34] and the relaxation time between each two consecutive pulses was set to 15 min.

Open-Circuit Voltage versus State-of-Charge Test
Similar to most of the Li-ion battery parameters, the OCV is strongly dependent on the operating conditions (i.e., temperature, SOC, etc.). Therefore, the goal of this test was to measure the OCV-SOC characteristic of the LFP/C battery at different temperatures for both charging and discharging conditions; the OCV was measured in steps of 5% SOC following the methodology described in [21].
For instance, the OCV-SOC characteristic of the tested battery, measured at 25 °C, is presented in Figure 5. For the measured characteristic, a hysteresis effect is observed, which is more pronounced for the 10%-40% SOC interval. This behavior of the OCV is intrinsic for Li-ion batteries based on active materials, which perform a two-phase transition lithium insertion and extraction process, as is the case of the LFP material [41]. The hysteresis effect is mainly influenced by the previous history (e.g., charging/discharging C-rate and relaxation time) of the LFP/C battery and can be modelled following various approaches reported in literatures [42,43].

Hybrid Pulse Power Characterization Test
The impedance represents a very important performance parameter of the Li-ion batteries since it defines the power capability of the battery. Moreover, it is the parameter which describes the dynamic behavior of Li-ion batteries. As most of the Li-ion battery parameters, the impedance shows a non-linear behavior depending on the operating conditions such as SOC, temperature, C-rate [44,45]. Furthermore, the battery impedance is highly dependent on the battery's state-of-health (SOH) [44]; however, this is out of the scope of the present work. In order to account for all the aforementioned dependences, the impedance of the LFP/C battery was measured following a modified version of the hybrid pulse power characterization (HPPC) test. The measurement profile applied to the LFP/C battery is presented in Figure 6 and consists of consecutive charging and discharging DC current pulses of different C-rates (i.e., 0.1C, 0.5C, 1C, 2C, 3C and 4C); the pulse duration was set to 18 s [34] and the relaxation time between each two consecutive pulses was set to 15 min.  Furthermore, in order to consider the dependence of the impedance on SOC and temperature, the test profile presented in Figure 6 was applied to the 5%-95% SOC interval (considering a 5% SOC resolution) and for the same temperatures (i.e., 15 • C, 25 • C, 35 • C and 45 • C) considered during measurement of the other performance parameters.
The typical voltage response of a Li-ion battery to a discharging current pulse is illustrated in Figure 7; a similar voltage behavior is obtained for a charging current pulse. As presented in the ISO 12405-1:2011 Standard [34], the charging and discharging impedance of the battery could be calculated for various pulse lengths (i.e., 0.1 s, 2 s and 18 s). Furthermore, the impedances corresponding to these pulse lengths are caused by/related to different electrochemical processes inside the battery [44]. Nevertheless, in this work, impedance values determined for a current length of 18 s are presented. Based on the current and voltage measurements, the impedance of the Li-ion battery is computed based on Ohm's law, as: where R i represents the impedance of the battery, ∆V represents the change in the voltage due to the applied current pulse, ∆I represents the change in the current, V 1 represents the voltage measured after 18 s, V 0 represents the voltage measured just before applying the current pulse, and I 1 represents the amplitude of the applied current pulse. Furthermore, in order to consider the dependence of the impedance on SOC and temperature, the test profile presented in Figure 6 was applied to the 5%-95% SOC interval (considering a 5% SOC resolution) and for the same temperatures (i.e., 15 °C, 25 °C, 35 °C and 45 °C) considered during measurement of the other performance parameters.
The typical voltage response of a Li-ion battery to a discharging current pulse is illustrated in Figure 7; a similar voltage behavior is obtained for a charging current pulse. As presented in the ISO 12405-1:2011 Standard [34], the charging and discharging impedance of the battery could be calculated for various pulse lengths (i.e., 0.1 s, 2 s and 18 s). Furthermore, the impedances corresponding to these pulse lengths are caused by/related to different electrochemical processes inside the battery [44]. Nevertheless, in this work, impedance values determined for a current length of 18 s are presented. Based on the current and voltage measurements, the impedance of the Li-ion battery is computed based on Ohm's law, as: where Ri represents the impedance of the battery, ∆V represents the change in the voltage due to the applied current pulse, ∆I represents the change in the current, V1 represents the voltage measured after 18 s, V0 represents the voltage measured just before applying the current pulse, and I1 represents the amplitude of the applied current pulse. The voltage response of the LFP battery to the current profile measured at 50% SOC is illustrated in Figure 8.   The voltage response of the LFP battery to the current profile measured at 50% SOC is illustrated in Figure 8. Furthermore, in order to consider the dependence of the impedance on SOC and temperature, the test profile presented in Figure 6 was applied to the 5%-95% SOC interval (considering a 5% SOC resolution) and for the same temperatures (i.e., 15 °C, 25 °C, 35 °C and 45 °C) considered during measurement of the other performance parameters.
The typical voltage response of a Li-ion battery to a discharging current pulse is illustrated in Figure 7; a similar voltage behavior is obtained for a charging current pulse. As presented in the ISO 12405-1:2011 Standard [34], the charging and discharging impedance of the battery could be calculated for various pulse lengths (i.e., 0.1 s, 2 s and 18 s). Furthermore, the impedances corresponding to these pulse lengths are caused by/related to different electrochemical processes inside the battery [44]. Nevertheless, in this work, impedance values determined for a current length of 18 s are presented. Based on the current and voltage measurements, the impedance of the Li-ion battery is computed based on Ohm's law, as: where Ri represents the impedance of the battery, ∆V represents the change in the voltage due to the applied current pulse, ∆I represents the change in the current, V1 represents the voltage measured after 18 s, V0 represents the voltage measured just before applying the current pulse, and I1 represents the amplitude of the applied current pulse. The voltage response of the LFP battery to the current profile measured at 50% SOC is illustrated in Figure 8.  0.1 C-rate 0.5 C-rate 1 C-rate 2 C-rate 3 C-rate 4 C-rate Based on battery's voltage response (similar to those presented in Figure 8), and considering Equation (2), the charging and discharging impedance of the tested LFP/C battery was determined for all the considered conditions; for example, Figure 9 presents the measured impedance of the LFP/C during pulse charging and discharging with 1C-rate (i.e., 2.5 A) at T = 25 • C. Both impedance characteristics are showing a parabolic dependence on the SOC, following a nearly flat region between 20% and 90% SOC. Based on battery's voltage response (similar to those presented in Figure 8), and considering Equation (2), the charging and discharging impedance of the tested LFP/C battery was determined for all the considered conditions; for example, Figure 9 presents the measured impedance of the LFP/C during pulse charging and discharging with 1C-rate (i.e., 2.5 A) at T = 25 °C. Both impedance characteristics are showing a parabolic dependence on the SOC, following a nearly flat region between 20% and 90% SOC.

Electrochemical Impedance Spectroscopy Test
EIS has developed as a reliable method for characterization and modelling the performance behavior of the Li-ion batteries. The small-signal impedance of the Li-ion batteries is determined by applying a small sinusoidal current (galavantostatic mode) or voltage (potentiostatic mode) and measuring the amplitude and phase shift of the output voltage or current, respectively. This procedure is repeated for a sweep of frequencies, and thus the battery impedance spectrum is obtained.
In this work, the EIS measurements were performed in galvanostatic mode and for the frequency range 10 kHz-10 mHz. Furthermore, all the EIS measurements were performed without superimposed DC current; thus, the influence of the current on the small-signal AC impedance was not determined. As for the other considered performance parameters, the small-signal impedance was measured at different temperatures and for the entire SOC interval (0%-100% SOC with 5% SOC resolution). The dependence of the LFP/C battery's impedance spectra on the SOC and temperature is presented in Figures 10 and 11, respectively.

Electrochemical Impedance Spectroscopy Test
EIS has developed as a reliable method for characterization and modelling the performance behavior of the Li-ion batteries. The small-signal impedance of the Li-ion batteries is determined by applying a small sinusoidal current (galavantostatic mode) or voltage (potentiostatic mode) and measuring the amplitude and phase shift of the output voltage or current, respectively. This procedure is repeated for a sweep of frequencies, and thus the battery impedance spectrum is obtained.
In this work, the EIS measurements were performed in galvanostatic mode and for the frequency range 10 kHz-10 mHz. Furthermore, all the EIS measurements were performed without superimposed DC current; thus, the influence of the current on the small-signal AC impedance was not determined. As for the other considered performance parameters, the small-signal impedance was measured at different temperatures and for the entire SOC interval (0%-100% SOC with 5% SOC resolution). The dependence of the LFP/C battery's impedance spectra on the SOC and temperature is presented in Figures 10 and 11, respectively. Based on battery's voltage response (similar to those presented in Figure 8), and considering Equation (2), the charging and discharging impedance of the tested LFP/C battery was determined for all the considered conditions; for example, Figure 9 presents the measured impedance of the LFP/C during pulse charging and discharging with 1C-rate (i.e., 2.5 A) at T = 25 °C. Both impedance characteristics are showing a parabolic dependence on the SOC, following a nearly flat region between 20% and 90% SOC.

Electrochemical Impedance Spectroscopy Test
EIS has developed as a reliable method for characterization and modelling the performance behavior of the Li-ion batteries. The small-signal impedance of the Li-ion batteries is determined by applying a small sinusoidal current (galavantostatic mode) or voltage (potentiostatic mode) and measuring the amplitude and phase shift of the output voltage or current, respectively. This procedure is repeated for a sweep of frequencies, and thus the battery impedance spectrum is obtained.
In this work, the EIS measurements were performed in galvanostatic mode and for the frequency range 10 kHz-10 mHz. Furthermore, all the EIS measurements were performed without superimposed DC current; thus, the influence of the current on the small-signal AC impedance was not determined. As for the other considered performance parameters, the small-signal impedance was measured at different temperatures and for the entire SOC interval (0%-100% SOC with 5% SOC resolution). The dependence of the LFP/C battery's impedance spectra on the SOC and temperature is presented in Figures 10 and 11, respectively.

Lithium-Ion Battery Performance Model
Li-ion battery performance models can be classified into three main categories: electrochemical, electrical, and mathematical models, depending on the development and parameterization approach. These modelling approaches are characterized by different degrees of complexity, different accuracy levels and could be used for different purposes (e.g., battery design, system simulation). Nevertheless, this work has focused on the development and parameterization of the LFP/C battery performance model using an electrical modelling approach.

Equivalent Electrical Circuit Model
Li-ion batteries' electrical models generally use EECs, composed of simple (e.g., resistors, capacitors, inductors) or more complex elements (e.g., constant phase elements, Warburg elements, etc.) and a voltage source to express the dynamic behavior of batteries. The basic configuration of an EEC-based performance model is presented in Figure 12. The structure of the EEC can have different configurations depending on various aspects, such as: model accuracy, model computation time, Liion battery chemistry, etc. According to the configuration shown in Figure 12, the voltage of the battery is computed according to Equation (3).
where Vbat represents the voltage of the battery, VOC represents the OCV, and VEEC represents the voltage drop across the EEC, which is used to model the dynamics of the battery. The voltage drop across EEC's impedance is obtained as the sum of various over-voltages, which are caused by different processes occurring inside the battery with different time constants during charging or discharging [46]: where  Figure 11. Dependence of the impedance spectra on the temperature; SOC = 50%.

Lithium-Ion Battery Performance Model
Li-ion battery performance models can be classified into three main categories: electrochemical, electrical, and mathematical models, depending on the development and parameterization approach. These modelling approaches are characterized by different degrees of complexity, different accuracy levels and could be used for different purposes (e.g., battery design, system simulation). Nevertheless, this work has focused on the development and parameterization of the LFP/C battery performance model using an electrical modelling approach.

Equivalent Electrical Circuit Model
Li-ion batteries' electrical models generally use EECs, composed of simple (e.g., resistors, capacitors, inductors) or more complex elements (e.g., constant phase elements, Warburg elements, etc.) and a voltage source to express the dynamic behavior of batteries. The basic configuration of an EEC-based performance model is presented in Figure 12

Lithium-Ion Battery Performance Model
Li-ion battery performance models can be classified into three main categories: electrochemical, electrical, and mathematical models, depending on the development and parameterization approach. These modelling approaches are characterized by different degrees of complexity, different accuracy levels and could be used for different purposes (e.g., battery design, system simulation). Nevertheless, this work has focused on the development and parameterization of the LFP/C battery performance model using an electrical modelling approach.

Equivalent Electrical Circuit Model
Li-ion batteries' electrical models generally use EECs, composed of simple (e.g., resistors, capacitors, inductors) or more complex elements (e.g., constant phase elements, Warburg elements, etc.) and a voltage source to express the dynamic behavior of batteries. The basic configuration of an EEC-based performance model is presented in Figure 12. The structure of the EEC can have different configurations depending on various aspects, such as: model accuracy, model computation time, Liion battery chemistry, etc. According to the configuration shown in Figure 12, the voltage of the battery is computed according to Equation (3).
where Vbat represents the voltage of the battery, VOC represents the OCV, and VEEC represents the voltage drop across the EEC, which is used to model the dynamics of the battery. The voltage drop across EEC's impedance is obtained as the sum of various over-voltages, which are caused by different processes occurring inside the battery with different time constants during charging or discharging [46]: where Vohmic represents the ohmic over-voltage (caused by the resistance of the poles, current collectors, electrolytes, etc.), Vchtr represents the charge transfer over-voltage (caused by According to the configuration shown in Figure 12, the voltage of the battery is computed according to Equation (3).
where V bat represents the voltage of the battery, V OC represents the OCV, and V EEC represents the voltage drop across the EEC, which is used to model the dynamics of the battery. The voltage drop across EEC's impedance is obtained as the sum of various over-voltages, which are caused by different processes occurring inside the battery with different time constants during charging or discharging [46]: where V ohmic represents the ohmic over-voltage (caused by the resistance of the poles, current collectors, electrolytes, etc.), V chtr represents the charge transfer over-voltage (caused by electrochemical reaction at inner surfaces), and V diff represents the diffusion over-voltage (caused by a deficit or surplus of reactants at the reaction's location). Thus, to model the voltage behavior of the LFP/C battery, the contributions of the V OC and of the V EEC components have to be known. The contribution of the OCV to the battery voltage was obtained for various conditions (i.e., SOCs and temperatures) during the battery characterization, as presented in Section 3.4. To obtain the contribution of the V EEC to the LFP/C battery voltage, the EEC had to be parameterized.
The parameterization of the EEC was realized by curve fitting the impedance spectra of the LFP/C battery, which were measured at different conditions as presented in Section 3.6. The configuration of the EEC which was used for curve fitting, and consequently to express the dynamics of the LFP/C battery, is illustrated in Figure 13. The curve fitting procedure was based on a complex non-linear least-square (CNLS) algorithm, where both real and imaginary components of the measured impedance were fitted simultaneously using a least square minimization; an example of the impedance spectrum curve fitting results is shown in Figure 14. electrochemical reaction at inner surfaces), and Vdiff represents the diffusion over-voltage (caused by a deficit or surplus of reactants at the reaction's location). Thus, to model the voltage behavior of the LFP/C battery, the contributions of the VOC and of the VEEC components have to be known. The contribution of the OCV to the battery voltage was obtained for various conditions (i.e., SOCs and temperatures) during the battery characterization, as presented in Section 3.4. To obtain the contribution of the VEEC to the LFP/C battery voltage, the EEC had to be parameterized.
The parameterization of the EEC was realized by curve fitting the impedance spectra of the LFP/C battery, which were measured at different conditions as presented in Section 3.6. The configuration of the EEC which was used for curve fitting, and consequently to express the dynamics of the LFP/C battery, is illustrated in Figure 13. The curve fitting procedure was based on a complex non-linear least-square (CNLS) algorithm, where both real and imaginary components of the measured impedance were fitted simultaneously using a least square minimization; an example of the impedance spectrum curve fitting results is shown in Figure 14.  As presented in Figure 14, the EEC configuration based on a series inductance Ls, a series resistance Rs, and two ZARC elements (i.e., parallel connections of a resistor and a constant phase element) was able to accurately fit the measured impedance spectrum of the tested LFP/C battery. Therefore, this EEC configuration was used to fit all the measured impedance spectra (at the conditions discussed in Section 3.6) and thus to model the performance behavior of the LFP/C battery. The impedance of the selected EEC is given by: where R1, Q1, and N1 represent the resistance, the generalized capacitance and the depression factor of the first ZARC element, respectively; the same explanations are valid for the second ZARC of the EEC.
A main advantage of the EIS technique over other parameterization techniques is represented by its feature to relate EEC elements to the physicochemical process that occur inside the battery [28,30]. Thus, in the present work, the first term in Equation (5) is used to describe the inductive behavior at high frequencies, the second term is used to calculate the over-voltage caused by the ohmic resistance (Vohmic), the third term is used to calculate the over-voltage corresponding to the electrochemical reaction at inner surfaces), and Vdiff represents the diffusion over-voltage (caused by a deficit or surplus of reactants at the reaction's location). Thus, to model the voltage behavior of the LFP/C battery, the contributions of the VOC and of the VEEC components have to be known. The contribution of the OCV to the battery voltage was obtained for various conditions (i.e., SOCs and temperatures) during the battery characterization, as presented in Section 3.4. To obtain the contribution of the VEEC to the LFP/C battery voltage, the EEC had to be parameterized.
The parameterization of the EEC was realized by curve fitting the impedance spectra of the LFP/C battery, which were measured at different conditions as presented in Section 3.6. The configuration of the EEC which was used for curve fitting, and consequently to express the dynamics of the LFP/C battery, is illustrated in Figure 13. The curve fitting procedure was based on a complex non-linear least-square (CNLS) algorithm, where both real and imaginary components of the measured impedance were fitted simultaneously using a least square minimization; an example of the impedance spectrum curve fitting results is shown in Figure 14.  As presented in Figure 14, the EEC configuration based on a series inductance Ls, a series resistance Rs, and two ZARC elements (i.e., parallel connections of a resistor and a constant phase element) was able to accurately fit the measured impedance spectrum of the tested LFP/C battery. Therefore, this EEC configuration was used to fit all the measured impedance spectra (at the conditions discussed in Section 3.6) and thus to model the performance behavior of the LFP/C battery. The impedance of the selected EEC is given by: where R1, Q1, and N1 represent the resistance, the generalized capacitance and the depression factor of the first ZARC element, respectively; the same explanations are valid for the second ZARC of the EEC.
A main advantage of the EIS technique over other parameterization techniques is represented by its feature to relate EEC elements to the physicochemical process that occur inside the battery [28,30]. Thus, in the present work, the first term in Equation (5) is used to describe the inductive behavior at high frequencies, the second term is used to calculate the over-voltage caused by the ohmic resistance (Vohmic), the third term is used to calculate the over-voltage corresponding to the Measured spectrum Fitted spectrum Figure 14. Measured and fitted impedance spectrum using the EEC presented in Figure 13.
As presented in Figure 14, the EEC configuration based on a series inductance L s , a series resistance R s , and two ZARC elements (i.e., parallel connections of a resistor and a constant phase element) was able to accurately fit the measured impedance spectrum of the tested LFP/C battery. Therefore, this EEC configuration was used to fit all the measured impedance spectra (at the conditions discussed in Section 3.6) and thus to model the performance behavior of the LFP/C battery. The impedance of the selected EEC is given by: where R 1 , Q 1 , and N 1 represent the resistance, the generalized capacitance and the depression factor of the first ZARC element, respectively; the same explanations are valid for the second ZARC of the EEC. A main advantage of the EIS technique over other parameterization techniques is represented by its feature to relate EEC elements to the physicochemical process that occur inside the battery [28,30]. Thus, in the present work, the first term in Equation (5) is used to describe the inductive behavior at high frequencies, the second term is used to calculate the over-voltage caused by the ohmic resistance (V ohmic ), the third term is used to calculate the over-voltage corresponding to the charge-transfer process (V chtr ), while the last term is used to determine the over-voltage corresponding to the diffusion process (V diff ). Consequently, based on the aforementioned aspects, the battery voltage Equation (3) was rewritten as given in Equation (6). Furthermore, considering Equation (6), the configuration of the EEC, which was used to model the performance behavior of the LFP/C battery, is presented in Figure 15.
Batteries 2016, 2, 37 11 of 19 charge-transfer process (Vchtr), while the last term is used to determine the over-voltage corresponding to the diffusion process (Vdiff). Consequently, based on the aforementioned aspects, the battery voltage Equation (3) was rewritten as given in Equation (6). Furthermore, considering Equation (6), the configuration of the EEC, which was used to model the performance behavior of the LFP/C battery, is presented in Figure 15.

Dependence on Load Current
Besides being dependent on SOC and temperature-dependences obtained from the EIS measurements (see Section 3.6)-the parameters of the EEC are also dependent on the load current as illustrated in Figure 15. Nevertheless, because all the EIS measurements were performed without superimposed DC current (see Section 3.6), the effect of the load current on the parameters of the EEC was not determined. Among the parameters of the EEC, the charge transfer resistance R1 shows a highly non-linear dependence on the load current. This non-linear dependence of the charge-transfer resistance, and of the corresponding over-voltage, on the LFP/C battery current is presented in Figure  16.
Theoretically, the exponential dependence between the Li-ion battery current and the potential can be determined from the Butler-Volmer equation given in Equation (7) [47].
where i is the electrode current density, i0 is the exchange current density, α is the symmetry factor, n is the number of electrons, R is the universal gas constant (8.3144 J·K −1 ·mol −1 ), F is the Faraday constant (96,485.339 C·mol −1 ), η is the overpotential (i.e., the difference between the electrode potential and equilibrium potential).

Dependence on Load Current
Besides being dependent on SOC and temperature-dependences obtained from the EIS measurements (see Section 3.6)-the parameters of the EEC are also dependent on the load current as illustrated in Figure 15. Nevertheless, because all the EIS measurements were performed without superimposed DC current (see Section 3.6), the effect of the load current on the parameters of the EEC was not determined. Among the parameters of the EEC, the charge transfer resistance R 1 shows a highly non-linear dependence on the load current. This non-linear dependence of the charge-transfer resistance, and of the corresponding over-voltage, on the LFP/C battery current is presented in Figure 16.
charge-transfer process (Vchtr), while the last term is used to determine the over-voltage corresponding to the diffusion process (Vdiff). Consequently, based on the aforementioned aspects, the battery voltage Equation (3) was rewritten as given in Equation (6). Furthermore, considering Equation (6), the configuration of the EEC, which was used to model the performance behavior of the LFP/C battery, is presented in Figure 15.

Dependence on Load Current
Besides being dependent on SOC and temperature-dependences obtained from the EIS measurements (see Section 3.6)-the parameters of the EEC are also dependent on the load current as illustrated in Figure 15. Nevertheless, because all the EIS measurements were performed without superimposed DC current (see Section 3.6), the effect of the load current on the parameters of the EEC was not determined. Among the parameters of the EEC, the charge transfer resistance R1 shows a highly non-linear dependence on the load current. This non-linear dependence of the charge-transfer resistance, and of the corresponding over-voltage, on the LFP/C battery current is presented in Figure  16.
Theoretically, the exponential dependence between the Li-ion battery current and the potential can be determined from the Butler-Volmer equation given in Equation (7) [47].
where i is the electrode current density, i0 is the exchange current density, α is the symmetry factor, n is the number of electrons, R is the universal gas constant (8.3144 J·K −1 ·mol −1 ), F is the Faraday constant (96,485.339 C·mol −1 ), η is the overpotential (i.e., the difference between the electrode potential and equilibrium potential).  Theoretically, the exponential dependence between the Li-ion battery current and the potential can be determined from the Butler-Volmer equation given in Equation (7) [47].
where i is the electrode current density, i 0 is the exchange current density, α is the symmetry factor, n is the number of electrons, R is the universal gas constant (8.3144 J·K −1 ·mol −1 ), F is the Faraday constant (96,485.339 C·mol −1 ), η is the overpotential (i.e., the difference between the electrode potential and equilibrium potential).
However, information about various parameters given in Equation (7) were not available and require specialized laboratory equipment in order to determine them. Consequently, in order to determine the dependence of the LFP/C battery voltage on the load current, the following approach was followed: a comparison was performed between the battery voltage response measured over a period of 18 s in the laboratory (during the HPPC tests-see Section 3.5) and the simulated voltage response. The obtained difference was considered to describe the impact of the current on the voltage of the LFP/C battery. The aforementioned procedure was repeated for all SOCs, temperatures, and C-rates, which were considered during the battery impedance measurement by HPPC test. The obtained voltage correction factors, assigned to the effect of the current on the LFP/C battery voltage were stored as a 3D look-up table; for instance, Figure 17 presents the look-up table implementation of the voltage correction factors for the 25 • C dimension. However, information about various parameters given in Equation (7) were not available and require specialized laboratory equipment in order to determine them. Consequently, in order to determine the dependence of the LFP/C battery voltage on the load current, the following approach was followed: a comparison was performed between the battery voltage response measured over a period of 18 s in the laboratory (during the HPPC tests-see Section 3.5) and the simulated voltage response. The obtained difference was considered to describe the impact of the current on the voltage of the LFP/C battery. The aforementioned procedure was repeated for all SOCs, temperatures, and C-rates, which were considered during the battery impedance measurement by HPPC test. The obtained voltage correction factors, assigned to the effect of the current on the LFP/C battery voltage were stored as a 3D look-up table; for instance, Figure 17 presents the look-up table implementation of the voltage correction factors for the 25 °C dimension.

Model Implementation
Based on the voltage Equation (6) and considering the dependences of the performance parameters on the operating conditions, the performance model of the LFP/C battery was implemented according to the block diagram proposed in Figure 18. The inputs of the proposed performance model for the LFP/C battery are the applied load current, the initial SOC and the temperature, while the output of the model is the estimated battery voltage; additionally, the model can return information about the battery SOC at every moment. Furthermore, since no thermal model was developed for the LFP/C battery, the temperature of the battery is considered constant during the whole simulation.
The parameterization of the proposed performance model was performed based on the results obtained from the characterization test, which was presented in detail in Section 3. The parameters of the LFP/C battery were implemented as 2D look-up tables (i.e., capacity, OCV, EEC's parameters)

Model Implementation
Based on the voltage Equation (6) and considering the dependences of the performance parameters on the operating conditions, the performance model of the LFP/C battery was implemented according to the block diagram proposed in Figure 18. However, information about various parameters given in Equation (7) were not available and require specialized laboratory equipment in order to determine them. Consequently, in order to determine the dependence of the LFP/C battery voltage on the load current, the following approach was followed: a comparison was performed between the battery voltage response measured over a period of 18 s in the laboratory (during the HPPC tests-see Section 3.5) and the simulated voltage response. The obtained difference was considered to describe the impact of the current on the voltage of the LFP/C battery. The aforementioned procedure was repeated for all SOCs, temperatures, and C-rates, which were considered during the battery impedance measurement by HPPC test. The obtained voltage correction factors, assigned to the effect of the current on the LFP/C battery voltage were stored as a 3D look-up table; for instance, Figure 17 presents the look-up table implementation of the voltage correction factors for the 25 °C dimension.

Model Implementation
Based on the voltage Equation (6) and considering the dependences of the performance parameters on the operating conditions, the performance model of the LFP/C battery was implemented according to the block diagram proposed in Figure 18. The inputs of the proposed performance model for the LFP/C battery are the applied load current, the initial SOC and the temperature, while the output of the model is the estimated battery voltage; additionally, the model can return information about the battery SOC at every moment. Furthermore, since no thermal model was developed for the LFP/C battery, the temperature of the battery is considered constant during the whole simulation.
The parameterization of the proposed performance model was performed based on the results obtained from the characterization test, which was presented in detail in Section 3. The parameters of the LFP/C battery were implemented as 2D look-up tables (i.e., capacity, OCV, EEC's parameters) The inputs of the proposed performance model for the LFP/C battery are the applied load current, the initial SOC and the temperature, while the output of the model is the estimated battery voltage; additionally, the model can return information about the battery SOC at every moment. Furthermore, since no thermal model was developed for the LFP/C battery, the temperature of the battery is considered constant during the whole simulation.
The parameterization of the proposed performance model was performed based on the results obtained from the characterization test, which was presented in detail in Section 3. The parameters of the LFP/C battery were implemented as 2D look-up tables (i.e., capacity, OCV, EEC's parameters) and 3D-look-up table (i.e., Butler-Volmer correction voltage) in order to consider their dependences on the operating conditions (i.e., SOC, temperature, and load current). Because for similar operating conditions, the capacity of the Li-ion batteries is different between charging and discharging case (i.e., Coulombic efficiency different than 100%), two separate look-up tables were implemented for the charging and discharging capacity, respectively. Because of the hysteresis effect (see Figure 5), which is inherent to LFP/C battery chemistry, a similar approach was followed for implementing the OCV characteristic; however, for other Li-ion battery chemistries, a single look-up table is enough to model the OCV characteristic for both charging and discharging cases.
During the fitting process of the measured impedance spectra, it was found out that only R s , R 1 , R 2 , Q 1 , and Q 2 are dependent on the SOC and temperature at which the EIS measurements had been performed. Consequently, these parameters were implemented into the performance model as 2D look-up tables since their values are dependent on the SOC and temperature; Figure 19 exemplifies the look-up table implementation of the series resistance R s . The other parameters of the EEC, L s , N 1 , and N 2 were found to be constant, independent of the operating conditions: L s = 1.8 × 10 −8 H, N 1 = 0.5, and N 2 = 0.8. and 3D-look-up table (i.e., Butler-Volmer correction voltage) in order to consider their dependences on the operating conditions (i.e., SOC, temperature, and load current). Because for similar operating conditions, the capacity of the Li-ion batteries is different between charging and discharging case (i.e., Coulombic efficiency different than 100%), two separate look-up tables were implemented for the charging and discharging capacity, respectively. Because of the hysteresis effect (see Figure 5), which is inherent to LFP/C battery chemistry, a similar approach was followed for implementing the OCV characteristic; however, for other Li-ion battery chemistries, a single look-up table is enough to model the OCV characteristic for both charging and discharging cases.
During the fitting process of the measured impedance spectra, it was found out that only Rs, R1, R2, Q1, and Q2 are dependent on the SOC and temperature at which the EIS measurements had been performed. Consequently, these parameters were implemented into the performance model as 2D look-up tables since their values are dependent on the SOC and temperature; Figure 19 exemplifies the look-up table implementation of the series resistance Rs. The other parameters of the EEC, Ls, N1, and N2 were found to be constant, independent of the operating conditions: Ls = 1.8·× 10 −8 H, N1 = 0.5, and N2 = 0.8. Finally, the estimated voltage of the LFP/C battery was obtained by summing up the contributions from the OCV, the voltage drop across the EEC, and the Butler-Volmer correction voltage, as illustrated in Figure 18.

Validation of the Performance Model
In order to verify the accuracy of the developed performance model of the LFP/C battery and to validate the proposed characterization procedure, different load current profiles were considered. The accuracy of the developed performance model was evaluated by computing the coefficient of determination R 2 which was obtained by comparing the measured and estimated battery's voltage profiles, according to Equation (8).
where Vmodel is the estimated battery voltage, Vmeas represents the measured battery voltage, x represents the present observation, and n represents the total number of observations. Moreover, the mean ( ε ) and the maximum (εmax) deviation between the measured and estimated battery voltage profiles were computed according to Equations (9) and (10) for the three considered verification cases. Finally, the estimated voltage of the LFP/C battery was obtained by summing up the contributions from the OCV, the voltage drop across the EEC, and the Butler-Volmer correction voltage, as illustrated in Figure 18.

Validation of the Performance Model
In order to verify the accuracy of the developed performance model of the LFP/C battery and to validate the proposed characterization procedure, different load current profiles were considered. The accuracy of the developed performance model was evaluated by computing the coefficient of determination R 2 which was obtained by comparing the measured and estimated battery's voltage profiles, according to Equation (8). (8) where V model is the estimated battery voltage, V meas represents the measured battery voltage, x represents the present observation, and n represents the total number of observations. Moreover, the mean (ε) and the maximum (ε max ) deviation between the measured and estimated battery voltage profiles were computed according to Equations (9) and (10) for the three considered verification cases.

Pulse Discharge
The first validation of the developed performance model was carried out using a discharging current pulse profile. From a fully charged state, the LFP/C battery was discharged with 1C-rate (i.e., 2.5 A) in steps of 10% SOC; between two consecutive discharging current pulses, a relaxation period of 15 min was applied. The measured and estimated voltage profiles of the LFP/C for this verification case are illustrated in Figure 20. The developed performance model was able to estimate very accurately the voltage of the tested LFP/C battery, with a maximum error ε max = 3.66% and a mean error ε = 6.8 mV.

Pulse Discharge
The first validation of the developed performance model was carried out using a discharging current pulse profile. From a fully charged state, the LFP/C battery was discharged with 1C-rate (i.e., 2.5 A) in steps of 10% SOC; between two consecutive discharging current pulses, a relaxation period of 15 min was applied. The measured and estimated voltage profiles of the LFP/C for this verification case are illustrated in Figure 20. The developed performance model was able to estimate very accurately the voltage of the tested LFP/C battery, with a maximum error εmax = 3.66% and a mean error ε = 6.8 mV. The distribution of the voltage deviations, expressed as relative error values, is presented in Figure 21, showing that the model has a slight tendency of over-estimating the voltage behavior of the tested LFP/C battery; the main contribution of this tendency is represented by the voltage estimation error after the current pulse is cut off. The dynamic behavior of the performance model might be improved (thus the accuracy of the model will further increase) by considering an EEC with more ZARC elements (or R-C parallel networks) [24]; however, this will cause a high increase of the computation time of the performance model.

Dynamic Pulse Charging and Discharging
In real applications (e.g., e-mobility, grid support etc.), Li-ion batteries are subjected to complex charging-discharging profiles. Consequently, the developed performance model was verified for a The distribution of the voltage deviations, expressed as relative error values, is presented in Figure 21, showing that the model has a slight tendency of over-estimating the voltage behavior of the tested LFP/C battery; the main contribution of this tendency is represented by the voltage estimation error after the current pulse is cut off. The dynamic behavior of the performance model might be improved (thus the accuracy of the model will further increase) by considering an EEC with more ZARC elements (or R-C parallel networks) [24]; however, this will cause a high increase of the computation time of the performance model.

Pulse Discharge
The first validation of the developed performance model was carried out using a discharging current pulse profile. From a fully charged state, the LFP/C battery was discharged with 1C-rate (i.e., 2.5 A) in steps of 10% SOC; between two consecutive discharging current pulses, a relaxation period of 15 min was applied. The measured and estimated voltage profiles of the LFP/C for this verification case are illustrated in Figure 20. The developed performance model was able to estimate very accurately the voltage of the tested LFP/C battery, with a maximum error εmax = 3.66% and a mean error ε = 6.8 mV. The distribution of the voltage deviations, expressed as relative error values, is presented in Figure 21, showing that the model has a slight tendency of over-estimating the voltage behavior of the tested LFP/C battery; the main contribution of this tendency is represented by the voltage estimation error after the current pulse is cut off. The dynamic behavior of the performance model might be improved (thus the accuracy of the model will further increase) by considering an EEC with more ZARC elements (or R-C parallel networks) [24]; however, this will cause a high increase of the computation time of the performance model.

Dynamic Pulse Charging and Discharging
In real applications (e.g., e-mobility, grid support etc.), Li-ion batteries are subjected to complex charging-discharging profiles. Consequently, the developed performance model was verified for a

Dynamic Pulse Charging and Discharging
In real applications (e.g., e-mobility, grid support etc.), Li-ion batteries are subjected to complex charging-discharging profiles. Consequently, the developed performance model was verified for a current profile, which was composed of charging and discharging pulses of different C-rates, over the whole SOC interval; the considered current profile is presented in Figure 22. To allow for a non-biased verification, the C-rates, which were used for this dynamic profile, are different than the C-rates used to compute the Butler-Volmer correction voltage (see Sections 3.5 and 4.2).
Batteries 2016, 2, 37 15 of 19 current profile, which was composed of charging and discharging pulses of different C-rates, over the whole SOC interval; the considered current profile is presented in Figure 22. To allow for a nonbiased verification, the C-rates, which were used for this dynamic profile, are different than the Crates used to compute the Butler-Volmer correction voltage (see Sections 3.5 and 4.2). The measured and estimated voltage profiles obtained by applying the dynamic current profile, presented in Figure 22, are illustrated in Figure 23. For this case, the proposed model is able to predict the voltage of the LFP/C battery with a mean error of ε = 6.4 mV and maximum error of εmax = 4.89%.
Furthermore, the distribution of the relative error, presented in Figure 24, shows in this case a tendency of the developed model to under-estimate the voltage of the battery; as shown in the zoomed view of Figure 23, this under-estimation of the voltage was mainly caused by the inability of the model to react to sudden and large variation of C-rates (i.e., 3.5C-rate to 1C-rate).   The measured and estimated voltage profiles obtained by applying the dynamic current profile, presented in Figure 22, are illustrated in Figure 23. For this case, the proposed model is able to predict the voltage of the LFP/C battery with a mean error of ε = 6.4 mV and maximum error of ε max = 4.89%. Furthermore, the distribution of the relative error, presented in Figure 24, shows in this case a tendency of the developed model to under-estimate the voltage of the battery; as shown in the zoomed view of Figure 23, this under-estimation of the voltage was mainly caused by the inability of the model to react to sudden and large variation of C-rates (i.e., 3.5C-rate to 1C-rate).
Batteries 2016, 2, 37 15 of 19 current profile, which was composed of charging and discharging pulses of different C-rates, over the whole SOC interval; the considered current profile is presented in Figure 22. To allow for a nonbiased verification, the C-rates, which were used for this dynamic profile, are different than the Crates used to compute the Butler-Volmer correction voltage (see Sections 3.5 and 4.2). The measured and estimated voltage profiles obtained by applying the dynamic current profile, presented in Figure 22, are illustrated in Figure 23. For this case, the proposed model is able to predict the voltage of the LFP/C battery with a mean error of ε = 6.4 mV and maximum error of εmax = 4.89%.
Furthermore, the distribution of the relative error, presented in Figure 24, shows in this case a tendency of the developed model to under-estimate the voltage of the battery; as shown in the zoomed view of Figure 23, this under-estimation of the voltage was mainly caused by the inability of the model to react to sudden and large variation of C-rates (i.e., 3.5C-rate to 1C-rate).    current profile, which was composed of charging and discharging pulses of different C-rates, over the whole SOC interval; the considered current profile is presented in Figure 22. To allow for a nonbiased verification, the C-rates, which were used for this dynamic profile, are different than the Crates used to compute the Butler-Volmer correction voltage (see Sections 3.5 and 4.2). The measured and estimated voltage profiles obtained by applying the dynamic current profile, presented in Figure 22, are illustrated in Figure 23. For this case, the proposed model is able to predict the voltage of the LFP/C battery with a mean error of ε = 6.4 mV and maximum error of εmax = 4.89%.
Furthermore, the distribution of the relative error, presented in Figure 24, shows in this case a tendency of the developed model to under-estimate the voltage of the battery; as shown in the zoomed view of Figure 23, this under-estimation of the voltage was mainly caused by the inability of the model to react to sudden and large variation of C-rates (i.e., 3.5C-rate to 1C-rate).

Field-Measured Profile
The third verification of the developed performance model was performed using a realistic current profile, which was measured on field on LFP/C batteries providing primary frequency regulation (PFR) in the Danish energy market [5]. The considered current profile with a length of two hours is illustrated in Figure 25. The comparison between the estimated and measured voltage of the LFP/C battery, when the field-measured current profile was considered, is presented in Figure 26. The performance model of the LFP/C battery estimates generally with good accuracy (ε = 8.4 mV and ε max = 4.54%) the voltage of the battery. Nevertheless, the model has the tendency to over-estimate the voltage response of the LFP/C battery, when sudden and large variation of the load current have occurred (i.e., changes in the current from below 1C-rate to 4C-rate); this behavior is illustrated in Figure 27, where the distribution of the relative error obtained by comparing the measured and estimated voltage profiles is plotted.

Field-Measured Profile
The third verification of the developed performance model was performed using a realistic current profile, which was measured on field on LFP/C batteries providing primary frequency regulation (PFR) in the Danish energy market [5]. The considered current profile with a length of two hours is illustrated in Figure 25. The comparison between the estimated and measured voltage of the LFP/C battery, when the field-measured current profile was considered, is presented in Figure 26.
The performance model of the LFP/C battery estimates generally with good accuracy ( ε = 8.4 mV and εmax = 4.54%) the voltage of the battery. Nevertheless, the model has the tendency to over-estimate the voltage response of the LFP/C battery, when sudden and large variation of the load current have occurred (i.e., changes in the current from below 1C-rate to 4C-rate); this behavior is illustrated in Figure 27, where the distribution of the relative error obtained by comparing the measured and estimated voltage profiles is plotted.

Field-Measured Profile
The third verification of the developed performance model was performed using a realistic current profile, which was measured on field on LFP/C batteries providing primary frequency regulation (PFR) in the Danish energy market [5]. The considered current profile with a length of two hours is illustrated in Figure 25. The comparison between the estimated and measured voltage of the LFP/C battery, when the field-measured current profile was considered, is presented in Figure 26.
The performance model of the LFP/C battery estimates generally with good accuracy ( ε = 8.4 mV and εmax = 4.54%) the voltage of the battery. Nevertheless, the model has the tendency to over-estimate the voltage response of the LFP/C battery, when sudden and large variation of the load current have occurred (i.e., changes in the current from below 1C-rate to 4C-rate); this behavior is illustrated in Figure 27, where the distribution of the relative error obtained by comparing the measured and estimated voltage profiles is plotted.

Field-Measured Profile
The third verification of the developed performance model was performed using a realistic current profile, which was measured on field on LFP/C batteries providing primary frequency regulation (PFR) in the Danish energy market [5]. The considered current profile with a length of two hours is illustrated in Figure 25. The comparison between the estimated and measured voltage of the LFP/C battery, when the field-measured current profile was considered, is presented in Figure 26.
The performance model of the LFP/C battery estimates generally with good accuracy ( ε = 8.4 mV and εmax = 4.54%) the voltage of the battery. Nevertheless, the model has the tendency to over-estimate the voltage response of the LFP/C battery, when sudden and large variation of the load current have occurred (i.e., changes in the current from below 1C-rate to 4C-rate); this behavior is illustrated in Figure 27, where the distribution of the relative error obtained by comparing the measured and estimated voltage profiles is plotted.   The developed performance model is able to accurately estimate the voltage behavior of the tested LFP/C battery-regardless of the applied current profile-returning a maximum error lower than 5% and a mean error below 8.5 mV as summarized in Table 2.

Conclusions
A seven-step characterization methodology for measuring the performance characteristics of the Li-ion batteries was proposed in the first part of this work. For exemplification of the procedure, a commercially available 2.5 Ah LFP/C battery was used. The results obtained during the characterization test were used to develop and parameterize the performance model of the LFP/C battery.
In order to model the dynamic behavior of the battery, a novel hybrid procedure was proposed. Thus, the developed performance model combines information obtained from the electrochemical impedance spectroscopy (EIS) test (i.e., dependence of the impedance on SOC and temperature) and HPPC test (i.e., dependence on the load current), respectively. The proposed hybrid performance model was verified using different representative dynamic current profiles. The obtained results, maximum error lower than 5% and a mean error below 8.5 mV, have suggested that the developed model is able to estimate with high accuracy the voltage of the LFP/C battery, independent of the considered test conditions. Thus, it can be concluded that the proposed seven-step generalized characterization procedure for Li-ion batteries was validated.
Even though the proposed characterization procedure was illustrated for a certain Li-ion battery chemistry (i.e., LFP/C), it can be applied to any type of Li-ion battery chemistry with the amendment that the values of the test conditions (i.e., temperature and load current levels) have to be adjusted in order to match the manufacturer data-sheets.