Скорость изменения температуры формула

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

> = alpha Aleft( <– T>
ight),] где (Q) − количество теплоты, (A) − площадь поверхности тела, через которую передается тепло, (T) − температура тела, (<
>) − температура окружающей среды, (alpha) − коэффициент теплопередачи , зависящий от геометрии тела, состояния поверхности, режима теплопередачи и других факторов.

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

><
> = frac<<alpha A>>left( <– T>
ight) = kleft( <
– T>
ight).] Решение данного уравнение имеет вид: [Tleft( t
ight) =
+ left( < >
ight)>,] где () обозначает начальную температуру тела.

Таким образом, температура тела уменьшается экспоненциально по мере охлаждения, приближаясь к температуре окружающей среды. Скорость охлаждения зависит от параметра (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<

><
> = kleft( <– T>
ight).] По условию задачи, (
= > + eta t.) Следовательно, последнее уравнение можно записать в виде: [ <frac<
><
> = kleft( <> + eta t – T>
ight)>;; < ext<или>;;T’ + kT = k
> + keta t.> ] Мы получили линейное дифференциальное уравнение, которое можно решить, например, с помощью интегрирующего множителя: [uleft( t
ight) = >> = >.] Общее решение уравнения записывается в форме [ >left(
> + keta t>
ight)dt> + C>><<>>> > = <frac<
>int <>dt> + keta int <>tdt> + C>><<>>>.> ] Второй интеграл в числителе находится интегрированием по частям: [ <int <underbrace <>>_>>\ >>\ \