Take the 2-minute tour ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free, no registration required.

Yesterday I randomly started thinking about Newton's Law of Cooling. The problem I realized is that it assumes the ambient temperature stays constant over time, which is obviously not true. So what I tried doing was to modify the differential equation into a system of differential equations, and taking the heat capacity of each into account. $$ \frac{dT_1(t)}{dt} = C_2(T_2(t)-T_1(t))\\ \frac{dT_2(t)}{dt} = C_1(T_1(t)-T_2(t))\\ T_1(0) = T_{10}\\ T_2(0) = T_{20} $$ which had the following solution: $$ T_1(t)={\frac {C_{{2}}T_{{20}}+C_{{1}}T_{{10}}}{C_{{1}}+C_{{2}}}}+{\frac {C_{ {2}} \left( -T_{{20}}+T_{{10}} \right) {{\rm e}^{- \left( C_{{1}}+C_{{ 2}} \right) t}}}{C_{{1}}+C_{{2}}}}\\ T_2(t)={\frac {C_{{2}}T_{{20}} +C_{{1}}T_{{10}}}{C_{{1}}+C_{{2}}}}-{\frac {C_{{1}} \left( -T_{{20}}+T_{{10}} \right) {{\rm e}^{ -\left(C_{{1}} +C_{{2}} \right) t}}}{C_{{1}}+C_{{2}}}} $$

When I plotted the two equations out they seem to be right. They also follow the fact that at any point in time $C_1\Delta T_1=C_2\Delta T_2$. The end behaviour also seem to be correct. However, in cases where $C_1\rightarrow\infty$ you would expect it to behave similarly to Newton's Law of Cooling, but in reality $T_2$ drops its temperature in a very short period of time, which doesn't seem to be right. I tried looking this up but couldn't find much on this topic. If anyone can point out what I did wrong that would be great. I've had very little experience with thermodynamics but I do realize how limited these models are at describing real world scenarios. This was more for the fun of playing with differential equations.

share|improve this question
add comment

1 Answer

up vote 6 down vote accepted

The flaw in your reasoning seems to be that $C$ is not in fact heat capacity. In Newton's Law of Cooling, the proportionality constant would be related inversely to the heat capacity of the two heated liquids/gasses/materials, and directly to the heat conductance of the object separating the two materials. A material with a higher heat capacity would have a smaller temperature change for a given temperature difference, and a thin piece of metal separating the materials would result in a much larger $\frac{dT}{dt}$ than a thick piece of styrofoam would.

When bringing $C_1$ to $\infty$, you are actually decreasing the heat capacity and increasing the conductance, both of which would cause $T_2$ to drop quickly as you observed.

share|improve this answer
    
Ah thank you. I modified the equation to be $\frac{dT_1(t)}{dt} = \frac{1}{C_2}(T_2(t)-T_1(t))$ and it looks much better now. As $C_1$ approaches $\infty$ the reciprocal approaches zero and that term is effectively gone. Just wondering, are there any literature online that's related these types of modelling? –  Neo Dec 31 '12 at 23:43
    
I'm not familiar with any specific literature on the topic, although the general field is computational physics. I too am interested in modeling physics problems with a computer (right now, I'm attempting to model a hydrogen atom). –  Draksis Dec 31 '12 at 23:50
    
Yeah modelling is extremely fun. Unfortunately I'm in first year and there's not much you can do with high school level mathematics. Onwards to MIT OCW! –  Neo Jan 1 '13 at 0:11
    
I created a predator-prey simulation and many-body gravity simulator, both of which were interesting simulations at a high school level. Also, you can look at scicomp.stackexchange.com. –  Draksis Jan 1 '13 at 16:51
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.