Research article Special Issues

Enhancement of cone beam CT image registration by super-resolution pre-processing algorithm


  • In order to enhance cone-beam computed tomography (CBCT) image information and improve the registration accuracy for image-guided radiation therapy, we propose a super-resolution (SR) image enhancement method. This method uses super-resolution techniques to pre-process the CBCT prior to registration. Three rigid registration methods (rigid transformation, affine transformation, and similarity transformation) and a deep learning deformed registration (DLDR) method with and without SR were compared. The five evaluation indices, the mean squared error (MSE), mutual information, Pearson correlation coefficient (PCC), structural similarity index (SSIM), and PCC + SSIM, were used to validate the results of registration with SR. Moreover, the proposed method SR-DLDR was also compared with the VoxelMorph (VM) method. In rigid registration with SR, the registration accuracy improved by up to 6% in the PCC metric. In DLDR with SR, the registration accuracy was improved by up to 5% in PCC + SSIM. When taking the MSE as the loss function, the accuracy of SR-DLDR is equivalent to that of the VM method. In addition, when taking the SSIM as the loss function, the registration accuracy of SR-DLDR is 6% higher than that of VM. SR is a feasible method to be used in medical image registration for planning CT (pCT) and CBCT. The experimental results show that the SR algorithm can improve the accuracy and efficiency of CBCT image alignment regardless of which alignment algorithm is used.

    Citation: Liwei Deng, Yuanzhi Zhang, Jingjing Qi, Sijuan Huang, Xin Yang, Jing Wang. Enhancement of cone beam CT image registration by super-resolution pre-processing algorithm[J]. Mathematical Biosciences and Engineering, 2023, 20(3): 4403-4420. doi: 10.3934/mbe.2023204

    Related Papers:

    [1] Fahmi Astuti, Rima Feisy Azmi, Mohammad Arrafi Azhar, Fani Rahayu Hidayah Rayanisaputri, Muhammad Redo Ramadhan, Malik Anjelh Baqiya, Darminto . Employing Na2CO3 and NaCl as sources of sodium in NaFePO4 cathode: A comparative study on structure and electrochemical properties. AIMS Materials Science, 2024, 11(1): 102-113. doi: 10.3934/matersci.2024006
    [2] Christian M Julien, Alain Mauger, Ashraf E Abdel-Ghany, Ahmed M Hashem, Karim Zaghib . Smart materials for energy storage in Li-ion batteries. AIMS Materials Science, 2016, 3(1): 137-148. doi: 10.3934/matersci.2016.1.137
    [3] Alain Mauger, Haiming Xie, Christian M. Julien . Composite anodes for lithium-ion batteries: status and trends. AIMS Materials Science, 2016, 3(3): 1054-1106. doi: 10.3934/matersci.2016.3.1054
    [4] Christian M. Julien, Alain Mauger . Functional behavior of AlF3 coatings for high-performance cathode materials for lithium-ion batteries. AIMS Materials Science, 2019, 6(3): 406-440. doi: 10.3934/matersci.2019.3.406
    [5] Tomoya Kawabata, Shuji Aihara, Yukito Hagihara . Coalescence judgment criteria for the interaction between two close surface cracks by WES2805 and its safety margin for brittle fracture assessment. AIMS Materials Science, 2016, 3(4): 1665-1682. doi: 10.3934/matersci.2016.4.1665
    [6] Mordechai Perl, Tomer Saley . The cumulative detrimental impact of pressure and autofrettage on the fatigue life of an externally cracked modern tank gun barrel. AIMS Materials Science, 2019, 6(5): 833-851. doi: 10.3934/matersci.2019.5.833
    [7] Hiroki Nara, Keisuke Morita, Tokihiko Yokoshima, Daikichi Mukoyama, Toshiyuki Momma, Tetsuya Osaka . Electrochemical impedance spectroscopy analysis with a symmetric cell for LiCoO2 cathode degradation correlated with Co dissolution. AIMS Materials Science, 2016, 3(2): 448-459. doi: 10.3934/matersci.2016.2.448
    [8] Ojo Sunday Isaac Fayomi, Adedamola Sode, Itopa Godwin Akande, Abimbola Patricia Idowu Popoola, Oluranti Agboola . Improving the structural properties and corrosion behaviour of electroless deposited Ni-P-Zn coatings on mild steel for advanced processes. AIMS Materials Science, 2020, 7(4): 441-452. doi: 10.3934/matersci.2020.4.441
    [9] Christian M. Julien, Alain Mauger . In situ Raman analyses of electrode materials for Li-ion batteries. AIMS Materials Science, 2018, 5(4): 650-698. doi: 10.3934/matersci.2018.4.650
    [10] Elisa Padovano, Francesco Trevisan, Sara Biamino, Claudio Badini . Processing of hybrid laminates integrating ZrB2/SiC and SiC layers. AIMS Materials Science, 2020, 7(5): 552-564. doi: 10.3934/matersci.2020.5.552
  • In order to enhance cone-beam computed tomography (CBCT) image information and improve the registration accuracy for image-guided radiation therapy, we propose a super-resolution (SR) image enhancement method. This method uses super-resolution techniques to pre-process the CBCT prior to registration. Three rigid registration methods (rigid transformation, affine transformation, and similarity transformation) and a deep learning deformed registration (DLDR) method with and without SR were compared. The five evaluation indices, the mean squared error (MSE), mutual information, Pearson correlation coefficient (PCC), structural similarity index (SSIM), and PCC + SSIM, were used to validate the results of registration with SR. Moreover, the proposed method SR-DLDR was also compared with the VoxelMorph (VM) method. In rigid registration with SR, the registration accuracy improved by up to 6% in the PCC metric. In DLDR with SR, the registration accuracy was improved by up to 5% in PCC + SSIM. When taking the MSE as the loss function, the accuracy of SR-DLDR is equivalent to that of the VM method. In addition, when taking the SSIM as the loss function, the registration accuracy of SR-DLDR is 6% higher than that of VM. SR is a feasible method to be used in medical image registration for planning CT (pCT) and CBCT. The experimental results show that the SR algorithm can improve the accuracy and efficiency of CBCT image alignment regardless of which alignment algorithm is used.



    1. Introduction

    Batteries have long been a very important commodity for satisfying the world’s mobile energy storage needs. But more than ever before, advanced lithium-ion batteries are receiving an increased level of interest due to expanded power capabilities and the near limitless assortment of application opportunities, from vehicles and laptop computers to orbiting satellites and electric tools. As transportation markets are coming into play, the environmental impact of lithium-ion batteries and durability characteristics against climate variation are becoming very important in both the component material selection and design [1,2,3,4,5]. In 1997, the concept of LiFePO4 as a possible cathode for today’s Li-ion batteries came to life [6]. This material, which has an olivine structure, has been a focal point for much experimentation and discussion, as it boasts several attractive qualities including a relatively high theoretical capacity of 170 mAh/g, great structural stability, long cycle life, environmentally benign qualities, and most importantly, safety due to great thermal stability [6,7]. However, LiFePO4 has a very low electronic conductivity, a property essential for battery design, but it has been greatly improved through the successful application of conductive coatings or cation doping, making it comparable with other cathode materials [8,9,10,11]. In addition to these attributes, the constituent materials are abundant and widely available at a relatively low cost [12], making LiFePO4 extremely attractive for high-power vehicles, military and space flight operations [13,14], and everyday consumer electronics.

    Phase transformation from lithium-poor (FePO4) to lithium-rich (LiFePO4) during intercalation induces different strains in each lattice direction: a, b, and c [15,16,17,18,19]. Stresses induced by the phase transformation create flaws or cracks that may become stress concentration areas when further cycling occurs [20,21,22,23,24,25,26]. As a result, cracks are often observed along the phase boundary on the ac-plane [15,27]; structural failure due to fracturing of active material is a primary factor in battery degradation [7,22,27,28,29,30]. A combination of tensile and compressive stresses will induce microscopic crack propagation if critical values are reached. Over time, crack propagation will continue until the crack reaches a critical length, at which time complete fracture initiates. When this active electrode material is fractured away, the battery will effectively lose capacity and maximum power output diminishes. Thus, understanding the relationship between stresses and imperfections during lithium ion intercalation can be beneficial to better advance the current battery technology.

    Olivine LiFePO4 is a brittle material and linear elastic fracture mechanics could be applied to better understand mixed-mode fractures in lithium-ion battery materials. Extensive theoretical and computational studies have focused on the growth of interfacial cracks of bi-materials, and on quantifying the extension of a crack and the elastic field near a crack tip [31,32]; the energy release rates (G) and stress intensity factors (K1 and K2) of engineering materials are generally discussed in this regard. Hutchinson and Suo have provided analytical energy release rates and stress intensity factors for interfacial cracks in isotropic bi-materials under mixed-mode loading [33]. Qian and Sun have extended the mixed-mode study to monoclinic and orthogonal bi-materials [34,35]. Sequentially, the development and application of the finite element-based virtual crack closure technique (VCCT) facilitated the analyses on fracture mechanics for orthogonal (orthorhombic) bi-materials [36,37,38,39,40]. Raju has calculated energy release rates via the VCCT with higher order and singular finite elements for isotropic materials [41]. Agrawal and Karlsson utilized the VCCT to obtain mixed-mode energy release rates and stress intensity factors for isotropic materials [40]. For details about the overview of the VCCT, please refer to studies by Krueger et al., Rybicki, Dattaguru, Xie and Agrawal [36,37,38,39,40].

    The current study aims to identify an optimal combination of the particle- and crack-size of LiFePO4 for better battery performance. Because reducing the LiFePO4 electrode particle size could effectively decrease stress inside the particles, in addition to reducing the diffusion path during electrochemical cycling [42,43]. The current study is designed to provide insight into the sub-micro scale fracture mechanisms within the cathode material by accomplishing the following objectives to help aid in the understanding and development of a longer lasting and more powerful battery: (1) Perform a finite element simulation to evaluate fracture mechanisms inside a single LiFePO4 particle due to the phase transformation; and (2) establish an energy-based approach to estimate the propensity of crack extension for lithium-ion battery materials.

    We achieve these objectives by utilizing the VCCT to study fracture mechanics in LiFePO4 battery materials under mixed-mode loading. Pre-existing cracks at the interface of orthotropic bi-materials are incorporated. Various particle and crack sizes are considered to illustrate the effects on the cracking information in LiFePO4 due to phase transformation. Our aim is to present an approach towards improving battery structural design by providing a better understanding of the micromechanics, as well as a basis for future fatigue analyses, which may incorporate a variety of other battery chemistries and additional fatigue life parameters.

    2. Materials and Method

    2.1. Model development

    To study fracture mechanics in LiFePO4 particles, VCCT with the ANSYS finite element software (ANSYS, Inc. Canonsburg, Pennsylvania) is utilized to obtain the energy release rates and stress intensity factors. Two-phase plate-like LiFePO4 finite element models are generated based on the separate experimental observations [15,27,44] and different orthotropic material properties are adopted for both phases [28]. An interfacial crack is considered to run parallel to the bc-plane (along the c-axis) where phase boundaries are present (Figure 1). Although crack extension is possible in other planes, literature and experimentation have reported that cracks are observed in the bc-plane between phases while advancing in the c-axis [15]. In the current study, an assumption was made that the cracks will extend along their original direction. Moreover, the only experimental data that was available was for the surface energy and the prediction of crack turning or kinking was beyond the scope of the investigation. However, there exists a slight possibility that the shear mode fracture toughness may be very small in these materials and the crack path may deviate from a straight line. Therefore, we assume that a crack extends along the interface and remains parallel to the bc-plane [15,27,29] (Figure 1). Based on the experimental observation [15], previous simulations [45,46,47,48], and LiFePO4 lattice constants of the unit cell—a = 10.3 Å, b = 6.0 Å, and c = 4.7 Å [49]—three differently sized plate-like LiFePO4 finite element models on the ac-plane are generated as follows: (a) 500 nm × 225 nm, (b) 200 nm × 90 nm, and (c) 100 nm × 45 nm, respectively. In general, a continuum mechanics-based finite element method is valid when the model size is above 25 nm, as several nano-scaled studies presented by Zhao et al. [22,50,51]. Furthermore, lithium-ions diffuse into LiFePO4 particles along the b-axis during the phase transformation; the elastic model is reduced to 2D by assuming plane stress and neglecting stresses in the b-direction, as described by Hu et al. [29]. A total of 21 finite element simulations are conducted, wherein the initial crack is considered to start from the top face with variable lengths of crack sizes 0.05-0.8 L/d, in which L represents the crack length and d represents the particle size along the c-axis (Figure 1). The initial crack mouth opening is set at 0.5 nm in the a-axis in the range of experimentally observed crack dimensions [27]. Strains were applied to the LiFePO4 phase according to the misfit strains previously measured and reported in the literature: expansion occurs along the a-direction (εa = 5.03%) and contraction occurs along the c-direction (εc = −1.9%) (Figure 1a) [52]. Of note, εb = 4.5%). Two-dimensional quadrilateral elements are used in the finite element analysis. Such eight-node parabolic elements allow for more flexibility and improved accuracy in contrast to four-node elements That is, the 8 noded parabolic quadrilateral element uses quadratic functions for the displacements, and the simulation shows the exact analytic solution for pure bending dominated problems even with a coarse mesh with only one element in depth [53,54]. The finite element models are densely meshed around the crack tip to a size of 50 nm (Figure 1b).

    Figure 1. (a) A representative orthotropic bi-material with variable lengths of crack sizes 0.05-0.8 L/d. Expansion occurs along the a-direction and contraction occurs along the c-direction when particles undergo phase changes [52]. (b) Two-dimensional finite element models are used in the analysis, and models are densely meshed around the crack tip to a size of 50 nm. Scale bar = 0.03 μm.

    2.2. Finite element method-based virtual crack closure technique (VCCT)

    The VCCT is used along with the finite element modeling package ANSYS to determine the energy release rates and stress intensity factors under mixed-mode loading. The theory behind the VCCT is that the energy needed to separate or slide a crack surface is the same as the energy needed to close the surface back on itself [36]. The tendency of a crack to extend, the energy release rates, and the stress intensity factors are calculated through command options in ANSYS [55] (Figure 2a), and they are currently not available on the graphical user interface of ANSYS. Stress intensity factors have been derived analytically for orthotropic bi-materials, as detailed in published studies [34,56,57]. Energy release rates can be calculated with relative displacements (u, w), reaction forces (X, Z), and a unit thickness t = 1, as shown in Equation 1 and Figure 2 [39]. From the VCCT and a few commands in ANSYS, the Mode 1 (tension mode) and Mode 2 (sliding shear mode) strain energy release rates for a flaw of given size in a two-phase LiFePO4 particle could be determined (Figure 2b-c):

    @{G_1} = \frac{{ - \left[{{X_1}\left( {\Delta {u_1}} \right) + {X_2}\left( {\Delta {u_2}} \right)} \right]}}{{2\Delta \alpha }} \qquad {G_2} = \frac{{ - \left[{{Z_1}\left( {\Delta {w_1}} \right) + {Z_2}\left( {\Delta {w_2}} \right)} \right]}}{{2\Delta \alpha }}@ (1)

    Figure 2. (a) The VCCT for 2-D quadrilateral elements with a unit thickness t. Energy release rates are calculated based on the crack size (Δa), reaction forces (Z and X), relative displacements (u, w, and a), and via the FEM-based VCCT from the command-coding interface of ANSYS. (b)-(c) Anticipated Mode 1 (tension mode) and Mode 2 (sliding shear mode) fractures for a two-phase LiFePO4 battery material under mixed-mode loading conditions.

    2.3. Crack extension

    A crack will advance if the total energy release rate for Mode 1 and Mode 2, GT= G1 + G2 , is larger than approximately twice the surface energy of the particle, i.e., GT > 2 g [29,58]. A first-principle analysis by Wang et al. has reported a g value of 0.66 N/m for LiFePO4 in the (100) crack face orientation [29,58]. In the current study, we collect finite element results satisfying GT > 2g, and the difference between these two energies (GT−2g) is used to predict crack extension. That is, crack extension is eminent and a crack will advance until the surface energy from the newly-formed crack faces bring the particle back to be equal to or larger than the strain energy release rate (Figure 3). In the current work, we follow steps in “Linear Elastic Fracture Mechanics” detailed by Fischer-Cripps [59]: for a crack to extend, the rate of strain energy release per unit of crack extension must be at least equal to the rate of surface energy requirement [60]. Therefore, a 0.01-nm crack extension is set as a unit growth and the resulting crack growth (da) would be determined as follows: First, we calculate the difference between the GT and 2g. Second, the additional energy release rate for a unit of 0.01 nm crack extension is determined based on the crack surface areas. Third, the updated additional surface energy for a 0.01 nm crack extension is determined. Fourth, the crack extension is calculated since the relieved internal energy due to crack propagation should be balanced by the additional strain energy release rate, as shown in Figure 3.

    Figure 3. Crack extension calculation. The argument GT > 2γ is used to determine the required crack advancement (da) to reach mechanical equilibrium, where b is the particle size in the b-axis and g = 0.66 N/m [29,58].

    3. Results and Discussion

    The results of the finite element method-based VCCT indicate that normal stresses dominate as compared to shear stress throughout the particle (stress fields not shown). In particular, the Mode 1 fractures predominated for all models, as the energy release rate G1 is much larger than G2 for all particle sizes (Figure 4). The fracture mode is a result of applied boundary conditions (misfit strains due to phase transformation in the current study: εa = 5.03% and εc = −1.9%), anisotropic elastic moduli for LiFePO4 and FePO4 [28], and as a result, the energy release rates and the stress intensity factors are indicators of the state of energy or stress around the crack tip. Since the stiffness value of LiFePO4 from [28] along the a-axis is much smaller than in the other orthorhombic directions, the material is weaker in this direction, giving way to the splitting or opening Mode 1 type of fracture. Therefore, in the current study we focus on the discussion of Mode 1 fractures in LiFePO4 battery materials from the mixed-mode boundary condition.

    Figure 4. Energy release rates for two modes of fracture at crack tips with various L/d ratios. G1 is highly dependent on the crack size for three different particle sizes. It is observed that Mode 1 cracks are likely to occur and has higher tendency to propagate, especially when L/d ≈ 0.5-0.6.

    Several studies have shown that the energy release rate is proportional to the crack size since the energy release rate is the energy dissipated during fracture per unit of newly created fracture surface area [33]. The variation of energy release rates with respect to the crack-to-particle ratios for two fracture modes are shown in Figure 4. It is revealed that G1is highly dependent on the crack size for all three particle sizes. Comparing to the 2g = 1.32 N/m, smaller particles (100 nm × 45 nm) are able to better accommodate initial flaws without the crack propagation, i.e., L/d ≤ 0.2. In contrast, particles with a larger size (500 nm × 225 nm) could only accommodate initial flaws as small as L/d ≤0.03 by a linear extrapolation. Therefore, it is the L/d point that needs to be considered to determine if crack propagation occurs during the lithiation/delithiation process. Below L/d = 0.5 and 0.6, the energy release rates increase with the increased crack lengths and then begin to level off. It is concluded that crack propagation acts to relieve internal stresses due to misfit strain during the phase transformation of LiFePO4 materials, and it is particle-size and initial flaw-size dependent. Crack propagation due to mixed-mode loading is possible in other planes; however, experimentation has reported that cracks are observed between phases while propagating in the c-axis [15]. The olivine crystal structure of LiFePO4 and the Pauling's third rule support this observation: (100) planes are linked through stable, shared corner bonds between FeO6 octahedra and PO4 tetrahedra, while (010) planes are connected by weaker, edge-sharing bonds between FeO6 octahedra [27,52,61]. Lower energy release rates are observed in the Mode 2 fracture and exhibit less dependency on particle sizes. However, G1 and G2 both contribute to the GT for the crack extension. In contrast to the study by Hu et al. [29], a maximum value of the normalized energy release rate was identified; however, it is not clear which fracture mode is considered in their study. The approach in their study was to inversely look for a critical particle size d by solving boundary value problems, and the normalized energy release rate is a function of one component of the anisotropic stiffness matrix (c11) and the volume misfit along the c-axis (ec). The decreasing normalized energy release rate after L/d = 0.5 from Hu et al. could be attributed by several simplifications [29], such as the choice of the normalization where the other anisotropic effects of volume expansions were not incorporated.

    Stress intensity factors (K) are also obtained from the command-coding interface of ANSYS via the VCCT (Figure 2a), and the particle-size-dependent relationship of stress intensity factors and L/d is shown in Figure 5a. For larger particles (500 nm × 225 nm), it is observed that K1reaches a relative maximum value of 1.4 MPa-m1/2. For a particle size of 200 nm × 90 nm, K1 reaches a relative maximum value of 1.0 MPa-m1/2. For smaller particles (100 nm × 45 nm), it is observed that K1 reaches a relative maximum value of 0.6 MPa-m1/2 (Figure 5a). Our results indicate that K1is particle-size dependent, which is consistent with the observations made by Krstic and Khaund [62].

    Figure 5. (a) Mode 1 stress intensity factors at crack tips with various L/d ratios. K1 is highly dependent on the particle and crack sizes. (b) Computational model predictions of crack advancement (Mode 1) for three different particle sizes. Note that these simulations are not based on Paris’ law, but rather from systematic analyses of independent VCCT models (i.e., 16 models are solved to generate the 16 data points plotted above). It is observed that smaller particles exhibit faster crack propagation (i.e., a higher value for the exponent of ΔK1 in the best-fit equations appearing in the plot above).

    Further, to reveal the relationship between the crack advancement (da) and the change in the Mode 1 stress intensity factor for LiFePO4, the estimation is depicted in Figure 5b. The da is determined via the method described in Figure 3, and ΔK1 is determined from ranges of the stress intensity factor of Mode 1. It is observed that da vs. ΔK1 curves are particle-size dependent, wherein smaller particles exhibit faster crack extension (i.e., a higher value for the exponent of ΔK1 in the best-fit equations appearing in Figure 5b): the crack growth is proportional to the DK1 to the power of 4.6229 for particle size of 100 nm × 45 nm, and for particle size of 500 nm x 225 nm, the crack growth is related to the DK1 to the power of 1.3849 (Figure 5b). Detailed calculated data from the FEM-based VCCT are listed in Table 1.

    Table 1. Results from the FEM-based VCCT of lithium-ion battery materials for each particle size with various crack-to-particle ratios: Mode 1 stress intensity factor (K1), total energy release rate (GT), twice of the surface energies (2γ), the energy difference between them (ΔE = GT - 2γ), the required additional strain energy release rate per 0.01 nm crack growth (ER), the required additional strain energy release rate per 0.01-nm crack growth (Es), and the required crack advancement (da).
    a × b × c = 500 nm × 300 nm × 225 nm
    L/dKI (MPa-m1/2)GT (N/m)2γ (N/m)ΔE (N/m)ER per 0.01 nm (Joules)ES per 0.01 nm (Joules)da (m)
    0.050.6595.4621.324.1422.485×10-173.960×10-186.275×10-11
    0.10.8335.9221.324.6022.761×10-173.960×10-186.972×10-11
    0.21.1368.7441.327.4244.454×10-173.960×10-181.125×10-11
    0.31.31011.0271.329.7075.824×10-173.960×10-181.471×10-11
    0.41.40312.3731.3211.0536.632×10-173.960×10-181.675×10-11
    0.61.45613.0531.3211.7337.040×10-173.960×10-181.778×10-11
    a × b × c = 200 nm × 120 nm × 90 nm
    L/dKI (MPa-m1/2)GT (N/m)2γ (N/m)ΔE (N/m)ER per 0.01 nm (Joules)ES per 0.01 nm (Joules)da (m)
    0.050.3540.0681.32-1.252---
    0.10.5310.0791.32-1.241---
    0.20.6743.4811.322.1615.185×10-181.584×10-183.273×10-11
    0.30.7794.4001.323.0807.393×10-181.584×10-184.667×10-11
    0.40.8364.8501.323.5308.472×10-181.584×10-185.348×10-11
    0.50.8615.1761.323.8569.254×10-181.584×10-185.842×10-11
    0.60.9025.2201.323.9009.361×10-181.584×10-185.910×10-11
    0.70.9035.2111.323.8919.337×10-181.584×10-185.895×10-11
    0.80.9265.3051.323.9859.564×10-181.584×10-186.038×10-11
    a × b × c = 100 nm × 60 nm × 45 nm
    L/dKI (MPa-m1/2)GT (N/m)2γ (N/m)ΔE (N/m)ER per 0.01 nm (Joules)ES per 0.01 nm (Joules)da (m)
    0.050.2970.0291.32-1.291--
    0.10.3750.0421.32-1.278--
    0.20.5101.7401.320.4205.044×10-197.920×10-196.369×10-12
    0.30.5892.2001.320.8801.055×10-197.920×10-191.332×10-12
    0.40.6292.4701.321.1501.380×10-197.920×10-191.742×10-12
    0.60.6522.6071.321.2871.545×10-197.920×10-191.951×10-12
     | Show Table
    DownLoad: CSV

    Though the particle sizes chosen for the study retain the same aspect ratio, it is believed the surface energy has greater effects for smaller particles [63], where the same surface energy for three particle-size models are adopted, i.e., g = 0.66 N/m [58]. That is, for a particle with a smaller size, the newly created surface area has a stronger effect since crack propagation in a smaller particle acts to relieve more internal stresses. As a result, less strain energy release rates and stress-intensity factors are found in smaller particles (Figures 4 and 5a). In addition, a convergence test was performed in our finite element analysis. It has been observed that models beyond L/d = 0.6 provide similar strain energy release rates and stress intensity factors, suggesting that L/d = 0.6 is the critical value for both computational simulations and battery materials design.

    In general, an empirical Paris’ Law is used to describe crack propagation in engineering materials. However, fatigue experiments are not currently available for nanosized LiFePO4 particles. The current study specifically aims to provide a method that estimates crack advancement based on the effects of particle and crack sizes under mixed-mode loading. The current study does not aim to estimate material parameters from the Paris’ Law for LiFePO4 particles. Our approach is in contrast to that of Deshpande et al. [64], who utilized Paris’ Law and incorporated material parameters for LiFePO4 to predict battery life (i.e., numbers of cycles). In their study, an isotropic core-shell model incorporating chemical and mechanical degradation was implemented, considering a surface crack that grows with each charge-discharge cycle [64].

    The rate of discharge will determine the resultant stresses within the material, which will in turn alter the energy release rates and stress intensity factors at a crack in the particle. Therefore a current-rate-dependent crack propagation analysis for LiFePO4 is currently being developed in our research group. In general, larger stresses over a shorter period of time may induce cracking more readily than for slower discharging rates; this is due to the limited volumetric expansion rate being unable to satisfying the rate of lithium diffusivity during fast discharging, similar to rush hour traffic in limited highway space. Zhu et al. [20] have reported current density dependent fracture mechanic characteristics for another battery chemistry, LiMn2O4, wherein isotropic spherical and ellipsoidal particles are considered. They have concluded that crack propagation occurs at the center of the particle even during relatively low current density discharging [20].

    Molecular dynamic (MD) simulation has been a powerful tool to numerically solve the classical equations of motion for a system at a smaller-scale. It has been widely used to understand the behavior of liquid electrolytes for lithium-ion batteries to improve the performance of portable electronics devices, electric vehicles, and hybrid electric vehicles. However, such considerations represent limitations of our study, as we did not have the computational tools in the present study. We have chosen to conducta continuum mechanics-based finite element analyses based available experimental and computational results from the literatures in which particles and models are on the same size order of ours. In this case, we could use the available surface energy values from the literature to calculate the crack extension, as shown in Table 1. Moreover, since our current study mainly focuses on interface-diffusion limited lithium-ion battery materials, an improved model by using MD could be developed to better predict interface diffusion and failure. Another limitation in the current study lies on the small element size at the crack tip in which a MD simulation is regarded as a better choice. To this matter, the logical approach in the desire to create truly multiple scale simulations that exist at disparate length and timescales has been to couple MD and FE. Unfortunately, the coupling of these methods is not straightforward since the major problem in multi-scale simulations is that of pathological wave reflection, which occurs at the interface between the MD and FE regions. We did not have the coupled MD and FE computational tools in the present study. Therefore, the idea of meshing the finite element region down to a small-scale was one of our first attempts to eliminate spurious wave reflections, which has been generally reported at the MD-FEM interface. To this end, we are in the process by adopting MD simulations to improve our computational technique.

    4. Conclusions

    To the best of the authors’ knowledge, this is the first study utilizing the finite element method-based VCCT to obtain cracking information in LiFePO4 battery materials. Various lengths of interfacial cracks in orthotropic bi-materials are incorporated. A method is established to calculate the amount of crack propagation for lithium-ion battery materials (Figure 3 and Figure 5b). It is observed that energy release rates, stress intensity factors, and crack propagation are all particle- and crack-size-dependent (Figures 4-5). Therefore, one potential future consideration is to reduce the LiFePO4 electrode particle size to allow initial flaws during material processing (Figure 4) and to effectively decrease stress inside the particles (Figure 5a) in addition to reducing the diffusion path for better battery performance [42,43].

    The current study provides an approach towards improving battery structural design by providing a better understanding of the micromechanics, as well as a basis for future fatigue analyses, which may incorporate a variety of other battery chemistries and additional fatigue life parameters. Thus, the current study has identified factors that play an important role in inducing fracture, and the study also provides insight into factors that can be optimized to minimize any detrimental capacity loss in a battery. This research will hopefully elucidate a relationship between micro-mechanics and battery usage so as to help in the design of a higher performance, longer lasting battery for the future.

    Conflict of Interest

    All authors declare no conflict of interest in this paper.



    [1] A. K. Jain, Y. Zhong, S. Lakshmanan, Object matching using deformable templates, IEEE Trans. Pattern. Anal. Mach. Intell., 18 (1996), 267–278. https://doi.org/10.1109/34.485555 doi: 10.1109/34.485555
    [2] D. P. Kingma, M, Welling, Auto-Encoding variational bayes, preprint, arXiv: 1312/6114.
    [3] A. Ahmadi, I. Patras, Unsupervised convolutional neural networks for motion estimation, in IEEE International Conference on Image Processing, IEEE, (2016), 1629–1633. https://doi.org/10.1109/ICIP.2016.7532634
    [4] S. Joshi, B. Davis, M. Jomier, G. Gerig, Unbiased diffeomorphic atlas construction for computational anatomy, NeuroImage, 23 (2004), S151–S160. https://doi.org/10.1016/j.neuroimage.2004.07.068
    [5] M. A. Viergevera, J. B. Antoine-Maintzb, S. Kleinc, K. Murphyd, M. Staringe, J. P. W. Pluimfa, A survey of medical image registration-under review, Med. Image Anal., 33 (2016), 140–144. https://doi.org/10.1016/j.media.2016.06.030 doi: 10.1016/j.media.2016.06.030
    [6] Y. Fu, N. M. Brown, S. U. Saeed, A. Casamitjana, Y. Hu, DeepReg: a deep learning toolkit for medical image registration, J Open Source Softw., 5 (2020), 2705. https://doi.org/10.21105/joss.02705 doi: 10.21105/joss.02705
    [7] D. FAIM, A ConvNet Method for Unsupervised 3D Medical Image Registration, in International Workshop on Machine Learning in Medical Imaging, Springer, (2019), 646–654. https://doi.org/10.1007/978-3-030-32692-0_74
    [8] Y. Fu, Y. Lei, T. Wang, W. J. Curran, T. Liu, X. Yang, Deep learning in medical image registration: A review, Phys. Med. Biol., 65 (2019), 20TR01. https://doi.org/10.1088/1361-6560/ab843e doi: 10.1088/1361-6560/ab843e
    [9] Y. Fu, S. Liu, H. H. Li, D. Yang, Automatic and hierarchical segmentation of the human skeleton in CT images, Phys. Med. Biol., 62 (2017), 2812–2833. https://doi.org/10.1088/1361-6560/aa6055 doi: 10.1088/1361-6560/aa6055
    [10] X. Yang, N. Wu, G. Cheng, Z. Zhou, D. S.Yu, J. J. Beitler, et al., Automated segmentation of the parotid gland based on atlas registration and machine learning: a longitudinal MRI study in head-and-neck radiation therapy, Int. J. Radiat. Oncol. Biol. Phys., 90 (2014), 1225–1233. https://doi.org/10.1016/j.ijrobp.2014.08.350 doi: 10.1016/j.ijrobp.2014.08.350
    [11] X. Huang, Y. Zhang, J. Wang, A biomechanical modeling guided simultaneous motion estimation and image reconstruction technique (SMEIR-Bio) for 4D-CBCT reconstruction, Phys. Med. Biol., 63 (2018), 045002. https://doi.org/10.1088/1361-6560/aaa730 doi: 10.1088/1361-6560/aaa730
    [12] J. Boda-Heggemann, F. Lohr, F. Wenz, M. Flentje, M. Guckenberger, Cone-Beam CT-Based IGRT, Strahlenther Onkol., 187 (2011), 284–291. https://doi.org/10.1007/s00066-011-2236-4
    [13] X. Zhen, X. Gu, H. Yan, L. Zhou, X. Jia, S. B. Jiang, CT to Cone-beam CT deformable registration with simultaneous intensity correction, Phys. Med. Biol., 57 (2012), 6807–6826. https://doi.org/10.1088/0031-9155/57/21/6807 doi: 10.1088/0031-9155/57/21/6807
    [14] C. Veiga, J. Mcclelland, S. Moinuddin, A. Lourenço, K. Ricketts, J. Annkah, et al., Toward adaptive radiotherapy for head and neck patients: Feasibility study on using CT-to-CBCT deformable registration for "dose of the day" calculations, Med. Phys., 41 (2014), 031703. https://doi.org/10.1118/1.4864240 doi: 10.1118/1.4864240
    [15] C. Veiga, G. Janssens, C. L. Teng, T. Baudier, B. K. Teo, First clinical investigation of cone beam computed tomography and deformable registration for adaptive proton therapy for lung cancer, Int. J. Radiat. Oncol. Biol. Phys., 95 (2016), 549–559. https://doi.org/10.1016/j.ijrobp.2016.01.055 doi: 10.1016/j.ijrobp.2016.01.055
    [16] B. Zhou, Z. Augenfeld, J. Chapiro, S. K. Zhou, C. Liu, J. S. Duncan, Anatomy-guided multimodal registration by learning segmentation without ground truth: Application to intraprocedural CBCT/MR liver segmentation and registration, Med Image Anal, 71 (2021), 102041. https://doi.org/10.1016/j.media.2021.102041 doi: 10.1016/j.media.2021.102041
    [17] B. Zhou, C. Liu, J. S. Duncan, Anatomy-Constrained contrastive learning for synthetic segmentation without ground-truth, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, (2021), 47–56. https://doi.org/10.1007/978-3-030-87193-2_5
    [18] X. Liang, Y. Jiang, Y. Xie, T. Niu, Quantitative cone-beam CT imaging in radiotherapy: Parallel computation and comprehensive evaluation on the TrueBeam system, IEEE Access, 7 (2019), 66226–66233. https://doi.org/10.1109/ACCESS.2019.2902168 doi: 10.1109/ACCESS.2019.2902168
    [19] K. Srinivasan, M. Mohammadi, J. Shepherd, Investigation of effect of reconstruction filters on cone-beam computed tomography image quality, Australas. Phys. Eng. Sci. Med., 37 (2014), 607–614. https://doi.org/10.1007/s13246-014-0291-8 doi: 10.1007/s13246-014-0291-8
    [20] L. Ouyang, T. Solberg, J. Wang, Noise reduction in low‐dose cone beam CT by incorporating prior volumetric image information, Med. Phys., 39 (2012), 2569–2577. https://doi.org/10.1118/1.3702592 doi: 10.1118/1.3702592
    [21] W. Kang, M. Patwari, Low Dose Helical CBCT denoising by using domain filtering with deep reinforcement learning, preprint, arXiv: 2104/00889.
    [22] J. Zheng, D. Zhang, K. Huang, Y. Sun, Cone-Beam computed tomography image pretreatment and segmentation, in 2018 11th International Symposium on Computational Intelligence and Design (ISCID), IEEE, (2018), 25–28. https://doi.org/10.1109/ISCID.2018.00012
    [23] X. Liang, L. Chen, D. Nguyen, Z. Zhou, X. Gu, et al., Generating synthesized computed tomography (CT) from cone-beam computed tomography (CBCT) using CycleGAN for adaptive radiation therapy, Phys. Med. Biol., 64 (2019), 125002. https://doi.org/10.1088/1361-6560/ab22f9 doi: 10.1088/1361-6560/ab22f9
    [24] D. C. Hansen, G. Landry, F. Kamp, M. Li, C. Belka, K. Parodi, et al., ScatterNet: A convolutional neural network for cone‐beam CT intensity correction, Med. Phys., 45 (2019), 4916–4926. https://doi.org/10.1002/mp.13175 doi: 10.1002/mp.13175
    [25] G. Wu, M. Kim, Q. Wang, Y. Gao, D. Shen, Unsupervised deep feature learning for deformable registration of MR brain images, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, (2013), 649–656. https://doi.org/10.1007/978-3-642-40763-5_80
    [26] M. Simonovsky, B. Gutiérrez-Becker, D. Mateus, N. Navab, N. Komodakis, A deep metric for multimodal registration, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, (2016). https://doi.org/10.1007/978-3-319-46726-9_2
    [27] C. Hao, D. Qi, N. Dong, J. Z. Cheng, P. A. Heng, Automatic fetal ultrasound standard plane detection using knowledge transferred recurrent neural networks, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, (2015). https://doi.org/10.1007/978-3-319-24553-9_62
    [28] S. Miao, Z. J. Wang, L. Rui, A CNN regression approach for Real-Time 2D/3D registration, IEEE Trans. Med. Imaging, 35 (2016), 1352–1363. https://doi.org/10.1109/TMI.2016.2521800 doi: 10.1109/TMI.2016.2521800
    [29] X. Yang, R. Kwitt, M. Niethammer, Fast rredictive image registration, in Deep Learning and Data Labeling for Medical Applications, Springer, (2016), 48–57. https://doi.org/10.1007/978-3-319-46976-8_6
    [30] D. Detone, T. Malisiewicz, A. Rabinovich, Deep image homography estimation, preprint, arXiv: 1606/03798.
    [31] T. Carvalho, E. Rezende, M. Alves, F. Balieiro, R. B. Sovat, Exposing computer generated images by eye's region classification via transfer learning of VGG19 CNN, in 16th IEEE International Conference On Machine Learning And Applications, IEEE, (2017). https://doi.org/10.1109/ICMLA.2017.00-47
    [32] J. Krebs, T. Mansi, H. DelingetteZ. Li, F. Ghesu, S. Miao, et al., Robust non-rigid registration through agent-based action learning, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, (2017), 344–352. https://doi.org/10.1007/978-3-319-66182-7_40
    [33] J. Zhang, C. Wang, S. Liu, L. Jia, J. Sun, Content-Aware unsupervised deep homography estimation, in European Conference on Computer Vision, Springer, (2020), 653–669. https://doi.org/10.1007/978-3-030-58452-8_38
    [34] M. Bevilacqua, A. Roumy, C. Guillemot, M. L. A. Morel, Super-resolution using neighbor embedding of back-projection residuals, in 2013 18th International Conference on Digital Signal Processing (DSP), 2013.07, IEEE, (2013), 1–8. https://doi.org/10.1109/ICDSP.2013.6622796
    [35] M. Bevilacqua, A. Roumy, C. Guillemot, M. L. A. Morel, Low-Complexity single-image super-resolution based on nonnegative neighbor embedding, in Proceedings of the British Machine Vision Conference 2012, (2012), 1–10. https://doi.org/10.5244/C.26.135
    [36] R. Timofte, V. De, L. V. Gool, Anchored Neighborhood Regression for Fast Example-Based Super-Resolution, in IEEE International Conference on Computer Vision, IEEE, (2014), 1920–1927. https://doi.org/10.1109/ICCV.2013.241
    [37] J. Kim, J. K. Lee, K. M. Lee, Accurate image super-resolution using very deep convolutional networks, in 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), IEEE, (2016), 1646–1654. https://doi.org/10.1109/Cvpr.2016.182
    [38] K. Simonyan, A. Zisserman, Very deep convolutional networks for large-scale image recognition, preprint, arXiv: 1409/1556.
    [39] C. Dong, C. C. Loy, K. He, X. Tang, Image super-resolution using deep convolutional networks, IEEE Trans. Pattern. Anal. Mach. Intell., 38 (2016), 295–307. https://doi.org/10.1109/TPAMI.2015.2439281 doi: 10.1109/TPAMI.2015.2439281
    [40] J. Chen, Y. Li, Y. Du, E. C. Frey, Generating anthropomorphic phantoms using fully unsupervised deformable image registration with convolutional neural networks, Med. Phys., 47 (2020). https://doi.org/10.1002/mp.14545
    [41] R. Sandkühler, C. Jud, S. Andermatt, P. C. Cattin, AirLab: Autograd image registration laboratory, preprint, arXiv: 1806/09907.
    [42] O. Ronneberger, P. Fischer, T. Brox, U-Net: Convolutional networks for biomedical image segmentation, in International Conference on Medical Image Computing and Computer-Assisted Intervention, Springer, (2015), 234–241. https://doi.org/10.1007/978-3-319-24574-4_28
    [43] D. Rueckert, L. I. Sonoda, C. Hayes, D. Hill, M. O. Leach, Nonrigid registration using free-form deformations: application to breast MR images, IEEE Trans. Med. Imaging, 18 (2012), 712–721. https://doi.org/10.1109/42.796284 doi: 10.1109/42.796284
    [44] T. Fechter, D. Baltas, One shot learning for deformable medical image registration and periodic motion tracking, IEEE Trans. Med. Imaging, 39 (2019), 12. https://doi.org/10.1109/TMI.2020.2972616 doi: 10.1109/TMI.2020.2972616
    [45] A. V. Dalca, M. Rakic, J. Guttag, M. R. Sabuncu, Learning conditional deformable templates with convolutional networks, in NIPS'19: Proceedings of the 33rd International Conference on Neural Information Processing Systems, Curran Associates Inc., (2019), 806–818.
    [46] P. Viola, W. Iii, Alignment by maximization of mutual information, Int. J. Comput. Vision, 24 (2008), 137v154. https://doi.org/10.1023/A:1007958904918 doi: 10.1023/A:1007958904918
    [47] Z. S. Saad, D. R. Glen, G. Chen, M. S. Beauchamp, R. Desai, R. W. Cox, A new method for improving functional-to-structural MRI alignment using local Pearson correlation, Neuroimage, 44 (2009), 839–848. https://doi.org/10.1016/j.neuroimage.2008.09.037 doi: 10.1016/j.neuroimage.2008.09.037
    [48] Z. Shen, X. Han, Z. Xu, M. Niethammer, Networks for joint affine and Non-parametric image registration, in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition, IEEE, (2019), 4224–4233. https://doi.org/10.1109/CVPR.2019.00435
    [49] M. Staring, S. Klein, J. Pluim, A rigidity penalty term for nonrigid registration, Med. Phys., 34 (2007), 4098–4108. https://doi.org/10.1118/1.2776236 doi: 10.1118/1.2776236
    [50] L. Hui, P. Du, W. Zhao, L. Zhang, H. Sun, Image registration based on corner detection and affine transformation, in International Congress on Image & Signal Processing, IEEE, (2010), 2184–2188. https://doi.org/10.1109/CISP.2010.5647722
    [51] C. Papazov, D. Burschka, Deformable 3D shape registration based on local similarity transforms, in Computer Graphics Forum, Wiley Online Library, (2011), 1493–1502. https://doi.org/10.1111/j.1467-8659.2011.02023.x
  • This article has been cited by:

    1. Sangwook Kim, Hsiao-Ying Shadow Huang, Mechanical stresses at the cathode–electrolyte interface in lithium-ion batteries, 2016, 31, 0884-2914, 3506, 10.1557/jmr.2016.373
    2. Susmita Sarkar, Rajankumar L. Patel, Xinhua Liang, Jonghyun Park, Unveiling the Role of CeO2 Atomic Layer Deposition Coatings on LiMn2O4 Cathode Materials: An Experimental and Theoretical Study, 2017, 9, 1944-8244, 30599, 10.1021/acsami.7b06988
    3. S. Esmizadeh, H. Haftbaradaran, F. Mossaiby, An investigation of the critical conditions leading to deintercalation induced fracture in two-phase elastic electrode particles using a moving interphase core-shell model, 2019, 74, 09977538, 96, 10.1016/j.euromechsol.2018.10.019
    4. H. Haftbaradaran, A. Maddahian, F. Mossaiby, A fracture mechanics study of the phase separating planar electrodes: Phase field modeling and analytical results, 2017, 350, 03787753, 127, 10.1016/j.jpowsour.2017.03.073
    5. Netanel Shpigel, Mikhael D. Levi, Xiaopeng Cheng, Tianci Cao, Rui Wu, Tyler S. Mathis, Yuefei Zhang, Doron Aurbach, Yury Gogotsi, Diffusion-Induced Transient Stresses in Li-Battery Electrodes Imaged by Electrochemical Quartz Crystal Microbalance with Dissipation Monitoring and Environmental Scanning Electron Microscopy, 2019, 4, 2380-8195, 1907, 10.1021/acsenergylett.9b00884
    6. Nan Luo, Yong Lin, Jian Guo, Emanuele Quattrocchi, Huaijiu Deng, Jian Dong, Francesco Ciucci, Filippo Boi, Chunfeng Hu, Salvatore Grasso, Spark Plasma Sintering of LiFePO4: AC Field Suppressing Lithium Migration, 2021, 14, 1996-1944, 2826, 10.3390/ma14112826
    7. D Clerici, F Mocera, F Pistorio, Analysis of fracture behaviour in active materials for lithium ion batteries, 2022, 1214, 1757-8981, 012018, 10.1088/1757-899X/1214/1/012018
    8. M. Nikpour, B. A. Mazzeo, D. R. Wheeler, A Model for Investigating Sources of Li-Ion Battery Electrode Heterogeneity: Part II. Active Material Size, Shape, Orientation, and Stiffness, 2021, 168, 0013-4651, 120518, 10.1149/1945-7111/ac3c1f
    9. Iheke Micheal Nwachukwu, Assumpta Chinwe Nwanya, A.B.C. Ekwealor, Fabian I. Ezema, Recent progress in Mn and Fe-rich cathode materials used in Li-ion batteries, 2022, 54, 2352152X, 105248, 10.1016/j.est.2022.105248
    10. Chihyun Nam, Bonho Koo, Juwon Kim, Jinkyu Chung, Jaejung Song, Danwon Lee, Sungjae Seo, Munsoo Song, Seyeon Shin, Namdong Kim, Markus Weigand, Jian Wang, Jongwoo Lim, Dynamic Lithium Transport Pathway via Crack Formation in Phase-Separating Battery Particles, 2025, 1936-0851, 10.1021/acsnano.4c15960
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4028) PDF downloads(240) Cited by(1)

Article outline

Figures and Tables

Figures(6)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog