CC BY
170
УДК 629.7:532.59:539.3
Вестник СибГАУ Т. 16, № 2. С. 386-394
ЗАТУХАНИЕ ВЫНУЖДЕННЫХ КОЛЕБАНИЙ И ГИДРОУДАРНЫХ ВОЛН В ТУРБУЛЕНТНОМ ПОТОКЕ ПРОТЯЖЁННЫХ ТРАКТОВ СИСТЕМ ПИТАНИЯ ЖИДКОСТНЫХ РАКЕТНЫХ ДВИГАТЕЛЕЙ
С. А. Барабанов 1, В. И. Бирюков 2*
1 ООО «Энергоавтоматика» Российская Федерация, 127206, г. Москва, Чуксин тупик, 9 2 Московский авиационный институт (национальный исследовательский университет) Российская Федерация,125993, г. Москва, ГСП-3, Волоколамское шоссе, 4
*E-mail: [email protected]
Волновые процессы в трубопроводных системах являются причиной неравномерности расходов подаваемых через них жидкостей. При открытии или закрытии регулирующей арматуры системы питания жидкостных ракетных двигателей возникают одиночные ударные волны, которые являются возмущениями давления в широком спектре частот и могут приводить к возбуждению неустойчивости рабочего процесса в газогенераторах или камерах сгорания. Целью работы является экспериментальное исследование процессов затухания вынужденных колебаний и одиночных гидроударных волн в турбулентных потоках слабосжимаемой капельной жидкости в гладких гидравлических магистралях круглого сечения. Поскольку процесс затухания ударных волн заметен только в протяженных трактах с отношением длины к диаметру трубопровода более нескольких тысяч, то методология проведения данной работы была построена на возбуждении гидроударных волн в магистральной транспортной системе, осуществляющей перекачку углеводородных горючих и регистрацию давления датчиками, расставленными равномерно вдоль магистрали. При исследованиях создавались возмущения с амплитудой более чем в 2 раза большей уровня постоянного гидродинамического шума. Одиночные гидроударные волны возбуждались путем проведения отбора жидкости. При этом в ответ на открытие крана (начало отбора) генерировалась волна разрежения, а на закрытие - волна сжатия. Проведённый спектральный анализ сигнала показал, что основная часть спектра гидродинамических колебаний находится в частотном диапазоне до 10 Гц. Более высокие частоты, генерируемые работой насосных агрегатов, быстро затухают. В результате экспериментов была получена аппроксимационная зависимость декремента затухания одиночных волн сжатия и разрежения в виде функции от чисел Рейнольдса. При переходных режимах течения жидкости от ламинарного к турбулентному динамические изменения происходят внутри ламинарного слоя и показатель затухания ударной волны стремится к постоянной величине. На основании сделанных предположений экспериментальные точки аппроксимированы экспоненциальной кривой вида y = ^0 + A [1 -exp(-Bx)] .
По результатам проведенной обработки были получены значения для коэффициентов и выполнено сравнение расчетных данных с результатами экспериментов. Показано удовлетворительное согласие результатов испытаний на модельной установке и на протяженных магистралях. Полученная экспериментальная зависимость, описывающая затухание гидроударных волн в турбулентных потоках, используется при проведении инженерных расчетов систем питания жидкостных ракетных двигателей на кафедре ракетных двигателей в Московском авиационном институте.
Ключевые слова: гидроударная волна, спектр частот, декремент затухания, резонанс, турбулентность, число Рейнольдса.
Vestnik SibGAU Vol. 16, No. 2, P. 386-394
ATTENUATION OF THE COMPELLED FLUCTUATIONS AND HYDROSHOCK WAVES IN THE TURBULENT STREAM OF EXTENDED PATHS OF POWER SUPPLY SYSTEMS
OF LIQUID ROCKET ENGINES
S. A. Barabanov1, V. I. Biryukov2*
1 JSC "Power automatics" 9, Chuksin deadlock, Moscow, 127206, Russian Federation 2Moscow Aviation Institute (National Research University) 4, Volokolamskoe shosse, GSP-3, Moscow, 125993, Russian Federation *E-mail: [email protected]
Wave processes in pipeline systems are at the bottom of non-uniformity of the expenses submitted liquids through them. At opening or closing of regulating armature of a power supply system of liquid rocket engines there are single shock waves which are pressure indignations in a wide spectrum offrequencies and can lead to excitation of instability of working process in gas generators or combustion chambers. The work purpose is the experimental research ofproc-esses of attenuation of the compelled fluctuations and single hydroshock waves in turbulent streams poorly compressed a drop liquid in smooth hydraulic highways of round section. As process of attenuation of shock waves is swept up only in extended paths with the relation of length to diameter of the pipeline more than several thousand the methodology of carrying out of the given work has been constructed on excitation of hydroshock waves in the main transport system which is carrying out swapping hydrocarbonic combustible and registration of pressure by gauges, placed in regular intervals along a highway. At researches indignations with amplitude more than in 2 times of level of constant hydrody-namic noise were created. Single hydroshock waves were raised by carrying out the selection of liquid. Thus in reply to crane opening (the selection beginning) the wave pressure reduction, and on closing - a compression wave was generated. The carried out spectral analysis of a signal has shown that the basic part of a spectrum of hydrodynamic fluctuations is in a frequency range to 10 Hz. Higher frequencies generated by work ofpump units, quickly fade. As a result of experiments approximation dependence of decrement of attenuation of single waves of compression and pressure reduction in the form offunction from Reynolds's numbers has been received. At transitive modes of a current of a liquid from laminar to turbulent dynamic changes occur in a laminar layer and the indicator of attenuation of a shock wave aspires to a constant. On the basis of the made assumptions experimental points are approximated exponential by a kind curve y = + A [i -exp (-Bx)] . By results of the spent processing values for factors have been received and
comparison of the settlement data with results of experiments is executed. The satisfactory consent of results of tests to modeling installation and extended highways is shown. The received experimental dependence describing attenuation of hydroshock waves in turbulent streams, is used at carrying out of engineering calculations of power supply systems of liquid rocket engines on chair "Rocket engines" at the Moscow aviation institute.
Keywords: a hydroshock wave, a spectrum of frequencies, decrement of attenuation, a resonance, turbulence, Reynolds's number.
Введение. Разработанные в 1980-х гг. маршевые жидкостные ракетные двигатели (ЖРД) РД 170 и РД 171 для ракетоносителей (РН) «Энергия» и «Зенит» на компонентах керосин и жидкий кислород с высоким уровнем тяги (806,4 тс тяги в пустоте) характеризуются значительными расходом компонентов топлива (2393 кг/с) и высоким давлением в камерах сгорания (24,5 МПа). Неизбежное возбуждение ударных волн при открытии или закрытии отсечных клапанов во время проведения испытаний ЖРД сопряжено с повреждением измерительной и управляющей арматуры протяженных магистралей заправки. Одним из основных средств борьбы с гидроударными волнами является установка в магистралях демпферов и других устройств, снижающих энергию ударных возмущений [1-7]. Однако всегда остается открытым вопрос влияния волновых возмущений в магистралях на режим горения в камере ЖРД и, как следствие, возможное возбуждение неустойчивости горения [5; 8].
Из экспериментов на магистралях стендовой системы и теоретически определено, что в случае отсутствия потока ламинарная модель движения жидкости вполне удовлетворительно описывает акустические характеристики тракта в широком диапазоне частот вплоть до пятого резонанса продольных колебаний [1; 2; 9; 10]. Для трактов относительно малой длины (Ь/О < 200-300) эффектом затухания волновых возмущений при расчётах процессов, протекающих в них, можно пренебречь. Данных об экспериментальных исследованиях распространения волновых возмущений в протяженных трактах (где отношение диаметра к длине 1: 1000 и более) в известных публикациях нет. В то же время влияние единичных гидро-
ударных волн и волновых процессов на характеристики турбулентных течений в протяженных трактах являются важным аспектом при построении физико-математических моделей систем питания мощных жидкостных ракетных двигателей [1; 2; 6; 7; 11], а также при моделировании течений и режимов работы магистральных трубопроводов и других задач машиностроения.
Исследование затухания колебаний гидроакустической волны в протяжённых круглых трактах осложнено трудоёмкостью создания экспериментальной установки необходимого масштаба. Именно по этой причине отсутствуют экспериментальные работы, посвященные исследованию процессов затухания волновых возмущений при их распространении в трубах с большим соотношением длины к диаметру. Можно разработать физико-математическую модель распространения волн в турбулентном жидкостном потоке, основываясь на известных классических работах Н. Е. Жуковского [12], Г. Шлихтинга по теории пограничного слоя [13], И. А. Чарного по турбулентным течениям в трубах [14]. Данные для коэффициента трения при ламинарном и турбулентном режимах течения жидкости в гладких каналах можно получить из работ Блазиуса, Халанда, Свами [15] и др. Однако достоверность результатов численных расчетов может быть подтверждена только эмпирическим путем. Учитывая практическую направленность работы, целью настоящего исследования является создание инженерной методики расчета затухания ударных возмущений и вынужденных колебаний в жидкостных турбулентных потоках.
Квазистационарное решение. Определение параметров гидроударной волны и показателя её затухания
дал Н. Е. Жуковский в своей работе «О гидравлическом ударе в водопроводных трубах» [12]. Аналитическое решение для коэффициента затухания волны в круглом тракте, полученное из предположения квазистационарности течения, было опубликовано в работе И. А. Чарного в 1975 г. Он, решая задачу распространения волны, использовал метод уравнений характеристик [14]. Представив уравнения для стационарного течения и возмущенного потока, он получил соотношение, связывающее изменение скорости u с пройденной волной расстояния х:
d 2 5u
4aD
-i
du
XquQ
, -X|-Xu|
(1)
где а - скорость звука в жидкости; Б - внутренний диаметр круглого тракта; X
ния стационарного и возмущенного течения; И1 скорость распространения волны возмущения в сечении х = 0; и - скорость распространения волны возмущения при х > 0; и0 - стационарная скорость течения до внесения возмущения в сечении х.
Коэффициент трения X появляется при переходе от двумерного течения к одномерному, когда уже нельзя связать напряжения вязкого трения с изменением профиля скорости. В этом случае гидравлические потери давления из-за трения связывают с изменением скорости эмпирическими зависимостями типа уравнения Блазиуса. В то же время известно, что влияние сил трения жидкости о стенки гидравлического тракта зависит от частоты [12]. Используя безразмерные вариации скорости и давления в виде отношения возмущений скорости и давления 5и' и др' к их средним
й ди' д др'
стационарным значениям ди = —, др =-, а также
Рср
безразмерных времени t —
ср t
(l / a)
_ х
х = 1, уравнения движения жидкости в малых отклонениях можно представить в виде
pua Ô5u Ô5p Ap „ ---— +——+——5u — 0;
Рср Ôt ÔX Рср
Ô5p + pua Ô5u о Ôt Рср Ôx
(2)
5p — 5p ов'ш ; 5u — 5ù0emt,
(3)
где др0 и дй0 - амплитуды вынужденных колебаний давления и скорости, и получено решение в виде
+ 1 ю2 -|5м — 0.
(4)
0, X - коэффициенты трем -
dx2 ^ pua j
Член в скобках формулы (4) является волновым
,2-2 — Ap числом k — ю -ira——. При относительно малом pua
гидравлическом сопротивлении и достаточно боль-
„ _ Ap _2 шои частоте выполняется условие ю-^ ю , тогда
pua
при решении уравнения (4) вторым слагаемым в волновом числе k можно пренебречь. В полученном решении возмущения давления и скорости при их распространении вдоль тракта зависят от координаты и времени
5p — PAp к_ - a exp i ( t - kx )- B exp i (юT + kx )]. (5)
Уравнение (5) описывает распространение волн постояннои формы от источника в разные стороны трубы. Первое слагаемое в квадратных скобках описывает распространение возмущения в прямом направлении, а второе - в обратном. Постоянные A и B зависят от исходного сигнала возмущения. На практике мы чаще встречаемся с эффектом затухания возмущения, т. е. когда влиянием сил вязкого трения пренебрегать нельзя. В этом случае обе составляющие волнового числа являются соизмеримыми величинами, а решение уравнения (4) приобретает несколько другую форму, в котороИ появляется дополнительный множитель exp [-Im(k) x ]. Для гармонических колебаний, распространяющихся в прямом направлении, возмущения давления записываются следующим образом:
5p — A exp i (ю t - kx )exp [im (k )x ].
(6)
и координаты
Мнимая часть комплексного волнового числа является коэффициентом затухания 1т (к) = "Л. Значение коэффициента А получим из граничных условий при х = 0, так как в этом случае возмущения равны колебаниям, возбуждаемым источником: А = др0.
В экспериментах показатель (декремент) затухания определяем как величину, обратную расстоянию, на котором амплитуда волны уменьшается в е раз:
1 5p0 1 5ÙQ q — — ln—- — — ln—-.
x 5p x 5u
(7)
где Ар - гидравлические потери давления из-за вязкого трения в тракте.
В работе [2] рассмотрен пример распространения возмущений вдоль магистрали. Решая данное уравнение для гармонических колебаний частоты ю (безраз-
„ — ю1
мерный параметр частоты ю = —):
а
Для гармонических колебаний в работе [2] аналитически было получено значение декремента затухания
q — ю
/ \2
1 1 + 4 ^-Ap 1 -1
2 \ ^pua j
(8)
Задачу определения коэффициента затухания для нестационарных течений в гидравлической магистрали можно решить аналитически лишь при специальных допущениях [14]. В общем случае необходимо выполнить точную математическую модель в частных производных и решить полученную систему уравнений численными методами.
В данной работе мы экспериментально исследовали затухание возмущений в узком диапазоне частот от 1 до 10 Гц и определяли зависимость коэффициента затухания от параметров течения в длинной магистрали.
Эксперимент. Волновые колебания в трубе - это распространяющиеся возмущения давления и скорости, которые связаны между собой известной формулой Н. Е. Жуковского. Поскольку измерение скорости потока технологически сложнее измерения давления, то в экспериментах регистрировались стационарные значения давления и пульсации давления равномерно расположенными вдоль трубы датчиками. При длинах трактов в несколько десятков километров, применительно к задачам транспортировки углеводородных топлив по магистральным трубопроводам, точки измерения контролируемых параметров устанавливались на расстоянии 10-20 км в зависимости от длины исследуемого трубопровода. Аналоговый сигнал датчика давления подавался на вход контроллера сбора, обработки информации и преобразовывался АЦП контроллера в цифровой сигнал. Оцифрованный сигнал фиксировался в оперативной памяти устройства, упаковывался и передавался на сервер базы данных, где были собраны в единую базу данных сигналы от всех точек измерений.
При анализе процессов распространения одиночной гидроударной волны в турбулентном потоке тракта мы столкнулись с задачей выделения полезного сигнала (волнового фронта) на фоне постоянного гидродинамического шума. Сам по себе гидродинамический шум представляет собой колебания различных частот, создаваемых работой насосных агрегатов и прочего технологического оборудования. Обрабатывая сигнал, в котором содержится только гидроди-
намический шум, можно наблюдать за процессами затухания волновых возмущений в широком спектре частот. В экспериментах инициировалась одиночная ударная волна с амплитудой более чем в 2 раза большей уровня постоянного гидродинамического шума. Гидроударные волны создавались путем отбора жидкости. При этом в ответ на открытие крана (начало отбора) генерировалась волна разрежения, а на закрытие - волна сжатия (рис. 1).
Эксперименты проводились на магистральных трубопроводах круглого сечения, осуществляющих постоянную перекачку углеводородных топлив. Физические характеристики перекачиваемой жидкости и параметры течения по каждому объекту представлены в таблице.
Проведённый спектральный анализ сигнала, измеренный преобразователем избыточного давления ТЖИУ 406, установленным на расстоянии 4 -104 Бу ,
показал, что основная часть спектра гидродинамических колебаний находится в частотном диапазоне до 10 Гц (рис. 2). Высокие частоты, генерируемые работой насосных агрегатов, успевают затухать на расстоянии 4 • 104 .
Измеряя амплитуду единичных гидроударных волн, а также и общий уровень гидродинамического шума, состоящего из колебаний частотой до 10 Гц, были получены идентичные картины затухания. Для каждого экспериментального образца протяженного тракта проводилась серия экспериментов по определению снижения пика давления волновых возмущений. Для этого строились кривые зависимости величины амплитуды гидроударной волны 5р от расстояния (рис. 3).
Рис. 1. Измеренные сигналы преобразователей давления, приведённые к отклонениям вблизи 0 вычитанием постоянной составляющей
Характеристики исследуемых объектов
№ п/п Длина, км Диаметр, м Плотность, кг/м3 Эффективная вязкость, сСт Скорость течения, м/с Скорость звука, км/с Число Рейнольдса (Яе),105
1 290 1 872,3 21 0,442 1,040 0,21
2.1 23 0,8 864,8 13,4 0,665 1,118 0,39
2.2 1,744 1,103 1,04
3 170 0,5 848 1 1,381 1,161 7,21
4 40 0,5 839,7 2 2,235 1,123 5,47
5 227 0,5 844 6 0,849 1,150 0,85
6 259 1 879 18 1,75 1,060 0,95
Рис. 2. Типичный спектр гидродинамического шума в тракте, записанный преобразователем
Рис. 3. Процесс затухания одиночной гидроударной волны
Следуя формуле (7), для каждого эксперимента набор точек был аппроксимирован при помощи стандартных функций математического программного обеспечения функцией вида /(х) = Аехр(Ьх) . Значения коэффициента Ь являются показателем (декрементом) затухания гидроударной волны. Очевидно, что затухание волновых процессов в турбулентном потоке происходит в основном благодаря силам вязкого трения. Отношение сил трения к силам инерции потока характеризуется числом Рейнольдса. В этой связи представляется целесообразным определять величину декремента затухания одиночной волны в зависимости от числа Рейнольдса:
ц = / (Re ) = f
Ч V
(9)
где Dy = const - внутренний диаметр исследуемого тракта; u0 - скорость течения принималась постоянной и равной среднему значению скорости; v - вязкость перекачиваемого продукта.
Вычислялась эффективная вязкость по формуле Блазиуса, исходя из уравнения потерь давления на трение на исследуемом участке на стационарном режиме, как
Ар = 0,3164 Re
l
Pu0
D,, 2
(10)
На основе обработанных данных всех экспериментов был построен график зависимости коэффициента затухания п от безразмерного параметра Яе. Результаты каждой серии экспериментов были обобщены и получено одно среднее значение коэффициента затухания п для каждой серии. Полученные экспериментальные значения коэффициента затухания приведены на рис. 4.
При малых числах Рейнольдса Яе << 105 течение в пограничном слое становится ламинарным [3] и функция затухания стремится к постоянной величине п0, которая может быть вычислена при помощи формулы Стокса для закона трения. Постоянная величина коэффициента волнового затухания в ламинарном потоке равна
Ло
16v
(11)
При распространении возмущений в потоках при больших числах Рейнольдса динамические изменения происходят внутри ламинарного подслоя. Дальнейшее увеличение числа Рейнольдса не приводит к уменьшению ламинарного подслоя и, следовательно, не даёт изменений в коэффициенте затухания, и в результате он остаётся постоянным и стремится к постоянной величине. Таким образом, искомая функция имеет асимптоту п = а при больших числах Рейнольдса. На основании сделанных предположений экспериментальные точки были аппроксимированы кривой вида
у = Л0 +А [1 -ехр (-Бх)].
По результатам проведенной обработки были получены следующие значения для коэффициентов:
А = 2 -105; Б = 10-5.
Запишем полученную зависимость в следующем виде:
- / _ V
1 - ехр
Л
1
Л0 + — 0 р
Re
Re
кр у
(12)
где — = 47,32 ± 0,32 = 47,3 км - расстояние, на котором сигнал гидроударной волны затухает в е раз в турбулентном пограничном слое, а коэффициент в зависит от толщины ламинарного подслоя, который в свою очередь является функцией от шероховатости стенки магистрали; Яекр = (0,897 ± 0,022) ■ 105 ~ 105 -
критическое значение числа Рейнольдса, при котором толщина ламинарного подслоя вблизи поверхности трубы становится менее 1 % от диаметра тракта, при этом дальнейший рост потерь из-за вязкого трения постепенно уменьшается, а коэффициент затухания
1
выходит на постоянную величину —.
Рис. 4. Зависимость коэффициента затухания от безразмерного числа Рейнольдса
Предложенная эмпирическая зависимость позволяет проводить простую оценку затухания колебательных процессов и одиночных гидроударных волн, распространяющихся вдоль протяженного тракта круглого сечения.
Проведённые на следующем шаге исследования сравнения результатов полученной зависимости с теоретическими расчетами распространения возмущений, полученных вычислительными методами, показали, что наибольшая сходимость эксперимента и теории наблюдается в области малых значений числа Рейнольдса. В области больших значений Re результаты теоретических расчетов расходятся с экспериментом: затухание за счет вязкого трения в эксперименте больше рассчитанных значений.
Колебательный процесс состоит из «разгонных» волн сжатия и «тормозящих» волн разрежения. Проведенные серии экспериментов содержали исследования распространения и затуханиях волн обоих типов. Результаты экспериментов (рис. 5) подтвердили теоретические предположения, полученные У. Р. Лийвом в 1983 г. [15]: экспериментальные значения коэффициента затухания волны разрежения значительно меньше коэффициента затухания волны сжатия.
Например, для трубопровода № 3, результаты которого приведены на рис. 5, эта разница составила 0,009 1/км, т. е. затухание волн при торможении потока почти в 2 раза больше чем при разгоне.
На рис. 6, а приведено изображение распределения скорости потока по сечению магистрали в случае отсутствия волны давления. При распространении волны сжатия, т. е. при разгоне потока, профиль течения на фронте волны вытягивается (рис. 6, б).
В этом случае при u1 > u0 > 0 либо при u1 < u0 < 0 проинтегрированное уравнение (1) в предположении равенства значения Я = Я0 = const:
л2 =
Я0 u0 , 11 —■—'- +—ln
2aD„ x
1 +
4P0
2pcu0
Я0 Ы x
2aD„
(13)
При торможении потока (рис. 6, в) при щ > и0 > 0, а также и при щ < и0 < 0 , решение того же уравнения принимает вид
т Л2
Я0 u0 , 11 -1-L + — ln
2aD„ x
1 -
ДР0
2pcu0
Я0 |u01
2aD^,
1 - e у
(14)
Из сравнения полученных формул видно, что л22 <Л2, т. е. при одинаковых по модулю величинах скорости потока коэффициент затухания при разгоне потока будет больше, чем при торможении. При разгоне потока градиент скорости вблизи стенок трубопровода увеличивается, следовательно, потери на вязкое трение увеличиваются (рис. 6, б). При торможении потока профиль скорости становится таким, как изображенный на рис. 6, в, градиент скорости уменьшается, а с ним и потери на трение также уменьшаются.
Заключение. В результате проведенных экспериментальных исследований по затуханию ударных гидравлических волн на протяженных магистралях получены эмпирические соотношения показателя затухания в виде функциональной зависимости от критерия Рейнольдса.
Определено, что гармонические составляющие широкополостного спектра импульсного возмущения в турбулентном жидкостном потоке с частотой более 10 Гц при распространении в магистрали быстро затухают и несущественны. Увеличение показателя затухания низкочастотных составляющих спектра возмущения существенным образом определяется вязкостью и гидравлическими потерями и обратно пропорционально связано с волновыми свойствами среды и диаметром трубопровода. Найденные соотношения, описывающие затухание гидроударных волн в турбулентных потоках, используются при проведении инженерных расчетов систем питания жидкостных ракетных двигателей на кафедре ракетных двигателей в Московском авиационном институте.
Рис. 5. Различие коэффициентов затухания гидроударной волны сжатия и разрежения
Рис. 6. Изменение профиля скорости на фронте волны
Полученные зависимости будут полезны для оценок воздействия гидроударных волн на агрегаты, арматуру и трубопроводы при проведении приемосдаточных испытаний (КТИ) и других видов испытаний жидкостных ракетных двигателей на предприятиях, занимающихся их разработкой и изготовлением.
Библиографические ссылки
1. Гликман Б. Ф. Автоматическое регулирование жидкостных ракетных двигателей. М. : Машиностроение, 1974. 296 с.
2. Гликман Б. Ф. Нестационарные течения в пневмогидравлических цепях. М. : Машиностроение, 1979. 253 с.
3. Ламб Г. Гидродинамика. М. : Гостехиздат, 1947. 929 с.
4. Пухов В. А., Чучеров А. И. Стендовые огневые испытания ЖРД. М. : Машиностроение, 1971. 543 с.
5. Неустойчивость горения в ЖРД / под ред. Д. Т. Харрье, Ф. Г. Рирдона. М. : Мир, 1975. 870 с.
6. Махин В. А., Присняков В. Ф., Белик Н. П. Динамика жидкостных ракетных двигателей. М. : Машиностроение, 1974. 833 с.
7. Мошкин Е. К. Нестационарные режимы работы ЖРД. М. : Машиностроение, 1970. 336 с.
8. Бирюков В . И., Козлов А. А. Внутрикамерная низкочастотная неустойчивость рабочего процесса в дросселируемых ЖРД. М. : МАИ-ПРИНТ, 2007. 384 с.
9. Бирюков В. И., Мосолов С. В. Акустика газовых трактов жидкостных ракетных двигателей. М. : Изд-во МАИ, 2013. 164 с.
10. Экспериментальное исследование акустических характеристик столба жидкости в непроточных трубах / Е. Д. Барбашов [и др.] // Акустический журнал. 1996. Т. 42, № 4. С. 478-488.
11. Беляев Е. Н., Черваков В. В. Математическое моделирование ЖРД. М. : МАИ_ПРИНТ, 2009. 280 с.
12. Жуковский Н. Е. О гидравлическом ударе в водопроводных трубах. М. ; Л. : Гостехиздат, 1949. 103 с.
13. Шлихтинг Г. Теория пограничного слоя. М. : Наука, 1974. 713 с.
14. Чарный И. А. Неустановившееся движение реальной жидкости в трубах. 2-е изд. перераб. и доп. М. : Недра, 1975. С. 186-191.
15. Лийв У. Р. Теоретические и экспериментальные основы расчета напорного ускоренного движения жидкости в цилиндрических трубах: автореф. дис. ... д-ра техн. наук : 05.14.09. Л., 1983. 48 с.
References
1. Glikman B. F. Avtomaticheskoe regulirovanie zhidkostnykh raketnykh dvigateley [Automatic control of liquid rocket engines]. Moscow, Mashinostroenie Publ., 1974, 296 p.
2. Glikman B. F. Nestatsionarnye techeniya v pnevmogidravlicheskikh tsepyakh [Non-stationary currents in pneumatichydraulic chains]. Moscow, Mashinostroenie Publ., 1979, 253 p.
3. Lamb G. Gidrodinamika [Hydrodynamics]. Moscow, Gostekhizdat Publ., 1947, 929 p.
4. Pukhov V. A., Chucherov A. I. Stendovye ognevye ispytaniya ZhRD [Bench fire tests LPR]. Moscow, Mashinostroenie Publ., 1971, 543 p.
5. Neustoychivost' goreniya v ZhRD [Liquid Propel-land Rocket Combustion Instability]. Editor D. T. Xarrje, ass. edit. F. H. Reardon]. Moscow, Mir Publ., 1975, 870 p.
6. Makhin V. A., Prisnyakov V. F., Belik N. P. Di-namika zhidkostnykh raketnykh dvigateley [Dynamics of liquid rocket engines]. Moscow, Mashinostroenie Publ., 1974, 833 p.
7. Moshkin E. K. Nestatsionarnye rezhimy raboty ZhRD [Non-stationary operating modes LPR]. Moscow, Mashinostroenie Publ., 1970, 336 p.
BecmHUK Cu6rAY. TOM 16, № 2
8. Biryukov V. I., Kozlov A. A. Vnutrikamernaya nizkochastotnaya neustoychivost' rabochego protsessa v drosseliruemykh ZhRD [Intrachamber low-frequency instability of working process in throttled LPR], Moscow, MAI-PRINT Publ., 2007, 384 p.
9. Biryukov V. I., Mosolov S. V. Akustika gazovykh traktov zhidkostnykh raketnykh dvigateley [Acoustics of gas paths of liquid rocket engines], Moscow, Izd-vo MAI Publ., 2013, 164 p.
10. Barbashov E. D., Glikman B. F., Kazakov A. A., Morozov S. A. [Experimental research of acoustic characteristics of a column of a liquid in landlocked pipes], Akusticheskiy zhurnal, 1996, Vol. 42, No.4, P. 478-488 (In Russ.).
11. Belyaev E. N., Chervakov V. V. Mate-maticheskoe modelirovanie ZhRD [Mathematical modeling LPR], Moscow, MAI_PRINT Publ., 2009, 280 p.
12. Zhukovskiy N. E. O gidravlicheskom udare v vodoprovodnykh trubakh [About hydraulic blow in water pipes]. M.-L., Gostekhizdat Publ., 1949, 103 p.
13. Shlikhting G. Teoriya pogranichnogo sloya [The interface theory]. Moscow, Nauka Publ., 1974, 713 p.
14. Charnyy I. A. Neustanovivsheesya dvizhenie real'noy zhidkosti v trubakh [Unsteady movement of a real liquid in pipes]. Moscow, Nedra Publ., 1975, P. 186-191.
15. Liyv U. R. Teoreticheskie i eksperimental'nye osnovy rascheta napornogo uskorennogo dvizheniya zhidkosti v tsilindricheskikh trubakh. Avtoref. dis. d-ra tekhn. nauk. [Theoretical and experimental bases of calculation of the pressure head accelerated movement of a liquid in cylindrical pipes. Dr. techn. sci. diss. author's abstr.], Leningrad, 1983, 48 p. (In Russ.).
© Eapa6aHOB C. A., EHPMKOB B. H., 2015
