Bridging micro nature with macro behaviors for granular thermal mechanics

The connection between micro-level characteristics and macroscopic properties in granular heat transfer and mechanics is fundamental and crucial. This study proposes a novel discrete element approach incorporating granular heat transfer, contact bonding, and granular stress tensor models to investigate the mechanical and thermal responses of continuum media composed of constituent spheres. Eight benchmark tests were devised to bridge the long-standing gap between micro and macro properties in granular materials. Through these tests, the numerical solutions obtained from discrete element modeling match well with existing analytical or finite element solutions derived from continuum-based theory. This validation underscores the rationality and reliability of the granular heat transfer model, contact bonding model, and granular stress tensor model. Moreover, the study highlights the consistency between continuum-based theory and discontinuum-based theory. A minor distinction between continuum-based models and discrete element models emerges near the boundaries due to variations in the specification of boundary conditions. This discrepancy can be clarified by Saint-Venant's Principle, thus validating the accuracy of the microscale heat transfer and mechanics theory for granular materials. Five mono-disperse packing structures, including simple cubic (SC), body-centered cubic (BCC), facecentered cubic (FCC), hexagonal close packing (HCP), and random packing (Random), were further analyzed to examine their influence on heat transfer performance. Numerical results reveal that higher coordination numbers and solid volume fractions correspond to higher apparent thermal conductivity of granular assemblies, thus elucidating the connection between micro packing configurations and macroscopic heat transfer properties. The apparent thermal conductivity for different crystal configurations follows the sequence: HCP & efDot; FCC > BCC & efDot; Random > SC. To improve the accuracy and physical relevance of the proposed model, the effect of particle contact area needs to be further incorporated into the granular heat transfer model.