Comparative Study of Surface Temperature Behavior of Commercial Li-Ion Pouch Cells of Different Chemistries and Capacities by Infrared Thermography

The non-uniform surface temperature distribution of a battery cell results from complex reactions inside the cell and makes efficient thermal management a challenging task. This experimental work attempts to determine the evolution of surface temperature distribution of three pouch type commercial cells: Nickel Manganese Cobalt oxide (NMC)-based 20 Ah cell, Lithium Iron Phosphate (LFP) 14 Ah, and Lithium Titanate Oxide (LTO) 5 Ah battery cell by using contact thermistor and infrared (IR) thermography. High current (up to 100 A) continuous charge/discharge and high current (80 A) micro pulse cycling profile were applied on the cells. It was found that thermistor based temperature profile varied cell to cell, especially the LTO cell. Among the investigated cells, the NMC cell shows highest temperature rise and the LTO cell the lowest rise. IR (Infrared) images revealed the spatial distribution of surface temperature, in particular the location of the hottest region varies depending not only on the geometrical and material properties of the cell, but also the type of loads applied on the cells. Finally, a modeling perspective of the cell temperature non-uniformity is also discussed. OPEN ACCESS Energies 2015, 8 8176


Introduction
As an alternative to gasoline fuel, rechargeable energy storage systems (RESS) have garnered the attention of automotive industries since the last decade.A growing number of electric vehicles (EV) and hybrid electric vehicles (HEV) in the present market depicts the rapid growing demand [1][2][3].Being of superior performance in terms of high power and energy density, EV and HEV manufacturers largely depends on lithium-ion batteries as the rechargeable energy source nowadays.However, automotive application requires batteries with high performance as well as with a high level of safety.For instance, heat generated within the cell is one of the major concerns for safety and performance of the battery cell.Total heat generation within the cell can generally be attributed to the contributions of reversible heat (entropic heat components related to the electrochemical reactions) and irreversible heat (ohmic and polarization resistance heat components).On the one hand, high cell temperature, even locally, can have a detrimental influence on the cell performance.For instance, cycle life of the cell can be heavily influenced leading to accelerated ageing of the cell [4,5].On the other hand, excessive and uncontrolled heat can lead to thermal runway and ultimately may lead to explosion of the cell [6].Therefore, monitoring and optimizing the operating temperature of the cell during operation are of crucial importance.
It has been well established that modeling of the thermal behavior of battery cells can play a vital role in cell temperature monitoring and can also provide scopes for the development of battery cells for better thermal management and improved thermal safety of, not only a single battery cell, but also a battery pack consisting several cells [7,8].In reality, a precise validation process of these thermal models can be complex and extremely difficult.Validation through comparing cell surface temperature obtained from model and physical measurements is, comparatively, a less complex and faster method of model validation [9].However, due to the complexity of electrochemical reactions (e.g., effect of side reactions) and uneven current density distribution inside the cell, the spatial distribution of cell surface temperature can be non-uniform [10][11][12][13][14][15], in particular, for large format pouch cells for automotive applications.Moreover, the pattern of the heat distribution, e.g., location of maximum and average temperature, can be different depending on the type of applied load profile [15].This particular behavior of temperature non-uniformity makes the thermal management of a battery pack more challenging in order to maintain a uniform temperature inside the battery pack during operation.Improvements have to be made in battery cells design, and the location oriented active/passive cooling system of battery pack has to be incorporated to overcome this challenge [16][17][18].Now, among the different methods of measuring cell surface temperature, the single point (minute area) measurement (e.g., contact thermocouple, thermistor, etc.) is not adequate to measure the spatially non-uniform temperature distribution.In this case, thermal infrared (IR) imaging is a potential tool that can be used to observe and measure this spatially non-uniform temperature distribution with reasonable high accuracy.Additionally, battery cell surface temperature distribution information obtained as such can be used for monitoring surface temperature locations of interest (maximum and/or average temperature) through single point measurements for a large number of cells during performance-characterization testing.Despite the high importance of non-uniform battery cell surface temperature measurement, to the authors' knowledge, very few works have been published dedicated to this topic using IR thermography [12,14,19].Veth et al. have analyzed surface temperature evolution through IR thermography in their work [14].However, the work was performed on an NMC cell only and the thermal behavior was studied only under a continuous charge and discharge load profile.
In this experimental work, non-uniform surface temperature evolution of three commercial pouch cells, of different chemistries and capacities, were studied through both infrared thermography and contact thermistor measurements.Continuous charge and discharge current, ranging from 10 A up to 100 A, and pulse current up to 80 A, were applied to determine the evolution of surface temperature distribution.

Experimental Procedures
The chemistries of three lithium-ion commercial pouch type battery cells that were used are: Lithium Nickel Manganese Cobalt oxide (NMC), Lithium Iron Phosphate (LFP), and Lithium Titanate Oxide (LTO).The nominal capacities of these cells are 20 Ah, 14 Ah, and 5 Ah, respectively.Two types of load profiles were applied on the cells.One type consisted of several consecutive continuous complete charge (constant current, CC) and discharge (CC) cycles with currents of 10 A, 20 A, 40 A, 60A, 80 A, up to 100 A. Sufficient rest time (1-3 h depending on the current rate) was incorporated between successive operations (charge/discharge) and cycles in order to allow the cell surface temperature to return to ambient temperature.The complete charge and discharge were performed according to the conditions presented in Table 1.On the other type of load profile, high-current micro-pulse cycling was incorporated.The micro-pulse cycling load profile consisted of ~1000 cycles of charge and discharge current pulses.Each of these cycles sequentially consisted of a charge pulse for 2 s at 80 A, 2 s rest (zero current), a discharge pulse for 2 s at 80 A, and 2 s rest.All the tests were performed by using an ACT 0550 (80 channels) battery tester (PEC ® , Leuven, Belgium).Surface temperature was logged by using an NTC 5K thermistor (EPCOS ® , Munich, Germany) with a measuring range of 0-100 °C and a tolerance of 1%.Simultaneously, IR images were captured by a Ti25 thermal imager (FLUKE ® , Everett, WA, USA) at regular time intervals.In order to achieve accurate results from the IR thermography, cells were placed in a semi-closed (to observe behavior at natural heat transfer) and dark environment (to avoid visible light interference).Moreover, the cell surfaces were painted uniformly with a dull black paint in order to minimize the effect of the visible light reflections.The IR camera was calibrated within the range 0-120 °C, with a maximum error of 2% by the manufacturer.The thermal sensitivity, noise-equivalent temperature difference (NETD) was ≤0.09 °C at a 30 °C target temperature (≤90 mK).For the calibration of the experiment, emissivity correction was performed by measuring a known uniform surface temperature painted with the same dull black paint as used on the cell surface.Emissivity was set to ~0.99.A comparative test was carried out to compare the measurement taken by the thermistor and the IR camera during a continuous discharge at 100 A of a NMC cell.The cell surface was divided into 16 equal rectangular portions and two points were selected at two points of intersection of the dividing lines, as shown in Figure 1a,b, which shows the temperature profile of the corresponding two positions on the cell surface.It was found that, throughout the process, the maximum difference between the measurements made by the thermistor and the IR camera was ~1 °C.By means of this comparative test measurement, the accuracy of the IR temperature measurement has been quantified and validated over the complete measurement range.

Cell Characteristics
Internal resistance values of the three cells have been measured by means of the standardized hybrid pulse power capability (HPPC) test.The pulse currents for this HPPC test were adapted to be the same as the currents used for surface temperature measurements, in particular 20, 40, 60, 80, and 100 A. The pulses were performed at three different SoC levels, namely at an 80%, 50%, and 20% SoC level.Figure 2 shows the internal resistance trend for all three cells.From these graphs, we can see that in the case of discharging, and for all three chemistries, higher internal resistance is observed at a low state of charge of the cell.Further, the lowest internal resistance is observed at a high state of charge, in combination with the higher current values.In the case of the LTO chemistry charging pulses, the lowest internal resistance is found at the mid value of the state of charge.

Continuous Charge and Discharge
Figure 3 shows the surface temperature profile of the 20 Ah NMC cell obtained by a thermistor placed near the positive tab during complete charge at 20 A and complete discharge at 10 A, 20 A, 40 A, 60 A, 80 A and 100 A respectively.The temperature profile according to this thermistor measurement during charging showed an initial rise until the state of charge (SoC) level reaches ca.50%.Following, the temperature remains fairly steady until it reaches a 100% SoC level (4.2 V).Discharge at different rates showed, comparatively, a different trend in the temperature profile.For instance, during discharge at 10 A, the temperature rose until a ca.70% SoC level was reached and was followed by a drop until it rose sharply again after it reached a ca.30% SoC level.In order to explain the phenomenon of temperature drop, one may consider the relative dominance of reversible and irreversible heat contributions.It was found that endothermic entropy change is the result of phase change in the electrode material at a certain SoC level range [20].Both of the electrodes may undergo volume and phase changes.At the cathode, these changes are based on the ratio of lithium and other elements (e.g., Cobalt) at the cathode.Additionally, at the anode, these changes are based on the ratio of lithium and carbon during intercalation/deintercalation of lithium [20][21][22].However, at a higher current rate, the contribution of irreversible heat (i.e., polarization resistance heat and ohmic resistance heat) becomes dominant.Therefore, temperature drop due to entropy change became comparatively less significant at 20 A discharge.At higher current of 40 A and 60 A, this effect became trivial, shown by a slight change in the steepness (slope) of the temperature profile.At 80 A and 100 A the effect is negligible.
The IR images of Figure 4 depict the spatial distribution of the cell surface temperature during a discharge at 100 A. Figure 4a shows the schematic of the relative size and orientation of the tabs of the NMC cell.The black border represents the extruded pouch of the cell.The following IR images were represented according to the orientation of the schematic.Now, it is clear from the temperature distribution pattern of the IR images that, initially, the most heated regions were at the adjacent areas of the tabs of the cell and slightly higher near the positive tab.This can be attributed to the comparatively higher resistance of the aluminum positive tab and current collector (compared with the copper negative tab).However, this finding is inconsistent with the findings of Veth et al., who observed that the maximum temperature was initially near the negative tab [14].Apparently, this inconsistency appeared because of the variation in commercial battery cell design (e.g., surface area of the tabs and current collectors, etc.) [16].Nonetheless, with the progression of the discharge towards ending, the temperature distribution became more spatially uniform over the whole surface of the cell (Figure 4c,d), with the most heated region located in the center region of the cell.This observation of the shifting of hottest point is in congruence with the observation of Li et al. [23], who investigated cell temperature by inserting sensors inside the cell.This observation can be attributed to the faster depletion of the active species near the positive area due to the higher current density.Similar patterns were observed with charge and discharge at other current rates.It is important to mention here that individual IR images have different temperature scales.Here, the principal objective was to observe the spatial non-uniform temperature distribution of a cell at a certain state.By incorporating a different temperature scale, it was possible to observe and show the distinguishable non-uniformity of surface temperature distribution at a particular stage of operation.Figure 4d shows the temperature distribution of heat dissipation at rest (no applied current).This pattern suggests that the heat dissipation rate was higher at the upper half region of the cell (the half that contains the tabs), apparently because of higher heat transfer through the tabs.
The LFP cell consisted of a similar graphite based anode, therefore, the thermistor (placed at adjacent area of the positive tab) based surface temperature profile shows similar trend as the trend of the NMC cell (see Figure 5), particularly during discharge at all current rates.However, the temperature profile during charge showed a slight variation.At the beginning of the charge process, a small drop in temperature was observed, followed by a sharp rise until the SoC reached a level of ca.50%-60%.After that, again, a small drop appeared, followed by a rise until the end of charge.
In order to explain this variation from the NMC cell, the following factors can be considered.(i) It was found that the entropy of an LFP cathode does not change significantly during the charge/discharge process [21], therefore, a major portion of the entropy related heat is contributed by the graphite anode; (ii) an LFP cell operates within a lower voltage window compared to NMC; (iii) LFP shows a two phase system of LiFePO4 and FePO4 during Li + intercalation/deintercalation [24].Additionally, it was also shown that the distribution of these two phases, which vary across the electrode, are at a certain SoC level [25].Additionally, the geometrical attributes of the cell also had to be considered.The LFP cell is slightly smaller than the NMC cell with comparatively wider tabs (see Figure 6a).Although, despite of this geometrical variation, spatial temperature distribution showed similar evolution of the contour, as shown in Figure 6a-c.The LTO anode based cells have a distinct difference compared to the other two cells (i.e., the NMC and LFP cells).Titanate based anode has a higher electrical conductivity and a bigger surface area, thus, it shows lower impedance compared to graphite anode (see Figure 2 and Table 1).One of the reasons behind this variation is the distinction of solid-electrolyte interface (SEI) layer formation and the thickness characteristics between graphite and tianate anodes [3,26].Thus, the net heat generation in the LTO cells is much lower, which makes the LTO cell a potential candidate for a cell with higher thermal safety, even during very fast charging (within less than 30 min).It was also found that the entropy change is insignificant in an LTO based commercial full cell [22].These properties were reflected in the surface temperature profile of the LTO cell (see Figure 7).Discharge at different rates showed a sharp rise without any noticeable drop in the surface temperature (Figure 7a).However, charge at a low current (10 A) showed a change of slope in the rising temperature profile, approximately at the middle of the charge process, which further diminished gradually during charging at high current rates (20-100 A) (Figure 7b).Additionally, the LTO battery cell had different geometrical attributes.For instance, much wider tabs were located at the opposite side of the cell (see Figure 8a).Despite of the difference in material and geometrical properties, the spatial temperature distribution of the cell surface showed similar contours to the NMC and LFP cells (see Figure 8b,c).A summary of surface temperature rise of all the three cells is presented in Figure 9. Here, ∆T represent the temperature rise for each process (charge/discharge) based on the measurement by the thermistor.In addition, it was calculated as follows: where, = Temperature at the end of the process (charge/discharge), = Temperature at the start of the process (charge/discharge).

Micro Pulse Cycling
Thermistor (placed at the adjacent area of the positive tab) based surface temperature profile of the three cells showed a similar trend during high current (80 A) micro pulse cycling (see Figure 10).The inset of Figure 10 shows one complete cycle of ~1000 micro pulses of charge and discharge cycles.At the beginning, the temperature followed a sharp rise.Approximately after ~500 cycles or ~1 h, the temperature reached a steady state condition (i.e., variation <1 °C per 5 min).Micro pulse cycling was performed on all the cells, which were at an 80% SoC level.In previous work, we showed that micro pulse cycling performed at different SoC levels showed similar surface temperature profiles.Additionally, the temperature rise did not change significantly at different SoC levels [15].According to the IR images (Figure 11a-c) of NMC during micro pulse cycling, it was observed that, as the cell surface temperature proceeded towards a steady state, the spatial temperature distribution over the surface became non-uniform.In addition, the upper half portion of the cell, which is near the tabs, was comparatively hotter than the other half.It is also visible from the IR images that the hottest point is mostly located near the positive tab of the cell.Local high current density at the adjacent areas of the tab, for a very short time, can be attributed to this localization of the hottest region.However, evolution of the surface temperature contour in the case of the LFP cell showed a different pattern compared to the contour of the NMC cell (Figure 11e-f).In this case, the location of the hottest region was observed around the center region of the cell throughout the cycling.This evolution pattern is more similar to the pattern during continuous charge/discharge of the LFP cell (Figure 6).Relative spatial distribution of two phases (LiFePO4/FePO4) over the electrode, which changes during the charge/discharge process, might be related to this surface temperature distribution evolution [25].In the case of the LTO cell, at the beginning of the cycling, the hottest region appeared at the upper half of the cell (the half that contains the positive tab) (Figure 11g,h).However, at a steady state, the hottest region shifted towards the lower half of the cell (Figure 11i).It is important to mention here, that, at the steady state, the highest temperature difference over the complete surface was less than 2 °C.Variation in geometrical attributes (size and position of the tabs) and material properties (LTO based anode) can be considered in order to explain this variation compared to other to cells.Further, individual IR images contains different temperature scales in order to show the distinguishable non-uniformity of the cell surface temperature at a certain stage of operation.

Modeling Perspective
A large number of works have been published on the modeling of Li-ion batteries [8,[27][28][29][30][31][32][33].These models were created to predict electrical, electrochemical, thermal and mechanical behavior of the battery system, either independently or in combination (coupled).According to the literature, the scale of these models may vary from the microscopic level (molecular level quantum behavior based) to the mesoscopic level (electrochemical process and transport phenomenon based) or to the macroscopic level (based on the distribution of current and voltage on electrodes) to blackbox system level (empirical, data fit based).A comprehensive review in this regard can be found in Reference [8].
Thermal models generally describe the thermal behavior in terms of heating of the battery and are built on the energy conservation of the system as follows: where, T-temperature of the system domain, ρ-density of the system component, Cp-specific heat capacity, K-thermal conductivity and Q-volumetric heat source.Depending on the level of the modeling, Equation ( 2) can be applied on the separate domain of the battery.For instance, for macroscopic level the battery can be divided in three domains: electrode, tab and casing.Thus, Equation (2) for the three domains will be: At the boundaries of the domains heat transfer to the surroundings by radiation and convection, where, hconv-convection coefficient, hrad-radiation coefficient and Tamb-ambient temperature.Now, the variation between the existing modeling principal lies in the determination of the volumetric heat Q.Q can be formulated as: where, QEch and QJ are electrochemical heat and joule heating, respectively.For modeling at the system level, considering the system as a blackbox, the right side of Equation ( 6) can be formulated as: for the electrode domain ( 7) , for the tab domain (8) where, εe, εt are the volume ratio factors of electrodes and tabs respectively, R is the internal resistance of the cell, is the resistance of the tab, Vocv and I are the charge/discharge currents.
To model the non-uniformity of the heat across the area of the cell a thermal equivalent circuit model can be assumed in the different direction of the cell.The circuit generally consists of several thermal resistance and thermal capacitance network.Through electrical characterization, relative amount of heat generation in different direction can be can be modeled.For instance, in Reference [9], a first order Cauer model was considered and found in fair agreement between the modeled and real thermal behaviors.Due to the simplicity and low number of parameters, this modeling is efficient in terms of computational power.As shown previously, the LTO cell has a lesser degree of temperature non-uniformity compared to the NMC cell and the LFP cell, a thermal equivalent circuit model can be very effective in this case.However, as the model is built based on characterization data at a sets of test conditions, the model accuracy can be low in dynamic conditions, especially in the case of a cell with severe temperature non-uniformity, as the NMC and LFP cells presented above.
In order to overcome the accuracy limitations, the model can be extended to macroscopic level, taking into account the potential gradient over the surface of the electrodes [29,30,33].In this model formulation, charge balance can be applied to a pair of electrode as: where, σj is the electrical conductivity of the current collector σ, Vj is the potential distribution on the current collector, J is the transvers current density perpendicular to the electrode surface, n is the unit normal vector, δj is the thickness of the current collector and the subscript j corresponds to the positive and negative electrode.In this case, Equations ( 6) and ( 7) become as follows: where, δpe, δs and δne are the thickness of positive electrode, separator and negative electrode respectively.The transverse current density can be expressed in terms of open circuit voltage and voltage of the electrodes [29] as: Here, Vocv and R can be determined experimentally.With this model, the thermal behavior of any pouch type cell can be modeled with very high accuracy during constant charge or discharge.However, against dynamic current profile such as micro cycle as presented previously, this model may show less accuracy.In order to overcome this limitation, the non-faradic component of the transverse current can be taken into account.As suggested by Taheri et al. [34], the transvers current density can be expressed as: where, in the faradic component, Y and U are fitting parameters and expressed in terms of depth of discharge (DoD).Additionally, in the non-faradic component, Cdl is the lumped double layer capacitance.
As presented above, a macroscopic level modeling, assumption was made that electrochemical reactions are homogenous over the electrode area.At the mesoscopic level, modeling the inhomogeneity of electrochemical reaction parameters, such as charge balance in the liquid and solid phases, can be incorporated to improve the accuracy of the model [30].For instance, in order to model the shifting of the hottest point at the end of a continuous discharge process, a mesoscopic electrochemical-thermal model can be incorporated for higher accuracy.However, this kind of modeling is computationally challenging.Moreover, the parameters needed for the model are very difficult to determine.

Conclusions
Battery cell surface temperature distribution under high current continuous charge and discharge up to 100 A, as well as under high current micro-pulse cycling at 80 A, were studied for three different commercial pouch type cells: NMC cathode based 20 Ah cell, LFP cathode based 14 Ah cell, and LTO anode based 5 Ah cells, by using both contact thermistor measurements and infrared thermography.Variations in temperature profiles between the different types were observed during continuous charge/discharge.However, the spatial distribution patterns of the surface temperature were observed to be similar in the case of these three types of cells.At the beginning of the continuous charge/discharge test, the hottest region appeared near the positive tab but as the process progressed towards the end, the temperature distribution became more uniform with the hottest region located at the center region of the cell's surface.Comprised of a titanate based anode, the LTO cell showed a very low temperature rise, which proved it as the best option from the thermal safety point of view.On the other hand, during micro pulse cycling, a similar trend of the thermistor based temperature profiles was observed.According to this trend, the temperature rose sharply until a common point in time and became steady without any significant rise until the end of the cycling.The location of the hottest region was observed near the positive tab throughout the micro pulse cycling in the case of the NMC cell.The pattern of the surface temperature distribution throughout the micro pulse cycling in the case of the LFP cell was similar to the pattern observed during continuous charge and discharge, showing the hottest region in the center area of the cell.While in the case of the LTO cell, although the hottest region was located at the upper half (towards positive tab) at the beginning of cycling, it shifted towards the area adjacent to the negative tab of the cell at the steady state.
Non-uniformity of the surface temperature cannot be attributed to a single independent factor.It is a result of a combined influence of electrochemical properties, such as the inhomogeneity of the reactions across the electrode, of electrical properties, such as localization of current density and SoC, and of physical properties, such as the cell size and the geometry, and locations of the tabs.In most cases those factors are interdependent.This makes it very cumbersome to point out the most influencing factor.In order to model the non-uniformity of the cell temperature, several modeling methods have been adopted by researchers.These methods have been discussed in this work.Moreover, the effectiveness and limitations of these methods were concisely discussed from the point of view of the thermal behavior of the investigated cells.
Nevertheless, as temperature is one of the major factors in battery cell ageing, non-uniformity of the temperature distribution may lead to imbalanced ageing, and, ultimately, may lead to the degradation and imbalance of cell performance.In order to avoid such undesirable behaviors, improvements needs to be made at both the cell level and pack level.At the cell level, material properties and physical properties can be improved, whereas at the pack level, improvements in cooling architecture can be incorporated.

Figure 1 .
Figure 1.(a) Two locations on the cell surface subjected to measurement comparison; (b) Temperature profile measured by the IR camera and the thermistor at the corresponding locations of the cell surface during a discharge at 100 A of an NMC, 20 Ah cell.

Figure 3 .
Figure 3. Thermistor (placed near the positive tab of the cell) based temperature profile (red line) of NMC cell and corresponding load profile (blue line) during continuous charge at 20 A and discharge at 10 A, 20 A, 40 A, 60 A, 80 A, and 100 A. Discharge is denoted by negative current.

Figure 4 .
Figure 4. (a) Schematic of the NMC 20 Ah cell.IR images during discharge at 100 A along with the location of respective maximum temperature; (b) after ~5 s of discharge or ~1% of total discharge time (DT); (c) after 5 min or ~40% DT; (d) after 10 min or ~90% DT; (e) during rest (zero current).Different temperature scales were used here to show distinguishable non-uniformity of surface temperature distribution.

Figure 5 .
Figure 5. Thermistor (placed near the positive tab of the cell) based temperature profile (red line) of an LFP cell and the corresponding load profile (blue line) during a continuous charge at 10 A and discharge at 10 A, 20 A, 40 A, 60 A, 80 A and 100 A. Discharge is denoted by negative current.

Figure 6 .
Figure 6.(a) Schematic of the LFP 14 Ah cell.IR images during discharge at 100 A along with the location of respective maximum temperature; (b) after ~4 s of discharge or ~1% of discharge time (DT); (c) after 7 min or ~90% DT.

Figure 7 .
Figure 7. Thermistor (placed near the positive tab of the cell) based temperature profile (red line) of LTO cell and corresponding load profile (blue line) during (a) continuous charge at 10 A and discharge at 10 A, 20 A, 40 A, 60 A, 80 A and 100 A, respectively; (b) continuous charge at 10 A, 20 A, 40 A, 60 A, 80 A and 100 A and discharge at 10 A, respectively.Discharge is denoted by negative current.

Figure 8 .
Figure 8.(a) Schematic of the LTO 5 Ah cell.IR images during discharge at 100 A along with the location of respective maximum temperature; (b) after ~3 s of discharge or ~1% of discharge time (DT); (c) after 3 min or ~90% DT.

Figure 9 .
Figure 9. Temperature rise evolution in the NMC cell (during discharge), the LFP cell (during discharge) and the LTO cell (during both charge and discharge) at different current rates.

Figure 10 .
Figure 10.Thermistor based surface temperature profile of NMC, LFP, and LTO cells at 80% SoC during micro pulse cycling (~1000 cycles) at 80 A. Inset shows one single micro pulse cycle.

Figure 11 .
Figure 11.IR images during micro pulse cycling at 80 A of NMC, LFP and LTO cells at an 80% SoC at different points of time along with the location of respective maximum temperature.(a) NMC after 2 min; (b) NMC after 7 min; (c) NMC after 50 min; (d) LFP after 2 min; (e) LFP after 7 min; (f) LFP after 50 min; (g) LTO after 2 min; (h) LTO after 7 min; (i) LTO after 50 min.Different temperature scales were used here to show the distinguishable non-uniformity of surface temperature distribution.

Table 1 .
Mechanical and electrical properties of the NMC, LFP and LTO Cells.