W e present numerical calculations of the effective thermal conductivity of unidirectional brous composite materials with an interfacial thermal resist ance bet ween the cont inuous and dispersed components. W e rst develop a continuous variat ional formulation of the heat conduction problem and then apply the method of homogenization to derive the appropriat e cell problem. Next , we solve the latter using the nite element met hod and present and validate effective conductivity results for ordered composit es. Finally, to highlight the wide applicability of our approach, we also present a posteriori contact conduct ance results for a disordered composite whose microstructure is based on a real optical micrograph.
INTROD UCTIONThe determination of the macroscopic behavior of composite materials is of fundamental and practical importance, particularly now, as new composites are being developed continually for the automotive and aerospace industries, electronic packaging, thermal insulation, and many other applications [1^4]. Typical microstructures of composites consist of unidirectional long ¢bers (mono¢laments), short chopped ¢bers (whiskers), or particles or voids dispersed in a solid matrix. One of the most important macroscopic thermal properties of composite materials is the effective thermal conductivity, which depends on the conductivities, relative amounts, geometries, and spatial distributions of the constituent phases. A recent comprehensive review of heat conduction in composites is presented by Furma·ski [3].Most manufacturing processes of composite materials do not ensure a perfect thermal contact between the constituent phases [1, 2, 5^8]. Therefore, in practice, the effective conductivity also depends on the interfacial thermal resistance, or contact resistance, between the continuous phase (the matrix) and the dispersed phase.