УДК 536.77
ТЕРМОГИДРОАКУСТИЧЕСКАЯ УСТОЙЧИВОСТЬ АВТОКОЛЕБАТЕЛЬНЫХ ПРОЦЕССОВ В АППАРАТАХ ПУЛЬСИРУЮЩЕГО ГОРЕНИЯ
В.И. Быченок1, Н.В. Мозговой2
Кафедра гидравлики и теплотехники, ТГТУ (1); кафедра теплотехники и промышленной теплоэнергетики, Воронежский государственный технический университет (2)
Представлена членом редколлегии профессором В.И. Коноваловым
Ключевые слова и фразы: аппараты пульсирующего горения; метод Ляпунова; метод Лагранжа; термогидроакустическая устойчивость.
Аннотация: Рассматривается проблема расчета термогидроакустической устойчивости автоколебательного возмущенного движения реагирующего газа в камерах сгорания теплоэнергетических установок, которое возникает в результате воздействия на рабочее тело акустических и тепловых возмущений. Для определения границ устойчивости анализируется возможность нахождения для таких систем функции Ляпунова и функции Лагранжа. Проводится сравнительный анализ применения их к процессам в энергетических установках. Показано, что в аппаратах пульсирующего горения с аэродинамическим клапаном для определения границ термогидроакустической устойчивости с успехом можно использовать метод Лагранжа. На основе этого метода получено уравнение для расчета резонансной частоты. Уравнение позволяет также решать задачу оптимизации геометрических размеров.
Традиционные требования, предъявляемые к теплоэнергетическим установкам, можно сконцентрировать в ключевых словах: экономичность, экология, простота и надежность. Всем этим требованиям в полной мере соответствуют теплоэнергетические установки, в которых реализуется механизм пульсирующего горения топлива [1]. Однако расчет их конструктивных размеров неразрывно связан с условиями устойчивости колебаний при горении. С физической точки зрения такой вид колебаний описывается системой нелинейных уравнений, а сами системы называются нелинейными. Дело в том, что «линейность» редко бывает свойством, присущим самой системе, в большинстве случаев линейность есть результат упрощения реальной системы. Класс нелинейных систем бесконечно шире и многообразнее, чем узкая область искусственно построенных линейных систем. Отметим некоторые существенные отличия в свойствах нелинейных и линейных систем. В линейных системах строго периодические колебания возможны только в форме вынужденных колебаний, возникающих от действия внешних периодических возмущающих сил. В то же время в нелинейных системах возможны периодические устойчивые свободные колебания даже при наличии сопротивления. При этом потери энергии автоматически компенсируются поступлениями ее из не колебательного источника. Дозировка поступления энергии во времени и по
величине регулируется самой колеблющейся системой. Такие системы называют автоколебательными. В нелинейных системах частота большей частью зависит от амплитуды колебаний, а в линейных системах такой зависимости нет. Поэтому простая линейная трактовка задачи о колебаниях не только не дает возможности раскрыть важные свойства колебательной системы, но и может существенно исказить реальную картину периодического процесса.
В то же время некоторые свойства нелинейных систем позволяют использовать для решения задачи ее устойчивости математический аппарат линейных систем. Так, например, для исследования малых колебаний упругих систем около положения устойчивого равновесия используют систему линейных уравнений. При этом пренебрегают членами, нелинейными относительно возмущений. Некоторые автоколебательные системы так же можно исследовать методами, характерными для линейных систем, если известна математическая зависимость механизма обратной связи. Так как в линейных системах собственная частота не зависит от начальных условий, и, в частности, от амплитуды, то ее изменение возможно только путем существенных изменений конструкции системы или перераспределением в ней масс. Это свойство позволяет использовать собственную частоту автоколебательной системы в качестве определяющего параметра, который можно задавать при проектировании аппарата пульсирующего горения на заданную тепловую мощность. Следовательно, трактовка периодических процессов в линейном приближении требует проведения подробного предварительного анализа механизма колебаний и введения обоснованных ограничений, чтобы математическое описание отражало действительную физическую картину поведения системы.
Так как горение является основным процессом всех энергетических установок, то исследование устойчивости этого процесса является основополагающей задачей. Решению этой задачи посвящено большое количество экспериментальных и теоретических работ. Наиболее подробный теоретический анализ проведен в работе К.И. Артамонова [2], где доказана принципиальная возможность линеаризация дифференциальных уравнений для возмущений. Доказательство базируется на том, что возмущения давления являются результатом возникновения акустических волн, а по амплитуде такие колебания можно считать малыми. В этом случае устойчивость описывается на основе фундаментальной теории Ляпунова для сплошной среды. Функция Ляпунова 2 в области устойчивого возмущенного движения обладает следующими свойствами
В работе [2] показано, что для исследования процесса горения на устойчивость в качестве функции Ляпунова можно использовать полную энергию потока реагирующего газа Е и воспользоваться нестационарной составляющей ее субстанциональной производной
где 8р, 8^, 8Q - возмущения давления, скорости и теплового потока; р - термический коэффициент расширения; Е - контрольная поверхность; V - исследуемый объем; р - плотность; Ср - изобарная теплоемкость.
Рассмотрим физический смысл величин, входящих в уравнение (2). Под энергией понимается выражение следующего вида
Z > 0, Z < 0 .
(1)
*=1
V
5м +S
( 2 ^ -w
dV,
(3)
т.е. рассматривается сумма внутренней и кинетическои энергии в возмущенном состоянии объема. Последний интеграл в (2) представляет собой критерий Релея, определяющий наиболее вероятные условия возбуждения акустических колебаний за счет периодического подвода тепла к объему. Наилучшие условия для возникновения колебаний будут при несовпадении фаз колебаний давления и скорости теплообмена. При наличии фазового сдвига ф между колебаниями, которые подчиняются гармоническому закону 5p = Ap cos ют и 5Q = Aq cos( ют -ф), интегрирование по периоду колебаний дает приращение акустической энергии равное
T г ApAQ 5*ак = -cosф j P------dV .
-cp
(4)
Очевидно, что усиление акустических колебаний возможно, если фазовый
сдвиг между 5р и 5Q удовлетворяет условию | ф | < п /2 . Последний интеграл в (2) можно записать в другом виде
jp— 5pdV = j—
-cp
PT d_(5s) d т
V p
5pdV.
(5)
Производная от вариации энтропии в (5) равна вариации от ее производной
ё(5^) 5(ds)
т.е. ----------------------------------------=-. В общем случае объем, в котором происходит горение,
ё т ё т
представляет собой открытую термодинамическую систему. Поэтому изменение энтропии в нем складывается из производства энтропии за счет внутренних процессов и потока энтропии через контрольную поверхность
ds d т
ds
( ds + 1 •
d т )e ^ d т/;
(6)
Единственный самопроизвольный процесс, который обычно учитывают при горении, - это образование продуктов сгорания вследствие химических реакций. В большинстве случаев этот процесс можно считать равновесным, поэтому производство энтропии (ё^ / ёт)е = 0. В то же время для учета реальных физических явлений объем, в котором происходит горение, надо считать открытой системой. Это означает, что изменение энтропии в ней вызывают колебания скорости подвода массы вещества. В связи с этим автор работы [2] отмечает, что колебательный подвод массы и колебательный подвод тепла эквивалентны с точки зрения возбуждения акустических колебаний, т.к.
( „ \
д_
5т
д_
5т
Q
pCpT
jpT -5QTSpdV = j—5
* Г\Г I J r
-CpT
pr,
cp
V p
ds d т
5pdV Ф 0. (7)
В зависимости от реальных факторов знак потока энтропии может быть < 0, > 0 или = 0. Поэтому априори знак интеграла (7) не определяется. Акустические волны, описываемые уравнением (7) часто называют «энтропийными волнами».
Физический смысл первого интеграла в правой части уравнения (2) подробно раскрыт в работе Б.В. Раушенбаха [3]. Поэтому здесь приведем только основные моменты.
Рассматривается канал постоянного сечения, в некотором объеме которого происходит горение. Зона горения ограничена двумя неподвижными сечениями: входное - 1 и выходное - 2. Пусть в одномерном потоке возникли возмущения, в результате которых давление и скорость получают отклонения от их стационарных значений. Величины отклонений определяются на основе гармонического закона 5р = |5р| 5ш(ют) и 5^ = |5^| $т(ют + ф), где ф - сдвиг фаз между возмущениями скорости и давления. Очевидно, что |5р| = Ар и |5^| = Лк являются амплитудами. Тогда текущие параметры возмущенного потока будут р = р0 + 5р и V = ^0 + 5^. Для анализа возмущенного движения произведениерм> разлагают в ряд
Интегрирование (8) за период колебания T = 2я/ю дает нам среднее значение
Поток энергии р0^0 не связан с колебаниями, т.к. соответствует стационарному течению (или фону). Учитывая гармонический закон для возмущений 5^ и 5р, выполним операцию интегрирования (9) за промежуток времени т = Т и получим
Выражение (10) Б.В. Раушенбах называет потоком акустической энергии через сечение и вводит обозначение
Отсюда следует первый вывод, что если в рассматриваемом сечении устанавливается узел давления (5р = 0) или узел скорости (5w = 0), то поток акустической энергии равен нулю. Второй важный вывод состоит в том, что для установившихся колебаний поток акустической энергии вдоль канала имеет постоянное значение A = const. Свойство A = const справедливо вне зоны горения, т.е. перед сечением 1 и после сечения 2 даже в том случае, когда на концах канала происходит рассеивание акустической энергии. Изменение А возможно внутри зоны горения вследствие протекающих там процессов. Подставляя (11) в первый интеграл в правой части уравнения (2), для одномерного потока получим
В зависимости от сдвига фаз ф величина А может быть как положительной, так и отрицательной. Если акустическая волна движется вправо вдоль оси х, то А >0, если влево, то А <0. Обозначим поток акустической энергии пересекающей
рм> = p0w0 + p0 5w + w0 5p + 5w5p.
(8)
произведения (pw)
0
(10)
(11)
0
2
(12)
сечение 1 А1, а сечение 2 - А2 . Тогда суммарное количество акустической энергии в зоне горения будет Ах = Л2 - А1. Введем правило знаков. Если Ах > 0, то область генерирует акустическую энергию, и ее поток направлен из области в окружающую среду, если Ах < 0 - область поглощает акустическую энергию и поток акустической энергии направлен внутрь области. Соответственно, если Ах = 0, то зона горения является нейтральной по отношению к акустической энергии. При наличии потока массы газа надо учитывать потери акустической энергии. Пусть эти потери сосредоточены на концах канала а и в. С учетом правила знаков потери на конце а обозначим -Ка, а на конце в - +Кв. Суммарные потери акустической энергии будут К = Кв - Яа. Тогда из энергетических соображений ясно, что при установившихся колебаниях Ах = К. Если Ах > К, то колебания в системе возбуждаются, а при Ах < К - колебания гасятся. Очевидно, что условие Ах = К соответствует границе устойчивости.
Здесь надо заметить, что в установках стационарного горения ставится задача, в результате решения которой определяются геометрические размеры зоны горения, обеспечивающие выполнение условия Ах < К. Для надежной работы камер пульсирующего горения необходимо условие Ах = К, а точнее Ах-К < |е|, где е - величина, заданная заранее. Именно из последнего условия формулируется задача устойчивости работы пульсирующего горения. При этом возникает задача вычисления величин Ка и Ке. Эту задачу можно решить, если известны конструктивные размеры трубы и методы расчета диссипативных эффектов на ее концах. Подставим (3), (7) и (12) в (2) и получим выражение, описывающее устойчивость нестационарного движения, вызванное акустическими возмущениями
д_
дт
5м + 5
( 2 ^ pw
dV =-JAdS + J-^5^—1 dV. (13)
1 vCp ^dT 1
Уравнение (13) представляет собой конкретный вид функции Ляпунова, которая может применяться при исследованиях термогидроакустической устойчивости процесса горения. Не трудно заметить, что при этом используется система отсчета Эйлера. Из анализа уравнения (13) можно сделать следующие выводы.
1. Для определения знака функции Ляпунова по уравнению (13) необходимо решать линеаризованную систему дифференциальных уравнений в частных производных.
2. Для учета диссипативных эффектов при анализе процесса горения на устойчивость необходимо знать граничные условия, которые определяются конкретной конструкцией.
3. Для определения «энтропийных» волн необходимо математическое описание кинетики химических реакций при горении и учитывать наличие потока эн-тропи, который так же зависит от конструктивных параметров зоны горения.
Перечисленные условия подчеркивают сложность решения задачи термогидроакустической устойчивости процесса горения. В самом общем виде уравнение (13) имеет большое значение при проведении теоретического анализа процесса горения на устойчивость. Для решения задач термогидроакустической устойчивости в конкретных теплоэнергетических установках предлагается использовать метод Лагранжа. Обоснование применения метода Лагранжа к исследованиям на устойчивость неравновесных консервативных потоков приведено в работе [4]. Практическое использование метода Лагранжа приводится в работе [5] при расчете регулярных аэротермоакустических автоколебаний в дозвуковых газотурбо-машинах.
Источник =(/ в о Регулирующее Емкостное Инерционное
энергии устройство звено (У, M1) звено (V2, M2)
(воздух + Gj. = const (клапан + (продукты (газ в резонанс-
горючее) зона горения) сгорания) ной трубе)
і Ь і ь
Sl S2 ,
поверхности разрыва
обратная связь
Рис. 1 Составные элементы автоколебательной системы АПГ
В настоящей работе предлагается использовать метод Лагранжа для исследования термогидроакустических колебаний в аппаратах пульсирующего горения (АПГ) с аэродинамическим клапаном, физическая модель рабочего процесса которого предложена в работе [1]. Согласно этой модели рабочее тело в АПГ представляет собой автоколебательную систему с обратной связью. Схема приведена на рис. 1
Введем основные допущения, соответствующие этой схеме.
1. Исследуется на устойчивость возмущенное движение потока идеального газа, образующего колебательную систему, состоящую из емкостного звена с объемом У в камере сгорания и инерционного звена с объемом У2 в резонансной трубе. Рассматриваемые объемы заполнены сплошной средой, которая имеет соответствующие массы М\ и М2.
2. Движение объемов установившееся и неразрывное в гидродинамическом понимании. Рассматривается одномерное движение сосредоточенных масс вдоль оси х со скоростями w-[ и ^2, соответственно. В акустических терминах такое движение называют фоном. Указанные на рис. 1 поверхности разрывов Е1 и Е2 имеют чисто математический смысл, который вводится для того, чтобы применить метод разрывных решений [6]. Под разрывом понимается мгновенное изменение состава и свойств потока в указанных сечениях проточного канала. Такой прием часто называют математической идеализацией.
3. В результате периодических возмущений в зоне горения в потоке газа в резонансной трубе возникает стоячая акустическая волна. При этом длина волны больше длины резонансной трубы, а точнее 1тр = 1/4 X. Следовательно, и длина емкостного звена намного меньше длины акустической волны.
4. Акустические возмущения в резонансной трубе накладываются на фон и вызывают смещение поверхности разрыва X со скоростью Б = ёх2 /й%, что приводит к деформации объема У, который заполнен реагирующими продуктами сгорания.
5. В результате деформации объема У в нем изменяется давление на величину 5р. Это, в свою очередь, приводит к изменению расхода воздуха Ова, поступающего в зону горения через аэродинамический клапан. Но изменение расхода воздуха приводит к отклонениям коэффициента избытка воздуха, участвующего в процессе горения на 5а. Следовательно, 5а(т) = Овд(х) / Ов°. В результате происходит изменение количества тепловой энергии, поступающей с потоком в объем У через неподвижную поверхность разрыва Таким способом осуществляется обратная связь, а система становится автоколебательной.
Энергетическое взаимодействие элементов автоколебательной системы можно показать на основе закона превращения энергии, записав его для возмущений в виде последовательности взаимодействия процессов по схеме на рис. 2.
3
1
гг
ди(т) + Д(PV)(т) = 80(т) - SL'(t)
4
Рис. 2 Схема обратной связи автоколебательной системы
В методе Лагранжа наблюдение ведется за фиксированной частицей или за сплошной средой в целом и прослеживается изменение ее параметров. При этом независимыми переменными, помимо времени, должны являться признаки, позволяющие отличать одну среду от другой. Это может быть масса, плотность, температура среды и т.д. Всегда можно найти такую систему отсчета, по отношению к которой пространство является однородным и изотропным, а время -однородным, т.е. д/ дт = ё / ёт . Такая система называется инерциальной и функция Лагранжа будет иметь следующий вид Ь = К(г, г)-П(г), где К -кинетическая и П - потенциальная энергии являются функциями обобщенных координат г и обобщенных скоростей г. При известной функции Лагранжа уравнения возмущенного движения в обобщенных координатах г , г принимают вид
Следовательно, возникает задача определения вида функции Лагранжа для конкретной автоколебательной системы.
В рассматриваемом случае кинетической энергией в возмущенном движении обладает масса газа в резонансной трубе (инерционное звено). Потенциальная энергия силового поля сосредоточена в продуктах сгорания (емкостное звено). Согласно принятой схеме взаимодействия на рис. 2, мы имеем два силовых поля: поле механической энергии и поле тепловой энергии. В связи с тем, что по определению тепловая энергия не связана с деформацией объема, то возникает вопрос определения действующих сил теплового поля. Из второго закона термодинамики известно, что параметром тепловой энергии является энтропия S, при этом существует связь между возрастанием энтропии и скоростью процесса. Если £ = £ (а), то с изменением а во времени связано изменение энтропии [7], тогда
Б =а(дБ/да)о . Следовательно, производная (дБ/да)0, которая вычисляется в установившемся течении (индекс «0»), рассматривается как «сила», вызывающая изменение а. Эта связь указывает на природу тепловых или «энтропийных» волн. При этом а = а (т). Для открытой системы такая трактовка будет означать, что поток энтропии вызывает изменение коэффициента избытка воздуха а, но справедливо и обратное утверждение: вариация 5а - вызывает изменение потока энтропии. Зона горения представляет собой открытую термодинамическую систему, поэтому необходимо доказательство, что действующие на нее силовые поля консервативны.
(14)
-> nk
Продукты
сгорания
Hnc
Рис. 3 Расчетная схема процессов в зоне горения
В пульсирующем режиме горения расход воздуха, поступающего в зону горения (область А на рис. 3) изменяется в соответствии с законом акустических возмущений, т.е. G/ = Ge(x). При этом, в связи с особенностью механизма течения потока через аэродинамический клапан [8], часть воздуха G'E постоянно участвует в поддержании процесса горения, а вторая его составляющая G"E = /(т) частично участвует в горении, а другая его часть выбрасывается из области А назад и создает реактивную силу. За счет этого обеспечивается установившееся и направленное в сторону +х движение потока в канале АПГ. Такой механизм позволяет считать, что открытая реагирующая зона горения является консервативной системой, т.к. в установившемся потоке при наличии в нем акустических возмущений за полный период а = const . В этом случае установившемуся потоку соответствует равновесное состояние продуктов сгорания. Для таких процессов в работе [9] получено уравнение, которое при равновесном протекании химических реакций имеет вид
dG = - SdT + Vdp + у p Td а .
(15)
Согласно Ландау [10] при переходе тела из одного состояния в другое, тело совершает максимальную работу Rmax, для обратного перехода требуется затрата внешним источником минимальной работы Rmm. Очевидно, что в потенциальном силовом поле эти работы совпадают, но имеют разные знаки. Для элементарной работы сил внешнего силового поля в работе [10] приводится выражение
dRmin = (T - Tj) dS - ( p - P0) dV . (16)
В равновесном состоянии (с индексом 0) T = T0 и р = р0 и работа d Rmin = 0.
Следовательно, в системе отсчета Лагранжа обобщенными координатами
являются S и V. Принимаем в уравнении (15) в состоянии равновесия открытой системы Т = const, и р = const. Тогда dG = уpjdа, а производная
(dG / dа)p t = уp t имеет смысл внутренней силы системы. Соответственно,
обобщенная сила внешнего поля, действующего на открытую систему, имеет противоположный знак или
| . (17)
p ,т
(уp,T )сист (уp,T )поля f
На основании (17) можно считать, что потенциальной энергией внешнего силового поля для открытой термодинамической системы в системе отсчета Лагранжа является функция состояния П = О = и + р¥- ТБ(а), а обобщенными координатами при равновесном протекании реакций будут 5 V и 5 Б(а).
Прежде чем написать функцию Лагранжа, с помощью которой будем анализировать устойчивость возмущенного движения, необходимо отметить следующее. Рассматриваемая схема автоколебательной системы на рис. 1 предполагает, что система с объемом V1 взаимодействует с системой объемом V2. В таком случае говорят, что система 2 движется в заданном внешнем поле, создаваемом системой 1. Для такого движения в работе [11] приводится функция Лагранжа в следующем виде
L = К2 (Г2 , Г2 ) - П (Г2 , Г (т))
(18)
или в наших координатах
L = К2( SVl) - П (5Vl, а (т)),
(l9)
а уравнение возмущенного движения, вызванного акустической волной, будет
d_
d т
дІ2 V д5 V 'l у
дІ2
д5V1
(i = 1)
(2О)
или
d_
d т
дК 2
дб V1
дП
+--------(i = 1).
д5V
(21)
Уравнение (21) показывает, что мы исследуем устойчивость возмущенного движения относительно возмущений одной независимой обобщенной координаты
5 VI. В этом случае для кинетической энергии объема У при / = £1 = £2 получаем
1 М2 ( d5V Y 1 M2 • 2
К =------2|------1 I =---------^( 5V1)2.
2 f2 I dт Jo 2 f
(22)
Потенциальная энергии в возмущенном состоянии для малых отклонений равна
П = По +l
і дП
Vд5V1
1 у0,5а
5V1 + f—1 5а +
V^aJo^V
д5V1
(5Vl )2 + 2
1 у
0,5 а
( д2 П ^ д5Vlд5a
5V[5a +
О
д5а2
(5а)2
0,5V
(2З)
По определению потенциального силового поля первые вариации от потенциальной энергии в равновесном состоянии (с индексом «0») дают нам обобщенные силы, которые равны нулю. Будем считать, что 5V-l является независимой обобщенной координатой, а коэффициент избытка воздуха есть функция от величины деформации объема, т.е. 5а =/(5V2). Функциональная зависимость а от 5V2 устанавливается на основе закона обратной связи (см. рис. 2). Используя принятое выражение для потенциальной энергии П = О, имеем
П - По = 1 5V1)2 - 2
2 I M,
Фі
да
RT +
^д2 HT (а) да2
0,5V
Л
(5а)2
д (RT )i
да
0,5V
5V,5a +
(24)
Продифференцируем (24) по 5V, и получим
д( RT V
^гт 2 2
дП Pi a
д (5Vi) M,
5V -
дрі
да
RT+
0,5V
да
0,5V
Sa+
д2 Hг (a)
да2
д(5а)
д (5Vi)
Sa = 0. (25)
Для двухкомпонентных топлив
( д2 Hт (a) ^ да2
= (xO )2 h GB0;
д(5а) 1 д(Од ) /кл д(р^)в „
а выражение ----------=------------=------------представляет собой обратную
Р д(5¥і) ов0 д(5¥і) ов0 д(5^) Р Р '
связь между возмущениями 5У1 и 5а. После подстановки (22) и (25) в (21) находим
M2
f
2 •• Pi al2
5V1 +±^- 5V1 -2 1 M, 1
Ч
Фі
да
RT+
0,5V
д( RT )i
да
0,5V
Sa+(xO )2 h^G0 ^ Sa = 0. (2б) д (5V2)
Приведем полученное выражение к классическому виду. Для этого поделим его на коэффициент, стоящий перед кинетической энергией и перенесем в правую часть члены, не содержащие обобщенной координаты 5 VI. В результате уравнение нестационарного движения в системе отсчета Лагранжа примет вид
5Vl +paJ- 5v = - f
У,^----------иу, —----------
P2V1l2 P2V2
дрі
да
RT +
5V
д( RT )i
да
5V
- (xO )2 KG | 5a . (27)
д (8V2)L
Для сокращения записей формул введем обозначения, которые обычно используются в теории колебаний
2 р1а1
Ш =
, h = р2у2 р2У2
Sp1
да
RT +
5У
д ( RT )1
да
ЪУ
- (х0 )2 hG ЦО)
д (ЪУ1)
Да
(28)
В (28) произведена замена 5а = /(т) = Аа sin (ют + ф), где Аа - амплитуда колебаний а. Правая часть (27) отражает вынужденные колебания, генератором которых является зона горения. Следовательно, выражение (27) представляет собой уравнение вынужденных колебаний системы с одной степенью свободы в системе отсчета Лагранжа. Подставим (28) в (27), тогда
5у +ю 5у = hsin(ют + ф) = hsinРт .
(29)
Решение уравнений такого вида известно из теории колебаний материальных точек и представляется в виде суммы общего решения однородного уравнения
5У1 +ю2 5 У = 0
(30)
и частного решения уравнения
5у +ю2 5у = к 8ш( ют + ф). (31)
Из решения однородного уравнения мы находим частоту собственных колебаний ю. Частота вынужденных колебаний р = (ю + ф) зависит только от внешних условий. Известно, что для надежной работы АПГ в пульсирующем режиме, необходимым условием является возникновение резонанса. Эти условия создаются за счет определенных геометрических размеров проточных каналов АПГ. В этом случае собственная частота свободных колебаний системы совпадает с частотой вынужденных колебаний ю = р. Уравнение (27) не учитывает наличие в периодических процессах диссипативных явлений (теплопотери, трение и др.). Поэтому из уравнения (27) невозможно определить сдвиг фаз ф между вынужденными и свободными колебаниями системы. Учет диссипативных эффектов в уравнении (27) в принципе возможен, но в данной работе не рассматривается. Проведенный авторами анализ показывает, что учет трения дает сдвиг фаз между возмущениями ф = п /2. Вычисление производных по а приведено в работе [12].
Выводы
Проведенный сравнительный анализ двух методов расчета термогидроакустической устойчивости позволяет сделать следующие выводы.
1. Метод Ляпунова носит более общий характер, т.к. система уравнений справедлива для решения задач термогидроакустической устойчивости в широком классе теплотехнических устройств, в том числе в трехмерной постановке.
2. Применение метода Ляпунова может вызвать затруднения при решении линеаризованной системы возмущенного движения в частных производных, осо-
бенно в трехмерной постановке. Так как система в частных производных интегрируется численными методами, то для достижения желаемой точности приходится увеличивать число точек в исследуемом диапазоне изменения параметров. Однако, как известно из публикаций, увеличение числа точек не всегда приводит к повышению точности.
3. В систему уравнений в методе Ляпунова геометрические размеры аппарата не входят, поэтому при интегрировании они должны задаваться в виде граничных условий на основе экспериментальны данных. Влияние «энтропийных» волн учитывают на основе эмпирических зависимостей для конкретного вида топлива и условий его сжигания, или вводят эмпирические коэффициенты в уравнение сохранения энергии в виде коэффициентов полноты сгорания.
4. Метод Лагранжа дает нам дифференциальные уравнения возмущенного движения в субстанциальных производных, для решения которых разработан соответствующий математический аппарат.
5. Возмущенное движение, вызванное термоакустическими волнами, в случае одной независимой обобщенной координаты описывается одним уравнением в полных производных. При этом даже одно уравнение включает в себя процессы, связанные с механизмом горения. Это позволяет одновременно учитывать влияние, как акустических возмущений, так и тепловых («энтропийные» волны).
6. В уравнение возмущенного движения, записанное в системе Лагранжа, входят геометрические размеры теплотехнического аппарата, что позволяет решать задачи оптимизации на этапе его проектирования.
7. Сложность применения метода Лагранжа к решению задач термогидроакустической устойчивости заключается в нахождении потенциальной энергии силового поля, обеспечивающего устойчивый режим колебаний.
Список литературы
1. Быченок В.И. Рабочий процесс аппаратов пульсирующего горения для сушки и термовлажностной обработки материалов / В.И. Быченок // Современные энергосберегающие тепловые технологии (сушка и термовлажностная обработка материалов): Тр. междунар. науч.-практ. конф. - Москва: МГАУ, 2002. - Т. 4. -С. 184 - 187.
2. Артамонов К.И. Термогидроакустическая устойчивость / К.И. Артамонов. - М.: Машиностроение, 1982. - 260 с.
3. Раушенбах Б.В. Вибрационное горение / Б.В. Раушенбах. - М.: Физмат-гиз, 1961. - 500 с.
4. Extended Conservation Laws from Hamilton's Principle Nonequilibrium Dissipative Fluids / Stanislaw Sieniutycz // Periodica Polytechnica. Ser. Phys. and Nucl. Sci. - 1994. - Vol. 2. No. 1-2. - Pp. 61 - 83.
5. Козырев В.Т. Регулярные автоколебания аэротермоакустической системы дозвуковых газотурбомашин (АКАТАС ГТМ) / В.Т. Козырев // Известия вузов. Машиностроение. - 1991. - № 4-6. - С. 62 - 70.
6. Седов Л.И. Механика сплошной среды / Л.И. Седов. - М.: Наука,1976. -Т. 1. - 535 с.
7. Беккер Р. Теория теплоты / Беккер Р. - М.: Энергия, 1974. - 504 с.
8. Быченок В.И. Расчет геометрических размеров аэродинамического клапана аппарата пульсирующего горения (АПГ) / В.И. Быченок, В. А. Русин // Труды ТГТУ. - Тамбов, 2004. - Вып.15. - С. 97 - 102.
9. Быченок В.И. Применение методов термодинамики для анализа динамических процессов в камерах пульсирующего горения / В.И. Быченок // Вестник ТГТУ. - 1998. - Т 4. № 4. - С. 495 - 503.
10. Ландау Л.Д. Теоретическая физика. Т.У. Статистическая физика: / Л.Д. Ландау, Е.М. Лифшиц. - М.: Наука, 1976. Ч. 1. - 584 с.
11. Ландау Л. Д. Теоретическая физика. Т. I. Механика / Л.Д. Ландау, Е.М. Лифшиц. - М.: Наука, 1973. - 208 с.
12. Быченок В.И. Частные производные термодинамических функций продуктов сгорания по химическому составу топлива и использование их для экстраполяции / В.И. Быченок, Р.А. Мухамедзянов // Труды КАИ. - Казань, 1973. -Вып. 153. - С. 54 - 59.
Thermal Hydroacoustic Resistance of Autooscillation Processes in Pulsating Combustion Apparatuses
V.I. Bychenock1, N.V. Mozgovoy2
Department «Hydraulics and Heat Engineering», TSTU (1);
Department «Heat Engineering and Industrial Heat Power Engineering, Vorinezh State Technical University (2)
Key words and phrases: thermal-hydroacoustic resistance; Lyapunov’s method; Lagrangian method; pulsating combustion apparatuses.
Abstract: The problem of calculating thermal-hydroacoustic resistance of au-tooscilating perturbed motion of reacting gas in combustion chambers of heat engineering plants is considered. This motion is the result of influence of acoustic and heat perturbation on working body. In order to determine the borders of resistance the possibility of finding Lyapunov’s function and Lagrangian function is analyzed. Comparative analysis of their application to the processes in power engineering plants is carried out. It is shown that Lagrangian method can be successfully used for determining borders of thermal-hydroacoustic resistance in pulsating combustion apparatuses with aerodynamic valve. The equation for calculating resonance frequency is obtained on the basis of this method. The equation enables to solve the problem of optimization of geometric sizes as well.
Thermohydroakustische Stabilitat der Autoschwankunsprozesse in den Apparaten des pulsierenden Brennens
Zusammenfassung: Es wird das Problem der Rechnung der thermohydroakus-tischen Stabilitat der autoschwingenden emporten Bewegung des reagierenden Gases in den Kammern der Verbrennung der warmeenergetischen Anlagen betrachtet. Diese Bewegung entsteht als Ergebnis der Einwirkung auf den Arbeitskorper der akustischen und thermischen Emporungen. Fur die Bestimmung der Stabilitatsgrenzen wird die Moglichkeit des Auffindens der Ljapunow-Funktion und der Lagrange-Funktion fur solche Systeme analysiert. Es wird die Vergleichanalyse ihrer Anwendung zu den Pro-
zessen in den energetischen Anlagen durchgefuhrt. Es ist aufgezeigt, daft die Lagrange-Methode in den Apparaten des pulsierenden Brennens mit dem aerodynamischen Ventil fur die Bestimmung der Grenzen der thermohydroakustischen Stabilitat mit Erfolg be-nutzt werden kann. Auf Grund dieser Methode ist die Gleichung fur die Rechnung der Resonanzfrequenz erhalten.
Stabilite thermohydroacoustique des processus d’auto-oscillation dans les appareils de la combustion pulsatoire
Resume: Est envisage le probleme du calcul de la stabilite
thermohydroacoustique du mouvement perturbe d’auto-oscillation du gaz reagissant dans les chambres de combustion des installations thermoenergetiques qui surgit au cours de l’action sur un corps des perturbations acoustiques et thermiques. Pour la definission des limites de la stabilite est analysee la possibilite de la recherche pour ces systemes des fonctions de Liapounov et de Lagrange. Est realisee l’analyse comparative de leur application pour les processus dans les installations energetiques. Est montre que dans les appareils de la combustion pulsatoire avec une soupape pour la definission des limites de la stabilite thermohydroacoustique on pourrait utiliser la methode de Lagrange. A la base de cette methode est regue une equation pour le calcul de la frequence de resonance. L’equation permet aussi de resoudre le probleme de l’optimisation des dimensions geometriques.