Multiscale coupling of surface temperature with solid diffusion in large lithium-ion pouch cells

by   Jie Lin, et al.

Untangling the relationship between reactions, mass transfer, and temperature within lithium-ion batteries enables control approaches that mitigate thermal hot spots and slow degradation. Here, we develop an efficient physics-based pouch-cell model to simulate lock-in thermography experiments, which synchronously record the applied current, cell voltage, and surface-temperature distribution. Prior modelling efforts have been confounded by experimental temperature profiles whose characteristics suggest anisotropic heat conduction. Accounting for a multiscale coupling between heat flow and solid-state diffusion rationalizes this surface-temperature nonuniformity. We extend an earlier streamlined model based on the popular Doyle–Fuller–Newman theory, augmented by a local heat balance. The reduced-order model is exploited to parametrize and simulate commercial 20 Ah lithium iron phosphate (LFP) cells at currents up to 80 A. This work highlights how microscopic intercalation processes produce distinctive macroscopic heat signatures in large-format cells, as well as how heat signatures can be exploited to fingerprint material properties.


page 3

page 7

page 11

page 12

page 24

page 25

page 26

page 29


Design of pin-fin heat sink for the platform inertial navigation system by surrogate assisted techniques

In this study, in order to reduce the local high temperature of the plat...

Multiscale seamless-domain method for nonperiodic fields: linear heat conduction analysis

A multiscale numerical solver called the seamless-domain method (SDM) is...

Numerical Simulations of Frost Growth Using Mixture Model on Surfaces with Different Wettability

Frost growth on cold surfaces is a transient process with coupled heat a...

An Object-Oriented Library for Heat Transfer Modelling and Simulation in Open Cell Foams

Metallic open cell foams have multiple applications in industry, e. g. a...

Simultaneous recovery of surface heat flux and thickness of a solid structure by ultrasonic measurements

This paper is concerned with a practical inverse problem of simultaneous...

Simulating Ising and Potts models at critical and cold temperatures using auxiliary Gaussian variables

Ising and Potts models are an important class of discrete probability di...

1 Introduction

Specific energy, cycle life, safety, and cost of lithium-ion batteries have all substantially improved in the past decade Gibb2021:rise; Chu2016;Path; Nitta:2015Review; Ciez:2020cost, but challenges remain for high-power applications. The heat generation that accompanies charge or discharge Bernardi:1985Ge generally increases both a cell’s mean temperature and the extent of its temperature non-uniformity, which can significantly degrade performance and cycle life Ma:2018Review; Finegan:2015In; Zhu:2019Fast; Wang:2020Under. In worst-case scenarios, non-uniform heating of lithium-ion batteries at high power can facilitate catastrophic thermal runaway Feng:2018Review; Liu:2018Th. Outside of the effect on degradation, local variations of temperature within a battery cell can significantly impact characteristics such as equilibrium voltage and apparent internal resistance Yang:2019Asy; Yang:2021The. Temperature measurements contain rich information about the physicochemical processes that govern battery behaviour. If the microscopic origins of a battery’s thermal footprint are understood in sufficient detail, it may be possible to diagnose complex microscopic information from macroscopic temperature measurements. Temperature distributions can be monitored transiently via non-invasive thermal imaging techniques Robinson:2016Review; Finegan:2017Ch; Howie:2020pa; Wang:2020Th and implanted sensors Huang:2020Op.

Battery models are essential tools for exploring how different physical mechanisms contribute to measured behaviour Li:2020Mul; Mistry:2018Ele

. Experimentally validated models are further useful for estimating material properties and optimising cell designs

Deng2018:Safe. Physics-based electrochemical simulations based on porous-electrode theory, such as the Doyle–Fuller–Newman model Doyle:1993Mo; Fuller:1994Re; Fuller:1994Si, are well established. Almost all of the numerous porous-electrode-theory investigations in the literature focus only on the electrochemical response of a single layer in the ‘through-plane’ direction perpendicular to the current collectors Doyle:1996Ca; Cai2011:Ma; Yang:2017Mo; Yao2019:Quan; Finegan2020:spa; Lu2020:3D. The ‘in-plane’ distribution of current can be equally or perhaps even more important, however, especially in large-format cells.

This paper shows that the surface-temperature distribution across a large-format cell is an effective probe for battery diagnostics. We investigate the causes of non-uniform in-plane temperature distributions during battery charging and discharging, and explore various electrochemical processes that may rationalize them. Through use of a judiciously designed test rig that minimizes heating due to external wiring and tab contacts Howie:2020pa, lock-in thermography Robinson:2015De of commercial large-format 20 Ah LFP/Graphite pouch cells from A123 Systems is performed while the cells undergo square-wave cycling or constant-current discharge. The results of these experiments are simulated using a new computationally efficient three-dimensional battery model.

We formally derive the streamlined model previously proposed by Chu et al. Howie:2020pa from a three-dimensional version of the Doyle–Fuller–Newman model that is extended with a local heat balance. This process reveals several natural routes to produce reduced-order models that account rigorously for additional phenomena, while retaining the high computational efficiency and parsimonious parameter set of the streamlined model.

Excellent agreement between simulations and experiments is obtained by extending the streamlined model to include solid-phase diffusion dynamics. A single lumped diffusion time for both electrodes suffices to match experimental temperature profiles. The three-dimensional aspect of the electrode model differs substantially from typical approaches based on porous-electrode theory, and provides new insights into the dominant physical mechanisms that result in non-uniform temperature distributions. As well as reducing the large set of unknown material properties involved in Doyle–Fuller–Newman theory to a set of just a few observable parameters, the order reduction makes our model computationally efficient enough to support inverse-modelling algorithms based on iterative solutions of full finite-element simulations. We demonstrate, by extracting parameter values from measured cell data, that the streamlined model with solid diffusion can accurately estimate solid-phase diffusivity, as well as key material properties such as electrolyte conductivity, interfacial exchange-current density, specific heat capacity, and cell-reaction entropy, among others.

2 Temperature non-uniformity and solid-state diffusion

Lock-in thermography experiments were conducted to measure the surface temperatures of pouch cells synchronously with their voltage output under given applied currents. Experimental data was gathered using the test rig depicted in Fig. 1, for which the experimental setup and procedures were established by Chu et al. Howie:2020pa.

Figure 1: An infrared image of the instantaneous surface-temperature distribution across a 20 Ah pouch cell during lock-in thermography. Inset: Photograph of the experimental test rig.

Here we report data from two sources: square-wave-excitation cycling experiments, of which the data sets that cycled around a 30% state of charge (SOC) were reported earlier Howie:2020pa, but the others were not; and full-cell discharge experiments, performed specially for this report. All experiments used 20 Ah pouch cells from A123 Systems, which have a lithium iron phosphate (LFP) positive electrode and a graphite negative electrode.

In all cases, lock-in thermography was performed using cells initially equilibrated at ambient temperature. Before each cycling experiment, the cell was discharged from 100% SOC using Coulomb counting to a predetermined initial SOC of 30%, 50%, or 70%. Cells were cycled galvanostatically, alternating between charge and discharge periods of equal length for the 2500 s duration of the experiment, although the first charge step was performed over a half-period to keep the cell’s time-averaged SOC centred at its initial value. The applied current was set at 2C or 4C, with periods of 50 s or 100 s for one charge/discharge cycle. The cell voltage and ambient temperature were recorded at 1.0 Hz.

Thermograms of the cell surface were captured via the thermal imaging camera, allowing online visual monitoring of cell temperature. A physics-based battery model was solved using COMSOL Multiphysics software to simulate the electrical and thermal responses to the square-wave current excitation (see Methods section for more details). For all the square-wave cycling tests, parameter estimation based on a previously introduced streamlined model Howie:2020pa could capture the hot-spot, cold-spot, and surface-average temperatures accurately, but failed to predict the correct horizontal temperature distribution on the surface — that is, the distribution across the largest surface of the pouch cells, in the direction perpendicular to the tabs (see Figure S2 of the Supplementary Material).

In hopes of improving fits of the horizontal temperature distribution, several reduced-order models were derived, each based on the Doyle–Fuller–Newman (DFN) model with an added local heat balance. These produced a variety of extensions to the streamlined model of Chu et al. Howie:2020pa, as described in Supplemental Note 2. We found that the assumption of linear kinetics in place of nonlinear Butler–Volmer kinetics did not affect observed surface-temperature distributions (see Figure S3). The inclusion of solid-state diffusion, however, did have a fairly large macroscopic effect. In particular, solid-state diffusivity was found to be the only parameter that had observable impact on the horizontal temperature distribution. This finding suggests that there is a close coupling between thermal transport and solid-state diffusion in the electrode particles, and further, that consideration of this coupling is necessary to account for the horizontal temperature variation in large-format pouch cells.

Figure 2 presents simulation results that demonstrate how the spatial concavity of the horizontal temperature distribution (along a horizontal axis through the hot spot) varies with the solid-state diffusion time constant (cf. Table 1).

Figure 2: Coupling of horizontal temperature concavity through the hot spot with solid-state diffusion time. Insets: simulated surface-temperature profiles after  s of square-wave cycling with 80 A applied-current amplitude and 100 s period with a 30% initial state of charge (4C–100s@30%), at various diffusion times .
Parameter Symbol This work Streamlined model
30% 50% 70% 30% 50% 70%
Local reaction current density ()
Reaction activation energy ()
Effective ionic conductivity ()
Temperature coefficient of ()
OCP gradient ()
Diffusion time ()
Entropy change ()
Heat transfer velocity ()
Effective thermal conductivity ()
Table 1: Parameters yielded by model fits based on the new reduced-order model (This work) and the model of Chu et al. Howie:2020pa (Streamlined model), extracted from experimental lock-in thermography data under 4C–100s square-wave applied currents, with differing initial states of charge.

Details of how this concavity was estimated are provided in Supplemental Note 1 and Figure S4. The horizontal concavity flattens as the diffusion time constant decreases, ranging from 168 to 35 . At   — an unrealistically fast value Dees:2020Ap; Chang:2021Im — the hot spot lies very close to the cell’s top edge, and the horizontal temperature variation is minimal across the vast majority of the cell surface. This behaviour qualitatively agrees with results presented by Chu et al. Howie:2020pa, whose model derives from the assumption that solid-state diffusivity is infinitely large, as explained in Supplemental Note 2. If the diffusion time constant increases (i.e., the diffusion coefficient decreases), then horizontal temperature non-uniformity increases. At diffusion time constants above , significant extra heat generation occurs, causing an even larger temperature difference between the central vertical axis of the cell surface and its left or right edges, with generally higher absolute temperatures everywhere.

Typical experimental results for a cell cycled from an initial 30% SOC under a square-wave current having 4C amplitude and period—abbreviated henceforth as 4C–100s@30%—are shown in Fig. 3 and Supplemental Video 1.

Figure 3: Experimental and simulation results for a 20 Ah LFP pouch cell under 4C–100s square-wave cycling around an initial 30% state of charge. (a) Cell voltage, temperatures and horizontal temperature concavity. (b) Experimental measurements of surface temperature after 100 s, 500 s, 1000 s, and 2500 s of 4C-100 s cycling. (c) Model training results based on the experimental data shown in panel (a) and the streamlined model by Chu et al. Howie:2020pa. (d) Model training results based on the experimental data from (a) and the model proposed in this work. The battery tabs (not shown) are on the top edge of the thermal images in (b)-(d).

Before , the maximum cell temperature occurs in the areas directly adjacent to the tabs, and a temperature gradient develops primarily in the vertical direction. Near  , a small ‘hot spot’ — a point maximum of temperature on the cell surface — appears near the top edge of the cell and begins to move downward along the vertical centre line. Between   and  , a domain of higher temperature around the hot spot gradually expands, eventually achieving a relatively stable shape and location when  . The average cell temperature also arrives at a state in which it fluctuates around a relatively constant elevated value after about  . The mean surface temperature in this ‘periodic steady state’ is determined by the balance of Joule heating and convection from the cell surface, and the temperature fluctuations arise from entropic effects associated with the cell reaction. In the periodic steady state, the hot spot generally swells and contracts when applied currents have opposite signs; the hot spot expands when the reaction entropy effect is exergonic, and contracts when it is endergonic.

The strong and relatively isolated correlation between solid-phase diffusion and horizontal temperature concavity (cf. Fig. 2) justifies the estimation of solid-state diffusivity with temperature measurements. Taking the concavity of the temperature distributions into account, model parameters were estimated by fitting 4C–100s@30%, 4C–100s@50%, and 4C–100s@70% data using maximum likelihood estimation, as detailed in the Methods section. The parameters resulting from this estimation process are listed in Table 1, and the fits for 4C–100s@30% are also plotted on Fig. 3(a). Corresponding to the experimental thermograms, Figs. 3(c) and (d) show simulation results for 4C–100s@30% using the prior streamlined model Howie:2020pa and the proposed model, respectively, both based on the best-fit properties provided in Table 1.

It is clear that the model incorporating solid-state diffusion captures more features of the surface-temperature distribution. This is because solid-state diffusion introduces an additional time constant to the system. Two key relaxation times are visible in Fig. 3(a) besides : a reaction relaxation time , and a thermal time constant . A dimensional analysis shows that


where is the gas constant, , the ambient temperature, and , Faraday’s constant; is the rated capacity of the cell per unit of superficial electrode area () and , the cell thickness (); the remaining parameters are defined in Table 1. The short-time relaxation of cell voltage is controlled by . This is an electrical time constant, which arises from the fact that the slope of the cell’s open circuit potential (OCP) with respect to its SOC acts like an admittance, and the exchange-current density acts like a conductance. Thus the interface relaxes somewhat like a parallel RC circuit, with the OCP slope providing the (inverse) capacitance, and interfacial charge-transfer kinetics providing the resistance. Notably, this time constant has not been considered in most asymptotic analyses of the Doyle–Fuller–Newman model Moura:2017Bat; Sulzer:2019Fast. The thermal relaxation time is much longer than the other time constants. It controls the relaxation of the voltage envelope and the average cell temperature. As mentioned earlier, the solid-state diffusion time constant controls how the horizontal temperature concavity relaxes.

Parameter values in Table 1 are similar to those found by fitting with the streamlined model Howie:2020pa, with three notable exceptions—ionic conductivity, reaction current density, and thermal conductivity. Since work by Chu et al. Howie:2020pa ignored the concentration polarization caused by solid diffusion, the resulting voltage drop could only be attributed to poor effective ionic conductivity in the electrolyte, which had to be underestimated by ca. 50% to fit the cell-voltage response. Reduced ionic conductivity increases the amount of Joule heating, causing an overestimation of exchange current density in order to lower interfacial resistance and match the temperature. The change in thermal conductivity owes in part to the inclusion of hot-spot position in the cost function used during parameter optimization, as discussed in the Methods section. Separate simulations showed that changing the electrode’s effective thermal conductivity tunes the vertical position of the hot spot on the cell surface. This change led to fitted thermal conductivities of the order of , placing results in good agreement with independent thermal-characterisation tests undertaken on similar electrode materials Zhang:2014Si; Werner:2017Th.

Among all the battery properties, one expects on a theoretical basis that the volumetric exchange-current density , OCP gradient , diffusion time , and entropy change may vary with cell SOC, while the other material properties should be nearly independent of it. In the range of SOC studied here, however, the variations of fitted and with SOC were also minimal. The apparent constancy of in this range is qualitatively confirmed by the horizontal temperature-profile concavity data shown in Figure S5, which are nearly invariant with respect to cell SOC at a given current density. For A123 20 Ah LFP pouch cells, it appears that a single, SOC-independent diffusion coefficient and exchange-current density suffice to describe measured voltage and temperature behaviour up to 4C.

For validation purposes, simulations using the fitted parameters were also compared with experimental results at conditions not used for fitting, with different C-rates and cycling periods, specifically 4C–50s and 2C–100s cycling, at each SOC. Figure S5 provides the experimental and simulated voltage and temperature responses for two validation tests at 30% SOC (2C–100s@30% and 4C–50s@30%), and two 4C–100s parameter-estimation tests at 50% and 70% SOC. Measured and predicted thermal images under various square-wave cycling profiles at are plotted in Figure S6. The root-mean-square errors comparing simulations with experiments are for temperature and for voltage.

3 Cross-scale effects of non-uniform temperature

The model parameterized above can, with minimal modification, be expanded to full discharge simulations that retain most of parameters estimated from the square-wave-excitation cycling tests. As mentioned before, most parameters vary negligibly with SOC. Full discharges were simulated by leaving every parameter constant apart from the local OCP gradient and the entropy change, which were replaced by local functions of SOC gathered from either full-cell measurements (OCP) or manufacturer-supplied data (entropy change). A description of the model parameterization is available in Supplemental Note 4; a complete set of parameters is given in Table S2; details of how OCP and entropy were handled are discussed in the Methods section.

The maximum C-rate explored in this work is 4C, a fairly challenging test of the accuracy and generality of our approach. The LFP pouch cell was discharged at constant current from 100% SOC to 0% SOC, after first charging the cell to with a ‘CC-CV’ protocol (in which a 4C current was applied until the voltage reached , at which it was held until current decayed to C/100) and then resting for an additional hour. The cell voltage, current, and surface-temperature distribution were measured at . Figure 4(a) and (b) show the battery voltage and temperature response during the 4C discharge.

Figure 4: Model parameters extracted from square-wave-current perturbations of a 20 Ah LFP pouch cell extrapolate to simulate complete constant-current discharge accurately. (a) Cell voltage and maximum, surface-average, and minimum temperatures during full discharge at 80 A (a 4C rate). (b) Experimental thermograms showing the surface-temperature distribution at various times. (c) Simulation results. The battery tabs (not shown) are on the top edge of the thermal images in panels (b) and (c).

Initially ( ), similarly to the square wave cycling tests, the region close to the battery’s tabs rapidly warms up due to the higher current density there, forming a hot spot near the cell’s top edge. Between   and  , the size of the hot spot gradually expands vertically; the point of maximum temperature detaches from the top edge at  . As the discharge continues, the hot spot moves down towards the centre of the cell and grows in size. The centre of the hot spot shifts downward substantially, passing through ca.   at s and   at  .

The average surface temperature increases by during the 4C discharge. At the end of discharge, the surface-temperature distribution spans between the maximum temperature, at the hot spot, and the minimum temperature, at the bottom edge. Simulated transient surface-temperature fields are shown in Figure 4(c) and a real-time comparison between the thermography test and model simulation is provided in Supplemental Video 2; the cell temperature rise and distribution are well matched between the model and experiment throughout the discharge.

Figure 5(a) shows spatial and temporal variations of solid-phase, liquid-phase, and reaction current density through-plane (i.e. normal to the - plane shown in Fig. 4) through the in-plane (-) location of the hot spot.

Figure 5: Fine-grained simulation results for 4C constant-current discharge. (a) Spatiotemporal variations of solid-phase, liquid-phase, and reaction current along an axis through the hot spot at horizontal position and vertical position . Labels PE, Sep and NE indicate positive-electrode, separator, and negative-electrode domains, respectively. (b) Local reaction current density in central cross-sections of the LFP and graphite electrodes through the blue dash-dotted lines in panel (a). The contour lines are spaced at intervals of . (c) In-plane total, Joule, reaction, and entropic heat generation in the LFP and graphite electrodes at and . Each location in these surface plots shows the integral of local heat generation across the entire electrode thickness; the contour lines are spaced at intervals of .

Generally, the magnitude of liquid-phase current density increases toward the separator, while the solid-phase current drops. Reaction current is associated with the concavity of these profiles. Extrema of the reaction current are observed in both electrodes at all times. Electrochemical reactions in the electrodes are initially favoured at the electrode/separator interfaces, so the extrema originate near the separator in both electrodes. As discharge progresses, the extrema become peaks, and the reaction front moves over time toward the current collectors. The presence of peaks in reaction current owes to the solid-phase diffusion limitation, a phenomenon familiar from one-dimensional DFN models in the literature Fuller:1994Si. Diffusion also controls the rate at which the reaction-current peaks move towards the current collectors. Because the graphite electrode is thinner, the extremum of reaction current has reached the current collector at , whereas the LFP reaction-current distribution still exhibits a peak. But because the LFP electrode is thicker, the peak in reaction current is narrower, and has a higher magnitude; the separator side of the LFP is fully discharged just after , and can no longer sustain reaction current.

The current densities also vary substantially in-plane direction. Figure 5(b) shows in-plane reaction-current distributions along cross-sections cut through the instantaneous extrema of reaction current (these positions are indicated with blue dash-dotted lines in Figure 5(a)). Peak-to-valley variation relative to the average ranges from 2% to 20%, and differs between the two electrodes at a given instant, despite the solid-phase diffusion times being the same. When the extremum in reaction current resides at the separator or current collector, the reaction current has relatively uniform gradient top-to-bottom. Otherwise, the -location of the peak, cf. Fig. 5(a), is also associated with a peak in the - plane. Again, the peak is sharper in the thicker electrode.

Figure 5(c) presents the in-plane distributions of instantaneous heat generation, integrated across the thickness of the whole cell after and of 4C discharge. The total heat generated is broken down into contributions from Joule heat (arising from bulk resistances in Eqn. (7) of the Methods section), reaction heat (interfacial resistances), and entropic heat. The local maxima in total through-plane heat generation are similar to the hot-spot locations in Fig. 4. The total heat generation generally increases as the cell is discharged. Joule heating accounts for about half of the total: at

it is relatively uniformly distributed in-plane; at

its distribution is dominated by an increase from the top of the cell to the bottom. Reaction heat and entropic heat, on the other hand, both come to relatively sharp peaks in the - plane. They interfere destructively at , where entropic heat is negative, and constructively at . Nevertheless the sum of the two is always positive; reaction heat and entropy together dominate the placement of the peak in total heat generation. Although the reaction-current distributions in Fig. 5(b) are difficult to parse, it does appear that the average position of the extreme reaction currents correlates with the maximum in total heat generation. Thus it can be concluded that the position of the hot spot is largely controlled by nonuniform in-plane reaction current.

4 Conclusions

Multiple processes within a battery have distinctive heat signatures, including electrochemical reactions, interfacial kinetics, and Joule heating. Thus, in large lithium-ion pouch cells, surface temperature can be used as an effective probe to provide microscopic understanding. Using a combination of lock-in thermography and physics-based modelling, we characterised several material properties by inverse modelling of experimental tests with square-wave applied currents. Significant in-plane temperature non-uniformity was observed, and was attributed to a balance of the heterogeneous distributions of local charge state, interfacial Joule heating, and ohmic heat generation. Importantly, solid-phase diffusion of intercalated lithium was found to have a distinctive macroscopic heat signature, leading to concavity in the temperature distribution along an axis parallel to the current collectors but perpendicular to the battery’s tabs.

We showed that many inferred properties could be assumed constant over a very wide range of states of charge. A model that kept most parameters fixed but included a lookup table of pseudo-OCP data accurately simulated the cell voltage and thermal response during galvanostatic full discharges at 4C. Joule heating and reaction heat are generally comparable in magnitude; the position of the maximum temperature on the cell surface is controlled primarily by the microscopic reaction distribution and solid-phase diffusion. Since large-format lithium-ion cells are favoured for high-energy-density packs, detailed knowledge about non-uniform thermal states is critical to the understanding of battery performance and cycle life, particularly during fast charging or discharging.


Electrochemical testing

All tests reported were performed using commercial LFP pouch cells (AMP20M1HD-A, A123 Systems). The equilibrium open-circuit potential (OCP) of the cells was measured with a pseudo-OCP approach at C/25 applied current between (identified as 100% SOC) and (0% SOC). The charge capacity was determined via Coulomb counting under a CC-CV protocol at 1C until the current decayed to C/100. When a cell was set to a certain initial SOC, it was first charged to 100% SOC with a CC-CV protocol, then discharged to the required setpoint by Coulomb counting at 1C. All these tests were conducted using an automated battery test system (Series 4000, Maccor Inc.). In the lock-in thermography tests, the cells placed in the test rig were charged and discharged using a high-power bipolar power supply (BOP 10-100MG, KEPCO Inc.). Square-wave-excitation cycling data were gathered based on initial SOCs of 30%, 50% or 70%, with applied current at 2C or 4C and periods of or Howie:2020pa. Full-discharge experiments were conducted from 100% SOC to 0% SOC with applied current at 0.2C, 1C, 2C and 4C.

Pouch cell disassembly

To determine the correct physical dimensions and layer structure to be used in finite-element modelling of the LFP pouch cell, a cell was disassembled in a glove box after being fully discharged to . Photographs of single layers of the cathode and anode extracted from the pouch cell are shown in Figure S1. The geometric parameters of the battery and its components were measured and are summarized in Table S1.

Lock-in thermography and thermal image analysis

Lock-in thermography was conducted following the method of Chu et al. Howie:2020pa using a thermal imaging camera (A35sc, FLIR Systems). Spatiotemporal temperature data was gathered by image processing, including reference-temperature data averaged over a piece of black felt, labeled as the ‘ambient’ rectangle in Fig. 1. The ambient-spot temperature from thermography was calibrated against measurements from a thermocouple ( accuracy, Type T, Omega Engineering) placed behind a piece of felt at that location, also labelled in the Fig. 1 inset. The resulting baseline temperature was subtracted from each pixel of the image during data processing.

Battery model

Transient surface-temperature and voltage profiles were simulated using a reduced-order continuum model, which was derived from an extended Doyle–Fuller–Newman model that incorporates a local energy balance, as detailed in Supplemental Note 2. This 3D model was applied throughout the pouch interior. As described in section 2 above, a simplified model, which neglects diffusion in the electrolyte but retains solid-state diffusion in the electrodes, was deemed sufficient to fit the data.

Within the macroscopic volume spanned by a given electrode, this reduced-order model considers ohmic charge balances in the solid and liquid phases, respectively, such that


in which and are respectively the current density and electrical potential in phase (with subscript ‘l’ designating liquid and ‘s’, solid), is the pore surface area per unit electrode volume, and is the current density across the pore surface, defined such that anodic currents are positive. The effective electronic conductivity of the solid is , and the effective ionic conductivity of the liquid, . Ionic conductivity was taken to vary with absolute temperature according to


where is the conductivity at reference temperature and expresses its linear variation.

Temperature is distributed across the electrode domains according to a macroscopic thermal energy balance, derived under the assumption that the interpenetrating liquid and solid phases that make up the electrode sit at equal temperatures:


Here is the effective local volumetric heat capacity, is the effective thermal conductivity, and is the reaction entropy of the electrode half-reaction.

The applied current () and the electric ground () are defined at the positive and negative terminals (copper bars, cf. Fig. S9), respectively. The component of normal to interfaces between the current collector and the electrodes was taken to vanish; similarly, the components of normal to interfaces between the anode and cathode were taken to vanish. The outer edges of the pouch were taken to be electrically insulating. To bound the thermal portion of the problem, Newton’s law of cooling was adopted at the exterior surfaces of the cell:


in which is the heat transfer coefficient and is the ambient temperature.

Active particles within each electrode are taken to exist across a microscopic radial dimension at each point within the electrode. Within these spherical particles of radius , the concentration of intercalated lithium is taken to satisfy Fick’s law


in which is the solid-phase diffusivity. This microscopic mass balance is coupled to the macroscopic problem through a boundary condition


where stands for Faraday’s constant, and is required to be finite. Thus the local macrosopic interfacial current density completely specifies the diffusion dynamics within the solid particles.

As justified in Supplemental Note 2, the interfacial reaction currents were taken here to follow linear kinetics,


where is the exchange current density and , the gas constant. The temperature dependence of exchange-current density was modeled as guo2010single


in which is the exchange-current density at and is an Arrhenius parameter. Charge transfer is driven by the surface overpotential between the liquid and solid phases, which breaks down as


in which is the electrode’s equilibrium OCP relative to a reference electrode of a given kind.

During square-wave cycling tests, the perturbation in SOC is small, allowing it to be linearized following Chu et al. Howie:2020pa. Whereas the equation of Chu is based on the average SOC within the particle, the diffusion limitation requires that this be replaced by the effective SOC at the particle surface, , defined as


Here is the instantaneous surface concentration at within the solid and is the maximum lithium concentration the solid particles can accept. Taking account of the diffusion limitation’s effect on the surface concentration of intercalated lithium, one finds that the OCP satisfies


in which is the average fractional state of charge of the whole pouch cell. The constant parameter is included to describe possible OCP hysteresis during slow charge or discharge srinivasan2006existence. Here, a nonzero value of was included to model LFP cathodes; was taken to be zero for graphite anodes. Here, is the OCP at , represents the OCP gradient with respect to fractional SOC at , and is the reaction entropy at and , which appears because OCP satisfies the Maxwell relation


Generally we assume that the reaction entropy depends weakly on temperature, and is therefore a function of SOC only.

The same numerical model was used to simulate full discharges, except the linearized OCP curves from equation 15 were replaced with full-cell experimental pseudo-OCP discharge data in the nonlinear form


where the functions and come from experiments. (Note that no hysteresis term is present here because only discharges were modelled.) These data were used assuming the anode as a reference potential: thus the model used as the equilibrium voltage in LFP, and assumed the graphite OCP to be ground (0 V). In the model, reversible heating was computed from full-cell data by equally apportioning the reaction entropy between the two electrodes. This could be refined if reference-electrode measurements were available, but practically, the thinness of the cell normal to the electrodes means that an unequal distribution of reversible heat is difficult to discern. The full-cell OCP and data are provided in Figures S7 and S8, respectively. A detailed description of the model parameterization is available in Supplemental Note 4, and a complete set of model parameters is given in Table S2.

When solving the model, computational speed can be improved by a scaling analysis of the governing equations. By applying the scaling argument put forward by Chu et al. Howie:2020pa, the multi-layer internal geometry was homogenised, allowing the electrochemical model to be solved across a single representative layer. The scaling procedure and details of the simulated geometry are described in Supplemental Note 3.

Parameter estimation

Inverse modeling was based on iterative solution of the transient battery model. The parameter vector


was identified using a nonlinear least-squares fitting algorithm that minimised the error between measurements and simulations (respectively denoted with superscripts ‘exp’ and ‘sim’) at each time-step of voltage, ; maximum, minimum, and surface-averaged temperatures , , and , respectively; concavity of the temperature distribution ; and the position of the hot spot in the plane, . The objective function for the minimisation was expressed as a sum over all entries in the time series, as


where , , , and denote the ranges of , , , and in the experiment. The width and length of the pouch cell are and . The temperature concavity was calculated at the hot spot location by fitting the horizontal temperature profile through the hot spot (see Figure S3) with a quadratic polynomial, as described in Supplemental Note 1.


The authors gratefully acknowledge funding from the EPSRC Translational Energy Storage Diagnostics (TrEnDs) project (EP/R020973/1), the Faraday Institution Multiscale Modelling project (subaward FIRG003 under grant EP/P003532/1) and the STFC Futures Early Career Award. The authors also thank Andrew Wang for his help with pouch-cell disassembly.