The effect of initial temperature on the modes of passivation of pyrophoric na-nopowder (macrokinetic approach)

Cover Page

Abstract


In the work, thermal modes of passivation of pyrophoric nanopowders at lower initial temperatures of the layer of the nanopowder are studied with analytical and numerical methods in the context of a macrokinetic approach. It is shown that the decrease in the initial temperature of the layer leads to a transition from a layer-by-layer mode of passivation to the volume one. That allows controlling the acceptable level of a warming up at passivation even at the high concentration of oxidizer in gas. Analytical expressions for delimitation of the layer-by-layer and volume passivation modes are obtained; these are confirmed by the results of numerical calculations.


Full Text

В литературе выделяют два основных фактора, которые учитываются при описании наносистем. Это изменение термодинамического состояния наносистем по сравнению с классическим состоянием и появление квантово-размерных эффектов с уменьшением характеристических размеров структурных элементов [1]. Наносистемы представляют собой макроскопические ансамбли малых частиц с размерами до нескольких нанометров. Исследования термодинамики малых частиц показывают, что размер частиц является переменной, определяющей вместе с другими термодинамическими переменными состояние системы [2], т. е. классические термодинамические представления могут быть применимы к ансамблям наночастиц. Совокупность экспериментальных данных, полученных в [3, 4], свидетельствует о том, что для засыпок из нанопорошков железа существуют критические условия теплового воспламенения, связанные с изменением соотношения между тепловыделением и теплоотводом. Следовательно, теоретические подходы классической макроскопической теории теплового взрыва могут быть применимы для описания явлений воспламенения в макроскопических объектах, состоящих из наночастиц.

Как показали результаты наших экспериментальных исследований [3, 5], одним из перспективных направлений интенсификации процесса пассивации нанопорошков является проведение пассивации при температурах ниже нуля. Существует ряд публикаций, посвященных проблемам получения и стабилизации нанопорошков, однако эти работы имеют преимущественно прикладной характер [6, 7]. В данной работе в рамках макрокинетического подхода аналитическими и численными методами исследовано влияние начальной температуры пирофорного нанопорошка на режимы пассивации.

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

 

Рис. 1. Модель процесса. 1– газ над слоем порошка, 2 – пассивированный нанопорошок, 3 – исходный нанопорошок, 4 – подложка.

 

Предполагается, что химическое взаимодействие пористого вещества с окислителем происходит за счёт диффузионного подвода окислителя из внешней среды (через верхнюю поверхность слоя), а продукты взаимодействия конденсированные. Кроме того, считается, что температуры пористого вещества и газа в порах равны (однотемпературная модель). Считается, что отвод тепла в окружающую среду происходит только через нижнюю границу слоя и описывается законом Ньютона. Предполагается, что при достижении некоторой степени превращения пористого вещества (которая может быть определена на основании экспериментальных исследований) на поверхности частиц образуется защитный слой, который препятствует дальнейшему окислению частиц. В принятых в теории горения безразмерных переменных [8, 9] математическая модель, описывающая такой процесс, имеет вид

θτ=2θξ2+f(θ,η), ητ=Le2ηξ2+γfθ,η,  

η1τ=νγfθ,ηfθ,η=1ηexpθ1+βθ при η<1 è η1<η1,k,0 при η=1 или η1=η1,k.(1)

 

Начальные и граничные условия:

τ=0:θ=θin,η=ηn,η1=0;τ>0:ξ=0:θξ=0,η=0;ξ=Lξ:θξ=Bi θinθ,ηξ=0.  

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

 T0t=ccRT02EQρcρga0k0expERT0, x=λctccρc.

Соответственно, безразмерные переменные и параметры выражаются следующим образом:

θ=ERT02TT0,τ=tt, ξ=xx, Lξ=lxx,β=RT0E,

γ=ccRT02EQρcρga0, η=a0aa0, η1=ρc0ρcρc0,

θin=ERT02TinT0,  Bi=αxλc,  Le=Defaс,

ac=λcccρc, ν=ν¯ρga0ρc.

Здесь размерные величины: t — время; x — координата; a0, a — исходная и текущая концентрация газообразного окислителя; ρc0,ρc – исходное и текущее содержание конденсированного вещества в единице объёма; ρg — плотность газовой фазы; cc – теплоёмкость конденсированного вещества; R — универсальная газовая постоянная; E и k0 — энергия активации и предэкспоненциальный множитель гетерогенной реакции, отнесенный к единице объёма пористой среды; Q — тепловой эффект гетерогенной реакции, отнесённый к единице массы окислителя, λc, ac  — коэффициенты теплопроводности и температуропроводности конденсированной фазы; α — коэффициент теплоотдачи; Def — эффективный коэффициент диффузии окислителя; ν¯ – количество конденсированного вещества, реагирующего с одним граммом окислителя, t=tch — характерное время химической реакции, x — характерная длина, lx — толщина слоя засыпки. Безразмерные переменные и параметры: τ, ξ — время и координата; θ — температура; θin – начальная температура засыпки, равная температуре окружающей среды; η — глубина превращения окислителя; η1 — глубина превращения конденсированного вещества; η1,k — максимальная степень превращения конденсированного вещества; ηn — начальная глубина превращения окислителя внутри засыпки (введённая для удобства численного счёта и не влияющая на его результат); Le — аналог параметра Льюиса для окислителя; β, γ, ν – параметры, Bi – число Био; Lξ – толщина слоя засыпки.

В качестве температуры T0 удобно выбрать комнатную температуру, так как при этой температуре проводится пассивация в обычных условиях (при малых концентрациях окислителя в газе). Тогда пассивации при пониженных температурах соответствует отрицательное значение θin.

В сформулированной задаче существует два характерных размера: это толщина слоя Lξ и глубина проникновения реакции в пористое тело Lcem при температуре θin. В соответствии с [8] глубина проникновения реакции в пористое тело с точностью до множителя порядка единицы равна

Lcem=Leγ  1ηξ=0nexpθin     Leγexpθin. (2)

При записи (2) пренебрегли слабо влияющей на условия воспламенения зависимостью скорости реакции от параметра  [8, 9]. Как следует из (2), при понижении начальной температуры слоя θin экспоненциально возрастает Lcem и может стать равным толщине слоя. Данный факт является необходимым условием проведения процесса пассивации в объёмном режиме.

Если процесс пассивации лимитируется скоростью подвода окислителя, то этот процесс происходит во фронтальном режиме [10]. В этом режиме зависимость координаты фронта ξfr и его скорости ufr от времени τ определяются выражениями

ξfr=2Leνη1,kτ, ufr=Leν2η1,kτ. (3)

 

Для температуры фронта θfr в [10] было получено следующее выражение:

θfr=η1,kγν1+2πη1,kLeν1. (4)

Отметим, что анализ режимов пассивации проводился в [10] при нулевой начальной температуре. Режимы пассивации при отрицательных начальных температурах (θin<0) требуют отдельного рассмотрения, так как изменение начальной температуры приводит к возникновению качественно новых режимов пассивации.

Анализ системы (1) показал, что с помощью замены переменных

θn=θθin, τn=τexpθin, ξn=ξexpθin2

задача о пассивации при температуре, равной θin, сводится к задаче о пассивации при температуре, равной нулю, что позволяет воспользоваться некоторыми результатами, полученными ранее [10]. При такой замене переменных задачу (1) надо рассматривать при новых значениях начальных и граничных условий:

τn=0:θinn=θinθin=  0,  η=ηn,η1=0, (5)

τn>0:ξn=0:θnξn=0,η=0,    ξn=Lξexpθin2:θnξn=expθin2Bi θn,ηξn=0. (6)

 

Как видно из (5), (6), в новых переменных имеет место экспоненциальная зависимость размеров образца Lξexp(θin/2) и эффективного значения критерия Био Bin=expθin2Bi от начальной температуры слоя θin, причем понижение температуры ведёт к уменьшению размеров слоя и увеличению Bin, так как θin<0.

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

Lcem=Leγ. (7)

В новых переменных, если процесс пассивации лимитируется скоростью подвода окислителя и происходит во фронтальном режиме, температура во фронте определяется выражением (4). В старых переменных температура диффузионного режима θfr линейно зависит от θin, т. е.

θfr=θfrθin=0+θin. (8)

Как следует из (6), при понижении начальной температуры размер образца в новых переменных Lξexpθin2 экспоненциально убывает и может стать равным глубине проникновения реакции в пористое тело слоя Lcem. Данное равенство является необходимым условием проведения процесса пассивации в объёмном режиме. Достаточным условием проведения процесса пассивации в объёмном режиме является отсутствие значительных разогревов. Для выполнения неравенства Lξexpθin2<Lcem в процессе химического взаимодействия необходимо, чтобы разогревы порошка были невелики. В соответствии с теорией теплового взрыва для отсутствия значительных разогревов необходимо, чтобы размер образца (при больших значениях критерия Био) был меньше единицы [8]. Поэтому равенство Lξexpθin2=1 позволяет определить значение температуры, до которой нужно охладить образец, чтобы разогрев в процессе химического взаимодействия не превышал один характеристический интервал:

θin=2ln1/Lξ. (9)

Итак, как показывает анализ, значение температуры, ниже которой реализуется объёмный режим пассивации и разогрев, не превышает один характеристический интервал (в размерных переменных — нескольких градусов) и зависит от толщины слоя нанопорошка. Следовательно, в образцах разной толщины при одной и той же отрицательной температуре могут реализоваться качественно различные режимы пассивации с отличающейся величиной разогрева.

 

Рис. 2. Зависимости на нижней границе слоя от времени: глубины превращения нанопорошка η1 (1), окислителя η (2) и максимальной температуры в слое maxθ (3). Диффузионный (фронтальный) режим пассивации. Начальная температура слоя

 

Рис. 3. Зависимости на нижней границе слоя от времени: глубины превращения нанопорошка η1 (1), окислителя η (2) и максимальной температуры в слое maxθ (3), скорос- ти тепловыделения f(Lζ) (4). Объёмный режим пассивации. Начальная температура слоя

 

Численные расчёты позволили проследить динамику процесса пассивации при различных начальных температурах, а также подтвердили правильность выводов теоретического анализа (см. рис. 2, 3). Расчёты проводились при следующих значениях определяющих параметров: β=0,01; γ=0,2; Bi=20; Le=2; Lξ=15, ν=0,06; η1,k=0,1; ηn=0,8; θign=5. При высоких начальных температурах, когда (Lcem<Lξ) реализуется послойный режим пассивации с большим максимальным разогревом. Процесс химического реагирования нанопорошка с окислителем, содержащимся в газе, происходит в тонком слое засыпки. Как видно из рис. 2, по мере приближения к нижней границе слоя температура в зоне реакции (кривая 3) уменьшается из-за отвода тепла в подложку, что приводит к увеличению размеров зоны реакции. Превращение нанопорошка на нижней границе слоя (кривые 1, 2) начинается с большой задержкой и происходит при температурах, близких к начальной температуре. Следовательно, в послойном режиме реагирования различные слои нанопорошка реагируют при разной температуре, что в реальном эксперименте может приводить к различной толщине окисного слоя на наночастицах в засыпке.

В другом предельном случае (низкая начальная температура, Lcem>Lξ) реализуется объёмный режим пассивации (рис. 3). В этом режиме максимальный разогрев при отсутствии теплового воспламенения не превышает один характеристический интервал (кривая 3), а газообразный окислитель находится в каждой точке объёма образца. Особенно наглядно объёмный характер пассивации проявляется при наблюдении за динамикой изменения концентрации окислителя (кривая 2) и степени превращения нанопорошка на нижней границе слоя (кривая 1). Процесс превращения начинается практически в момент контакта слоя порошка с газом, содержащим окислитель, концентрация окислителя близка к концентрации окислителя в окружающем газе и скорость превращения (кривая 4) почти постоянна. Следовательно, в объёмном режиме различные слои нанопорошка реагируют при одинаковой температуре, что в эксперименте может обеспечить одинаковую толщину окисного слоя на наночастицах при пассивации.

Результаты теоретического анализа качественно согласуются с экспериментальными данными, полученными нами при исследовании динамики разогрева образца в процессе воспламенения и пассивации нанопорошков железа и никеля при изменении начальной температуры [3, 11].

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

Исследование выполнено при финансовой поддержке Российского научного фонда (проект № 16–13–00013).

About the authors

B. S. Seplyarskii

Merzhanov Institute of Structural Macrokinetics and Materials Science

Author for correspondence.
Email: seplb1@mail.ru

Russian Federation, Chernogolovka Moscow region

N. M. Rubtsov

Merzhanov Institute of Structural Macrokinetics and Materials Science

Email: nmrubtss@mail.ru

Russian Federation, Chernogolovka Moscow region

T. P. Ivleva

Merzhanov Institute of Structural Macrokinetics and Materials Science

Email: nmrubtss@mail.ru

Russian Federation, Chernogolovka Moscow region

M. I. Alymov

Merzhanov Institute of Structural Macrokinetics and Materials Science

Email: nmrubtss@mail.ru

Russian Federation, Chernogolovka Moscow region

References

  1. Суздалев И.П. Нанотехнология: физико-химия нанокластеров, наноструктур и наноматериалов. М.: КомКнига, 2006. 592 с.
  2. Азаренков Н.А., Веревкин А.А., Ковтун Г.П., Литовченко С.В. Нанотехнологии и наноматериалы. Харьков: Изд-во Казах. нац. ун-т им. аль-Фараби, 2009. 68 с.
  3. Alymov M.I., Rubtsov N.M., Seplyarskii B.S., Zelensky V.A., Ankudinov A.B. // Mendeleev Communs. 2017. № 27. P. 482–484.
  4. Hosokawa M., Nogi K., Naito M., Yokoyama T.Nanoparticle Technology Нandbook. L.: Elsevier, 2007. 644 p.
  5. Rubtsov N.M., Seplyarskii B.S., Alymov M.I. Ignition and Wave Processes in Combustion of Solids. B.: Springer Intern. Publ., 2017. 284 p.
  6. Meziani M.J., Bunker C.E., Lu F., Li H., Wang W., Guliants E.A., Quinn R.A., Sun Y.-P. // ACS Appl. Mater. Interfaces. 2009. V. 1. № 3. P. 703–709.
  7. Nagarajan R., Hatton T.A. Nanoparticles: Synthesis, Stabilization, Passivation, and Functionalization. Wash. (D.C.): Amer. Chem. Soc., 2008. 451 p.
  8. Франк-Каменецкий Д.А. Диффузия и теплопередача в химической кинетике. М.: Наука, 1987. 491 с.
  9. Мержанов А.Г, Барзыкин В.В., Абрамов В.Г. // Хим. физика. 1996. Т. 15. № 6. С. 3–44.
  10. Seplyarsky B.S., Ivleva T.P., Alymov M.I. // Nanotechnol. in Russia. 2017. V. 12. № 11/12. P. 583–588.
  11. Alymov M.I., Rubtsov N.M., Seplyarskii B.S., Ze-lensky V.A., Ankudinov A.B. // Mendeleev Communs. 2017. № 27. Р. 631–633.

Supplementary files

Supplementary Files Action
1.
Fig. 1. Process Model. 1– gas above the powder layer, 2 — passivated nanopowder, 3 — initial nanopowder, 4 — substrate.

Download (53KB) Indexing metadata
2.
Fig. 2. Dependencies at the lower boundary of the layer on time: the depth of transformation of the nanopowder η1 (1), oxidant η (2), and the maximum temperature in the layer maxθ (3). Diffusion (frontal) passivation mode. Initial layer temperature

Download (48KB) Indexing metadata
3.
Fig. 3. Dependencies at the lower boundary of the layer on time: the conversion depth of the nanopowder η1 (1), oxidant η (2) and maximum temperature in the maxθ layer (3), heat release rate f (Lζ) (4). Volume passivation mode. Initial layer temperature

Download (45KB) Indexing metadata

Statistics

Views

Abstract - 238

PDF (Russian) - 148

PlumX


Copyright (c) 2019 Russian Academy of Sciences

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies