В конце (17) века британский ученый Исаак Ньютон изучал охлаждение тел. Эксперименты показали, что скорость охлаждения примерно пропорциональна разнице температур между нагретым телом и окружающей средой. Этот факт можно записать в виде дифференциального уравнения: [frac<><

Поскольку (Q = CT,) где (C) − теплоемкость тела, то дифференциальное уравнение можно записать как [frac<

Таким образом, температура тела уменьшается экспоненциально по мере охлаждения, приближаясь к температуре окружающей среды. Скорость охлаждения зависит от параметра (k = largefrac<<alpha A>>
ormalsize) ( коэффициента теплопроводности ). С увеличением коэффициента (k) (например, вследствие увеличения площади поверхности), тело будет охлаждаться быстрее (рисунок (1).)

Мы решим задачу сначала для случая произвольной температуры окружающей среды, а затем вычислим конечную температуру тела при температуре среды (0^circ.)

Пусть начальная температура нагретого тела составляет ( = 200^circ.) Последующее изменение температуры описывается формулой: [ + left( <– >
ight)> > = <+ left( <200^circ – >
ight)>.> ] В конце первого часа тело охладилось до (100^circ.) Следовательно, можно записать следующее соотношение: [
ight) = 100^circ = + left( <200^circ – >
ight)>,>;; <Rightarrow 100^circ = + left( <200^circ – >
ight)>.> ] Спустя (2) часа температура тела становится равной (X) градусов: [X = + left( <200^circ – >
ight)>.] Таким образом, мы получили систему двух уравнений с тремя неизвестными: (,) (k) и (X:) [left< egin 100 = + left( <200 – >
ight)>\ X = + left( <200 – >
ight)> end
ight..] Мы не можем однозначно определить температуру тела (X) через (2) часа из данной системы. Однако можно вывести зависимость (X) от температуры окружающей среды (.) Выразим функцию (>) из первого уравнения: [> = frac<<100 – >><<200 – >>.] Тогда [> = <left( <>>
ight)^2> = <left( <frac<<100 – >><<200 – >>>
ight)^2>.] Следовательно, зависимость (Xleft( <>
ight)) имеет вид: [ >
ight) = + left( <200 – >
ight)<left( <frac<<100 – >><<200 – >>>
ight)^2> > = <+ frac<<<<left( <100 – >
ight)>^2>>><<200 – >>.> ] Если, например, положить температуру окружающей среды равной нулю, то температура тела через два часа будет составлять [ = 0>
ight) = 0 + frac <<<<left( <100 – 0>
ight)>^2>>><<200 – 0>> > = <frac<<10000>><<200>> = 50^circ.> ] Зависимость температуры тела (X) от температуры окружающей среды в данной задаче показана выше на рисунке (2.)

Прежде всего, отметим разницу со случаем когда тело охлаждается в среде, температура которой постоянна. В этом случае температура тела формально будет бесконечно долго приближаться к температуре окружающей среды. В нашей же задаче температура среды линейно возрастает. Поэтому, рано или поздно температура тела станет равной температуре среды, то есть задача имеет решение. Будем считать также, что соблюдается квазистационарный режим, т.е. все переходные процессы в системе быстро затухают.

В таком случае процесс можно описать дифференциальным уравнением: [frac<