Файл: Численное моделирование процессов энергоразделения в потоках сжимаемого газа.pdf

80 0.0 0.2 0.4 0.6 0.8 1.0
Λ
0
{Λ
0
????????
2 4
6 8
10
????
????
Re
˚˚
10 3
10 4
10 5
10 6
Рисунок 2.31 — Изменение коэффициента аналогии Рейнольдса для турбулентного пограничного слоя при наличии положительного продольного градиента давления потерям давления торможения. Следовательно, использование конвективного охлаждения в сверхзвуковом потоке при наличии продольного градиента дав­
ления представляется нецелесообразным с точки зрения увеличения давления торможения.
Испарительное охлаждение
Ещё одним способом охлаждения потока является испарительное охла­
ждение: впрыск в газовый поток испаряющейся жидкости.
В работе [
89
] Шапиро и Хоторном было предложено уравнение описы­
вающее изменение давления торможения в канале постоянного сечения при впрыске испаряющейся жидкости:

81
торможения за счёт испарительного охлаждения (термопрессия) было названо авторами работы [
89
] аэротермопрессором (АТП). Схема устройства и измене­
ние основных параметров показаны на рис.
2.32
Работа [
89
] послужила началом крупному теоретико-экспериментальному исследованию процессов испарения жидкости в высокоскоростном высокотем­
пературном потоке газа. Целью этого исследования было изучение условий повышения давления торможения. Серия экспериментов и расчётов была проведена в Массачусетском технологическом институте (МТИ) в период
1952–1956 гг [
90
].
Одно из первых экспериментальных исследований было проведено в рабо­
те [
91
]. Стоит отметить, что предварительные исследования показали, что
– эффективность АТП является функцией диаметра испарительной ча­
сти;
– АТП наиболее эффективно будет работать при сверхзвуковых скоро­
стях.
Другими словами, аэротермопрессор должен быть крупноразмерным и высокоскоростным устройством. Однако, несмотря на это в работе [
91
] был
Воздух
Вода
I
Сопротивление
II
Испарение
III
Трение
100 %
????
ω
????
M
????
˚
????
????
????
????
????
????
????
????
????
˚
????
????
????
????
????
Рисунок 2.32 — Качественное изменение параметров в аэротермопрессоре

82
исследован малоразмерный АТП (????
ℎ
« 50
мм, ????{????
ℎ
« 35
). Такой выбор пре­
следовал следующие цели:
– оперативно и при минимальных затратах исследовать процессы проис­
текающие в АТП и режимы его работы;
– получить экспериментальные данные для валидации теоретических рас­
чётов;
– путём экстраполяции данных оценить возможные характеристики ра­
боты крупноразмерного АТП;
– накопить информацию для создания среднеразмерного АТП.
Малоразмерный АТП состоял из конического сопла, системы впрыска воды, испарительной части (цилиндрической трубы постоянного сечения) диа­
метром ????
ℎ
“ 53.975
мм и длиной ???? “ 1.86 м (????{????
ℎ
“ 34.47
) и диффузора с углом раскрытия α “ 6
˝
и степенью расширения ???? “ 7.97.
В эксперименте измерялось статическое давление по длине испаритель­
ной части АТП. Так же использовались специально разработанные зонды для отбора пара.
В работе были проведены экспериментальные исследования 20 различных систем впрыска. Наилучшей была признана система осевого подвода воды с малой скоростью по ходу высокоскоростного газового потока. Система состояла из семи трубок с наружним диаметром ????
????????????
“ 3.175
мм и внутренним ????
????????
“
2.795
мм. Одна трубка располагалась в центре, а остальные шесть равномерно по окружности диаметром ???? “ 31.75 мм.
Влияние впрыска воды на распределение статического давления показано на рис.
2.33
. Измерения были проведены при приблизительно постоянном чис­
ле Маха на входе M
0
« 0.5
, но для трёх различных значений относительного расхода воды Ω
0
“
????
????
{
????
и двух значений температуры торможения на входе ????
˚
0
Кривые A0 и B0 (см. рис.
2.33
) соответствуют течению без впрыска воды при ????
˚
0
“ 21
и 560 °C, соответственно. Так как режим течения дозвуковой,
то статическое давление падает вдоль канала. Кривая A0 лежит выше кривой
B0, так как, во-первых, начальное число Маха M
0
для A0 немного ниже и,
во-вторых, вследствие меньшего числа Рейнольдса для «холодного» течения,
коэффициент трения для A0 также ниже.
Кривые A1 и A2 (см. рис.
2.33
) соответствуют распределениям давления при различных уровнях впрыска воды при низкой температуре. При этих режи­
мах происходит незначительное испарение жидкости и падение давления здесь

83
????
˚
0
,°C
M
0
Ω
0
‚
‚
‚
A0 21 0.465 0.000
‚
‚
‚
A1 21 0.475 0.104
‚
‚
‚
A2 21 0.475 0.201
‚
‚
‚
B0 560 0.500 0.000
‚
‚
‚
B1 560 0.485 0.105
‚
‚
‚
B2 560 0.495 0.209 0
5 10 15 20 25 30 35
????{????
ℎ
0.5 0.6 0.7 0.8 0.9
????{????
˚
0
Рисунок 2.33 — Влияние впрыска воды на распределение статического давления в канале постоянного сечения при наличии трения. Данные работы [
91
]
больше, так как сила лобового сопротивления капель добавляется к сопротив­
лению сил трения.
Сравнение кривых B1 и B2 с кривыми A1 и A2 (см. рис.
2.33
) показы­
вает, что впрыск жидкости при высокой начальной температуре приводит к снижению потери давления. При начальной температуре ????
˚
0
“ 560
°C испа­
рение жидкости весьма ощутимо и существенно снижает потери давления в сравнении со случаем без впрыска воды.
На рис.
2.34
показано изменение статического давления по длине кана­
ла при различных значениях начального относительного расхода воды Ω
0
. Все кривые соответствуют начальному числу Маха M
0
“ 0.63
и начальной темпера­
туре торможения ????
˚
0
“ 560
°C. Так как течение без впрыска воды при M
0
“ 0.63
соответствует критическому, то кривая для этого случая не показана.
Все кривые на рис.
2.32
и
2.34
можно условно разделить на следующие участки, в соответствии с протекающими процессами:
а) Резкое падение давления вследствие впрыска воды, т.к. на этом участке доминируют силы сопротивления капель.

84
????
˚
0
, °C ????
˚
0
, атм M
0 560 1.01 0.63 0
5 10 15 20 25 30 35
????{????
ℎ
0.55 0.60 0.65 0.70 0.75 0.80
????{????
˚
0
Ω
0 0.10 0.21 0.31 0.42
Рисунок 2.34 — Влияние относительного расхода воды на изменение статического давления по длине АТП. ????
ℎ
“ 53.975
мм; ????
˚
0
“ 560
°C;
????
˚
0
“ 1.01
атм; M
0
“ 0.63
. Сплошные линии — расчёт, штриховые линии — эксперимент [
91
]
б) Последующий быстрый рост давления, соответствующий быстрому ис­
парению жидкости, происходящему в «турбулентном» режиме, т.к.
относительная скорость (между каплями жидкости и газом) значи­
тельна. На этом участке эффект испарения доминирует над трением и сопротивлением капель.
в) После того как капли жидкости были разогнаны газовым потоком и относительная скорость стала близка к нулю, испарение перешло в «ламинарный» режим, т.к. коэффициент теплообмена существенно уменьшился, а кроме того существенно упал перепад температур меж­
ду каплями и газом. На данном этапе превалирующим становится сила трения и давление опять падает.
Совместный эффект сил трения, сопротивления капель воды и испарения при высоких температурах, как показывают результаты работы [
91
], может привести к непрерывному переходу от дозвукового потока к сверхзвуковому в канале постоянного сечения. Такой интересный феномен проиллюстрирован на рис.
2.35
. Все кривые на рисунке соответствуют ????
˚
0
“ 560
°C и Ω
0
“ 0.21

85
Стоит отметить, что кривые на рис.
2.35
сходны по своему поведению с рис.
2.34
В эксперименте число Маха не было измерено, однако из результатов рас­
чёта (см. рис.
2.35
б) видно, что на начальном участке число Маха возрастает,
когда сила сопротивления капель превалирует, затем, по мере того, как всё
большую роль играет испарение, число Маха уменьшается и наконец вновь возрастает с того момента, как сила трения становится главенствующей. Сле­
довательно, число Маха проходит через максимум когда испарение становится доминирующим. Если число Маха на входе постоянно увеличивать, то макси­
мальное число Маха тоже будет увеличиваться, пока не достигнет значения единицы. Такое «критическое» начальное число Маха, зафиксированное в экс­
перименте было равно M
0
????????
“ 0.78
. Значение массового расхода в этом случае достигает своего максимума и поток становится «запертым». Дальнейшее увели­
чение начального числа Маха за счёт изменения противодавления невозможно.
Для «запертого» течения (кривая, с M
0
“ 0.779
на рис.
2.35
) давление сначала падает, а число Маха возрастает, так как поток дозвуковой и сила сопротивления капель доминирует. Критические условия достигаются, когда число Маха становится равным единице, на расстоянии ????{????
ℎ
« 0.4
от сечения впрыска. В этих условиях испарение жидкости становится наиболее значимым,
давление продолжает падать, а число Маха продолжает возрастать, т.е. проис­
ходит непрерывный переход в сверхзвуковую область течения. На расстоянии
????{????
ℎ
« 10
силы трения начинают превалировать над всеми остальными эф­
фектами и, соответственно, в сверхзвуковом потоке давление начинает расти,
а число Маха снижаться. На расстоянии ????{????
ℎ
« 28
происходит скачок уплот­
нения и поток переходит в дозвуковую область течения. По-прежнему сила трения является основной, а значит давление начинает падать, а число Маха расти и достигает единицы в конце испарительной части. Течение в диффузоре становится сверхзвуковым, а затем скачком переходит в дозвуковую область.
Положение скачка уплотнения определяется противодавлением. По мере увели­
чения противодавления скачок уплотнения будет продвигаться вверх по потоку.
Характерные распределения давления вдоль канала показаны на рис.
2.35
а.
Для оценки эффективности работы АТП было предложено использовать следующий комплекс [
91
]:
????
˚
2
´ ????
˚
0
????
˚
0
M
2 0
,
(2.61)

86 0.3 0.4 0.5 0.6 0.7 0.8 0.9
????{????
˚
0
a)
Эксп. [
91
]
Расчёт
0 1
0.4 0.6 0.8 1.0 1.2
M
5 10 15 20 25 30 35
????{????
ℎ
M
0
“ 0.48
M
0
“ 0.63
M
0
“ 0.73
M
0
“ 0.779
M
????????????
0
“ 0.78
пере ло м
масштаб а
б)
Рисунок 2.35 — Влияние начального числа Маха на распределение статического давления (а) и числа Маха (б) по длине АТП. Непрерывный переход через скорость звука. ????
˚
0
“ 560
°C; ????
˚
0
“ 1.01
атм;
Ω
0
“ 0.21
.Сплошные линии — расчёт, штриховые линии — эксперимент [
91
]
где ????
˚
0
— давление торможения на входе в сопло, ????
˚
2
— давление торможения на выходе из диффузора и M
0
— начальное значение числа Маха, на входе в испа­
рительную секцию. Не трудно заметить, что показатель эффективности (
2.61
)
есть ни что иное как нормированная левая часть уравнения (
2.60
).
Значения коэффициента эффективности приведены на рис.
2.36
. На гра­
фике показаны максимальные значения коэффициента эффективности АТП,
наблюдаемые в эксперименте [
91
] при варьировании относительного расхода воды.

87 0.4 0.5 0.6 0.7 0.8
M
0
´0.6
´0.5
´0.4
´0.3
´0.2
´0.1 0.0
ˆ ????
˚
2
´ ????
˚
0
????
˚
0
M
2 0
˙
????????????
????
˚
0
,°C ????
˚
0
, атм
Ω
0 1200 0.98 ˜ 1.01 0.15 ˜ 0.18 1500 0.98 ˜ 1.01 0.23 ˜ 0.25 1800 0.98 ˜ 1.01 0.30 ˜ 0.32 1500 0.72 0.22 ˜ 0.24
Рисунок 2.36 — Эффективность малоразмерного аэротермопрессора. Данные работы [
91
]
Как и ожидалось, для малоразмерного АТП не удалось получить суммар­
ного роста давления торможения, так как малый размер аппарата влечёт за собой неблагоприятный эффект доминирующего трения. Тем не менее, ожида­
емое увеличение эффективности при росте числа Маха и уровня температуры подтвердилось.
В работе [
92
] исследовался среднеразмерный АТП состоящий из сужа­
ющегося сопла, испарительной части постоянного сечения, диаметром ????
ℎ
“
276.4
мм, длиной ???? “ 2.172 м и выходного диффузора с углом раскры­
тия 5
˝
. Максимальное увеличение давления торможения, зафиксированное в эксперименте p????
˚
2
´ ????
˚
0
q
{
????
˚
0
« 5 %
. Кроме того, исследовался канал с переменной площадью поперечного сечения. Площадь варьировалась специальной встав­
кой. Однако, дополнительные потери трения, вызванные вставкой, превысили положительный эффект изменения площади.
Экспериментальные исследования работы [
93
] проводились на той же установке, что и [
92
], однако геометрия испарительной части была модифици­