Файл: Смольский Б.М. Нестационарный теплообмен.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 22.07.2024

Просмотров: 139

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

 

1

■t,-

= А exp (— p? Foj)

 

 

 

h

 

 

 

 

 

1 - . . t-Ч .. Jo

А exp ( - t f Fon),

 

где'

^it

'La

 

 

 

 

 

 

 

 

 

 

 

2 (sin px — pt cos px)

r2 sin рх —

 

А =

тп

 

px — sin px COS p1

 

 

 

 

r\kx

 

Объединяя оба полученных выражения, находим

 

 

 

 

ig

^Жер

tl

 

 

 

 

Іжг*.

t II

(2.14)

 

М-х =

 

 

 

 

(Fon — Foi) lg e

 

 

Y

 

 

 

 

 

 

 

 

В выражение (2. 14) подставляем найденные из опы­ та значения температуры, измеренные в сплошном шаре на одном и том же расстоянии от центра для двух мо­ ментов времени. По найденному значению рі находим значение критерия Ві (рис. 8) [55], а по нему — коэф­ фициент теплообмена а.

Рис. 8. Корни характеристического уравнения tg р = —

р [55]

43


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

На рис. 9 приведены' результаты расчета коэффици­ ента теплообмена по первому способу для радиусов

/■'= 45 мм.

Рис. 9. Зависимость коэффициента теплообмена от времени, найден­ ная по первому способу: 1г' —45 лш; 2 — г"=35 мм

Второй способ. Здесь также делается предположение о постоянстве коэффициента теплообмена а.

Запишем выражение для температуры стенки шара для одного и того же произвольного момента времени Fo, но для двух значений радиуса г' и г" (соответственно tr' и tr ):

 

tr.

 

/•» Sin

p

 

t

=

В

г рх

exp (— Fo),

tn

 

 

 

^Жор

tf'

 

Г2 S in

р х

 

=

В

r Pi

exp (— p2Fo),

где

 

 

 

 

 

 

 

 

 

 

B =

2 (sin pi — px cos px)

 

 

 

Pl sin Pi COS Pi

44


Объединяя эти два выражения, получим

 

г

^ ---- — . (2.15)

J! sІП

г'

Жсо

tf'

В этом случае значения коэффициента ці находятся

путем графического решения

уравнения

(2.15). Затем,

используя рис. 8, находим значения критерия Ві и коэф­ фициент теплообмена.

Расчеты с помощью как первого, так и второго ме­ тода показали, что коэффициент теплообмена является в наших условиях переменной величиной, зависящей от времени. Он изменялся в 1,8—3,2 раза при использова­ нии первого метода для разных значений радиуса и в 3 раза при использовании второго метода. Полученные аб­ солютные значения коэффициента теплообмена нельзя было принять за действительные значения а, так как предложенная методика позволяет лишь оценить харак­ тер процесса теплообмена (квазистационарный или не­ стационарный). В случае квазистационарного процесса теплообмена метод дает действительное значение коэф­ фициента теплообмена.

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

3. МЕТОДЫ ЭКСПЕРИМЕНТАЛЬНОГО ОПРЕДЕЛЕНИЯ НЕСТАЦИОНАРНЫХ ГРАНИЧНЫХ УСЛОВИЙ

Определение теплового потока и коэффициента теп­ лообмена при малом значении критерия Ві. В случае стационарного режима теплообмена производные в при­ веденных формулах (2.17) и (2.18) являются посто­ янными величинами. При постоянном коэффициенте теплообмена из формулы (2.17) для шара как частный случай получается следующая зависимость для а:

а = И *

іп-І- .

Зт

t0

В этом случае зависимость ln t от т выражается прямой линией.

45

Для тонкостенных полых шаров, а также для корот­ ких цилиндров с большим коэффициентом теплопровод­ ности в работах [53, 105, 117] был использован метод определения граничных условий (qw, а), основанный на условии малости критерия Ві. В этом случае количество подведенного к телу тепла за интервал времени dx рав­ но теплоемкости тела, умноженной на изменение темпе­ ратуры тела dt. Температура в этом случае одинакова для всех точек тела. Таким образом, запишем

Qdx = GCpdt — а (т) (іЖооt) Fdx.

(2.16)

Из уравнения (2.16) можно получить следующие формулы коэффициента теплообмена:

для тонкостенного полого шара, цилиндра и пла­ стины

а (т) = — рсб d ln (/Жсо — t)

(2.17)

для сплошного шара

 

dx

 

рcR

d ln (іЖао — t)

 

а (т) = —

 

3

dx

 

 

 

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

Формулы теплового потока имеют вид: для вырезки из пластины

-

9 = р С6 —

(2.18)

для сплошного шара

dx

 

рcR

dt

 

_

'

^

3

dx

46

)


Этот метод был использован в работе [43] для опреде­ ления коэффициента теплообмена при поперечном об­ текании тонкостенного полого цилиндра. Температура стенки цилиндра в этой работе измерялась термометром сопротивления, а производная температуры по времени была заменена отношением конечных приращений. В ре­ зультате получаем конечную расчетную формулу

бср

АRcy

(2.19)

К «“ Кот

Ат

 

Аналогичная методика для полого цилиндра была ис­

пользована в работах [45, 46],

для полой сферы — в

работе [34], для сферических

частиц — в работах

[75,97]'.

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

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

Авторы описанных в литературе методов прибегают к некоторым приемам, упрощающим эту задачу.

Метод последовательных интервалов [24]. Заменив произвольную монотонную зависимость /теплового пото­ ка на поверхности пластины от времени ступенчатой кривой, авторы записывали решения задач для последо-

47


вательных интервалов времени. Для первого интервала времени формулировка задачи имеет вид

dt

=

дЧ

( 2. 20)

-----

а ------ ,

дт

 

дх2

 

dt(R,

т)

Д_=0

 

дх X

ді(0, т) ==0

дх

Температурное поле пластины при граничном условии 2-го рода для интервала времени 0 т хх определяется выражением [55]

Ѳ(х,

x) = t — 10

q-tR

I" пт

R2 — 3x2

 

 

 

6R2

 

 

 

 

 

 

 

+

V

( - i ) ,i+i

 

COS LI,,

' n R-

(2.21)

К

 

 

n~1

 

R

 

 

 

 

 

 

 

 

При Fo >

0,5

с погрешностью до 0,5%

можно пренебречь

членалш ряда.

Таким образом, для второго и последующих

интервалов при Fo>0,5 начальным условием является па­ раболическое распределение температуры по х:

Ѳ (х, tj) =

q1-

R

AFo, —

1 '

+

Ql X2.

(2.22)

 

 

X

 

 

6

 

2XR

 

Решением для п-го интервала будет

 

 

 

 

Ңх, хп) t0 =

^

г п~

] Яі Foi + Яп

^AFon

 

 

+

 

-

G

 

 

 

 

6

 

X*

\1

 

 

 

 

 

 

1

 

 

 

(2.23).

 

 

 

2R2 )_

 

 

 

 

Из решения для температуры на п-м

интервале времени

можно получить выражение для теплового потока

 

 

[t{x,

т) — *0]

 

І=П—1

 

 

------\ > , F

o ;

 

qn = ________________ _____ __________ .

(2.24)

 

F o

l

x

2

 

 

 

----- - +

2б2

 

 

 

 

 

6

 

 

 

48