Abstract
The thermal conductivities and elastic properties of carbon nanotubes (CNTs) are estimated by using the double-inclusion model, which is based on rigorous elasticity approach. The model regards a CNT as one inclusion (the inner cylindrical void) embedded in the other (the outer coaxial single-crystal graphite shell). The concept of homogenization is employed, and vital microstructural parameters, such as CNT diameter, length, and aspect ratio, are included in the present model. The relationship between microstructure and thermal conductivities and elastic stiffness of CNTs is quantitatively characterized. Our analytical results, benchmarked by experimental data, show that the thermal conductivities and elastic stiffness of CNTs are strongly dependent on the diameter of CNT with little dependence on the length of CNT.
1. Introduction
The present work employs and modifies the double-inclusion model [1] to estimate, from a microstructural point of view, the thermal conductivities and elastic stiffness of carbon nanotubes (CNTs). Several models based on molecular dynamics simulation are developed for estimating the thermal conductive and elastic properties of CNTs [2–7] yet are not very computationally efficient. Some of the existing models are based on beam theories [8, 9] and can only estimate the Young moduli of CNTs. Robertson et al. [10] and Lu [11] employ energy methods to estimate the elastic properties of CNTs. Finite deformation is considered in Gao and Li [12] where the formulation is aimed at single-walled CNTs. A multiscale shear-lag-based approach later is developed for estimating Young’s moduli of CNTs as well as stresses in CNT-based composites [13]. Chantrenne and Barrat [14] derive analytical expressions for the thermal conductivities of a single graphene and a CNT as a function of their characteristic lengths. Their analytical estimates show good agreement for single graphene but not for CNTs.
The reported experimental data of thermal conductivities and Young’s and shear moduli of CNTs from literature are summarized in Tables 1 and 2, respectively. It is noted that the reported data are significantly different from each other as the thermal and elastic properties of CNTs are functions of the microstructure of CNTs [15–18]. Qian et al. [19] also collect many experimental data on the Young modulus of CNTs, and they range from 320 GPa to 5.0 TPa. It is the objective of the present work to develop an analytical model to delineate the relationship between the microstructure and thermal conductive and elastic stiffness of CNTs.
2. The Model
2.1. Mori-Tanaka Theorem
Consider an infinite medium , embodying a subregion , which includes an inclusion undergoing a heat flux-free transformation temperature gradient . The matrix is denoted by , and , see Figure 1. In view of the concept of Eshelby [20], Mori and Tanaka [21], and Hatta and Taya [22], one may find the disturbance temperature gradient at a typical location, , in due to a transformation temperature gradient in as follows: where the case of ellipsoidal and are considered. Eshelby’s matrix [] is a matrix and is a function of the geometry of inclusion . In particular, the components of [] for a circular cylindrical are only functions of the aspect ratio of the cylinder [22]. It can be shown in a straightforward manner that the average temperature gradient of , , is given by where the angle brackets represent an average quantity, is the average temperature gradient field of the inclusion, is the average temperature gradient field of the matrix, and is the inclusion volume fraction. The average disturbance temperature gradient in region due to is given by where is Eshelby’s matrix of inclusion . Rearranging (3) gives Substituting (1) and (4) into (2) gives the disturbance temperature gradient averaged over due to a uniform transformation temperature gradient prescribed in ,
Now consider a uniform transformation temperature gradient in . The heat flux and temperature gradient fields can be obtained by superposing two fields—one is due to and distributed over , and the other is due to and distributed over . The average disturbance temperature gradient in due to in is given by The average disturbance temperature gradient in due to prescribed in is given by The first term on the right-hand side of (7) is in essence identical to (3) except being replaced by . Thus, (4) can be used to give Substituting (8) and (6) into (2), one obtains the disturbance temperature gradient averaged over due to transformation temperature gradient prescribed in as follows:
2.2. Extended Mori-Tanaka Theorem
Consider an infinite homogeneous solid with a thermal conductivity . The solid contains two ellipsoidal inclusions, and . Let and () undergo distinct heat flux-free transformation. The uniform transformation temperature gradient in the absence of the surrounding medium is and in and , respectively. The average of the disturbance temperature gradient resulting from is given by (1) and (5), and the volume average of the corresponding heat flux field is given bywhere is the identity matrix.
Similarly, the temperature gradient resulted from is given by (6) and (9), and the corresponding heat flux fields are given by Combining (1), (6), (10a), and (11a) gives the average temperature gradient and heat flux in due to and ,The combination of (5), (9), (10b) and (11b) gives the average temperature gradient and heat flux in due to and ,The following relations can be derived in a straightforward manner:Substituting (12a) and (12b), (13a) and (13b) into (14a) and (14b) gives
2.3. Thermal Conductivities and Elastic Constants of Double Inclusion
Consider an ellipsoidal matrix , containing an ellipsoidal inhomogeneity . The matrix is embedded in an infinite medium . The thermal conductivity of , , and is given by respectively; see Figure 2. To solve the double-inhomogeneity problem, one may replace the inhomogeneity, , and the matrix, , by a reference material with thermal conductivity , which is the thermal conductivity of . Then one could prescribe transformation temperature gradient and in and , respectively, to compensate the material mismatch, which is essentially Eshelby’s concept of equivalent inclusion [20]. When the infinity domain is subjected to far-field temperature gradient , the heat flux and temperature gradient fields of the double inclusion () are in general not uniform, and their averages are, with the help of (12a), (12b), and (15a) and (15b), given bySince the temperature gradient and the heat flux fields must be preserved after homogenization, the following constraint conditions must be satisfied:Solving (19a) and (19b) simultaneously, one findswhere Substituting (17a) and (17b) and (18a) and (18b) into (14a) and (14b), one findsThe thermal conductivities of the double inclusion are defined by The combination of (20a) and (20b), (23a) and (23b), and (24) gives A reasonable choice for the reference material for homogenization is the matrix material. Thus, . Letting in (25) gives In the case of CNT, the inner inclusion is a void, and the outer M is a graphene. Therefore, the [] in (26) may be neglected as the thermal conductivity of air is much less than that of the graphene, or in the case of a vacuum. Therefore, the thermal conductivities of a CNT can be estimated by Specifically,
In view of the analogy between heat conduction and linear elasticity, the elastic constants , which is a matrix, of a CNT can be estimated by where [] is the stiffness matrix of a single-crystal graphite, [] is the identity matrix, and the components of Eshelby’s matrix for anisotropic elastic inclusions are given by Mura [23]. In particular, the shear moduli of a CNT can be estimated by
3. Results
The following elastic properties of a single-crystal graphite are used in our calculations [24]: As it is indicated in Figure 3, where the diameter and thickness of CNT are assumed to be 1.2 nm and 0.34 nm, respectively, the length of CNT does not affect its axial and transverse Young’s moduli significantly. The relationship between axial Young’s modulus and the diameter of CNT is shown in Figure 4, where the length and thickness of CNT are assumed to be 1 μm and 0.34 nm, respectively. The predicted axial Young’s modulus increases with a decrease in CNT diameter, which is consistent with the available experimental data, yet is in contrast with what is predicted by Li and Chou [18]. It is noted that axial Young’s modulus drops significantly when the CNT diameter is greater than 1.2 nm. The transverse Young modulus of CNT exhibits a similar relation with CNT diameter. The relationship between CNT diameter and thickness, for a fixed axial Young’s modulus of 1 TPa, is shown in Figure 5, which agrees with Thostenson et al. [25].
The in-plane and out-of-plane thermal conductivities of graphene are assumed to be equal to 2,000 and 10 W/m-K, respectively [26] in the present calculations. The thermal conductivity of air varies from 0.02 to 0.08 W/m-K as the temperature varies from −55°C to 1,000°C and is negligible. As it is indicated in Figure 6, where the diameter of CNT is assumed to be 1.2 nm and 1.4 nm and the CNT thickness is assumed to be 0.34 nm, the length of CNT does not affect its axial thermal conductivity significantly, which is consistent to the findings of Yang et al. [27] and Hone et al. [28]. The relationship between CNT diameter and its thermal conductivity for a fixed length of 100 nm is shown in Figure 7, where the thickness of CNT is assumed to be 0.34 nm, and the axial thermal conductivity of CNT decreases dramatically with an increase in CNT diameter when the CNT diameter is less than 10 nm. When the CNT diameter is greater than 10 nm, the decrease in axial thermal conductivity of CNT becomes gradual. Figure 8 shows the relationship between CNT aspect ratio and its axial thermal conductivity. The axial thermal conductivity increases with an increase in the aspect ratio of CNT [17]. The relationship between the axial thermal conductivity of MWCNT and its number of layers is shown in Figure 9, where the MWCNT outer diameter, length, and thickness of each layer are assumed to be 14 nm, 100 nm, and 0.34 nm, respectively.
4. Conclusion
In the present work, the double-inclusion model [1], which is based on elasticity solutions [20, 21] and is originally presented in the context of linear elasticity, is employed and extended to predict the thermal conductivities and elastic constants of CNTs. Though the theorem proposed by Mori and Tanaka [21] is not restricted to ellipsoidal inclusions, the present work regards a CNT as a circular cylindrical shell containing a concentric circular cylindrical void, that is, two coaxial circular cylindrical inclusions—one is embedded in the other. The outer shell is composed of a single-crystal graphite sheet.
Our results indicate that (1) the length of CNT does not significantly affect its thermal conductivities, axial and transverse Young’s and shear moduli; (2) with an increase in CNT diameter, axial and transverse Young’s moduli and thermal conductivities as well as shear moduli decrease. This may explain why the reported experimental data in the literature (see Tables 1 and 2) are scattered when the thermal conductivity and elastic properties of a CNT are considered as functions of its microstructure, in particular, its diameter. The results derived from the present model agree quite well with the available experimental measurement and observation.
Acknowledgment
This work is supported by the National Science Council in Taiwan under Grant no. NSC95-2221-E-155-016 to Yuan Ze University, Taiwan.