книги из ГПНТБ / Зысина-Моложен, Л. М. Теплообмен в турбомашинах
.pdfпограничного слоя и капля имеет форму сфероида. Тогда, пре небрегая испарением капли, дифференциальное уравнение ее движения в пограничном слое можно записать в виде
|
|
|
т Ч |
™*кР" Ч |
- w " n f |
(VI. 19) |
|
|
|
dt |
|
|
|
где т — масса капли; dK— диаметр |
капли. |
|
||||
Примем |
сх = |
24/R, где |
|
|
|
|
|
|
|
D _ d K Ч -w 'n ) р ' |
|
||
|
|
|
|
ц" |
|
|
Решая |
(VI. 19), получим |
|
|
|
||
|
|
|
= 1 |
e~At/dl t |
(VI.20) |
|
где А = |
18ц 7р'; |
t — время полного торможения капли нормаль |
||||
ным к |
поверхности паровым |
потоком. |
|
Рис. 96. Торможение капли у поверхности нормальным паровым потоком
На рис. 96 приведены результаты расчетов, выполненных с помощью зависимости (VI.20), показывающие, как влияет диа метр капли на время t (в расчетах принималось, что начальная скорость капли равна половине средней скорости пара).
Перечисленные обстоятельства позволяют сделать вывод о том, что интенсивность сепарации на несмоченную стенку при прочих равных условиях ниже, чем на жидкую пленку, и должна зави сеть от дисперсности потока. Значительные трудности, связанные с получением монодисперсных потоков, объясняют отсутствие в настоящее время экспериментальных работ, посвященных деталь ному изучению этого вопроса.
На основании визуальных наблюдений отмечается что при температуре поверхности Tw ^ 220° С капли, выпадающие на стенку, принимают форму сфероида и отделяются от стенки устой-
232
чивой паровой пленкой. Интенсивность теплоотдачи при этом невысокая. Для оценки времени испарения сфероидальной капли, достигшей нагретой поверхности, может быть использована полу ченная в работах С. С. Кутателадзе и В. М. Боришанского зави симость
dl (?' — V") I0'5 |
= 2,5е-35.43 + 0,18, |
(VI.21) |
а |
|
|
где
Р= г X' ’
36.Теплоотдача в переходном режимес (Тw Тs) Г
Унос влаги из пленки в результате срыва жидкости с гребней
длинных волн при RnjI > Rnfl (см. п. 35), парообразование внутри пленки и интенсивное испарение влаги с поверхности создают предпосылки для механического разрушения пленки. Пленка вырождается в пульсирующие жгуты, площадь смоченной поверх ности постепенно уменьшается, пока не исчезнет полностью [54]. Переход к кризису в данном случае носит не резкий харак тер, изменение температуры происходит постепенно. Подобные режимы, наиболее характерные для условий, возникающих в про точной части и системах охлаждения турбомашин, являются одно временно -наиболее сложными для аналитического изучения. Здесь рассмотренные выше элементарные процессы чередуются во времени и находятся в сложном взаимодействии. Привести полную систему уравнений, описывающих данный процесс, не представляется возможным; речь может идти только о прибли женной постановке задачи.
Приняв допущение, что при малом влагосодержании основного потока процесс обтекания можно описывать уравнением однофаз ного потока — уравнением Навье-Стокса, получим при обычных упрощениях, соответствующих рассматриваемой задаче, в ка честве единственного критерия, характеризующего аэродинами ческую обстановку процесса, критерий Рейнольдса R " = w"b/v (здесь b — характерный размер обтекаемого тела или текущая координата вдоль обвода поверхности в направлении дви жения).
Если теперь абстрагироваться от сносящего действия основного потока и рассматривать процесс теплообмена между каплями, движущимися параллельно поверхности нагрева в пограничном слое, то изменение теплового потока в элементе пограничного
слоя можно записать в виде |
|
dQ = X " ^ - d yd F . |
(VI.22) |
233
Из данного уравнения при пренебрежимо малом торможении капель на участке dtу получается уравнение
Г |
(** - **)] dy, |
(VI.23) |
где GK— количество жидкости в данном элементе. Отсюда можно получить еще один определяющий критерий
I. — к |
(VI.24) |
К = с (Tw — Т s) |
Для случая, когда несущим взвешенные капли потоком яв ляется неконденсируемый газ, можно принять, что температура капель вне пограничного слоя будет равна температуре газа Т 0.
Если температуру и давление потока в процессе обтекания принять постоянными, то отношение (VI.24) может быть связано с температурным фактором (см. гл. III) следующим образом:
К = |
const |
С |
(VI.25) |
То (Ч>-1) |
тр.— Г ’ |
||
где С — некоторая константа, зависящая от |
температуры и |
||
давления потока. |
|
|
|
Пародинамическую сторону процесса можно характеризовать количеством влаги GBbIn, выпадающей на поверхность, и режимом
теплообмена при испарении этой влаги. |
На основании (VI. 1) |
в число определяющих процесс критериев кроме числа R следует |
|
ввести относительное влагосодержание |
|
G = |
<V I-26) |
Однако формула (VI. 1) не учитывает процесса испарения, про тиводействующего процессу выпадения влаги. Уравнение тепло
обмена при парообразовании позволяет, как известно |
[105], |
|||||||
получить |
еще |
два |
основных критерия: |
|
|
|
||
|
|
R* |
|
я |
Рг = — |
|
(VI.27) |
|
|
|
ry’v' |
|
|||||
Таким образом, интенсивность процесса теплообмена является |
||||||||
функцией |
следующих определяющих критериев: |
|
||||||
|
|
Nu = |
N |
(R", Pr", G, К, R*, у "If, |
/)• |
(VI.28) |
||
Здесь f — некоторый |
параметр, определяющий |
геометрию ка |
||||||
нала. |
J |
jrp |
|
|
|
|
|
|
Если |
принимаются постоянными, уравнение (VI.28) |
|||||||
|
и |
|||||||
примет вид |
|
|
|
|
|
|
||
|
|
Nu = |
~ |
= Z(R"t Pr", G, |
ф, ~q, f), |
|
(VI.29) |
234
где q — безразмерный тепловой поток, который может норми роваться параметром q0 = const.
Полагая, что влияние чисел R" и Рг" при малых значениях О будет таким же, как и в однофазном газовом потоке, представим
(VI.29) в виде
Nu == X (R, Рг) Y (G, ф, q, /). |
(VI.30) |
Введем коэффициент интенсификации теплоотдачи N и получим следующую форму записи для (VI.30):
|
|
|
|
|
|
|
^ |
= |
|
= |
|
|
Я* /)■ |
|
|
|
(у 1-31) |
||
Обратимся теперь к имеющимся экспериментальным данным. |
|||||||||||||||||||
В работах КАИ |
[50 ] |
был исследован процесс теплообмена между |
|||||||||||||||||
стенкой |
полой |
турбинной |
N |
|
|
|
|
|
|
|
|
|
|||||||
лопатки |
и |
охлаждающей |
Г ! |
|
|
|
|
|
|
|
|||||||||
жидкостью. Рабочим телом |
|
|
|
|
|
|
|
|
|||||||||||
служили воздушно-жидко |
|
|
|
|
|
|
|
|
|
|
|||||||||
стные |
смеси. |
Увеличение |
12 |
|
|
|
|
|
|
|
|
|
|||||||
коэффициента теплообмена |
|
|
|
|
|
|
|
|
|
||||||||||
при |
|
температуре |
стенки |
|
-а‘- |
|
|
|
|
|
|
|
|||||||
450—480° С было двукрат |
|
|
|
|
|
|
|
|
|||||||||||
ным для воды, а для керо |
|
I о |
|
|
|
|
|
|
|
||||||||||
сина |
составило |
|
5—10%. |
д и |
|
|
|
|
|
|
|
|
|||||||
Полученная |
|
обобщающая |
|
|
|
|
|
|
|
|
|
|
|||||||
зависимость для воздушно |
|
|
|
|
|
|
|
|
|
|
|||||||||
водяной |
смеси имеет |
вид: |
|
|
|
V |
|
|
|
|
|
|
|||||||
|
|
|
N = 1 + |
|
|
|
|
|
Т |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
. Аv |
|
|
|
|
|
||||||
I |
л |
*7г* '"'0,8 |
. —-1 |
65гл —0,08 |
|
|
|
|
|
|
|
|
|||||||
( |
0,750 |
лр |
|
|
Ксм |
, |
|
|
|
ч |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
(VI.32) |
|
|
|
|
|
|
|
|
|
|
||
где |
RCM— число Рейнольд |
о |
|
|
1Л |
|
|
|
18 |
|
|
||||||||
са для |
смеси. |
|
|
|
|
1,0 |
|
|
|
|
|
|
|||||||
В |
ЛПИ |
|
проводились |
Рис. 97. Зависимость |
интенсивности тепло |
||||||||||||||
исследования |
теплообмена |
отдачи от температуры поверхности для |
|||||||||||||||||
плоских |
гладкостенных |
и |
лопаточных |
радиаторов |
при |
G = |
8,55% и |
||||||||||||
оребренных |
|
радиаторов |
|
|
|
q = |
3,6- 10ь |
Вт/м2: |
|
||||||||||
для лопаток с двухкон |
Л — влажный |
пар; |
A, |
Ot |
-Ь — воздушно- |
||||||||||||||
водяная |
смесь (д |
и |
А |
— решетка |
радиаторов |
||||||||||||||
турной системой охлажде |
п р и / — 3,5; О — решетка |
радиаторов при / = |
|||||||||||||||||
ния |
|
[54]. |
Опыты прово |
== 2,6; # |
и + |
— одиночный |
радиатор |
соответст |
|||||||||||
|
|
|
венно при / |
= |
2,6 и / = 1) |
|
|||||||||||||
дились |
при |
|
давлении, |
|
|
|
|
|
|
|
смеси |
и влаж |
|||||||
близком к атмосферному, на воздушно-водяной |
|||||||||||||||||||
ном |
|
паре |
в |
диапазоне |
изменения |
режимных |
параметров: |
||||||||||||
q = |
(0,5-5-4,8) 10* Вт/м2; |
0 |
= 0-0,14; |
ф5 |
= |
T J T S = 1,05-2,5; |
|||||||||||||
R" = |
(3—10) 104. |
Коэффициент оребрения f составлял 2,6 и 3,5. |
|||||||||||||||||
Исследования |
|
сопровождались |
|
визуальными |
наблюдениями |
||||||||||||||
поверхности |
|
теплообмена. |
Обнаружена |
сильная |
зависимость |
235
интенсивности теплоотдачи N от режимных параметров. Наиболь
шее значение N = |
10-ь35 наблюдалось, как и следовало ожидать, |
в случаях, когда |
вся модель была покрыта жидкой пленкой |
(рис. 97). Область существования этих режимов ограничивается значением г|з5 1,15. При увеличении тепловой нагрузки пленка жидкости разрывается на отдельные жгуты, которые с видимой при визуальном наблюдении частотой перемещаются по поверх ности нагрева в направлении, перпендикулярном течению потока.
На участках поверхности модели, находящихся в аэродинами
ческой тени, |
вместо жгутов |
образовывались мокрые пятна, мед |
|||||||||||||||
|
|
|
|
|
|
|
ленно |
изменяющиеся |
очер |
||||||||
|
|
|
|
|
|
|
тания * которых |
позволяют |
|||||||||
|
|
|
|
|
|
|
предположить, |
что |
характер |
||||||||
|
|
|
|
|
|
|
процесса здесь не такой, как |
||||||||||
|
|
|
|
|
|
|
в жгутах |
|
жидкости. |
|
Кипе |
||||||
|
|
|
|
|
|
|
ния жидкости в пятнах не |
||||||||||
|
|
|
|
|
|
|
наблюдалось. Область сущест |
||||||||||
|
|
|
|
|
|
|
вования |
подобного |
режима |
||||||||
|
|
|
|
|
|
|
лежит |
в |
интервале |
1,15 < |
|||||||
|
|
|
|
|
|
|
< |
< |
2,0 |
(рис. |
97). |
Зна |
|||||
|
|
|
|
|
|
|
чения |
коэффициента |
интен |
||||||||
Рис. 98. Влияние концентрации |
влаги |
сификации |
теплоотдачи |
по |
|||||||||||||
мере |
возрастания |
ф5 |
умень |
||||||||||||||
в потоке |
на теплоотдачу |
лопаточных |
шаются, |
и при |
|
я» 2 имеем |
|||||||||||
радиаторов |
при |
a|)s = |
2 |
(кривая |
|
1) и |
|
||||||||||
t|’s — 1,3 (кривые 2—5) и следующих |
зна |
N —>I. |
|
Отметим, |
что |
|
при |
||||||||||
чениях q (в |
Вт/ма): |
|
|
определении N |
теплоотдача |
||||||||||||
/ — 4,0- 10s; |
2 — |
1,9-10“; |
3 — 3,0-Ю 5; |
4 — |
|||||||||||||
как |
в опытах |
на |
влажном |
||||||||||||||
— 4,0- 10s; 5 — расчет по |
формуле (VI.32) |
||||||||||||||||
|
|
|
|
|
|
|
паре, |
так |
и |
в опытах |
на |
||||||
воздушно-водяной смеси относилась к коэффициенту |
теплообмена |
для сухого насыщенного пара, результаты с достаточной точностью совпали. Этим подтверждается выдвинутое выше предположение о наличии двухфазного (жидкостно-парового) пограничного слоя для общего случая газожидкостного потока.
При > 2 жидкости на поверхности не наблюдалось и N =
= 1,0-ь 1,2.
Во всей области переходного режима теплоотдача примерно линейно зависела от удельного теплового потока. Увеличение
влажности потока в интервале 0 < G <С 6% приводило к увели чению теплоотдачи, при больших значениях G теплообмен стаби лизировался (рис. 98).
Интересно отметить особенность переходного режима на оребренных поверхностях. Здесь пленка при ^ 1,2 полностью исчезала на межреберной части стенки. Движение жидкости при больших значениях % осуществлялось только по поверхности ребер, т. е. пленка оттеснялась в область меньших тепловых потоков. Этим, очевидно, объясняется то обстоятельство, что
236
значения коэффициента интенсификации для оребренных поверх ностей ниже, чем для гладких.
Экспериментальные данные по средней и локальной теплоот
даче |
с точностью ±20% описываются корреляцией: |
|
||
|
|
N = 1 + |
Yq 0°'27я|Г4,65 |
(VI.33) |
|
|
Y = 92,5f~l-\ |
||
|
|
|
||
где / — коэффициент |
оребрения (/ = 1,0 ±3,5). |
|
||
В |
практических |
расчетах |
применение соотношения (VI.33) |
связано с необходимостью искать один из определяющих крите риев косвенным путем.
Заслуживают внимания данные работы [1 ], где получено
следующее выражение |
для критических тепловых нагрузок в труб |
|||||
ках диаметром d = |
5* |
10~3 м щ длиной/ — 0,240±0,425 м при G = |
||||
= 0,1 ±0,5 |
и р ---- |
(1,96 ±9,8) |
10БПа: |
|
|
|
9кр |
Дi |
@пл + Gb |
exp |
0.3164R'—°-25дг,кр |
. (VI.34) |
|
л dx,кр |
2d |
Здесь At — изменение энтальпии при парообразовании; хкр — значение координаты, отсчитываемой в направлении потока, при которой наступает кризис; GB= G' — Опл; GnjI — расход жидкости в пленке для сечения на входе в канал.
Явления кризиса, изучавшиеся в работе [1], носили явно выраженный характер. Там же предлагается корреляция для определения расхода жидкости в пленке
: 0,985 |
0,44 lg [ - b ( J i ^ L ) 2 Ш4]. |
(VI.35) |
Формулу (VI.35) можно рекомендовать для практического использования при расчетах, например, трактов охлаждения лопаток. Однако распространение данной зависимости на более короткие трубы не рекомендуется. Необходимы специальные исследования областей существования переходного режима и режима с явно выраженным кризисом.
37. Теплообмен при течении влажнопаровых потоков в межцилиндровых пространствах корпусов паровых турбин
Приведенные в п. 33 данные о теплообмене в моделях полостей корпусов паровых турбин могут быть использованы при анализе течения перегретого пара, когда последний можно рассматривать как однородный газ. Такое рассмотрение допустимо примени тельно к условиям ЦВД и ЦСД мощных паровых турбин, однако при исследовании течения пара в полостях корпуса ЦНД этих турбин нельзя абстрагироваться от влажности пара и дисперсной структуры влажнопарового потока. Аналогичное явление имеет
237
место и при течении в условиях влажнопаровых турбин АЭС. Эти турбины благодаря сравнительно низкому рабочему давлению обладают известными преимуществами перед турбинами сверх критических параметров и по мнению ряда специалистов являются предпочтительными при работе на переменных режимах [141].
Для выработки обоснованных рациональных графиков пуска и остановки машины необходима правильная оценка температур ных напряжений и тепловых расширений ее элементов. В усло виях переменных режимов работы результаты расчетов темпера турных полей в элементах ротора и статора турбины будут су щественно зависеть от точности задания граничных условий теплообмена.
В потоке влажного пара, как уже отмечалось выше, интенсив ность теплообмена в значительной мере определяется характером контактирования влаги с поверхностью и характером ее выпадения на поверхность.
В настоящее время практически отсутствуют эксперименталь ные и теоретические данные по теплообмену во влажнопаровых потоках, которые могли бы быть использованы для определения коэффициентов теплоотдачи в межцилиндровых полостях корпусов современных влажнопаровых турбин. Пограничный слой во влаж нопаровых потоках не изучен, имеются только первые опытные исследования, которые относятся к течению без теплообмена [38]. По этим данным режим течения пленки (ламинарный, волновой, турбулентный) существенно зависит от сепарации влаги и скорости течения пара. Наличие теплообмена, как уже говорилось выше, еще более усложняет задачу, при этом возможно возникновение большого количества различных режимов течения, которые вкратце были описаны в предыдущих параграфах. Для межцилиндровых полостей современных влажнопаровых турбин из всех возможных режимов течения интерес представляют дисперсный, когда пленка на поверхности не возникает, и дисперсно-кольцевой, когда на поверхности имеется пленка, а в остальной части — дисперсный поток пара. При пусковых режимах, когда температура стенки ниже температуры насыщения, возникает дисперсно-кольцевой режим течения. При сбросах нагрузки, сопровождающихся резким снижением параметров пара, возможно возникновение таких режимов, когда пленка конденсата вся испарится или будет отде лена от поверхности паровой прослойкой. В этом случае интенсив ность теплоотдачи может значительно уменьшиться. При боль ших скоростях движения пара и при наличии на стенке тонкой пленки может произойти вышеописанное подавление кипения, когда теплообмен будет осуществляться только за счет конвек тивного переноса в пленке, а испарение происходить только с по верхности пленки.
Обширный экспериментальный материал, имеющийся в ли тературе по теплообмену во влажнопаровых потоках примени тельно к задачам парогенераторстроения, не может быть исполь
238
зован для рассматриваемых задач вследствие значительно боль ших скоростей движения пара и существенно меньших степеней влажности, характерных для турбин.
Выполненная в ЦКТИ работа по исследованию теплоотдачи пластины, обтекаемой влажнопаровым потоком, в условиях, близких к условиям течения пара в межцилиндровых простран ствах влажнопаровых ступеней турбин, показала, что в этом случае возникают другие закономерности. Опыты проводились с плоской пластиной, находящейся в пародинамической трубе, и с пластиной, являющейся одной из стенок плоской модели, геометрически подобной продольному сечению выхлопной части ЦВД турбины АЭС. Исследования проводились при следующем изменении параметров: влажность у = 0 —30%; давление р 0 =
=(0,5—2,0) 105 Па; R" = (1,5-20) 10»; К = 80-2700.
Режимы течения в опытах соответствовали пусковым, темпе
ратура стенки Tw всегда была ниже температуры насыщения Ts. Средние коэффициенты теплоотдачи в зависимости от режимных параметров колебались в широком диапазоне: а ^ (2—200) X X 103 Вт/(м2-К). Для оценки влияния скорости движения пара на интенсивность теплоотдачи обычно вводится параметр N = = а/ан, где ан — средний коэффициент теплоотдачи в неподвиж
ном паре. |
Значение а н можно определить с помощью приведенного |
в работе |
[131 1 соотношения |
з„=о,72 (vi.зб)
Для медленно движущегося пара в работе [1051 приводится фор
мула, позволяющая определить N при турбулентном |
движении |
|
в пленке: |
|
|
Л 7 - 1 + 0 ,0 1 3 ( £ ) М ^ |
. |
(VI.37) |
Сопоставление опытных данных ЦКТИ со значениями, вы численными по этим формулам, показывает, что параметры, вхо дящие в формулы (VI.36) и (VI.37), не учитывают всех особен ностей теплообмена в быстро движущемся влажнопаровом по токе. Как видно из рис. 99, расчетные данные (сплошная линия), полученные по формуле (VI.37), отличаются от опытных данных в несколько раз, причем в разные стороны. Анализ этих опытных данных позволил обнаружить следующую зависимость пара
метра N от чисел Рейнольдса Rlny" d/v" и К = г/(сАТ):
N = -¥-=■ I 0,0018/С ( ^ y ,5(Rl)°'4. |
(VI.38) |
На рис. 100 приведено сопоставление значений, вычисленных по этой формуле (сплошная линия), с теми же опытными данными,
239
N
10 20 JO SO 100 2 0 0 / у " ) 0-5 w"
\ f ) ((fv')'l3
Рис. 99. Сопоставление опытных данных по теплоотдаче с расчетными, полу' ченными для медленно движущегося пара:
• — у = О (сухой насыщенный пар); О — У = 1 н-30% (влажный пар)
Рис. 100. Зависимость N от определяющих параметров:
Ф - У = 0 ; 0 - У = 1 +30%
которые приведены на рис. 99. Как видно, в этом случае расхож дение экспериментальных и расчетных данных практически не превышает ±25%, причем опытные точки для различных давле ний и различных влажностей практически совпадают.
Эти данные могут быть использованы для оценки теплоотдачи в камерах при пусковых режимах, однако в натурных условиях процессы будут более сложными: в криволинейных каналах при наличии градиента давления будет возникать дополнитель ная сепарация влаги, возможно возникновение отрывов погра ничного слоя, поток может быть пространственным и т. д. Все эти явления скажутся на величине коэффициентов теплоотдачи. При режимах остановки турбины характер теплообмена будет отличаться от вышеописанного, так как в этом случае темпера тура стенки будет выше температуры потока и вместо конден сации возможно возникновение испарения или кипения. Эти процессы в условиях турбомашин в настоящее время не иссле дованы.
16 Л. М. Зысина-Моложен ДР •