Научная статья на тему 'ВЛИЯНИЕ НЕОДНОРОДНОСТИ ТЕМПЕРАТУРНОГО ПОЛЯ НА ПЕРЕНОС ПАРОВ КРЕМНИЯ ОТ ЗЕРКАЛА РАСПЛАВА К ИЗДЕЛИЮ В ХОДЕ ВЫСОКОТЕМПЕРАТУРНОГО СИЛИЦИРОВАНИЯ'

ВЛИЯНИЕ НЕОДНОРОДНОСТИ ТЕМПЕРАТУРНОГО ПОЛЯ НА ПЕРЕНОС ПАРОВ КРЕМНИЯ ОТ ЗЕРКАЛА РАСПЛАВА К ИЗДЕЛИЮ В ХОДЕ ВЫСОКОТЕМПЕРАТУРНОГО СИЛИЦИРОВАНИЯ Текст научной статьи по специальности «Физика»

CC BY
27
6
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
УГЛЕРОДНОЕ ВОЛОКНО / ПАРЫ КРЕМНИЯ / ДИФФУЗИОННЫЙ ТРАНСПОРТ / СРЕДНИЙ ВАКУУМ

Аннотация научной статьи по физике, автор научной работы — Демин В. А., Зинурова В. Э.

Рассматривается модифицированная физико-математическая модель высокотемпературного переноса паров кремния в разреженной несущей среде в условиях среднего вакуума от зеркала расплава до поверхности поглощающего образца пористой среды. В роли пористой среды выступает углеродное волокно искусственного происхождения. Ранее уже были изучены процессы переноса газообразного кремния от зеркала расплава к поверхности изделия в случае однородного распределения температуры в пространстве. Однако в реальном технологическом процессе в общем случае это условие не выполняется. Таким образом, возникает задача исследования влияния этого фактора на интенсивность массопереноса паров кремния в случае навязанного условиями эксперимента неоднородного распределения температуры. Показано, что физико-математическая модель переноса по-прежнему допускает описание процесса с помощью одного нелинейного дифференциального уравнения в частных производных для концентрации.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

INFLUENCE OF THE TEMPERATURE FIELD INHOMOGENEITY ON THE TRANSFER OF SILICON VAPOR FROM THE MELT MIRROR TO THE PRODUCT DURING HIGH-TEMPERATURE SILICONIZATION

We study a modified physical and mathematical model of high-temperature transfer of silicon vapor from the melt mirror to the surface of an absorbing sample of a porous medium that occurs in a rarefied carrier medium under medium vacuum conditions. Carbon fiber of artificial origin acts as a porous medium. The process of gaseous silicon transfer from the melt mirror to the surface of the product was previously studied in the case of a uniform temperature distribution in space. However, in a real technological process this condition is generally not fulfilled. Thus, this phenomenon needs to be studied in the case of an inhomogeneous temperature distribution. It is shown that, as was the case in the previous study, the physical and mathematical model of transfer allows the description of the process with the help of one nonlinear partial differential equation for the concentration.

Текст научной работы на тему «ВЛИЯНИЕ НЕОДНОРОДНОСТИ ТЕМПЕРАТУРНОГО ПОЛЯ НА ПЕРЕНОС ПАРОВ КРЕМНИЯ ОТ ЗЕРКАЛА РАСПЛАВА К ИЗДЕЛИЮ В ХОДЕ ВЫСОКОТЕМПЕРАТУРНОГО СИЛИЦИРОВАНИЯ»

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА

2023

• ФИЗИКА •

Вып. 2

УДК 533.72; 533.59 PACS 51.20.+4 81.90+

Влияние неоднородности температурного поля на перенос паров кремния от зеркала расплава к изделию в ходе высокотемпературного силицирования

В. А. Демин1*, В. Э. Зинурова1'2

1 Пермский государственный национальный исследовательский университет, Пермь

2 Институт механики сплошных сред УрО РАН, Пермь *[email protected]

Рассматривается модифицированная физико-математическая модель высокотемпературного переноса паров кремния в разреженной несущей среде в условиях среднего вакуума от зеркала расплава до поверхности поглощающего образца пористой среды. В роли пористой среды выступает углеродное волокно искусственного происхождения. Ранее уже были изучены процессы переноса газообразного кремния от зеркала расплава к поверхности изделия в случае однородного распределения температуры в пространстве. Однако в реальном технологическом процессе в общем случае это условие не выполняется. Таким образом, возникает задача исследования влияния этого фактора на интенсивность массопереноса паров кремния в случае навязанного условиями эксперимента неоднородного распределения температуры. Показано, что физико-математическая модель переноса по-прежнему допускает описание процесса с помощью одного нелинейного дифференциального уравнения в частных производных для концентрации.

Ключевые слова: углеродное волокно; пары кремния; диффузионный транспорт; средний вакуум

Поступила в редакцию 10.04.2023; после рецензии 29.05.2023; принята к опубликованию 01.06.2023

Influence of the temperature field inhomogeneity on the transfer of silicon vapor from the melt mirror to the product during high-temperature siliconization

V. A. Demin1*, V. E. Zinurova1,2

1 Perm State University, Perm, Russia

2 Institute of Continuous Media Mechanics UB RAS *[email protected]

We study a modified physical and mathematical model of high-temperature transfer of silicon vapor from the melt mirror to the surface of an absorbing sample of a porous medium that occurs in a rarefied carrier medium under medium vacuum conditions. Carbon fiber of artificial origin acts as a porous medium. The process of gaseous silicon transfer from the melt mirror to the surface of the

© Демин В. А., Зинурова В. Э., 2023

распространяется на условиях лицензии

Creative Commons Attribution 4.0 International (CC BY 4.0).

product was previously studied in the case of a uniform temperature distribution in space. However, in a real technological process this condition is generally not fulfilled. Thus, this phenomenon needs to be studied in the case of an inhomogeneous temperature distribution. It is shown that, as was the case in the previous study, the physical and mathematical model of transfer allows the description of the process with the help of one nonlinear partial differential equation for the concentration.

Keywords: carbon fiber; silicon vapor; diffusive transport; medium vacuum

Received 10 April 2023; revised 29 May 2023; accepted 01 June 2023 doi: 10.17072/1994-3598-2023-2-51-58

1. Введение

В качестве одного из перспективных активно развиваемых направлений современного материаловедения выступает совершенствование подходов и технологий производства высокотемпературных композиционных материалов. Типичным примером композиционного материала, который обладает высокими огнеупорными и антиокислительными свойствами, является кремний-углеродное волокно. Технология создания композиционного материала, которая состоит в насыщении объема исходной пористой матрицы парами кремния при высокой температуре, называется силицированием [1]. Однако актуальной проблемой является не только реализация самого процесса насыщения пористой среды связующим, но и доставка паров кремния от зеркала расплава к изделию.

В данной работе рассматривается модельный теоретический подход к описанию доставки паров к изделию, согласно которому газообразный кремний, как некий флюид, фильтруется через несущую среду за счет того, что атомы кремния рассеиваются на точечных, более тяжелых, но, тем не менее подвижных атомах аргона. Ввиду того, что данный процесс осуществляется в условиях среднего вакуума и при чрезвычайно высоких температурах выше точки плавления кремния, считается, что диффузия играет большую, но не решающую роль в обеспечении транспорта газообразного кремния от зеркала расплава к изделию.

2. Теоретическая модель переноса паров кремния

2.1. Система уравнений

Ранее в работе [2] в одномерной постановке были изучены процессы переноса газообразного кремния от зеркала расплава к поверхности изделия в случае однородного распределения температуры в пространстве. Однако в реальном технологическом процессе это условие, как следует из анализа динамики прогрева пространства внутри реторты, не выполняется. Схематически конфигу-

рация технологического процесса представлена на рис. 1. Пары кремния генерируются на правой границе при температуре Т2, переносятся за счет диффузии между атомами несущего газа и полностью поглощаются на левой границе. Таким образом, возникает задача исследования влияния неоднородности теплового поля в реторте в случае каких-либо реалистичных предположений о распределении температуры при математической формулировке задачи.

Сначала рассмотрим проблему в общей постановке. А именно, следуя [3], проанализируем уравнение Навье - Стокса, которое описывает движение вязких жидкостей и газов в пористой среде.

В гидродинамике пористых сред следует отличать поровую скорость V и скорость фильтрации и. Скорость фильтрации определяется через общий расход флюида и связана с поровой скоростью соотношением и = фу . Здесь ф - пористость материала.

• • •

Ti T Т 2

—♦ ♦ ♦ ♦♦

• • • —♦ ♦ ♦ ♦ ♦

1 • • 2

Рис. 1. Двухкомпонентный газ: кружки -атомы несущей среды (аргон), ромбики -атомы кремния; 1 - поверхность образца, 2 - зеркало расплава

Для поровой скорости в среде имеем исходное уравнение движения [3]:

Р/ (Ц + (^)У) = -Ур-Ко, (1)

где р - плотность прокачиваемого флюида, ц -динамическая вязкость, к- проницаемость среды, р - поле давления. В этом уравнении вязкое слагаемое определяется линейной зависимостью от ско-

рости фильтрации. Переходя в уравнении (1) полностью к скорости фильтрации, получаем

р} ф |тV* \ = -Ур -к.

(2)

Как следует из опыта, парциальное давление аргона отвечает условиям среднего вакуума и равно по порядку величины 100 Па. В то же время парциальное давление паров кремния на порядок меньше и составляет всего 10 Па.

2.2. Формула Козени-Кармана

Ввиду крайней разреженности несущей среды, в роли которой выступает аргон, ожидается, что проницаемость к по отношению к фильтрационному потоку паров кремния будет иметь аномально большое значение. Проследим, как выводится формула Козени - Кармана для обоснования ее дальнейшей применимости при оценке проницаемости

к = -

(1 -ф)2

Рис. 2. Модель пористой среды

В теории пористых сред расход через такой материал определяется законом Дарси [4,5]:

ф кБ Ар кБ Ар ру I г I

(3)

Здесь к - проницаемость среды, Б - общая площадь сечения, р - плотность жидкости, у и г - соответственно, кинематическая и динамическая вяз-

кости, Ар - перепад давления на торцах, I - длина массива.

В то же время через один каналец протекает флюид, расход которого определяется формулой Пуазейля [6]. Не вдаваясь в детали множителя перед формулой и предполагая, что по порядку величины длина криволинейного канала в пористом массиве равна длине образца, имеем

г4 Ар

Ч =—Г

г 1

(4)

где г - средний радиус трубочки.

Введем понятие пористости среды ф. Традиционно она определяется так:

V- = Б! = V = V =

N 'г 21 V

(5)

здесь к - константа Козени - Кармана.

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

А именно, пористость - это объем пор V' отнесенный к общему объему массива V. Здесь N ' - число трубочек. По порядку величины Б ' = N'r2 (общая площадь каналов в поперечном сечении). Пусть N - число омываемых препятствий. Заметим, что N и N' - в общем случае являются независимыми величинами. Если а - характерный размер препятствий, то их объем составляет

а ^ = V-V '.

(6)

Отсюда имеем

=v_-v;=V-фv = (1

3 3 3 * '

а

Общая площадь омываемых препятствий в поперечном сечении равна а2К Далее с учетом (5) получаем:

N = Б- = ^

г2 г 21

(8)

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

а 2N =

= N —1 = Ф .

а

(9)

Здесь произведение радиуса канальца на его длину дает по порядку величины площадь его боковой поверхности. Из формулы (9) выражаем радиус канальца:

г=

(1 -Ф)

(10)

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

а

а

Это отличает формулу (10) от моделей, в которых предполагается плотная упаковка препятствий. Тем не менее в пределе плотной упаковки препятствий реализуемый подход тоже дает значения проницаемости, хорошо согласующиеся с известными формулами в теории пористых сред [7].

Полный расход жидкости складывается из всех расходов через отдельные трубочки

д = чы ' = Г! ^Ры'.

V 1

(11)

Принимая во внимание формулу (3) и сравнивая ее с (11), получаем нетривиальное по смыслу соотношение

кЯ = г4Ы'. Отсюда выражаем проницаемость

к = -

г 4ы'

Я

(12)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Подставляем в (12) соотношение (8) и с учетом, что Я' = Ы'г2, получаем

Я 'г2

к=

Отношение площадей всех трубочек в поперечном сечении к общей площади образца дает пористость ф = Б'/Б, а г выражается формулой (10). В итоге приходим к окончательной формуле

к = I

ф3а2

(1 -Ф) (1 -{

(13)

как инерционное слагаемое, так и слагаемое, отвечающее за нестационарность процессов, пренебрежимо малы по сравнению с вязким членом

Г)У

>> Р/ф (иУ)и\

В результате из (2) приходим к выражению для скорости в виде известного закона Дарси [2]

и = —Урк. V

2.3. Учет неоднородности температурного поля

За счет испарения на зеркале расплава и поглощения на изделии имеется некий средний градиент плотности паров кремния. Так как давление идеального газа пропорционально плотности, это порождает градиент давления кремния, который может выступать в роли дополнительной движущей силы помимо диффузии. Согласно уравнению состояния идеального газа:

РкУ я = — тк,

где тк - масса кремния, ц - его молярная масса, Рк - парциальное давление паров кремния, Я -универсальная газовая постоянная, Т - абсолютная температура газа. Следовательно,

Я Т

Рк =-Т Рк.

Ц

(14)

Заметим, что константа Козени - Кармана к в формуле (13) пропорциональна квадрату характерного размера препятствия.

В нашем случае среда представляет собой систему из небольших твердых сферических центров, которые омываются гидродинамическим потоком. Из логики вывода (13) следует, что эта формула вполне применима к нашей модели.

Оценим величину каждого слагаемого в уравнении (1), принимая во внимание, что пористость такой среды в виду ее крайней разреженности близка к единице:

р/ф 2(0У)0

,.2 А12

■Рг — = 2.86 -10-4 — = 0.96 , И/ L 0.5

Л-1 10

р/ф и

^ 0,

Г)У

к

3.5-10-4 • 41 3.110-5

= 463 .

Из этих оценок видно, что вязкое слагаемое является доминирующим в этом уравнении. А именно,

По определению под массовой концентрацией необходимо понимать

С = -

Рк

Ра +Рк

Откуда

С

Рк = 1-С Ра .

Далее следует подставить данное выражение в формулу для парциального давления кремния (14)

Я С Т

Рк = —Г^РаТ . Цк 1 - С

Затем необходимо разложить множитель С/(1 - С) в ряд по малым С, после чего ограничимся в итоговом выражении первым неисчезающим членом:

С 1-С

С,

Рк -Р^С.

Ц

Рассмотрим ситуацию, когда температура неоднородна в пространстве между поверхностью изделия и зеркалом расплава. Для простоты будем

к

3

2

2

в дальнейшем решать одномерную задачу. В этом случае температура является функцией только одной переменной Т(х). Подставив этот результат в закон Дарси, получаем

¿к =--(Т VC + С VT).

гмк

Общее уравнение переноса примеси с учетом дополнительного конвективного транспорта газа как сплошной среды [8] записывается в виде

дС -

— + (VV)C = DАС . дг

Здесь V - макроскопическая скорость физически малого элемента газа, Б - коэффициент диффузии. Выразим в уравнении для переноса примеси сред-немассовую скорость:

V = -

Ра +Рк Ра крк Ы

гмк

(Т VC + С VT),

так как рк /ра ^ 1. Окончательно получаем уравнение вида

дС кркЯ дг гмк

(т(VC)2 + CVT •VC) = DАС . (15)

В отличие от работы [2] в уравнении (15) присутствует дополнительный член пропорциональный градиенту температуры.

В уравнение (15) входят такие нетривиальные параметры для разреженной газовой среды, как коэффициент диффузии С и коэффициент динамической вязкости Эти параметры требуют дополнительного обсуждения и могут быть вычислены с высокой степенью достоверности с помощью, в некотором смысле, точных формул молекулярно-кинетической теории газов [9]:

с = 3-

кТ \жкТ = 3 (кТ )

8 сг12 р^ 2/

32

-П2

; р^

жш0

г =

1' 2 3

3/2

(т0 кТ )

12

й

■Лг

Здесь 012 - эффективное сечение рассеяния для двух частиц, / - приведенная масса. Для двух примерно одинаковых по массе и размерам частиц имеем 012 = тгй2, /2 = те. Масса одного атома кремния равна т0 = //Ыа = 0.028 кг/Моль/ 6.02-1023 1/Моль = 4.7-10-26 кг. Из табличных данных известно, что диаметр атома кремния равен йя, = 0.23-10-9 м. Таким образом, имеются все воз-

можности для аккуратной оценки значений коэффициентов диффузии и динамической вязкости.

3. Программная реализация

При численном решении одномерной задачи с помощью метода конечных разностей была использована явная схема решения дифференциального уравнения в частных производных. Для реализации устойчивости схемы пространственные производные вычислялись как односторонние разности «назад». Число узлов в сетке принималось равным N = 85. Компьютерный код был реализован на языке программирования FORTRAN-90. В соответствии с явной схемой данные на следующем шаге по времени находятся непосредственно из данных на предыдущем шаге.

Материальные параметры среды, использованные в ходе вычислений, имеют следующие значения: проницаемость вычислялась по формуле (13) и оказалась равной к = 2-10-3 м2, динамическая вязкость п = 3.5-10-4 Пас, коэффициент диффузии Б = 0.7 м2/с, молярные массы кремния и аргона равны соответственно ]ик = 0.028 кг/моль, = 0.04 кг/моль.

Расстояние между изделием и поверхностью расплава менялось в реалистичных пределах. Ниже будут представлены результаты расчетов для Ь = 0.9 м, и 1.6 м. Температуры расплава Т2 и образца Т соответственно менялись в интервале 1720 -г-1770 К. Рассматривались два возможных варианта: Т2 > Т и Т2 < Т1. Предполагалось, что распределение температуры зависит от координаты по линейному закону.

4. Результаты расчетов

Как показывают расчеты, система выходит на установившийся режим фильтрации за очень короткое время. В результате вне зависимости от разности температур в расчетах время выхода примерно составляло ге = 1 с. Столь малое время установления объясняется высокой степенью разреженности среды. Как уже неоднократно указывалось выше, обсуждаемый технологический процесс происходит в условиях среднего вакуума. В ходе вычислений были получены распределения концентраций кремния для двух расстояний между поверхностью изделия и зеркалом расплава, а также величина потока кремния и его среднее значение на левой границе в зависимости от времени.

На рис. 3 и 4 представлены результаты расчета концентрации паров кремния в пространстве между образцом и зеркалом расплава. Зависимость, изображенная прямой штрихпунктирной линией 1, представляет собой стационарное решение при учете только диффузионного механизма переноса.

Р^а +РА РЛ

2

0.8-С

0.6-

0,4-

0,2-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

У* ¿2 / / / /

// / / / / ' 1

/ /

/ / /

20

40

60

Лг 100

0.8-С

0.6-

0,4-

0,2'

- ^......" / 7

? / /

// / / / / / ' 1

/

/

/

/ г

0 20 40 60

N 100

Рис. 3. Концентрация паров кремния между изделием и зеркалом расплава при L = 1.6 м; штрихпунктирная линия 1 соответствует расчетам с помощью классического уравнения диффузии [8] без учета конвективного переноса, семейство кривых 2 отвечает модифицированному уравнению (15), сплошная линия — Т1 = Т2, штриховая с короткой чертой — Т2 > Т1, штриховая с длинной чертой — Т1 > Т2

Семейство кривых 2 получено с учетом нелинейных слагаемых в уравнении (15). Видно, что при учете дополнительного макроскопического механизма переноса кремний начинает заполнять все пространство внутри реторты и его концентрация резко убывает только вблизи изделия.

Это приводит к возникновению достаточно высокого градиента концентрации кремния на левой границе, что обеспечивает наблюдаемый приток кремния в образец. Сплошная линия в семействе кривых 2 описывает распределение паров кремния в отсутствие градиента температуры. Кривая с коротким штрихом лежит выше сплошной линии и соответствует условию Т2 > 71, что вполне закономерно, т.к. более высокая температура расплава обеспечивает дополнительное приращение давления на правой границе. В противоположном случае, когда образец разогревается до более высокой температуры нежели расплав, наблюдается замедление массопереноса кремния от зеркала расплава к изделию.

Заметим, что относительные изменения концентрации в задаче велики, т.к. на правой границе за счет испарения имеет место ее определенное значение, а на левой границе предполагается полное поглощение паров кремния пористой средой. В то же время относительная разница температур в расчетах не может быть большой, т.к. она отражает реальные условия технологического процесса. В расчетах максимальное значение ДТ/Т составляло

Рис. 4. Концентрация паров кремния между изделием и зеркалом расплава при L = 0.9 м; штрихпунктирная линия 1 соответствует расчетам с помощью классического уравнения диффузии [8] без учета конвективного переноса, семейство кривых 2 отвечает модифицированному уравнению (15), сплошная линия — Т1 = Т2, штриховая с короткой чертой — Т2 > Т1, штриховая с длинной чертой — Т1 > Т2

примерно 0.03, поэтому эффект неоднородности температурного поля в рабочем пространстве печи не оказывает существенного влияния на распределение паров кремния. При увеличении разности температур этот эффект закономерно должен стать сильнее. Также при рассмотрении семейств кривых 2 на рис. 3 и 4 можно увидеть, что включение разности температур не изменяет качественно форму решения. В работе [2] для распределения концентрации в стационарном одномерном случае аналитически был получен логарифмический профиль. При включении разности температур на левой и правой границах расчетной области в задаче не удается получить аналитическое решение. Однако из численных расчетов видно, что профиль концентрации качественно не меняется.

Следует отметить, что в связи с добавлением фактора неоднородности температуры возникает необходимость немного изменить подход к определению управляющих параметров в задаче.

Как было показано в [2], при постоянной температуре и плотности кремния р = 1.87-10-5 кг/м3 стационарное распределение концентрации, а именно, крутизна графика, характеризуется всего одним знакоположительным управляющим параметром

¥ =

У/к ° кЯТ рк

0,8-

С

0,60,40,2-

1 ^

/ «

ч//

0 +

0

0,4

0,8

1,2 х [м] 1,6

10-

дх 8-

- 2

-

1 1

1 -1- 1 -,- 1

0,2

0,4

0,6

0,8 I [с] 1

Рис. 5. Семейства кривых концентрации паров кремния между изделием и зеркалом расплава для разных L: 1 - L = 0.9 м, 2 - 1.6 м, сплошная линия - Т1 = Т2, штриховая (с короткой чертой) - Т2 > Т1, штриховая (с длинной чертой) - Т1 > Т2

Параметр у является безразмерным и представляет собой отношение коэффициента диффузии к блоку параметров, которые в совокупности отвечают за формирование градиента давления в реторте вследствие испарения кремния. Таким образом, он отвечает за вклад конвективного слагаемого в общий транспорт паров кремния. Для приведенных выше значений материальных параметров среды имеем у ~ 0.012. Это означает, что конвективное слагаемое вносит существенный вклад в общий массоперенос кремния.

Теперь температура меняется от точки к точке, и механизмы массопереноса развязаны. Однако размерный комплекс

Г = = 0.32,

ЛМк

стоящий перед скобкой в левой части уравнения (15), по-прежнему отвечает за величину суммарного конвективного массопереноса в системе, и он не является малым параметром в задаче по сравнению с еще одним параметром в правой части уравнения (15) Б = 0.7 м2/с.

Помимо сравнения решений модифицированного уравнения (15) и обычного уравнения диффузии имеет смысл наложить друг на друга профили концентраций в реальном масштабе, полученные для разных расстояний между образцом и зеркалом расплава. Два семейства кривых для L = 0.9 м и 1.6 м представлены рис. 5. В результате можно видеть, что с ростом расстояния L массоперенос в системе происходит менее интенсивно, что визуально определяется меньшей крутизной графиков вблизи левой границы.

Рис. 6. Значение производной от концентрации кремния (линия 1), вычисленное исходя из простейшего уравнения диффузии, и аналогичная производная (кривая 2) при решении уравнения (15) в зависимости от времени

На рис. 6 изображен график производной концентрации кремния по поперечной координате х на границе изделия в зависимости от времени (кривая 2). Сначала эта производная имеет нулевое значение (при t < 0.08 с кремний еще не успевает дойти до левой границы), далее можно наблюдать ее резкий рост и быстрое приближение к асимптотическому значению.

Поток кремния в образец в каждый момент времени определяется тангенсом угла наклона графика концентрации на левой границе расчетной области. В свою очередь, значение производной, которое определяется из простейшего уравнения диффузии (линия 1), конечно, тоже как-то меняется с течением времени, но после ее выхода на асимптотическое значение (сплошной участок) данное значение заметно меньше, чем при учете дополнительного конвективного переноса паров кремния. Таким образом, расчеты показывают, что при добавлении конвективного массопереноса имеет место превышение потока кремния в образец больше чем на порядок по сравнению с ситуацией, когда рассматривается только диффузионный механизм переноса.

5. Заключение

Выполнено теоретическое исследование влияния неоднородности температурного поля на перенос газообразного кремния от зеркала расплава к изделию при реализации технологического процесса высокотемпературного силицирования в условиях среднего вакуума. Показано, что реалистичные значения перепадов температуры между зеркалом расплава и изделием незначительно вли-

2

6

4

2

0

0

яют на интенсивность переноса. Тем не менее, установлено, что превышение температуры расплава над температурой изделия усиливает массо-перенос. В противоположном случае имеет место ослабление потока кремния к изделию.

Список литературы

1. Демин В. А., Марышев Б. С., Меньшиков А. И. Численное моделирование процесса высокотемпературного силицирования углерод-углеродного пористого материала // Вестник Пермского университета. Физика. 2021. № 3. С. 56-62. DOI: 10.17072/1994-3598-2021-3-56-62

2. Демин В. А., Демина Т. В., Марышев Б. С. Физико-математическая модель переноса газообразного кремния в ходе высокотемпературного силицирования углеродных композитных материалов // Вестник Пермского университета. Физика. 2022. № 3. С. 48-55. DOI: 10.17072/1994-3598-2022-3-48-55

3. NieldD. A., Bejan A. Convection in Porous Media. New York: Springer, 2006. 654 p.

4. Лейбензон Л. С. Движение природных жидкостей и газов в пористой среде. М.: Гостехиздат. 1947. 244 с.

5. Баренблатт Г. И., Ентов В. М., Рыжик В. М. Движение жидкостей и газов в природных пластах. М.: Недра, 1984. 211 с.

6. Ландау Л. Д., Лифшиц Е. М. Механика сплошных сред. М.: Наука, 1953. 788 с.

7. Леонтьев Н. Е. Основы теории фильтрации: учебное пособие. М.: МАКС Пресс, 2017. 88 с.

8. Ландау Л. Д., Лифшиц Е. М. Курс теоретической физики. Т. 6. Гидродинамика. М.: Физма-тлит, 2001. 736 с.

9. Гиршфельдер Дж., Кертисс Ч., Берд Р. Молекулярная теория газов и жидкостей. М.: Изд-во ин. лит-ры. 1961. 929 с.

References

1. Demin V. A., Maryshev B. S., Menshikov A. I. Numerical modelling of a high temperature silicifi-cation process in the carbon-carbon porous medium. Bulletin of Perm University. Physics, 2021, no. 3, pp. 56-62. DOI: 10.17072/1994-3598-20213-56-62.

2. Demin V. A., Demina T. V., Maryshev B. S. Physical and mathematical model of gaseous silicon transfer during high-temperature siliconization of carbon composite materials. Bulletin of Perm University. Physics, 2022, no. 3, pp. 48-55. DOI: 10.17072/1994-3598-2022-3-48-55.

3. Nield D. A., Bejan A. Convection in Porous Media. New York: Springer, 2006. 654 p.

4. Leibenzon L. S. Dvizhenie prirodnykh zhidkostei i gazov v poristoi srede [The Movement of Natural Fluids and Gases in a Porous Medium]. Moscow: Gostekhizdat, 1947. 244 p. (In Russian).

5. Barenblatt G. I., Yentov V. M., Ryzhik V. M. Dvizhenie zhidkostei i gazov v prirodnykh plastakh [The Movement of Liquids and Gases in Natural Reservoirs]. Moscow: Nedra, 1984. 211 p. (In Russian)

6. Landau L. D., Lifshitz E. M. Mekhanika splosh-nykh sred [Continuum Mechanics]. Moscow: Nau-ka, 1953. 788 p. (In Russian)

7. Leontiev N. E. Osnovy teorii fil'tratsii [Fundamentals of the Theory of Filtration]: a study guide. Moscow: MAKS Press, 2017. 88 p. (In Russian)

8. Landau L. D., Lifshitz E. M. Fluid Mechanics: Course of Theoretical Physics, Vol. 6. Oxford: Butterworth-Heinemann, 1987. 554 p.

9. Hirschfelder J. O., Curtiss C. F., Bird R. B. Molecular Theory of Gases and Liquids. New York: Wiley, 1954. 1219 p.

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом: Демин В. А., Зинурова В. Э. Влияние неоднородности температурного поля на перенос паров кремния от зеркала расплава к изделию в ходе высокотемпературного силицирования // Вестник Пермского университета. Физика. 2023. № 2. С. 51-58. DOI: 10.17072/1994-3598-2023-2-51-58

Please cite this article in English as:

Demin V. A., Zinurova V. E. Influence of the temperature field inhomogeneity on the transfer of silicon vapor from the melt mirror to the product during high-temperature siliconization. Bulletin of Perm University. Physics, 2023, no. 2, pp. 51-58. doi: 10.17072/1994-3598-2023-2-51-58

Сведения об авторах

1. Виталий Анатольевич Демин, доктор физ.-мат. наук, профессор, заведующий кафедрой теоретической физики, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614990

2. Вероника Эдуардовна Зинурова, студентка физического факультета, Пермский государственный национальный исследовательский университет, ул. Букирева, 15, Пермь, 614990

Author information

1. Vitaly A. Demin, Doctor of Physical and Mathematical Sciences, Professor, Head of the Department of Theoretical Physics, Perm State University; 15, Bukireva st., Perm, 614990, Russia

2. Veronika E. Zinurova, Student of the Faculty of Physics, Perm State University; 15, Bukireva st., Perm, 614990, Russia

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.