CC BY
57
DOI: 10.15593/2224-9982/2015.41.01 УДК 662.2:658.567.5
М.В. Попова1, А.В. Литвинов1, С.Н. Козлов2
ОАО ФНПЦ «Алтай», Бийск, Россия 2 Бийский технологический институт (филиал) Алтайского государственного технического университета им. И.И. Ползунова, Бийск, Россия
СИСТЕМНАЯ ФОРМАЛИЗАЦИЯ И ИДЕНТИФИКАЦИЯ ПРОЦЕССОВ ПРИ ОГНЕВОЙ УТИЛИЗАЦИИ КРУПНОГАБАРИТНЫХ ЗАРЯДОВ РДТТ НА ОТКРЫТОМ СТЕНДЕ
Применительно к структуре информационно-управляющей системы «Экология» рассмотрены вопросы физико-математической формализации газодинамических процессов, протекающих в струе и облаке продуктов сгорания твердотопливных зарядов при их огневой утилизации на открытом стенде с применением водной экологической защиты. По мнению авторов, сложная картина и многофакторность процессов не позволяет достоверно формализовать их детерминированными методами, поэтому предлагается информационные потоки по параметрам процессов разделить по аналогии с радиотехникой на «сигнал» и «шум», используя для описания сигнала инженерные модели струи и облака продуктов сгорания, а все виды неопределенностей вынести в «шумовую» область, величину которой определять статистической обработкой результатов сжиганий зарядов и использовать эти данные при оценке экологического риска. Исходя из этого в статье предложены аналитические зависимости для определения параметров струи и облака продуктов сгорания, выбраны контрольные параметры и разработаны схемы идентификации этих зависимостей по результатам статических сжиганий заряда, исследована чувствительность выбранных моделей к управляющему воздействию - расходной характеристике водной экологической защиты по воде, даны определения экологических рисков. Показано, что выбранные аналитические модели струи и облака продуктов сгорания отвечают требованиям управляемости, идентифицируемости, адаптивности и могут использоваться в структуре информационно-управляющей системы «Экология».
Ключевые слова: экологическая безопасность, экологический риск, статическое сжигание, огневая утилизация, физико-математическая модель, многофакторные процессы, энергетическая установка, информационно-управляющая система, системный анализ, техническая система, проектирование информационных систем.
M.V. Popova1, A.V. Litvinov1, S.N. Kozlov2
OJSC "Altay", Biysk, Russian Federation 2 Biysk Technological Institute, Biysk, Russian Federation
SYSTEMIC FORMALIZATION AND IDENTIFICATION OF PROCESSES AT FIRING UTILIZATION OF LARGE-SIZE SOLID ROCKET MOTOR CHARGES AT AN OPEN TEST BENCH
Issues of physico-mathematical formalization of gas-dynamic processes occurring in jet and cloud of combustion products of solid propellant charges on their firing utilization at an open test bench equipped with water ecological protection have been considered as applied to the structure of the information-controlling system Ecology. In the authors' opinion the complicated picture and multi-factor processes do not make it possible to formalize them reliably by determined methods. Therefore, they suggested dividing data flows into process parameters as signal and noise, as accepted in radioengineering, using for signal description engineering models of a jet and a cloud of combustion products and transferring all kinds of uncertainties into the noise field, its value to be determined by static treatment of results of charge combustions for evaluating ecological risk. Based on this, analytical relationships to define parameters of a jet and a cloud of combustion products have been proposed in the paper. Control parameters have been selected and schemes to identify these relationships on the base of static firing results have been developed. Response of the models selected to the control effect - water discharge characteristic of water ecological protection has been studied and definitions have been given. The selected analytical models of combustion products jet and cloud are shown to meet the requirements of controllability, identifiability, adaptability and can be used in the structure of the information-controlling system Ecology.
Keywords: ecological safety, ecological risk, static firing, firing utilization, physic-mathematical model, multi-factor processes, power plant, information-controlling system, systems analysis, technical system, information systems design.
В статье [1] рассмотрена методология создания информационно-управляющей системы под задачи огневой утилизации твердотопливных зарядов крупногабаритных ракетных двигателей на открытом стенде при действии водной экологической защиты (ВЭЗ). По разработанной ФНПЦ «Алтай» технологии статических сжиганий (СтС) предусмотрено сжигание твердотопливного заряда в двигателе без соплового блока, при этом в процессе горения заряда струя продуктов сгорания (ПС) орошается водой из коллекторов (рис. 1), расположенных вдоль струи ПС, что позволяет осадить большую часть вредных компонентов (ВК) ПС в рабочей зоне стенда.
Структурно схема контроля и обеспечения экологической безопасности (ЭБ) стенда представлена как управляемая техническая система открытого типа с регулятором непрямого действия, в которой регулируемым параметром является концентрация ВК ПС на границе сани-тарно-защитной зоны (СЗЗ), а регулирующим элементом - расходная характеристика ВЭЗ по воде.
Применительно к рассматриваемой задаче такой подход, на наш взгляд, является единственно возможным, поскольку:
- объединяет все процессы, а также входные и выходные параметры в единую систему и позволяет рассматривать задачу контроля и обеспечения ЭБ в рамках теории технических и информационных систем;
- позволяет объединить детерминированные и стохастические методы для формирования и последующего параметрического исследования системы.
Действительно, основу блока управления информационно-управляющей системы (ИУС) составляют концептуальные логические модели, которые разрабатываются на основе физико-математических моделей процессов, протекающих в системе. При этом физико-математические модели процессов должны обеспечить основные характеристики ИУС -управляемость, идентифицируемость и адаптивность.
Под управляемостью понимается способность системы реагировать на изменение ее входных параметров. При этом качество и эффективность управления зависит от правильности выбора параметров, используемых для управления системой.
Система ВЭЗ
Рис. 1. Схема подачи воды в коллекторы
Под идентификацией понимается установление соответствия между результатами математического моделирования и результатами экспериментального определения контролируемых параметров системы. При этом, если вход X и выход У системы объединены оператором соответствия А: У = АХ, то задача идентификации состоит в определении оператора А по результатам сравнения расчетных и опытных данных.
Под адаптацией понимается способность системы к повышению качества управления в процессе испытаний или эксплуатации за счет уточнения, по мере набора статистических данных по параметрам оператора А, и повышения точности идентификации.
Используемые в ИУС «Экология» функциональные модели, разрабатываемые на основе детерминированных математических моделей регулятора и распределителя, должны обеспечивать все перечисленные свойства ИУС и включаться в систему только после проверки на соответствие этим требованиям. Поэтому выбор физико-математических моделей процессов играет ключевую роль в формировании блока управления и обеспечения эффективности ИУС.
Основными в системе ЭБ являются:
- газодинамические процессы в струе ПС, орошаемой водой;
- газотермодинамические процессы, определяющие состояние облака ПС, которое после орошения струи поднимается, расширяется и рассеивается в атмосфере.
В обоих случаях состояние и струи, и облака ПС обусловливается взаимодействием многих процессов, которые практически невозможно описать в рамках каких-либо существующих детерминированных подходов. Использование численных методов [2] не помогает: они позволяют построить многопараметрические модели, обеспечить любую требуемую вычислительную точность, но адекватность их действительно существующему описываемому моделью процессу зависит от достоверности исходных посылок, которые не всегда можно корректно сформулировать. Кроме того, высокая трудоемкость вычислений ограничивает их возможность для параметрических исследований системы и ее численной оптимизации. С учетом этого применительно к формированию блока управления ИУС целесообразно рассматривать инженерные упрощенные модели, позволяющие выделить основные действующие процессы, а второстепенные факторы и различные виды неоп-
ределенностей вынести в «шумовую зону», т.е., по аналогии с радиотехническими системами, разделить информационный поток по параметрам процессов на сигнал и шум. Это позволит формализовать концептуальные и логические модели блока управления ИУС и обеспечить все основные показатели этой системы: управляемость, идентифицируемость, адаптивность. Рассмотрим этот методический подход для каждой из подсистем ИУС - регулятора и распределителя.
1. Регулятор: орошаемая струя ПС.
Исходной базой для построения математической модели регулятора является структура струи ПС, образующейся при горении заряда, и газотермодинамические процессы, протекающие в струе при ее орошении водой.
Физика процессов в свободно истекающих высокоскоростных, высокотемпературных струях изучена достаточно хорошо в работах Г.Н. Абрамовича [3], И.П. Гинзбурга [4] и др. Однако при проведении СтС в двигателе без соплового блока с орошением струи ПС водой через коллекторы, расположенные вдоль струи, существенно меняется газодинамическая картина струи ПС, при этом необходимо учитывать следующие факторы:
1. Сжигание заряда в корпусе двигателя без соплового блока меняет газодинамическую картину течения ПС внутри заряда [5] - в начальный период горения критическое сечение для ПС располагается внутри канала заряда и по мере его разгорания перемещается в сторону выходного сечения корпуса. Этот процесс необходимо учитывать при расчетном определении массового расхода ПС по результатам замера давления в двигателе, которое на начальном участке разгара канала заряда снижается с 1,0 до 0,2 МПа и фиксируется на этом уровне при стабилизации критического сечения в выходном сечении корпуса.
Поэтому при разработке алгоритмов управления целесообразно рассматривать два участка горения заряда (рис. 2):
1) высоконапорный (ВНУ) - соответствующий периоду разгара внутреннего канала заряда, на котором давление снижается до минимального;
2) низконапорный (ННУ) - в течение которого критическое течение устанавливается в выходном сечении корпуса и давление в двигателе поддерживается на уровне 0,20-0,25 МПа.
Р. МПа
О 60 120 180 240 300
т, с
Рис. 2. Диаграммы давления (избыточного) в ЭУ-1, 2
2. Горение заряда в двигателе без соплового блока протекает при пониженном, относительно штатного режима, давлении в камере сгорания (0,2-1,0 МПа вместо 7-10 МПа в штатных двигателях), что приводит к значительному снижению полноты сгорания топлива и, соответственно, к снижению температуры ПС. Этот факт подтвержден экспериментально [6, 7] и должен учитываться при расчете температуры струи ПС.
3. Как показывает анализ, в течение всего времени горения заряда при СтС струя ПС истекает из двигателя в надкритическом, т.е. сверхзвуковом, режиме и переходит в дозвуковой режим течения через серию скачков уплотнения, которые образуют переходный участок - «бочку». Поскольку ввод воды в струю на сверхзвуковом участке мало эффективен при давлении подачи воды 0,5-1,0 МПа, то при разработке конструкции и схемы расположения коллекторов системы орошения необходимо учитывать газодинамическую структуру струи и проводить расчетное обоснование оптимального их расположения по длине струи.
4. Энергия свободно истекающих струй расходуется на работу расширения и увеличение кинетической энергии струи [8]. При орошении струи водой возникают затраты энергии на массообмен [9], что снижает энтальпию и кинетическую энергию струи, приводит к снижению ее температуры и скорости, однако при этом расходы энергии на работу расширения не исчезают. Эту особенность энергетического баланса струи ПС необходимо учитывать при разработке расчетного алгоритма.
5. При орошении струи возникает сложная картина взаимодействия воды и струи ПС. Условно можно нарисовать следующую виртуальную схему: при подаче воды в струю ПС через кольцевые коллекторы, расположенные по внешнему контуру струи, создается водяной «чулок» с тремя слоевыми зонами:
- в 1-й зоне вода проникает в струю, испаряется и взаимодействует с компонентами ПС. Толщина к этой зоны и, соответственно, степень орошения струи зависят от соотношения гидродинамических напоров струи ПС и воды и увеличиваются по мере торможения струи ПС;
- во 2-й зоне реализуется только поверхностный контакт между водой и струей ПС, вследствие этого вода, не проникая в струю ПС, нагревается и испаряется;
- в 3-й зоне вода не контактирует со струей ПС и полностью сливается на пол огневого двора.
Эти зоны можно рассматривать чисто условно, поскольку значительную роль в протекании процессов играет динамика истечения струи. Действительно, длина огневого двора составляет ~ 100 м. На этом участке скорость течения струи меняется от ~ 900 до ~ 100 м/с, а температура от 1900 до 400 К, т. е. при средней скорости течения струи ~ 500 м/с ПС преодолевают это расстояние за 0,2 с с градиентом по скорости течения струи ~ 4-10 м/с , по температуре ~ 7,5-10 К/с. Это должно приводить к неравновесности и незавершенности практически всех физико-химических процессов взаимодействия воды и ПС в струе [10].
С учетом перечисленных факторов построение базовой физико-математической модели регулятора включает в себя:
- расчет давления и температуры в камере сгорания двигателя;
- расчет газотермодинамических параметров орошаемой и неорошаемой струй по их длине;
- покомпонентный расчет массовых характеристик струи ПС при действии ВЭЗ.
Расчет внутрибаллистических параметров внутри двигателя проводится в следующем порядке:
1. С использованием формулы Бори, исходя из поверхности горения с учетом разгара канала заряда, рассчитывается диаграмма изменения давления во времени в камере двигателя:
Р СО =
( Я
V ^ ,
где Я(е) - площадь поверхности заряда в функции свода; р - плотность топлива; А - коэффициент истечения; ^кр - площадь критического сечения, определяемого с учетом соотношения площадей поверхности горения и проходного сечения канала заряда; V - показатель зависимости скорости горения от давления и = и1Ру.
Полученная расчетная функция Рк (т) сравнивается с осциллограммой давления, полученной при СтС.
2. Для диаграммы Рк (т) с использованием температурных коэффициентов Кт табл. 1 определяется температура Т (т.. ) в камере сгорания:
Т0. (т. , Р ) = Т0Кт (Рко ) . Температурный коэффициент
Таблица 1
Рк, кгс/см2 40 20 14 10 8 6 4 2 1
Кт = Т } т Т 1о 1,00 0,98 0,96 0,94 0,91 0,88 0,84 0,70 0,62
Величина коэффициента Кт определена по результатам работ [6, 7]. 3. Рассчитываются значения критической скорости звука по фор-
муле
а =
кр
\2kRT к +1
при этом ориентировочно могут быть приняты R = 291,30 кДж/кг-К, к = 1,14.. .1,20. Конкретно для каждого РДТТ эти величины уточняются по паспортным характеристикам топлива.
4. Определяются параметры ПС в выходном сечении для каждого .-го момента времени:
- скорость истечения Жа = Хакр;
- давление Ра = Ркп(ка);
1
- температура Та = т0х(Ха).
Здесь п(Ха ), т(Х а) - газодинамические функции [11]. 5. Массовый расход истекающих газов определяется по формуле
где А = ф2
к +1)
™ ПС = АРРкр,
е
, ф2 = 0,98.. .1,0, при этом ^кр для началь-
ного высоконапорного участка определяется исходя из геометрии разгорающегося канала заряда.
Полученные функции Ра (т), Та (т), Xа (т), Жа (т) служат исходными данными для расчета параметров струи ПС в зоне стенда.
Расчет параметров неорошаемой струи проводится с учетом следующих факторов.
Как показано ранее, при работе без соплового блока критическое сечение потока, в котором X = 1, устанавливается во внутреннем канале заряда и по мере его разгорания перемещается в сторону выходного сечения двигателя с диаметром ёа, где окончательно позиционируется. В любом случае в атмосферу выходит сверхзвуковая струя, которая через серию скачков уплотнения переходит в дозвуковую с давлением, равным рн, и далее изобарически расширяется с уменьшением скорости и температуры потока по длине струи вследствие затрат энергии на расширение. При этом диаметр струи увеличивается [3, 4].
Структура и форма начального участка определяются соотношением давлений на выходе из двигателя Ра и окружающей среды рн
: при п > 2 начальный участок состоит из одной «бочки», со-
( р!
I Рн )
держащей ударные волны, а при нерасчетностях, близких к единице, будет содержать несколько «бочек» [4], при этом
- максимальный диаметр «бочки» определяется по формуле
^^ = 1,244п°'6;
- длину начального газодинамического участка можно оценить по формуле
I «-а-
0 1,2 + 0,3 и
Понимание структуры начального участка истечения ПС важно для проведения системного анализа и расчета параметров струи в зоне расположения водных коллекторов ВЭЗ, которые целесообразно располагать на участке струи после ее перехода в дозвуковой режим течения.
При расчете неорошаемой струи скорость и температура на переходном участке условно принимаются постоянными и равными Шхп = = ш т = т
'' а 1хп 1 а-
В качестве примера в табл. 2 и на рис. 3 приведены форма и размеры переходных участков струи для высоконапорного и низконапорного участков.
Таблица 2
Размеры переходного участка струи ПС в зависимости от давления в камере
Рк, кг/см2 Ра, кг/см2 dmax, м 1, м
30 15,7 2,5 4,2
20 10,5 1,9 3,8
10 5,2 1,3 2,9
5 2,619 0,8 2,1
2 1,048 0,5 1,1
у
у
Рис. 3. Расположение и схемы переходного участка струи ПС: а - высоконапорный участок; б - низконапорный участок
Расчет изменения скорости и температуры ПС по оси струи после перехода ее в дозвуковой режим проводится по следующим зависимостям [3, 4]:
Ж. = Ж ^^
х
где Жх - скорость потока ПС на оси струи; х - расстояние от координаты диска Маха; ёа - радиус выходного сечения.
Т = Т +(Т - Т) 5,45 ' ёа ,
х н \ а н / х '
где Тх - температура потока ПС на оси струи; Тн - температура окружающей среды.
Приведенные зависимости получены аппроксимацией экспериментальных данных и учитывают работу изоэнтропического расширения потока [3, 12].
Параметры струи при орошении могут быть определены следующим образом [13]. При истечении расширяющейся струи ПС энергия потока содержит две составляющие (тепловую и кинетическую) и может быть записана в виде
2 2
Е = г + А— = С Т + А—, 2^ р 2
где г - энтальпия потока; А - тепловой эквивалент механической работы, А = кДж/кг ; Ср - теплоемкость ПС при постоянном давлении; 102
Т - температура; ч - скорость потока.
При орошении вода вводится в струю через последовательно расположенные вдоль струи кольцевые коллекторы, вследствие чего образуется термодинамическое тело переменной массы [9], при этом энергия потока расходуется дополнительно к работе расширения на массообмен между водой и ПС и на испарение воды, поступающей в струю ПС. Если принять допущение о мгновенности взаимодействия, то параметры струи ПС в каждом сечении струи, соответствующем расположению водяного коллектора, будут определяться следующими законами газодинамики свободных струй: - закон сохранения массы в виде
Ч(п+1)= ™г(п+1)= тгп + ™вп, (1)
где дагп - секундный массовый расход газов через сечение впрыска п-й ступени коллектора; шъп - секундный расход воды при впрыске в п-й
ступени коллектора; т^^) - секундный расход газов после впрыска
воды в п-й ступени коллектора;
- равнение количества движения в виде
тШх = тг(п+1)Шп+1 = (™гп + ™вп )Ш
п+1'
(2)
где Шх - скорость движения ПС на оси струи; шх - расход ПС на оси струи;
- закон сохранения энергии в виде
(
т
п+1)
С Т + А • ^г(п+1)
С рг( п+1)Тг( п+1)+ А 2 ^
\
(
= т„
СргпТгп + А
Ш
2Л
2 Я
-т„
С 1 шв Л
Срвп АГн + г +1 • А • ^
V Р н 2 2^,
Здесь в соответствии с зависимостью (1) из энтальпии газа перед впрыском воды вычитается энтальпия, необходимая для нагрева воды до температуры кипения и парообразования. В этой формуле приняты следующие обозначения: А - тепловой эквивалент механической работы, А = 1012 кДж/кг ; АТн - температура нагрева воды до точки кипения; г - скрытая теплота парообразования, г = 2200 кДж/кг.
Теплоемкость газов после впрыска воды определяется по формуле
С = С
рг (п+1) ргп
т„
- + С •
т„
т
п+1)
рв
т
п+1)
Тогда из формулы (2) определяется температура ПС после впрыска воды:
Т,
т„
(
п+1)
т г( п+1)Срг( п+1)
С
X
СргпТгп + А •
]
2 Я )
2 Л
т„
т
п+1)
"Срс( п+1) Х
срвАТн + г +1 А^
V рв н 2 2я )
При оценке массообмена между водой и ПС необходимо учитывать, что глубина проникновения воды в струю ПС зависит от соотношения гидравлических напоров струи воды и ПС; приближенно эта величина для каждого коллектора может быть определена по формуле
где Рв - давление подачи воды из коллектора; Рг - давление в струе ПС, Рг = Ркл(Х); М - число Маха для скорости потока в данном сечении; с1 - диаметр отверстий коллектора; к - коэффициент изоэнтропы;
С учетом этого коэффициент орошения струи для каждого /-го сечения впрыска воды
где Як - радиус коллектора по центрам расположения отверстий.
Тогда масса воды, израсходованная на орошение струи, приближенно
по всем / = 1...п коллекторам.
Приведенные зависимости позволяют определить все параметры струи ПС - Т, Ж, Р, исходя из расходной характеристики ВЭЗ Ов (х, т),
а также степень орошения струи по количеству воды т^р, проникшей в струю.
Эта величина, исходя из зависимости (3), должна увеличиваться по мере выгорания заряда, снижения давления в камере двигателя и по длине струи ПС по мере ее торможения.
(3)
ф= - коэффициент напора.
(я - Ь)2
Я 2
С учетом рассмотренной ранее схемы взаимодействия воды со струей ПС может быть определен баланс масс воды следующим образом.
Вода, подаваемая из коллекторов для орошения струи ПС, распределяется следующим образом:
- часть т™ - не входит в струю и напрямую сливается на пол;
- часть т^сп - испаряется при поверхностном контакте со струей ПС;
- часть т™ - входит внутрь струи, при этом часть ее т™'исп, захватив часть ПС, испаряется, а часть т™'ос, захватив ПС, осаждается
на пол огневого двора, т.е. в технологические стоки.
По длине струи соотношение между этими частями будет меняться.
Тогда баланс масс воды
т* = твсл + тисп + т.™ , тв = тв + тв + т ' + т ' .
в в в в в
Отсюда следует, что только часть подаваемой воды эффективно используется для осаждения ПС в рабочей зоне стенда. Если эта доля мала, то увеличение расхода воды неэффективно.
Этот фактор необходимо учитывать при оценке эффективности ВЭЗ и анализе чувствительности подсистемы регулятора при выборе управляющей функции.
Идентификация регулятора.
Для регулятора как подсистемы ИУС «Экология» должны быть определены контрольные параметры и порядок идентификации функциональной модели. При этом должны быть выбраны:
а) контрольные параметры регулятора, по которым проводится идентификация;
б) идентификационные характеристики для подгонки модели к идентифицируемому процессу.
Для идентификации могут быть использованы:
- результаты натурных СтС и модельных огневых испытаний;
- результаты расчетно-теоретических исследований процессов или их отдельных элементов на основе математически точных числен-
ных методов решения уравнений газовой динамики. Применительно к рассматриваемому вопросу в качестве таковой может использоваться методика ФНПЦ «Алтай» [2], позволяющая на основе численных методов проводить расчетные исследования термогазодинамических параметров струи, состоящей из смеси ПС, паров воды и атмосферного воздуха. Такая модель может рассматриваться как опорная модель в цилиндрической системе координат на основе укороченных уравнений Навье - Стокса. Методика позволяет исследовать взаимодействие ПС и паров воды и проводить параметрическую оценку газодинамической структуры орошаемой струи в любом ее сечении.
Функциональным назначением регулятора является осаждение ВК ПС в рабочей зоне стенда. Исходя из этого эффективность регулятора как подсистемы должна оцениваться по степени «очищения» струи ПС от ВК. Такую оценку позволяют проводить результаты аналитического контроля содержания Al2O3 и HCl в технологических стоках, собранных в бассейне-отстойнике. В свою очередь, степень осаждения ВК в рабочей зоне стенда коррелирована со степенью торможения струи ПС перед бассейном, вследствие чего перед бассейном уменьшаются ее скорость и температура Жк, Тк до конечных величин.
С учетом этого в качестве контрольных параметров регулятора целесообразно использовать:
- массу и долю осажденных в рабочей зоне стенда Al2O3 и HCl: т°С тОС
т°С т°С к ос - m°A к °С - тхв •
'ОА > '"ХВ > лОА _ £ > ЛХВ _ б 5 тОА тХВ
- температуру струи ПС на конечном участке струи перед бассейном Тк.
Использование температуры струи для контроля целесообразно, поскольку:
- существующие инструментальные средства позволяют измерять температуру с относительно небольшими погрешностями, а датчики могут быть размещены в нескольких точках по длине струи;
- температура функционально связана со скоростью потока, что позволяет косвенно оценивать и этот параметр струи ПС на ее конечном участке.
Эти же параметры используются в качестве контрольных и при идентификации функциональной модели, при этом температура Тк
может контролироваться как для орошаемой струи, так и для струи без орошения, что позволяет повысить информативность идентификации.
В качестве начального этапа проводится точечно-параметрическая идентификация функциональной модели регулятора в следующей последовательности (схема на рис. 4).
1. Идентификация функциональной модели неорошаемой струи ПС по конечной температуре. По результатам сравнения опытной и расчетной температур Тк в конце струи ПС корректируется начальная температура струи Т0 на величину Т0'.
Рис. 4. Схема идентификации
2. Исходя из откорректированной величины температуры ТО проводится расчет параметров струи ПС для заданной расходной характеристики Ов (х, т), определяется конечная температура Ткр' и сравнивается с опытным значением Т°п'. Принимается, что расхождение опытной и расчетной температур Т°п' < > Ткр' возникает вследствие
неопределенностей в оценке степени орошения струи. Подгонка модели проводится за счет корректировки глубины проникновения воды в струю ПС:
2P„
Фн
h(x, т) - d •
13 +-
Gm
л_
kM2 P( x)
путем изменения коэффициента напора фн. В результате определяются величины к'( х, т) и масса воды, вошедшая в струю ПС т™'.
3. Идентификация эффективности осаждения ВК проводится по балансу масс воды, пошедшей на орошение. В качестве исходных параметров используются:
- масса воды в бассейне т°С - по результатам замеров;
ОС
- масса оксида алюминия в технологических стоках тОА и доля
ОС
7 ОС т°А
его осаждения: к°А - —у2-
т,
- по результатам замеров в бассейне;
ОА
— масса воды, вошедшая в струю ПС, т™ .
Из структуры орошаемой струи следует, что часть воды, вошедшая в струю, испаряется, а другая часть обеспечивает осаждение ВК из струи: ту - т°С + т^. В частности, в работе [10] показано, что хлористый водород образует хлористую кислоту с выделением тепла, достаточного для ее закипания и выделения газообразной HCl даже на конечном участке струи при температуре 110-70 °С.
С учетом этого логично оценить массу воды, израсходованной на осаждение ВК, по степени осаждения Al2O3, образующего в конце струи твердые компоненты:
,ян,ОС в
т ' - т„
• к
ОС ОА ,
т°С
где к°С - т°А
где ^ОА - у
т,
ОА
Эти соотношения позволяют определить все компоненты массового баланса воды и ВК и их распределение между рабочей зоной стенда и атмосферой.
Чувствительность выхода ко входным параметрам определяет степень управляемости системы. Требуемые параметры расходной характеристики Gв(x, т) обеспечиваются:
- количеством и координатами размещения коллекторов;
- диаметром отверстий для воды;
- давлением подачи воды.
Кроме того, при анализе чувствительности необходимо учитывать, что количественно (а не функционально) степень зависимости между выходом и входом регулятора будет различной на высоко- и низконапорных участках горения заряда, что следует из зависимости (3).
Исходные алгоритмы, принятые для разработки функциональной модели регулятора, позволяют учесть все эти факторы, исследовать влияние на выходные параметры регулятора на основе численного моделирования зависимостей и проводить оптимизацию расходной характеристики Gв(x, т) путем соответствующего размещения коллекторов, выбора диаметров для подвода воды и давления подачи воды. На начальном этапе для такой оптимизации проводится однопараметриче-ское исследование влияния расхода воды на температуру и скорость струи на конечном этапе ее расширения после орошения и степень орошения струи.
Для иллюстрации возможностей принятой для ИУС «Экология» функциональной модели регулятора на рис. 5-9 приведены зависимости скоростей, температуры и степени орошения струи ПС от среднего секундного расхода воды.
Рис. 5. Изменение скорости струи ПС на ВНУ в зависимости от среднего расхода воды Ов(х, т): 1 - Ов = 0 т/с; 2 - Ов = 0,2 т/с; 3 - Ов = 0,4 т/с; 4 - Ов = 1,0 т/с; 5 - Ов = 2,0 т/с
поп
У 00 700
= 300'
2 100
§
-20
г 2 ;
х/
\ ^
\
: (7в(д,т) = \аг, т/с
! Участок размещения коллекторов! --1-1-1-'-1-I-
10
20
25
Рис. 6. Изменение скорости струи ПС на ННУ в зависимости от среднего расхода воды Ов(х, т): 1 - Ов = 0 т/с; 2 - Ов = 0,2 т/с; 3 - Ов = 0,4 т/с; 4 - Ов = 1,0 т/с; 5 - Ов = 2,0 т/с
0.14 0.13 0,12 0,1! 0.10 0,09 0,08 0.07 0.06
_____-2
3
-1-■-1-■- —-— 4 -1-1-
7.4 10
15
20
25
0,50 0,45 0,40 0,35 0.30 0,25 0,20 0,15 0.10 0,05
..............!......... — /
_2
т ......................... .... 2
- 4
-•- ---
7,4 10
15
20
25
Рис. 7. Глубина проникновения воды в струю ПС на ВНУ по длине струи в зависимости от среднего расхода воды Ов(х, т): 1 - Ов = 2,0 т/с; 2 - Ов = 1,0 т/с; 3 - Ов = 0,4 т/с; 4 - Ов = 0,2 т/с
Рис. 8. Глубина проникновения воды в струю ПС на ННУ по длине струи в зависимости от среднего расхода воды
Ов(х, т): 1 - Ов = 2,0 т/с; 2 - Ов = 1,0 т/с; 3 - Ов = 0,4 т/с; 4 - Ов = 0,2 т/с
Рис. 9. Степень орошения струи ПС в зависимости от среднего расхода воды Ов(х, т)
Из графиков, в частности, видно, что чувствительность всех параметров струи к изменению расхода воды существенно выше для низконапорного участка горения заряда, чем для высоконапорного.
Исследование зависимости степени осаждения ВК в рабочей зоне стенда от расхода воды представляет определенные трудности, поскольку проведение натурных СтС при различных расходах воды затруднительно вследствие их высокой стоимости. В этом случае может быть использована следующая схема оценки зависимости
тОС = /(Gв (х, т)) , где I - ВК и вода. Если определена степень осаждения при некотором расходе О0( х, т), например Gв(х, т) = 2 т/с, то эта
точка используется в качестве первой опорной точки для построения искомых зависимостей (рис. 10).
Вторую опорную точку дает СтС без орошения струи, при котором Gв(х, т) = 0 т/с и тОС = 0. Это дает две точки для зависимости
тОС (Ов (х, т)), которая между этими точками может быть (см. рис. 10):
- линейной;
- выпуклой;
- вогнутой.
Проведение СтС при одном из промежуточных значений Gв(х, т) < GB) (х, т) позиционирует вид этой функции.
Таким образом, выбранная функциональная модель регулятора отвечает требованиям системы по чувствительности и позволяет обеспечить ее идентификацию и управление по расходной характеристике ВЭЗ.
Адаптируемость. Воздействие большого числа второстепенных «шумовых» факторов приводит к разбросу выходных параметров регулятора Тк, Жк, тОС . Поэтому по мере проведения СтС и набора статистических данных проводится уточнение всех выходных параметров и их разбросов и, соответственно, уточнение зависимостей между расходной характеристикой ВЭЗ и выходными параметрами регулятора.
Распределитель: облако ПС при подъеме в атмосфере.
В рамках подсистемы распределителя рассматривается облако ПС, которое формируется в процессе горения заряда и представляет собой аэрозольный термик, содержащий ВК, не осажденные в рабочей зоне стенда, пары воды и хлорной кислоты, содержание которых зависит от эффективности орошения струи ПС. Оксид алюминия в составе облака существует в виде твердых частиц k-фазы, размеры которых при среднем диаметре dk = 5,5 мкм подчиняются гамма-распределению с параметрами [13, 14]:
f (d) = 0,578d°'68 • е"0'68.
При разработке функциональной модели облака ПС необходимо учитывать следующие факторы:
1. Двухрежимность горения заряда при постоянной расходной характеристике ВЭЗ Ов (x, т) = const обусловливает слоистую структуру облака ПС. Физически это вызывается резким увеличением степени орошения струи по мере понижения давления в камере сгорания и соответственно расхода ПС.
Слоистость структуры облака ПС хорошо видна из рис. 11.
Рис. 11. Вид облака ПС: 1 - начальный участок; 2 - осадки; 3 - расширение и рассеивание
2. При подъеме термика изменяются его масса, габариты, температура и компонентный состав, которые зависят от следующих процес-
сов: эжекции атмосферного воздуха в облако при его подъеме; изменения агрегатного состояния его компонентов, в том числе конденсации паров воды и хлористой кислоты; формирования и выпадения (седиментации) аэрозолей и твердых частиц, при этом в качестве ядер конденсации могут служить твердые частицы &-фазы оксида алюминия.
Однако часть частиц оксида алюминия вымывается при орошении струи, а при подъеме облака может происходить их коагуляция
[14]. Эти процессы не исследованы ни экспериментально, ни теоретически, вследствие чего достоверный учет этого фактора при анализе конденсации и седиментации компонентов ПС практически невозможен.
3. Одним из важных факторов, которые могут определять процессы седиментации ВК из облака, является пульсационный характер термодинамического состояния облака в процессе его подъема в атмосфере. Проведенные измерения температуры облака на основе метода инфракрасной сигнатуры с помощью регистратора тепловизионных изображений реального времени, а также определение траектории и формообразования облака с помощью двухракурсной съемки показали
[15], что изменения температуры облака и его размеров носят периодический пульсационный характер с периодом 6-8 с.
Очевидно, это связано с неравновесным термодинамическим состоянием облака и пульсационным характером процессов конденсации и седиментации. Эти процессы носят высокочастотный характер (6-8 с на цикл по сравнению с 300 с времени горения заряда) и, видимо, не могут повлиять на интегральные показатели зон рассеивания ВК [16].
Для оценки динамики подъема облака ПС разработан ряд математических моделей [14, 16, 17], из них при использовании приближенных оценок целесообразно выбрать методику [16], в которой на инженерном уровне учитываются фазовые переходы компонентов.
В общем случае скорость подъема термика определяется на основе решения уравнения
(М + П)(1 -Р0> = ^ -
а т
где М - масса облака ПС; П - присоединенная масса воздуха; р0 - степень расширения облака ПС; V - скорость подъема облака; ^А - сила Архимеда; - сила сопротивления среды.
Методика1 позволяет рассчитать скорость, массу, температуру облака с учетом эжекции и конденсации для использования в качестве опорных параметров на начальных этапах исследований.
4. Для оценки зон рассеивания ПС в атмосфере большое значение имеет конечная высота подъема облака, которая определяется как высота, на которой значение вертикальной скорости подъема облака мало по сравнению со скоростью ветра. На основе теории распространения турбулентных затопленных струй предложена следующая зависимость [13]:
Н -- 3,75 ™
и10
, 1 о ЯЛТ 1 + 1,3-2 —
и10Та
где Я, Ж, и10, ДТ - соответственно начальный радиус облака, начальная скорость подъема облака и скорость ветра на высоте 10 м, температурный перегрев.
Анализ показывает, что зависимости такого вида могут использоваться в расчетном алгоритме ИУС «Экология».
5. Рассеивание ВК в атмосфере. В настоящее время существует несколько расчетных методик для определения распространения компонентов выбросов в атмосфере. Методика2 предназначена для расчета рассеивания в атмосфере выбросов от промышленных предприятий. Она содержит большой набор расчетных формул и рекомендаций и использовалась при отработке функциональной модели регулятора.
Исходя из перечисленных выше факторов, определяющих динамику подъема и рассеивания облака ПС, может быть принят следующий подход к разработке функциональной модели облака ПС:
1. Исходя из диаграммы давления р(т) и расходной характеристики ВЭЗ выбираются временные участки для расчета параметров отдельных слоев облака ПС.
2. Для каждого расчетного участка, исходя из баланса масс воды и ВК, осажденных в рабочей зоне стенда, определяются массы и концентрации ВК и воды (пара) в составе облака ПС.
3. На основе расчетов и результатов экспериментального измерения определяются начальные параметры облака ПС: температура Т0, вертикальная скорость Ж, радиус Я0.
1 ОНД-86. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий. - Л.: Гидрометеоиздат, 1987. - 68 с.
2 Там же.
4. По начальным параметрам облака определяются конечная высота подъема каждого слоя облака, начальные концентрации ВК в облаке ПС, зоны рассеивания ВК.
5. Полученные высоты и содержание ВК в пределах этих зон сравниваются с ПДК и результатами экологического наземного контроля санитарно-защитной и селитебных зон.
При определении базовых параметров процессов используются следующие зависимости:
Баланс масс по воде, Л120з и HCl определен следующим образом:
ОС ОС ОС
- для технологических стоков тв , тОА, тХВ - по результатам
аналитического контроля;
- для облака ПС в целом
тва = тв - тОС,
тОл = тОл - т°А, тХв = т|в - т££, тПС = тз + т^ -(т°А + т££). Для каждого временного участка т,:
- определяется масса сгоревшего заряда как тПС = = | AP(т)/гкр (x)dт , где ^кр(т) - площадь критического сечения заряда;
A - коэффициент согласования;
- определяются коэффициенты масс Кт. =
тПС,.
тз
- определяется содержание Л1203 и HCl в i-й части облака:
тОл, = тОл • Kт,, тХв,- = тХв • Кт,;
- определяются секундные массовые расходы Л1203 и HCl:
тОл,
тяэл, =-
т,
- а ^ХВ,-тХВ,- =-1
т,
и объемные расходы УОА- , УХВ- .
Высота подъема ,-й части облака находится по формуле
ЛЖЯ
Н, =■
и
в + (Щ/р - Та)
Т и 2
(4)
где Ж, Я, Т0 - начальные вертикальная скорость подъема, радиус и температура облака ПС; и - скорость ветра; Л, В, С - коэффициенты согласования, которые определяются по экспериментальным данным на предварительных СтС; Та - температура атмосферы.
Расчет зон рассеивания и концентраций ВК определяется для каждого слоя облака ПС с учетом рекомендаций [18] по зависимостям.
1. Максимальная концентрация ст для каждого ,-го ВК:
АМЕтпц
Ст "
Н 23ЦАУ'
Коэффициенты, входящие в эту зависимость, определяют следующие факторы: Л - коэффициент, зависящий от температурной стратификации атмосферы; М - масса вредного вещества, выбрасываемого в атмосферу в единицу времени, г/с; ^ - безразмерный коэффициент, учитывающий скорость оседания вредных веществ в атмосферном воздухе; т и п - коэффициенты, учитывающие условия выхода газовоздушной смеси из устья источника выброса; Н - высота источника выброса над уровнем земли, м; п - безразмерный коэффициент, учитывающий влияние рельефа местности, в случае ровной или слабопересеченной местности с перепадом высот, не превышающим 50 м на 1 км, п = 1; ДТ - разность между температурой выбрасываемой газовоздушной смеси Тг и температурой окружающего атмосферного воздуха Тв, °С; У1 - расход газовоздушной смеси, м3/с.
Принимаются в соответствии с рекомендациями3.
3 ОНД-86. Методика расчета концентраций в атмосферном воздухе вредных веществ, содержащихся в выбросах предприятий.
2. Расстояние хт от источника выбросов, на котором приземная концентрация с (мг/м3) при неблагоприятных метеорологических условиях достигает максимального значения ст, определяется по формуле
5 - F
хт =-dH,
т 4
где d - коэффициент согласования.
3. Максимальное значение приземной концентрации вредного вещества сти (мг/м3) при неблагоприятных метеорологических условиях и скорости ветра u
с = rc
ти т 1
где r - безразмерная величина, определяемая в зависимости от отношения и/ит по формулам
r = 0,67(u / ит ) +1,67(u / ит )2 -1,34(и / ит )3
при и/ит < 1,
r = 3(u / ит )
2(u / ит )2 - (и / и„ ) + 2
при и/ит > 1.
Расстояние от источника выброса хти (м), на котором при скорости ветра и и неблагоприятных метеорологических условиях приземная концентрация вредных веществ достигает максимального значения сти (мг/м3), определяется по формуле
Хти = pхm,
где р - безразмерный коэффициент, определяемый в зависимости от отношения и/ит по формулам
p = 3 при и/ит < 0,25; p = 8,43(1 - и / ит )3 +1 при 0,25< и/ит < 1; p = 0,32и / ит + 0,68 при и/ит > 1.
Для идентификации расчетных зависимостей параметров облака и зон рассеивания целесообразно использовать:
- количество вошедших в состав облака ПС, ВК и воды, определяемое расчетом по данным контроля состава технологических стоков в бассейне-отстойнике;
- высоты подъема слоев облака для разных участков горения заряда в зависимости от степени орошения струи; при проведении СтС эти высоты определяются средствами инструментально-приборного контроля на основе разработанной методики, при проектировании - по зависимости вида (4);
- зоны рассеяния и концентрация ВК в приземном слое атмосферы на границе СЗЗ, которые:
а) могут быть определены расчетным путем, исходя из состава и высоты подъема облака. При этом проведение расчета позволяет определить зоны рассеяния ВК, величину и координаты максимальных концентраций в зависимости от скорости ветра при проведении СтС;
б) контролируются в процессе проведения СтС с помощью передвижных лабораторий мониторинга.
Особенность идентификации параметров облака ПС состоит в том, что если состав облака ПС может быть расчитан достаточно точно по результатам аналитического контроля бассейна отстойника, то замеры высоты содержат значительные неопределенности за счет условности понятия «высота», поскольку протяженность отдельных слоев облака на высоте составляет 50-70 м, а замеры содержания ВК в приземном слое атмосферы вообще не могут быть использованы для количественной идентификации необходимых функциональных зависимостей, так как зависят от степени попадания контролирующей лаборатории в зону «следа» облака. Тем не менее результаты экологического контроля можно использовать при проведении факторного анализа и для сравнительной оценки расчета и множества контрольных замеров, полученных при СтС.
Состав и содержание компонентов облака ПС определяются по балансу масс с учетом эффективности осаждения ПС, ВК и воды в рабочей зоне стенда.
Высота подъема облака может определяться по зависимости вида (4), в которой начальные параметры облака ПС Я, Ж и Т0 замеряются при проведении СтС, а скорость ветра и температура воздуха также
контролируются. Неудобство этой формулы состоит в том, что высота подъема Н значительно зависит от скорости ветра, причем при и ^ 0 Н и формула не идентифицируется при и = 0, т.е. для условий штиля.
Кроме того, в рекомендуемом в работе [13] виде данная формула не дает достаточно хорошей сходимости для высот подъема облака при СтС вследствие действия ряда факторов, зависящих от водосодержа-ния облака, в частности:
- изменяются, по сравнению с «сухими» СтС, габариты облака, его плавучесть и коэффициент аэродинамического сопротивления;
- внешний и внутренний тепломассообмен дополняется тепловыделением при конденсации водяного пара и т.д.
Тем не менее три степени свободы формулы (7) по коэффициентам А, В, С позволяют использовать ее для идентификации с учетом переменных АТ, и и Ж.
В частности, для параметров облака, приведенных в табл. 3, эта зависимость может быть записана в виде
Н = 4,5
(1 + и)
1 + 3,3
gR АТ
(1 + и)2 Т
Параметры облака ПС при и = 1 м/с
(5)
Таблица 3
Параметр Вариант СтС
с ВЭЗ без ВЭЗ
Я, м 6 3,25
Ж, м/с 6 10
АТ, К 30 140
Т, К 303 343
Н, м 500 900
Результаты СтС зарядов с массой от 6,0 до 50,0 т подтверждают расходимость этой зависимости от опыта ниже 10 %.
Концентрация и зоны рассеивания ВК. Анализ показывает, что расчетные зависимости для оценки концентраций ВК в приземном слое атмосферы могут использоваться и для оценки экологических рисков при проведении СтС. При этом величина коэффициентов, входящих в эту зависимость, составляет: А = 300, Е = 3,0 для А1203; Е = 2,5
для HCl; m = 1,1, n = 1, n = 1,5^10-2, d = 5,17 для AI2O3; d = 4,67 для HCl. Все другие параметры - M, V, AT, H - определяются из баланса масс и других параметров облака.
Чувствительность к входным параметрам. На основе приведенных выше зависимостей проведено сквозное моделирование параметров струи и облака ПС для исследования влияния изменений расходной характеристики ВЭЗ GB(x, т) на рассеивание ВК. Результаты представлены на графиках рис. 12, 13 и показывают следующее:
- концентрации ВК с, зависят от расходной характеристики GB(x, т), при этом при среднем расходе GB(x, т) = 2,0 т/с достигают минимума. Однако при снижении расхода воды до 1,5 т/с концентрации как Al2O3, так и HCl увеличиваются незначительно и не превышают ПДК;
- зависимость концентрации от расхода воды не является монотонной и должна исследоваться для каждого двигателя на всем диапазоне изменения расходной характеристики;
- увеличение степени орошения струи ПС сдвигает координату максимума концентраций cmi в сторону рабочей зоны стенда и наоборот, что должно учитываться при выборе расходной характеристики и назначении границ СЗЗ;
- координата максимума концентрации xm для HCl расположена дальше от стенда, чем для Al2O3, что необходимо учитывать при определении границ СЗЗ;
- при увеличении среднего расхода воды масса выбросов ВК в атмосферу снижается (см. ниже).
Рис. 12. Максимальная концентрация HCl Рис. 13. Максимальная концентрация на ВНУ и ННУ в зависимости от среднего Al2O3 на ВНУ и ННУ в зависимости расхода воды Gв(x, т) от среднего расхода воды G^x, т)
Суммарные выбросы компонентов ПС в атмосферу
Ов(х, т), т/с 0 0,2 0,4 1 2
тПС, т 11,58 10,59 9,59 6,60 3,95
тол, т 9,37 8,62 7,86 5,61 3,23
тхВ, т 2,22 1,97 1,73 0,99 0,71
Адаптируемость. При проектировании ИУС все контрольные параметры облака ПС определяются расчетным путем и затем первично идентифицируются по результатам установочных СтС. По мере проведения СтС и набора статистических данных проводится уточнение всех контрольных параметров и определение величины их рассеяния, что позволяет снижать экологические риски за счет корректировки расходной характеристики по уточненным данным.
Выходные параметры облака ПС обусловливают величину экологического риска Я0, определяемого по каждому /-му ВК как Я = Вер (е. > [е, ]), где [е,] - допустимая величина концентрации /-го ВК.
Эти оценки могут быть проведены:
- по результатам численного моделирования параметров струи и облака ПС;
- по результатам экологического контроля при проведении СтС.
В обоих случаях концентрация с рассматривается как случайная
величина, распределенная по некоторому заданному закону.
При допущении о нормальности статистического распределения с величина риска может быть определена по зависимости
Я = Ф
V. -[е,] ^
где Ф - табулированная функция Лапласа; те ., а е . - математическое
ожидание и среднеквадратическое отклонение замеренных концентраций с .
При расчетной оценке величину те . (Ов (х, т)) принимают как точечный результат расчета, а а е. - по результатам анализа ошибок расчета [18] и по статистическим данным, полученным при СтС двигателей-аналогов.
В табл. 4 приведены результаты расчетных оценок экологического риска по HCl в зависимости от степени орошения струи.
Таблица 4
Влияние расхода воды на экологический риск по HCl
Расход воды, т/с 0,0 0,2 0,4 1 2
Величина риска ВНУ 0,50 0,84 0,99 0,16 0,023
ННУ 10-3 3 • 10-3 3 10-3 10-4 10-6
Из табл. 4 видно, что снижение суммарного расхода воды ниже 1,0 т/с на высоконапорном участке значительно увеличивает экологический риск. Вместе с тем для низконапорного участка имеются значительные резервы для такого снижения.
Заключение
На основе инженерных физико-математических моделей предложена методология описания термогазодинамических процессов, происходящих при огневой утилизации зарядов крупногабаритных РДТТ на открытом стенде с использованием водной экологической защиты.
Предложенные модели позволяют:
- контролировать основные параметры процессов, определяющих формирование и рассеивание облака ПС в атмосфере и концентрацию ВК в приземном слое атмосферы;
- выбирать конструкцию и расходную характеристику ВЭЗ по воде с учетом внутрибаллистических параметров заряда, особенностей его горения при проведении СтС РДТТ без выходного блока;
- проводить идентификацию и последовательное уточнение моделей по результатам СтС и оценивать экологические риски с учетом физических неопределенностей процессов.
Библиографический список
1. Попова М.В., Литвинов А.В., Козлов С.Н. Информационная структура контроля и обеспечения экологической безопасности при огневой утилизации крупногабаритных зарядов РДТТ на открытом стенде // Вестник Пермского национального исследовательского политехнического университета. Аэрокосмическая техника. - 2013. -№ 35. - С. 5-20.
2. Шевырин А.Ю., Лебедев А.С. Оценка параметров струи продуктов сгорания твердого топлива крупногабаритных ракетных двигателей при их ликвидации методом открытого сжигания // Высокоэнергетические материалы: демилитаризация, антитерроризм и гражданское применение: материалы IV Междунар. конф. «HEMS - 2008» (Белокуриха, 3-5 сентября 2008 г.). - Бийск, 2008. - С. 189.
3. Абрамович Г.Н. Прикладная газовая динамика. - М.: Гостехиз-дат, 1953. - 736 с.
4. Гинзбург И.П. Прикладная гидрогазодинамика - Л.: Изд-во Ленинград. гос. ун-та, 1958. - 338 с.
5. Особенности процессов при запуске энергетических установок без выходных блоков / В.И. Марьяш, В.С. Аверин, А.А. Назаров [и др.] // Материалы всерос. науч.-практ. конф. ФГУП ФНПЦ «Алтай» (Бийск, 26-28 сентября 2001 г.). - Бийск, 2003. - С. 28-33.
6. Гребенкин В.И., Павлов А.Д., Шмачков Е.А. Итоги экспериментальных бессопловых сжиганий РДТТ и перспективы исследований в обеспечение отработки технологии их утилизации // Материалы все-рос. науч.-практ. конф. ФГУП ФНПЦ «Алтай» (Бийск, 26-28 сентября 2001 г.). - Бийск, 2003. - С. 21-24.
7. Определение температурных параметров факела продуктов сгорания и его воздействия на элементы конструкции двигателя и стенда при бессопловом сжигании / В. И. Гребенкин, А. Д. Павлов, Е.А. Шмачков [и др.] // Материалы всерос. науч.-практ. конф. ФГУП ФНПЦ «Алтай» (Бийск, 26-28 сентября 2001 г.). - Бийск, 2003. -С. 57-60.
8. Нащокин В.В. Техническая термодинамика и теплопередача: учеб. пособие для неэнергетич. специальностей вузов. - М.: Высшая школа, 1975. - 496 с.
9. Мамонтов М.А. Вопросы термодинамики тела переменной массы. - М.: Оборонгиз, 1962. - 102 с.
10. Жигульский А.И., Крюков В.Г., Литвинов А.В. Интенсификация процессов абсорбции при водном орошении струи продуктов сгорания РДТТ // Высокоэнергетические материалы: демилитаризация, антитерроризм и гражданское применение: материалы IV Междунар. конф. «HEMS - 2008» (Белокуриха, 3-5 сентября 2008 г.). - Бийск, 2008. - С. 193.
11. Диментова А. А. Таблицы газодинамических функций: справ. пособие. - М.: Машиностроение, 1966. - 138 с.
12. Газодинамические параметры струи при внешнем орошении ее водой / С.Н. Козлов, М.В. Попова, А.В. Литвинов, А.В. Скворцов // Материалы всерос. науч.-техн. конф. «Информационные технологии в науке, экономике и образовании» (Бийск, 22-23 сентября 2011 г.). -Бийск, 2011. - С. 36-41.
13. Математическое моделирование горения и взрыва высокоэнергетических систем: монография / под ред. И.М. Васенина. - Томск: Изд-во Том. ун-та, 2006. - 322 с.
14. Полуэмпирическая модель процессов, происходящих при образовании и распространении высокодисперсных аэрозолей окислов алюминия при открытом сжигании крупногабаритных РДТТ / С.Э. Пащенко, А.Е. Осоченко, В.Е. Зарко [и др.] // Материалы всерос. науч.-практ. конф. ФГУП ФНПЦ «Алтай» (Бийск, 26-28 сентября 2001 г.). - Бийск, 2003. - С. 83-97.
15. Динамика подъема облака продуктов сгорания при огневой утилизации зарядов твердотопливных энергоустановок на открытом стенде. Процессы и их идентификация / В.П. Лушев, А.В. Литвинов, Н.Ю. Демидов, С.Н. Козлов, В.В. Рейно // Оптика атмосферы и океана. -Новосибирск, 2014. - Т. 27, № 10. - С. 917-924.
16. Берлянд М.Е. Прогноз и регулирование загрязнения атмосферы. - Л.: Гидрометеоиздат, 1985. - 272 с.
17. Воздействие на окружающую среду кратковременных выбросов большой мощности: учеб. пособие / В.М. Суслонов, Н.Г. Максимович, В.Н. Иванов, В. А. Шкляев; Перм. гос. ун-т. - Пермь, 2005. - 126 с.
18. Ситуационный анализ экологических рисков при сжигании твердотопливных зарядов энергетических установок на открытом стенде / М.В. Попова, А.В. Литвинов, С.Н. Козлов, В.П. Лушев // Известия Том. политехн. ун-та. - 2013. - Т. 322, № 1. - С. 157-161.
References
1. Popova M.V., Litvinov A.V., Kozlov S.N. Informatsionnaya struk-tura kontrolya i obespecheniya ekologicheskoy bezopasnosti pri ognevoy utilizatsii krupnogabaritnykh zaryadov RDTT na otkrytom stende [Information structure of control and provision of ecological safety at the firing utilization of large-size SRM charges at an open bench]. Vestnik Permskogo
natsionalnogo issledovatelskogo politekhnicheskogo universiteta. Aerokos-micheskaya tekhnika, 2013, no. 35, pp. 5-20.
2. Shevyrin A.Yu., Lebedev A.S. Otsenka parametrov strui produktov sgoraniya tverdogo topliva krupnogabaritnykh raketnykh dvigateley pri ikh likvidatsii metodom otkrytogo szhiganiya [Evaluation of jet parameters of solid propellant combustion products of large-size SRM charges on their utilization by open firing]. Materialy IV Mezhdunarodnoy konferentsii "HEMS - 2008" "Vysokoenergeticheskie materialy: demilitarizatsiya, anti-terrorizm i grazhdanskoe primenenie". Biysk, 2008, p. 189.
3. Abramovich G.N. Prikladnaya gazovaya dinamika [Applied gas dynamics]. Moscow: Gostekhizdat, 1953. 736 p.
4. Ginzburg I.P. Prikladnaya gidrogazodinamika [Applied hydro-gasdynamics]. Leningradskiy gosudarstvennyy universitet, 1958. 338 p.
5. Maryash V.I., Averin V.S., Nazarov A.A. [et al.]. Osobennosti prot-sessov pri zapuske energeticheskikh ustanovok bez vykhodnykh blokov [Features of processes on firing propulsions without nozzles]. Materialy vserossiyskoy nauchno-prakticheskoy konferentsii FGUP FNPTs "Altay". Biysk, 2003, p. 28-33.
6. Grebenkin V.I., Pavlov A.D., Shmachkov E.A. Itogi eksperimental-nykh bessoplovykh szhiganiy RDTT i perspektivy issledovaniy v obe-spechenii otrabotki tekhnologii ikh utilizatsii [Results of experimental firings of SRM without nozzles and research prospects of finishing their utilization technology]. Materialy vserossiyskoy nauchno-prakticheskoy konferentsii FGUP FNPTs "Altay". Biysk, 2003, p. 21-24.
7. Grebenkin V.I, Pavlov A.D., Shmachkov E.A. [et al.] Opredelenie temperaturnykh parametrov fakela produktov sgoraniya i ego vozdeystviya na elementy konstruktsii dvigatelya i stenda pri bessoplovom szhiganii [Determination of temperature parameters of combustion products plume and its influence on engine construction elements and test bench at nozzle-free firing]. Materialy Vserossiyskoy nauchno-prakticheskoy konferentsii FGUP FNPTs "Altay". Biysk, 2003, p. 57-60.
8. Nashchokin V.V. Tekhnicheskaya termodinamika i teploperedacha [Technical thermodynamics and heat transfer]. Moscow: Vysshaya shkola, 1975. 496 p.
9. Mamontov M.A. Voprosy termodinamiki tela peremennoy massy [Problems of thermodynamics of varying mass body]. Moscow: Oborongiz, 1962. 102 p.
10. Zhigulskiy A.I., Kryukov V.G., Litvinov A.V. Intensifikatsiya protsessov absorbtsii pri vodnom oroshenii strui produktov sgoraniya RDTT [Intensification of absorption processes on water irrigation of SRM combustion products jet]. Materialy IV Mezhdunarodnoy konferentsii "HEMS -2008" "Vysokoenergeticheskie materialy: demilitarizatsiya, antiterrorizm i grazhdanskoe primenenie". Biysk, 2008, p. 193.
11. Dimentova A.A. Tablitsy gazodinamicheskikh funktsiy [Tables of gasdynamic functions]. Moscow: Mashinostroenie, 1966. 138 p.
12. Kozlov S.N., Popova M.V., Litvinov A.V., Skvortsov A.V. Gazo-dinamicheskie parametry strui pri vneshnem oroshenii ee vodoy [Gas-dynamic jet parameters at its external irrigation]. Materialy Vserossiyskoy nauchno-tekhnicheskoy konferentsii "Informatsionnye tekhnologii v nauke, ekonomike i obrazovanie". Biysk, 2011, p. 36-41.
13. Matematicheskoe modelirovanie goreniya i vzryva vysokoener-geticheskikh sistem [Mathematical modeling of combustion and explosion of high energy systems]. Ed. I.M. Vasenin. Tomskiy universitet, 2006. 322 p.
14. Pashchenko S.E., Osochenko A.E., Zarko V.E. [et al.]. Poluem-piricheskaya model protsessov, proiskhodyashchikh pri obrazovanii i raspros-tranenii vysokodispersnykh aerozoley okislov alyuminiya pri otkrytom szhi-ganii krupnogabaritnykh RDTT [Semiempirical model of processes occurring on formation and propagation of highly dispersed aerosol of aluminum oxides at open firing of large-size SRM]. Materialy vserossiyskoy nauchno-prakticheskoy konferentsii FGUP FNPTs "Altay". Biysk, 2003, p. 83-97.
15. Lushev V.P., Litvinov A.V., Demidov N.Yu., Kozlov S.N., Reyno V.V. Dinamika podema oblaka produktov sgoraniya pri ognevoy uti-lizatsii zaryadov tverdotoplivnykh energoustanovok na otkrytom stende. Protsessy i ikh identifikatsiya [Dynamics of lifting combustion products cloud when firing utilization of propulsion solid propellant charges at an open test bench]. Optika atmosfery i okeana. Novosibirsk, 2014, vol. 27, no. 10, pp. 917-924.
16. Berlyand M.E. Prognoz i regulirovanie zagryazneniya atmosfery [Prediction and control of atmosphere pollution]. Leningrad: Gidrometeoiz-dat, 1985. 272 p.
17. Suslonov V.M., Maksimovich N.G., Ivanov V.N., Shklyaev V.A. Vozdeystvie na okruzhayushchuyu sredu kratkovremennykh vybrosov bol-shoy moshchnosti [Influence of short-term exhausts of high power on the environment]. Permskiy gosudarstvennyy universitet, 2005. 126 p.
18. Popova M.V., Litvinov A.V., Kozlov S.N., Lushev V.P. Situat-sionnyy analiz ekologicheskikh riskov pri szhiganii tverdotoplivnykh zary-adov energeticheskikh ustanovok na otkrytom stende [Situational analysis of ecological risks on firing of solid propellant propulsion plants at an open test bench]. Izvestiya Tomskogo politekhnicheskogo universiteta, 2013, vol. 322, no. 1, pp. 157-161.
Об авторах
Попова Мария Владиславовна (Бийск, Россия) - кандидат физико-математических наук, ведущий инженер ОАО ФНПЦ «Алтай» (659322, г. Бийск, ул. Социалистическая, д. 1, e-mail: popova.maria.v@gmail.com).
Литвинов Андрей Владимирович (Бийск, Россия) - кандидат технических наук, заместитель генерального директора - директор -главный конструктор по НИОКР ОАО ФНПЦ «Алтай» (659322,
г. Бийск, ул. Социалистическая, д. 1, e-mail: post@frpc.secna.ru).
Козлов Станислав Николаевич (Бийск, Россия) - кандидат технических наук, профессор кафедры «Ракетные двигатели и высокоэнергетические устройства автоматических систем» Бийского технологического института (филиала) ФГОУ ВПО «Алтайский государственный технический университет им. И.И. Ползунова» (659305, г. Бийск, ул. Трофимова,
д. 27, e-mail: magistrus@city.biisk.ru).
About the authors
Mariya V. Popova (Biysk, Russian Federation) - Ph. D. in Physical and Mathematical Sciences, Lead Engineer, OJSC "Altay" (1, Sotsya-listycheskaya st., Biysk, 659322, Russian Federation, e-mail: pop-ova.maria.v@gmail.com).
Andrey V. Litvinov (Biysk, Russian Federation) - Ph. D. in Technical Sciences, Deputy Director General - Director, OJSC "Altay" (1, Sot-syalistycheskaya st., Biysk, 659322, Russian Federation, e-mail: post@frpc.secna.ru).
Stanislav N. Kozlov (Biysk, Russian Federation) - Ph. D. in Technical Sciences, Professor, Department of Rocket Engines and High-energy Devices of Automatic Systems, Biysk Tecnological Institute (27, Trofimova st., Biysk, 659305, Russian Federation, e-mail: magistrus@city.biisk.ru).
Получено 2.03.2015
