Influence of boundary conditions on computation of the effective thermal conductivity of foams

Abstract

Accurate numerical simulation of the effective thermal conductivity (ETC) of 3D pore-scale foam models requires a judicious choice of boundary conditions, as the computational domains are often smaller than the representative volume element, giving rise to considerable edge effects. Within the finite element homogenization framework, a set of mixed boundary conditions are considered alongside the usual uniform and periodic boundary conditions. Validity criteria and order relations, demonstrated from entropy-based principles, are numerically verified on unit cell-based geometries, random virtual periodic foams, and non-periodic tomography-reconstructed foams of equivalent microstructure. A statistical treatment based on the integral range provides confidence intervals for the estimated ETC. For foam samples with random homogeneous porosity, the mixed boundary conditions are shown to fulfill the macrohomogeneity condition and thus provide thermodynamically valid ETC estimates. For periodic foams with irregular microstructure, the ETC is very slightly underestimated under the mixed boundary conditions. For non-periodic geometries, it is shown that periodic boundary conditions–commonly viewed as the reference–underestimate the ETC due to boundary geometry mismatch, while the mixed boundary conditions give a more accurate and precise estimate.