Estimation of Effective Thermal and Mechanical Properties of Particulate Thermal Interface Materials by a Random Network Model

[+] Author and Article Information
Pavan Kumar Vaitheeswaran

Department of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907
e-mail: pvaithee@purdue.edu

Ganesh Subbarayan

Department of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907
e-mail: gss@purdue.edu

1Corresponding author.

Contributed by the Electronic and Photonic Packaging Division of ASME for publication in the JOURNAL OF ELECTRONIC PACKAGING. Manuscript received September 20, 2017; final manuscript received January 12, 2018; published online May 9, 2018. Assoc. Editor: Reza Khiabani.

J. Electron. Packag 140(2), 020901 (May 09, 2018) (7 pages) Paper No: EP-17-1092; doi: 10.1115/1.4039136 History: Received September 20, 2017; Revised January 12, 2018

Particulate thermal interface materials (TIMs) are commonly used to transport heat from chip to heat sink. While high thermal conductance is achieved by large volume loadings of highly conducting particles in a compliant matrix, small volume loadings of stiff particles will ensure reduced thermal stresses in the brittle silicon device. Developing numerical models to estimate effective thermal and mechanical properties of TIM systems would help optimize TIM performance with respect to these conflicting requirements. Classical models, often based on single particle solutions or regular arrangement of particles, are insufficient as real-life TIM systems contain a distribution of particles at high volume fractions, where classical models are invalid. In our earlier work, a computationally efficient random network model (RNM) was developed to estimate the effective thermal conductivity of TIM systems (Kanuparthi et al., 2008, “An Efficient Network Model for Determining the Effective Thermal Conductivity of Particulate Thermal Interface Materials,” IEEE Trans. Compon. Packag. Technol., 31(3), pp. 611–621; Dan et al., 2009, “An Improved Network Model for Determining the Effective Thermal Conductivity of Particulate Thermal Interface Materials,” ASME Paper No. InterPACK2009-89116.) . This model is extended in this paper to estimate the effective elastic modulus of TIMs. Realistic microstructures are simulated and analyzed using the proposed method. Factors affecting the modulus (volume fraction and particle size distribution (PSD)) are also studied.

Copyright © 2018 by ASME
Your Session has timed out. Please sign back in to continue.


Maxwell, J. C. , 1873, Treatise on Electricity and Magnetism, Clarendon Press, Oxford, UK.
Rayleigh, L. , 1892, “On Influence of Obstacles Arranged in Rectangular Order Upon the Properties of a Medium,” Philos. Mag., 34(211), pp. 481–502. [CrossRef]
Bruggeman, D. , 1935, “Berechnung Verschiedener Physikalischer Konstanten Von Heterogenen Substanzen,” Ann. Phys. (Liepzig), 24(7), pp. 636–679. [CrossRef]
Kanuparthi, S. , Rayasam, M. , Subbarayan, G. , Sammakia, B. , Gowda, A. , and Tonapi, S. , 2009, “Hierarchical Field Composition for Simulations of Near-Percolation Thermal Transport in Particulate Materials,” Comput. Methods Appl. Mech. Eng., 198(5–8), pp. 657–668. [CrossRef]
Batchelor, G. K. , and O'Brien, R. W. , 1977, “Thermal or Electrical Conduction Through a Granular Material,” Proc. R. Soc. London, Ser. A, 355(1682), pp. 313–333.
Kanuparthi, S. , Subbarayan, G. , Siegmund, T. , and Sammakia, B. , 2008, “An Efficient Network Model for Determining the Effective Thermal Conductivity of Particulate Thermal Interface Materials,” IEEE Trans. Compon. Packag. Technol., 31(3), pp. 611–621.
Dan, B. , Kanuparthi, S. , Subbarayan, G. , and Sammakia, B. , 2009, “An Improved Network Model for Determining the Effective Thermal Conductivity of Particulate Thermal Interface Materials,” ASME Paper No. InterPACK2009-89116.
Fu, S.-Y. , Feng, X.-Q. , Lauke, B. , and Mai, Y.-W. , 2008, “Effects of Particle Size, Particle/Matrix Interface Adhesion and Particle Loading on Mechanical Properties of Particulate-Polymer Composites,” Composites, Part B, 39(6), pp. 933–961. [CrossRef]
Einstein, A. , 1905, “Ueber Die Von Der Molekularkinetischen Fluessigkeiten Suspendierten Teilchen,” Ann. Phys. (Leipzig), 322(8), pp. 549–560. [CrossRef]
Guth, E. , 1945, “Theory of Filler Reinforcement,” J. Appl. Phys., 16(1), pp. 20–25. [CrossRef]
Kerner, E. H. , 1956, “The Elastic and Thermoelastic Properties of Composite Media,” Proc. Phys. Soc. B, 69(8), pp. 808–813. [CrossRef]
Torquato, S. , 2002, Random Heterogeneous Materials, Springer-Verlag, New York. [CrossRef]
Hashin, Z. , 1962, “The Elastic Moduli of Heterogeneous Materials,” ASME J. Appl. Mech., 29(1), pp. 143–150. [CrossRef]
Hashin, Z. , and Shtrikman, S. , 1963, “A Variational Approach to the Theory of the Elastic Behaviour of Multiphase Materials,” J. Mech. Phys. Solids, 11(2), pp. 127–140. [CrossRef]
Budiansky, B. , 1965, “On the Elastic Moduli of Some Heterogeneous Materials,” J. Mech. Phys. Solids, 13(4), pp. 223–227. [CrossRef]
Hill, R. , 1965, “A Self-Consistent Mechanics of Composite Materials,” J. Mech. Phys. Solids, 13(4), pp. 213–222. [CrossRef]
Boucher, S. , 1974, “On the Effective Moduli of Isotropic Two-Phase Elastic Composites,” J. Compos. Mater., 8(1), pp. 82–89. [CrossRef]
McLaughlin, R. , 1977, “A Study of the Differential Scheme for Composite Materials,” Int. J. Eng. Sci., 15(4), pp. 237–244. [CrossRef]
Eshelby, J. D. , 1957, “The Determination of the Elastic Field of an Ellipsoidal Inclusion, and Related Problems,” Proc. R. Soc. A, 241(1226), pp. 376–396. [CrossRef]
Halpin, J. C. , 1969, “Stiffness and Expansion Estimates for Oriented Short Fiber Composites,” J. Compos. Mater., 3(4), pp. 732–734. [CrossRef]
Mooney, M. , 1951, “The Viscosity of a Concentrated Suspension of Spherical Particles,” J. Colloid Sci., 6(2), pp. 162–170. [CrossRef]
Santio, E. , and Muller, E. A. , 2002, “Dense Packing of Binary and Polydisperse Hard Spheres,” Mol. Phys., 100(15), pp. 2461–2469. [CrossRef]
Zhang, X. , 2004, “Constructive Modeling Strategies and Implementation Frameworks for Optimal Hierarchical Synthesis,” Ph.D. thesis, Purdue University, West Lafayette, IN. https://docs.lib.purdue.edu/dissertations/AAI3252937/
Rayasam, M. , Srinivasan, V. , and Subbarayan, G. , 2007, “CAD Inspired Hierarchical Partition of Unity Constructions for NURBS-Based, Meshless Design, Analysis and Optimization,” Int. J. Numer. Methods Eng., 72(12), pp. 1452–1489. [CrossRef]
Smith, J. C. , 1976, “Experimental Values for the Elastic Constants of a Particulate-Filled Glassy Polymer,” J. Res. Natl. Bur. Stand., Sect. A, 80(1), pp. 45–49. http://nvlpubs.nist.gov/nistpubs/jres/80A/jresv80An1p45_A1b.pdf
Nie, S. , and Basaran, C. , 2005, “A Micromechanical Model for Effective Elastic Properties of Particulate Composites With Imperfect Interfacial Bonds,” Int. J. Solids Struct., 42(14), pp. 4179–4191. [CrossRef]
Braem, M. , Doren, V. E. V. , Lambrechts, P. , and Vanherle, G. , 1987, “Determination of Young's Modulus of Dental Composites,” J. Mater. Sci., 22(6), pp. 2037–2042. [CrossRef]


Grahic Jump Location
Fig. 1

RVE of microstructure (left) with 58% volume fraction of alumina particles in silicone matrix with (right) heat flux field observed at the midsection (y = 0.5) using HPFC

Grahic Jump Location
Fig. 2

Equivalent network of conductors for a given microstructure

Grahic Jump Location
Fig. 3

Interparticle conductance as a series combination of gap conductance and particle conductances

Grahic Jump Location
Fig. 4

Effective conductance for the interaction between a particle and a wall

Grahic Jump Location
Fig. 5

Comparing the results of the random network model with those obtained using full-field simulations

Grahic Jump Location
Fig. 6

Heat flux transport within particle extended to entire hemisphere instead of just a cylindrical region

Grahic Jump Location
Fig. 7

Force transfer across the matrix as an effective gap stiffness

Grahic Jump Location
Fig. 8

Particle divided into infinitesimal stiffnesses in series

Grahic Jump Location
Fig. 9

Illustration of the drop-fall-shake procedure used to generate realistic random microstructures from a given particle size distribution

Grahic Jump Location
Fig. 10

Random network model predictions against experimental measurements. Accuracy improves with increase in particle volume fractions.

Grahic Jump Location
Fig. 11

The effective elastic modulus of the system increases exponentially with filler fraction

Grahic Jump Location
Fig. 12

Different log-normal particle size distributions generated and analyzed

Grahic Jump Location
Fig. 13

The effective elastic modulus shows a slight increase with polydispersity at high volume fractions



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In