Numerical Evaluation Of Effective Thermal Properties For Materials With Variable Porosity

Open access


In this paper a computational homogenization technique is applied to thermal analyses in porous materials. A volume fraction of pores on the microstructural level is the key factor that changes the macroscopic thermal properties. Thus, the distribution of thermal fields at the macroscopic level is analysed through the incorporation of the microstructural response on the representative volume element (RVE) assuming a uniform distribution of pores. For the numerical analysis the scaled boundary finite element method (SBFEM) is introduced to compute the thermal response of RVE. The SBFEM combines the main advantages of the finite element method (FEM) and the boundary element method (BEM). In this method, only the boundary is discretized with elements leading to the reduction of spatial dimension by one, similarly as in the BEM. It reduces computational efforts in the mesh generation and CPU time. The proposed method is used to study square RVE with a circular and elliptic pore under the thermal load. Dimensions of the pore are varied to obtain different volume fractions of matrix material. Numerical results for effective thermal conductivities obtained via SBFEM modelling show an excellent agreement with the finite element analysis using commercial software COMSOL Multiphysics.

[1] Abbasbandy, S. and Shirzadi, A. (2010), “A meshless method for two-dimensional diffusion equation with an integral condition”, Eng Anal Bound Elem 34, 1031-1037.

[2] Abreu, A. I., Canelas, A. and Mansur, W. J. (2013), “A CQM-based BEM for transient heat conduction problems in homogenous materials and FGMs”, Applied Mathematical Modelling 37, 776-792.

[3] Atluri, S. N. (2004), The meshless method (MLPG) for domain and BIE discretizations, Tech Science Press, Forsyth, USA.

[4] Bathe, K. J. (1976), Numerical Methods in Finite Element Analysis, Prentice-Hall, New Jersey, USA.

[5] Benveniste, Y. (1987), “A new approach to the application of Mori-Tanaka’s theory in composite materials”, Mechanics of Materials 6 (2), 147-157.

[6] Brebbia, C. A., Telles, J. C. and Wrobel, L. C. (1984), Boundary Elements Techniques: Theory and Applications in Engineering, Springer-Verlag, Berlin-New York.

[7] Cai, Z., Mandel, J. and McCormick, S. (1991), “The finite volume element method for diffusion equations on general triangulations”, SIAM Journal on Numerical Analysis 28, 392-403.

[8] Carslaw, HS. and Jaeger, J. C. (1959), Conduction of Heat in Solids, Oxford University Press, Vivian Ridler, Walton Street, Oxford, second edition.

[9] Gruescu, C., Giraud, A., Homand, F., Kondo, D. and Do, D. P. (2007), “Effective thermal conductivity of partially saturated porous rocks”, International Journal of Solids and Structures 44, 811-818.

[10] Hill, R. (1963), “Elastic properties of reinforced solids: some theoretical principles”, Journal of the Mechanics and Physics of Solids 11 (5), 357-372.

[11] Christensen, R. M. (1993), “Effective properties of composite materials containing voids”. Proceedings of the Royal Society A, 440 (1909), 461-473.

[12] Jasiuk, I., Chen, C. and Thorp, M. F. (1994), “Elastic moduli of two dimensional materials with polygonal and elliptical holes”, Applied Mechanics Review 47, 1S, 18-28.

[13] Juncu, Gh. (2008), “Unsteady conjugate forced convection heat/mass transfer from a finite flat plate”, International Journal of Thermal Sciences 47, 972-984.

[14] Li, C., Ooi, E.T, Song, C. and Natarajan, S. (2015), “SBFEM for fracture analysis of piezoelectric composites under thermal load”, Int J Solids Struct, 52, 114-129.

[15] Mori, T. and Tanaka, K. (1973), “Average Stress in the Matrix and Average Elastic Energy of Materials with Misfitting Inclusions”, Acta metallica 21, 571-574.

[16] Ochiai, Y. (2001), “Steady heat conduction analysis in orthotropic bodies by triple-reciprocity BEM”, Computer Modeling in Engineering and Sciences 2(4), 435-446.

[17] Ozisik, M. N. (1994), Finite Difference Methods in Heat Transfer, CRC Press Inc, Boca Raton, USA.

[18] Sladek, J., Sladek, V., Solek, P., Tan, C. L. and Zhang, Ch. (2009), “Two- and three-dimensional transient thermoelastic analysis by the MLPG method”, CMES: Computer Modeling in Engineering & Sciences 47, 61-95.

[19] Sladek, J., Sladek, V., Zhang, Ch. and Tan, C. L. (2006), “Meshless local Petrov-Galerkin method for linear coupled thermoelastic analysis”, CMES: Computer Modeling in Engineering & Sciences 16, 57-68.

[20] Sladek, J., Sladek, V., and Zhang, Ch. (2003), “Transient heat conduction analysis in functionally graded materials by the meshless local boundary integral equation method”, Comput Mat Sci 28, 494-504.

[21] Sladek, J., Stanak, P., Han, Z. D., Sladek, V. and Atluri, S. N. (2013), “Applications of the MLPG Method in Engineering & Sciences: A Review”, CMES: Computer Modeling in Engineering & Sciences 92(5), 423-475.

[22] Song, C. and Wolf, J. (1997), “The scaled boundary finite-element method- alias consistent infinitesimal finite-element cell method- for elastodynamics”, Comput Methods Appl Mech Engng 147, 329-355.

[23] Song, C., (2006), “Analysis of singular stress fields at multi-material corners under thermal loading”, Int J Numer Methods Eng 65(5), 620-652.

Journal Information


All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 227 227 9
PDF Downloads 67 67 0