 The smooth particle hydrodynamics, SPH, method was used to predict the effective thermal conductivity of porous media. The accuracy of the simulations was verified against ideal homogenous and multi-layered media. A concise method was also developed to generate unsaturated media using the SPH framework. The results of the simulation showed that the reciprocal of the normalized thermal conductivity and the reciprocal of the saturation are linearly related, which is consistent with the empirical formula proposed by Coat and Conrad. This article was authored by Bing Bai, Yan Wang, Ding Yu Rao, and others.