Анализ характера возмущенного процесса на основании уравнения (31.4) возможен двояким образом. Во-иер — вых, путем рассмотрения числешшх примеров и, во-вторых, путем теоретического исследования частных случаев, допускающих аналитическое рассмотрение. В настоящем параграфе будет использована первая из названных возможностей.
Найдем по аналогии с расчетами, приведенными в предыдущей главе, границу устойчивости при наличии потерь на концах трубы, т. о. будем искать такие значения 6Х и ЬЕ, при которых колебательная система нейтральна. Будем считать режим течения и положение плоскости подвода тепла заданными (т. е. заданными п, Мг, М2, li и /2), а в качестве переменного параметра выберем частоту колебаний оз. Так как по условию v=0 (рассматривается граница устойчивости), то тем самым заданы (для каждого ©) все входящие в уравнение (31.4) величины, кроме коэффициентов Сх, Съ С з и С4, зависящих (через аи, ап, я21, а23) от искомых ЬХ и ЬЕ.
Сравнив формулы (17.1) и (31.1), легко видеть, что прп ЬЕ=ухр0 и bX=y. zv0 (где уj и у* — комплексные коэффициенты, которые математически выражают существование некоторого механизма обратпой связи):
Рассмотрим для определенности случаи возбуждения колебаний при ЬЕ = 0. При сделанном предположении //i=0 и уравнение (31.4) становится линейным относительно уг соотношением, из которого последняя величина легко определяется.
Результаты такого рода расчета приводятся ниже. Расчет был проведен для двух случаев. В первом анализировалось влияние частоты колебаний па положение границы устойчивости, во втором — влияние изменения
Относительной длины трубы ~ ^точнее, обратной величи — L
Ны d ) Па ТОї же фактор. В обоих случаях расчет велся для
ЛД=0,15; п= 2,5; т^ 1; ^-^=0,5.
Ь первой серии расчетов, имевшей целью проиллюстрировать влияние частоты колебаний на положение границы устойчивости, величина ± была принята равной 0,2, т. е. рассматривалась сравнительно короткая труба, длина которой равнялась 5 .калибрам. Результаты расчета прп — ведены на рис. 61. Здесь по оси абсцисс отложена вещественная, а но оси ординат — мнимая часть у2. Комплексный коэффициент у2 дает возможность судить об амплитудно-фазовом соотношении между ЬХ и v0. В предыдущей
Рін Gi. Іїлііннітс частоты колебаний со на ноложепио границы устойчивости.
Главе п ранее было получено чрезвычайно простое условно самовозбуждения колебапин ігріт отсутствии потерь:
Угол между фазами vQ и б Л" должен быть менее Поскольку при bX~y2v„ фазы 6.Y u v0 совпадают при положительных (вещественных) значениях у2, границей устойчивости в рассмотренном ранее примере является мнимая ось, а областью неустойчивости — вся правая полуплоскость на рис. 61. Как видно из расположения кривых у2 на рис. 61, прп учете потерь область неустойчивости уменьшается — границы устойчивости смещаются вправо от мнимой оси и смещаются тем сильнее, чем больше частота колебаний (на диаграмме «неустойчивая» сторона границы устойчивости заштрихована). При этом условие возбуждения усложняется — теперь существенна пе только фаза ЬХ, но и амплитуда.
Наименьшая относительная амплитуда величины ЬХ (гс условно принято за единицу), прп которой возможно возбуждение, получается тогда, когда фазы ЬХ п и() совпадают (у2 — действительная величина). На рис. 61 приведены лишь две ветви уг (соответствующие первой н второй гармопнкам). Как видно пз диаграммы, увеличение частоты (переход от первой ко второй гармонике) нлечет за собою увеличение абсолютных величин т. е. увеличение относительных амплитуд ЬХ, потребных для возбуждения системы.
Этот результат легко поддается объяснению. Оцепим влияние частоты колебаний на потери акустической эпергии в окружающем трубу пространстве. С этой целью воспользуемся выражением (30.7), которое позволяет вычислить поток акустической эпергии, движущейся через концевое сечение трубы. При этом учтем формулу (30.0), позволяющую выразить | £ |2 через вещественную и мнимую части нмпеданца z. Тогда
Полученная формула указывает, в частности, что знак А (направление потока акустической энергии) свя — зап со знаком коэффициента активного сопротивления г. Эта формула может быть упрощепа следующим образом. Заметим, что выражение, стоящее в прямых скобках, близко к единице, так как для рассматриваемых задач величины 7-их малы по сравнению с единицей. Это видно, например, нз приближенной формулы для z (30.8) и*особо подчеркивалось в конце § 30.
Таким образом, для достаточно «длинных» труб и для низких гармоник
(32.2)
Приближенное равенство (32.2) показывает, что поток акустической энергии через концевое сечение Трубы пропорционален квадрату частоты. Следовательно, умеиь — шеиие вероятности возбуждения по мере перехода к высшим гармоникам вполне естественно и ход кривых на рис. 61 надо считать типичным. Более высокие частоты дали бы границы устойчивости, еще более сдвипутые вправо.
Иа практике величина г/2 пе может увеличиваться беспредельпо. Реально существующие механизмы обратпой связи обусловливают существование верхнего предела этой величины. Если, в частности, верхний предел у. г больше значений, требуемых для возбуждения первой гармоники, ио меньше значений, требуемых для возбуждения второй (точка А па рис. 61), то возбуждение колебаний основного тона возможно, а всех более высоких гармоник — невозможно. Рассмотренный пример хорошо объясняет то обстоятельство, что в опытах всегда наблюдается возбуждение колебании сравнительно низ — коп частоты.
Оценим порядок у2. Из экспериментов известно, что при развитом процессе вибрационного горения амплитуда колебаний давления бр имеет порядок перепада давления между холодным н горячим газом, т. е. порядок поднора, вызванного тепловым сопротивлением. Этот перепад может быть выражен следующей формулой:
Таким образом,
Переходя к безразмерным переменным и приняв, что в первом приближении —:= — = /?,-, получим:
При развитых колебаниях амплитуда колебаний скорости имеет порядок скорости установившегося течения холодного газа и, следовательно,
Порядок величины ЬХ^р.,—рj можно считать совпадающим с порядком величии рх її р2. Тогда
І г’і I
Более вероятна меньшая величина у. г, поскольку разность давлений слева и справа от зоны горения обычно составляет не более х12—V3 полной амплитуды колебаний б р.
В рассмотренном здесь численном примере (гс3 = 6,25; Л/j =0,15) порядок абсолютпой величины у2 будет, таким образом, но более |y2| = 0,25. Сравнивая это значение с абсциссами точек пересечения кривых у2 с вещественной осыо (рис. 61), нетрудно прийти к выводу, что в данном случае возможно возбуждение лишь первой гармоники системы и то при особеппо благоприятных обстоятельствах. Это впдпо из того, что окружность у21=0,25, захватывает лшнь весьма незначительную часть области неустойчивости (прп выбранных на рпс. 61 масштабах эта окружность перешла в эллипс). Такой результат связан с тем, что принятое в расчете отпошение -^ = 0,2 слишком велико. По сути, здесь рассматривается слишком «короткая» труба (£.•= 5 с?), для которой относительная величина потерь будет большой.
Это впдно, например, из формулы (32.2), по которой потери растут не только пропорционально квадрату частоты, но и пропорциональпо квадрату отноніепня — j — . Следовательно, по мере уменьшения ~ количество пзлучасмой
В окружающее пространство акустической энергии должно падать п колебательная система стаповится все более и более склонной к возбуждению. Этим обстоятельством объясняется тот известный экспериментаторам факт, что внбрацпонпое горение можно наблюдать лшнь в достаточно «длинных» трубах, т. е. в трубах с малым значением параметра.
Чтобы придать этим качественным соображениям наглядность, па рнс. 62 приведены границы устойчивости для вторых гармоник тон жо колебательной системы, что и рассматривавшаяся выше, по для трех значений параметра — j-: 0,2; 0,00 и 0. На приведенной диаграмме хоро-
‘I
Шо видно, что г уменьшением отношении — треоуемая
Устойчивости,
Для возбуждения системы величина | //„ ] уменьшается d
И стремится к нулю вместе С — J — .
Экспериментальным подтверждением иолученпых здесь теоретических выводов может служить наблюдение, сделанное Киркби и Уилером1). Указанные авторы пишут, что при диаметре трубы 100 мм возбуждение вибрационного горения в пей становилось возможпым лишь при общей длине трубы, превышающей 1400 мм.
Если вернуться к численному примеру и продолжать считать верхпей границей | у2 | значения 0,25, то из рис. 62
!) К і г k b у, Wheeler, Jonrn. of Ни — Llheni. Soc., 1931, стр. 847.
Сдодуот, что при Л с/ ^-^- = 0,06^ вторая гармоника М. ожрт возбудиться, в то время как для -]г—0,2 это невозможно.
Приведенный выше пример демпфирующего влияния увеличения частоты колебаний позволяет сделать наглядными те рассуждения о возможпостн препебреженпя начальными условиями, которые приводились в главе II.
Пусть колебательная система имеет те же параметры, что были использованы в численном примере § 23. Единственное отличие будет заключаться в том, что потерями на излучение акустической энергии во внешшою среду пренебрегать не будем. На рпс. 33 приведены величины о) п v, полученные для разных гармоник без учета потерь па излучение. Если учесть эти потери, то v п со для первой гармопикп вряд лп изменятся сколько-нибудь заметным образом. Вторая и третья гармоники и без учета потерь на излучение задемнфировапы достаточно сильно (v < 0), учет потерь на пзлучепне лишь усилит это демпфирование. Что касается четвертой гармоники, то без учета потерь на нзлученио она, как и первая гармоника, почти нейтральна. Одпако учет потерь резко изменит это положение, так как частота этой гармоники в 4,5 раза превосходит частоту основпого тона, а потери на излучение возрастут, следовательно, в 20 раз. Это соображение позволяет считать все гармоники выше третьей задемтг — фировапными зпачителыю сильнее первой. Будем поэтому в последующих расчетах учитывать лишь три первые гармоники.
Пусть характер процесса в зоне горепия, на основании которого были построены кривые рис. 33, слабо изменится. Пусть указанное слабое изменение произошло так, что первая гармоппка стала нейтральной (v — 0). Что касается второй и третьей гармоник, то сохраним соответствующие им значения v, хотя учет потерь скорее всего увеличит декременты затухания.
Пред ни 11 і, і,чі у,. .1То в начальный момепт т=0 в трубе происходили с.-иысное колебательное движение воздушных масс, составленное нз стоячих волн первых трех гармоник, причЄ5г амплитуды второй н третьей гармопик вдвое превышали амплитуды стоячей волны осповпого тона, а фазы всех трех гармоник совпадали. Тогда эпюра мгновенного значения возмущения р имела бы вид, изображенный в верхней части рпс. 03. Вид этой же эпюры в последующие моменты времени т15 т2, т3 дан на том
6’3. Постепенное исчезновение влігшшя начальных условий.
Же рисунке ниже. Нетрудно видеть, что уже через время tad2 сохраняются лишь колебания основного тона. Это время соответствует приблизительно 6—7 колебаниям основпого то на. Если учесть, что продольные акустические колебания, с которыми приходится пметь дело па практике, находятся в пределах 20—1000 герц, то очевидно, что весь эффект, связанный с начальными условиями, исчезнет за доли секунды-
Приведенная здесь грубая оценка имела целью дать наглядное представление о быстроте, с которой сглаживается влияние начальных условий в рассматриваемых системах.