УДК 538.953:538.975
И. Н. Алиев, С. В. Резник, С. О. Юрченко
О ФРАКТОННОЙ МОДЕЛИ ТЕПЛОВЫХ СВОЙСТВ НАНОСТРУКТУР
Рассмотрена модель теплоемкости самоподобных структур на основе фрактонного представления локальных движений групп атомов. Получены определяющие соотношения для свободной энергии Гельмгольца, тепловой энергии и теплоемкости системы в зависимости от температуры, аналогичные соотношениям теории теплоемкости кристаллов Дебая. Анализируются предельные случаи полученных соотношений, а также условия теплового равновесия однородных наноструктур фрактальной размерности (с нулевой связностью). Найдена единственная стабильная структура для однородного фрактала с нулевой связностью.
Исследованию теплофизических свойств материалов в структурных состояниях с дробной размерностью уделяется значительно меньшее внимание, чем изучению фрактальной механики деформирования и разрушения [1-3]. Возможно это связано с тем, что в обычных условиях и температурных режимах эксплуатации тепловые свойства материалов изменяются относительно слабо. Однако трудно переоценить важность теоретических моделей и анализа свойств и структурных превращений материалов, которые неизбежны уже с наномасштабного уровня для широкого круга применений пористых и нанокристалли-ческих материалов.
Примером конструкций, для которых важно прогнозирование тепловых и структурных свойств материала, являются теплозащитные конструкции аэрокосмических систем. В последнее время в связи с их разработкой возрос интерес к углерод-керамическим композиционным материалам (УККМ) [4]. Вместе с тем для некоторых карбидов были обнаружены фрактальные модификации структуры [5, б], которые приведены ниже:
Исходный материал............ poly-SiC poly-TiC poly-Mo2C 6H-SiC
Фрактальная размерность D.... 1,97-2,97 1,64-2,57 2,36-2,70 2,34-2,71 Настоящая работа посвящена анализу тепловых свойств кристаллических тел фрактальной структуры с учетом различной размерности и связности. Физическая сущность теплоемкости предполагается аналогичной фононной, а электронная составляющая при этом не рассматривается.
Решение задачи о теплоемкости фрактального кристалла, приведенное в работе [7], по ряду причин можно считать некорректным. Во-первых, рассмотрены фрактальные структуры только с нулевой связностью, что справедливо лишь для достаточно плотных структур; во-вторых, при исследовании условий стационарности состояния
структуры не учитывается изменение характеристическом температуры Дебая, что приводит впоследствии к ошибочному выводу о стационарности структуры с размерностью 2,18.
В самом общем случае можно рассмотреть нанокластер фрактальной размерности, а затем перейти к относительно непрерывным и достаточно протяженным структурам. В предположении, что кластер можно считать изотропным, в теле возможны собственные локальные движения, описываемые мнимыми частицами — фрактонами [1,8]. Частота фрактонов связана с абсолютной величиной волнового вектора к нелинейным дисперсионным соотношением
мо- 2 + в . ( £ 4'/2
ш = Лко, а =-, Л и ио —
2 ' 0 V 2п/
где и0 — скорость низкочастотных (фононных колебаний); £ — размер нанокластера; в — индекс связности фрактальной структуры кластера. Примеры структур с различной связностью приведены на рис. 1.
Число собственных состояний в спектре волн с абсолютной величиной волнового вектора от к до к + д,к есть
2 В В-о
4пк з Б ш о
V-г д,к = V—--д,ш. (1)
(2п)3 2п2а л В
Последние выводы могут быть получены и из рассуждений, связанных с масштабной инвариантностью фрактальных структур и процессов фрактонных колебаний в них. Характерно, что атомы фрактального нанокластера могут иметь менее трех степеней свободы. Размерность же фрактального кластера остается постоянной, по крайней мере, в течение малого промежутка времени. Составляющие части фрактала движутся во фрактальном пространстве, не нарушая геометрии кластера, т.е. не пересекаясь друг с другом. При высоких температурах возбуждены все БИа^ колебаний, поэтому целесообразо построить интерполяционную теорию, в которой на всем протяжении спектра колебаний частоты распределены по закону (1), который в действительности справедлив для малых частот. Спектр начинается с частоты
Рис. 1. Вид структур с разной фрактальной связностью [8]
ш = 0 и обрывается на некоторой частоте шт; последняя может быть определена условием равенства полного числа колебаний предельному значению БИа^, откуда
а
Шт
= Л
Таким образом, распределение частот в данной модели задается формулой
D-a
■>2,
, , NaD2 РШ a
g (Ш) = ----—.
О ^
Шт
Свободная энергия системы определяется соотношением [9]
F = Neo + kT ^ ln (1 - e-ñw/kT) =
а 2 Wm
= Neo + kTNDV I шln (1 - e-nw/kT) dш. ошт o
Можно ввести обобщенную фрактонную дебаевскую температуру тела в, тогда
T \ а
©
T
F = Neo + NaD2vkT^-j Jz°-1 ln(1 - e-z)dz, ^
" o (2)
Ншт D
в = ——, a = —. k a
Вводя обобщенную функцию Дебая порядка a и интегрируя (2) по частям, для энергии Гельмгольца получаем выражение
F = Neo + vRTAD lní 1 - e-© j - dJ T i />
x
n . . D í za ,
Da(x) = — --dz.
aw xa J ez - 1
o
dF д (F\
Тогда для энергии E = F = T^t = -T2 dT\T) имеем
E = Neo + vDRTD^T^j , где R — универсальная газовая постоянная. Отсюда для теплоемкости
Рис. 2. Температурная зависимость теплоемкости фрактальных кластеров различной размерности при нулевой связности (в = 0):
1, 2, 3 — Б = 1; 2; 3 соответственно
справедлив ряд соотношений:
C = vDRt
Da (x) =
xa I
ez — 1
dz;
(3)
2 D
а =
2 + в'
На рис.2 показан вид зависимостей (3) теплоемкости от безразмерной температуры при параметре связности в = 0.
Уменьшение размерности системы сопровождается уменьшением ее теплоемкости при повышенных температурах, что, однако, становится неверным для некоторых структур со связностью в > 1. Последнее обстоятельство проиллюстрировано на рис.3, где фрактал с размерностью 2,5 обладает большей теплоемкостью, чем таковая для 3Б-кристалла.
При пониженных температурах теплоемкость фракталов с меньшей размерностью становится выше, чем у фракталов с большей размерностью. Увеличение связности должно приводить к увеличению теплоемкости. Вместе с тем, изменение фрактальной размерности и связности приводит к изменению обобщенной температуры Дебая. Эволюция фрактальных кластеров может осуществляться двумя путями — либо изменением фрактальной размерности системы, либо изменением связности фрактального ансамбля.
В случае однородной фрактальной структуры с нулевой связностью функция свободной энергии принимает вид
x
а
z
Рис. 3. Зависимость теплоемкости фрактальной структуры (О = 2,5) от температуры в случае различной связности:
1,2,3 — 0 = — 1, 0, 1 соответственно
D T
F = Nsq + vRTD2 (T) j zD-1 ln (1 - e-z) dz =
©
D-1
0 1
= Neo + vRTD2 J zD-1 ln ^i - e-dz; o
@ = Пщ ^ 2n2 Nav D
Обобщенная температура Дебая выражается через обычную характеристическую температуру Дебая для 3Б-кристаллов как
1 3-D 1
2 лт , A q / \ п / ть\ п 3
f?(3 - е-Ш D (f)D e<D
и зависит от структуры, так как связана с размерностью; при скорости звука и0 ~ 103 м/с имеем
1
8(3-Р) Л 3-Л
з"
е-тео, TT - 10 d — е0D
Характерный вид функции Т приведен на рис.4. Даже малое уменьшение фрактальной размерности структуры приводит к существенному увеличению обобщенной температуры Дебая.
При характерных значениях фрактальной размерности для кластеров карбидных материалов обобщенная температура Дебая может отличаться от обычной характеристической температуры Дебая в несколько десятков раз.
Рис. 4. Вид функции безразмерной обобщенной температуры Дебая в зависимости от фрактальной размерности (в = 103 ф
Таким образом, температуры порядка тысяч градусов для наноструктур с фрактальной размерностью D < 3 оказываются "низкими". В этом случае верхний предел в выражении для свободной энергии можно заменить на бесконечность, откуда
п
D
F = Neo + vRTD^-J Г (D + 1) Z (D + 1).
где r(z) — гамма-функция Эйлера; Z(z) - дзета-функция Римана.
В состоянии термодинамического равновесия тепловая энергия имеет максимальное значение
Е ^ max.
Подстановкой оценочных выражений для обобщенной температуры Дебая легко получить характерный вид зависимости тепловой энергии однородного фрактально-структурного материала от его фрактальной размерности, который показан на рис. 5.
В работе [7] рассматривалась аналогичная задача о теплоемкости фрактальной структуры и ее состояниях равновесия. Однако при
Рис. 5. Зависимость тепловой энергии от фрактальной размерности для однородной структуры
выводе условий равновесия не учитывалось изменение обобщенной температуры Дебая при уменьшении фрактальной размерности, которое ощутимо сказывается на тепловой энергии системы. Стабильность структур с размерностью меньше 3 при этом оказывается невозможной.
Проведенный анализ особенностей тепловых свойств фрактальных структур кристаллических тел позволяет сделать вывод о том, что единственной термодинамически устойчивой структурой, находящейся в тепловом равновесии при характерных значениях параметров для твердых тел, в случае однородного фрактального кластера является 3Б-кристаллическая решетка. По всей видимости, устойчивые структуры с размерностью, отличной от 3, следует искать в случае ненулевой связности.
Отдельные результаты работы получены при финансовой поддержке РФФИ по гранту РФФИ 06-08-01516а.
СПИСОК ЛИТЕРАТУРЫ
1. Зосимов В.В., Лямшев Л. М. Фракталы в волновых процессах // УФН. 1995.-№ 4.-С. 361-402.
2. Бобылев С. В., М о р о з о в Н. Ф., Овидько И. А. Испускание дислокации порами в нанокристаллических металлах // Физика твердого тела. -2007. - № 6. - С. 1044-1049.
3. Веттегрень В. И., Башкарев А. Я., Б а р у с о в А. В. и др. Температурная зависимость прочности углеродного волокна и трехмерно армированного углерод-углеродного композита // Журнал технической физики. - 2008. -№ 1.-С. 63-67.
4. Материалы и покрытия в экстремальных условиях. Взгляд в будущее: В 3 т. / Под ред. С В. Резника. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2002.
5. Сморгонская Э. А., Кютт Р. Н., Гордеев С. К. О фрактальном характере структуры нанопористого углерода, полученного из карбидных материалов // ФТТ. - 2000. - № 6. - С. 1141-1146.
6. Ч м е л ь А. Е., С е м е н о в А. Д., С м и р н о в А. Н. Проявление в рама-новском спектре фрактальной геометрии трещины в стекле // ФТТ. - 1999. -№ 6. - С. 1030-1034.
7. Рехвиашвили С. Ш. К вопросу о теплоемкости нанокристаллических веществ // ЖТФ. - 2004. - № 11. - С. 65-69.
8. З е л е н ы й Л. М., М и л о в а н о в А. В. Фрактальная топология и странная кинетика // УФН. - 2004. - № 8. - С. 809-852.
9. Л а н д а у Л. Д., Л и ф ш и ц Е. М. Статистическая физика. - М.: Физматлит,
2002.- 616 с. ООШОПП7
Статья поступила в редакцию 29.10.2007
Исмаил Новрузович Алиев родился в 1945 г., окончил в 1969 г. Московский инженерно-физический институт (МИФИ). Д-р физ.-мат. наук, профессор кафедры "Физика" МГТУ им. Н.Э. Баумана. Автор более 60 научных работ в различных областях физики.
I.N. Aliev (b. 1945) graduated from the Moscow Engineering and Physical Institute in 1969. D. Sc. (Phys.-Math.), professor of "Physics" department of the Bauman Moscow State Technical University. Author of over 60 publications in the various fields of physics.
Сергей Васильевич Резник родился в 1947 г., окончил МВТУ им. Н.Э. Баумана в 1971 г. Д-р техн. наук, профессор, заместитель зав. кафедрой "Ракетно-космические композитные конструкции" МГТУ им. Н.Э. Баумана. Автор более 200 научных работ в области теплового проектирования, производства и испытания объектов ракетно-космической техники.
S.V. Reznik (b. 1947) graduated from the Bauman Moscow Higher Technical School in 1971. D. Sc. (Eng.), professor, deputy head of "Rocket and Space Composite Structures" department of the Bauman Moscow State Technical University. Author of more than 200 publications in the field of thermal design, manufacturing and testing of objects of rocket and space technology.
Станислав Олегович Юрченко родился в 1985 г. Бакалавр техники и технологии по направлению "Авиа- и ракетостроение", магистрант кафедры "Ракетно-космические композитные конструкции" МГТУ им. Н.Э. Баумана.
S.O. Yurchenko (b. 1985) - bachelor in engineering and technology for areas of aircraft and rocket building, undergraduate of of "Rocket and Space Composite Constructions" of the Bauman Moscow State Technical University.
В издательстве МГТУ им. Н.Э. Баумана в 2008 г. вышла в свет книга
Бром А.Е., Колобов А.А., Омельченко И.Н.
Интегрированная логистическая поддержка жизненного цикла наукоемкой продукции: Учебник / А.Е. Бром, А.А. Колобов, И.Н. Омельченко; Под ред. А.А. Колобова. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. - 296 с.
Изложены основы концепции логистической поддержки жизненного цикла наукоемкой продукции, обеспечивающей эксплуатационную надежность и качество послепродажного этапа обслуживания изделий, т.е. ключевые факторы конкурентоспособности современной сложной техники. Рассмотрены основы теории надежности и методы продления эксплуатационного ресурса сложных технических объектов, методы математической формализации потоковых процессов в логистических системах и функционального моделирования с использованием стандарта IDEF0. Описаны основные элементы структуры системы интегрированнной логистической поддержки, приведены образцы задач логистического анализа, проанализированы функции системы материально-технического обеспечения, стратегии эксплуатации и формы технического обслуживания и ремонта сложной техники.
Содержание учебника соответствует курсу лекций, который авторы читают в МГТУ им. Н.Э. Баумана.
Для студентов, аспирантов и преподавателей технических и экономических специальностей высших технических учебных заведений, а также для специалистов и руководителей промышленных наукоемких предприятий, научно-производственных объединений.
По вопросам приобретения обращаться по тел. (499) 263-60-45; e-mail: press@bmstu.ru