автореферат диссертации по авиационной и ракетно-космической технике, 05.07.05, диссертация на тему:Экспериментально-теоретическое исследование и разработка электрофизического метода диагностики ракетных двигателей

УДК 629.7.036.54

На правах рукописи

Рудинский Александр Викторович

Экспериментально-теоретическое исследование и разработка электрофизического метода диагностики ракетных двигателей

Специальность 05.07.05 - Тепловые, электроракетные двигатели и энергоустановки летательных аппаратов

Автореферат диссертации на соискание ученой степени кандидата технических наук

2 9 Ш 2015

Москва-2015

005570949

005570949

Работа выполнена в федеральном государственном бюджетном образовательном учреждении высшего профессионального образования «Московский государственный технический университет имени Н. Э. Баумана».

Научный руководитель: Д.т.н., профессор Ягодников Д. А.

Официальные оппоненты: д.т.н. Семенов В. И.;

(НПО «Энергомаш им. академика В. П. Глушко)

д.т.н. Головин Ю. М.

(ГНЦ ФГУП «Исследовательский центр им. М. В. Келдыша»)

Ведущая организация: ОАО Ракетно-космическая корпорация

«Энергия» им. С. П. Королёва (РКК «Энергия»)

Защита состоится «14» октября 2015 г. в 14 часов 00 минут на заседании диссертационного совета Д 212.141.08 при Московском государственном техническом университете имени Н. Э. Баумана по адресу: 105005, г. Москва, Лефортовская наб., д. 1, факультет «Энергомашиностроение».

С диссертацией можно ознакомиться в библиотеке и на сайте www.bmstu.ru Московского государственного технического университета имени Н. Э. Баумана.

Ваш отзыв на автореферат в двух экземплярах, заверенный печатью организации, просим направлять по адресу: 105005, г. Москва, 2-я Бауманская ул., д. 5, МГТУим. Н. Э. Баумана, ученому секретарю диссертационного Совета Д 212.141.08

Автореферат разослан «_»_2015 г.

Ученый секретарь диссертационного Sfít л/

совета д.т.н., доцент ¿¡М/'У/А^' В. В. Перевезенцев

ОБЩАЯ ХАКАКТЕРИСТИКА РАБОТЫ

Актуальность темы. В настоящее время средства аварийной защиты (САЗ) жидкостной ракетной двигательной установки (ЖРДУ) для любого летательного аппарата (ЛА) являются её неотъемлемой частью конструкции и обеспечивают её функционирование с заданной вероятностью безотказной работы 0,9995...0,9998.

В связи с вышесказанным появляется необходимость в создании таких активных САЗ, которые бы позволили осуществлять диагностику рабочего процесса в ЖРДУ в реальном времени с целью эффективного и своевременного обнаружения развития аномальных явлений. Существующие системы

диагностики, например, в ЖРДУ основываются на получении информации о работоспособности узлов двигателя, в частности турбонасосного агрегата, регистрируя частоту вращения и перемещения ротора, температуру генераторного газа и др. Описанные системы косвенно отражают рабочие процессы в двигателе и не могут дать первичную информацию о последних с необходимым быстродействием. Высокие значения температуры газового потока (3100.. .3900 К) в камерах сгорания ЖРД и сложность конструкции также накладывают определенные трудности при реализации стандартных методов и систем диагностики рабочих процессов, в частности зондовых.

Среди современных методов диагностики внутрикамерных физико-химических процессов, протекающих в газовых трактах энергосиловых установок, представляются наиболее перспективными бесконтактные методы, основанные, например, на регистрации электрофизических характеристик продуктов сгорания (ПС) топлив, которые содержат положительно и отрицательно заряженные ионы, электронный газ. Часть электрического заряда выносится из камеры двигателя основным потоком ПС. Реализуемое при этом течение можно назвать электрогазодинамическим, т.к. характеризуется малым значением плотности электрического тока (10 мА) и большим электрическим потенциалом (-100 кВ).

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

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

Цель диссертационной работы заключается в экспериментально-теоретическом обосновании использования электрофизических характеристик ПС ЖРД на углеводородном топливе для разработки бесконтактных эффективных средств диагностики и аварийной защиты.

Предмет исследований. В диссертации экспериментально и теоретически исследовалось влияние изменения режимных параметров модельной камеры ЖРД на электрофизические характеристики ПС углеводородного топлива.

Объект исследований. В качестве объектов исследований в диссертации рассматриваются камеры модельного и натурных ЖРД на углеводородном топливе и первичные измерительные преобразователи - датчики магнитного и электрического поля.

Научная новизна работы.

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

2. Получены расчетные зависимости от основных режимных параметров работы модельного ЖРД и камеры маршевого ЖРД РД-171М электрических потенциалов элементов конструкции сопла двигателя и амплитуды напряженности магнитного поля в диапазоне изменения коэффициента избытка окислителя а,,*= 0,5...0,85 для углеводородного топлива.

3. Определена близкая к линейной зависимость напряженности магнитного поля продуктов сгорания углеводородного топлива в области критического сечения от давления в камере сгорания, которая может быть рекомендована в качестве задающего параметра системы контроля и диагностики ЖРД.

4. С помощью частотного анализа колебаний напряженности магнитного и электрического полей впервые получены спектры сигналов и их частотно-амплитудные характеристики в зависимости от времени испытания и текущего давления в камере сгорания.

5. Разработана математическая модель и методика расчетного анализа параметров быстродействия системы аварийной защиты ЖРД на основе контроля характеристик магнитного поля продуктов сгорания при скачкообразном изменении давления в КС и площади критического сечения сопла.

Практическая значимость результатов заключается в следующем.

1. В разработке математической модели, описывающей генерацию и протекание электромагнитных процессов в ПС ЖРД при горении топлив с учетом геометрических особенностей проточного тракта камеры ЖРД и состава топлива, а также режимных параметров двигателя.

2. В выявлении зависимостей электрофизических характеристик ПС от режимных параметров ЖРД для назначения диагностических признаков протекания рабочих процессов.

3. Разработана структурная схема схемы системы аварийной защиты ЖРД на основе регистрации амплитуды напряженности магнитного поля, обладающей большим быстродействием (-0,02 с) по сравнению с традиционными средствами (-0,04...0,05 с).

Основные положения диссертации, выносимые на защиту.

1. Осесимметричная математическая модель, алгоритм и результаты численного анализа электрофизических характеристик модельного и маршевого ЖРД на углеводородном топливе (УВТ).

2. Методика огневых испытаний модельного ЖРД с регистрацией электрофизических характеристик ПС углеводородного топлива.

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

4. Результаты анализа возможности использования САЗ на основе регистрации амплитуды напряженности магнитного поля ПС. Достоверность полученных результатов подтверждаются.

1. Использованием фундаментальных уравнений газодинамики, магнитогидродинамики, теории горения.

2. Удовлетворительным согласованием результатов численных расчетов и проведенных в работе экспериментов, а также их удовлетворительным согласованием с экспериментальными данными других авторов.

3. Применением аттестованных средств измерения и регистрации, обеспечивающих необходимую точность для теплотехнических экспериментов. Апробация работы проведена на Всероссийской научно-технической конференции «Студенческая весна 2010», ФГБОУ ВПО МГТУ им. Н. Э. Баумана, (Москва. 2010); на Всероссийской научно-технической конференции «Ракетно-космические двигательные установки», ФГБОУ ВПО МГТУ им. Н. Э. Баумана, (Москва. 2010); на «XXXVI Академических Королёвских научных чтениях по космонавтике», ФГБОУ ВПО МГТУ им. Н. Э. Баумана, (Москва. 2012); на Всероссийской научно-технической конференции «Ракетно-космические двигательные установки», ФГБОУ ВПО МГТУ им. Н. Э. Баумана, (Москва. 2013); на XX научно-технической конференции молодых ученых и специалистов РКК «Энергия», (Москва. 2014)

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, выводов и списка литературы; содержит 142 страницы, 75 рисунков и 19 таблиц. Список литературных источников состоит из 83 позиции.

Основное содержание работы

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

Первая глава посвящена обзору экспериментально-теоретических исследований в области электрофизических методов диагностики рабочих процессов в ракетных и реактивных двигателях. Сформулированы задачи исследования. Приводится обзор математических моделей, используемых авторами для описания генерации заряженных частиц (ионов и электронов) и процессов разделения зарядов в пламени горелок, в камерах сгорания модельных и натурных ЖРД и ВРД. Проведен анализ выбранных граничных условий для решения задач распределения электрического тока, коэффициента электропроводности ПС и потенциала элементов конструкций камеры модельных и натурных ЖРД и ВРД. Приводятся выводы о необходимости учета взаимного движения ионов и электронов на участках камеры с наибольшей концентрацией заряженных частиц (камера сгорания) с целью повышения точности в определении характеристик магнитных полей. Выполнен анализ влияния коэффициента избытка окислителя на потенциал элементов конструкции как в энергетических силовых установках, так и в модельных ЖРД, что позволило выявить диагностические признаки состояния ЖРД. Проведен сравнительный анализ данных различных авторов по экспериментальному определению тока, суммарного электрического заряда в потоке продуктов сгорания и потенциала конструкции ракетных двигателей на различных типах топлив. Количественные расхождения в результатах, приводимых в обзоре, свидетельствуют о различных методах градуировки измерительных преобразователей (электростатических зондов) и условиях проведения экспериментов. Однако характер зависимостей электрофизических свойств от режимных параметров всех исследуемых двигательных установок идентичен. В заключении главы сформулированы задачи экспериментально-теоретического исследования.

Во второй главе описывается осесимметричная математическая модель электрофизических процессов в камере ЖРД. Предлагаемая модель включает в себя электрогазодинамический расчет истечения ПС углеводородного топлива (УВТ) из сверхзвукового конического сопла ЖРД и состоит из уравнений газовой динамики, в которые включены уравнения электродинамики и уравнения сохранения для электрически заряженных компонент (ионов и электронов). Применительно к осесимметричной геометрии камеры ЖРД выбрана двумерная модель течения. Ввиду малого отношения электрогазодинамической силы к инерционным силам отдельно исследовалась общая газодинамическая система уравнений. Затем с помощью найденных распределений газодинамических параметров определялись электрические токи, поля и концентрации электрически заряженных компонент газовой смеси на основе уравнений, описывающих только электрические явления (магнитогидродинамики, электростатики, уравнений Максвелла и т.д.).

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

ионизированной реактивной струи возникновение собственного электромагнитного поля протекает по следующим основным стадиям.

1. Генерация заряженных частиц (ионов и электронов) за счет химических реакций во фронте горения в сечении КС, где уровни концентраций электронов обеспечивают электрическую проводимость ПС а.

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

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

4. Создаваемое за счет такого движения заряженных частиц электрическое

—►

поле напряженностью Е и потенциалом фпс порождает магнитное поле, векторы напряженности Н которого расположены в радиальной плоскости и представляют собой окружности, соосные потоку ПС.

Последовательность подготовки исходных данных для решения газодинамической задачи истечения ионизированных ПС и электрофизической задачи, описывающей собственное электромагнитное поле ПС, схематично представлена на Рис. 1.

Рис. 1.

Схема расчета характеристик собственного электромагнитного поля ПС УВТ

Моделирование осесимметричного течения газа в сопле осуществлялось в программе «Ansys Fluent» с помощью системы уравнений, состоящей из уравнения: неразрывности (1), сохранения количества движения в проекциях на координатные оси г (2) и х (3) в форме Навье-Стокса, энергии (4) и уравнение состояния (5). Диссипация энергии и касательные напряжения определялись в

рамках стандартной двухпараметрической (к-е)-модели турбулентности уравнений (6) и (7).

из

d(pwx) | d(pwr) | pwr дх дт г

= О

i э_

г дх

. 3wx ах

+

г дх

(rpwxwx) + ¿I (rpwrwx) = - g + [гц (2

(rpwxwr) + £(rpwrwr) = -2 +1¿[r,(5 + Э] + [r,(2^) -

2 /3wx 3wr ууЛ1 _ w, 2 ti ídwx dwr wr\

3 Ux + дт + r JJ r^ 3 r \ öx 5r r /

ö(pwxH0) , a(pwrH0) i pwrH0 _ a2 /к, \ a2 /kt \ , iaH0r

~lr + Г 0xJ är2^ P = РКцТ

a(Pkwx) a r/ ut^ski awr

3x

r dr

: + тх

ax'

(O (2)

(3)

(4)

(5)

(6)

awr E2

Г~зГ~С2Ер1?

(7)

дх dr LV сткУ a(p£wx) _ a r/ jv\£e| ax 3rLV ajdrl где p - плотность, wr, wx - проекции скорости на оси г и х, р - давление, Т -температура, р. - вязкость, Но - энтальпия, kt - турбулентная теплопроводность, Ср — теплоемкость, R^ - газовая постоянная, к - турбулентная кинетическая энергия, е - скорость диссипации турбулентной энергии, - турбулентная вязкость, оЕ, Ok — числа Прандтля для кие соответственно, CiE, СгЕ - константы модели турбулентности, тхг - тензор напряжений.

Собственное электромагнитное поле ПС в модуле «Ansys Fluent MHD» описывалось уравнениями Максвела методом электрических потенциалов

ЗЕГ ЗЕХ „ (8)

дх Зг дЕг Ег

+-- + — :

дх Эг г

ЗЕ

4тг

-q.

(9)

заряд в объеме ПС;

здесь Е - напряженность электрического поля; ц во - диэлектрическая проницаемость вакуума.

Для модельной камеры ЖРД (Рис. 2, а) была построена конечно-элементная (КЭ) модель в виде структурированной сетки (Рис. 2, б).

Рис. 2.

Проточной часть (а) и осесимметричная КЭ модель камеры ЖРД (б) Моделирование рабочего тела проводилсь путем расчета термодинамических свойств ПС этилового спирта (СгНзСЩж)) с объемной концентрацией 75% и газооразного кислорода (СЬ) при а = 0,58.. .2. В результате термодинамического расчета горения указанной топливной пары в программе «Terra» были получены распределения теплоемкостей Ср в зависимости от

температуры ПС Т при исходном значении а, которые аппроксимировались полиномиальными зависимостями:

Cpi(T) = 1,939 10-7 Т3-7,539-10^Т2+0,98Т+1774, при а=0,58;

СР2(Т) = -1,875-10"I2-Ts+3,226-10"8T4-2,17-10j,-T3-0,713-T2-1143T+71650, при а=0,85;

Срз(Т) = 2,304• 10"7- Т3-9,377• 10"4 ■ Т2+1,515 ■ Т+б 12,807, при а=0,85.

Для исследования полученного магнитогазодинамического течения ионизированных ПС УВТ проводилось совместное решение систем уравнений (1) ~ (7) и (8) - (9)- При этом граничные условия дополнялись значениями электрической проводимости рабочего тела (ПС УВТ). В математической модели электрическая проводимость ПС УВТ определялась выражением

а = kCTqeNem, (10)

где - поправочный (кинетический) коэффициент; qe, - элементарный заряд; Ne — объемная концентрация электронов; це - подвижность электрона.

Подвижность электронов находилась из зависимости

Це me(ve>'

здесь ше - масса электрона; (ve) - средняя частота электронных столкновений, которая определялась выражением

(ve) = 3,64 • 10""6 • Ne(ln(A))1,s ,

где 1пЛ = In |l,24 • 107 j - кулоновский логарифм.

Плотность электрического тока уноса описывалась уравнением

¡у = ЧеВД, (11)

где Nj, Wj - концентрация и скорость положительно заряженных ионов.

Ток уноса, обусловленный объемным зарядом, выносимым из внутрикамерного пространства, определялся выражением

Iy = /S qe-NH30+-WH30+ds, (12)

где s — элемент площади слоя, занятый зарядом.

При методе, представленном на схеме (Рис. 1) концентрации заряженных частиц в ПС УВТ (газообразный кислород и 75% водный раствор этилового спирта), определялись с помощью термодинамических расчетов по программе «Terra» в зависимости от коэффициента избытка окислителя а.

По результатам термодинамических расчетов в программе «Terra» для указанной топливной пары носителями положительных зарядов, преобладающими в ПС, являются ионы НзО+. Носителями отрицательного электрического заряда принимались свободные электроны е".

Граничными условиями для решения системы уравнений (8) - (9) являлись потенциалы электростатического поля в КС и на срезе сопла.

Значение электрического потенциала во входном сечении камеры сгорания определялось выражением

\Edr=- J ^ ' }-RDdrMNr-NZ\RL

здесь NeK , Nh3o+ - концентрации электронов и ионов в КС, rd = - радиус Дебая.

Значение электрического потенциала в сечении среза сопла задавалось, исходя из условия выноса положительного заряда истекающей струей ПС из внутрикамерного пространства:

ьо

где , — концентрация ионов и электронов на срезе сопла, Та - статическая

температура на срезе сопла.

Расчет напряженности собственного магнитного поля, создаваемого положительными ионами НзО+, движущимися со скорость основного потока ПС, в месте установки датчика МП определелялся из уравнения

„ _ <^„30+ _ деМНзо+Ук„г (13)

М — -2~ — —7Г1-2 "Н30+1

здесь — объемный суммарный положительный заряд, Мн30+ - скорость ионов (равная скорости газа); элементарный заряд; Гдмп — расстояние от датчика МП до линии, вдоль которой движутся заряженные частицы (ионы); Ук — объем КС.

Описание граничных условий (ГУ) и области их приложения применительно к камере модельного ЖРД, приведено в Таблице 1.

Таблица 1.

Граничные условия магнитогазодинамической математической модели

Тип ГУ Описание ГУ

Inlet (вход в КС) Расход топлива тт = 0,5 кг/с; температура Тк = 3240 К (а = 0,85); давление рк = 3 МПа; параметры турбулентности; электрический потенциал <ркс1 = 0,268 В

Axis (Ось камеры) Градиент скорости, температуры и давления по нормали к границе в*, 0 а*. _ 0 аг _ 0 аР 0 дг ' дх ' дг ' дг

Wall (Стенка камеры) Адиабатная стенка, тепловой поток равен 0 (г.у. второго рода). Условие изолированной стенки, электрический ток равен 0 (¡у=0 А/м2)

Outlet (Вых. сечение сопла) Значения на границе находятся путем экстраполяции из центра ячейки: \уг = 0; ра = 0,099 МПа; Та = 293 К. Электрический потенциал при условии выноса положительного заряда (ра = 0,01 В

Тестирование математической модели электрофизических процессов в камере сгорания модельного ЖРД позволило в развитие метода использовать её для моделирования и исследования электрофизических особенностей рабочего процесса и продуктов сгорания кислород-керосинового топлива маршевого ЖРД РД-171М конструкции НПО «Энергомаш» имени академика В.П. Глушко. 8

Распределения электрофизических характеристик на оси камеры модельного и натурного ЖРД представлены на Рис. 3.

, А/м2

, А/м2

% В ,444

0,8 0,4

0 0,02 0,0^ ст-103, Ом'м-'

(¡1=0,851 а=0,5!

ПГ7

Рис. 3.

Распределение электрофизических характеристик ПС по длине камеры модельного ЖРД при а=0,85 (а) и РД-171М при а=0,78 (б): 1 - плотность электрического тока уноса ^ 2 - потенциал ф; 3 - электропроводность а По результатам полученных распределений (Рис. 3) и оценочных расчетов уровень тока уноса на срезе сопла, определявшийся по выражению (12), для модельного ЖРД составил I/ =117 мА при а - 0,85. Для натурного ЖРД I/ 4170 мА при а = 0,78. Плотность тока уноса, найденная из зависимости (12), для модельного ЖРД в критическом сечении (хкр=48,5 мм) составила ]уч> = 0,26 А/м2, в сечении среза (ха=98 мм) = 0,042 А/м2 при а = 0,85. Для натурного ЖРД в критическом сечении (х,ф=537 мм) Зу45 = 2 А/м2, в сечении среза (ха=2241 мм) = 6,188 10"3 А/м2 при а = 0,78. Электропроводность ПС в камере сгорания модельного ЖРД а = 0,675- Ю-3 Ом-'м"1 при а = 0,85 и а = 0,986-10"3 Ом-'м"1 при а = 0,78 для натурного ЖРД. Значение потенциала ПС в камере сгорания модельного ЖРД фксо = 268 мВ (проверить) при а = 0,85 и ф ксо = 444 мВ при а = 0,78 для натурного ЖРД. Уровни напряженностей собственного магнитного поля ПС в месте установки датчика магнитного поля (на внешней стенке в области критического сечения сопла), найденные из выражения (13), для модельного ЖРД составили Н1 =0,818 А/м при сп = 0,58; Нг = 8,679 А/м при а2 = 0,85 и Нз = 0,019 А/м при а = 2, что соответствует температурам потока в КС Тк = 2803; 3240; 2445 К.

Из расчетных распределений на Рис. 3 для камер модельного и натурного ЖРД следует, что плотность тока уноса, определяющая напряженность собственного магнитного поля ПС, увеличивается по мере движения газового

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

Наибольшее значение уровня тока 1у натурного двигателя (РД-171М) значительно превышает уровень тока модельного ЖРД (например, на срезе сопла I/ = 4170 мА и 117 мА соответственно). Это объясняется значительными различиями в расходах топлива (598,25 кг/с для РД-171М и «0,5 кг/с для модельного ЖРД ). Однако рекомбинация заряженных частиц на срезе натурного ЖРД с большим расходом топлива также интенсифицируется. Вследствие этого расчетные значение уровня тока уноса на срезе сопла натурного двигателя больше модельного ЖРД почти в 18 раз.

Для верификации разработанной математической модели дополнительно проводилось сравнение с опубликованными экспериментальными данными других авторов, например, по определению электропроводности ПС в камере сгорания РД (Рис. 5, а) (Пинчук В. А., Грибакин В. А. и др.) и электрического потенциала сопла ЖРД (Рис. 5, б) (Ягодников Д. А., Воронецкий A.B., Пушкин Н. М.). Экспериментальные данные по определению потенциала сопла ЖРД были получены при параметрах: рк = 1 МПа, dKp =10 мм, da = 18 мм, горючее -этиловый спирт, окислитель - кислород (газ). Расчетные значения потенциала сопла (Рис. 5, б) определялись из условия дрейфа электронов в стенку, что придавало элементу конструкции суммарный отрицательный заряд.

а-П^Ом-'-м-1

-oir

• • > а)

• •

2 1

b ▲ •

v :

0,8

1,2

1,4 а.

Рис. 5.

Зависимости ст(а) в КС РД (а): 1 - сечение №1; 2- сечение №2 и <рс(а) ЖРД (б): 1 - рк=1 МПа; 2 - рк=0,5 МПа. Значки - эксперимент Максимальные расчетные значения электропроводности ПС в сечениях камеры РД№1 и №2 составили =4,49-10"3 Ом"' м"' и а№2 = 4,62-1 (Г3 Ом-1-м-1. Этим значениям соответствует соотношение компонентов близкое к 1 (а ~ 1), что обусловлено наибольшей температурой ПС 'ГК№1 = 2102 К (сечение №1), ТК№2 = 2107 К (сечение №2) и подтверждает наибольшие расчетные значения концентраций электронов в этой области N¡¡№1 = 5,5-1013 м"3 и N^2 = 5,8-1013 м"3.

Максимальное значение <рс = —210 мВ соответствует области а ~ 1 (Рис. 5, б, кривая 1), при котором достигается максимальная температура (Тг = 3614 К), что обусловливает максимальную расчетную концентрацию заряженных частиц вследствие термической ионизации (ЪГг = 2,569-1019 м"3; 1Чнзо+ = 9,835-1018 м"3).

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

В третьей главе приведены описания экспериментальной установки, систем измерения и регистрации, а также методика проведения огневых стендовых испытаний, основными задачами которых являются экспериментальное подтверждение прогнозируемых электрофизических характеристик ПС ЖРД и их изменение в процессе варьирования режимных параметров (давления в КС и соотношения компонентов топлива газообразный кислород и 75% - ый раствор этанола).

Конструкция камеры модельного ЖРД состоит из нескольких секций, закрепленных вертикально на стапеле (Рис. 6). Выбор данного двигателя обусловлен простотой его конструкции, практически неограниченным ресурсом работы, безопасностью эксплуатации и возможностью реализации различных режимов истечения из сопла. В серии экспериментов по определению электромагнитных свойств ПС в процессе разгара элементов проточного тракта двигателя на присоединительный фланец последней секции КС 1 монтировался

неохлаждаемый вкладыш 4, из композитных материалов типа стеклопластика П5-7, или углепластика П5-13.

Для регистрации

электромагнитного поля выхлопной струи ПС ЖРД использовались датчики магнитного 2 и электрического 3 полей (ДМП и ДЭП), которые позволяли в процессе эксперимента регистрировать

переменные составляющие

магнитного и электрического полей в частотных диапазонах 10 Гц...20 кГц. Сигнал с первичных преобразователей поступал на промежуточные усилители.

Информация регистрировалась и обрабатывалась с помощью средств микропроцессорной и вычислительной техники с использованием аппаратно-программного комплекса «ПИРС».

Пространственное положение ДМП и ДЭП по отношению к потоку ПС выбиралось таким, чтобы регистрируемый ими сигнал был наибольшим. Это

Рис. 6.

Конструкция модельного ЖРД с охлаждаемым (а) и неохлаждаемым вкладышем критического сечения (б)

достигалось с учетом данных предварительных градуировок датчиков, а также на основе электрофизической картины истечения ПС из сопла.

Осциллограммы на Рис. 7 были получены при моделировании разгара критического сечения и сверхзвуковой части сопла. Уменьшение давления в КС в диапазоне рк = 2,2... 1,5 МПа происходило за время х = 4... 12 с, что сопровождалось увеличением площади критического сечения в интервале 2,4-10-4...3-1<Г1м2.

рю МПа

,0006 - 0,0004 0,0002

Т, С

Рис. 7.

Осциллограммы пуска: 1 - рк, 2 - Р,ф, 3 — Н при разгаре критического сечения сопла модельного ЖРД, материал сопла П5-13 Из осциллограммы пуска (Рис. 7) следует, что по мере уменьшения давления в КС происходит монотонное уменьшение амплитуды напряженности магнитного поля, регистрируемой измерительным преобразователем магнитного поля, с Н = 40 А/м (рк = 2,2 МПа) до Н = 25 А/м (рк= 1,4 МПа).

Анализ результатов серии ОСИ при различных давлениях в КС позволил установить зависимость амплитуды напряженности магнитного поля от давления в КС (Рис. 8), которая в диапазоне рк = 2,5...3,5 МПа аппроксимировалась линейной функцией Н(рк) = 27,66-рк-25,86.

н, А/м

3 Рк, МПа

Рис. 8.

Зависимость амплитуды Н от давления в КС при разгаре критического сечения модельного ЖРД. Значки-эксперимент

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

характеристики электромагнитного поля истекающей струи продуктов сгорания, пропорционально

произведению значений диаметра среза сопла, скорости потока и концентрацией заряженных частиц в этом сечении, т.е. 1у ~ с1а-ша-Ые,1. Тогда при постоянном диаметре среза сопла

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

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

При испытании с разгорающимся критическим сечением интенсивно изменялся спектр амплитуд электрического поля по времени (Рис. 9).

Рис. 9.

Трехмерные графики амплитудных спектров сигналов с ДМП (а) и ДЭП (б) при I разгаре сопла (время анализа 5.. .6 с). Материал вкладыша П5-13

Из Рис. 9, а наглядно видно, что максимальные амплитудные спектры сигнала магнитного поля сосредоточены в ограниченном диапазоне частот 2900...3200 Гц. В отличие от амплитудного спектра магнитного поля, максимальные амплитуды спектральных плотностей мощности сигнала электрического поля находятся в более широком диапазоне 10...900 Гц низкочастотной области. Так, в процессе разгара критического сечения максимум плотности мощности смещается в коротковолновую область спектра,

где соответствует частоте £пах = 9241 Гц, а в случае постоянного давления в КС Гтах = 41,6 Гц (Рис. 10).

Спектральные характеристики сигнала ДМП при изменении давления в КС рк изменяются менее существенно и плотности мощности (Сспм) сигнала соответствует частоте Гшах ~ 3 кГц (Рис. 11). Однако плотность мощности сигнала ДМП при уменьшении давления в КС уменьшается почти в 10 раз при уменьшении давления, вызванного разгаром критического сечения.

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

вспм'Ю7, В2/Гц

1 1 1

1 и у ш !

Рис. 10.

Спектр плотности мощности сигнала ДЭП при постоянном значении рк(1)и рк=уаг (2). Материал вставки критического сечения П5-13

С другой стороны, продукты горения материала вкладыша критического сечения (углепластика П5 -13 ), электризуясь, непосредственно

изменяют параметры магнитного и электрического поля.

В четвертой главе рассмотрена структурная схема системы аварийной защиты (САЗ) ЖРД на основе регистрации амплитуды

напряженности магнитного поля (МП), обладающей большим быстродействием по сравнению с традиционными средствами (40... 50 мс). Выполнено моделирование системы камера ЖРД - первичный измерительный преобразователь магнитного поля высокотемпературного потока ПС углеводородного топлива и определено время инерционности (быстродействия) при скачкообразном изменении давления в камере сгорания (КС) и диаметра критического сечения сопла. В результате параметрического анализа установлено, что на конечное быстродействие такой системы главным образом влияет объем КС диагностируемого двигателя, конструктивные особенности измерительного преобразователя и параметры колебательного контура. Временные характеристики переходного процесса предлагаемой САЗ позволяют рекомендовать электрофизические характеристики как источник информации о работоспособности ЖРД.

Описание динамического процесса возникновения сигнала в измерительном преобразователе МП в результате изменения давления в камере двигателя проводилось с помощью составления структурной схемы «двигатель -измерительный преобразователь МП» в программе «МВТУ 3.5». Схема включала в себя задающее воздействие, которое описывалось ступенчатым законом изменения площади критического сечения, передаточную функцию камеры двигателя и измерительного преобразователя МП. Результатами программных вычислений являлись временные графики, описывающие переходные процессы в звеньях системы (Рис. 12).

По результатам проведенных расчетов можно заключить, что ступенчатое увеличение площади критического сечения вызывает монотонное падение давления в КС за время 14,2 мс. Длительность переходного процесса в измерительном преобразователе магнитного поля составило 17,6 мс, что свидетельствует о слабой инерционности последнего при мгновенном

Оспм'Ю7, В2/Гц

0 1 2 3 4 5 6 7 8 9 £ кГц

Рис. 11. Спектр плотности мощности сигнала ДМП при постоянном значении рк (1) и при падении рк в результате разгара критического сечения. Материал П5-13

. F ^ a) _ изменении площади критического

м»-1 "" —I_ сечения и подтверждает возможность

В)

к использования САЗ на основе

_ -Ц"2"' контроля электрофизических

х,с 105 характеристик ПС. — [дд Основные результаты и выводы

1. Разработана осесимметричная

. математическая модель электрофизических процессов в камере ЖРД на углеводородном Рис. 12. топливе, учитывающая особенности

Изменение давления (б), напряжения К01П7ра камеры

в измерительном преобразователе и термодинамических свойств

МП (в) и амплитуды напряженности ионшированных продуктов старания при МП (г) при увеличении площади различныхначальныхрежимныхшрамехрах критического сечения (а) и геометрии проточного тракта ЖРД

2. Получены расчетные зависимости электрических потенциалов ПС и конструкции сопла двигателя, уровни напряженносгей магнитного поля и электрической проводимости ПС в диапазоне изменения коэффициента избытка окислителя а»=0,58.. .2 и давления в КС Рк=3 МПа для топлива кислород-этанол. Для топлива кислород-керосин камеры маршевого ЖРД РД-171М в диапазонах а»=0,5.. .2 и р*=24,5 МПа.

3. Разработана методика проведения эксперимента, система регистрации и обработки результатов измерений характеристик электромагнит! юга поля при огневых испытаниях модельного ЖРД на топливе кислород-этанол, при моделировании аномальных режимов работы и процессов разгара сопла.

4. В диапазоне изменения давления в камере сгорания р* = 1,5...3,5 МПа экспериментально определена близкая к линейной зависимость напряженности магнитного поля от давления в камере сгорания, которая может бьпь использована как входной сигнал и диагностический признак для системы аварийной зашщы ЖРД

5. В результате частотного анализа сигнала, полученного с измерительного преобразователя магнитного поля, установлено, что в процессе уменьшения давления в диапазоне р, = 2^.. .1,5 МПа, вызванного разгаром критического сечения сопла, максимум плотности мощности сигнала магнитного поля, соответствующий частоте 3 кГц, уменьшается с 8-10*В2сдо lO^B^c.

6. Разработана структурная схема системы аварийной зашипи двигателя на основе регистрации амплитуды напряженности магнитного поля, обладающей большим быстродействием (17,6 мс) по сравнению с традиционными средствами (40.. .50 мс).

7. Выполнено моделирование системы «камера двигателя - измерительный преобразователь магнитного поля» и определено время инерционности (17,6 мс) при скачкообразном юменении давления в камере сгорания и площади критического сечения сопла. В результате параметрического анализа установлено, что на быстродействие такой системы главным образом влияет объем камеры сгорания диагностируемого двигателя, конструктивные особенности измерительного преобразователя и параметры колебательного контура.

Список публикаций

1. Ягодников Д. А., Рудинский А. В. Моделирование электрофизических процессов в газовых трактах энергосиловых установок // Вестник МГТУ имени Н. Э. Баумана. Сер. Машиностроение. 2010. №4(81). С.39-46. (0,2 пл./ 0,1 пл.).

2. Ягодников Д. А., Бобров А. Н., Рудинский А. В. Частотный анализ электрофизических характеристик рабочего процесса ЖРД на углеводородном топливе // Наука и образование. МГТУ им. Н. Э. Баумана. Электрон, журн. 2011. № 11. С. 1-26. http://technomag.edu.ru/doc/250245.html. (1,625 пл./ 0,5 пл.).

3. Ягодников Д. А., Рудинский А. В. Моделирование неравномерного течения и параметров электромагнитного поля в камере жидкостного ракетного двигателя при наличии пристеночного слоя // Наука и образование. МГТУ им. Н. Э. Баумана. Электрон, журн. 2015. № 4. С. 141-154. http:www.technomag.edu.ru/issue/763210.html. (0,875 пл./ 0,44 пл.).

4. Ягодников Д. А., Рудинский А. В. Моделирование электрофизических процессов в камере жидкостного ракетного двигателя // Ракетно-космические двигательные установки: Сборник материалов Всероссийской науч.-тех. конф. М., 2010, С. 51-52. (0,125 пл./ 0,07 пл.).

5. Ягодников Д. А., Рудинский А. В. Электрофизические процессы в ЖРД // Студенческая научная весна-2010: Сб. тез. докладов. НТК. М., 2010.Т.10. ч.З. С. 219-220. (0,125 пл./ 0,07 пл.).

6. Ягодников Д. А., Рудинский А. В. Исследование характеристик магнитных элементов для диагностики процессов разгара теплонапряженных элементов конструкции жидкостных ракетных двигателей // XXXVI Академические чтения по космонавтике 2012, Секция 3. С. 76-77. (0,127 пл./ 0,07 пл.).

7. Ягодников Д. А., Рудинский А. В. Исследование частотных характеристик электромагнитного поля сверхзвуковой струи ЖРД // Вестник МГТУ имени Н. Э. Баумана. Сер. Машиностроение. Спецвыпуск №7. Теория и практика энергетического и транспортного машиностроения. 2012. С.12-18. (0,44 пл./ 0,22 шт.).

8. Ягодников Д. А., Рудинский А. В. Оценка параметров быстродействия САЗ ЖРД по электромагнитным свойствам продуктов сгорания // Ракетно-космические двигательные установки: сб. материалов Всероссийской научно-технической конференции. 2013. С. 88-89. (0,125 пл./ 0,07 пл.).

9. Ягодников Д. А., Рудинский А. В. Экспериментально-теоретическое определение параметров быстродействия САЗ ЖРД по электромагнитным свойствам продуктов сгорания // Труды МГТУ имени Н.Э. Баумана. Теория и практика современного ракетного двигателестроения. 2013. №607. С.24-31. (0,5 пл./ 0,25 пл.).

10. Рудинский A.B. Экспериментально-теоретическое исследование магнитометрического метода контроля и диагностики рабочего процесса ЖРД // XX научно-технической конференции молодых ученых и специалистов РКК «Энергия». Тезисы докладов. 2014. С.424-425. (0,125 пл.).

Подписано к печати 16.06.15. Заказ № 367 Объем 1 печ.л. Тираж 100 экз. Типография МГТУ им. Н.Э. Баумана 105005, Москва, 2-я Бауманская ул., д.5,стр.1 (499) 263-62-01