Файл: Вайсман М.Д. Режимы и способы пуска блоков сверхкритического давления учеб. пособие.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, имеем:
6« ( 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 |
u° |
Внесем выражение Ап в (38"); тогда формула, описывающая в рассматриваемых условиях распределения по радиусу и времени избыточных температур стенок цилиндра, прини мает вид
О |
____ -Л* (*„/?,.)_____ |
ЯА\, |
(48) |
||
^г (А/і^?ц) |
•/(г [hnRв) |
Ѵ'о ( V ) <?' |
|||
Обозначим 7?н/^в — ш и |
/гл/?в = |
ц„. |
Тогда akn4 — р?і2 • атIRB2. |
||
Как известно, критерий Фурье |
Fo = ат//2 (/ — характер |
||||
ный размер), |
следовательно, |
ц„2■ах/Р— pn2Fo. |
Произведя |
замену в (48), получим выражение для относительной избы точной температуры в такой форме:
4б
П о |
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 |