Issue 
EPJ Appl. Metamat.
Volume 8, 2021



Article Number  16  
Number of page(s)  7  
DOI  https://doi.org/10.1051/epjam/2021010  
Published online  23 August 2021 
https://doi.org/10.1051/epjam/2021010
Research Article
Estimating the effective conductivity for ellipseinclusion model with Kapitza thermal resistance
Faculty of Mechanical Engineering, Hanoi University of Industry, 298 Cau Dien Street, Bac Tu Liem District, Hanoi, Vietnam
^{*} email: luatnv1980@gmail.com
Received:
4
March
2021
Accepted:
14
July
2021
Published online: 23 August 2021
The ellipse assemblage model with imperfect interface has quite complex microstructure, that can be considered an extension of the circle assemblage model with imperfect interfaces. The paper introduces an approximate method for computing the effective conductivity of isotropic composites with imperfect interfaces in twodimensional space. Based on the coatedellipse assemblage model and the equivalent inclusion approximation, one can determine the effective thermal conductivity of the composites. The polarization approximation is given in an explicit form (PEK) and this method will be applied to calculate the effective conductivity of the composite with Kapitza thermal resistance model. The PEK result will have compared with the Fast Fourier Transform (FFT) simulation and Hashinstrikman bounds (HS).
Key words: Effective conductivity / homogenization / imperfect interface / elliptical Kapitza thermal resistance
© V.L. Nguyen, Published by EDP Sciences, 2021
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Due to limited available information about the materials as well as the complexity of microstructure, the determination of macroscopic moduli of composite materials is not simple. Especially composite materials with imperfect contact surfaces, this model needs to have its own approach solution, that will be presented in this work. Determination of macroscopic properties of composite materials from the properties of the component materials and their microstructure (homogenization) is the subject of many studies in recent decades. There are only a few cases where macroscopic properties are determined exactly. The solution of the problem of homogenization can be divided into the following directions: The first one is to construct the upper and lower bounds on the effective properties based on variational principles [1–4], the second one is to develop the approximate methods using Eshelby's solution [5] such as [6–10]. Besides, with the fast development of computer science, the numerical methods [10–15] have been applied in mechanics of materials and effective in solving complex problems. The objective of this research part is to present the methods (approximateanalytic and numerical simulation), which determine macroscopic properties of composite materials with complex microstructures. With complex inclusion such as coatedinclusions or those with imperfect interfaces, the fields inside and around the inclusion are complex, making more difficult the development of approximate models. Benveniste et al. [16] has developed two models of thin interface in the problem of thermal conductivity and elasticity, the model allows determining the fields of adjacent phases without these fields of interphase itself. Le Quang et al. [17] have been derived the closedform expressions for the effective conductivity of anisotropic multiphase composites with imperfect interfaces by using the Eshelby's formalism and dilute solution results. Yvonnet et al. [11] have used XFEM method to the thermal conductivity of the material with imperfect interfaces. Monchiet [15] has developed algorithms based on XFFT to simulate discontinuous phenomenon of temperature between the phases (phenomenon Kapitza). Nguyen [12] has developed a simple approach to approximate macroscopic properties of the material in this case by replacing the coatedellipse inclusion by a homogeneous equivalent inclusion. Many authors have mentioned the phenomenon of imperfect interfaces but most applications of these studies are based on the solution of circleinclusion models. This work focuses on the conductivity of ellipse assemblage model with imperfect interfaces, in which the solution for this problem is based on polarization approximation (PEK) with elliptical Kapitza thermal resistance.
2 The ellipse inclusion with interfacial thermal resistance
We consider the twocomponent coated ellipsoid assemblage with thin coating (volume fractions v_{2} ≪ 1, v_{1} + v_{2} = 1), in which the ellipse made material1 is the inclusion the ellipse shells made material2 is the matrix. In the imperfect interfaces model with the infinitely thin coating thickness h, . Because elliptical inclusions of phase 1 coated by the infinitely thin shells, we can consider a uniform thickness interphase.
We have (1)where a_{1}, a_{2}, a_{3} are the principal radii of the ellipsoid in the direction of the axes x_{1}, x_{2}, x_{3}.
In the case of two component matrix composites with inclusions of ellipsoidal forms, in which v_{I}, c_{I} are the volume fraction and conductivity of the inclusion; v_{M}, c_{M} are the volume fraction and conductivity of the matrix. The expression of the effective conductivity for the dilute suspension of randomly oriented ellipsoids (hence the effective medium is isotropic) can be presented as (v_{I} ≪ 1, v_{I} + v_{M} = 1) [5](2)where D is the respective inclusionfunction. In the case of twodimensional space (d = 2) with the same aspect ratio r = a_{2}/a_{1} in a continuous matrix, one has(3)
On the other side, we used the polarization approximation for twocomponent, the effective conductivity C^{eff} has form(4)where the reference parameter c_{*} should be determined from a reference dilute solution result, or reference effective conductivity of the composite at certain finite volume proportions of the components. Note that, function P_{C}(c_{*}) is a positive monotonousincreasing function of the positive parameter c_{*}.
The polarization approximation P_{C}(c_{*}) would obey HashinStrikman bounds [1]:(5)the respective reference parameter c_{*} should lie within the limits(6)
The polarization approximation (4) using dilute solution reference (2) or equalizing (4) and (2), one finds the reference parameter c_{*} in the twodimensional space (d = 2)(7)
Substituting (7) into (4) after some manipulation, one obtains(8)
The application (8) for coatedellipse inclusions, consider the ellipse made of material1 are embedded in the ellipse shell of material2, in which the ellipse made material1 is the inclusion with conductivity c_{I} = c_{1}, the ellipse shell made material2 is the matrix with conductivity c_{M} = c_{2}.
In the context of thermal conduction with the imperfect interface models, the most widely used ones are the Kapitza's interface thermal resistance model, called also lowly conducting (LC) interface model, and the highly conducting (HC) interface model. In this work, we focus on the lowly conducting imperfect interface model [17,18], it is assumed that(9)where α is called the Kapitza thermal resistance.
Based on (9), (8) and (3), one derives the asymptotic expression of the effective conductivity c_{1L} of the assemblage of ellipse inclusions of phase1 coated by the infinitely thin shells (phase2) with thermal resistance α (Fig. 1)(10)
Substituting (1) and (9) into (10) (in the limit h → 0, v_{1} = 1) we come to(11)
Especially, when r = 1, a_{1} = a_{2} = R, d = 2, we derive the asymptotic expression of the effective conductivity c_{1L} of the assemblage of circular inclusions, that has been presented in [19](12)
For the coated ellipse assemblage model, if these coated inclusions are embedded in the elliptic shells of the matrix phase M, one obtains the effective conductivity of the elliptic assemblage with interfacial thermal resistance α (in the limit h → 0) between the inclusion 1 and matrix M components by replacing c_{EI} for c_{I} in (4), (7) and (8) (Fig. 2). The twophase coated inclusions with imperfect interface are substituted by equivalent onephase inclusions c_{EI} = c_{1L} (12), from (7) we have(13)
The polarization approximation has the form (PEK)(14)
The polarization approximation (14) contains all the components of the coatedellipse assemblage model with imperfect interface such as c_{1}, c_{M}, α, a_{1}, a_{2}. So, one can apply the simple polarization approximation (14) to obtain an estimate for the effective conductivity of the original elliptic inclusion with imperfect interface.
Fig. 1
A single Coated ellipsoid. 
Fig. 2
Unit cell: Coated ellipse (left), equivalent inclusion (right). 
3 FFT simulation for composite materials with elliptical Kapitza thermal resistance
Transverse isotropic unidirectional fiberreinforced composite is well presented by polydispersed coated ellipse cylinders, whose phase boundaries are cylindrical ellipse surfaces, with generators (in x_{3}direction) orthogonal to the plane of isotropy (Fig. 1). Consider a composite reinforced by fibers aligned in the direction x_{3}. The ellipse inclusions are periodically distributed in the directions x_{1} and x_{2}. The contact between inclusion and matrix is the imperfect interface of Kapitza's type. This problem can be modelled as a composite of threecomponent: Inclusions (v_{1}, c_{1}), matrix (v_{M}, c_{M}) and interphase with volume fraction v_{2}, conductivity c_{2} = h/α (Fig. 2). By substitution scheme, one can consider the twocomponent of coatedellipse with perfect interface in which the conductivity of equivalent inclusion is determined by (11) (15)
The main content of the FFT method is based on the known equations, equilibrium conditions and the Fourier transform of the gradient E(x), flux field J(x), the second order local conductivity tensor of the heterogeneous medium C(x) to establish LippmannSchwinger integral equations for heterogeneous material problem and use of the Green periodic operator.
The algorithm for determining the effective conductivity of twophase periodic materials based on the Fourier transform method (FFT) has been introduced in the literature [12,14,19]. In the next section, this method will be applied to calculate the effective conductivity of the composite with Kapitza thermal resistance.
4 Application and comparison
For numerical FFT illustrations, consider periodic models: Square model (Fig. 3). We take c_{M} = 1, the inclusion with c_{1} > c_{M}, c_{1} = 10 or c_{1} = 100, h=0.0001, c_{2} = h/α, thermal resistance α = 0.05; α = 0.2; α = 1. Consider a unit cell having the dimension L = 1 along each space directions containing inclusion with the dimensionless, ratio r = a_{2}/a_{1} varies from 0.1 to 1, with nonoverlapping condition of ellipse inclusions, we assumed that a_{1} = 0.25 for square model (corresponding volume fraction from 0.0393 to 0.3927).
The results of calculations are reported in Tables 1–3 and Figures 4–9. From all tables and figures, we can see that: the FFT simulation and PEK approximation are coincided and always obey HS bounds (those includes upper bound denoted by HSU and lower bound denoted HSL), one also sees the dependence of the effective conductivity on the Kapitza thermal resistance. When α = 1, volume fraction of inclusion increases but effective conductivity decreases. α = 0.2 one observes a slight decrease then increase. α = 0.05 two curves increase. In all cases, when conductivity of the inclusion increases 10 times then effective conductivity increases only slightly. This shows the dominance of the thermal resistance coefficient for the macroscopic conductivity.
The FFT results in a special case r =1, a_{1} = a_{2} = R which were compared with results of Monchiet [15] for square model with one circle inclusion in the unit cell (Fig. 2), the unit cell having the dimension L = 1 along each space directions containing inclusion with the dimensionless radius R varies from 0.1 to 0.5. The conductivity of the equivalent inclusion is determined from (13). The result of comparison is reported in Table 4, those are completely close.
Fig. 3
Unit cell with square ellipse inclusion (left), equivalent inclusion (right). 
Comparison of results (C^{eff}) of FFT and PEK for square model, α = 1, a_{1} = 0.25.
Comparison of results (C^{eff}) of FFT and PEK for square model, α = 0.2, a_{1} = 0.25.
Comparison of results (C^{eff}) of FFT and PEK for square model, α = 0.05, a_{1} = 0.25.
Fig. 4
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 1, . 
Fig. 5
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 1, . 
Fig. 6
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.2, . 
Fig. 7
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.2, . 
Fig. 8
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.05, . 
Fig. 9
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.05, . 
Comparison of results (C^{eff}) of FFT, PEK and Monchiet's results for square model with one circle inclusion in the unit cell α = 1, a_{1} = a_{2} = R.
5 Conclusion
The paper has presented a solution for effective conductivity of the coatedellipse assemblage model with imperfect interfaces, the effective conductivity of isotropic multiphase composites with imperfect interfaces can be determined with two different approaches. The first approach is based on the model of coatedellipse inclusion with thin coating, the PEK expression for the effective conductivity have been derived by using the Eshelby's formalism and dilute solution results combined with polarization approximation. The second approach is based on the FFT simulation combining equivalent inclusion approach. Using different approaches, polarization approximation and FFT simulation have been proposed to estimate the effective conductivity of isotropic material composites with imperfect interfaces of Kapitza's type while still ensuring accuracy. FFT and PEK give almost same results, the estimation obtained by the FFT and PEK scheme can be viewed an excellent approximation for the effective thermal conductivity of the square models. The FFT results in a special case r = 1, a_{1} = a_{2} = R which were compared with results of Monchiet, the results of the comparison are completely close. Especially with the simple expression PEK for coated ellipse aggregate model that can be applied to good results for imperfect interface contact.
References
 Z. Hashin, S. Shtrikman, A variational approach to the theory of the elastic behaviour of multiphase materials, J. Mech. Phys. Solids 11, 127–140 (1963) [CrossRef] [MathSciNet] [Google Scholar]
 R. Hill, Theory of mechanical properties of fibrestrengthened materials: I. Elastic behavior, J. Mech. Phys. Solids 12, 199–212 (1964) [CrossRef] [MathSciNet] [Google Scholar]
 J.R. Willis, Variational and related methods for the overall properties of composite materials, in Advances in Applied Mechanics (Academic Press, New York, 1981) [Google Scholar]
 D.C. Pham, L.D. Vu, V.L. Nguyen, Bounds on the ranges of the conductive and elastic properties of randomly inhomogeneous materials, Philos. Mag. 93, 2229–2249 (2013) [CrossRef] [Google Scholar]
 J.D. Eshelby, The determination of the elastic field of an ellipsoidal inclusion, and related problems, Proc. Roy. Soc. London A 241, 376–396 (1957) [CrossRef] [MathSciNet] [Google Scholar]
 S. Torquato, Random Heterogeneous Materials (SpringerVerlag, New York, 2002) [CrossRef] [Google Scholar]
 R.A. Hill, Selfconsistent mechanics of composite materials, J. Mech. Phys. Solids 13, 213–222 (1965) [CrossRef] [Google Scholar]
 T. Mori, K. Tanaka, Average stress in the matrix and average elastic energy of materials with misfitting inclusions, Acta Metall. 21, 571–574 (1973) [CrossRef] [Google Scholar]
 T.K. Nguyen, D.C. Pham, V.L. Nguyen, Conduction in 2D and 3D dimensional sphericallysymmetric anisotropiccoating inclusion composites, Int. J. Eng. Sci. 154, 103352 (2020) [CrossRef] [Google Scholar]
 J. Michel, H. Moulinec, P. Suquet, Effective properties of composite materials with periodic microstructure: a computational approach, Comput. Methods Appl. Mech. Eng. 172, 109–143 (1999) [CrossRef] [Google Scholar]
 J. Yvonnet, Q.C. He, Q.Z. Zhu, J.F. Shao, A general and efficient computational procedure for modeling the Kapitza thermal resistance based on XFEM, Comput. Mater. Sci. 50, 1220–1224 (2011) [CrossRef] [Google Scholar]
 V.L. Nguyen, FFT and Equivalentinclusion approach for effective conductivity of threephase composites with random coatedellipse inclusion, Eng. Res. Express 3, 025014 (2021) [CrossRef] [Google Scholar]
 S. NematNasser, M. Hori, Micromechanics: Overall Properties of Heterogeneous Materials (NorthHolland, Amsterdam, 1999) [Google Scholar]
 N. Van Luat, N. Trung Kien, FFTsimulations and multicoated inclusion model for macroscopic conductivity of 2D suspensions of compound inclusions, Vietnam J. Mech. 37, 169–176 (2015) [CrossRef] [Google Scholar]
 V. Monchiet, FFT based iterative schemes for composites conductors with nonoverlapping fibers and Kapitza interface resistance, Int. J. Solids Struct. 135, 14–25 (2018) [CrossRef] [Google Scholar]
 Y. Benveniste, Model of thin interphases and the effective medium approximation in composite media with curvilinearly anisotropic coated inclusions, Int. J. Eng. Sci. 72, 140–154 (2013) [CrossRef] [Google Scholar]
 H. Le Quang, D.C. Pham, G. Bonnet, Q.C. He, Estimations of the effective conductivity of anisotropic multiphase composites with imperfect interfaces, Int. J. Heat Mass Transfer 58, 175–187 (2013) [CrossRef] [Google Scholar]
 H. Le Quang, T.L. Phan, G. Bonnet, Effective thermal conductivity of periodic composites with highly conducting imperfect interfaces, Int. J. Therm. Sci. 50, 1428–1444 (2011) [CrossRef] [Google Scholar]
 N. Trung Kien, N. Van Luat, P. Duc Chinh, Effective conductivity of isotropic composite with Kapitza thermal resistance, Vietnam J. Mech. 40, 377–386 (2018) [CrossRef] [Google Scholar]
Cite this article as: VanLuat Nguyen, Estimating the effective conductivity for ellipseinclusion model with Kapitza thermal resistance, EPJ Appl. Metamat. 8, 16 (2021)
All Tables
Comparison of results (C^{eff}) of FFT and PEK for square model, α = 1, a_{1} = 0.25.
Comparison of results (C^{eff}) of FFT and PEK for square model, α = 0.2, a_{1} = 0.25.
Comparison of results (C^{eff}) of FFT and PEK for square model, α = 0.05, a_{1} = 0.25.
Comparison of results (C^{eff}) of FFT, PEK and Monchiet's results for square model with one circle inclusion in the unit cell α = 1, a_{1} = a_{2} = R.
All Figures
Fig. 1
A single Coated ellipsoid. 

In the text 
Fig. 2
Unit cell: Coated ellipse (left), equivalent inclusion (right). 

In the text 
Fig. 3
Unit cell with square ellipse inclusion (left), equivalent inclusion (right). 

In the text 
Fig. 4
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 1, . 

In the text 
Fig. 5
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 1, . 

In the text 
Fig. 6
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.2, . 

In the text 
Fig. 7
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.2, . 

In the text 
Fig. 8
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.05, . 

In the text 
Fig. 9
Effective conductivity of square model with imperfect interface of ellipse Kapitza's type, α = 0.05, . 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.