НАУЧНОЕ ИЗДАНИЕ МЕТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. КБМ 1994-0408
электронный научно-технический журнал
Метод расчета нестационарного нагрева наноструктур
# 09, сентябрь 2013
DOI: 10.7463/0913.0617255
Воробьев Д. А., Хвесюк В. И.
УДК 536.2.022
Россия, МГТУ им. Н.Э. Баумана den vorobyow@mail.ru khves@power.bmstu.ru
Введение
Проблема расчета теплопереноса в многослойных наноструктурах связана с развитием ряда направлений нанотехнологий, где эти структуры нашли применение, в частности, микроэлектроники, основанной на структурированных полупроводниковых микроструктурах.
Одним из видов многослойных упорядоченных структур являются бинарные структуры - плоские образования с толщиной одного слоя 10-50 нм. Количество слоев в структуре лежит в пределах от сотен до тысяч. Слои выполнены из двух различных материалов AlAs и GaAs и чередуются друг с другом так, что у каждого слоя оба соседних выполнены из другого материала.
Толщина слоя в таких системах обычно значительно меньше длины свободного пробега переносчиков теплоты (электронов и фононов). В таком случае закон Фурье не применим для описания распространения теплоты в одном слое. Использование уравнения теплопроводности Фурье-Кирхгофа исключается.
Для описания процесса теплопереноса в слоистых структурах предлагается следующая модель. Внутри каждого слоя температура одинакова по всей его ширине. На границах слоев имею место конечные скачки температуры, впервые экспериментально обнаруженные П.Л. Капицей в 1941 году на границе жидкого гелия и материала сосуда, в котором находился гелий. Таким образом, в
многослойных наноструктурах картина распределения температуры значительно отличается от таковой для массивных сплошных тел.
Перенос теплоты в слоистых наноструктурах определяется не процессами рассеяния электронов и фононов внутри решетки, а процессами рассеяния частиц на интерфейсах (границах слоев). Иными словами, процесс теплопереноса определяется тепловым сопротивлением Рк (обратная величина - тепловая проводимость ) Капицы. Вместо закона Фурье следует использовать формулу, связывающую удельный тепловой поток д с разностью температур соседних слоев (ТГТ2), впервые предложенную И.М. Халатниковым для случая контакта жидкого гелия с твердым телом [1]:
Позднее было доказано, что сопротивление Капицы имеет место между любыми парами различных материалов [5].
Процесс теплопереноса считается одномерным, то есть распространение теплоты рассматривается только в направлении перпендикулярном поверхности слоев.
Модель теплопереноса в слоистой структуре Рассматривается одномерный нестационарный процесс распространения тепла в теплоизолированной упорядоченной многослойной наноструктуре, состоящей из 1000 чередующихся слоев материалов Л1Лб и ОаЛБ. В таблице 1 для справки приведены теплофизические характеристики рассматриваемых материалов (теплоемкость С, плотность р, скорость звука у8, средняя величина сопротивления Капицы при контакте материалов рк).
Таблица 1. Теплофизические характеристики материалов [4].
Т епло физическая характеристика Материал
АІАб ОаЛБ
С, Дж/кгК 424 327
р, кг/м 3730 5320
/с £ 6258 5242
у -, м2К/Вт 0,625-10-9
В начальный момент времени первые 50 слоев имеют температуру 900 К, а все последующие слои - 300 К (имитация импульсного нагрева, [2]). На рисунке 1 представлен тепловой баланс отдельно взятого і-го слоя.
Рис. 1. Тепловой баланс 1-го слоя.
В соответствии с рисунком 1 для расчета переноса теплоты в наноструктуре используется следующее соотношение:
где 1 - номер слоя в многослойной структуре;
АТ, К - изменение температуры слоя за время т;
И, м - толщина слоя.
Левая часть уравнения (2) есть изменение энергии 1-го слоя, а правая -подведенная и отведенная теплоты за время т, соответственно. Данное соотношение позволяет найти распределение температур в наноструктурах при известных значениях сопротивления Капицы, тогда как формула Халатникова используется для определения сопротивления Капицы.
Условие теплоизоляции внешних границ наноструктуры получено обнулением проводимостей на этих границах.
В расчетах использована неявная конечно-разностная схема, позволяющая избежать ограничений на временной шаг имеющих место при использовании явной схемы расчета. Необходимо отметить, что в случае явной схемы расчета максимальный временной шаг имеет физическое ограничение, связанное со скоростью распространения тепловой волны в материалах многослойной наноструктуры. Максимальное значение шага может быть определено как:
л,
В случае одинаковых толщин слоев скорость звука берется для материала с меньшим ее значением.
Численное исследование процесса теплопереноса
Согласно (2) теплоперенос в многослойных структурах зависит от свойств материалов слоев (теплоемкости С), их толщин И и проводимости интерфейсов между слоями. В последнем случае важно отметить, что сопротивление интерфейсов Рк не одинаково при переходе от слоя из материала А1А8 к слою из материала ОаАБ и от ОаАБ к А1Аб. На рисунке 2 видны «малые» и «большие» скачки температуры на границах слоев, обусловленные различной величиной сопротивлений интерфейсов.
354---------------V-
3001-------------------- ^ - —-----------
О 4х1(Гб 8x10”6 1.2хЮ_ Э 1.6x10“5 2х10-5
X, м
Рис. 2. Характерное распределение температуры по слоям наноструктуры при
ЬА1А8=ЬоаА8=50 нм, 1=10 С.
В дальнейших расчетах для удобства используется отношение
проводимостей:
где акА1Аз - проводимость интерфейса при переходе от А1Аб к ОаАв;
- проводимость интерфейса при переходе от ОаАБ к А1Аб;
На рисунке 3 показаны профили температур, полученные при различных временах расчета.
900
780
660
Т.К.
540
420
300
1 о т—I II 8 с
■'ч 1
\У = 0_?с
1: = Ч 0 6С 1=10 5С
8x10 6 1.6x10 5 2.4x10 5 3.2x10 5 4x10'
Рис. 3. Распределение температуры в многослойной структуре при различных
временах расчета ^ c (ЬА1А<=Ь^А<5=50 нм).
На рисунке 4 представлены распределения температур при различной толщине слоев для случая одинаковых толщин слоев AlAs и GaAs, т.е.
Из графиков видно, что с увеличением толщины слоя процесс теплопереноса замедляется. Хотя ранее было сказано, что в многослойных наноструктурах толщины слоев имеют величину 10 - 50 нм, расчет был также произведен при больших толщинах слоя.
Рис. 4. Распределение температуры в многослойной структуре при толщинах слоев от 10 до 70 нм с шагом 10 нм. На рисунке цветом выделены распределения температуры при крайних значениях толщины слоя. 1=10 с.
Исследование теплопереноса при различных толщинах слоев AlAs и GaAs (ЬА1А#Ь5аА<0 показало, что в случае ЬА1А!5=20 нм и ^^^=40 нм, а также в аналогичных случаях, полученное распределение температуры с высокой точностью совпадает с распределением для ЬА1А<5=Ь5аА<5=(20+40)/2=30 нм.
В расчетах также было получено, что соотношение между скоростями теплопереноса сохраняется во времени при различной толщине слоя.
Все выше описанные расчеты были произведены при значении отношения проводимостей к=5. Расчеты, проведенные для значений к в диапазоне от 0,1 до 10, показали незначительное различие температурных профилей при различных значениях к. Данный факт объясняется неизменностью суммарного термического сопротивления интерфейсов структуры, т.е.:
Несмотря на то, что для пары A1As/GaAs термическое сопротивление Капицы является известной величиной, в целях исследования влияния последнего на процесс теплопереноса были произведены расчеты при различных величинах сопротивления Ръ, результаты которых показаны на рисунке 5.
660
570
Т,К
480
390
300
- _9м*К - рь= 1-9 10 У—— Ч\\ К Вг
к = 0.1 10 2 V
_ Р Вт
6x10'
-6
1.2x10 3 1.8x10 3 2.4x10 3
ЗхЮ
Рис. 5. Распределение температуры в многослойной структуре при изменении
термического сопротивления интерфейсов от 0,1*10-9 до 1,9-10-9 м2К/Вт с шагом
10 2
0,2-10 м К/Вт. На рисунке цветом выделены распределения температуры при
7
крайних значениях термического сопротивления. 1=10 с, к=2.
Расчеты, проведенные при к=5 и к=8, выявили, что тенденция, показная на рис. 5 имеет место и в этих случаях.
Проведенные расчеты основаны на идеализированном случае, когда не наблюдается отклонений проводимостей интерфейсов, т.е. все переходы от АІЛб к ОаЛБ имеют одинаковую проводимость акЛіЛ8, и все переходы от ОЛЛб к АІЛб имеют одинаковую проводимость ак0аЛі5. На рисунке 6 показно распределение температуры при наличии отклонения проводимости произвольно взятого интерфейса.
£00
'00
Г к 600
<00
400
ЭОС
\
ч,
1 /
2 //
Л У
1x10"6 2х 10"6 3x10"6 4x10"* 5x10"*
к. к
Рис. 6. Распределения температуры при отклонении проводимости интерфейса:
1 - проводимость в 100 раз меньше номинальной;
2 - проводимость в 10 раз меньше номинальной.
На графиках рис. 6 виден значительный скачок температуры на интерфейсе с низкой проводимостью. Величина скачка закономерно увеличивается при уменьшении проводимости (увеличении сопротивления) интерфейса, причем эта зависимость носит нелинейный характер.
Выводы
В работе рассмотрен нестационарный процесс теплопереноса в теплоизолированной многослойной упорядоченной наноструктуре А1Ав/ОаАв с заданным начальным распределение температуры. Показаны распределения температуры в различные моменты времени. Установлено замедление процесса теплопереноса при увеличении толщины слоя. Аналогичная тенденция выявлена при увеличении отношения проводимостей Капицы и сопротивления интерфейсов.
Список литературы
1. Халатников И.М. Введение в теорию сверхтекучести. М.: Наука, 1965.
С. 128-136.
2. Попович А.Ф., Ральченко В.Г., Власов А.В., Цыганков П.А. Температуропроводность нанокомпозитных многослойных пленок Ti-Al // Тонкие пленки в электронике: материалы XXII Междунар. симп. (Москва, МГТУ им.
Н.Э. Баумана, 9-11 сентября 2009 г.). М.: ЦНИТИ «Техномаш», 2009. С. 469-473.
3. Samvedi V., Tomar V. The role of interface thermal boundary resistance in the overall thermal conductivity of Si-Ge multilayered structures // Nanotechnology. 2009. Vol. 20, no. 36. Art. no. 365701. DOI: 10.1088/0957-4484/20/36/365701
4. Sadao Adashi. Properties of Semiconductors Alloys: Group-IV, III-V and II-VI Semiconductors. John Wiley & Sons, 2009. 400 p. DOI: 10.1002/9780470744383
5. Cahill D.G., Ford W.K., Goodson K.E., Mahan G.D., Majumdar A., Maris H.J., Merlin R., Phillpot S.R. Nanoscale thermal transport // Journal of Applied Physics.
2003. Vol. 93, no. 2. P. 793 (26 pages). http://dx.doi.org/10.1063/1.1524305
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Calculation method for non-stationary heating of nano-structures
# 09, September 2013 DOI: 10.7463/0913.0617255 Vorob’ev D.A., Hvesyuk V.I.
Bauman Moscow State Technical University, 105005, Moscow, Russian Federation
den vorobyow@mail.ru khves@power.bmstu.ru
A non-stationary process of heat transfer in an insulated ordered multilayer nanostructure with specified initial temperature distribution was considered in the article. The structure consists of 1000 alternating layers of AlAs and GaAs. A heat transfer process model based on heat balance of a single layer was described in this paper. This model was used for heat transfer computations in ordered binary multilayer nanostructures for the first time. A process of heat transfer was investigated at different time moments, with different thicknesses of layers and heat resistances of interfaces. A case of deviation of periodicity in the structure was considered.
Publications with keywords: nanostructure, heat transfer, Kapitza resistance, Kapitza conductivity, binary structure
Publications with words: nanostructure, heat transfer, Kapitza resistance, Kapitza conductivity, binary structure
References
1. Khalatnikov I.M. Vvedenie v teoriyu sverkhtekuchesti [Introduction to the theory of superfluidity]. Moscow, Nauka, 1965. S. 128-136.
2. Popovich A.F., Ral'chenko V.G., Vlasov A.V., Tsygankov P.A. Temperaturoprovodnost' nanokompozitnykh mnogosloynykh plenok Ti-Al [Thermal diffusivity of nanocomposite multilayer films Ti-Al]. Tonkieplenki v elektronike: materialy 22Mezhdunar. simp. [Thin Films in Electronics: proc. of the 22nd Intern. Symp.], Moscow, Bauman MSTU, 9-11 September 2009. Moscow, TsNITI “Tekhnomash” Publ., 2009, pp. 469-473.
3. Samvedi V., Tomar V. The role of interface thermal boundary resistance in the overall thermal
conductivity of Si-Ge multilayered structures. Nanotechnology, 2009, vol. 20, no. 36, art. no. 365701. DOI: 10.1088/0957-4484/20/36/365701
4. Sadao Adashi. Properties of Semiconductors Alloys: Group-IV, III-V and II-VI Semiconductors. John Wiley & Sons, 2009. 400 p. DOI: 10.1002/9780470744383
5. Cahill D.G., Ford W.K., Goodson K.E., Mahan G.D., Majumdar A., Maris H.J., Merlin R., Phillpot S.R. Nanoscale thermal transport. Journal of Applied Physics, 2003, vol. 93, no. 2, p. 793 (26 pages). http://dx.doi.org/10.1063/L1524305