Electrochemical Change Induced by Spherical Indentation in Lithium-Ion Batteries

: Since the Young’s modulus of the separator is weaker than that of the other materials inside a lithium-ion battery, local deformation may cause blockage or rupture of the separator, resulting in internal short-circuit or other disasters. This study collects the stress–strain relationship of various materials within the battery, and combines the mechanical model with the electrochemical model through the coupled relationship between the volumetric strain and the volume fractions of solid and liquid phases. From a two-dimensional electrochemical simulation of a spherical indentation on a layer-structured battery, it is found that there is local negative value of the side reaction overpotential on the negative electrode adjacent to the separator after the battery is deformed. A higher strain will cause a decrease in the negative overpotential, leading to a more serious deposition of lithium during the charge process. The deformation-dependent overpotential is evaluated and the lithium deposition is then quantiﬁed. Moreover, the issue of the separator thickness is explored. We ﬁnd that under an indentation, the thickness does not affect the charging voltage, while a thinner layer will reduce the separator porosity and thus lower the overpotential and increase the chance of lithium deposition.


Introduction
Lithium-ion batteries (LIBs) became one of the most reliant energy sources due to their advantages of high energy density and long cycle times. As LIBs quickly penetrated into our daily life, their safety became a major concern for consumers [1][2][3][4].
In post-mortem analysis of the spontaneous combustion of mobile phones, was found that the main cause of the incidents is the internal short circuit (ISC) [5], which can result from deformation of the batteries [6]. Two major factors were noted to account for the deformation: (i) the rupture of the separator under compression [7], and (ii) the penetration of lithium dendrites through the separator [8][9][10].
The compression of the separator further leads to deposition of lithium metal at the interface between the negative electrode and the separator [11] and batteries generally suffer from uneven stress caused by the early winding process [12]. When a LIB is charged at a high current or at a low temperature, in particular, a high amount of lithium metal will be accumulated on the negative electrode [13][14][15]. In addition to deposition, lithium metal could also gradually oxidize into lithium ions during the discharging process, called Li stripping [16].
A basic indicator of the lithium deposition is the overpotential between the electrode and the electrolyte [17], and it can be obtained in the electrochemical model by solving concentration diffusion, electric field, and reaction kinetics. The relevant reaction kinetics of the lithium metal deposition reported in previous studies are summarized in Table 1. α a,2 = 0.5, α c,2 = 0.5 Same charge transfer coefficients. Overestimation of Li metal oxidation. [11] α a,2 = 0.33, α c,2 = 0.67 Transfer coefficient is measured. Interfacial resistance is attributable to degradation. [18] Tafel Li metal will not re-oxidize to lithium ions. Li metal generated continuously. [9] Modified Butler-Volmer Reduced-order model with plating equations is presented.
Li metal deposition when the overpotential is less than zero. [19] Studies also show that ISC can be induced by the damage of the materials inside the batteries through external force, which gives rise to a new trend in the mechanics of the battery, focusing on the decrease in the electric resistance between the positive and negative electrodes in pressurized batteries [6], the relationship between the battery force and deformation in pouch batteries [20,21], the relationship between the mechanical properties of the separator and the strain rate [22], or the mechanical behavior of graphite anode via the discrete element method [21,23].
The progressive failure of LIBs under external force was also examined by the spherical indentation test [24,25]. Using this test, Zhou et al. analyzed the short-circuit behavior of lithium-ion pouch cells from experiments [26], and Li et al. perform a simulation of the coupling of four sub-models, i.e., the 2D mechanical model, the 2D short-circuit model, the 2D thermal model, and the 1D battery model, to predict the short-circuit process under mechanical indentation [27,28]. Since these battery models are constrained to one dimension, there remains a lack of research about the evolution of the local electrochemical quantities in the direction parallel to the battery layers.
In the present study, we establish a coupled model related to mechanics and electrochemistry to understand the electrochemical response of a deformable LIB under a spherical indentation before a short circuit occurs. Based on this model, we explore the possible lithium deposition caused by different degrees of deformation, and quantify the maximum deformation that the battery can withstand. The current paper is organized as follows. We first introduce the mechanical and electrochemical models, which are used to analyze the two-dimensional electro-chemical change induced by non-uniform deformation. Next, an electrochemical simulation is carried out to thoroughly discuss the time evolution of various field variables within the battery, such as volume fractions of solid and liquid phases, overpotential, current density, and state of charge when the battery is subjected to a spherical indentation. Finally, the conclusions of this study are presented.

Methodology
This study uses a Python package to carry out the analysis of the mechanical model, while the electrochemical behavior of the LIB is delineated by the Doyle-Fuller-Newman (DFN) electrochemical model [29][30][31], which is also called the pseudo 2-dimensional (P2D) model. Through the strain-porosity relation, the two models can be coupled so that the deformation-induced local lithium deposition within the battery can be explicated.

Mechanical Analysis
We assume that the battery is stacked with different materials along the z-direction in the order of copper (Cu) current collector, anode active material (Ano), separator (Sep), cathode (Cat) active material, and aluminum current collector (Al). Since it was reported that the Poisson ratio of a LIB is close to 0 [21], the battery can be greatly simplified by a one-dimensional mechanical model.
As shown in Figure 1, the materials in the spherical indentation problem are subjected to the stress of σ zz , and thus the stress field is only related to the r-direction. The total deformation in the z-direction is evaluated by the superposition of the deformation of each layer, which is written as where H and W represent the thickness of the battery before and after deformation. ∆h i is the deformation of each layer, and it can be expressed by in which h i and ε i indicate the thickness and the strain of each layer whose index i (from 1 to 5) refers to the layer associated with the positive current collector, positive active material, separator, negative active material, and negative current collector, respectively. As shown in Figure 1, the materials in the spherical indentation problem are subjected to the stress of σ zz , and thus the stress field is only related to the r-direction. The total deformation in the z-direction is evaluated by the superposition of the deformation of each layer, which is written as where H and W represent the thickness of the battery before and after deformation. Δh i is the deformation of each layer, and it can be expressed by in which h i and ε i indicate the thickness and the strain of each layer whose index i (from 1 to 5) refers to the layer associated with the positive current collector, positive active material, separator, negative active material, and negative current collector, respectively. Along with the observation that the Poisson ratio is close to zero when a LIB is subjected to compression, a small deformation of the battery (see Figure 1) allows us to assume that the strain ε i only takes the z-direction and it is the same as the volumetric strain ε V,i [20]. Under the quasi-static condition, the stress-strain relationship takes the form ( ) where σ zz,i is the stress tensor and g layer,i (σ zz,i ) is the strain function of the layer i. The strain function for the material associated with each layer, displayed in Figure 2, can be determined from the uniform compression experiment [21,23]. It can be seen from Figure  2a that when subjected to the same stress (up to 80 MPa), the copper and the aluminum are basically not deformed. However, the deformation of the negative electrode is the largest when the stress is small (below 15 MPa). With the increase in stress, the deformation of the separator rises sharply and it becomes the largest among the five materials. Figure 1. Schematic of the spherical indentation and a simplified one-dimensional mechanical LIB model. The thickness of each layer is taken from a former study [7].
Along with the observation that the Poisson ratio is close to zero when a LIB is subjected to compression, a small deformation of the battery (see Figure 1) allows us to assume that the strain ε i only takes the z-direction and it is the same as the volumetric strain ε V,i [20]. Under the quasi-static condition, the stress-strain relationship takes the form where σ zz,i is the stress tensor and g layer,i (σ zz,i ) is the strain function of the layer i. The strain function for the material associated with each layer, displayed in Figure 2, can be determined from the uniform compression experiment [21,23]. It can be seen from Figure 2a that when subjected to the same stress (up to 80 MPa), the copper and the aluminum are basically not deformed. However, the deformation of the negative electrode is the largest when the stress is small (below 15 MPa). With the increase in stress, the deformation of the separator rises sharply and it becomes the largest among the five materials. On account of the force balance, the stress of each layer must be the same, which means σ zz,i (r) = σ zz,j (r). Since each layer produces different strains under the same stress, the overall strain performance of the battery is the summation of the deformation of each layer under the same stress. Accordingly, with the help of the layer thickness information in Figure 1, we obtain the stress-strain relation of the whole battery, shown in Figure 2b.
To explore the deformation of the battery under non-uniform compression, we formulate the stress in Figure 2b in terms of the total strain as where the total strain can be expressed as Here, W(r) is specified by the geometric shape. Taking the sphere of a radius R pressed into the battery with the distance D as an example, we could have Substitution of Equations (5) and (6) into Equation (4) enables us to evaluate the stress field σ zz (r). Figure 2c displays the stress σ zz (r) in the case of R = 10 mm and D = 0.46 mm; the figure exhibits that the stress decreases to zero when r is around 2.9 mm. Having the stress distribution, Equation (3) helps to acquire the strain of each layer, shown in Figure 2d. It can be seen that, around the force center (r = 0), the separator layer is most strained, followed by the anode layer. However, within the range of 2.3 < r < 3.0, the simulation indicates that the anode owns a higher strain than the separator. On account of the force balance, the stress of each layer must be the same, which means σ zz,i (r) = σ zz,j (r). Since each layer produces different strains under the same stress, the overall strain performance of the battery is the summation of the deformation of each layer under the same stress. Accordingly, with the help of the layer thickness information in Regarding the deformation of each layer, the position y i,n for the i-type material of the nth layer (n indicates the number of a unit of layers) is given as where the bottom position of anode at the unit n = 1 is set to be zero. The deformation shape of the battery with 20 layered units is exemplified in Figure 2e, which shows that the uppermost area has a large curvature. Theoretically, this area shares a shear strain, which, however, is not considered in the current study for the sake of modeling simplicity.
To have a closer look, the first layered unit is enlarged and shown in Figure 2f. As can be seen, this unit contains two layers of each of the five materials, and it is observed that the deformation of the separator is larger than that of the other materials. This is in line with the results indicated in Figure 2d.
In the above discussion, we reviewed the mechanical model used to investigate the stress-strain relation of layer-structured batteries. In our study, the model will be coupled with an electrochemical method to examine the internal reactions of batteries when they are deformed. The procedure is as follows: As soon as the values of R and D are specified, we can have the total battery strain, shown in Equation (5), and then in turn the total stress can be obtained from Equation (4). Given that the stress of each layer is equal to the total stress, the constitutive relation in Equation (3) further enables us to determine the strain of each layer. Specifically, the coupling of the strain of each layer to the volume fraction will be performed.

Strain-Induced Change in Porosity
In order to characterize the effect of deformation on the battery, we convert the strain field into the volume fraction (or porosity) field, which is one of the main field variables in the electromechanical modeling. Since the inside of a LIB is a porous material (denoting anode, separator, and cathode here) and can be treated as a mixture of a solid phase and a liquid phase, when the mixture is compressed, the fluid would flow freely relative to the solid phase.
As shown in Figure 3, the orange area signifies the solid portion of the porous material, while Vs refers to the volume of the solid phase. The blue part symbolizes the fluid portion of the porous material, whose volume is represented by Vl. The control volume, featured by the red box, can be defined as the volume spanned by the solid skeleton. When the mixture is compressed, denoted by the three arrows, the fluid flows out of the control volume, and the fluid inside the control volume decreases. Inspection of Figure 3 reveals the geometrical relations before and after the compression: where e s and e l indicate the volume fractions of solid and fluid, and the quantities with prime above denote those quantities after the deformation. Given that there is a volume change ∆V in the control volume, we could have V' = V + ∆V = V(1 + ε V ). Furthermore, we assume that the volume change is caused by the fluid flowing out of the control volume, which implies V' s = V s and V' l = V l + ∆V. Hence, from Equation (8), the volume fractions of the two phases after deformation can be obtained as based on which the mechanical and electrochemical models are coupled.
based on which the mechanical and electrochemical models are coupled. Considering the battery structure in Figure 2e, the layer thickness in Figure 1, and the case of R = 10 mm and D = 0.46 mm, we obtain the distribution of the volume fractions for the solid active material in Figure 4a and those for the electrolyte in Figure 4b. As shown in Figure 2a, when the battery is compressed by a stress of more than 20 MPa, the separator will deform greatly, resulting in a rapid decrease in the volume fraction of electrolyte within the separator. This accounts for the large e s ′ and the small e l ′ in the vicinity of r = 0 in Figure 4a,b. An additional simulation is carried out by increasing the thickness of the separator from 13 to 23 μm, while the thickness of the other materials remains unchanged. We find in Figure 4c that this increase results in a rise of the volume fraction of electrolyte in the vicinity of r = 0, which indicates that a great reduction in the thickness of the separator would significantly affect the mechanical properties of the battery. Considering the battery structure in Figure 2e, the layer thickness in Figure 1, and the case of R = 10 mm and D = 0.46 mm, we obtain the distribution of the volume fractions for the solid active material in Figure 4a and those for the electrolyte in Figure 4b. As shown in Figure 2a, when the battery is compressed by a stress of more than 20 MPa, the separator will deform greatly, resulting in a rapid decrease in the volume fraction of electrolyte within the separator. This accounts for the large e s ' and the small e l ' in the vicinity of r = 0 in Figure 4a,b. An additional simulation is carried out by increasing the thickness of the separator from 13 to 23 µm, while the thickness of the other materials remains unchanged. We find in Figure 4c that this increase results in a rise of the volume fraction of electrolyte in the vicinity of r = 0, which indicates that a great reduction in the thickness of the separator would significantly affect the mechanical properties of the battery.

Electrochemical Simulation
The DFN model, developed by Newman's group, is the most widely used electrochemical model [29][30][31]. It usually treats the temperature as a constant parameter

Electrochemical Simulation
The DFN model, developed by Newman's group, is the most widely used electrochemical model [29][30][31]. It usually treats the temperature as a constant parameter and simplifies the jelly roll of the cell to a one-dimensional geometry that is composed of five segments, which correspond to the thickness of the two current collectors, of two electrodes, and of a separator. This model is dominated by several partial differential equations, characterizing the interfacial reaction kinetics at the electrode surfaces, the charge conservation in the solid and in the liquid phases, the solid phase diffusion in the electrodes, and the mass transfer process in the electrolyte phase. These equations are listed in Table 2. Table 2. Equations of the P2D model of a LIB.

Descriptions Equations
Interfacial reaction kinetics In the current study, the axis-symmetrical deformation of the layer-structural battery is explored. Thus, we extend the accustomed one-dimensional electrochemical analysis to a two-dimensional simulation, in which the zand r-dependent field variables are considered.
In view of Equation (9), the mechanical-electrochemical coupling can be introduced into the set of the governing equations through the volume fraction-dependent diffusion coefficient and conductivity [11], that is D l = D l,bulk e l γ , D s = D s,bulk e s γ , κ l = κ l,bulk e l γ , κ s = κ s,bulk e s γ , where γ stands for the Bruggeman exponent, whose value is given in Table 3. This exponent is used to take the tortuosity factors of porous media into account. The interfacial reaction kinetics in Table 2 is the Butler-Volmer equation, where the charge transfer coefficients α a,1 and α c,1 are both set to be 0.5 in this study. This equation expresses the reaction flux j 1 , which comes from the difference between reduction and oxidation reaction rates of lithium at the electrode-electrolyte interface. In addition to the flux flowing into or out of the two electrodes, certain side reactions, such as lithium plating, may occur in the negative electrode. When this reduction reaction of lithium ions takes place, they may return to lithium metal in the electrolyte, and the reduction flux j 2 can be written as where the lithium plating exchange current density i 0,2 is set to be a constant of 0.001 A/m 2 [9]. In this study, the two parameters α a,2 and α c,2 are given as 0.33 and 0.67, respectively [19], and α a,2 < α c,2 indicates that the reducing power of lithium ions is stronger than the oxidizing power of lithium. Since lithium ions are directly reduced to lithium metal, there is no interface impedance and the reference potential (Li vs. Li + /Li) is zero. Hence, the SR overpotential η 2 for this side reaction (SR) process can take the form: According to Equation (12), when η 2 < 0, there is a reduction tendency for the lithium ions so that the lithium plating on the surface of negative electrode particles is more liable to occur.

Battery Strain of 8%
A higher strain case is first studied to highlight the effect of deformation on the deposition amount of lithium metal. We consider a rigid sphere of radius R = 6 mm, and it presses the battery into an indentation depth of D = 0.368 mm, giving rise to 8% compression with respect to the thickness of the battery (or a strain value of 8% along r = 0). It is found in Figure 5a that the volume fraction of the electrolyte in the separator greatly decreases around r = 0, which indicates that the separator suffered from a large deformation. When r > 2 mm, the deformation effect significantly diminished, and thus all of the volume fractions of the electrolyte in the three materials approach their original values (see Table 3). Now the deformed battery is charged with a constant current of 1C and then discharged at 1C, it is expected that lithium will be deposited at the junction of the separator, and the negative electrode in the charge/discharge process, and the deposition amount, can be evaluated in terms of the molar number per unit volume as:  Figure 5b exhibits the lithium metal production field at the ending point of charging, i.e., t = 3600 s. We observe that more lithium deposition occurs at the center of the indentation and that there is a large gradient of the production field along the vertical direction. Figure 5c shows the electrolyte potential ϕ l at t = 3600 s, and the colors refer to the Now the deformed battery is charged with a constant current of 1C and then discharged at 1C, it is expected that lithium will be deposited at the junction of the separator, and the negative electrode in the charge/discharge process, and the deposition amount, can be evaluated in terms of the molar number per unit volume as:  Figure 5b exhibits the lithium metal production field at the ending point of charging, i.e., t = 3600 s. We observe that more lithium deposition occurs at the center of the indentation and that there is a large gradient of the production field along the vertical direction. Figure 5c shows the electrolyte potential φ l at t = 3600 s, and the colors refer to the potential magnitude of the electrolyte. It can be found that the highest potential appears at the top of the positive electrode (where it is in contact with the current collector) around the center of the indentation, while the lowest potential appears at the current collector of the negative electrode. On the other hand, the solid electrode potential field φ s at t = 3600 s is exhibited in Figure 5d. Since the solid electrode has good electrical conductivity, the potential of the positive electrode is kept at about 4.2 V, while the negative electrode is kept at a potential close to 0 V.
During the charging process, lithium ions will move from the positive to the negative electrode. The deformation due to the indentation results in a nonuniform distribution of the concentration of lithium ions, which is zand r-dependent. This concentration can be readily expressed in terms of the local state of charge (SOC) on the surface of solid electrode particles, defined as where C s,surf , C s,max , and C s,min are the local solid-phase surface, the maximum, and the minimum concentrations of lithium in the solid electrodes, respectively. The time evolution of the distribution of the SOC cs field is simulated and displayed in Figure 6. It is observed that before charging, the SOC cs in the positive electrode is close to 0.95 while nearly 0.1 in the negative electrode. In the middle stage of charging (t = 1800 s, see Figure 6b), the SOC cs rises quickly in the negative electrode near the separator under the compression. At the end of charging (t = 3600 s, see Figure 6c), there is a large gradient of the SOC cs around the compression area of the negative electrode, and this area, which is near the current collector and under the indentation, is difficult to charge.
Batteries 2022, 8, x FOR PEER REVIEW 10 of 1 evolution of the distribution of the SOC cs field is simulated and displayed in Figure 6. It i observed that before charging, the SOC cs in the positive electrode is close to 0.95 whil nearly 0.1 in the negative electrode. In the middle stage of charging (t = 1800 s, see Figure  6b), the SOC cs rises quickly in the negative electrode near the separator under th compression. At the end of charging (t = 3600 s, see Figure 6c), there is a large gradient o the SOC cs around the compression area of the negative electrode, and this area, which i near the current collector and under the indentation, is difficult to charge.

Charge with CCCV
Here, the battery is charged at a constant current (CC) for 3600 s to 4.2 V, and is then further charged at a constant voltage (CV) of 4.2 V. As shown in Figure 7a, the curren drops rapidly during the CV process. At t = 5000 s, the SOC cs distribution inside the battery is shown in Figure 7b, and we observe that around r = 0 this SOC cs distribution in the positive electrode is more uniform when compared with the distribution at the end of CC charging (t = 3600 s) in Figure 6c. In addition, it can be seen that the size of the domain with a relatively low value of SOC cs in the negative electrode is reduced, indicating tha lithium ions diffuse to the compression zone of the negative electrode. This diffusion wil continue and the SOC cs will gradually become uniform along the r-direction, which can

Charge with CCCV
Here, the battery is charged at a constant current (CC) for 3600 s to 4.2 V, and is then further charged at a constant voltage (CV) of 4.2 V. As shown in Figure 7a, the current drops rapidly during the CV process. At t = 5000 s, the SOC cs distribution inside the battery is shown in Figure 7b, and we observe that around r = 0 this SOC cs distribution in the positive electrode is more uniform when compared with the distribution at the end of CC charging (t = 3600 s) in Figure 6c. In addition, it can be seen that the size of the domain with a relatively low value of SOC cs in the negative electrode is reduced, indicating that lithium ions diffuse to the compression zone of the negative electrode. This diffusion will continue and the SOC cs will gradually become uniform along the r-direction, which can be seen in Figure 7c when the CV continues and the charging time extends to 15,000 s. Moreover, at the junction of the negative electrode and the separator, the time evolutions of the SR overpotential for the Li-metal deposition at the point of indentation (r = 0) and at the point deviating from the indentation (r = 4 mm) are demonstrated in Figure 7d. It is found that the SR overpotential at r = 0 is less than 0 when charging for about 1300 s. Even if this overpotential still keeps negative within 3600 s to about 6000 s in the CV process, the deposition of lithium ions is less likely to accumulate due to the small current in this process.  Figure 7d. It is found that the SR overpotential at r = 0 is less than 0 when charging for about 1300 s. Even if this overpotential still keeps negative within 3600 s to about 6000 s in the CV process, the deposition of lithium ions is less likely to accumulate due to the small current in this process.

Deformation and the SR Overpotential
Lithium metal deposition will not only result in battery aging, but will also increase the chance of internal short circuits. To determine whether lithium metal would deposit under various geometry of spherical indentation, i.e., change the radius of sphere R and

Deformation and the SR Overpotential
Lithium metal deposition will not only result in battery aging, but will also increase the chance of internal short circuits. To determine whether lithium metal would deposit under various geometry of spherical indentation, i.e., change the radius of sphere R and alter the strain of the battery, the SR overpontential of the negative electrode adjacent to the separator with r = 0 is evaluated. By treating the zero overpotential as the division, it can be seen from Figure 8a that when subjected to a small radius of the indentation sphere, the battery deformation under a strain of about 2.5% does not cause lithium metal deposition; however, when a larger radius is used, only about 2.2% of the strain can be tolerated. It can also be shown that changing the internal design of the battery would affect the ability of the battery to resist external forces. To explore the issue, we consider different choices of separator thickness and evaluate the terminal voltage when a sphere with a radius of 6 mm is pressed into a 6% relative displacement. From the electrochemical simulation in Figure 8b, it is found that changing the thickness of the separator from 8, 13, to 18 mm has little effect on the charging voltage. However, with the same compression, the mechanical analysis in Figure 8c indicates that thinner separators have greater strain, resulting in a lower volume fraction of the electrolyte (or porosity). Figure 8d illustrates that this drop further leads to a more negative value of the SR overpotential at the intersection of r = 0 and the junction of the negative electrode and the separator. Hence, it is concluded that a separator with a thinner thickness would increase the risk of the lithium deposition.

Conclusions
The mechanical and electrochemical response of a LIB subjected to a spherical indentation is numerically simulated. Two points can be summarized regarding the small compression on the electrochemistry of a LIB. First, it is found that the indentation causes local strain inside the battery, which induces the decrease of the volume fraction of the electrolyte in the separator and electrodes. This decrease makes the SR overpotential at the junction of the negative electrode and the separator become negative, resulting in lithium deposition during the CC charging process. Second, the map of the overpotential at the position under the indentation and at the junction of the anode and the separator is It can also be shown that changing the internal design of the battery would affect the ability of the battery to resist external forces. To explore the issue, we consider different choices of separator thickness and evaluate the terminal voltage when a sphere with a radius of 6 mm is pressed into a 6% relative displacement. From the electrochemical simulation in Figure 8b, it is found that changing the thickness of the separator from 8, 13, to 18 mm has little effect on the charging voltage. However, with the same compression, the mechanical analysis in Figure 8c indicates that thinner separators have greater strain, resulting in a lower volume fraction of the electrolyte (or porosity). Figure 8d illustrates that this drop further leads to a more negative value of the SR overpotential at the intersection of r = 0 and the junction of the negative electrode and the separator. Hence, it is concluded that a separator with a thinner thickness would increase the risk of the lithium deposition.

Conclusions
The mechanical and electrochemical response of a LIB subjected to a spherical indentation is numerically simulated. Two points can be summarized regarding the small compression on the electrochemistry of a LIB. First, it is found that the indentation causes local strain inside the battery, which induces the decrease of the volume fraction of the electrolyte in the separator and electrodes. This decrease makes the SR overpotential at the junction of the negative electrode and the separator become negative, resulting in lithium deposition during the CC charging process. Second, the map of the overpotential at the position under the indentation and at the junction of the anode and the separator is plotted, indicating that a strain higher than about 2% will result in a negative overpotential, and thus have a negative impact on a LIB.
The negative overpotential is an indicator of the lithium deposition. The current study shows that the overpotential at the junction of the negative electrode and the separator could become negative, which is compatible with the simulation result by Cannarella and Arnold [11], who show that the local electrode potential can decline to negative values along the separator/negative electrode interface in the charging process when there are defects (closed pores) in the separator. Since the reduction in porosity caused by the spherical indentation is analogous to the effect of closed pores in the separator, we believe that our analytical result is qualitatively correct.
Finally, it is noticed that the shear-induced volume change, which can be observed in porous media [32,33], is not explored in our mechanical modeling for the sake of simplicity. Whether the effect is significant in a battery system is an interesting topic to explore in future research.