Heat conduction in a semi-infinite medium with a spherical inhomogeneity and time-periodic boundary temperature

We solve the problem of heat conduction in a homogeneous media below a planar boundary subjected to time-periodic temperature (of frequency omega), in the presence of a spherical inhomogeneity (of radius R), whose center is at distance d > R from the boundary. In the absence of the sphere, the well known one dimensional solution can be regarded as an oscillating thermal boundary layer of displacement thickness delta = root 2 alpha/omega, where alpha is the heat diffusivity. The general solution depends on four dimensionless parameters: d/R, delta/R, the heat conductivity ratio K and the heat capacity ratio C. An analytical solution is derived as an infinite series of Bessel functions, which converges quickly. The results are illustrated and analyzed for a given accuracy and for a few values of the governing parameters. The general solution can be simplified considerably for asymptotic values of the parameters. A first approximation, obtained for R/d << 1, pertains to an unbounded domain. A further approximate solution, for R/delta << 1, while K and C are fixed, can be regarded as pertaining to a quasi-steady regime, and is similar in structure to Maxwell's solution for steady state. However, its accuracy deteriorates for kappa << 1, and a solution, coined as the insulated sphere approximation, is derived for this case. Comparison with the exact solution shows that these approximations are accurate for a wide range of parameter values. Besides providing insight, they can be employed for solving in a simple manner more complex problems, e.g. effective properties of a heterogeneous medium made of an ensemble of spherical inclusions. (C) 2011 Elsevier Ltd. All rights reserved.