
книги / Тепловые и гидродинамические процессы в колеблющихся потоках
..pdfБ. М. ГАЛИЦЕЙСКИЙ.
Ю.А. РЫЖОВ,
Е.В. ЯКУШ
ТЕПЛОВЫЕ И ГИДРОДИНА МИЧЕСКИЕ ПРОЦЕССЫ В КОЛЕБЛЮ
ЩИХСЯ
ПОТОКАХ
Под редакцией профессора В. К . КОШ КИНА
М о с к в а
«М А Ш И Н О С ТРО ЕН И Е»
УДК 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