Приборы и техника эксперимента, 2019, № 1, стр. 94-98
УСТАНОВКА ДЛЯ БЫСТРОГО И БЕЗОПАСНОГО ПРИГОТОВЛЕНИЯ РАБОЧИХ СМЕСЕЙ В РЕАКТОРАХ ИМПУЛЬСНЫХ ХИМИЧЕСКИХ HF(DF)-ЛАЗЕРОВ
В. Я. Агроскин, Б. Г. Бравый, Г. К. Васильев, В. И. Гурьев, С. А. Каштанов, Е. Ф. Макаров *, С. А. Сотниченко, Ю. А. Чернышев **
Институт проблем химической физики РАН
142432 Московской обл., Черноголовка, ул. Академика Семенова, 1, Россия
* E-mail: makarov@icp.ac.ru
** E-mail: chern@icp.ac.ru
Поступила в редакцию 21.03.2018
После доработки 21.03.2018
Принята к публикации 15.04.2018
Аннотация
Описана система для приготовления рабочих смесей в реакторах импульсных химических лазеров на основе реакций F2 + H2(D2). Разработана конструкция компактного капиллярного смесителя с достигнутой в экспериментах пропускной способностью до 2 л · атм/с. Предложена схема напуска газов с использованием расходных шайб с отверстиями, работающими в критическом режиме, которая обеспечивает постоянство расхода смешиваемых газов в течение всего напуска. Теоретически проанализирована и экспериментально осуществлена описанная технология приготовления рабочих смесей импульсных химических HF(DF)-лазеров.
Многолетний опыт работы с мощными импульсными химическими лазерами на реакции F2 + H2(D2) показал актуальность проблемы быстрого и безопасного приготовления рабочих смесей в активных объемах лазеров. Фтор – один из самых реакционно-способных химических элементов. Реакция F2 + H2(D2) протекает по разветвленно-цепному механизму, вследствие чего приготовление смесей при одновременном напуске реагентов в реактор может сопровождаться вхождением их в область самовоспламенения. Как показано в работе [1], смесь стабилизируется образующимся при напуске фтористым водородом, однако, с другой стороны, фтористый водород может негативно влиять на характеристики генерации. Поэтому добавочно возникает проблема приготовления смесей с минимальным количеством фтористого водорода.
Эти и другие возникающие при смешении проблемы, обусловленные разветвленно-цепным механизмом реакции фторирования водорода, детально обсуждаются в работе [1]. Существует еще одна причина воспламенения приготавливаемой смеси, обусловленная возможным наличием застойных зон в напускном тракте. Застойные зоны являются очагами воспламенения и стабилизаторами горения.
Данная работа посвящена описанию разработанных нами эффективного компактного смесителя и технологии быстрого и безопасного приготовления фтороводородных смесей, обеспечивающей постоянную скорость напуска реагентов в реактор. Предложенные решения подходят в случае приготовления любых газовых смесей.
КОНСТРУКЦИЯ СМЕСИТЕЛЯ
Для безопасного приготовления смесей использовали технологию, в основе которой лежат следующие главные моменты: 1) быстрота смешения реагентов обеспечивается за счет разделения их на множество перемежающихся струй; 2) устраняется возможность образования застойных зон; 3) скорость движения смеси по напускному тракту всюду превышает нормальную скорость горения. Подходящим оказался смеситель на основе структуры из чередующихся отверстий для окислителя и горючего. Характерный размер структуры (расстояние между соседними отверстиями) выбирался достаточно малым, чтобы исключить возможность возникновения застойных зон путем образования отрывного циркуляционного течения, способствующего стабилизации пламени за счет переноса тепла и активных центров к свежей смеси обратными токами.
Изготовленный из нержавеющей стали смеситель содержал 48 капилляров (из нержавеющей стали) с внешним диаметром 0.1 cм и внутренним 0.04 cм, плотно упакованных в цилиндрический канал с внутренним диаметром 0.8 cм. Конструкция смесителя схематически представлена на рис. 1.
Через капилляры 3 (на рис. 1 показаны 2 капилляра из 48), используя ввод 5, подавали “предсмесь”, содержащую H2 (или D2), а через промежутки между капиллярами, используя ввод 6, – “предсмесь”, содержащую F2. Обычно использовали “предсмеси” с составом 0.1H2 (или D2) + 0.9He и 0.3F2 + 0.06O2 + 0.64He, дающие при равном расходе требуемую рабочую смесь 0.15F2 + 0.05H2 + + 0.03O2. Кислород добавляется для стабилизации смесей.
Концы капилляров впаивались с помощью свинцово-оловянного припоя в круглую латунную пластину 4 с 48 отверстиями. Пластина вакуумно плотно устанавливалась между корпусом 1 и фланцем 2. Эта пластина разделяла полости, через которые подавались “предсмеси” в смеситель.
Конструкция смесителя обеспечивала в экспериментах суммарную пропускную способность до 2 л · атм/с. Для обычных химических HF(DF)-лазеров с энергией в импульсе 100–150 Дж и рабочим объемом 5–10 л, используемых в наших экспериментах, время готовности лазера к работе при этом составляло 3–5 с. Наибольший габаритный размер смесителя не превышал 10 см.
ТЕОРЕТИЧЕСКАЯ МОДЕЛЬ РАБОТЫ КАПИЛЛЯРНОГО СМЕСИТЕЛЯ
Рассмотрим случай истечения газа через капилляр в среду с возрастающим по мере напуска газа давлением. При движении сжимаемого газа в капилляре с трением связь между безразмерными скоростями потока во входном 1 и выходном 2 сечениях задается соотношением [2, с. 249]:
(1)
$1{/\lambda }_{1}^{2} - 1{/\lambda }_{2}^{2} - \ln ({\lambda }_{2}^{2}{/\lambda }_{1}^{2}) = {\chi },$Для капилляра с круглым сечением при сравнительно небольших числах Рейнольдса (Re < 2300), когда течение носит ламинарный характер, коэффициент ζ = 64/Re. Число Рейнольдса записывается как Re = сd/ν, где с – средняя расходная скорость движения газа, ν – кинематическая вязкость газа. Перепишем число Рейнольдса в следующем виде: Re = cdρf/(νρf), где ρ – плотность газа, f = πd2/4 – площадь сечения капилляра. Учитывая, что массовый расход газа через капилляр m = cρf = M/N (M – массовый расход газа через N капилляров), νρ = η (η – динамическая вязкость газа), для числа Рейнольдса получим выражение Re = 4M/πηdN. Тогда
Массовые расходы газа в сечениях 1 и 2 с учетом всех капилляров записываются M1 = c1ρ1F и M2 = c2ρ2F, где F = fN – сумма площадей отверстий всех капилляров. Представим выражения для массовых расходов в виде $({{c}_{1}}{\text{/}}с _{1}^{*})с _{1}^{*}({{\rho }_{1}}{\text{/}}{{\rho }_{{01}}}){{\rho }_{{01}}}F$ и $({{c}_{2}}{\text{/}}с _{2}^{*})с _{2}^{*}({{\rho }_{2}}{\text{/}}{{\rho }_{{02}}}){{\rho }_{{02}}}F$. Здесь и далее индекс “0” будем относить к параметрам заторможенного газа. По определению ${{c}_{1}}{\text{/}}с _{1}^{*} = {{\lambda }_{1}}$, ${{c}_{2}}{\text{/}}с _{2}^{*} = {{\lambda }_{2}}$. Из [2, с. 58–59] $\frac{{{{{\rho }}_{1}}}}{{{{{\rho }}_{{01}}}}} = {{\left( {1 - \frac{{k - 1}}{{k + 1}}{\lambda }_{1}^{2}} \right)}^{{\frac{1}{{k - 1}}}}}$, $\frac{{{{{\rho }}_{2}}}}{{{{{\rho }}_{{02}}}}} = {{\left( {1 - \frac{{k - 1}}{{k + 1}}{\lambda }_{2}^{2}} \right)}^{{\frac{1}{{k - 1}}}}}$, $c_{1}^{*}$ = = $\sqrt {\frac{{2k}}{{k\, + \,1}}\frac{{R{{T}_{{01}}}}}{{\mu }}} $, $c_{2}^{*} = \sqrt {\frac{{2k}}{{k + 1}}\frac{{R{{T}_{{02}}}}}{{\mu }}} $.
Из уравнения состояния ρ01 = P01μ/RT01, ρ02 = = P02μ/RT02, где μ – молекулярный вес газа, P01, P02 и T01, T02 – давление и температура торможения в сечениях 1 и 2 соответственно. При течении без энергообмена (изоэнтропийное течение) T01 = = T02 = T0. С учетом сказанного можно получить:
(4)
${{M}_{1}} = {{{\lambda }}_{1}}{{\left( {1 - \frac{{k - 1}}{{k + 1}}{\lambda }_{1}^{2}} \right)}^{{\frac{1}{{k - 1}}}}}\sqrt {\frac{{2k}}{{k + 1}}\frac{{\mu }}{{R{{T}_{0}}}}} {{P}_{{01}}}F,$(5)
${{M}_{2}} = {{\lambda }_{2}}{{\left( {1 - \frac{{k - 1}}{{k + 1}}\lambda _{2}^{2}} \right)}^{{\frac{1}{{k - 1}}}}}\sqrt {\frac{{2k}}{{k + 1}}\frac{{\mu }}{{R{{T}_{0}}}}} {{P}_{{02}}}F.$Статическое давление P2 в сечении 2, выраженное в долях давления полного торможения P02, записывается [2, с. 59] в виде
(6)
$\frac{{{{P}_{2}}}}{{{{P}_{{02}}}}} = {{\left( {1 - \frac{{k - 1}}{{k + 1}}{\lambda }_{2}^{2}} \right)}^{{\frac{k}{{k - 1}}}}}.$С учетом (6) соотношение (5) перепишется как
(7)
${{M}_{2}} = {{{\lambda }}_{2}}{{\left( {1 - \frac{{k - 1}}{{k + 1}}{\lambda }_{2}^{2}} \right)}^{{ - 1}}}\sqrt {\frac{{2k}}{{k + 1}}\frac{{\mu }}{{R{{T}_{0}}}}} \,{{P}_{2}}F.$При истечении через капилляр скорость газа не может быть больше местной скорости звука, следовательно, расширение в капилляре осуществляется до давлений, не меньших критического давления $Р _{2}^{*}$. При Р2 = $Р _{2}^{*}$ скорость λ2 = 1. Если давление газа в реакторе P3 меньше критического, то расход газа через капилляр определяется величиной $Р _{2}^{*}$ и не зависит от давления в реакторе. Если давление в реакторе становится больше критического, скорость истечения становится дозвуковой и давление P2 на срезе капилляра сравнивается с давлением в реакторе P3.
Подставив в формулу (1) выражение для χ(3), можем получить
(8)
$M = \frac{{32{\pi \eta }NL}}{{(1{/\lambda }_{1}^{2} - 1{/\lambda }_{2}^{2} - \ln ({\lambda }_{2}^{2}{/\lambda }_{1}^{2}))(k + 1)}}.$Массовые расходы газа M1 и M2 в сечениях 1 и 2 в силу неразрывности течения равны между собой. При выполнении равенства M1 = M2 = M система уравнений (4), (7), (8) дает связь между пятью параметрами – P01, Р2, λ1, λ2, M, и если задать два из них, то остальные три однозначно определяются этой системой. Удобно вместо массового расхода газа M [г/c] использовать объемный расход газа Q [л · атм/с]. Связь между ними дается соотношением Q = (M/μ)V0, где V0 – объем идеального газа, равный при стандартных условиях 22.4 л · атм/моль. На рис. 2 приведены зависимости объемного расхода “предсмеси”, содержащей водород, от давления в реакторе P3 при различных величинах давления на входе капилляров. Использованные в расчетах значения параметров:
По расчетам при Q = 3.73 л · атм/с (соответствует значению P01 = 2 атм на рис. 2) число Рейнольдса Re = 2100 и при дальнейшем увеличении давления P01 газовый поток перестает быть ламинарным. Горизонтальные участки кривых на рис. 2 отвечают ситуации, когда давление в реакторе P3 меньше критического давления $Р _{2}^{*}$. При P3 > $Р _{2}^{*}$ наблюдается падение расхода газа через капилляры с увеличением давления в реакторе для всех значений P01, причем относительная величина падения растет с уменьшением давления на входе капилляров.
ТЕХНОЛОГИЯ ПРИГОТОВЛЕНИЯ РАБОЧЕЙ СМЕСИ
На практике желательно иметь постоянную скорость напуска газа в реактор вне зависимости от возрастающего по мере напуска давления в нем. Для решения этой задачи в напускной тракт перед капиллярами помещали шайбу с отверстием, работающим в критическом режиме. Шайбу изготавливали из алюминиевой фольги толщиной 0.5 мм. Отверстие в фольге пробивали стальной конической иглой. Размер отверстия контролировали с помощью микроскопа отсчетного типа с увеличением 24×.
Уравнение расхода газа через критическое сечение записывается в виде M* = ρ*c*F* = (ρ*/ρ0)c*F*ρ0.
Учитывая, что $\frac{{{\rho *}}}{{{{{\rho }}_{0}}}} = {{\left( {\frac{2}{{k + 1}}} \right)}^{{\frac{1}{{k - 1}}}}}$ [2, с. 62], ρ0 = = P0μ/RT0, ${{c}^{*}} = \sqrt {\frac{{2k}}{{k + 1}}\frac{{R{{T}_{0}}}}{{\mu }}} $, получаем
(9)
${{M}^{*}} = {{\left( {\frac{2}{{k + 1}}} \right)}^{{\frac{{k + 1}}{{2(k - 1)}}}}}\sqrt {\frac{{k{\mu }}}{{R{{T}_{0}}}}} {{P}_{0}}{{F}^{*}},$Отметим, что проведенный расчет касался изоэнтропийного течения. Реально картина течения газа сложнее. Изоэнтропийным можно считать течение до и после капилляра, течение газа в капилляре происходит с энергообменом в силу малости поперечного размера капилляра. В предельном случае бесконечно быстрого энергообмена течение в капилляре будет изотермическим с температурой газа, равной температуре стенки. Можно показать, что в этом случае уравнение (1) преобразуется в уравнение
(10)
$((k + 1){\text{/}}2k)(1{/\lambda }_{1}^{2} - 1{/\lambda }_{2}^{2}) - \ln ({\lambda }_{2}^{2}{/\lambda }_{1}^{2}) = {\chi }{\text{.}}$.Отличие в результатах расчетов при замене уравнения (1) на уравнение (10) не превышает нескольких процентов.
В описываемом смесителе, как указывалось, 48 капилляров с внешним диаметром 0.1 см упакованы в цилиндр с внутренним диаметром 0.8 см. Между капиллярами при их упаковке образуются каналы, отличающиеся, вообще говоря, как поперечными размерами, так и формой вследствие неидеальности упаковки. В каналы с наибольшими размерами вставляли проволоку подходящего диаметра. Общее число каналов в конечном итоге составило 64. Через эти каналы подавали “предсмесь” с фтором. Из-за неодинаковости каналов строгий расчет для “предсмеси” с фтором затруднителен. Однако с учетом вышеизложенного можно обойтись вообще без расчета: при заданном Q л · атм/с нужно просто подобрать начальное давление P0 и площадь отверстия F* такими, чтобы они обеспечивали нужный расход и его постоянство в течение напуска, что нетрудно.
Объемные расходы Q обеих “предсмесей” выбирали одинаковыми. В этом случае суммарный расход газа через смеситель равен 2Q и средняя скорость в соединительной трубке с внутренним диаметром D = 0.8 см составляет с = (2Q/P3)/(πD2/4). Скорость потока минимальна в конце напуска, когда P3 = 1 атм, и при Q = 1 л · атм/с, например, она составит 40 м/с, что больше скорости распространения пламени [3].
В напускной тракт созданной установки помимо смесителя и шайб с отверстиями, задающими критический расход “предсмесей”, входили быстродействующие электроклапаны для включения и отключения газовых потоков и резервуары с “предсмесями”. Давление газа в резервуаре задает начальное давлению P0 перед соответствующим отверстием. Объем резервуара выбирался таким, чтобы в ходе напуска величина P0 падала несильно.
Обычно падение давления не превышало 10–15%. Соответственно скорости напуска смеси в начале и в конце отличались на 10–15%, что допустимо. При необходимости эту разницу можно уменьшить, например, увеличив объем расходных резервуаров или используя устройства, поддерживающие в них постоянное давление. После напуска давление в резервуарах восстанавливалось до требуемой величины за счет “предсмесей”, заранее приготовленных и хранимых в стальных баллонах объемом 10 л при давлении 15–25 атм. Эксплуатационные параметры созданной установки сохраняли стабильность в течение длительного времени работы, обеспечивая воспроизводимость характеристик лазера от опыта к опыту.
ЗАКЛЮЧЕНИЕ
Подытожим сказанное. Разработана конструкция компактного смесителя, обеспечивающего эффективное смешение газов, с суммарной производительностью до 2 л · атм/с. Смеситель обеспечивает готовность лазера к работе в течение нескольких секунд. При разработке смесителя особое внимание уделено проблеме безопасного смешения высокореакционных газов. Предложена схема напуска газов с использованием расходных шайб с отверстиями, работающими в критическом режиме, которая обеспечивает постоянство расхода смешиваемых газов в течение всего напуска.
Работа выполнена при финансовой поддержке госзадания 01201361840.
Список литературы
Васильев Г.К., Макаров Е.Ф., Чернышев Ю.А. // Физика горения и взрыва. 1983. Т. 39. № 3. С. 9.
Дойч М.Е., Зарянкин А.Е. Гидрогазодинамика. М.: Энергоиздат, 1984.
Бурцев В.В., Великанов С.Д., Фролов Ю.Н. // Квантовая электроника. 1995. Т. 22. № 2. С. 123.
Дополнительные материалы отсутствуют.
Инструменты
Приборы и техника эксперимента