Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

8458

.pdf
Скачиваний:
0
Добавлен:
25.11.2023
Размер:
1.63 Mб
Скачать

d 2 T

 

q V

0 .

(3.48)

dx 2

λ

 

 

 

Будем рассматривать лишь правую половину пластины, поскольку поле температуры симметрично относительно плоскости х=0. Тогда граничные условия имеют вид

x 0;

dT

 

0;

 

 

 

 

 

 

 

 

 

dx x 0

(3.49)

 

 

dT

x δ;

 

α Tx δ Tж .

λ

 

 

 

 

 

dx x δ

После интегрирования (3.48) получим:

dT

 

q V x

 

C1

;

(3.50)

dx

λ

 

 

 

 

 

 

Т

q V x 2

 

C1x C 2 .

(3.51)

 

 

 

 

 

 

 

 

 

Постоянные интегрирования С1 и С2 найдем из граничных условий (3.49). При х=0 из уравнения (3.50) получаем С1=0. Подставив уравнения (3.50) и (3.51) в (3.49) при х= , получим

 

 

q

V

δ

 

 

 

q

V

δ 2

λ

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

Откуда имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

C2

Tж

q V δ

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

C2 Tж .

q V δ 2 .

Подставив значения постоянных С1 и С2 в выражение (3.51), получим уравнение температурного поля

 

q V δ

 

q V δ

2

 

x

2

 

 

 

 

 

 

 

 

 

Т Т ж

 

 

 

 

1

 

 

 

 

.

(3.52)

 

 

 

 

 

 

α

 

 

 

 

δ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из соотношения (3.50) имеем, что плотность теплового потока изменяется по толщине стенки линейно

 

T

q V x.

 

q -λ

 

 

 

x

 

 

Плотность теплового потока на поверхности пластины равна

 

q α Tc Tж

q V δ.

(3.53)

Общий поток тепла, отданный обеими поверхностями стенки, составит

Q q 2F q V δ2F.

Если в уравнении (3.52) положить , то граничные условия третьего рода вырождаются в граничные условия первого рода, ибо при получим Тжс. Тогда из уравнения (3.52) имеем

31

 

q V δ

2

 

x

2

 

 

 

 

 

 

 

 

T TC

 

 

1

 

 

 

 

.

(3.54)

 

 

 

 

 

 

 

 

δ

 

 

 

 

 

 

 

 

 

 

 

 

 

Из уравнений (3.52) и (3.54) можно определить максимальную температуру в центре пластины при х=0.

Для цилиндра, как и для пластины, задача одномерна и симметрична. Уравнение (3.48) при этом имеет вид

 

d 2 T

 

 

1 dT

 

q V

 

 

0.

(3.55)

 

dr 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r dr

 

 

λ

 

 

 

 

 

 

Граничные условия:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при r 0

dT

 

 

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr

r 0

 

 

 

 

 

 

 

 

 

 

при r r0

dT

 

 

 

 

α

TC Tж .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.56)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dr r r

 

 

 

 

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

После интегрирования (3.55), используя подстановку dT/dr=u, получим

 

 

dT

 

q V r

 

C1

 

;

(3.57)

 

 

 

 

 

 

 

 

 

 

 

 

dr

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

T

 

q V r

2

 

C1lnr C2 .

(3.58)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Постоянные С1 и С2 определяются из граничных условий (3.56). После их подстановки в (3.57) и (3.58) получим окончательно

dT q V r ; dr 2λ

 

q

r

2

 

 

 

r

2

 

 

q

r

T

 

V 0

 

 

 

 

 

 

 

 

V 0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На поверхности цилиндра плотность теплового соответственно равны

 

(3.59)

Tж .

(3.60)

потока и

поток тепла

q α T

T

 

q V r0

,

 

(3.61)

 

 

C

ж

2

 

 

 

 

 

 

 

 

Q qF

q V r0

2ππ l q

πr 2l.

(3.62)

 

 

2

0

 

V 0

 

 

 

 

 

 

 

Если в полученных решениях полагать = , то граничные условия третьего рода вырождаются в граничные условия первого рода Тжс, и из (3.60) получим выражение для поля температуры.

3.9. Теплопередача через ребристую стенку

Наличие ребер на стенке увеличивает поверхность ее соприкосновения с теплоносителем и тем самым уменьшает термическое сопротивление

32

теплопередаче. Поэтому наличие ребер может использоваться как средство интенсификации процесса теплопередачи или как средство снижения температуры стенки. Пусть на рис. 3.2 поверхность стенки со стороны жидкости с температурой Тж2 имеет в наличии ребра. При Тж1 Тж2 температура ребра, равная у его основания температуре поверхности между ребрами Тс2, будет уменьшаться к его концу. Температуру среды Тж2 можно считать неизменной для всей поверхности ребра. Поэтому участки поверхности ребра, удаленные от основания, будут отдавать меньше тепла, чем участки, расположенные вблизи основания ребра. Отношение количества тепла, передаваемого поверхностью ребер в окружающую среду Qp, к теплу, которое эта поверхность могла бы передать при постоянной температуре на поверхности ребер, равной температуре у их основания, Q p, называется коэффициентом эффективности ребер

ηp Q p .

Q p

Коэффициент эффективности ребер всегда меньше единицы. Чем резче меняется температура вдоль ребра, тем меньше коэффициент эффективности. Для коротких ребер, выполненных из высокотеплопроводных материалов, коэффициент эффективности близок к единице. Определим тепловой поток через плоскую стенку, площадь гладкой поверхности которой F1, а площадь ребристой поверхности – F2. Площадь F2 складывается из площади боковой поверхности ребер Fp и площади межреберных участков Fм. При стационарном режиме поток тепла от горячей среды к стенке, через стенку и от стенки к холодной среде выразится формулами

 

Q α1 (Tж1

Tc1) F1

(3.63)

 

Q

λ

(T

 

T

)F

(3.64)

 

 

 

 

 

δ

c1

 

c2

1

 

 

 

 

 

 

 

 

 

Q Qp Qм

 

(3.65)

Так как

 

 

 

ηp α2 (Tc2

Tж2 )Fp

Qp ηp Qp

Qм α 2 (Tc2 Tж2 )Fм ,

то, подставив это в (3.65) и исключив из (3.63), (3.64) и (3.65) температуры Тс1 и Тс2, получим

Q

 

 

 

(Tж1 Т ж2 )

 

 

 

 

.

(3.66)

 

1

 

δ

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α F

λF

α

2

(F

η

p

Fp)

 

 

 

 

 

 

 

 

1 1

 

1

 

 

м

 

 

 

 

 

Если проанализировать эту формулу при

α 2 α1, то можно обнаружить,

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

33

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

3.10. Температурное поле и коэффициент эффективности ребра постоянного поперечного сечения

Рассмотрим на рис. 3.5 ребро постоянного поперечного сечения, для которого изменением температуры по поперечному сечению ребра можно пренебречь по сравнению с изменениями температуры по длине ребра.

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент

теплоотдачи

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

θ0

 

 

 

 

 

 

 

 

 

 

температура

 

окружающей

среды

Тж

 

 

θ

 

 

 

θL

 

постоянны для всей

поверхности

ребра.

 

 

 

 

 

 

 

 

 

 

 

 

 

Площадь

 

поперечного

сечения ребра

f,

 

 

 

 

 

 

 

 

 

 

 

Tж

 

 

 

 

 

 

 

 

 

 

 

 

 

 

периметр

 

сечения

u,

длина

ребра

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

x

dx

 

 

 

коэффициент теплопроводности материала

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

также

 

постоянны.

При стационарном

 

 

 

Qx

 

 

 

 

 

 

Qx+dx

 

 

 

 

 

 

 

 

 

 

 

 

тепловом

 

режиме

поток

тепла

 

Qx,

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

 

 

 

dQ

 

 

 

подводимый в

сечении ребра

х частично

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

теплопроводностью передается вдоль ребра

 

 

 

 

 

 

 

Рис. 3.5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Qx+dx частично рассеивается в окружающую

 

 

 

 

 

 

 

 

 

 

 

 

 

 

среду dQ.

 

Q x

Q x dx

dQ .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(а)

 

 

 

 

 

 

 

 

 

 

 

Вводя избыточную температуру =T-Tж, согласно закону Фурье имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q x λ

d

f

 

 

 

 

 

(б)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q x dx

λ

 

 

 

 

 

 

dx f

 

 

 

 

(в)

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

Но по закону Ньютона-Рихмана

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ α u dx .

 

 

 

 

 

(г)

 

Подставив соотношения (б), (в) и (г) в формулу (а), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2

 

α u

m

2 ,

 

 

 

 

 

(3.67)

 

 

 

 

 

 

 

 

 

 

 

 

dx 2

λ f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где m

 

 

α u

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ f

Решением дифференциального уравнения (3.67) является функция

34

 

C

e mx

C

2

e mx .

(3.68)

 

1

 

 

 

 

 

Константы интегрирования С1 и С2

можно определить из

граничных

условий:

 

 

 

 

 

 

 

 

при х=0;

= 0,

 

 

 

 

 

 

x= ;

d

 

 

 

 

λ

 

 

 

α ,

 

 

 

 

 

dx x

 

 

 

 

 

 

 

 

 

где l, – избыточная температура и коэффициент теплоотдачи на торце

ребра.

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

при x= ;

d

0.

 

 

 

 

 

dx x

 

Определим константы интегрирования, пренебрегая теплоотдачей с торца ребра. Подстановка граничных условий в уравнение (3.68) дает

ddx

0 C1 C2 ,

mC1em mC 2 e m 0.

x

Из совместного решения этих уравнений определяются С1 и С2

C1

 

0

e m

 

 

;

 

 

 

 

C2

 

 

 

0

em

.

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2 ch( m )

 

ch (m )

 

 

 

 

 

 

 

 

 

Гиперболический косинус выражается формулой

 

 

 

 

ch (m )

 

em

e m

 

.

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив С1 и С2

 

 

в формулу (3.68), получим уравнение для

температурного поля в ребре

 

 

 

ch m x

 

 

 

 

 

 

 

 

0

 

.

 

 

 

 

(3.69)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch (m )

 

 

 

 

 

 

 

Избыточная температура на конце ребра при x = равна

 

 

 

 

 

 

 

 

 

 

 

0

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ch (m )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Весь рассеиваемый ребром поток тепла проходит путем теплопроводности

через сечение его основания. Поэтому

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

Q p λ

 

 

 

f.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx x 0

 

 

 

 

 

Продифференцировав (3.69) и подставив в эту формулу, получим

 

 

 

 

 

 

 

 

 

 

 

Qp 0

 

 

 

α u λ f th (m ).

 

 

(3.70)

Эта формула не учитывает теплоотдачу торца ребра. Если принять, что коэффициент теплоотдачи на торце ребра такой же, как на его боковой

35

поверхности, то теплоотдачу торца ребра можно учесть удлинением ребра на величину f/u и вместо длины ребра в расчетах использовать эффективную длину эф= +f/u. Тогда коэффициент эффективности ребра постоянного по длине произвольного поперечного сечения равен

ηp

Qp

 

0

α u λ f

th (m эф )

th (m эф )

(3.71)

 

 

 

 

.

Q1p

α 0 u эф

 

 

 

 

m эф

 

4. Нестационарная теплопроводность

4.1. Общие положения

При нестационарном процессе теплопроводности или теплопередачи температура в каждой точке тела изменяется во времени. Если при стационарном режиме поток тепла, подводимый к телу, равен потоку тепла отводимому от тела, то при нестационарном режиме переноса тепла эти потоки не равны. Разность этих потоков тепла в соответствии с законом сохранения энергии расходуется на изменение энтальпии тела, что связано с изменением температуры во времени. Поэтому, в отличие от задач стационарной теплопроводности или теплопередачи, при решении уравнений (3.1) или (3.2) задают еще и начальные условия. Т.е. для того чтобы рассчитать, каким будет температурное поле в теле для некоторого момента времени, необходимо задать, каким оно было в некоторый начальный момент времени.

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

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

36

времени, но ровно через период температура в любой точке тела принимает тоже значение.

Для определения температурного поля в телах классической формы (плоская неограниченная по длине и ширине стенка, неограниченный цилиндр, шар) используют аналитические методы решения. В телах сложной формы температурное поле определяют в результате численного интегрирования уравнений (3.1) и (3.2).

4.2. Охлаждение (нагревание) неограниченной пластины

Рассмотрим расчетные зависимости, полученные аналитическим методом [1,2,3], на примере охлаждения плоской стенки, размеры которой по длине и ширине настолько велики, что теплообменом с торцов можно пренебречь. При одинаковых условиях теплообмена на обеих поверхностях температурное поле будет одномерным и симметричным относительно середины стенки. Поэтому на рис. 4.1 ее толщина обозначена через 2 .

Т

 

 

 

 

 

 

 

 

Будем

полагать,

 

что

температура

 

 

окружающей пластину жидкости Тж=const и

 

 

коэффициент теплоотдачи = const. С учетом

 

Т0

принятых

ранее обозначений для

избыточных

 

 

температур

Т Т ж

дифференциальное

0

Тж

уравнение (3.2) будет иметь вид

 

 

 

х

а

2 .

 

 

 

 

(4.1)

 

 

 

τ

х 2

 

 

 

 

Начальные условия:

 

 

 

Рис. 4.1

 

при τ 0 0 Т0 Т ж .

 

(4.2)

При этом граничные условия в середине и на поверхности пластины соответственно можно записать так:

при х=0

 

 

0;

 

(4.3)

 

 

 

 

 

х

х 0

 

 

 

при х=

 

 

 

α

х δ .

 

 

 

 

(4.4)

 

 

 

х

х δ

λ

 

 

Решение этой задачи позволяет выразить безразмерную температуру в виде

ряда

 

 

 

 

n

 

 

2 sin μ n

 

 

 

 

 

 

 

 

 

 

θ

 

 

 

 

 

 

cos (μ n X) exp ( μ n2 Fo).

(4.5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n 1 μ n

sin μ n cos μ n

 

 

Здесь

θ

 

 

 

Т Тж

;

Х

х

;

 

Fo

 

a τ

число Фурье.

 

 

 

 

δ

 

 

δ2

 

 

 

0

 

 

Т Т

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Числа n определяются из трансцендентного уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

ctg μ

μ

,

 

(4.6)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Bi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

37

 

 

 

 

 

где Bi

α δ

число Био.

λ

 

 

Так как 1, 2,….., n представляет собой ряд возрастающих чисел, то чем больше , тем меньше роль последующего члена ряда в (4.5) по сравнению с предыдущим. Кроме того, чем больше число Fo, тем быстрее будут убывать члены знакопеременного (из-за cos n X) ряда (4.5) с увеличением номера n. Исследования показали, что уже при Fo 0,3 ряд (4.5) становится настолько быстросходящимся, что распределение температуры достаточно точно, с ошибкой не более 1 %, описывается первым членом ряда (4.5)

θ

 

2 sin μ1

cos (μ

 

X) exp ( μ 2

Fo).

(4.7)

 

 

1

 

 

μ1 sin μ1 cos μ1

 

 

1

 

 

 

 

 

 

 

 

 

 

Как видно из

(4.7), при фиксированной

безразмерной

координате Х

натуральный логарифм безразмерной температуры линейно зависит от безразмерного времени Fo. Для Х=0 и Х=1 решение (4.7) в [1] представлено в виде номограмм для различных значений чисел Био, так как 1 зависит в соответствии с (4.6) только от числа Био. Эти номограммы можно использовать в практических расчетах.

При Bi (практически при Bi 100) из (4.6) имеем ctg =0 и n

(2n 1)π , где n=1, 2, …… В этом случае температура поверхности пластины

2

сразу при Fo 0 становится равной температуре окружающей среды, а безразмерная температура при Х=1 =0. Иначе говоря, процесс отвода теплоты с поверхности тела происходит существенно интенсивнее, чем выравнивание температуры в теле. В этом случае процесс охлаждения (нагревания) определяется размерами тела и его физическими свойствами.

При Bi 0 (практически при Bi 0,1) ctg = и n =(n-1) . В этом случае все члены ряда (4.5), кроме первого, равны нулю. Но при 1 0

 

 

2 sin μ1

 

 

 

1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μ1

sin μ1 cos μ1 μ 0

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

С другой стороны, при 1

0 ctg μ1

cos μ1

 

 

1

 

μ1

;

μ12 Bi .

sin μ1

 

μ1

Bi

 

 

 

 

 

 

 

 

 

 

Учитывая это, уравнение (4.5) можно переписать в виде

 

θ cos

 

X exp Bi Fo .

 

 

 

(4.8)

Bi

 

 

 

 

 

 

cos

 

X 1

 

 

При малых числах

Bi

значение

Bi

и

температура на

поверхности пластины мало отличается от температуры в центре. Это указывает на то, что температура по толщине пластины практически распределяется равномерно. Иначе говоря, процесс выравнивания температуры в теле происходит существенно интенсивнее, чем отвод теплоты с поверхности тела, и задача охлаждения (нагревания) становится внешней.

38

Количество отданной (полученной) стенкой теплоты на любой момент времени 1 можно определить из следующих соображений. Количество отданной (полученной) стенкой теплоты с обеих ее поверхностей за время от=0 до = должно равняться (в изобарном процессе) изменению энтальпии пластины за период ее полного охлаждения (нагревания)

 

 

 

 

Qп 2δfρc (T0 Tж ) ,

 

 

 

 

 

 

(4.9)

где с – удельная теплоемкость;

 

 

 

 

 

 

 

f – площадь боковой поверхности;

 

 

 

 

 

 

 

– плотность.

 

 

 

 

 

 

 

Тогда за любой промежуток времени от =0

до = 1.,

или, что то же, от

Fo=0 до Fo= Fo1, энтальпия пластины изменится на

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

T

 

 

 

 

 

 

 

 

1

ж

 

 

 

 

 

 

 

 

 

Q Qп Q1 2δfρc (T0 Tж ) 1

 

 

 

 

,

 

 

 

 

 

 

 

 

 

T0 Tж

или

 

 

 

 

 

 

 

 

 

 

 

Q Qп (1

 

1 ),

 

 

 

 

 

 

(4.10)

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где θ1 (

Т

1 Т ж )/(Т0 Т ж ) средняя безразмерная

 

 

температура по толщине

пластины в момент времени 1.

Из соотношений (4.9) и (4.10) следует, что расчет количества теплоты, отданного (полученного) пластиной, сводится к нахождению средней безразмерной температуры в интересующий нас момент времени. Средняя безразмерная температура для слоя пластины от ее центра до плоскости Х равна

1 X

θ X 0 θ dX

в соответствии с теоремой о среднем.

Если в это выражение подставим под знак интеграла значение из уравнения (4.5) и проинтегрируем в пределах от нуля до единицы, то получим

 

 

 

 

2 sin

2

μ n

 

 

 

θ

 

 

 

exp ( μ n2 Fo).

(4.11)

 

μ n2

μ n sin μ n cos μ n

 

 

n 1

 

 

Подставив в уравнение (4.10) вычисленное по формуле (4.11) значение средней температуры пластины для интересующего нас момента времени Fo, получим количество теплоты, отданное пластиной в окружающую среду за этот промежуток времени.

4.3. Регулярный режим

Как показано выше, начиная с некоторого момента времени 1 охлаждения пластины, начальные условия начинают играть Второстепенную роль. Процесс полностью определяется только условиями охлаждения на границе пластины и среды, физическими свойствами материала тела и его геометрической формой и размерами. Математически это означает, что температурное поле в пластине описывается первым членом ряда (4.5)

39

 

 

 

0

2 sin μ

1

 

x

 

 

 

 

 

 

 

 

cos μ1

 

exp

μ12

 

 

A1 U1 exp m1 τ .

(4.12)

μ1

sin μ1 cos μ1

 

 

 

 

δ

 

 

 

δ 2

 

В этом уравнении А1-постоянный коэффициент, не зависящий ни от

координат, ни от времени, так как 1

определяется из соотношения (4.6).

Множитель

Un

является

функцией

только

координаты х. Для тел

других

геометрических форм температурное поле в стадии регулярного режима также будет описываться уравнением вида (4.12). Специфика геометрической формы учитывается различным видом множителей А1 и U1. Логарифмируя последнее уравнение и опуская индексы, получим

ln ln AU mτ .

(4.13)

После дифференцирования обеих частей уравнения (4.13) имеем

 

1

 

 

 

 

τ m const.

(4.14)

 

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

ln 1 ln 2

m const.

(4.15)

τ 2 τ1

 

 

Зависимость темпа охлаждения от физических свойств материала тела, его геометрической формы и размеров, а также условий теплообмена на поверхности тела можно найти из теплового баланса. При отводе от тела объемом V тепла dQ изменение энтальпии тела составит

 

 

 

 

 

 

dQ c ρ V

 

V

dτ ,

(4.16)

 

 

 

 

τ

 

где V – средняя по объему избыточная температура.

За тот же промежуток времени эта теплота должна быть отведена с

поверхности тела в окружающую среду за счет теплоотдачи

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dQ α F F dτ .

 

 

 

(4.17)

Приравнивая выражения (4.16) и (4.17), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

V

m

 

 

F

 

α F

Ψ

α F

.

(4.18)

 

 

 

 

τ

 

 

 

 

 

 

V

 

V

 

cVρ

 

C

 

Из этого уравнения следует, что относительная скорость охлаждения или темп охлаждения тела при конечном значении коэффициента теплоотдачи прямо пропорциональна коэффициенту теплоотдачи, поверхности тела и обратно пропорциональна его теплоемкости (первая теорема Кондратьева [1]).

 

 

 

 

 

 

 

В этом уравнении множитель Ψ

F

называется коэффициентом

 

 

 

 

 

V

 

 

40

 

 

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]