## Post-mortem body cooling in variable environments

In a recent post I noted that the standard post-mortem body cooling method used to estimate time of death (TOD) does not take into account a varying environmental temperature (i.e. it assumes a constant ambient temperature when applying Newton’s Law of Cooling).

Consider the differential equation

.

Suppose we take the difference between the body temperature and the ambient temperature, , as being and the time, , and being . Then

.

Since, as I noted, the ambient temperature is usually considered constant, the

term is eliminated and the usual solution is given.

So I figured I’d take some actual data from the NOAA/NWS website, throw it into Excel, and get a rough estimate of the function . The 24-hour data I chose from Sanford, Maine (nearest data to my house) fit a sixth-order polynomial with an value of 0.977. Since I then had a pretty good function for I simply plugged it (and ) in to the above to obtain a first-order non-homogenous differential equation of the form

.

Note that there are ‘s in there but not intuitively distributed. In terms of the variables in this situation, we have

.

In any case, I then plotted this for a 24-hour period assuming the body started out at 98.6º. Unfortunately, the data for is just nuts. In short, as one of my students noted, it is too strongly coupled to the ambient temperature for this instance. The question is: what’s wrong with the above and what is the correct way to handle this? I really think the solution, if found, will demonstrate that in some instances treating the ambient temperature as a constant results in an incorrect TOD.

## Comment (obtuse, impolite, or otherwise "troll"-like comments may be deleted)