Расчет температурного поля и массопереноса углерода при выращивании монокристаллов алмаза в расплаве металлов

С. А. Ивахненко, О. А. Заневский, А. А. Будяк, И. С. Белоусов

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

Основным методом выращивания монокристаллов алмаза на затравке является метод Стронга— Венторфа, разработанный в 1970—1971 гг. [1, 2]. Перепад температуры между источником углерода и затравкой обычно составляет от нескольких единиц до нескольких десятков градусов; величина градиентов температуры при этом изменяется от 2 до 20 град/мм [3, 4]. Перепад температуры между источником углерода и затравкой, величины осевых и радиальных градиентов температуры обычно подбираются экспериментально, что достаточно трудоемко и требует много времени. Основным материалом нагревателя является графит. Резистивные свойства различных марок графитов отличаются незначительно, поэтому возникают проблемы при задании конфигурации нагревательной системы реакционной ячейки, если требуется изменить величину осевого или радиального градиента температуры или их соотношения.

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

Применение дисперсно-композиционных материалов для формирования резистивной системы нагрева ячеек позволяет значительно упростить процесс задания распределения температуры в реакционном объеме. Особенно эффективно использование методов компьютерного моделирования потенциальных и тепловых полей. Эти методы подробно изложены ранее [6, 7]; они позволяют путем поэтапного моделирования изучить распределение температуры в ячейке и выполнить расчет полей концентраций углерода.

Рассмотрим этот подход на примере используемой нами ростовой ячейки для аппарата высокого давления типа тороид (диаметр полости высокого давления составляет 40 мм).

Схема ячейки для выращивания монокристаллов алмаза на затравке представлена на рис. 1. Нагрев ячейки осуществляется электрическим током через токоподводы 7, 75, трубчатый нагреватель 8, конфигурационные и нагревательные диски 3—5, и 11—13. В стационарном режиме в реакционной ячейке устанавливается тепловое поле с перепадом температуры между источником углерода и кристаллом-затравкой 10— 40 °С. Рост алмаза осуществляется на кристалле-затравке путем диффузии углерода через слой металла-растворителя. Поскольку растворимость углерода прямо пропорциональна температуре, то скорость

Рис. 1. Электрическая и тепловая схемы реакционной ячейки (1/2 часть осевого сечения):

1, 15 — токоподводы; 2, 14 — электрофокусы; 3, 5, 11, 13 — теплоразводящие диски; 4, 12 — верхний и нижний нагревательные диски; 7, 10 — электро- и теплоизоляционные втулки; 8 — трубчатый нагреватель; 6 — источник углерода; 9 — металл-растворитель

роста монокристалла алмаза зависит от величины вышеуказанного перепада температуры.

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

где Цх, у, z) — в зависимости от типа решаемой задачи коэффициент электропроводности, или теплопроводности, или диффузии; U (х, у, z) — потенциальная функция (электрическое напряжение, температура или концентрация углерода); W(x, у, z) — удельная мощность источников электрического тока, тепла или углерода.

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

Представление решения уравнений электропроводности, теплопроводности или уравнения диффузии в виде ряда или интеграла практически невозможно ввиду сложности конфигурации ростовой ячейки и задания граничных условий, а также большого количества составляющих элементов. Для расчета теплового поля ростовой ячейки применялся метод конечных элементов в виде метода конечных разностей [7]. Если затравка будет помещена в центре подложки, то диффузионную задачу можно решать тоже как двумерную (осесимметричную); при расположении еще одной затравки на периферии подложки, диффузионную задачу уже нужно решать как трехмерную; в этом случае количество уравнений возрастает на порядок. Для решения системы уравнений нами был использован метод Гаусса—Зейделя [8].

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

Расчет температурного поля в ячейке роста и распределения температуры в сплаве-растворителе позволил решить задачу диффузионного массопереноса углерода; при этом в качестве граничных условий брали значения растворимости углерода на нижней и верхней поверхностях металла-растворителя при температурах, определенных на предыдущем этапе решения задачи. Плотность диффузионного потока / можно рассчитать как:

где D — коэффициент диффузии в расплаве; vMe-c ~~ мольный объем расплава, насыщенного углеродом, по отношению к алмазу; Ха — растворимость алмаза в расплаве металл-углерод; А Т — перепад температуры между источником углерода и затравкой;

— безразмерная концентрация углерода; С — концентрация углерода; Cmin и Стах — минимальное и максимальное значение концентрации углерода, соответственно, п — нормаль к поверхности равной концентрации.

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

На рис. 2 представлено распределение температуры в ростовой ячейке, обусловленное различными соотношениями длин и диаметров резистивных элементов. При температуре 1420 °С на затравке максимальная температура источника углерода составляет 1530 °С, а максимальная температура на затравочной плоскости равна 1482 °С. Варьируя размеры системы нагрева, массоперенос углерода можно изменить в направлении от центра к периферии подложки, на которой располагаются затравочные кристаллы; расчет указанным методом это хорошо демонстрирует (рис. 2, б, в).

Поделитесь этой записью или добавьте в закладки