Analysis of diffusion delay in a layered medium. Application to heat measurements from muscle

An analysis is presented of diffusional delays in one-dimensional heat flow through a medium consisting of several layers of different materials. The model specifically addresses the measurement of heat production by muscle, but diffusion of solute or conduction of charge through a layered medium will obey the same equations. The model consists of a semi-infinite medium, the muscle, in which heat production is spacially uniform but time varying. The heat diffuses through layers of solution and insulation to the center of the thermal element where heat flow is zero. Using Laplace transforms, transfer functions are derived for the temperature change in the center of the thermopile as a function of the temperature at any interface between differing materials or as a function of heat production in the muscle. From these transfer functions, approximate analytical expressions are derived for the time constants which scale the early and late changes in the central temperature. We find that the earliest temperature changes are limited by the diffusivities of the materials, whereas the approach to steady state depends on the total heat capacity of the system and the diffusivity of muscle. Hill (1937) analyzed a similar geometry by modeling the layered medium as a homogeneous system with an equivalent half thickness. We show that his analysis was accurate for the materials in his system. In general, however, and specifically with regard to modern thermopiles, a homogeneous approximation will lead to significant errors. We compare responses of different thermopiles to establish the limits of time resolution in muscle heat records and to correct them for diffusional delays. Using numerical techniques, we invert the Laplace transforms and show the time course of the temperature changes recorded by different instruments in response to different patterns of heat production.