Теплофизика высоких температур, 2020, T. 58, № 1, стр. 123-127
Интенсификация процесса очистки микропористого оксида алюминия от примеси железа с помощью продувки аргоном
Е. И. Школьников 1, П. П. Иванов 1, *
1 Объединенный институт высоких температур РАН
Москва, Россия
* E-mail: peter-p-ivanov@yandex.ru
Поступила в редакцию 25.02.2019
После доработки 12.08.2019
Принята к публикации 22.10.2019
Аннотация
Предложена модель интенсификации процесса высокотемпературной вакуумной очистки микропористого оксида алюминия от примеси железа путем продувки слоя частиц аргоном. В характерном для этого процесса диапазоне температур давление насыщающих паров над чистым железом имеет порядок 0.1 мм рт. ст., а для примеси с концентрацией несколько частиц на миллион цифра пропорционально меньше, и это уже создает проблемы с вакуумированием камеры. Использование аргона как несущего газа разрушает жесткую связь между парциальным давлением паров примеси и давлением в вакуумной камере. Расчеты показывают, что при умеренном вакууме 10–5 мм рт. ст. поток аргона с удельным массовым расходом порядка 10–7 кг/(м2 с) обеспечивает весьма большую глубину очистки.
ВВЕДЕНИЕ
Вакуумная очистка оксида алюминия от примеси железа основана на том, что железо испаряется быстрее оксида алюминия. Если взять оксид с примесью железа 400 мг/кг и привести эту систему в равновесие при температурах 1600 и 1750°С, то состав газовой фазы, по данным ИВТАНТЕРМО [1], будет таким, как показано в таблице. Пары железа превалируют с отрывом минимум на два порядка. Очистка происходит путем удаления паровой фазы вакуумным насосом.
С помощью данных из базы ИВТАНТЕРМО и программы расчета равновесного состава можно найти зависимость давления насыщающих паров железа над чистым веществом от температуры. В диапазоне температур 1200–2000°С эта зависимость аппроксимируется формулой
(1)
$\begin{gathered} p = {\text{exp}}(--8.2575{{x}^{2}} + 40.009x--40.769), \\ x = {T \mathord{\left/ {\vphantom {T {1000}}} \right. \kern-0em} {1000}} \\ \end{gathered} $Давление паров малой примеси в [2] рекомендуется определять как произведение давления над чистым веществом на мольную долю примеси. В любом случае это давление мало, и требуется оценить длину свободного пробега атомов железа, чтобы понять характер течения паров через слой гранул оксида алюминия. Если считать атомы железа сферами диаметром $d,$ то длина их свободного пробега, согласно молекулярной теории, определяется формулой [3]
В (2) $p$ и $T$ – давление и температура, $k$ – постоянная Больцмана. Диаметр атома железа можно принять на уровне $d$ = 0.3 × 10–9 м (например, [4]). В этом случае минимально возможное значение длины свободного пробега (при $T$ = 1200 + + 273, $p$ = 52 × 0.001 – учет мольной доли примеси) составит чуть меньше 1 м, т.е. намного больше эффективного радиуса пор между частицами подвергаемого очистке оксида. Это значит, что столкновениями атомов железа между собой можно пренебречь по сравнению с их столкновениями с эффективными стенками пор, т.е. это кнудсеновское течение. Для атомов азота (приблизительно того же диаметра) при тысячекратном превышении давления над давлением умеренного вакуума (10–5 Тор) длина свободного пробега составляет 38 мм, что на порядок больше размера гранул.МАТЕМАТИЧЕСКАЯ МОДЕЛЬ
Одномерное уравнение сохранения импульса для кнудсеновского течения по оси Z – от дна к поверхности слоя частиц очищаемого оксида алюминия – имеет вид [5]
Предположим, что в случае, когда поток ${{J}_{Z}}$ представляет собой поток смеси аргона и паров железа $({{J}_{Z}} = J_{Z}^{{{\text{Fe}}}} + J_{Z}^{{{\text{Ar}}}}),$ уравнение сохранения импульса в виде (3) применимо, если под давлением понимать суммарное давление ${{p}^{\Sigma }} = {{p}_{{{\text{Fe}}}}} + {{p}_{{{\text{Ar}}}}}{\text{:}}$
(5)
$\frac{{d{{p}^{\Sigma }}}}{{dz}} = - \frac{{3{{{v}}_{T}}}}{{8R_{Z}^{p}}}({{J}_{Z}} + {{J}_{{{\text{Ar}}}}}).$(6)
${{p}_{{{\text{Fe}}}}} = \frac{{{{J}_{Z}}}}{{{{J}_{Z}} + {{J}_{{{\text{Ar}}}}}\frac{{{{M}_{{{\text{Fe}}}}}}}{{{{M}_{{{\text{Ar}}}}}}}}}{{p}^{\Sigma }}.$Тепловую скорость в уравнении (5) можно вычислять по среднему молекулярному весу. Вопрос осреднения не принципиален при наличии в знаменателе экспериментально идентифицируемого параметра эффективного радиуса зазора между гранулами $R_{Z}^{p}.$
Система обыкновенных дифференциальных уравнений (4), (5) относительно фазовых координат ${{J}_{Z}}$ и ${{p}^{\Sigma }}$ описывает процесс выхода паров железа, поступающих через поверхность гранул в пространство между ними (поток ${{J}_{Y}}$), из толщи слоя гранул на поверхность в вакуумируемое пространство.
В правую часть (4) входит поток ${{J}_{Y}},$ представляющий собой функцию ${{J}_{Y}} = f({{p}_{{{\text{Fe}}}}}),$ за которой стоит решение сложной задачи об испарении примеси из частицы. Простейший аналитический подход к решению этой задачи для сферических микропористых гранул оксида алюминия, сложенных из сферических микрочастиц, предложен в [6]:
(8)
$J_{{{\text{Fe}}}}^{{\text{1}}} = {{c}_{1}}\rho _{{{\text{Fe}}}}^{{\text{1}}} - {{c}_{2}}{{p}_{{{\text{Fe}}}}},$(9)
$J_{{{\text{Fe}}}}^{{\text{0}}} = {{c}_{1}}\rho _{{{\text{Fe}}}}^{{\text{0}}} - {{c}_{2}}p_{{{\text{Fe}}}}^{0},$Здесь $\rho _{{{\text{Fe}}}}^{{\text{0}}},$ $\rho _{{{\text{Fe}}}}^{{\text{1}}}$ – плотность примеси в центре и на периферии пористой частицы (кг/м3); $J_{{{\text{Fe}}}}^{0},$ $J_{{{\text{Fe}}}}^{{\text{1}}}$ – поток пара с поверхности микрочастицы, находящейся в центре и на периферии пористой гранулы; $p_{{{\text{Fe}}}}^{{\text{0}}}$ – давление пара в центре пористой гранулы. Уравнения (7)–(10) позволяют в неявном виде определить зависимость ${{J}_{Y}} = f({{p}_{{{\text{Fe}}}}},\rho _{{{\text{Fe}}}}^{{\text{0}}},\rho _{{{\text{Fe}}}}^{{\text{1}}}),$ а значения $J_{{{\text{Fe}}}}^{{\text{0}}}$ и $J_{{{\text{Fe}}}}^{{\text{1}}}$ являются побочным результатом, используемым для нахождения временнóго шага при решении уравнения
(11)
$\begin{gathered} \frac{d}{{dt}}\left( {\rho _{{{\text{Fe}}}}^{i}\frac{4}{3}\pi {{R}^{3}}} \right) = - J_{{{\text{Fe}}}}^{i}4\pi {{R}^{2}} \to \\ \to \frac{{d\rho _{{{\text{Fe}}}}^{i}}}{{dt}} = - \frac{3}{R}J_{{{\text{Fe}}}}^{i},\,\,\,\,i = 0,1. \\ \end{gathered} $Константы задачи в уравнениях (7)–(10) равны
(12)
$\begin{gathered} {{c}_{1}} = \frac{1}{{\frac{{\rho {\text{*}}}}{{J{\text{*}}}} + \frac{R}{{2D}}}},\,\,\,\,{{c}_{2}} = \frac{{\rho {\text{*}}}}{{J{\text{*}}{{{v}}_{T}}}}{{c}_{1}},\,\,\,\,{{c}_{3}} = \frac{{3{{{v}}_{T}}{{R}_{y}}}}{{16R_{Y}^{p}}}, \\ {{c}_{4}} = \frac{{{{R}_{y}}N_{S}^{0}}}{6}, \\ \end{gathered} $(13)
$J{\text{*}} = \frac{{\rho {\text{*}}}}{{{{\rho }_{{{\text{A}}{{{\text{l}}}_{{\text{2}}}}{{{\text{O}}}_{{\text{3}}}}}}}}}\frac{{p_{{{\text{Fe}}}}^{*}}}{{{{{v}}_{T}}}}.$ИЛЛЮСТРАЦИЯ РЕЗУЛЬТАТОВ ИНТЕГРИРОВАНИЯ
Расчет процесса очистки гранулированного оксида алюминия от примеси железа проведен с помощью системы уравнений (4)–(13) при следующих исходных данных: температура – 1750°С; вакуум (pvac) – 1 × 10–5 Тор; толщина слоя гранул – 4 см; радиус пористых гранул – 1 мм; радиус микрочастиц – 0.8 мкм; начальная концентрация примеси железа – 450 ppm; длительность процесса сублимационной очистки – 120 мин.
На рис. 1 приведены кривые зависимости остаточной средней концентрации примеси железа от длительности процесса при различных значениях потока аргона. Кривая 1 рассчитана в отсутствие потока аргона. Кривая 2 соответствует потоку аргона 1 × 10–9 кг/(м2 с). Это начало влияния потока аргона на процесс. Кривые 1 и 2 практически совпадают. Кривые 3 и 4, соответствующие потокам аргона 1 × 10–6 и 1 × 10–5 кг/(м2 с), демонстрируют положительное влияние потока аргона на процесс очистки.
Процесс очистки во времени иллюстрируется рис. 2 и 3. На рис. 2 показаны кривые зависимости плотности примеси в центре гранул от координаты z по толщине слоя через каждые 12 мин процесса для двух случаев: без продувки и с потоком аргона 1 × 10–5 кг/(м2 с). Обращает на себя внимание тот факт, что благодаря продувке аргоном процесс очистки идет не только с задней, но и с передней поверхности слоя.
На рис. 3 показаны кривые зависимости давления (в единицах, кратных давлению вакуумной камеры) от координаты z по толщине слоя через каждые 12 мин процесса. На рис. 3а и 3б – это давление паров примеси, на рис. 3в – суммарное давление (аргона и паров примеси), на рис. 3г – суммарное давление и давление паров примеси в начальный момент времени.
ПАРАМЕТРИЧЕСКИЕ ИССЛЕДОВАНИЯ
На рис. 4 приведена зависимость остаточной средней концентрации примеси железа от длительности процесса при различных значениях параметров процесса: температуры, давления в вакуумной камере и потока аргона. Вариация температуры от 1600 до 1750°С разделяет кривые на группы медленного (1–4) и быстрого (5–8) изменения концентрации. При температуре 1600°С кривые 1 и 2, соответствующие расчету без аргона при двух значениях давления – 2 × 10–5 и 1 × 10–5 Тор, показывают сильную зависимость процесса от глубины вакуума. Кривые 3 и 4, соответствующие расчету с потоком аргона 1 × 10–5 кг/(м2 с) при тех же значениях давления вакуума, практически совпадают. Аналогичная картина наблюдается и при температуре 1750°С. Кривые 5 и 6, соответствующие расчету без аргона при двух значениях давления – 1 × 10–5 и 1 × 10–6 Тор, показывают зависимость процесса от глубины вакуума, в то время как кривые 7 и 8, соответствующие расчету с потоком аргона 1 × 10–5 кг/(м2 с) при тех же значениях давления вакуума, практически совпадают.
ЗАКЛЮЧЕНИЕ
Предложена модель интенсификации процесса высокотемпературной вакуумной очистки микропористого оксида алюминия от примеси железа с помощью продувки слоя частиц аргоном.
Показано, что давление в вакуумной камере складывается из давления аргона и пара железа, а парциальное давление пара железа зависит только от температуры и концентрации примеси, благодаря чему уровень вакуума перестает быть ограничителем глубины очистки.
Зависимость остаточной средней концентрации примеси железа от потока аргона при прочих равных условиях не имеет экстремума, чем больше аргона, тем меньше примеси. Но в целях минимизации затрат целесообразно ограничиться минимумом производной этой функции, который в рассмотренных условиях имеет место при величине потока аргона 3.443 × 10–7 кг/(м2 с).
Список литературы
Белов Г.В., Иориш В.С., Юнгман В.С. Моделирование равновесных состояний термодинамических систем с использованием ИВТАНТЕРМО для Windows // ТВТ. 2000. Т. 38. № 2. С. 209.
Powell A., Van Den Avyle J., Damkroger B., Szekely J., Pal U. Analysis of Multicomponent Evaporation in Electron Beam Melting and Refining of Titanium Alloys // Metall. Mater. Trans. 1997. V. 38B. P. 1227.
Чепмен С., Каулинг Т. Математическая теория неоднородных газов. М.: Изд-во иностр. лит., 1960. 510 с.
Бекман И.Н. Мембраны в медицине. http://profbeckman.narod.ru/MedMemb.files/medmemb5.pdf
Welty J.R., Wicks C.E., Wilson R.E., Rorrer G.L. Fundamentals of Momentum, Heat and Mass Transfer. 5th ed. Hoboken: John Wiley and Sons, 2008.
Иванов П.П., Жук А.З., Власкин М.С. Высокотемпературная вакуумная очистка микропористого оксида алюминия от примеси железа // Вестн. ОИВТ РАН. 2018. Т. 1. С. 61.
Doremus R.H. Diffusion in Alumina // J. Appl. Phys. 2006. V. 100. № 10. P. 101 301.
Дополнительные материалы отсутствуют.
Инструменты
Теплофизика высоких температур