УДК 623.52
ОЦЕНКА СТОЙКОСТИ ТЕРМОЗАЩИТНЫХ ПОКРЫТИЙ ПРИ ИМПУЛЬСНОМ ВОЗДЕЙСТВИИ ВЫСОКОЭРОЗИОННЫХ ПРОДУКТОВ ГОРЕНИЯ ТОПЛИВНЫХ ЭЛЕМЕНТОВ
Г.В. Лепеш1, Е.Н. Моисеев2, Д.Ю. Латышев3, М.В. Басова4
1 Санкт-Петербургский государственный экономический университет (СПбГЭУ), Россия, 191023, Санкт-Петербург, наб. канала Грибоедова, д. 30-32, литер А.
1-4АО «Центральный научно-исследовательский институт материалов» им. Д.И. Менделеева,
Россия, 191014, Санкт-Петербург, Парадная ул., д. 8.
Статья посвящена программно-методическому и экспериментальному обеспечению исследований стойкости термозащитных покрытий в жестких условиях термосилового нагружения импульсным потоком топливных газов, при горении топливных элементов в условиях высоких давлений и температур.
Ключевые слова: термозащитные покрытия, разгарно-эрозионное действие, газодинамическая импульсная сопловая установка,
ASSESSMENT OF THE RESISTANCE OF THERMAL PROTECTIVE COATINGS UNDER THE PULSED EFFECT OF HIGHLY CORROSIVE COMBUSTION PRODUCTS OF FUEL CELLS
G.V. Lepesh, E.N. Moiseev, D.Y. Latyshev, M.V. Basova
St. Petersburg State University of Economics (SPbSEU), Russia, 191023, St. Petersburg, nab. Griboedov Canal,. 30-32, letter A.
Central Research Institute of Materials JSC named after D.I. Mendeleev, Russia, 191014, St. Petersburg, 8 Paradnaya str.
The article is devoted to the software, methodological and experimental support of research on the resistance of thermal protective coatings under harsh conditions of thermal force loading by a pulsed flow of fuel gases, during combustion of fuel cells under high pressures and temperatures.
Keywords: thermoprotective coatings, heat-erosion effect, gas-dynamic pulse nozzle installation.
Введение
Термозащитные покрытия (ТЗП) применяются для защиты поверхностей высоконагру-женных деталей от повреждений высокими тепловыми и газодинамическими потоками продуктов горения (ПГ) топливных элементов (ТЭ), вызывающими такие явления как: эрозия, деформация или плавление, являющиеся результатом их разгарно-эрозионного действия (РЭД). ТЗП должны обладать высокими температурами плавления и предотвращать передачу тепла на подложку. Они могут состоять из специальных теплостойких материалов, таких как керамика, керамические композиты или металлические сплавы, или быть металлическими [1,2], способными выдерживать силовые и термосиловые нагрузки. Современные ТЗП (Сг, Ц ВК-8, М, NbN, Ta, TaW,
ЭП-131, ТВ2, Мо, Ци др.) в целом удовлетворяют требованиям, обеспечивающим выполнение функции защиты по эрозионной стойкости в условиях газодинамического действия ПГ при высоких давлениях и скоростях потока. Они обладают высокой твердостью и механической прочностью. Им присущи также свойства высокой температуры плавления (выше 1700 °С) и химической нейтральности (инертности) к компонентам ПГ. Однако, для обеспечения стойкости таких покрытий в условиях высоких давлений и температур их коэффициенты линейного термического расширения должны быть близки к соответствующим значениям для защищаемой, как правило стальной подложки, что является проблемой [1,2].
EDN FAKNQU
1 Лепеш Григорий Васильевич - доктор технических наук, профессор, заведующий кафедрой Безопасность населения и территорий от ЧС, СПбГЭУ, тел.: +7 (921) 751-28-29, e-mail: [email protected], ORCID: 0000-0002-4160-3292, Scopus: 57215412255;
2Моисеев Евгений Николаевич - заместитель начальника лаборатории наноматериалов и карбидных композитов, тел.: +7 904 640-75-92, e-mail: [email protected];
ъЛатышев Денис Юрьевич - начальник отдела «Живучести», тел.: +7 921 318-57-38, e-mail: [email protected];
4Басова Мария Владимировна -начальник сектора НИОКР отдела «Живучести», тел.: +7 911 959-08-10, e-mail: basovam@ cniim.spb.ru.
Проблемой является и технология нанесения покрытия, которая должна обеспечивать многослойную конструкцию, состоящую, как правило, из слоев разной толщины, обладающих различными механическими и термическими свойствами. При этом должна быть обеспечена адгезионная прочность покрытия с основой на уровне механической прочности основного металла.
В случае воздействия коротких импульсных нагружений большой интенсивности [1,2], измеряемых сотыми долями секунды, в качестве защитных покрытий могут применяться металлические составы и композиции. Несмотря на их высокую теплопроводность в относительно короткое время с их применением, можно обеспечить тепловой барьер, необходимый, чтобы на стальной поверхности подложки не были получены критические температуры, приводящие к разрушению металла. Как правило металлические ТЗП формируют непосредственно на поверхности уже готового изделия, применяя при этом технологии, обеспечивающие необходимую адгезию в тонком (несколько микрон) поверхностном слое и сохраняющие неизменными свойства материала подложки. Многие из этих технологий были подвергнуты испытаниям в АО «ЦНИИМ».
К перспективным в данном направлении технологиям относят:
- гальваническое хромирование;
- комбинированное покрытие: гальванический хром + лазерное упрочнение на незащищенной поверхности в меньшей степени подверженной воздействию ПГ;
- вакуумная ионно-плазменная технология нанесения покрытия Сг, Ж, ВК-8, Ш, МЬЫ, Та, ТаЫ (магнетронное напыление, термическое электроннолучевое, дуговое напыление);
- вакуумное ионно-плазменное магнетрон-ное напыление импульсами высокой мощности (И1Р1М8);
- сварка взрывом (Та-2,5 Ж, Та, Со, стеллит);
- лазерное и плазменное упрочнение поверхностей и др.
Постановка задачи исследования
Применение технологий нанесения ТЗП, обеспечивающего работу изделия в экстремальных условиях ведет к необходимости применения методов имитационного моделирования физико-химических процессов, происходящих при действии ПГ. При этом существенную практическую значимость представляют исследования, проводимые на натурных образцах, изготовленных из
материалов идентичных материалу изделий с ТЗП и с реальными образцами ТЭ, которые должны применяться, как для подтверждения адекватности и верификации математических моделей, так и для сравнительной оценки мероприятий по защите от РЭД.
Многообразие перспективных технологий нанесения защитных покрытий, разработка перспективных ТЭ, сложность механизма воздействия продуктов ТЭ и неизученность их воздействия на поверхности и их защитные покрытия, ставят сегодня задачу широкомасштабных экспериментальных исследований.
Проведение полномасштабных экспериментальных исследований на образцах изделий, работающих в экстремальных условиях не целесообразно. Поэтому прибегают к исследованиям РЭД на баллистических установках, которые будут способны формировать подобный воздействию в реальных условиях импульс ГП на образцах с защитным покрытием. Наиболее экономически целесообразным методом, обеспечивающим такие исследования, является метод, основанный на использовании газодинамической импульсной сопловой установки (ГИСУ) [3]. Метод позволяет формировать импульс ГП в широком диапазоне давлений за счет варьирования диаметром сопла, массой ТЭ и объемом зарядной каморы. Основной проблемой, решаемой при применении ГИСУ, является обеспечение адекватности РЭД ПГ процессам, происходящим в реальных условиях.
Прогноз стойкости защитного покрытия может быть обеспечен использованием расчетно-экспериментальной методики, основанной на воздействии на него потока ПГ, приводящего к эрозионному изнашиванию.
Особенности устройства и действия газодинамической импульсной сопловой установки
Схема ГИСУ приведена на рисунке 1. Схема размещения испытуемых образцов в сопловом блоке ГИСУ приведена на рисунке 2. Основными требованиями, предъявляемыми к ГИСУ, являются следующие [3]:
- обеспечение запаса прочности при рабочих давлениях;
- обеспечение возможности регистрации основных эксплуатационных параметров (рабочее давление, температура на поверхности образца);
- соблюдение эргономических требований и безопасность проведения работ.
Подробное описание устройства ГИСУ (Рис.1) приведено в работе [3]. Инициирование воспламенения заряда 10 производится электрическим импульсом малой мощности, подведенным через канал трубы. При сгорании ТЭ 10 поток ПГ проходит через направляющую втулку 5 и омывает поверхность испытуемых образцов во втулке 5. Здесь же поток ПГ омывает поверхность с установленным безинерционным датчиком температуры. В процессе функционирования ГИСУ регистрируется изменение во времени давления в каморе 3 и температуры поверхности во втулке 5 (поверхности образцов). Сема измерительного комплекса приведена на рисунке 3.
7 3 10 1 5 9
Рисунок 1 - Экспериментальная газодинамическая импульсная сопловая установка: 1 - корпус; 2 -затвор; 3 - камора; 4 - обтюратор; 5 - втулка с образцами; 6 - труба; 7 - кожух; 8 -датчик давления; 9 -датчик температуры; 10 - ТЭ; 11 - запальная линия
Основу измерительного комплекса составляет модуль USB-3000, который позволил построить мобильную систему сбора и обработки цифровой информации, получаемой от датчика давления и термопар (т/п) через комутатор SCC-68. В качестве измерительного устройства использовался персональный компьютер. Высокое быстродействие модуля USB-3000 позволило исследовать быстропротекающие процессы. Достоверность измерения давления обеспечивается применением аттестованной штатной приборной базы. Достоверность измерения температуры поверхности образцов обеспечивается тарировкой специальной безинерционной термопары, изготавливаемой из материала, аналогичного защитному покрытию. Верхние пределы измерений для
давления - 1000 МПа и для температуры - 1700 К. Для обеспечения моделирования температурных, силовых и временных параметров воздействия ПГ в каморе производится воспламенение и сгорание ТЭ сходных с натурными по химическому и компонентному составу. Для обеспечения подобия процесса применяются специально изготовленные более тонкие ТЭ.
Рисунок 2 - Схема размещения образцов:
1 - сопловой блок; 2 - сопло; 3 - испытуемые образцы; 4 - втулка крепления образцов
Имитационное моделирование процесса теплового воздействия
Методика обеспечения баллистического подобия процесса основана на подборе объема каморы, массы ТЭ и толщины горящего свода элементов, обеспечивающих:
- наибольшее давление ПГ ртах в каморе 3 установки (Рис.1) - на уровне штатной моделируемой системы;
- длительность нарастания импульса давления тк- на уровне штатной моделируемой системы;
- 5и - площадь импульса теплового воздействия - на уровне штатной моделируемой системы.
В основе методики баллистического проектирования лежат основные зависимости горения топливного заряда в постоянном объеме с истечением ПГ через небольшое отверстие [3] при условии, что максимум давления ПГ совпадает с концом горения.
Рисунок 3 - Схема измерительного комплекса
Таблица 1 - Результаты подбора параметров процесса
Номер кривой Топливный элемент 4/1
^0, дм3 со, кг Ък, с
1 1 0,49 0,024
2 1,8 0,87 0,042
Топливный элемен 7/7
3 1,0 0,52 0,025
4 1,8 0,89 0,044
Решение задачи анализа температурно-си-лового воздействия при параметрах заряжания, представленных в таблице 1. приведено на графиках на рисунке 4
Величина температурного импульса обеспечивается подбором объема каморы Шо, за счет установки специальных недеформируемых втулок и массы пороха о. Варьирование производится в пределах (0,5 ^ 2,0) дм3 (Рис. 4). Однако вследствие сложности моделирования процесса истечения в геометрически сложном пространственном канале и в процессе горения заряда ТЭ окончательный подбор заряда (таблица 1) производится экспериментально. Как правило, расчетные значения давлений ртах в каморе отличаются в большую сторону, причем тем больше, чем ниже требуемое расчетное значение.
о
0,01
0,02
Время,с
Рисунок 4 - Импульсы давления (в соответствии с таблицей 1)
В качестве модели температурно-сило-вого воздействия ПГ на поверхность образцов применяется система уравнений, приведенная в работе [3].
В результате решения задачи моделирования температурно-силового воздействия ПГ были
получены кривые давления и температуры ПГ в зависимости от времени.
На рисунке 5 представлены кривые импульса давления ПГ в каморе и на образцах.
Таблица 2 - Параметры имитационного процесса
№ п/п Наименование параметра Значение
1. Диаметр сопла, мм 10
2. Объем каморы, дм3 2,1
3. Масса ТЭ, кг 0,8
4. Плотность заряжания 0,38
5. Теплота горения топлива, ккал/кг 840
6. Конечный импульс топлива, КПахс 500
Очевидно, что кривая скорости потока ПГ эквидестантна температуре ПГ и значительно (почти в два раза) превышает скорость в штатных системах.
и 5.00Е+08 0>
Ч 4.00Е+08 3,00Е+08 2.00Е-Ю8 1.00Е+08 О.ООЕ+ОО
0,01 0,015 0,02
Время,с
Рисунок 5 - Импульс давления: 1 - в каморе, 2 - на образцах, 3 - штатный, 4,5 - с уменьшенным критическим сечением (диаметр 10 мм)
Математическая модель анализа теплового состояния поверхности основана на законе теплопередачи Ньютона-Рихмана, определяемом уравнением
= а(Т1-Тсг)Р. (1)
с учетом изменяющихся во времени коэффициентом теплоотдачи« и температуры внутренней поверхности образцов Тст.
Осредненные по поверхности теплообмена значения коэффициента теплоотдачи при горении топлива в постоянном объеме определяются в соответствии с данными специальных
0,025
опытов Мюраура и Серебрякова по следующей формуле [3]
/ ч 0,65
а* = Стс/Ж0,55(^/Л) , (2)
где Ш - свободный объем каморы, в котором происходит горение топлива; йф/ п
ш / = " _ приток массы газа в результате горения топлива.
Осредненные по поверхности теплообмена значения коэффициента теплоотдачи в канале ГИСУ малых диаметров определяются введением дополнительного члена в формулу (2)
а = а* + Стс/Шо,55"т1«0]0,65, (3)
который учитывает теплообмен при движении ПГ в канале. Здесь коэффициент пт и показатели степени устанавливаются экспериментально, так чтобы получить температуру нагреваемых элементов ГИСУ.
Пример результатов решения задачи теплообмена приведен в виде графиков на рисунках 6 - 8. Эти результаты могут быть приняты в качестве критериев оценки адекватности расчетного метода и предметом согласования модельных и натурных параметров изнашивания образцов при применении мероприятий по защите изнашиваемой поверхности.
160 000 140 000 т 120 000 100 000 80 000 60 000 40 000 20 000 0
й «
н о о
о
о «
__ПО ФОРМУЛЕ 2 - - ПО ФОРМУЛЕ 3
1
/
У
/
0,01 0,02 0,03 0,04 0,05 Время процесса, с
Рисунок 6 - Коэффициент теплоотдачи
Таким образом, сравнение полученных результатов расчетов с заданными вариантами показывает возможность моделирования темпера-турно-силового воздействия на испытуемые образцы на разработанной экспериментальной ГИСУ.
0
Рисунок 7- Тепловой импульс: 1 - в каморе, 2 - на образцах, 3 - штатный, 4,5 - с уменьшенным критическим сечением (диаметр 10 мм)
Температура образца определяется из решения уравнения теплопроводности с учетом изменения температуропроводности материала образца от температуры а(Т). Для решения задачи необходимо предварительно определить значения коэффициентов теплоотдачи на внутренней аг и наружной анж поверхностях образца, а также температуры потоков изнутри Тг и снаружи Тнж. Подобные задачи можно решать с использованием универсальных программных продуктов [4].
тирующими коэффициентами, характеризующими результат теплового и динамического воздействия С^, времени теплового действия Ст, давления ПГ Ср, скорости потока Cv, массы топлива Сш (удельной тепловой энергии на поверхности образцов) в соответствии с условиями физического подобия (Fo = idem; Bi = idem). (Здесь обозначено % - глубина изношенного слоя.)
Такие соотношения могут быть получены на основании рассмотрения уравнений:
- теплопроводности
- изнашивания:
срТсТ= a(Tt- ТсТ) +Л—1 х=$;
- уравнения баланса энергии при функционировании ГИСУ;
- зависимости коэффициента теплоотдачи от давления и скорости ПГ и др.
Основным критерием теплового подобия является равенство температуры защищаемой поверхности КС и образцов на глубине адгезионного слоя. Как правило, именно разрушение этого слоя приводит к разрушению защитного покрытия. Причем верхним критерием является достижение температуры в этом слое критической точки аустенитного превращения Ас\ в стали.
Отдельной задачей исследования является определение динамического эффекта от действия скоростного напора ПГ, скорость которого значительно выше, чем скорость в реальных системах и практически не зависит от давления, поскольку образцы располагаются в критическом сечении сопла.
0,00Е+00 5,ООЕ-О3 1,0ОЕ-О2
Время, с
1,50Е-0
Рисунок 8 - Температура поверхности защитного покрытия образца
Традиционным подходом при переносе результатов модельных испытаний на полномасштабные ГИУ является - основанный на определении соотношения между масштабными резуль-
Результаты проведения стендовых испытаний эрозионной стойкости ТЗП
На рисунках 9-12 приведены внешний вид образцов с ТЗП, прошедших испытания на ГИСУ (рис.1, 2). Результаты проведенных испытаний являются сравнительными, а также могут применяться для параметризации моделей прогнозирования стойкости защитных покрытий применительно к конкретной гетерогенной системе потока ПГ и материала защитного покрытия.
Испытанию подвергались образцы с тантал-вольфрамовым (TaW10) покрытием, нанесенным методом ионно-плазменного вакуумного магнетронного осаждения, с хромовым гальваническим покрытием» и др..
Рисунок 9 - Внешний вид образцов перед проведением испытаний
Рисунок 1 0- Внешний вид образца с хромовым гальваническим покрытием после 1-го нагружения
Рисунок 11 - Внешний вид образца с покрытием TaW10 после 4-х нагружений
Рисунок 12 - Внешний вид образца с покрытием TaW10 после 7-ми нагружений
На всех образцах видны следы схода ТЗП после испытаний. Измерение площади и характера схода ТЗП позволило в результате расчетных и экспериментальных исследований получить сравнительную характеристику стойкости различных видов покрытий и технологий их нанесения, а также дать характеристику их применимости по отношению к условиям эксплуатации.
Заключение
1. Разработана газодинамическая импульсная сопловая установка, предназначенная для испытания стойкости защитных покрытий. Проведены экспериментальные исследования, обеспечивающие соответствующие реальным основные
характеристики температурно-силового воздействия на поверхности образцов с защитными покрытиями.
2. Разработано программно-методическое обеспечение и проведен анализ характеристик газодинамического процесса. Установлено, что экспериментальная установка позволяет моделировать баллистический процесс, происходящий в газодинамических импульсных устройствах в широком диапазоне температурно-силового воздействия ПГ, позволяющем проводить в жестких условиях испытания эрозионной стойкости образцов и их ТЗП.
3. Проведены испытания образцов, изготовленных из легированной стали с нанесенными различными способами ТЗП.
Литература
1. Лепеш, Г. В. Оценка напряженно-деформированного состояния термозащитного покрытия канала трубы, нагруженной высокотемпературным силовым импульсом / Г. В. Лепеш, М. В. Басова // Технико-технологические проблемы сервиса. - 2023. - № 4(66). -С. 31-37. - ББМ ^РТЖ.
2. Лепеш, Г. В. Оценка теплового состояния стального цилиндра с термозащитным покрытием, нагруженного высокотемпературным тепловым импульсом / Г. В. Лепеш // Технико-технологические проблемы сервиса. - 2023. - № 3(65). - С. 33-39. - ББМ ХБЛИР!
3. Лепеш, Г.В. Разработка и обоснование метода экспериментального исследования стойкости антиэрозионных защитных покрытий газодинамических импульсных устройств. /Г.В. Лепеш, Д.Ю.Латышев, М.С.Черкасов // Технико-технологические проблемы сервиса. -2015. - №2(28), С.59- 66.
4. Лепеш, Г.В. Моделирование процесса тепломассо-переноса в программной среде Ansys/Fluent при дифференцированном отоплении подземного перехода./ Г.В. Лепеш, Т.В. Потемкина, Г.А. Спроге // Технико-технологические проблемы сервиса. -2015 № 4(34) с. 41-49.
5. Лепеш Г.В. Моделирование процесса нагружения трубы внутренним давлением с перемещающимся с высокой скоростью фронтом нагружения. // Сб.докла-дов X межд. конф. по мягким вычислениям и измерениям. 25-27 июня 2007 г. СПб.: ЛЭТИ. С.152 - 157.