Impact of the Temperature in the Evaluation of Battery Performances During Long-Term Cycling — Characterisation and Modelling

This paper presents the results regarding the thermal characterisation and modelling of high energy lithium-ion battery cells at both room (25 ◦C) and cycling (35 ◦C) temperatures. In this work two types of Nickel Manganese Cobalt (NMC) batteries are studied: a fresh (or uncycled) and an aged (or cycled) battery cells. The ageing of the studied NMC battery cells is achieved by means of accelerated ageing tests (i.e., repetition of numerous charge and discharge cycles) at 35 ◦C cycling temperature. Temperature at the surface of the battery cells is characterised, with a set of three discharge current rates 0.3C (i.e., 6 A), 1C (i.e., 20 A) and 2C (i.e., 40 A), and the evolutions at three different locations on the surface of the battery cells namely, at the top, in the center and at the bottom regions are measured. In addition, temperature and ageing dependent electrochemical-thermal modelling of the uncycled and cycled battery cells is also successfully accomplished in case of both room and cycling temperatures. Numerical simulations were carried out in case of high 2C constant current rate, and the assessment of the modelling accuracy by comparison of the predicted battery cells voltage and temperature with respect to the experimental data is further presented. With this paper, thermal performances of battery cells prior and after long-term cycling are evaluated at the cycling temperature, next to the ambient temperature. Hence, thermal characterisation and modelling results are more closely reflecting that encountered by the battery cells in real cycling conditions, so that their performances are believed in this way to be more objectively evaluated.


Introduction
Whether in mobile phones, in portable computers, for starting the car engine or for powering the satellite that sends radio communication signals down to earth [1][2][3], there is no wonder anymore that batteries have become essential components in our daily lives.Particularly, in the transportation sector as part of the daily mobility, electric vehicles (EVs) including hybrid electric vehicles (HEVs), plug-in hybrid electric vehicles (PHEVs) and pure battery electric vehicles (BEVs) are predicted to dominate the vehicle market [4,5].In recent years, Lithium Nickel Manganese Cobalt (NMC) for use in positive electrodes eventually offered improved cycle life, thermal stability and energy density capabilities for lithium-ion batteries [6][7][8][9].In practice, limitations in the lifetime of lithium-ion batteries still reveal often to be an issue in battery applications as in automotive [10][11][12][13].
Throughout their lifetime, batteries are experiencing several kinds of degradation mechanism.When cycled at high temperature, degradation mechanisms leading to capacity loss are observed through the rise of the positive electrode impedance especially at low lithium content.This may be due to the loss of conductive carbon or to a non-ionically conductive organic film at the surface of the electrode [14].In addition, the electrodes expansion associated to the lithium insertion and extraction processes or the external stresses on the battery stack, are inducing mechanical mechanisms such as the viscoelastic creep behavior of the electrochemically inactive separator sheets.These mechanisms are altering the lithium-ion battery cells due to the pore closure effect, which is affecting the ion transport and is leading on the long term to internal resistance increase and capacity fade [15].At the level of the porous electrodes (composed of solid active material particles and liquid electrolyte), diffusion-induced stresses are causing degradations in the active material particles, as demonstrated by non-local damage evolution modelling [16].Life-prognostic methods for the predictions of damages on lithium-ion batteries caused by mechanical fatigue, can further be examined by physics-based modelling including the growth of solid electrolyte interface (SEI) at the carbonaceous negative electrode for instance, as single source of aging [17].
Further in case of advanced energy storage as Lithium metal (LMB) batteries, suitable for high energy density applications, robust SEI is required for improved cycling performances by preventing from the occurrence of dendrite growth.Factors such as the nature of the electrolyte, the temperature or the current density are influencing the SEI formation.Although, the introduction of pre-formed protective layer or modifications in electrode design or in electrolyte composition and additives represent efficient methods to modify the SEI layer [18].Further in solid-state batteries, the expansion induced by the intercalation of Lithium ions in the electrodes active material particles is constrained by the solid-state electrolyte (SE), therefore leading to mechanical degradation.Mechanical degradation of the solid electrolyte (SE) is responsible for worsen transport properties, decreased capacity, and decrease in electrodes active material strength.Electro-chemo-mechanical finite element modelling are developed for this purpose for the study of the corresponding damage evolutions [19,20].In view of this, long-term cycling of battery cells and their decreased performances are extensively studied in the literature [21][22][23][24][25].
For battery-powered vehicles, the operating temperature and the uniformity of the temperature distribution are two important factors at the scale of battery modules in which batteries are embedded [26].For lithium-ion batteries the best operating temperatures are in the range between 25 • C and 40 • C [27].Adiabatic tests in accelerated rate calorimeters demonstrate batteries thermal behaviour during charge and discharge, in abuse conditions such as in case of defective cooling systems in battery modules [28,29].
To avoid the occurrence of an excessive temperature rise [30] leading to the degradation in their lifetime or even worse to the thermal runaway phenomenon, local temperature of lithium-ion batteries should be carefully monitored.
So far although, research is commonly focusing on the evaluation of battery performances at room temperature rather than at the temperature at which they are cycled.However, battery performances evaluated in this manner obviously do not reflect those happening in real cycling conditions under which batteries are used and are further ageing.This is specifically believed to lead to biased and unfair assessments as regards batteries real capabilities.For accelerated ageing tests carried out at a temperature different from the room temperature, a lack of knowledge is hence recognized regarding the impact of the temperature chosen for the evaluation of batteries performances during long-term cycling.Therefore, with the research work achieved and presented in this paper, the impact of the temperature in the assessment of lithium-ion battery cells performances during long-term cycling is aimed to be investigated through characterisations together with numerical modelling for its further prediction.

Studied Batteries
In this research work, high energy density (174 Wh/kg) lithium-ion battery cells (EIG C020, Chuncheongnam-do, South Korea) are studied.These batteries have been theoretically designed to be optimised for Electric (EV) and Plug-In Hybrid (PHEV) applications with a 20 Ah nominal capacity and Nickel Manganese Cobalt-LiNi 0.4 Mn 0.4 Co 0.2 O 2 (NMC) positive electrodes.The battery cell specifications and characteristics are listed in Table 1.
Table 1.Main properties of the high energy lithium polymer EIG C020 battery cells [31].

Long-Term Cycling
For the ageing study presented in this research work, accelerated ageing tests have been performed by cycling of a high energy lithium-ion battery cell at 35 • C .The cycling conditions are defined by the following specifications:

•
The State-of-Charge (SoC) (in %) which is defined by the ratio between the battery cell actual capacity and its nominal capacity, as expressed in Equation ( 1):

•
The Depth-of-Discharge (DoD) (in %) which reflects the percentage of battery capacity that has been discharged expressed as a percentage of its actual maximum discharge capacity and is linked to the SoC in Equation (2):

•
The middle State-of-Charge (mid-SoC) which is the SoC around which the battery cell is cycled.
• Nominal Capacity which is the capacity expressed in Ampere-hours (Ah), specified by the manufacturer at specific conditions.In this research work, the nominal capacity of the studied battery cells is equal to 20 Ah.

•
The Charge/Discharge current rates expressed in terms of C-rate where the C-rate defines the rate at which a battery is discharged or charged relatively to its nominal capacity.

•
The Full Equivalent Cycles (FEC) (in case the cycling is performed with a DoD inferior to 100%), defines the amount of cycles needed to equal the same amount of Charge/Discharge capacity observed in case the cell would have been cycled with 100% DoD.

•
The Capacity Fade (∆C) (in %) indicates the loss of capacity of the battery observed in time by comparison of its actual capacity with its measured initial capacity, as expressed in Equation (3): • The State of Health (SoH) (in %) translates the ageing of the battery cell and is defined with the Equation (4): In this research work, the cycling conditions were defined by 80% Depth of Discharge (DoD), which is defined as 50% mid-State of Charge (SoC) and 1C/1C (i.e., 20 A/20 A) current for the constant charge and discharge currents.
The reference amount of cycles is equal to 100 FEC where each cycle is defined by a repetition of one charge and one discharge sequences.The research work presented in this paper is based on one uncycled (or fresh) cell (i.e., cell 1) and one cycled (or aged) cell (i.e., cell 2) that was cycled up to 900 FEC.Table 2 summarizes the testing of the studied cells and their purpose.Cells 1 and 2 were aimed for thermal characterisations of fresh and aged battery cells, respectively.Capacity tests have been performed on both battery cells at the Begin Of Life (BoL) and as part of the Check-Ups (CUs) every 100 FEC on the cycled cell to monitor its SoH evolution and its capacity fade due to the accelerated ageing.This test is defined by a 3 times repetition of charge and discharge sequences at C/3 (6.666 A) current rate.In this paper, this test is presented at room temperature.
The evolution of the discharge capacity and the State-of-Health (SoH) of the cell 2 with the number of cycles is illustrated in Figure 1.The capacity loss observed in Figure 1 translates for the ageing of the cells caused by their cycling.The SoH evolution of the cells is based on the recording of their capacity at 25 • C every 100 FEC.
At the level of the electrodes, the measured Equilibrium Potential of the uncycled and cycled NMC positive electrodes are shown in in Figure 2. The Equilibrium Potential of the electrodes dependent on the stoichiometry are in agreement with [22,32,33].Stoichiometry is defined by the amount of Li + intercalated in an electrode over the maximum amount of Li + that can be stored that electrode.As in Figure 2, the stoichiometry varies between 0 and 1 representing the x in Li x Ni 0.4 Mn 0.4 Co 0.2 .The highest Equilibrium Potential value corresponds to the lowest stoichiometry (0) of the electrode.With increasing stoichiometry, the Equilibrium Potential of the uncycled and cycled electrodes are decreasing down to 3 V during the lithiation of the electrodes.Reversely, as the delithiation of the electrodes proceeds the Equilibrium Potential of the electrodes are increasing up to 4.15 V.After long-term cycling, both loss of active material and loss of cyclable lithium are affecting the stoichiometries in the battery cells electrodes [34].NMC electrodes exhibit lower potential near the maximum stoichiometry in the end of the lithiation process in accordance with [22,32,33].This reflects for a faster lithiation process of cycled electrodes due to ageing.Reversely, higher potential near the minimum stoichiometry is observed for cycled electrodes due to their faster delithiation process.This behaviour leads to an increase in the hysteresis (difference between the Equilibrium Potential associated to the lithiation and the delithiation processes) of the cycled electrodes compared to uncycled electrodes.In Figure 3 (left), the measured discharge capacity of the uncycled and cycled NMC electrodes are illustrated.The cycled electrode exhibit a discharge capacity equal to 5.67 mAh when the initially uncycled electrode show a capacity equal to 5.87 mAh.Conformingly to what expected, the cycled NMC electrode is discharging faster, leading to a lower discharge capacity than prior long-term cycling.The analysis of the incremental discharge capacity (dQ/dV) of both electrodes in Figure 3 (right) confirms the impact of the long-term cycling on the performances of the initially uncycled electrodes.Indeed, a decrease in the intensity of the peak in the dQ/dV curve is observed for the cycled electrode as observed in [22,32].
At the level of the negative electrodes, in Figure 4, the measured Equilibrium Potential of the uncycled and cycled Graphite electrodes are presented.As previously observed for the positive electrode and in agreement with the literature [33,[35][36][37][38], the Equilibrium Potential of both uncycled and cycled electrodes are dependent on their stoichiometry.After long-term cycling, Graphite electrodes exhibit lower potential near the maximum stoichiometry in the end of the delithiation process.This reflects for a faster lithiation process of cycled electrodes due to ageing.Reversely, higher potential near the minimum stoichiometry is observed for cycled electrodes due to their faster delithiation process.These observations are in line with the behaviour of NMC electrodes measured after long-term cycling.Although, the hysteresis associated to the uncycled and cycled electrodes during their lithiation and delithiation processes is more marked for Graphite electrodes compared to NMC electrodes.This indicates for more significant ageing shown by the Graphite electrodes.In Figure 5 (left), the measured discharge capacity of the uncycled and cycled Graphite electrodes are illustrated.The cycled electrode exhibit a discharge capacity equal to 7.07 mAh when the initially uncycled electrode show a capacity equal to 7.34 mAh.Similarly to what observed for the NMC electrodes, the cycled Graphite electrode is discharging faster, leading to a lower discharge capacity than prior long-term cycling.Although, the difference in the discharge capacity of the uncycled and cycled electrodes is more strong for Graphite electrodes compared to NMC electrodes.This also indicates for more significant ageing shown by the Graphite electrodes.The incremental discharge capacity (dQ/dV) curves of both electrodes, as shown in Figure 5 (right), are exhibiting three peaks (two majors and one minor) in accordance with the literature [35][36][37].The analysis of the incremental discharge capacity (dQ/dV) of both electrodes in Figure 5 (right) confirms the impact of the long-term cycling on the performances of the initially uncycled electrodes.Indeed, the peaks in the dQ/dV curve of the cycled electrode exhibit lower intensity than observed for the uncycled electrode.In addition, the incremental capacity curve is seen to shift towards lower potentials.This might be explained by the faster lithiation process of the cycled compared to the uncycled electrodes as a consequence of the long-term cycling (in agreement with what suggested by Figure 4).

Thermal Characterisation
Thermal characterisation is performed at the scale of the whole fresh and aged pouch battery cells, by recording their temperature at three different locations namely at the top, middle (or center) and bottom regions as demonstrated in Figure 6.
Thermal behaviour of the battery cells basically results from the overall heat generation in the different layers inside the cells, as from the rate at which this heat can be transfered to the surroundings.It is assumed that the temperature distribution is independent of the z-direction, along the thickness of the battery cells.The thickness of the battery cells in the z-direction is fundamentally small compared to its two other dimensions (width and length) in x and y axes.Hence, the heat is propagating quickly by conduction along the thickness of the studied cells, so that the temperature distribution can be seen as uniform along the z-direction.

Numerical Simulation
In this work, numerical simulations are carried out according to an ageing and temperature dependent electrochemical-thermal modelling of the battery cells.The electrochemical simulations are performed with a pseudo two dimensional (P2D) modelling [39,40].Wherein the electrodes are assumed to be porous materials composed of spherical particles immersed in the electrolyte.As part of physico-chemical models, electrochemical models are well documented in the literature [41][42][43].These models are capable of accurately predicting the performance of lithium-ion batteries, which allows the battery manufacturers to optimise and improve the cell design and materials properties chosen for their products.With its wide set of physical parameters, porous electrochemical models are useful tools for design purposes of battery cells since they allow predicting and understanding the physical processes occurring inside the battery during its lifetime.
For the electrochemical modelling the internal structure of the battery is considered to be composed of multiple layers (electrodes, separator, current collectors) which is reduced (using a superposition approach) to a single unit layer made of negative current collector, negative electrode, separator, positive electrode and positive current collector.
The coupling between the electrochemical (EC) and the thermal (T) parts of the developed temperature and ageing dependent electrochemical-thermal (ECT) modelling of the studied battery cells, is illustrated in Figure 7.The heat equation which describes the heat transfer phenomenon occurring between the cell and its environment, for the coupling between both electrochemical and thermal modelling parts, is given in Equation ( 5): where ρ represents the density, c p the heat capacity, T the temperature, k the thermal conductivity, and Q represents the heat generation.
The heat generation Q inside the cell is derived from Equation ( 6): is the effective value of the electrical conductivity of the electrolyte phase and c e represents the concentration of the electrolyte [44].The first term in Equation ( 6) defines the reaction heat or the irreversible heat generation, the second term defines the reversible heat generation or active heat generation, while the third one (in brackets) defines the ohmic heat generation [45,46].By definition, the reaction and the reversible heat generations are caused by the negative and the positive electrodes, while the ohmic heat generation comes from the electrolyte and both negative and positive electrodes.
The resulting thermo-physical properties of the whole battery cell after computation, according to Equations (A1, A2, A3 and A4) and Table A4, are summarised in Table 3.
Table 3. Thermo-physical properties of the studied Nickel Manganese Cobalt(NMC)/Graphite type of battery cell.

Thermal Parameter Value
Unit Tables A1 and A2 list the electrochemical parameters values associated to the uncycled and cycled porous electrodes, respectively, while Table A3 lists the electrochemical parameters values associated to the electrolyte and the separator of the studied battery cells.
Figure 8 illustrates the different steps carried out in this research work as part of the electrochemical-thermal modelling of the fresh and aged battery cells at both room and cycling temperature.First, the domains and the geometry related to the battery cells are defined, and further associated to their respective materials properties.After discretisation by meshing of the domains, the physics and the boundary conditions associated to each domain are defined for the numerical solving of the governing equations.Last, as output of the numerical simulations, the predicted voltage and temperature of the battery cells are compared against experimental data for validation of the electrochemical-thermal modelling.
The commercial finite element software COMSOL Multiphysics R 4.3b was used for solving of the governing electrochemical and thermal differential equations.The temperature and ageing dependent modelling of the studied battery cells was achieved via the COMSOL R (Stockholm, Sweden) LiveLink TM for MATLAB R interface.A fully coupled MUMPS direct solver was opted as linear solver with a minimum step tolerance equal to 1 × 10 −4 .The mesh resolution for the discretisation of the 3D battery cells and the 1D equivalent electrochemical layer domains, consist of a total of 2163 tetrahedral elements of minimum size equal to 1.72 × 10 −2 and 158 elements of 1.57 × 10 −6 minimum size, respectively.

Thermal Response-Room and Cycling Temperatures
The thermal response in terms of temperature increase of the cell 1 measured during the discharge process under various current rates 6 A, 20 A, 40 A (corresponding to 0.3C, 1C, 2C) at 25 • C and 35 • C at the top, middle and bottom regions are presented in Figure 10.The measurements were performed in a climate chamber (Clima Temperatur Systeme GmbH (CTS), Hechingen, Germany) in which the cell is subjected to forced convection.The highest temperature rises are established in the middle and bottom regions of the cell, with up to 25% increase.Although for low current rates similar evolutions are observed for both studied temperatures in the three regions, at high current rates more significant temperature rises are observed at room in comparison to the the cycling temperature.The differences in the temperature increases at the three top, middle and bottom regions indicate a non-uniform temperature distribution on the surface of the fresh cell.The temperature distribution originates from the non-uniform potential developed in the electrodes sheets leading to non-uniform current density.The heat generation being proportional to the square root of the current, is hence non-uniform and leads in turn to non-uniform temperature distribution at the surface of the battery cell.Further, differences in the temperature increase of an uncycled battery cell under room compared to the cycling temperature in Figure 10 are around 10%.

Thermal Response-Room and Cycling Temperatures
The thermal response in terms of temperature increase of the cell 2 measured during the discharge process under various current rates 6 A, 20 A, 40 A (corresponding to 0.3C, 1C, 2C) at 25 • C and 35 • C at the top, middle and bottom regions are presented in Figure 12.The measurements were performed in a climate chamber (Clima Temperatur Systeme GmbH (CTS), Hechingen, Germany) in which the cell is subjected to forced convection.Similarly to the fresh cell, the difference observed in the temperature increases at the three top, middle and bottom regions on the surface of the aged battery cell indicate for a non-uniform temperature distribution.The highest temperature rises are detected in the middle and bottom regions of the cell, with up to more than 30% increase under room temperature.The evolution of the temperature in the three regions is similar at low current rates, although more significant temperature rises are observed at high current rates under room in comparison to the the cycling temperature.These differences in temperature increases initially observed in Figure 10, are seen now in Figure 12 to be even more prominent due to the ageing of the battery cell caused by long-term cycling.
These results stress the importance of a good thermal monitoring of cycled battery cells, especially in real applications when embedded in battery modules and packs.Thermal runaway phenomenon at the scale of a battery pack is in practice often initiated by the thermal failure of a single cell, which in turns causes other cells to fail.Within a pack composed of aged cells, this phenomenon is even more likely to occur and hence critical.The strong difference up to 15% (Figure 12) in the temperature increase of a cycled battery cell under room compared to the cycling temperature, emphasises for the issues of possible mismatch in the thermal behaviour of battery cells at the level of a battery pack not being initially appropriately conditioned.

Cells Voltage and Thermal Responses Comparisons
The comparison between the voltage response of the fresh and aged cells at 25 • C and 35 • C are illustrated in Figure 13, respectively.A decrease in capacity of the aged cell while discharged is observed at both room and cycling temperature.Considering the long-term cycling under 35 • C , the capacity fade observed can be attributed to the side reactions inducing further an important loss of usable Li + and active material particles participating in the electrochemical lithiation/delithiation reactions inside the cell.The comparison between the temperature increases of the fresh and aged cells at 25 • C and 35 • C are illustrated in Figures 14 and 15, respectively.These figures highlight the impact of the choice of the temperature for assessing the thermal responses of fresh and aged cells respectively, when discharged with various currents.At both room and cycling temperatures, the aged cell exhibits more pronounced temperature rises compared to the fresh cell.After long-term cycling, the voltage of the aged cell tends to be lowered and the internal resistance tends to increase, which leads to greater heat generation and hence greater temperature rises.Further, Figures 14 and 15 reveal that when assessed at the same temperature (respectively at room or cycling), the evolutions of the temperature at the surface of the battery cell prior in contrast to after long-term cycling are although not significantly differing.These differences are less in comparison to the differences observed in the thermal behaviour of the cells when assessed at different temperatures (either room/cycling or cycling/room).At both room and cycling temperatures, the maximum difference between the fresh and aged cells surface temperature evolutions being round 5%.In the long term, results presented in Sections 3.1 and 3.2 suggest to carry out in a systematic manner, the evaluation of batteries performances at the cycling temperature.According to this, in real applications as in automotive, it can be recommended to use the track of the temperature within the battery pack (e.g., mean) recorded by the battery management system, to define the temperature at which battery cells performances should be evaluated.In addition, room temperature can be used if needed for the comparison of battery cells performances, of same type that are though not operating at the same temperature.

Voltage Predictions
As outcome of the numerical simulations, the predicted voltage response of the cell 1 (or uncycled cell) at room and cycling temperatures are illustrated in Figures 16 and 17 during high constant current discharge (2C i.e., 40 A).From the comparison with the experimental measurements, the absolute value of the voltage deviations are overall less than 2%.Particularly, at the beginning and at the end of the discharge curves, highest voltage deviations can be found due to steeper slopes presented in the voltage behaviour of the battery cell.The good agreement in the simulated and the experimental voltages observed in Figures 16 and 17 translates for the accurate electrochemical characterisation as for the effective implementation of the parameters estimation in the framework of the coupled electrochemical-thermal modelling at both room and cycling temperatures.
The predicted voltage response of cell 2 (or cycled cell) at room and cycling temperatures are illustrated in Figures 18 and 19 during high constant current discharge (2C i.e., 40 A).From the comparison with the experimental measurements, the absolute value of the voltage deviations are overall less than 2% and 3%, respectively.Particularly, at the beginning and at the end of the discharge curve, highest voltage deviations can be found due to steeper slopes presented in the voltage behaviour of the battery cell.The good agreement in the simulated and the experimental voltages observed in Figures 18 and 19 translates for the accurate electrochemical characterisation as for the effective implementation of the parameters estimation in the framework of the coupled electrochemical-thermal modelling at both room and cycling temperatures.The slightly longer simulation time observed for the voltage curves in Figure 18, highlight the effect of the temperature on the internal battery cell electrochemical processes and further the non-negligible role of the temperature in the assessment of the battery performances during long-term cycling.

Temperature Predictions
As outcome of the numerical simulations, the predicted temperature behaviour of cell 1 (or uncycled cell) at room and cycling temperatures are illustrated in Figures 20 and 21 during high constant current discharge (2C i.e., 40 A).From the comparison with the experimental measurements at three different locations (Top, Middle and Bottom regions), the absolute value of the temperature deviations are overall less than 1 and 2 • C, respectively.The good agreement in the simulated and the experimental data observed in Figures 20 and 21 translates for the accurate electrochemical-thermal characterisation and for the effective temperature predictions delivered by the coupled electrochemical-thermal modelling of the studied battery cells at room temperature.The temperature differences in the three regions indicate for non-uniform temperature distribution taking place at the surface of the battery cell.Further, these results emphasise the need for appropriate thermal management solutions that would on the long term decrease the battery ageing, by preventing the development of temperature differences at the surface of the battery cell (i.e., hot and cold regions).
The predicted temperature behaviour of a cycled cell at room and cycling temperatures are illustrated in Figures 22 and 23 during high constant current discharge (2C i.e., 40 A).Predictions of the temperature distribution under 35 • C are illustrated in Figure 24.With these results, the higher temperature observed at the surface of the cycled battery cell, may be attributed to the rise in internal resistance, and hence in the heat generation inside the cell due to the ageing of the cell caused by long-term cycling.From the comparison with the experimental measurements at three different locations (Top, Middle and Bottom regions), the absolute value of the temperature deviations are overall less than 1.5 and 2 • C, respectively.In view of the results regarding cell 2 (or cycled battery) cell, less temperature uniformity is observed in Figure 23 compared to the case of an uncycled cell in Figure 21.The latter observation suggest for differences in the internal resistance of the cell in the different regions.Further thermal investigations might confirm this behaviour.The good agreement in the simulated and the experimental data observed in Figures 22 and 23 translates for the accurate electrochemical-thermal characterisation and for the effective temperature predictions delivered by the coupled electrochemical-thermal modelling of the studied battery cells at room temperature.

Conclusions
With this paper, thermal characterisation of both uncycled and cycled battery cells, is originally performed and investigated at the cycling temperature (or 35 • C), beyond the commonly studied 25 • C room temperature.By considering the cycling temperature, thermal characterisations of the battery cells are seen to more closely reflect that encountered in real cycling conditions, leading to a more objective assessment of the battery cells performances.
Experimental investigations of the temperature at the surface of the uncycled and cycled battery cells at various current rates revealed the impact of ageing, with higher temperature increases observed at the cycled battery cell surface.Temperature evolution at three different locations on the surface of the fresh battery cell namely, top, center and bottom regions also indicate for non-uniformity of the temperature distribution, which phenomenon is further accentuated for cycled battery cells that experienced long-term cycling.
Additionally, temperature and ageing dependent electrochemical-thermal modelling of the uncycled and cycled battery cells at both room and cycling temperatures were successfully achieved for high constant discharge current rates (2C i.e., 40 A), with absolute deviations less than 3% (voltage predictions) and 2 • C (temperature predictions).
In the long term, results presented with this paper suggest to carry out in a systematic manner, the evaluation of batteries performances at the cycling temperature.According to this, in real applications as in automotive, it can be recommended to use the track of the temperature within the battery pack (e.g., mean) recorded by the battery management system, to define the temperature at which battery cells performances should be evaluated.In addition, room temperature can be used if needed for the comparison of battery cells performances, of same type that are though not operating at the same temperature.With the subscripts (a),(b),(c) in Table A2 denoting for the origin of the parameters as follows: (a) experimentally determined; (b) literature and (c) estimated.

Parameters
In Tables A1 and A2, for the experimental determination of the parameters associated to the superscript (a), both NMC and Graphite electrodes sheets were recovered from the uncycled and cycled battery cells according to the procedure detailed in [47].
The electrodes thickness and the current collectors thickness were measured by means of a thickness gauge (i.e., elcometer).
Solid phase electrode potential of both Graphite and NMC electrodes were measured by charging first 18 mm diameter electrodes punches (in EL cells containers (Hamburg, Germany)) to their upper cut-off potential using a CC-CV method (C/30-0.5 C) (as referred in [47]) and discharging them subsequently to their lower cut-off voltage with C/30 (ca.250 µA).
Maximum lithium ions concentration in both electrodes were determined from XRD (X-ray Diffraction) experiments as referred by the procedure in [47].
Radius of active material (AM) particles in both electrodes were determined from the analysis of SEM (Scanning Electron Microscopy) micrographs as referred in [47].
Li + diffusion coefficient in the solid phase in both electrodes were determined from GITT (Galvanostatic Intermittent Titration Technique) according to the procedure detailed in [47].
The activation Energy associated to the Li + diffusion coefficient in the electrodes was determined during the discharge of the electrodes at each 10% discharge step, by means of an incremental temperature profile between 20 • C and 30 • C including 5 steps, each of 2.5 • C increment.
The volume fraction of active material (AM) and electrolyte in both electrodes, the electrochemical reaction rate constant associated to both electrodes, as well as the electronic conductivity of the NMC electrode were estimated.The estimation of these parameters was performed using a non linear least square based algorithm, which was implemented via MATLAB R as referred in [50].
The remaining parameters were referred directly to the literature as indicated by the corresponding reference(s) next to the superscript (b).
Table A3 lists the electrochemical parameters values associated to the electrolyte and the separator of the studied battery cells.

Parameters Value
Thickness separator 24 × 10 −6 m (a) Porosity separator 0.42 (b) , [51] Li + diffusion coefficient in the electrolyte 1.5 × 10 −10 m 2 /s (b) , [51] Li + initial concentration in the electrolyte 2000 mol/m 3 (c) Transference number of electrolyte for Li + diffusion 0.38 (b) , [52] With the subscripts (a),(b),(c) in Table A3 denoting for the origin of the parameters as follows: (a) experimentally determined; (b) literature and (c) estimated.The separator thickness was measured by means of a thickness gauge (i.e., elcometer).The Li + initial concentration in the electrolyte was estimated using a non linear least square based algorithm which was implemented via MATLAB R as referred in [50].
The thermo-physical properties of the electrodes of the studied battery cells are listed in Table A4.The reported thermo-physical electrodes properties are obtained for similar electrodes in the literature [49].The thermo-physical properties of the whole batteries are computed from those of their constitutive elements (electrodes, current collectors and separator).
The thermal capacity of each studied battery cell is defined as follows: The internal density of each studied battery cell is defined as follows:

Figure 1 .Figure 2 .
Figure 1.Evolution of the discharge capacity and the State-of-Health (SoH) of the cell 2 with the number of cycles.

Figure 6 .
Figure 6.Battery cell schematic and regions for disposition of temperature sensors.

Figure 7 .
Figure 7. Representation of the coupling between the Electrochemical and the Thermal modelling parts.EC: Electrochemical; ECT: electrochemical-thermal.

s
∇φ s ∇φ s + σ e f f e ∇φ e ∇φ e + κ e f f D ∇ln(c e )∇φ e ]ohmic heat(6) where j Li is the transfer current due to the intercalation or deintercalation of lithium ions, φ s is the potential of the solid phase, φ e is the potential of the electrolyte phase, U is the open circuit voltage, T is the temperature, σ e f f s is the effective value of the electrical conductivity of the solid phase, σ e f f e

3. 1 . 1 .
Voltage Response-Room and Cycling TemperaturesThe voltage response of the fresh battery cell (or cell 1) measured during the discharge processes under various current rates 6 A, 20 A, 40 A (corresponding to 0.3C, 1C, 2C) at 25 • C and 35 • C are presented in Figure9.The voltage responses in Figure9highlight an increase in capacity of the fresh cell while discharged at the cycling temperature.The energy barrier required for the electrochemical reaction linked to the insertion and extraction processes of Li + in the electrodes is reduced as the temperature increases.As a result, a greater amount of Li + can be displaced, resulting in greater discharge capacities exhibited by the battery cells.

Figure 9 .
Figure 9. Discharge capacity at 25 • C and 35 • C of the cell 1 at different current rates.

Figure 10 .
Figure 10.Temperature increase of the cell 1 at 25 • C and 35 • C under various current rates.
3.2.Experimental Characterisation-Aged Cell3.2.1.Voltage Response-Room and Cycling Temperatures The voltage response of the aged battery cell (or cell 2) measured during the discharge process under various current rates 6 A, 20 A, 40 A (corresponding to 0.3C, 1C, 2C) at 25 • C and 35 • C are presented in Figure 11.Similarly to what observed for the fresh cell, the voltage responses of the aged cell in Figure 11 highlight an increase in capacity while discharged at the cycling temperature.The energy barrier required for the electrochemical reaction to occur is reduced as the temperature increase.With a lower energy barrier associated to the insertion and extraction processes of Li + in the electrodes, a greater amount of Li + can be displaced, resulting in greater discharge capacities of the battery cells.

Figure 11 .
Figure 11.Discharge capacity at 25 • C and 35 • C of the cell 2 at various current rates.

Figure 12 .
Figure 12.Temperature increase of the cell 2 at 25 • C and 35 • C under various current rates.

Figure 13 .
Figure 13.Discharge capacity at 25 • C and 35 • C of the fresh and aged cells at various current rates.

Figure 14 .
Figure 14.Temperature increase at 25 • C for fresh and aged cells.

Figure 15 .
Figure 15.Temperature increase at 35 • C for fresh and aged cells.

Figure 16 .
Figure 16.(Upper) Uncycled cell voltage during high constant current discharge at 25 • C room temperature; (Lower) Absolute value of the voltage deviation between experimental and simulations data.

Figure 17 .
Figure 17.(Upper) Uncycled cell voltage during high constant current discharge at 35 • C cycling temperature; (Lower) Absolute value of the voltage deviation between experimental and simulations data.

Figure 18 .
Figure 18.(Upper) Cycled cell voltage during high constant current discharge at 25 • C room temperature; (Lower) Absolute value of the voltage deviation between experimental and simulations data.

Figure 19 .
Figure 19.(Upper) Cycled cell voltage during high constant current discharge at 35 • C cycling temperature; (Lower) Absolute value of the voltage deviation between experimental and simulations data.

Figure 20 .
Figure 20.(Upper) Uncycled cell surface temperature at three different locations (Top, Middle, Bottom) during high constant current discharge at 25 • C room temperature; (Lower) Corresponding absolute values of the temperature deviation between experimental and simulations data.

Figure 21 .
Figure 21.(Upper) Uncycled cell surface temperature at three different locations (Top, Middle, Bottom) during high constant current discharge at 35 • C cycling temperature; (Lower) Corresponding absolute value of the temperature deviation between experimental and simulations data.

Figure 22 .
Figure 22. (Upper) Cycled cell surface temperature at three different locations (Top, Middle, Bottom) during high constant current discharge at 25 • C room temperature; (Lower) Corresponding absolute values of the temperature deviation between experimental and simulations data.

Figure 23 .
Figure 23.(Upper) Cycled cell surface temperature at three different locations (Top, Middle, Bottom) during high constant current discharge at 35 • C cycling temperature; (Lower) Corresponding absolute value of the temperature deviation between experimental and simulations data.

Figure 24 .
Figure 24.Predicted temperature distribution at the surface of a cycled cell in the end of a high constant current discharge under 35 • C .

Table 2 .
Cells testings and purposes.

Table A2
lists the electrochemical parameters values associated to the cycled porous electrodes.

Table A4 .
Thermo-physical properties of the electrodes, separators, current collectors of the NMC/ Graphite type of battery cell.