Next Article in Journal
Reinforcement Learning-Based Delay-Aware Path Exploration of Parallelized Service Function Chains
Next Article in Special Issue
Numerical Investigation on a Diffuser-Augmented Horizontal Axis Tidal Stream Turbine with the Entropy Production Theory
Previous Article in Journal
A Note on Type-Symmetries in Finite Games
Previous Article in Special Issue
Transient Dynamic Response Analysis of Two-Dimensional Saturated Soil with Singular Boundary Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Diffusion-Induced Nonlinear Local Volume Change on the Structural Stability of NMC Cathode Materials of Lithium-Ion Batteries

1
Department of Mechanical, Robotics and Energy Engineering, Dongguk University, Seoul 04620, Republic of Korea
2
Gas Processing Center, College of Engineering, Qatar University, Doha P.O. Box 2713, Qatar
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(24), 4697; https://doi.org/10.3390/math10244697
Submission received: 16 November 2022 / Revised: 8 December 2022 / Accepted: 9 December 2022 / Published: 11 December 2022

Abstract

:
Electrochemical stress induced by the charging/discharging of electrode materials strongly affects the lifetime of lithium-ion batteries (LIBs) by regulating mechanical failures. Electrochemical stress is caused by a change in the local volume of the active materials associated with the lithium-ion concentration. The local volume change of certain active materials, such as nickel-rich LiNixMnyCozO2 (NMC), varies nonlinearly with the lithium content, which has not been considered in the stress calculations in previous studies. In this paper, the influence of nonlinear local volume change on the mechanical response of NMC-active materials is investigated numerically. The goal is achieved by using a concentration-dependent partial molar volume calculated from the previously obtained local volume change experimental results. A two-dimensional axisymmetric model was developed to perform finite element simulations by fully coupling lithium diffusion and stress generation at a single particle level. The numerical results demonstrate that (1) the global volume change of the particle evolves nonlinearly, (2) the stress response correlates with the rate of change of the active particle’s volume, and (3) stress–concentration coupling strongly affects the concentration levels inside the particle. We believe this is the first simulation study that highlights the effect of a concentration-dependent partial molar volume on diffusion-induced stresses in NMC materials. The proposed model provides insight into the design of next-generation NMC electrode materials to achieve better structural stability by reducing mechanical cracking issues.

1. Introduction

With the growing demand for energy storage devices, lithium-ion batteries (LIBs) are gaining more interest due to their higher capacity and longer cycle life [1]. LIB applications range from small medical and portable electronic devices to electric vehicles (EVs); however, their capacity decreases over time [2]. One of the main goals of the development of the next-generation batteries is to increase their efficiency by limiting the capacity-fading mechanisms [3]. Layered cathode materials, such as lithium nickel manganese cobalt oxide LiNixMnyCozO2 (NMC), are the most promising active materials with lower costs and higher energy density [4,5]. However, the battery capacity quickly decreases due to various mechanical failures caused by diffusion-induced stress (DIS) when lithium is introduced into (lithiation) or extracted from active materials (delithiation) [6,7].
A lot of effort has gone into the development of charging/discharging models to understand the mechanisms underlying DIS generation. One of the pioneering works are the studies of Christensen and Newman [8]. They investigated the lithiation-induced stresses and predicted mechanical failure in a single spherical graphite anode. A thermal analogy was then employed to calculate the DIS using finite element method (FEM) simulations [9]. Later, the DIS model was further updated by including the phase transition effects [10,11], grain boundaries [12,13,14], charging/discharging rates [15], material properties [16,17], active material morphology [18], surface stresses [19,20], yielding and plastic deformations of the active material [21,22,23,24], and solid–electrolyte interface [25,26,27]. Furthermore, several researchers have developed various stress-regulated charging/discharging strategies to reduce stresses using such DIS models [28,29,30]. Going a step further, the dynamic growth of failures in LIB model systems under the influence of DIS has also been simulated [31,32,33]. However, in order to reduce the simulation complexity, the mentioned studies have mainly used constant material properties, whereas most of the active material properties depend on the lithium concentration. Therefore, to accurately predict the charging/discharging behavior of the material, precise property values should be used when performing simulations.
In the last few years, many researchers have devoted themselves to determining the effect of concentration-dependent material properties on the charging and discharging processes of lithium-ion batteries. The researchers included the lithium concentration-dependent elastic properties especially Young’s modulus [34,35,36,37,38,39,40,41], yield stress [24,42,43], Poisson’s ratio [44], coefficient of chemical expansion or partial molar volume [45,46], toughness [47], and lithium diffusion coefficient [48,49] to DIS simulation models. For example, Deshpande et al. [50] used a simple cylindrical electrode particle to find that Young’s modulus variation has a significant effect on the evolution of DIS. Zhang et al. [51] examined the effect of the concentration-dependent elastic modulus on lithium diffusion and DIS generation using composition-gradient LCO cathode particles. They found that lithiation-induced stiffening regulates DIS in composition-gradient electrodes. Hong et al. [44] investigated the effects of concentration-dependent diffusivity, Young’s modulus, and Poisson’s ratio on stress evolution during lithiation of the Sn particle and concluded that the change in material properties with lithium content significantly alters mechanical failure modes. Cai and Guo [48] examined the effect of changing the diffusion coefficient and elastic modulus hardening with the lithium concentration on DIS in an anisotropic graphite anode particle. They found that diffusivity as a function of concentration increased concentration gradients and thus enhanced DIS because the change in the volume of the active material induced by the lithium concentration change is the main source of DIS. However, the effects of the concentration-dependent chemical expansion coefficient have rarely been investigated in the past. The study of nonlinear volume change and its effects on DIS is essential to determine strategies for avoiding mechanical failures in high-energy-density cathode materials such as NMC.
During charging, both the lattice parameters a and c of the NMC unit cell change strongly, which ultimately leads to a large change in the unit cell volume [52]. This volume change evolves nonlinearly with the lithium content. Distinct changes in the lattice structure are the main sources of stress generation that lead to various modes of mechanical cracks [53]. The disintegration and cracks induced by lithium diffusion in NMC materials disrupt the ionic and electronic conduction pathways [54]. Consequently, structural degradation accelerates the capacity decrease of NMC cathodes. However, the effects of the nonlinear volume change of NMC materials in stress diffusion problems have not yet received much attention.
In this paper, we focus on the effects of diffusion-induced nonlinear volume change on the evolution of lithium concentration distributions and the related stress development. To achieve the goal, we first calculated the partial molar volume (chemical expansion coefficient) based on experimentally obtained volume change values given in reference [52]. Next, we simulated the charging (delithiation) phenomenon in a single NMC particle using the concentration-dependent partial molar volume under galvanostatic charging conditions. We then compared the results of the constant and variable (concentration-dependent) partial molar volumes to evaluate the effects of the nonlinear volume change on the chemomechanical response of the NMC-active material. The main objectives of this paper are as follows:
1.
To construct a fully coupled finite element chemomechanical model in order to investigate the effects of a concentration-dependent partial molar volume on the mechanical response of NMC particle.
2.
To study and compare the simulation results for the constant and concentration-dependent partial molar volume of the NMC particle during the charging (delithiation) process.
3.
To investigate the simultaneous effects of the particle size or the charging rate and the concentration-dependent partial molar volume on the mechanical response of the NMC-active material.
The rest of the paper is organized as follows: we first developed a fully coupled chemomechanical model that considered the local effects of a nonlinear volume change on the distribution of stress and lithium concentration. We then validated the application of the concentration-dependent partial molar volume by comparing the finite element simulation results with the experimentally obtained volume changes. Afterwards, we compared the stress and lithium concentration results for the constant and variable (concentration-dependent) partial molar volumes for the delithiation (charging) process. We then evaluated the effects of the particle size and charge rate. Finally, the article is concluded in the last section.

2. Methodology

For numerical simulations, an isolated active particle with a radius R s was considered, as shown in Figure 1. The axisymmetric model system and its spatial discretization for a representative case of lithiation are illustrated in Figure 1a,b, respectively. We employed full coupling between lithiation/delithiation kinetics and mechanics of diffusion-induced deformations. The computational modeling of this fully coupled chemomechanical model is presented in the following subsections.

2.1. Fully Coupled Chemomechanical Model

2.1.1. Modeling of Lithium Ion Diffusion in an Active Particle

The diffusion of lithium ions inside the active particle was assumed to follow Fick’s law [55]:
c t + J = 0
where c is the lithium ion molar concentration; t is the lithiation/delithiation time, and J is the lithium ion flux inside the active material. The flux J is defined as a function of chemomechanical potential as [9]:
J = M   c   μ
where M is the mobility factor given by M = D / R T , and μ is the stress-dependent chemical potential, which can be defined by hydrostatic stress ( σ h ) as [9]:
μ = μ o + R T ln X Ω σ h
where μ o is the reference state potential; R is the universal gas constant; T is the absolute temperature; X is the molar fraction, and Ω is the partial molar volume. Combining Equations (2) and (3) leads to:
J = D ( c c R T ( Ω σ h ) )
Finally, putting Equation (4) to Equation (1) gives the final form of the partial differential equation for solving lithium concentrations in a coupled stress–concentration manner [56]:
c t D 2 c + D R T c ( Ω σ h ) + D c R T 2 ( Ω σ h ) = 0
The third and fourth terms on the left-hand side of Equation (5) result from the stress–concentration coupling. Thus, without stress–concentration coupling, the partial differential equation above reduces to:
c t D 2 c = 0

2.1.2. Modeling of Diffusion-Induced Stress

The DIS in the active particle was solved by the following partial differential equation for mechanical equilibrium [20]:
  σ + F b = 0
where σ is the Cauchy stress tensor, and F b is the body fore. In this work, no body force was assumed, so F b = 0 . For elastic deformations during the charging/discharging process, the stress–strain relationship is governed by Hook’s law as [33]:
σ = C : ε e
where C is the fourth-order stiffness matrix, and ε e is the elastic strain. The total strain caused by the elastic and diffusion-induced deformations is given as [33]:
ε t = ε e + ε d = 1 2 ( ( u ) T + u )
where ε t and ε d are the total and diffusion-induced strains, and u denotes the displacement field. Thus, the elastic strain is calculated by subtracting the diffusion-induced strain from the total strain as [9,14]:
ε e = ε t ε d
For the DIS calculations, the diffusion strain was calculated using the thermal analogy [8,9,57] as:
ε d = 1 3 c ˜ Ω I
where c ˜ is the concentration difference between the current and the initial state, and I is the identity matrix. In another method, using directly the volume change ( Δ V / V o ) obtained from the experiments, the diffusion-induced strain can be represented by the following diffusion-induced deformation gradient F d :
F d = ( 1 + Δ V V o ) 1 / 3 I

2.2. Material Properties and Numerical Simulations

To perform the finite element fully coupled chemomechanical simulations, the structural mechanics and transport of the dilute species modules of COMSOL Multiphysics (version 6.0) was employed in this work. The partial differential equation of the mechanical equilibrium was solved in the structural mechanics’ module, while the partial differential equations of the mass balance of lithium ions inside the active material were solved using the transport of the dilute species module. The built-in parallel sparse direct solver (PARDISO) transient solver was used to solve the weak forms of partial differential equations in a fully coupled manner. The delithiation phenomenon was modeled by applying a constant negative flux to the exposed surfaces of the active particle. The partial differential equations were solved in time increments until the local minimum state of charge reached the lower limit of the respective NMC-active material. To explore the effects of the concentration-dependent volume change of the active material, a representative case of 2 µm particle size of LiNi1/3Mn1/3Co1/3O2 (NMC-111) at 1C charge/discharge rate was considered, unless otherwise specified. Furthermore, the following initial and boundary conditions were considered to solve the partial differential equations given in Equations (5) and (7).
Initial conditions: At the start of the simulation, it was assumed that the maximum lithium concentrations ( c m a x ) were homogeneously distributed throughout the NMC active particles as:
  c = c t o t a l   at   t = 0
where c t o t a l is the stoichiometric lithium concentration in the NMC material. To solve the partial differential equation of mechanical equilibrium Equation (7), initially no stress was considered in the particle as:
( σ ) t = 0 = 0
Boundary conditions: To solve the issue of lithium diffusion inside the particle, a constant lithium flux was applied to the particle surface, such as:
J n = i n F = C rate α ρ F × R s 3
where n , i n , F , and C rate denote the outward normal unit vector on the external surface of the particle, the current density, the Faraday’s constant, and the charging rate, respectively. The α , ρ , and R s are the specific capacity, density, and radius of the active particle, respectively.
To quantify the effects of the nonlinear volume change, during the simulation the total volume change of the active material was calculated as:
Δ V V o = ( d v ) t = t V t = 0 V t = 0
where Δ V is the change in particle volume between the current (at t = t ) and initial (at t = 0 ) delithiation state. To explore the mechanical response of active materials, the maximum first principal stress ( σ m a x ) was calculated and compared for various cases. Further, to investigate the change in lithium concentration behavior, the difference between the maximum and minimum states of charge ( Δ SOC ) was evaluated as:
Δ SOC = c max c min c t o t a l × 100 %
where c max and c min are local maximum and minimum lithium concentrations. Finally, to assess the change in charge storage, the normalized capacity ( Π ) of the active particle was calculated as follows:
Π = F c d v α ρ V o
where Π will change between 1 and 0 during delithiation.
The convergence of both the maximum stress ( σ m a x ) and change in state of charge ( Δ SOC ) parameters was confirmed by increasing the degrees of freedom by refining the mesh size. According to the mesh independence test, the simulations were carried out with 5673 mesh elements and 46,356 degrees of freedom.

Concentration-Dependent Partial Molar Volume

To calculate the concentration-dependent partial molar volume, the volume change data during the delithiation process of NMC materials were obtained from reference [52]. The concentration-dependent partial molar volume of each NMC structure (i.e., NMC-111, LiNi0.5Mn0.2Co0.3O2 (NMC-523), LiNi0.6Mn0.2Co0.2O2 (NMC-622), and LiNi0.8Mn0.1Co0.1O2 (NMC-811)) was calculated separately. The volumetric strain ( λ ) due to the volume change caused by the delithiation was calculated as:
λ = ( Δ V V o + 1 ) 1 / 3 1 = ( V x = x V x = x max V x = x max + 1 ) 1 / 3 1
where Δ V is the volume change that is calculated as the delithiation states increase. The V x = x is current, and V x = x max is the initial (with maximum lithium content) unit cell volume. Based on this strain value, the partial molar volume is given by:
Ω = 3 × λ ( x x max ) × c t o t a l
Using Equation (20), the concentration-dependent partial molar volume of each NMC was calculated and plotted in Figure 2. For simplicity, we proposed that other material properties, such as Young’s modulus and diffusion coefficient of the active material, were independent of the lithium-ion concentrations. Other material properties and simulation parameters used in this work are listed in Table 1.

3. Results and Discussions

3.1. Validations of Numerical Results

In this section, we validate the calculations of the concentration-dependent partial molar volume of NMC with different nickel contents; NMC-111, NMC-523, NMC-622, and NMC-811. If the distributions of lithium concentration inside the active particle are sufficiently homogenized, then the local volume change will equal the percentage of the total volume change of the active particle. Therefore, in Figure 3, we compare the results by plotting the global volume change obtained by simulating a small particle under a lower charge rate (i.e., R s = 2   μ m and C r a t e = 1 C ) against the normalized charge capacity and the local volume change of the respective NMC material against the lithium content obtained from [52]. The close volume change values obtained from the simulations and the experiments suggest that this concentration-dependent partial molar volume can be used for stress calculations of NMC-active materials.
Then, the simulation results obtained by directly using the local volume change as the intercalation strain (Equation (12)) and using the calculated concentration-dependent partial molar volume as the thermal analogy (Equation (11)) are compared in Figure 4. To avoid changes caused by stress–concentration coupling, only uncoupled simulations were performed in this section using Equation (6). All simulation results are similar either by using the intercalation strains or thermal analogy. This, therefore, proves the validity of the thermal analogy method and the use of concentration-dependent partial molar volume instead of the direct use of the local volume change as the volumetric strain in finite element simulations.

3.2. Effects of Variable Partial Molar Volume

Figure 5 compares the simulation results obtained for a representative case of a 2   µ m NMC-111 particle at a discharge rate of 1C with and without considering the nonlinear local volume change during the delithiation process. Figure 5a shows that the total volume of the active particle ( Δ V / V o ) changes nonlinearly when the variable (concentration-dependent) partial molar volume ( Ω v a r ) is considered. On the other hand, using a constant partial volume ( Ω c o n s t ), the change of Δ V / V o is linear. This indicates that the local volume change significantly affects the global volume change, suggesting that the local volume change will also increase the stress due to the surrounding constraints provided by the binder, other active particles, current collector, and separator.
The evaluation of the stress generation caused by lithium diffusion is essential for the analysis of mechanical failure. For this purpose, the local maximum of the first principal stress ( σ m a x ) was analyzed here. DIS developed because of the local mismatch of lithium concentration levels between the inner and outer regions of the active particle. This mismatch caused the regions to expand/contract based on the magnitude of the partial molar volume. Consequently, using Ω v a r , the rate of volume change induced by lithium diffusion is different at different concentration levels, and the associated stress increase is different along the radial direction. As a result, σ m a x evolves differently for Ω v a r , as illustrated in Figure 5b. The evolution of σ m a x reveals several smaller peaks before a larger peak toward the end of the delithiation process, while these peaks are not visible when Ω c o n s t is used. Instead, σ m a x increases gradually after achieving a distinct peak in the early stages of the delithiation process. Figure 5b also shows that the peak magnitude of σ m a x is almost three times higher for the Ω v a r case compared to the Ω c o n s t case indicating a higher probability of mechanical failures. In summary, the evolution and levels of diffusion-induced stress in the active particle are significantly affected by the use of concentration-dependent partial molar volume. The reason for the appearance of such smaller peaks in the evolution of σ m a x is given in the following paragraphs.
The maximum concentration difference in the active particle measures concentration gradients to some extent. Therefore, the maximum SOC difference ( Δ SOC ) was evaluated to express the change in the evolution of the concentration gradients. Figure 5c compares the evolution of Δ SOC using Ω c o n s t and Ω v a r . For Ω c o n s t , Δ SOC initially increases sharply, and after reaching a distinct peak, the increase becomes more gradual. In the initial stages of a particle’s delithiation, only the outermost regions undergo the deintercalation process, so, over time, the SOC difference between the surface and the core increases rapidly, but when the deintercalation front reaches the particle’s core, the core starts to delithiate. The change in Δ SOC drops, which is the second stage of a gradual increase in Δ SOC . In contrast, when Ω v a r is included, the change of Δ SOC becomes irregular. It is affected by either a local volume change or a stress change through stress–concentration coupling.
To explore the effects of the concentration-dependent partial molar volume on the capacity stored inside the active particle, the normalized capacity ( Π ) was compared for both the Ω c o n s t and Ω v a r cases, as shown in Figure 5c. The same capacity value was observed for both partial volume cases. This indicates that although the nonlinear local volume affects the particle’s response to stress, the total stored charge remains the same.
Uncoupled numerical simulations were performed to isolate the effects of stress–concentration coupling and concentration-dependent partial molar volume on the SOC differences. In Figure 6a, the Δ SOC evolution is again plotted for Ω c o n s t and Ω v a r , without considering the coupling between the stress and concentration. The Δ SOC is the same for both partial molar volumes, although changes are still present in the evolution of σ m a x , as shown in Figure 6b when Ω v a r is considered. Thus, this proves that (1) changes in Δ SOC are caused by stress–concentration coupling, meaning that inhomogeneous stress levels due to nonlinear volume changes affect the lithium concentration levels, and (2) σ m a x is not directly affected by the change in Δ SOC .
Further, to figure out the mechanism underlying the stress response caused by the concentration-dependent partial molar volume, we investigated the rate of change of Δ V / V o   over time. The evolution of the 1st temporal derivative of Δ V / V o was plotted with the evolution of σ m a x in Figure 7. The trend of σ m a x perfectly matches the trend of the 1st derivative of Δ V / V o . This proves that the change in maximum stress behavior is directly affected by the rate of change in Δ V / V o .
We then investigated the effects of the concentration-dependent partial molar volume on the distribution patterns of the lithium concentration and the associated stress generation. It is clear from the comparison that the distribution patterns of the local lithium concentration and the first principal stress are consistent as the partial molar volume changes from constant to a variable. As shown in Figure 8a,b, the maximum local concentration of lithium occurs at the center of the particle, while the minimum occurs at the surface of the particle. Similarly, the local maximum first principal stress occurs at the surface of the particle, and the minimum stress occurs at the center of the particle, as shown in Figure 8c,d. This indicates that although the use of concentration-dependent partial volume enhances stress levels, it does not change the location of maximum stress.

3.3. Effects of Particle Size

Usually, the particle size ( R s ) inside the electrode varies. Particle size plays an essential role in the stress increase caused by lithium diffusion. Therefore, in this section, the particle size was changed from R s = 1   µ m to R s = 30   µ m in order to study the effects of concentration-dependent partial molar volume on the stress increase. Figure 9 compares the simulation results for various particle sizes with and without the use of a diffusion-induced nonlinear local volume change. The red dots to the right of each 3D plot represent the peak value during the delithiation time evolution of the corresponding parameter. The absolute (positive) volume change values are plotted in Figure 9 for better visibility and comparison of trends.
In both cases, active particles with different R s shrink in the same manner. However, for Ω c o n s t , the volume change is linear, and for Ω v a r , the volume change is nonlinear, as illustrated in Figure 9a,d, respectively. Moreover, as R s increases, the peak values of volume change decrease, which is consistent for both cases of partial molar volume. The peak values of Δ V / V o decrease because the stop condition for larger particles is achieved earlier, and the maximum volume achieved decreases. This indicates that the nonlinear volume change will have minimal effects for larger particles.
The Δ SOC curves in Figure 9b,e show that as R s increases, the magnitude of Δ SOC increases significantly. Larger particles take longer to start delithiating the core, resulting in enhanced concentration gradients within the particle that remain for a longer period of time. As a consequence, Δ SOC increases. However, the trends and magnitude for both Ω c o n s t and Ω v a r remain almost the same, suggesting that the concentration-dependent partial molar volume has a minimal effect on Δ SOC with increasing particle size.
Meanwhile, as the concentration gradients increase, σ m a x also increases, as shown in Figure 9c,f. Although the evolution in time and magnitude of σ m a x for both the constants Ω c o n s t and Ω v a r are different, the trend in the rise of the peak value of σ m a x with increasing R s is the same. For larger particles, the difference between the peak values of σ m a x is due to a very small volume change. In summary, the trends in the evolution of Δ V / V o , Δ SOC , and σ m a x remain consistent, while R s increases for both the Ω c o n s t and Ω v a r cases. However, the propensity for mechanical fracture increases with increasing particle size. Meanwhile the effects of concentration-dependent partial molar volume on the stress rise decrease for larger particles.

3.4. Effects of Charge Rate

Fast charging is the most desirable aspect of EVs. However, fast charging significantly increases the capacity degradation of lithium-ion batteries. Therefore, to explore the effects of fast charging in combination with the use of a concentration-dependent partial molar volume, the C rate was varied from 0.5 C to 10 C with a difference of 0.5 C. Simulations were performed for R s = 2   µ m of the NMC-111-active particle.
Figure 10 shows that the evolution trends of Δ V / V o , Δ SOC , and σ m a x are consistent for different C rates , even when Ω v a r was used. As the stopping condition of the numerical simulations was achieved earlier with an increasing C rate , the peak magnitude of Δ V / V o decreases, as presented in Figure 10a,d for Ω c o n s t and Ω v a r , respectively. Although the trends are consistent, the evolution of Δ V / V o becomes nonlinear once the concentration-dependent partial molar volume is used.
Furthermore, increasing the C rate causes an increase in Δ SOC due to enhanced concentration gradients. Figure 10b,e show a slight increase in the Δ SOC values for Ω v a r compared to Ω c o n s t . This is caused by this stress–concentration coupling effect.
Moreover, small SOC differences for a lower C rate indirectly indicate a more uniform distribution of lithium inside the particles and thus exhibit reduced concentration gradients. Since the lithium diffusion-induced strain mismatch is the main cause of stress generation, with smaller SOC differences, stress increase is also reduced for the lower C rate . As the C rate increases, σ m a x in a given delithiation time increases. The increase in the stress raises the probability of fracture. Although trends in the evolution of σ m a x remain consistent for Ω c o n s t and Ω v a r , using Ω v a r , the stress rise is significantly higher, as shown in Figure 10c,f for Ω c o n s t and Ω v a r , respectively. Thus, the propensity for mechanical failure increases significantly when using a concentration-dependent partial molar volume under fast charging conditions.

4. Conclusions

In this paper, the influence of a concentration-dependently local volume change of active materials on the lithium concentration evolution and diffusion-induced stress generation was studied by performing finite element simulations using a fully coupled chemomechanical model at the particle level. The concentration-dependent partial molar volume was calculated based on the previously obtained local nonlinear volume change of the NMC-active materials. We first validated the calculations of the partial molar volume by comparing the results of the global volume change obtained by finite element simulations with the local volume change values that were previously received in experiments described in the literature. We then compared the chemical and mechanical response of the active material with and without considering the concentration-dependent partial molar volume.
The main findings of this work are given below:
1.
The local volume change induced by the concentration-dependent chemical expansion of the active material significantly alters the global volume change of the active particles, which suggests that the stress increase due to the surrounding materials in electrodes will be affected by the concentration-dependent partial molar volume of the active materials.
2.
The concentration-dependent partial molar volume significantly changes the stress evolution trends and SOC differences. The peak stress due to diffusion is almost three times greater for a variable partial molar volume. However, the accumulated capacity within the particle remains independent of the change in partial molar volume.
3.
The trends of the maximum diffusion-induced stress in the particle correlate with variations in the time rate of Δ V / V o change (total volume change of particle). Although changes in the partial molar volume affect the stress response of the active particles, the lithium concentration patterns and stress distributions remain the same.
4.
As the particle size increases, the propensity for mechanical failure increases with the use of the concentration-dependent partial molar volume. However, the effect of changing the partial molar volume on the stress rise decreases for larger particles.
5.
Faster charging with a concentration-dependent partial molar volume increases diffusion-induced stress levels compared to using a constant partial molar volume.

Author Contributions

Conceptualization, N.I. and S.L.; methodology, N.I.; software, N.I., J.C., C.L., D.M., Y.K., J.K. and M.K.; validation, N.I., H.M.U.A., J.C., C.L., D.M., Y.K., J.K. and M.K.; investigation, N.I., J.C., C.L., D.M., Y.K., J.K. and M.K.; resources, S.L.; writing—original draft preparation, N.I. and J.C.; writing—review and editing, N.I., H.M.U.A. and S.L.; visualization, N.I. and H.M.U.A.; supervision, S.L.; project administration, S.L.; funding acquisition, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Research Foundation of Korea grants, funded by the Ministry of Science and ICT (No. 2018R1A5A7023490 and No. 2022R1A2C1003003).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ding, Y.; Cano, Z.P.; Yu, A.; Lu, J.; Chen, Z. Automotive Li-Ion Batteries: Current Status and Future Perspectives. Electrochem. Energ. Rev. 2019, 2, 1–28. [Google Scholar] [CrossRef]
  2. Nie, L.; Chen, S.; Liu, W. Challenges and strategies of lithium-rich layered oxides for Li-ion batteries. Nano Res. 2022, 1–12. [Google Scholar] [CrossRef]
  3. Armand, M.; Tarascon, J.-M. Building better batteries. Nature 2008, 451, 652–657. [Google Scholar] [CrossRef] [PubMed]
  4. Lu, Z.; MacNeil, D.D.; Dahn, J.R. Layered cathode materials Li [Ni x Li (1/3− 2x/3) Mn (2/3− x/3)] O 2 for lithium-ion batteries. Electrochem. Solid-State Lett. 2001, 4, A191. [Google Scholar] [CrossRef]
  5. Yabuuchi, N.; Ohzuku, T. Novel lithium insertion material of LiCo1/3Ni1/3Mn1/3O2 for advanced lithium-ion batteries. J. Power Sources 2003, 119–121, 171–174. [Google Scholar] [CrossRef]
  6. Iqbal, N.; Choi, J.; Lee, C.; Khan, A.; Tanveer, M.; Lee, S. A Review on Modeling of Chemo-mechanical Behavior of Particle—Binder Systems in Lithium-Ion Batteries. Multiscale Sci. Eng. 2022, 4, 79–93. [Google Scholar] [CrossRef]
  7. Ali, Y.; Iqbal, N.; Lee, S. Inhomogeneous stress development at the multiparticle electrode of lithium-ion batteries. Int. J. Energy Res. 2021, 45, 14788–14803. [Google Scholar] [CrossRef]
  8. Christensen, J.; Newman, J. Stress generation and fracture in lithium insertion materials. J. Solid State Electrochem. 2006, 10, 293–319. [Google Scholar] [CrossRef]
  9. Zhang, X.; Shyy, W.; Marie Sastry, A. Numerical Simulation of Intercalation-Induced Stress in Li-Ion Battery Electrode Particles. J. Electrochem. Soc. 2007, 154, A910. [Google Scholar] [CrossRef]
  10. Deshpande, R.; Cheng, Y.-T.; Verbrugge, M.W.; Timmons, A. Diffusion Induced Stresses and Strain Energy in a Phase-Transforming Spherical Electrode Particle. J. Electrochem. Soc. 2011, 158, A718. [Google Scholar] [CrossRef]
  11. Park, J.; Lu, W.; Sastry, A.M. Numerical Simulation of Stress Evolution in Lithium Manganese Dioxide Particles due to Coupled Phase Transition and Intercalation. J. Electrochem. Soc. 2011, 158, A201. [Google Scholar] [CrossRef]
  12. Bai, Y.; Zhao, K.; Liu, Y.; Stein, P.; Xu, B.-X. A chemo-mechanical grain boundary model and its application to understand the damage of Li-ion battery materials. Scr. Mater. 2020, 183, 45–49. [Google Scholar] [CrossRef]
  13. Barai, P.; Higa, K.; Ngo, A.T.; Curtiss, L.A.; Srinivasan, V. Mechanical Stress Induced Current Focusing and Fracture in Grain Boundaries. J. Electrochem. Soc. 2019, 166, A1752–A1762. [Google Scholar] [CrossRef]
  14. Han, S.; Park, J.; Lu, W.; Sastry, A.M. Numerical study of grain boundary effect on Li+ effective diffusivity and intercalation-induced stresses in Li-ion battery active materials. J. Power Sources 2013, 240, 155–167. [Google Scholar] [CrossRef]
  15. Zhao, K.; Pharr, M.; Vlassak, J.J.; Suo, Z. Fracture of electrodes in lithium-ion batteries caused by fast charging. J. Appl. Phys. 2010, 108, 073517. [Google Scholar] [CrossRef]
  16. Lee, S.; Park, J.; Sastry, A.M.; Lu, W. Molecular Dynamics Simulations of SOC-Dependent Elasticity of LixMn2O4 Spinels in Li-Ion Batteries. J. Electrochem. Soc. 2013, 160, A968–A972. [Google Scholar] [CrossRef]
  17. Schmidt, A.P.; Bitzer, M.; Imre, Á.W.; Guzzella, L. Experiment-driven electrochemical modeling and systematic parameterization for a lithium-ion battery cell. J. Power Sources 2010, 195, 5071–5080. [Google Scholar] [CrossRef]
  18. Stein, P.; Xu, B. 3D Isogeometric Analysis of intercalation-induced stresses in Li-ion battery electrode particles. Comput. Methods Appl. Mech. Eng. 2014, 268, 225–244. [Google Scholar] [CrossRef]
  19. Singh, G.; Bhandakkar, T.K. Semianalytical study of the effect of realistic boundary conditions on diffusion induced stresses in cylindrical Lithium ion electrode-binder system. Int. J. Mech. Sci. 2019, 163, 105141. [Google Scholar] [CrossRef]
  20. Stein, P.; Moradabadi, A.; Diehm, M.; Xu, B.-X.; Albe, K. The influence of anisotropic surface stresses and bulk stresses on defect thermodynamics in LiCoO2 nanoparticles. Acta Mater. 2018, 159, 225–240. [Google Scholar] [CrossRef]
  21. Pharr, M.; Suo, Z.; Vlassak, J.J. Variation of stress with charging rate due to strain-rate sensitivity of silicon electrodes of Li-ion batteries. J. Power Sources 2014, 270, 569–575. [Google Scholar] [CrossRef]
  22. Ma, Z.S.; Xie, Z.C.; Wang, Y.; Zhang, P.P.; Pan, Y.; Zhou, Y.C.; Lu, C. Failure modes of hollow core–shell structural active materials during the lithiation–delithiation process. J. Power Sources 2015, 290, 114–122. [Google Scholar] [CrossRef]
  23. Li, J.; Fang, Q.; Wu, H.; Liu, Y.; Wen, P. Investigation into diffusion induced plastic deformation behavior in hollow lithium ion battery electrode revealed by analytical model and atomistic simulation. Electrochim. Acta 2015, 178, 597–607. [Google Scholar] [CrossRef]
  24. Wen, J.; Wei, Y.; Cheng, Y.-T. Stress evolution in elastic-plastic electrodes during electrochemical processes: A numerical method and its applications. J. Mech. Phys. Solids 2018, 116, 403–415. [Google Scholar] [CrossRef] [Green Version]
  25. Colclasure, A.M.; Smith, K.A.; Kee, R.J. Modeling detailed chemistry and transport for solid-electrolyte-interface (SEI) films in Li–ion batteries. Electrochim. Acta 2011, 58, 33–43. [Google Scholar] [CrossRef]
  26. Deng, J.; Wagner, G.J.; Muller, R.P. Phase Field Modeling of Solid Electrolyte Interface Formation in Lithium Ion Batteries. J. Electrochem. Soc. 2013, 160, A487–A496. [Google Scholar] [CrossRef]
  27. Liu, L.; Park, J.; Lin, X.; Sastry, A.M.; Lu, W. A thermal-electrochemical model that gives spatial-dependent growth of solid electrolyte interphase in a Li-ion battery. J. Power Sources 2014, 268, 482–490. [Google Scholar] [CrossRef]
  28. Gao, E.; Lu, B.; Zhao, Y.; Feng, J.; Song, Y.; Zhang, J. Stress-Regulated Protocols for Fast Charging and Long Cycle Life in Lithium-Ion Batteries: Modeling and Experiments. J. Electrochem. Soc. 2021, 168, 060549. [Google Scholar] [CrossRef]
  29. Iqbal, N.; Haq, I.U.; Lee, S. Chemo-mechanical model predicted critical SOCs for the mechanical stability of electrode materials in lithium-ion batteries. Int. J. Mech. Sci. 2022, 216, 107034. [Google Scholar] [CrossRef]
  30. Iqbal, N.; Lee, S. Stress-regulated pulse charging protocols via coupled electrochemical-mechanical model for the mechanical stability of electrode materials in lithium-ion batteries. J. Power Sources 2022, 536, 231376. [Google Scholar] [CrossRef]
  31. Iqbal, N.; Ali, Y.; Haq, I.U.; Lee, S. Progressive interface debonding in composite electrodes of Li-ion batteries via mixed-mode cohesive zone model: Effects of binder characteristics. Compos. Struct. 2021, 259, 113173. [Google Scholar] [CrossRef]
  32. Klinsmann, M.; Rosato, D.; Kamlah, M.; McMeeking, R.M. Modeling crack growth during Li extraction and insertion within the second half cycle. J. Power Sources 2016, 331, 32–42. [Google Scholar] [CrossRef]
  33. Xu, R.; Yang, Y.; Yin, F.; Liu, P.; Cloetens, P.; Liu, Y.; Lin, F.; Zhao, K. Heterogeneous damage in Li-ion batteries: Experimental analysis and theoretical modeling. J. Mech. Phys. Solids 2019, 129, 160–183. [Google Scholar] [CrossRef]
  34. Guo, Z.-S.; Zhang, T.; Zhu, J.; Wang, Y. Effects of hydrostatic pressure and modulus softening on electrode curvature and stress in a bilayer electrode plate. Comput. Mater. Sci. 2014, 94, 218–224. [Google Scholar] [CrossRef]
  35. He, Y.-L.; Hu, H.J.; Song, Y.-C.; Guo, Z.-S.; Liu, C.; Zhang, J.-Q. Effects of concentration-dependent elastic modulus on the diffusion of lithium ions and diffusion induced stress in layered battery electrodes. J. Power Sources 2014, 248, 517–523. [Google Scholar] [CrossRef]
  36. Qi, Y.; Hector, L.G.; James, C.; Kim, K.J. Lithium Concentration Dependent Elastic Properties of Battery Electrode Materials from First Principles Calculations. J. Electrochem. Soc. 2014, 161, F3010. [Google Scholar] [CrossRef] [Green Version]
  37. Sitinamaluwa, H.S.; Wang, M.C.; Will, G.; Senadeera, W.; Zhang, S.; Yan, C. Lithium concentration dependent structure and mechanics of amorphous silicon. J. Appl. Phys. 2016, 119, 245103. [Google Scholar] [CrossRef] [Green Version]
  38. Xie, H.; Kang, Y.; Song, H.; Zhang, Q. Real-time measurements and experimental analysis of material softening and total stresses of Si-composite electrode. J. Power Sources 2019, 424, 100–107. [Google Scholar] [CrossRef]
  39. Yang, B.; He, Y.-P.; Irsa, J.; Lundgren, C.A.; Ratchford, J.B.; Zhao, Y.-P. Effects of composition-dependent modulus, finite concentration and boundary constraint on Li-ion diffusion and stresses in a bilayer Cu-coated Si nano-anode. J. Power Sources 2012, 204, 168–176. [Google Scholar] [CrossRef]
  40. Zhang, T.; Guo, Z. Effects of electrode properties and fabricated pressure on Li ion diffusion and diffusion-induced stresses in cylindrical Li-ion batteries. Model. Simul. Mater. Sci. Eng. 2014, 22, 025016. [Google Scholar] [CrossRef]
  41. Zhang, Y.; Zhan, S.; Zhang, K.; Zheng, B.; Lyu, L. Buckling behavior of a wire-like electrode with a concentration-dependent elastic modulus based on a deformed configuration. Eur. J. Mech.-A/Solids 2021, 85, 104111. [Google Scholar] [CrossRef]
  42. Liu, Y.; Guo, K.; Wang, C.; Han, J.; Gao, H. Concentration dependent properties and plastic deformation facilitate instability of the solid-electrolyte interphase in Li-ion batteries. Int. J. Solids Struct. 2020, 198, 99–109. [Google Scholar] [CrossRef]
  43. Weng, L.; Xu, C.; Chen, B.; Zhou, J.; Cai, R.; Shi, Y. A comparative study on ratcheting deformation between negative Poisson’s ratio electrode and thin film electrode in Li-ion battery cyclic operation. Mech. Mater. 2020, 150, 103567. [Google Scholar] [CrossRef]
  44. Hong, C.S.; Qaiser, N.; Nam, H.G.; Han, S.M. Effect of Li concentration-dependent material properties on diffusion induced stresses of a Sn anode. Phys. Chem. Chem. Phys. 2019, 21, 9581–9589. [Google Scholar] [CrossRef]
  45. Dora, J.K.; Sengupta, A.; Ghosh, S.; Yedla, N.; Chakraborty, J. Stress evolution with concentration-dependent compositional expansion in a silicon lithium-ion battery anode particle. J. Solid State Electrochem. 2019, 23, 2331–2342. [Google Scholar] [CrossRef]
  46. Malavé, V.; Berger, J.R.; Martin, P.A. Concentration-Dependent Chemical Expansion in Lithium-Ion Battery Cathode Particles. J. Appl. Mech. 2014, 81, 091005. [Google Scholar] [CrossRef]
  47. Qian, Y.; Lu, B.; Bao, Y.; Zhao, Y.; Song, Y.; Zhang, J. Prelithiation design for suppressing delamination in lithium-ion battery electrodes. Appl. Math. Mech.-Engl. Ed. 2021, 42, 1703–1716. [Google Scholar] [CrossRef]
  48. Cai, X.; Guo, Z. Influence of Li Concentration-Dependent Diffusion Coefficient and Modulus Hardening on Diffusion-Induced Stresses in Anisotropic Anode Particles. J. Electrochem. Soc. 2021, 168, 010517. [Google Scholar] [CrossRef]
  49. Karami, A.; Nayebi, A. Kinematic hardening analysis of Li-ion battery with concentration-dependent material behaviours under cyclic charging and discharging. J. Power Sources 2020, 461, 228155. [Google Scholar] [CrossRef]
  50. Deshpande, R.; Qi, Y.; Cheng, Y.-T. Effects of Concentration-Dependent Elastic Modulus on Diffusion-Induced Stresses for Battery Applications. J. Electrochem. Soc. 2010, 157, A967–A971. [Google Scholar] [CrossRef]
  51. Zhang, K.; Li, Y.; Zheng, B. Effects of concentration-dependent elastic modulus on Li-ions diffusion and diffusion-induced stresses in spherical composition-gradient electrodes. J. Appl. Phys. 2015, 118, 105102. [Google Scholar] [CrossRef]
  52. de Biasi, L.; Kondrakov, A.O.; Geßwein, H.; Brezesinski, T.; Hartmann, P.; Janek, J. Between Scylla and Charybdis: Balancing Among Structural Stability and Energy Density of Layered NCM Cathode Materials for Advanced Lithium-Ion Batteries. J. Phys. Chem. C 2017, 121, 26163–26171. [Google Scholar] [CrossRef]
  53. Noh, H.-J.; Youn, S.; Yoon, C.S.; Sun, Y.-K. Comparison of the structural and electrochemical properties of layered Li [NixCoyMnz] O2 (x = 1/3, 0.5, 0.6, 0.7, 0.8 and 0.85) cathode material for lithium-ion batteries. J. Power Sources 2013, 233, 121–130. [Google Scholar] [CrossRef]
  54. Mao, Y.; Wang, X.; Xia, S.; Zhang, K.; Wei, C.; Bak, S.; Shadike, Z.; Liu, X.; Yang, Y.; Xu, R.; et al. High-Voltage Charging-Induced Strain, Heterogeneity, and Micro-Cracks in Secondary Particles of a Nickel-Rich Layered Cathode Material. Adv. Funct. Mater. 2019, 29, 1900247. [Google Scholar] [CrossRef]
  55. Iqbal, N.; Ali, Y.; Lee, S. Chemo-mechanical response of composite electrode systems with multiple binder connections. Electrochim. Acta 2020, 364, 137312. [Google Scholar] [CrossRef]
  56. Takahashi, K.; Higa, K.; Mair, S.; Chintapalli, M.; Balsara, N.; Srinivasan, V. Mechanical degradation of graphite/PVDF composite electrodes: A model-experimental study. J. Electrochem. Soc. 2016, 163, A385–A395. [Google Scholar] [CrossRef] [Green Version]
  57. Ryu, I.; Choi, J.W.; Cui, Y.; Nix, W.D. Size-dependent fracture of Si nanowire battery anodes. J. Mech. Phys. Solids 2011, 59, 1717–1730. [Google Scholar] [CrossRef]
  58. Sun, H.; Zhao, K. Electronic structure and comparative properties of LiNi x Mn y Co z O2 cathode materials. J. Phys. Chem. C 2017, 121, 6002–6010. [Google Scholar] [CrossRef]
  59. Wei, Y.; Zheng, J.; Cui, S.; Song, X.; Su, Y.; Deng, W.; Wu, Z.; Wang, X.; Wang, W.; Rao, M.; et al. Kinetics tuning of Li-ion diffusion in layered Li (Ni x Mn y Co z) O2. J. Am. Chem. Soc. 2015, 137, 8364–8367. [Google Scholar] [CrossRef]
  60. Huang, X.; Zhu, W.; Yao, J.; Bu, L.; Li, X.; Tian, K.; Lu, H.; Quan, C.; Xu, S.; Xu, K.; et al. Suppressing structural degradation of Ni-rich cathode materials towards improved cycling stability enabled by a Li2 MnO3 coating. J. Mater. Chem. A 2020, 8, 17429–17441. [Google Scholar] [CrossRef]
  61. Cheng, E.J.; Hong, K.; Taylor, N.J.; Choe, H.; Wolfenstine, J.; Sakamoto, J. Mechanical and physical properties of LiNi0. 33Mn0. 33Co0. 33O2 (NMC). J. Eur. Ceram. Soc. 2017, 37, 3213–3217. [Google Scholar] [CrossRef]
  62. Mistry, A.; Juarez-Robles, D.; Stein, M.; Smith, K.; Mukherjee, P.P. Analysis of long-range interaction in lithium-ion battery electrodes. J. Electrochem. Energy Convers. Storage 2016, 13, 031006. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic representation and (b) spatial discretization of an axisymmetric model system used to perform numerical simulations.
Figure 1. (a) Schematic representation and (b) spatial discretization of an axisymmetric model system used to perform numerical simulations.
Mathematics 10 04697 g001
Figure 2. Calculated partial molar volume as the function of normalized lithium concentrations or lithium contents.
Figure 2. Calculated partial molar volume as the function of normalized lithium concentrations or lithium contents.
Mathematics 10 04697 g002
Figure 3. Local and global volume changes obtained by experiments (Adapted with permission from Ref. [52]. Copyright 2017, American Chemical Society.) and finite element simulations, respectively. The local volume change is plotted against the lithium content, and the global volume change is plotted against the normalized capacity of the various NMC materials.
Figure 3. Local and global volume changes obtained by experiments (Adapted with permission from Ref. [52]. Copyright 2017, American Chemical Society.) and finite element simulations, respectively. The local volume change is plotted against the lithium content, and the global volume change is plotted against the normalized capacity of the various NMC materials.
Mathematics 10 04697 g003
Figure 4. Validation of the use of concentration-dependent partial molar volume by performing simulations by intercalation strain and thermal analogy methods. (a) change in the total volume of the active particle ( Δ V / V o ), (b) the maximum value of the local first principal stress ( σ m a x ), (c) the difference between maximum and minimum local charge states ( Δ SOC ), and (d) the normalized lithium charge capacity ( Π ) inside the active particle.
Figure 4. Validation of the use of concentration-dependent partial molar volume by performing simulations by intercalation strain and thermal analogy methods. (a) change in the total volume of the active particle ( Δ V / V o ), (b) the maximum value of the local first principal stress ( σ m a x ), (c) the difference between maximum and minimum local charge states ( Δ SOC ), and (d) the normalized lithium charge capacity ( Π ) inside the active particle.
Mathematics 10 04697 g004
Figure 5. Comparison of the results of the delithiation process simulation with and without considering the concentration-dependent partial molar volume. (a) change in the total volume of the active particle ( Δ V / V o ), (b) the maximum value of the local first principal stress ( σ m a x ), (c) the difference between the maximum and minimum local states of charge ( Δ SOC ), and (d) the normalized lithium charge capacity ( Π ) inside the active particle.
Figure 5. Comparison of the results of the delithiation process simulation with and without considering the concentration-dependent partial molar volume. (a) change in the total volume of the active particle ( Δ V / V o ), (b) the maximum value of the local first principal stress ( σ m a x ), (c) the difference between the maximum and minimum local states of charge ( Δ SOC ), and (d) the normalized lithium charge capacity ( Π ) inside the active particle.
Mathematics 10 04697 g005
Figure 6. Comparison of results for uncoupled simulations. (a) The difference between the maximum and minimum local state of charges ( Δ SOC ). (b) The maximum value of the local first principal stress ( σ m a x ).
Figure 6. Comparison of results for uncoupled simulations. (a) The difference between the maximum and minimum local state of charges ( Δ SOC ). (b) The maximum value of the local first principal stress ( σ m a x ).
Mathematics 10 04697 g006
Figure 7. The delithiation time evolution of the maximum stress ( σ m a x ) and the 1st derivative of the volume change ( Δ V / V o ).
Figure 7. The delithiation time evolution of the maximum stress ( σ m a x ) and the 1st derivative of the volume change ( Δ V / V o ).
Mathematics 10 04697 g007
Figure 8. (a) Representative contour plots for (a,b) lithium concentration distributions ( c / c t o t a l ) and (c,d) first principal stress at the end of the delithiation. The results (a,c) are based on a constant partial molar volume, and (b,d) were obtained using a concentration-dependent partial molar volume.
Figure 8. (a) Representative contour plots for (a,b) lithium concentration distributions ( c / c t o t a l ) and (c,d) first principal stress at the end of the delithiation. The results (a,c) are based on a constant partial molar volume, and (b,d) were obtained using a concentration-dependent partial molar volume.
Mathematics 10 04697 g008
Figure 9. Impact of particle size (Rs) on the transient chemomechanical response during the delithiation process. The results (ac) were obtained using a constant partial molar volume, and the results (df) were obtained using a variable partial molar volume. The red dots to the right of each graph denote the evolution of the peak values with respect to the particle size increase. (a,d) Evolution of the volume change ( Δ V / V o ) of particle size. (b,e) Evolution of the maximum difference in this state of charge ( Δ SOC ). (c,f) Evolution of the first principle maximum stress ( σ m a x ) in the active particle.
Figure 9. Impact of particle size (Rs) on the transient chemomechanical response during the delithiation process. The results (ac) were obtained using a constant partial molar volume, and the results (df) were obtained using a variable partial molar volume. The red dots to the right of each graph denote the evolution of the peak values with respect to the particle size increase. (a,d) Evolution of the volume change ( Δ V / V o ) of particle size. (b,e) Evolution of the maximum difference in this state of charge ( Δ SOC ). (c,f) Evolution of the first principle maximum stress ( σ m a x ) in the active particle.
Mathematics 10 04697 g009
Figure 10. Impact of the charging rate ( C rate ) on the transient chemomechanical response during the delithiation process. The results (ac) were obtained using a constant partial molar volume, and the results (df) were obtained using a variable partial molar volume. The red dots to the right of each graph denote the evolution of the respective peak values. (a,d) Evolution of volume change ( Δ V / V o ) of the particle size. (b,e) Evolution of the maximum difference in this state of charge ( Δ SOC ). (c,f) Evolution of the first principle maximum stress ( σ m a x ) in the active particle.
Figure 10. Impact of the charging rate ( C rate ) on the transient chemomechanical response during the delithiation process. The results (ac) were obtained using a constant partial molar volume, and the results (df) were obtained using a variable partial molar volume. The red dots to the right of each graph denote the evolution of the respective peak values. (a,d) Evolution of volume change ( Δ V / V o ) of the particle size. (b,e) Evolution of the maximum difference in this state of charge ( Δ SOC ). (c,f) Evolution of the first principle maximum stress ( σ m a x ) in the active particle.
Mathematics 10 04697 g010
Table 1. Summary of simulation parameters and material properties for NMC materials used for numerical simulations.
Table 1. Summary of simulation parameters and material properties for NMC materials used for numerical simulations.
ParametersSymbolsUnitsValues
NMC-111NMC-523NMC-622NMC-811
Young’s modulus 1 E GPa202.98191.79181.52194.4
Diffusion coefficient 2 D m2 s−13.39 × 10−153.89 × 10−157.5 × 10−154.0 × 10−14
Specific capacity 3 α mAh g−1188.75194.89203.18213.42
Stoichiometric lithium concentration 4 c t o t a l mol m−333,45234,54236,00937,825
Minimum SOC 5 SOC min %21191811
Maximum SOC 6 SOC max %94939390
Poisson’s ratio 7 ν -0.25
Density 8 ρ kg m−34750
Faraday’s constant F C mol−196,487
Absolute temperature T K300
Universal gas constant R J mol−1 K−18.314
1 Sun and Zhao [58]; 2 Wei et al. [59] and Huang et al. [60]; 3,5,6 de Biasi et al. [52]; 4 Calculated by c t o t a l = α ρ / F ; 7 Cheng et al. [61]; 8 Mistry et al. [62].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Iqbal, N.; Choi, J.; Lee, C.; Ayub, H.M.U.; Kim, J.; Kim, M.; Kim, Y.; Moon, D.; Lee, S. Effects of Diffusion-Induced Nonlinear Local Volume Change on the Structural Stability of NMC Cathode Materials of Lithium-Ion Batteries. Mathematics 2022, 10, 4697. https://doi.org/10.3390/math10244697

AMA Style

Iqbal N, Choi J, Lee C, Ayub HMU, Kim J, Kim M, Kim Y, Moon D, Lee S. Effects of Diffusion-Induced Nonlinear Local Volume Change on the Structural Stability of NMC Cathode Materials of Lithium-Ion Batteries. Mathematics. 2022; 10(24):4697. https://doi.org/10.3390/math10244697

Chicago/Turabian Style

Iqbal, Noman, Jinwoong Choi, Changkyu Lee, Hafiz Muhammad Uzair Ayub, Jinho Kim, Minseo Kim, Younggee Kim, Dongjae Moon, and Seungjun Lee. 2022. "Effects of Diffusion-Induced Nonlinear Local Volume Change on the Structural Stability of NMC Cathode Materials of Lithium-Ion Batteries" Mathematics 10, no. 24: 4697. https://doi.org/10.3390/math10244697

APA Style

Iqbal, N., Choi, J., Lee, C., Ayub, H. M. U., Kim, J., Kim, M., Kim, Y., Moon, D., & Lee, S. (2022). Effects of Diffusion-Induced Nonlinear Local Volume Change on the Structural Stability of NMC Cathode Materials of Lithium-Ion Batteries. Mathematics, 10(24), 4697. https://doi.org/10.3390/math10244697

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop