Research Papers: Conduction

Conduction in Jammed Systems of Tetrahedra

[+] Author and Article Information
Kyle C. Smith

e-mail: kyle.c.smith@gmail.com

Timothy S. Fisher

e-mail: tsfisher@purdue.edu
Birck Nanotechnology Center and School of Mechanical Engineering,
Purdue University,
West Lafayette, IN 47907

Contributed by the Heat Transfer Division of ASME for publication in the JOURNAL OF HEAT TRANSFER. Manuscript received June 13, 2012; final manuscript received April 1, 2013; published online June 27, 2013. Assoc. Editor: Jose L. Lage.

J. Heat Transfer 135(8), 081301 (Jun 27, 2013) (7 pages) Paper No: HT-12-1285; doi: 10.1115/1.4024276 History: Received June 13, 2012; Revised April 01, 2013

Control of transport processes in composite microstructures is critical to the development of high-performance functional materials for a variety of energy storage applications. The fundamental process of conduction and its control through the manipulation of granular composite attributes (e.g., grain shape) are the subject of this work. We show that athermally jammed packings of tetrahedra with ultrashort range order exhibit fundamentally different pathways for conduction than those in dense sphere packings. Highly resistive granular constrictions and few face–face contacts between grains result in short-range distortions from the mean temperature field. As a consequence, ‘granular’ or differential effective medium theory predicts the conductivity of this media within 10% at the jamming point; in contrast, strong enhancement of transport near interparticle contacts in packed-sphere composites results in conductivity divergence at the jamming onset. The results are expected to be particularly relevant to the development of nanomaterials, where nanoparticle building blocks can exhibit a variety of faceted shapes.

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


Shakouri, A., 2011, “Recent Developments in Semiconductor Thermoelectric Physics and Materials,” Annu. Rev. Mater. Res., 41, pp. 399–431. [CrossRef]
Smith, K. C., Mukherjee, P. P., and Fisher, T. S., 2012, “Columnar Order in Jammed LiFePO4 Cathodes: Ion Transport Catastrophe and Its Mitigation,” Phys. Chem. Chem. Phys., 14, pp. 7040–7050. [CrossRef]
Schlapbach, L., and Züttel, A., 2001, “Hydrogen-Storage Materials for Mobile Applications,” Nature, 414, pp. 353–358. [CrossRef] [PubMed]
Smith, K. C., and Fisher, T. S., 2012, “Models for Metal Hydride Particle Shape, Packing, and Heat Transfer,” Int. J. Hydrogen Energy., 37(18), pp. 13417–13428. [CrossRef]
Baxter, J., Bian, Z., Chen, G., Danielson, D., Dresselhaus, M. S., Fedorov, A. G., Fisher, T. S., Jones, C. W., Maginn, E., Kortshagen, U., Manthiram, A., Nozik, A., Rolison, D. R., Sands, T., Shi, L., Sholl, D., and Wu, Y., 2009, “Nanoscale Design to Enable the Revolution in Renewable Energy,” Energy Environ. Sci., 2, pp. 559–588. [CrossRef]
Maxwell, J. C., 1873, Treatise on Electricity and Magnetism, Clarendon Press, Oxford.
Rayleigh, L., 1892, “On the Influence of Obstacles Arranged in a Rectangular Order Upon the Properties of a Medium,” Philos. Mag., 34, pp. 481–502.
Torquato, S., 2002, Random Heterogeneous Materials, Spring-Verlag, New York, pp. 467–468.
Bruggeman, D. A. G., 1935, “Calculation of Various Physics Constants in Heterogenous Substances I. Dielectricity Constants and Conductivity of Mixed Bodies From Isotropic Substances,” Ann. Phys., 24(7), pp. 636–664. [CrossRef]
Yonezawa, F., and Cohen, M. H., 1983, “Granular Effective Medium Approximation,” J. Appl. Phys., 54, p. 2895. [CrossRef]
Sen, P. N., Scala, C., and Cohen, M. H., 1981, “A Self-Similar Model for Sedimentary-Rocks With Application to the Dielectric-Constant of Fused Glass-Beads,” Geophysics, 46(5), pp. 781–795. [CrossRef]
Hasselman, D. P. H., and Johnson, L. F., 1987, “Effective Thermal-Conductivity of Composites With Interfacial Thermal Barrier Resistance,” J. Compos. Mater., 21(6), pp. 508–515. [CrossRef]
Nan, C.-W., Birringer, R., Clarke, D. R., and Gleiter, H., 1997, “Effective Thermal Conductivity of Particulate Composites With Interfacial Thermal Resistance,” J. Appl. Phys., 81(10), pp. 6692–6699. [CrossRef]
Every, A. G., Tzou, Y., Hasselman, D. P. H., and Raj, R., 1992, “The Effect of Particle-Size on the Thermal-Conductivity of ZnS Diamond Composites,” Acta Metall. Mater., 40(1), pp. 123–129. [CrossRef]
Tsotsas, E., and Martin, H., 1987, “Thermal Conductivity of Packed Beds: A Review,” Chem. Eng. Process., 22, pp. 19–37. [CrossRef]
Bauer, R., and Schlünder, E. U., 1978, “Effective Radial Thermal Conductivity of Packings in Gas Flow—Part II: Thermal Conductivity of Packing Fraction Without Gas Flow,” Int. Chem. Eng., 18, pp. 189–204.
Cates, M. E., Wittmer, J. P., Bouchaud, J.-P., and Claudin, P., 1998, “Jamming, Force Chains, and Fragile Matter,” Phys. Rev. Lett., 81(9), pp. 1841–1844. [CrossRef]
Haji-Akbari, A., Engel, M., Keys, A. S., Zheng, X., Petschek, R. G., Palffy-Muhoray, P., and Glotzer, S. C., 2009, “Disordered, Quasicrystalline and Crystalline Phases of Densely Packed Tetrahedra,” Nature, 462, pp. 773–777. [CrossRef] [PubMed]
Smith, K. C., Alam, M., and Fisher, T. S., 2010, “Athermal Jamming of Soft Frictionless Platonic Solids,” Phys. Rev. E, 82(5), p. 051304. [CrossRef]
Smith, K. C., Fisher, T. S., and Alam, M., 2011, “Isostaticity of Constraints in Amorphous Jammed Systems of Soft Frictionless Platonic Solids,” Phys. Rev. E, 84, p. 030301. [CrossRef]
Ma, Y., Hao, Q., Poudel, B., Lan, Y., Yu, B., Wang, D., Chen, G., and Ren, Z., 2008, “Enhanced Thermoelectric Figure-of-Merit in p-Type Nanostructured Bismuth Antimony Tellurium Alloys Made From Elemental Chunks,” Nano Lett., 8(8), pp. 2580–2584. [CrossRef] [PubMed]
Zebarjadi, M., Esfarjani, K., Bian, Z., and Shakouri, A., 2011, “Low-Temperature Thermoelectric Power Factor Enhancement by Controlling Nanoparticle Size Distribution,” Nano Lett., 11(1), pp. 225–230. [CrossRef] [PubMed]
Jeng, M.-S., Yang, R., Song, D., and Chen, G., 2008, “Modeling the Thermal Conductivity and Phonon Transport in Nanoparticle Composites Using Monte Carlo Simulation,” ASME J. Heat Transfer, 130(4), p. 042410. [CrossRef]
Wang, X., Yang, Y., and Zhu, L., 2011, “Effect of Grain Sizes and Shapes on Phonon Thermal Conductivity of Bulk Thermoelectric Materials,” J. Appl. Phys., 110(2), p. 024312. [CrossRef]
Hsieh, T.-Y., Yang, J.-Y., and Hong, Z.-C., 2009, “Thermal Conductivity Modeling of Compacted Type Nanocomposites,” J. Appl. Phys., 106(2), p. 023528. [CrossRef]
Zhou, J., Li, X., Chen, G., and Yang, R., 2010, “Semiclassical Model for Thermoelectric Transport in Nanocomposites,” Phys. Rev. B, 82, p. 115308. [CrossRef]
Jaoshvili, A., Esakia, A., Porrati, M., and Chaikin, P. M., 2010, “Experiments on the Random Packing of Tetrahedral Dice,” Phys. Rev. Lett., 104(18), p. 185501. [CrossRef] [PubMed]
Keller, J. B., 1963, “Conductivity of a Medium Containing a Dense Array of Perfectly Conducting Spheres or Cylinders or Nonconducting Cylinders,” J. Appl. Phys., 34, pp. 991–993. [CrossRef]
Jacobson, M., and Völker, M., 2011, “N-Dimensional Sparse Arrays,” MATLAB Central File Exchange, http://www.mathworks.com/matlabcentral/fileexchange/29832, Retrieved April 29, 2013.
McPhedran, R. C., and McKenzie, D. R., 1978, “The Conductivity of Lattices of Spheres. I. The Simple Cubic Lattice,” Proc. R. Soc. London, Ser. A, 359, pp. 45–63. [CrossRef]
Jaluria, Y., and Torrance, K., 2003, Computational Heat Transfer (Series in Computational Methods in Mechanics and Thermal Sciences), Taylor & Francis, London.
Johansen, H., and Colella, P., 1998, “A Cartesian Grid Embedded Boundary Method for Poisson's Equation on Irregular Domains,” J. Comp. Phys., 147, pp. 60–85. [CrossRef]
Notay, Y., 2010, “An Aggregation-Based Algebraic Multigrid Method,” Electron. Trans. Numer. Anal., 37, pp. 123–146.
Napov, A., and Notay, Y., 2010, “An Algebraic Multigrid Method With Guarranteed Convergence Rate,” Technical Report No. GANMN 10-03.
Notay, Y., 2011, “Aggregation-Based Algebraic Multigrid for Convection-Diffusion Equations,” Technical Report No. GANMN 11-01.
O'Hern, C. S., Silbert, L. E., Liu, A. J., and Nagel, S. R., 2003, “Jamming at Zero Temperature and Zero Applied Stress: The Epitome of Disorder,” Phys. Rev. E, 68, p. 011306. [CrossRef]
Neudecker, M., Ulrich, S., Herminghaus, S., and Schröter, M., 2012, “Jammed Frictional Tetrahedra are Hyperstatic,” e-print arXiv:1202.6272.
Jiao, Y., and Torquato, S., 2011, “Maximally Random Jammed Packings of Platonic Solids: Hyperuniform Long-Range Correlations and Isostaticity,” Phys. Rev. E, 84, p. 041309. [CrossRef]
Choy, T. C., 1999, Effective Medium Theory: Principles and Applications, Clarendon Press, Oxford.


Grahic Jump Location
Fig. 1

Consolidation of 400 tetrahedra. Density is indicated below each image and color indicates the number of particles belonging to the face–face cluster of a given particle at the jamming threshold density of ϕJ = 0.634 for this system. Boundaries of the periodic supercell are indicated by black edges of the surrounding cube.

Grahic Jump Location
Fig. 2

Reconstructed interface geometry. (a) Phase of a given cell is determined by that of its centroid and is indicated by the centroid's color. The true and reconstructed interfaces are depicted in black and gray, respectively. (b) Sample cross-section of the three-level refined mesh for a system of 400 tetrahedra at ϕ = 0.611. The pore phase is indicated in white, while distinct grains are shaded.

Grahic Jump Location
Fig. 3

Octree refinement of κ¯ as a function of density for 400 tetrahedra with κg = 103. For each data point, the coarsest level of the mesh has a resolution of Δx=0.025(NVg/φ)1/3, where Δx is the length of coarse cube-shaped finite volume cells, N is the number of particles, and Vg is the volume of an individual grain. The position of the jamming point is marked by the red dashed line.

Grahic Jump Location
Fig. 4

Stencil employed for discretization of flux at faces between cells of dissimilar size

Grahic Jump Location
Fig. 5

Radial distribution function of the present systems for various system sizes. Experimental [27,37] and simulated [18,38] values of other contemporary works are shown as well. The system of 25 tetrahedra is excluded because of insufficient sampling. In brackets are listed the jamming threshold densities. Rmin is the minimum possible centroidal separation between contacting tetrahedra. The referenced prior simulation results are offset by an amount of two.

Grahic Jump Location
Fig. 6

Cluster number distribution for jammed systems of tetrahedra of varying system size

Grahic Jump Location
Fig. 7

(a) Normalized temperature field T˜ for 400 tetrahedra with κg = 103 at ϕ = 0.611. The field is induced by a homogenized temperature gradient (∇T)h applied to the system along the indicated direction. The length scale of a tetrahedron edge is indicated as a reference for the size of an individual grain. (b) Histogram of the deviation δT˜ from the mean normalized temperature field. The fitted normal distribution is shown along with its standard deviation σ and mean value μ.

Grahic Jump Location
Fig. 8

Mean effective conductivity for tetrahedra with κg = 103 as a function of (a) density and (b) density deviation, ϕJ − ϕ, where ϕJ is the jamming density for the respective system. The open triangular symbol represents a cubic polynomial extrapolation to the jamming point. 50:50 denotes the bidisperse mixture of tetrahedra described in the text. For simple cubic-packed spheres ϕJ = π/6, the density at which spheres just touch. GEMA refers to the granular effective medium approximation presented in Refs. [10,11], and MG-EMA refers to the Maxwell-Garnett effective medium approximation (see Ref. [39]). Data for simple cubic-packed spheres are from [30] at low density and [28] at high density.

Grahic Jump Location
Fig. 9

Deviatoric conductivity Δκ as a function of ϕ for κg = 103. The open triangular symbol represents a cubic polynomial extrapolation to the jamming point, and 50:50 denotes the bidisperse mixture of tetrahedra described in the text.

Grahic Jump Location
Fig. 10

Mean effective conductivity as a function of grain conductivity κg for (a) 100 tetrahedra at ϕ = 0.626 and (b) 400 tetrahedra at ϕ = 0.611. Error bars represent κ1 and κ3. Theoretical data are presented based on the GEMA in Refs. [10,11] and the Zehner–Bauer–Schlünder (ZBS) model [15,16].




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