Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Тепловые и гидродинамические процессы в колеблющихся потоках

..pdf
Скачиваний:
9
Добавлен:
12.11.2023
Размер:
13.64 Mб
Скачать

Б. М. ГАЛИЦЕЙСКИЙ.

Ю.А. РЫЖОВ,

Е.В. ЯКУШ

ТЕПЛОВЫЕ И ГИДРОДИНА­ МИЧЕСКИЕ ПРОЦЕССЫ В КОЛЕБЛЮ­

ЩИХСЯ

ПОТОКАХ

Под редакцией профессора В. К . КОШ КИНА

М о с к в а

«М А Ш И Н О С ТРО ЕН И Е»

УДК 621.1.0167

Рецензент д-р техн. наук Н. И. Мелик-Пашаев

Галицейский Б. М., Рыжов Ю. А., Якуш Е. В. Тепловые и гидродинамические процессы в колеблющихся потоках. М., «Машиностроение», 1977. 256 с.

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

Книга предназначена для научных работников и инженеров, конструкторов и расчетчиков научно-исследовательских и проектно-конструкторских организаций.

Табл. 4, ил. 137, список лит. 75 назв.

20406—018

18—77

 

Г 038 (01)—77"

© Издательство «Машиностроение», 1977 г.

ПРЕДИСЛОВИЕ

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

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

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

исходит его разрушение.

течения

я.

Неустойчивый процесс

газового потока возникает

не только в жидкостных

ракетных

двигателях. Подобные про­

цессы возможны в воздушно-реактивных и в плазменных дви­

гателях,

а

также

в магнитогидродинамических генераторах,

в ядерных силовых и энергетических установках.

 

энерге­

При проектировании и анализе надежности различных

тических

установок

важно знать, как под действием

колебаний

в них изменяются тепловые процессы.

колебаний

Другая

важная

проблема, связанная с влиянием

жидкости на теплообмен, возникает при проектировании

и раз­

работке топливных баков и систем подачи рабочего тела

в лета­

тельных аппаратах. Как известно, маневрирование летательного

аппарата (торможение, разгон, разворот по курсу)

вызывает

i*

з

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

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

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

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

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

Настоящая монография является первым опытом обобщения методов экспериментального и теоретического анализа тепло­

обмена и

гидродинамики в ^колеблющихся

потоках.

Гл. I,

II

написаны

Б. М. Галицейским, гл. III и V—Б. М. Га-

лицейским

и Ю. А.

Рыжовым, гл. IV

и разд. 1 гл. V —

Е. В. Якушем.

 

 

Авторы благодарны рецензенту проф. Н. И. Мелик-Пашаеву за ряд ценных указаний и критику. Все замечания по книге будут с благодарностью приняты авторами. Просьба присылать заме­ чания по адресу: Москва, Б-78, 1-й Басманный пер., д. 3, изд-во «Машиностроение».

УСЛОВНЫЕ ОБОЗНАЧЕНИЯ

ДА — амплитуда колебания тела или его поверхности;

а— скорость звука, коэффициент температуропроводности;

ср, су — теплоемкость при постоян­ ном давлении и объеме;

d — диаметр;

Е — кинетическая энергия; d3 — эквивалентный диаметр; F — площадь;

/ — частота колебаний;

g — ускорение

силы

тяжести;

i — энтальпия;

 

 

k — — -----показатель

адиабаты,

су

 

 

волновое число;

 

L, / — длина;

трения;

/тр — работа сил

т — коэффициент потерь;

п— номер резонансной гармо­ ники;

р— давление;

q — плотность теплового потока; qy — мощность внутренних источ­

ников тепла;

R — газовая постоянная; S — энтропия;

Т — температура; Т* — температура торможения;

t — время;

и— скорость потока, внутрен­ няя энергия;

V — объем;

v — поперечная скорость потока;

v* — динамическая

скорость;

W —- безразмерная

скорость по­

тока,

фазовая скорость;

х, у %г — координаты;

импеданс;

Z — акустический

а — коэффициент теплоотдачи;

Р — коэффициент

ослабления,

коэффициент

объемного рас­

ширения;

 

толщина

6 — толщина

стенки,

пограничного

слоя;

 

8 — безразмерная

амплитуда ко­

лебания

скорости,

безраз­

мерный

коэффициент турбу­

лентной

вязкости;

 

ц — коэффициент

объемной вяз­

кости;

 

 

температура;

0 — безразмерная

Л — длина

волны

колебаний;

X — коэффициент

теплопровод­

ности;

 

 

турбулент­

|ЯТ — коэффициент

ной

теплопроводности;

р, — коэффициент динамической

вязкости;

турбулентной

р,т — коэффициент

вязкости;

сопротивления

| — коэффициент

трению;

 

кинематиче­

v — коэффициент

 

ской вязкости;

 

vT — коэффициент

турбулентной

кинематической

вязкости;

П — периметр;

 

 

р — плотность;

 

 

Т — период колебаний;

т — касательное

напряжение;

Ф — сила трения,

отнесенная

кединице поверхности;

Ф— фаза колебаний;

со — круговая частота колебаний;

Bi =

-----критерий

Био;

 

Fo

-----критерий

Фурье;

Gr _

«РАП8

— критерий Грас-

 

V 2

 

 

 

 

Но-

гофа;

 

гомохронности;

■критерий

J

АА ^ о 2! 3 — критерий

ви­

 

брационного ускорения;

k =

----- относительный коэф­

 

фициент

теплоотдачи;

k\ =

i------относительный коэф-

 

ьо

 

сопротивления тре­

 

.фициент

 

ния;

 

 

 

 

М0 =

— -----критерий Маха-Маев-

 

а0

 

 

 

 

 

ского; _

*

 

 

Мд =

Ди

 

 

число

----------колебательное

 

ао

 

 

 

 

 

Маха;

 

 

 

 

Nu =

X — критерий

Нуссельта;

Ре =

аI

критерий

Пекле;

Рг ■

критерий

Прандтля;

 

Re:

критерий Рейнольдса;

 

Ш

 

 

 

 

 

—------- амплитудный крите­

 

рий Рейнольдса;

.

ReA:

ДиА

 

 

 

----------- амплитудно-колеба-

5

 

тельный

критерий

Рей­

 

нольдса;

 

 

Re^ =

----------

колебательный

кри­

 

терий

Рейнольдса;

 

Re*, =

----------Аи2

критерий Рейнольдса

 

вторичных течений;

 

Sh =

-----критерий Струхаля.

Символы: А — пульсационные со­ ставляющие; <з > — осреднение по периоду колебаний, черта сверху — осреднение по длине канала.

Индексы: оо — параметры на бес­ конечности пограничного слоя; f — параметры, средние по сечению канала;

W — параметры

на поверхности; 0 —

амплитудные

значения

параметров,

осредненные

по

времени

значения

критериев;

i — проекция

на

ось

Xi (xt у у

z);

s — параметры,

соответ­

ствующие

резонансным

колебаниям;

(*) — безразмерные параметры,

пара­

метры торможения; штрихи «'» — зна­ чения турбулентных пульсационных составляющих параметров потока; «(/)» — порядок приближенного реше­ ния.

Г Л А В А I

ОСОБЕННОСТИ ИССЛЕДОВАНИЯ ТЕПЛОВЫХ

ИГИДРОДИНАМИЧЕСКИХ ПРОЦЕССОВ

ВКОЛЕБЛЮЩИХСЯ ПОТОКАХ

1.КЛАССИФИКАЦИЯ КОЛЕБЛЮЩИХСЯ ПОТОКОВ

ИИХ ФОРМАЛЬНОЕ МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ

В практике встречаются различные колебания веще­ ственной среды, которые могут быть обусловлены как гидродина­

мическими, так

и тепловыми возмущениями.

К гидродинамическим колебаниям относятся колебания ско­

рости, давления,

плотности среды или

колебания самого тела;

к тепловым — колебания температуры

среды, температуры по­

верхности тела, теплового потока на поверхности тела или тепло­ выделения в среде.

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

Рассмотрим периодическое изменение произвольного пара­ метра (А) среды или тела с амплитудой колебаний ДА 0, круго­ вой частотой со и фазой колебаний <р. В случае простейшего си­ нусоидального колебания изменение рассматриваемой величины А

во времени определяется из выражения

 

А = ДА0sin (со< — <р),

(1)

где t — время.

Промежуток времени Т = 2зт , в течение которого параметр А

совершает одно колебание, называется периодом гармонических

колебаний,

а число полных колебаний в единицу времени / =

= Tjr =

называется частотой гармонических колебаний.

В общем случае фаза колебания <р может зависеть от простран­ ственных координат х и х 2, х3. В частном случае, если q> = ± k x (где k — постоянное число), то уравнение (1) описывает плоскую бегущую недеформированную волну. Моментальный снимок такой волны представляет синусоиду, например А = AA0 sin(Ax) при t = 0. Коэффициент пропорциональности между фазой и расстоя­ нием х называется волновым числом. При положительном k волна распространяется в сторону возрастающих значений х.

7

Пространственный период, т. е. такое расстояние Л, при котором sin [k (х + Л) ] = sin (kx), называется длиной волны. Длина волны связана с волновым числом соотношением kA = 2п и характеризует пространственную периодичность волны, подобно тому как период колебаний Т определяет ее временную периодич­ ность.

В большинстве случаев гармонические колебания удобно представлять в комплексной форме

А = ДЛ*ехр (to/) = ДА0ехр — t<p);

(2)

здесь ДЛ* = а + ib — комплексная амплитуда. При такой форме записи амплитуда и фаза колебаний будут соответственно равны

ДЛ0= | ДЛ | = |/ а 2+

Ь \ <р =

arg (ДЛ) =

arctg -j- = kx.

(3)

При распространении

волны в

реальных

средах амплитуда

ее колебаний под действием диссипативных эффектов (например, преобразование кинетической энергии под действием трения

втепло) будет изменяться как во времени, так и в пространстве.

Вэтом случае волновое число является величиной комплексной:

k = k 0 — фх.

(4)

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

В этом случае волна будет описываться комплексным выраже­

нием

 

 

 

 

 

Л = ДЛ0ехр (— Рд-х) ехр (— ik9x) ехр (Ш)

(5)

или зависимостью относительно действительных переменных

 

Л = ДЛ0ехр (— рх) sin (ait — ф).

(6)

Физически коэффициент ослабления р характеризует ослабле­

ние

амплитуды колебания волны при распространении

волны

на

единицу пути х.

 

 

 

 

 

Таким образом, в общем случае как амплитуда, так и фаза

зависят от характерных координат.

В точке (х + 6x4

t + bt)

с точностью до членов первого

порядка малости фаза

 

 

iff — i |ф (х) +

бх] + mt + тЫ.

 

 

Она совпадает в точке (х,

t) с исходной фазой (конечно, ампли­

туда при этом изменяется),

если

 

 

 

 

бх +

Ш =

0.

 

8

Следовательно, отсюда можно определить фазовую скорость (или скорость распространения волны)

W — =

со:

0)

(7)

дх k

т

St

 

В практике, как правило, колебания отличаются от синусои­ дальных и носят более сложный характер. В этом случае для формального математического описания периодических колеба­ ний используется их спектральное разложение, основанное на рядах Фурье. Согласно методу Фурье периодическую функцию f (0 периода Т можно разложить в ряд по отдельным гармоникам

п k—n

f(t)= -тр + ^

[a* cos (Ш) +

bksin (Ш)] =

^

ckexp (ikat),

k= 1

 

 

 

k= —П

 

 

 

 

 

 

(8)

где коэффициенты

ряда

Фурье

соответственно

равны:

 

т

 

 

 

 

 

ak = ^

1 f(t) cos {ha>t)dt {k = 0,

1, . .

n)\

 

о

 

 

 

 

 

 

т

 

 

 

 

 

bk = Y

\

sin(A»0 dt (k = 0,

1, . . ., n);

 

0

 

 

 

 

 

T

 

 

 

 

 

 

cA= Y " J / (т)ехр (— 1Ш) dt (k = 0, ± 1, . . . , ± n)\

о

 

 

 

 

 

 

Ck = Y (ak— ibk)\

c_k = Y (ak +

ibk).

Амплитуда и фаза 6-й гармоники соответственно будут равны

 

АЛ0* = | с* |;

Ф* = arg (ck).

 

 

Таким образом, колебания любого параметра среды характе­ ризуются амплитудой ДА о (зависящей от коэффициента ослабле­ ния Р), частотой со, длиной волны А или скоростью ее распро­ странения W (поскольку длина волны А пропорциональна ско-

рости распространения Л = WТ = — W ).

В зависимости от значения амплитуды, частоты и длины волны колеблющиеся потоки среды можно классифицировать следующим образом (рис. 1).

По характеру изменения амплитуды и частоты колебаний во времени их можно разделить на две группы: стационарные и не­ стационарные колебания. К нестационарным колебаниям следует отнести колебания, амплитуда и частота которых изменяется во

9

Рис. 1. Структурная схема классификации колеблющихся потоков

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

При анализе теплового взаимодействия колеблющегося потока среды с поверхностью рассматриваемого тела следует выделить две области возможных частот колебаний — низкочастотные и высокочастотные. К низкочастотным колебаниям можно отнести такие колебания, при которых температура на поверхности тела изменяется во времени, а к высокочастотным — колебания, при которых температура поверхности тела практически не реагирует на колебания среды вследствие тепловой инерционности мате­ риала тела. Важной характеристикой колеблющихся потоков является длина волны Л или скорость ее распространения W.

В зависимости от соотношения между характерным размером тела L и длиной волны колебаний (например, колебаний скорости) колеблющиеся потоки можно подразделить на две группы:

длинноволновые колебания (характерный размер тела меньше чем половина длины волны колебаний L < Л/2);

коротковолновые колебания (характерный размер тела больше, чем половина длины волны колебаний L ^ L/2).

Последнюю группу колеблющихся потоков можно подразде­ лить на две подгруппы: резонансные и нерезонансные колебания.

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

10

Соседние файлы в папке книги