Файл: Вайсман М.Д. Режимы и способы пуска блоков сверхкритического давления учеб. пособие.pdf

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

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

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

Добавлен: 28.07.2024

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

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

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

нечное множество зависимостей для 0:

О, = C,J0 (k,r) е-°** + D, Y0 (k,r) е~ак- ;

0, = C2J0(А,г) е~ак-х-I- D,Yti (V ) в " ak* ;

Каждая

такая

зависимость, служит частным

интегралом

уравнения (32").

 

 

 

Общий интеграл уравнения теплопроводности строят по

принципу наложения — как

сумму бесконечного

множества

частных

интегралов:

 

 

 

0 = 2

ІСЛ а д

+ D J a( V ) ] <?“ °Ѵт-

(38)

 

п=1

 

 

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

Заметим, -что ни один-из частных интегралов сам по себе не удовлетворяет заданному начальному распределению тем­

пературы: при т =

0 каждое

Ѳ0і = С ^0(к{Г) + 1 \У0(/2;г)

есть

функция радиуса,

тогда как

по условию Ѳо = const.

Алге­

браическое суммирование множества частных решений позво­ ляет близко подойти к заданному начальному условию.

Перейдем к определению постоянных по краевым усло­ виям. Из граничного условия на внутренней поверхности имеем

CJ0{kRB) + D Y 0(kRB) = Q.

(39)

Условие іііа внешней поверхности дает

CJ0'\ k R tt) + DY0'(kRH) = 0.

В теории бесселевых функций доказывается [19], что произ­ водная от любой функции Бесселя порядка ѵ равна одно­ именной функции порядка ѵ+1, взятой с обратным знаком:

Л/ а д = - kJ>m « ) и к / а д = — к г, а д ,

(/і и Y, — бесселевы функции первого порядка, первого и второго рода).

Таким образом, из гранитного условия на внешней поверх­ ности цилиндра вытекает

су, а д + D r , а д = о.

(40)

Зависимости (39) и (40) есть линейные, относительно С и D однородные уравнения. Неизвестные таких уравнений могут

V 43


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

Л (АЛ„)

v 0 ( k R B) = 0.

 

(АЯ„)

Y\ (kR„)

 

Отсюда

 

 

Л № , ) Yi (bRn) -

h (kR„) T0 (АЯ.) = 0.

(4.1)

Из уравнения (41), называемого характеристическим уравне­ нием, определяются значения /г. Функции J и Y периодиче­ ские, и уравнение (41) имеет бесконечное множество веще­ ственных корней /г„, каждый из которых отвечает какомулибо частному решению вида (37).

В общем решении (38) число неизвестных постоянных можно уменьшить на единицу. Из (39) следует

 

 

г ) _ _/->

-А, (*/?„)

 

Подставив (38),

получим:

ги(М?вГ

 

 

 

ѳ =

у

С "

[ Л ( k nr ) Y n ( k n R B)

~

 

л*^1То [knRa\

 

- A ( W Y

ü(knr ) ] e - t,kY \

(38')

Для краткости обозначим

 

 

у ^ щ = Ап И Л ( Ѵ ) П ( ^ ^ в ) - Л ( А „ / ? в) Г 0( Ѵ ) = і / 0( Ѵ ) -

Тогда

 

0 = ^ А , Ѵ

о М

е - ак" \

(38")

 

п — 1

 

 

 

Постоянную

определим из

начальных

условий; при

т = 0 имеем

 

 

 

 

 

І

A „ V 0{ k nr ) .

(42)

 

я.-l

 

 

 

Поскольку распределение температур вначальный момент представлено рядом, членами которого являются комбина­ ции бесселевых функций, предполагается, что ряд можно ин­ тегрировать почленно. Умножим левую и правую части (42) на rV0(kmr)dr и проинтегрируем от RB до Яи. Так* как Ѳ0 = = const, имеем:

( r V 0{ k mr ) d r =

V А „

f r V 0( k nr ) V , ( k mr ) d r .

(43)

К

« - 1

к

 

44


Не приводя здесь доказательства, укажем, что при т ф п все интегралы правой части (43) обращаются в пуль; конечные значения интеграл имеет лишь при т — п. Тогда

R

R

 

% { г Л (V ) dr = Ап

f г vü\k„r) dr.

(430

К

«В

 

Из формул интегрирования функций Бесселя и произвольных' линейных комбинаций этих функций [19] следует:

? r V *(КП <*г = *\Я [ ИЛ V ) +-V7 (Кг)] =

к к

= - f [ ѵ0* (КЪ) + ѴгЧККд) - ^ [ Уо4KK*)+Vr(W]

(44)

Здесь

 

 

 

 

 

 

 

Ю (V ) = •/, (Кг) Уо ( К Ъ ) - J o (КЪ) у, (Kr)-

(45)

При r = /?H [см. (41)]

 

 

 

 

vx(KR») =

Ji (KRH) YO(КЪ) -

Jo (КЪ) К, (КЪ) = 0.

 

Кроме того, из определения функции Ѵ0(knr) имеем

 

Ѵ0 (КЪ) =

Л (КЪ) Yo (КЪ) -

Jo (КЪ) Y0(КЪ) = 0.

 

Таким

образом,

соотношение (44)

принимает вид

 

/

г Но2 (Kr) dr = 4

[Я,,21/0 (КЪ) -

RB21/, (КЪ)]

(440

«в

 

 

 

 

определения Ѵ0 (knr).

Преобразуем правую

часть (440-

Из

следует

Л (КЪ) = Jo (КЪ) у о (КЪ) - Jo (КЪ) Y0 (КЪ)-

Из характеристического уравнения (41) вытекает

Jp [kgRв) К {knRu)

Уо(КЪ) = (kt/Rn)

Подставив в выражение Ѵй(/г71/?„) значение Yo(KRo), за­ пишем

Ъ (КЪ) —.

Joі*„/?в)

(Jo ( К Ъ )

Ух (КЪ) Ji(knRn) Yo(knRii)].

JI (^л^?||1

Одно из свойств бесселевых функций выражается (см. напри­ мер [20]) соотношением

J\ (z)Y0( z )~ Jo (z )Y ,(z )^^ .

45


Следовательно,

 

~

Лі [fe/lR^

V o ( W

= -

 

Ji [fcfiRri)

Одновременно

 

 

 

^ (А Л ) = j, (А Л ) Yo ( W - J A W Гі (А Л ) =

 

2

(46)

Подстановка в (44') значений V0(knRn) и V\(knRв) дает

J г V* (k„r) dr = ^

.

(44")

VD

 

 

Интеграл в левой части (43'), согласно формулам интегриро­ вания функций Бесселя [19], выражается следующим об­ разом:

,

[(А Л )

V, (knRn) - (А Л ) I/, (А Л )]. (47)

I г Ѵа(V ) dr = ^

На основании равенств

 

 

1 /,(А Л )= 0

и і/,(А Л )= ;^ 5 ? 7

имеем

 

 

 

S

rV ü(k„r)dr =

«ft«2 ‘

 

 

 

Подставив значения интегралов в (43') и решая -относитель­ но Ап, получим:

Л21*«/?!,)

О,,.

Л" * У,* (АЯЛ„) —У«3 І*я/?в1

Внесем выражение Ап в (38"); тогда формула, описывающая в рассматриваемых условиях распределения по радиусу и времени избыточных температур стенок цилиндра, прини­ мает вид

О

____ -Л* (*„/?,.)_____

ЯА\,

(48)

^г (А/і^?ц)

•/(г [hnRв)

Ѵ'о ( V ) <?'

Обозначим 7?н/^в — ш и

/гл/?в =

ц„.

Тогда akn4 — р?і2 • атIRB2.

Как известно, критерий Фурье

Fo = ат//2 (/ — характер­

ный размер),

следовательно,

ц„2■ах/Р— pn2Fo.

Произведя

замену в (48), получим выражение для относительной избы­ точной температуры в такой форме:


П о

T(г,т)- Тс ___X 1

•/г(і^'и)

 

Оо ~

Т0- тс - пZé Jt (в,,/«)- Л2Ы

 

 

п = 1

 

 

 

х Ѵа( ѵ'п' щ ) е

(48)

 

 

Применим полученную формулу для расчета полей избы­ точных температур в стенке паропровода высокого давления блоков 300 МВт. Размеры поперечного сечения паропровода 0 245/45 мм; следовательно, RB= ll,b мм, m = Rj,jRB= 1,58. Коэффициент температуропроводности стали 15Х1М1Ф а »

яа12-10_6 м2/с, поэтому критерий

FO ^ 2 - 1 0 - 3T.

урав­

Значения

первых

семи корней

характеристического

нения

(41) указаны

в таблице:

 

 

 

 

 

п

1

2

ч 3

4

 

5

6

7

 

2,46

8,035

13,525

18,96

24,374

29,78

35,208

Пусть в паропровод поступает среда, вызывающая резкое охлаждение внутренней поверхности стенки: Тс</ Т0 и б ]> 0. Результаты расчета относительных избыточных температур

Рие. 10. Распределение по радиусу относитель­ ных інзбыточных температур

представлены кривыми на рис. 10 и 11. На рис. 10 приведено распределение б- вдоль радиуса в различные отрезки времени, отсчитываемые от момента мгновенного охлаждения внутрен-

і

47