Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Физика композитов.Термодинамические и диссипативные свойства

.pdf
Скачиваний:
7
Добавлен:
12.11.2023
Размер:
15.34 Mб
Скачать

ионного механизма соответствующее среднее дается формулой (2.76), а для времени хик подобная процедура также элементарна, и мы на­ ходим, что

/ 1 \

д3с,2е?рг5

 

(2.80а)

\ти/

8io*4p?« <£ при

|£”

 

/ Л

= ^ Щ ц Г 1 прк т > е ш .

(2.806)

\ х1*/

648йр\а{си

 

 

Сравнив теперь (2.75) и (2.806), определяем диапазон температур,

при которых время xlp/Mn < Xj:

 

 

QXD< T < T \

 

 

(2.81)

где

 

 

 

 

Г* =6,45Г 8 \* * \

M £ LЛ2 рМ <£

(2.81a)

 

VGi )

V J \e x У

 

 

Для приведенных выше параметров температура Г* оказывается та­

кой: Т* = 4 (g "/G ,)2 (К). Если

g” = 1, a G, =0,02, то

 

1000 К.

Итак, для 0 1О= 2 0 0 К

при температурах

из диапазона

200 К <

< Т < 1000 К будет преобладать магнонный механизм.

 

 

Для случая структуры D + M общая формула для теплопроводности

(2.67) должна быть представлена в том же виде, т.е.

 

 

х ( Г , V ) = (1 - %’ f

+ X (1 -

) х 0, + V (1 -

) х ,о +

1.

 

но параметры здесь таковы:

 

 

 

 

Ы) _ ( II)

 

 

 

 

 

*00 = X *'00

= х,00

 

 

 

 

и,и

где Ход1*дается формулами (см. выражения (2.63а) и (2.66))

Ц (П )

_

Р ас]

т

а

 

*00

-

2 ~

ПРИ 1

 

 

 

 

о7С

1

 

 

(2.82)

 

 

 

 

 

 

 

 

 

\ 3

2

 

 

*оо

-

 

 

^ хо

При

Т < 0£).

\ QD J

у-

 

 

 

 

 

Напомним, что х0 есть время R/cs,

v

-

Y „(««')_ „(11,12)

(21,22)

*01 “

1* *01

“ *01

+ *01

Если бы не было магнитной фазы, то согласно (2.64), (2.65) и (2.66) для

81

температур Т > 0D мы получили бы формулы

(Ш_

4

3

2

7

 

д а | р р 1с;

 

°'

 

TCifi2 IDQ2 (d)

ПР" C°s ls’

. . ( I D

_

2 V a

4 C|6, c , P P ?

при c0s <cu.

x 01

-

r

 

7 ft2

 

^Ol^OD

Для T<QD

V „2 *01) - CST01 X01 “

\ Q D J

Присутствие же магнитной фазы обязывает ввести эффективную теп­ лопроводность согласно газокинетическому выражению

*01U2) = c f \ \ 0pim » { n 0t ) l d T ) ^ k 2dk/6 %2,

(2.83)

о

 

где согласно формуле (2.78) (после ее небольшой модификации) можно записать, что эффективное время есть

_ ( е / ) _

^ 01 k p h - p h ^ 0 \ k p h - m

(2.84)

T01Jt--------------------------

^01 k p h - p h

^ O lk p h - m

 

 

Времена T0]kph_ph и To\kP h -m нам известны (см. формулы (2.56) и (2.75)). Небольшая модификация формулы (2.75) дает

J ____

(g m )2( M * i ) 2^ /;2

(2.85)

 

\впраъJ2ex\ieH

"^ O lkph -m

 

Очевидно, зависимость (2.56) не изменится. Единственное отличие - это дополнительные латинские буквы (ph - ph) у нижних индексов.

Начнем со случая cs < с^. В этом случае обратное время релак­ сации

1

= В,к,

_

n2G l $ lXDT

 

£ ,=

( 2.86)

^01 k p h - p h

 

 

10368р'4р1а3а16с15

 

 

 

Тогда коэффициент теплопроводности XQ2) с учетом формул (2.84), (2.85) и (2.86) есть

 

 

 

 

 

*0

x2dx

х,(П.12)

i c W j L v - ' H

 

д Т J

 

(2.87)

01

 

^ ( r f t f O i . A f ,) 2 д Т ‘ J

(ех - \ )(* + #•)’

где х0 = nhcJaT, а

 

 

 

 

 

Кз

Л \2

\2

JleXCs

г =

0 1 Р

 

'01

 

 

 

 

 

648 к *01 )

МеМ 0

т

PPi^fq5,

82

Оценим г. Для параметров GQI порядка g0], 0Ою = IK, cs = cls, \ieH = 10К, T = 300 К, /] СХ = 1000 К = 10-13 эрг находим, что г = 10-4. Благодаря этой малости параметра г мы имеем право пренебречь им по сравнению с х (последующие оценки это только подтверждают). Итак, для высоких температур Т > 0О мы имеем право разложить подын­ тегральную функцию в ряд по степеням х, ограничиваясь при этом вто­ рым порядком. В самом деле,

J L / f

xdx

 

д

х°

 

xdx

 

 

= —

Т \

(\ + х / 2 )(х + г)

 

дТ

i

(ех —l)(x + r)

дТ

J

 

= х0 - г!п(х0 / г) г х0 = 0D / Г.

 

 

 

Значит,

 

 

 

 

 

 

 

 

„(11.12) - - 4рагс Л ^ , Н В п

при

T>QD

(2.88а)

01

 

3Jt(s” )2ft2(neM,)2r

( 11. 12) _

2к2Т \ 01

при Т

< 0D.

 

(2.886)

к01

45П5сг

 

 

 

 

 

 

 

 

 

В случае же, когда скорость звука в основной матрице превышает скорость звука в частицах примесной фазы (с5 > с^), обратное время

релаксации здесь (см. (2.56)) есть

 

1 _

1

У —___

G02,eJ,DT(3+d)

(in =

т

v3 •

,(П) ”

 

 

Хг

х 01*рЛ-рЛ

 

1728itp'Pio‘c?A3a - ‘0-

С01

 

Вполне аналогичные выкладки (правда, с немного другими интег­ ралами) приводят нас к такому промежуточному ответу:

(11,12) _ ЙС

д _ Л °

xdx

*01

6 к2

дТ В2г \J

(ех - 1)(х2 + г*2) ’

где новый параметр

 

r *2 _ I

** Л

108Й3рр,с,2,с У (р„Л/, )2(1- d f

gpi

 

701 У

j L ^ H 0 lw O + d)T2

При GQI порядка ggi* -Лех = 1000 К = 10-13 эрг, 0Ою = IK, \1еМх= 0,1 К, Т = 300 К, \LeH = 10 К, р = 3 г/см3 найдем, что по порядку величины г* = Ю-4. Ограничиваясь, так же как и выше, разложением подын­ тегральной функции вплоть до членов вторго порядка по х и пользуясь тем, что параметр г* значительно меньше, чем х0 = 0р/Т, найдем в ре­ зультате простых вычислений, что

х'оУ12’ = 576Р2Р|°6с11^ ( | - ^ ) 3 |п(,+ е / 2Г) при Г > е„.

(2.89)

* G d , 0 2D(3+<i)

Для Т < 0D здесь будет справедлива формула (2.886).

83

Таким образом, окончательно

 

 

V -2

 

х01 - 2

сех01 при Г < 0О,

(2.90а)

ДЛЯ CQJ > сь и при Т > 00

х288p2p,a6C|3,c^fi(i-<J)3

01

%TGlfil,De D(3 + d)

(2.906)

 

для c0s < cls и при Т > 00

 

Х 01 _

4ра3Су/|2«,Ц«Д9р

(2.90в)

З л (^;)2Л2(р,Л/,)2Т '

 

Переходим теперь к вычислению х]0. Для нее имеем

 

 

(11. 12) . (21. 22)

 

х10 “ х10

+ х10

 

Поскольку фаза D немагнитная, то xjo2) = xjg2) = 0. Значит,

 

хю = W

+ xg«

(2.91)

Вычислим эти два коэффициента. Что касается параметра x jo \ то его выражение нам известно (см. соотношение (2.64)). В самом деле,

*(11) -

а Ъг г т 1/4

(2.92)

а CJC1Jt 0

хю -

1927Cfl,6Ql3,4(rf)a3,4( r ) ’

 

Все параметры, фигурирующие в данной зависимости, приведены перед выражением (2.64).

Вычислим x g ” Время релаксации, описывающее этот параметр, связано с механизмом взаимодействия магнона фазы "1" с фононом фа­ зы "0". Для двухмагнонного процесса взаимодействия с фононами (ам­ плитуда (2.696)) интеграл столкновений имеет вид

L { /t } = 2 J t r ' X 1у< 2) I2 {[(1 + Л )/* , (1 + % ) - Л (1 +

Ч:

+Л, К , № « (* ) - е .(* ,)+ е ^ (* 2)]Д (к -к1+ к 2)+

+К1 + Л >Л, % " Л 0 + Л, )(1+ % )]S[em(A) - £„(*,)-

-Ем (*)]Д(к-к, - к 2)+[(1+л)(1+ Л , )ntl - л л , ( 1+

+«tj ) № . (*)+ея № ) - е „ (*2 )]Д(к + к, - кг)}.

Отсюда следует, что время релаксации есть

-(2ГК2)- = 2 ^ " 11 I Vl2) I2 { ((% )-(Л ,)№ /п (Л) - е т ( ^ ) +

Х10*

84

+ е рл (*2 )1 д ( к - к , - k 2 ) + (1+ ( / * , ) + ( n * 2 ))6[Em (k) -

-em(kx) - e ph(k)]A(k - k x- k 2)+((.fkl) - (n*2 ))6[Em (k)+

+Em(kx) - z ph(*2)]A(k+ k, - k2)}.

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

Анализ первых двух законов сохранения приводит к условию, что импульс к2 должен лежать в диапазоне hcJJjex - к2 ^ hcs/JXex + 2к, а соответствующее время релаксации будет тогда (после проведения несложных математических действий, связанных с использованием обычных свойств 6-функции) для первого процесса (новый индекс /) таким:

г(21)(2)/

(8Г)2(ц, м , ) У '

—1] ^ Т

x2exdx

12тсрс^, JXexk

\ hcs j

Л*> (ех -l)[ex+£mWIT - l ]

40*

где пределы интегрирования есть

х, = (2к + hcs / Jlexa})bcs / JXex,

XQ = (-2к + hcs / JXexa2 )hcs / JXex.

Для температур T>0rD/ JXex (правая часть неравенства соответст­ вует примерно 40К) можно приближенно написать, что

1

(g:, )2(n M i)2<"

при ак< 1.

(2.93а)

_(21)(2)/

Л.

40Jк

блрс]а12J \exh ' a

 

 

По порядку величины это время соответствует примерно 10-4- 10"5 с.

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

0 ^ к 2 ^

2 к - hcs/JXex.

 

 

Тогда обратное время этого процесса есть

 

1

 

2 г

 

x2dx

 

(g,")2(n«M,)4e‘”w - i]

 

_(21)(2)//

блрс^! JXexk

\ hcsj

(ех -\)[e~x+Emik)IT - 1]

40*

 

где х* = 2k-hcJJXexa2.

 

 

Если сравнить данное время с временем

(2)/, увидим, что оно

значительно больше. Это же относится и к времени т*о*(2)/// Самое маленькое время из этих трех - это время (2.93а).

Вычислим теперь время релаксации для двухфотонного процесса взаимодействия с одним магноном. Гамильтониан взаимодействия ха­

рактеризуется при этом амплитудой vp,^.

85

+{n(,Zm{k)-hcskx)^dky

После аналогичных выкладок можно получить, что время релакса­

ции для этого процесса будет

 

 

1

 

2,1:

- Ш /й с1 2г1 + /„

\ +

w > “

36T tp V A

l , l '

m( 1 J [1+\ М

+

х ш Г ‘/

3б7СР “

к0

 

 

(2.936)

Пределы интегрирования таковы: ICQ= em(k)/2hcs- к/2, = £m(k)/2hcs + + к/2. Область же изменения волнового вектора к лежит в интервалах

ak*zQD/ 2 Jlex + (Од / 47?„ - \ieH / 7lftC)1/2

0 ’& a k ^ e D/2JUx- ( e l lA J t„ - H ' H I J l'X)ln

 

 

При малых к, таких, что к < £m(k)/hcs, имеем

 

 

1

_ ( M W g ] y

е ^ ) 1

Л m W I T -1]

 

(2.93в)

Д21К1)

576лр2а3с3

%cs

[ е ^ к)11Т - \ \

 

' 10*

 

 

Искомое время релаксации будет, очевидно, определяться суммой

обратных времен (2.93а) и (2.93в). То есть

 

 

1___ { (Ц«М,)‘( а )

\2r\2

 

V

г Л т (к )/Т

-1]

 

 

Ь # )

г(21)

бтсрc3sa2J2exh2a}

576тф2я3с3 he,S /

[ее"*(*)/2Г-1]

с10*

 

 

 

 

 

 

(2.93)

2

для температур Г > 0 D//\ex>Следовательно, теплопроводность после

обезразмеривания всех переменных есть

 

__ 01)

. «С21) _

(11)

+

 

х 10 - х 10

+ х 10

х 10

 

гЗ

( т \ 5/2

 

М 2

dx

2 4 ^

д

 

 

 

 

Зтс2а ,Й 2

д Т

 

о

 

(е * - l ) [ ( 7 / / I e ) 3j c V - lX e * 72 - 1 ) + A i ] '

где время

 

 

 

 

 

 

576тф2Й4с?я3

 

«2

 

960n^2a5pc4

 

т

 

 

° "

 

"(я Г М * ,)2./,4/

 

 

а хй” дается формулой (2.92).

Здесь верхний предел заменен на поскольку Порядок величины параметра XQ колеблется вблизи значений 103—104. Поэтому, если пренебречь членом, пропорциональным х4 в знаменателе интегра­ ла, найдем

___ (ID

/

т \ ъп

 

Г(7 / 2)^(7 / 2),

х,о = х 10

 

 

Зл2а,й2Хо \J\ex )

где Г(7/2) = 15 Vrc/8, £(7/2) s 0,8 и Г(7/2) £(7/2) s 3 Vic/2. 86

И так, окончательно

5, J^2 т

(

т \ ъп

*10 =

{

J lex J

г к ъ ,2 а хН 2 \ 20

аъс с

т1/4

(2.94)

а

С5СЬ

Т 0

6 ^ 3 / 4 ,

, 3 / 4 /

\ 9 2 n a ? Q i " \ < 1 ) а > ' \ Т )

 

Что касается теплопроводности хц, то для нее вполне аналогично имеем

„ _

(11.12)

(21.22)

*

*11 -

*11

+ *11

но xff> - 0 (нет вклада от собственных магнонных взаимодействий - они участвуют лишь в установлении внутренней квазиравновесной тем­ пературы внутри магнонной подсистемы), а потому, комбинируя выра­ жения (2.90а), (2.90в) и (2.94), которые для фазы 1 дают

X i I —

7Г~]

 

ПРИ 7

< 0ш

 

 

V0ID J

а\

 

 

 

 

v

■ Plfllc ls

5 J 2x X

 

r

J

 

 

 

при e 1D < T < J lex,

(2.95)

Л11 —

_

I

T/о

 

 

 

 

87i l T

 

2 n 3 l \ h

l X i

\

J \ex J

 

где

 

 

 

 

 

 

 

 

1

(g,V,M|)2J|4CT

.2

9 6 0 ^ 4

 

T*

576TtpfS4c,7X ’

0

 

Jlx

 

для теплопроводности композита со структурой типа D + М с учетом всех полученных соотношений (2.82), (2.906), (2.94) и (2.95) при ц.«Л < <0iD < 7 < hex находим:

*о+м(7Л ) = ) 4 D+M +$ О -la )*O,O+M + $ 2*ID+M. (2.96a)

где коэффициенты

г.pacs

*D+M~%n2T

 

 

 

 

\ 3 / 2

27Q13Q4C|6JC5P P ^

 

 

15J|3„pg,V;

 

 

 

 

У

^^010010^

1

4У12<)Грр,Я ера3с,

з „ _ l / 4

 

о

**

__. О *

 

 

З я Л ^ м М * ,) 2

192тса16Q3' 4(<i)a3' 4(Г) ’

n

PlfllCb ■

f T

\ 3 / 2

^Л«Р1°1С1д

(2.966)

 

W + M

S n 2 T

J H Q 2D {g ; \ i eM , ) 2 \ 7\ex )

приведены с учетом явных выражении для параметрических зависи-

- * мостеи т0 и А-о-

87

Рис. 2.8. Схематическое поведение коэффициента теплопроводности струк­ тур типа D + М

В частности, формулы (2.96) показывают, что если то кривая будет смещаться влево вместе со всеми экстремальными точками (пунктир­ ная линия). Такая тенденция ясна и качественно тоже понятна, поскольку скорость звука в композите cs есть функция от £*, причем монотонно

убывающая с ростом Температура Т**(Е,*) должна целиком определяться магнитной примесной фазой

Как видно из полученного выражения, теплопроводность линейно растет вместе с магнитным полем (см. последнее слагаемое в R0i). Этот рост, безусловно, ограничен, ибо при увеличении Н время релаксации

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

Из соотношений (2.96) видно, что поведение х при температурах Jex> T > Во можно представить в следующем схематичном виде:

\ 3 / 2

0,

1J \ex )

снекоторыми функциями N^*) и Н-^'). Отсюда очень хорошо видно, что у у-о+м должна наблюдаться тенденция проходить через минимум в

некоторой окрестности температур Т* = (01Z>/iex)2/5ср(4*)» который (ми­ нимум) должен "гулять" влево или вправо в зависимости от концент­ рации.

При совсем высоких температурах (Т > Jex) магнитная структура становится парамагнитной и в соответствии с законом Кюри средняя намагниченность есть линейная функция от отношения HIT. Благодаря этому обстоятельству зависимость х также станет обратно пропорцио­ нальной температуре, но уже в области парамагнетизма и вплоть до температуры плавления.

88

Все сказанное наглядно иллюстрирует рис. 2.8, на котором показа­ но качественное поведение теплопроводности структуры типа D + М в широком интервале температур.

И еще. Обратим внимание на один весьма тонкий момент теории теплопроводности в структурах типа D + М. В нашем случае, когда речь шла о примесной мелкодисперсной магнитной фазе, говорилось о двух подсистемах - о фононной и магнонной. Диапазон температур при этом был выбран не экзотический (сверхнизкие температуры мы не затрагивали), а вполне реальный и экспериментально легко изучаемый. Главное предположение всей изложенной выше теории базировалось на условии, что времена релаксаций внутри магнонной подсистемы и внут­ ри фононной (обозначим их соответственно как время т(т-и,) и т(рл-рл>) являлись наименьшими временами, а взаимодействие между этими двумя подсистемами и определяло некоторую аддитивную часть всей теплопроводности структуры. Для диэлектрической фазы (D) основной матрицы мы говорили только о фононной подсистеме, подразумевая при этом, что установление равновесной температуры есть результат кон­ куренции двух механизмов - четырехфононного и трехфононного взаи­ модействия. Хотя надо добавить, что такие процессы, как пятифо­ нонный, шестифононный и т.д., существуют, но они во внимание не принимались. Если учесть эти механизмы, то окажется, что они вполне конкурируют с механизмом переброса. Действительно, многофононные процессы взаимодействия приводят к степенной зависимости време­ ни релаксации этих процессов от температуры в области Т < 0О, т.е. 1/х(Л'-рЛ) = S/I T 2лм^ где /у _ количество фононов, принимающих участие в

акте рассеяния. Для N = 6 мы видим, что зависимость обратного време­ ни релаксации есть slT11. Но такое поведение т практически неот­ личимо от времени релаксации, связанного с процессами переброса: экс­ поненциальная и степенная зависимости в области низких температур настолько близки, что их можно считать практически одинаковыми. И это если ограничиться лишь шестифононным механизмом. Так что ме­ ханизму Пайерлса (перебросы в область соседних зон Бриллюэна на постоянный вектор Ь) есть вполне достойная конкуренция, забывать о которой не следует.

2.4. СТРУКТУРА M +D

Теплопроводность подобных типов композитов может быть весьма легко описана, опираясь на результаты раздела 2.3. В самом деле, здесь

XAf+D =(1-1r ) 2* * + $ 'a - V ) * o .+ V a - V ) * io + |i*2*ii.

где

„(11.12)

. „(21.22)

»

400 = х00

+ х00

„(П.12)

.„(21.22)

.

Х01 = Х01

+ х01

„(ПЛ2)

+ „(21.22)

*

х10 = х 10

+ х 10

„(П) Х11 = Х11 •

(2.97а)

(2.976)

89

Все эти х мы сейчас и приведем. Действительно, согласно (2.95) с

* *

учетом явного вида параметров т0 и Х,0

(

т \ 3

при Т<TQd,0„,

XQO =

\

9 D J

а

х

- Р E£s_.

1 5 ^ р ас3

х 00 — О +

V^02D(S*neM„)3

 

8л

гт

U D

, ( 2 D

и

Далее, х01 = х01

+ х01

 

5л1т“

т. Л3' 2

Хгк|

 

 

01 “ 2к3,2а П Х 2 \ Jex J

где

т \ 312

при 0D< ^r^y e

па\Ъг г 1Тt 01 / 4

CJC J

192ла6еГ3/4(^*)а*3/4(Г)

(2.98)

(2.99)

** _ 5767cp2fi4c1Ja3

. * * 2 _ 968^дЙ2а15Р|С14

 

х “ / * А Л \2

/4 ’ а

Л0 -

,5 „4

*

(^ о М 'е ^ о )

^ех

 

^ е х ^

 

* _

^ 01еР1р Щ

(d

)

0 *(d*)~

d

+3

d * - c

а — ,

___2 _ 8

2 2

3 ’

i i l v“ ) —

, ,

, * . з ♦

« ~ c Q

s ' c \ s m

1728jcpfpa,VcJc,V

 

(1-</*)3

 

 

(1U2)

Коэффициент Xjo = х10 ’ , и, следовательно, с помощью (2.90), меняя в них формально индексы "1" и "0", местами имеем

f j

N3 cuR

ПРИ

_ _

(2.100a)

х ю =2

“V

^ < 0 1D.

v0lDJ

 

 

 

для Cl, > C0j и при T > 0!D

 

 

_ 288p?pflfc3c1Jf t( l - ^ ) 3

(2.1006)

Xin —

 

 

 

KTG№owQw (3 + d )

 

Д Л Я C U < Cgs И при T> 01D

 

 

х ш -

 

 

 

( 2 . 1 0 0 B )

Зк(«ы)2Й2(ЦгМ„)2Г

 

 

Наконец, x n = х01*, и с помощью (2.82) имеем

 

.3

 

 

 

(2. ЮОг)

х п = Г е

ПРИ

T > 9 1D>

8n2T

 

 

 

 

x, i —

" V

при

r < 0 1D.

 

V0iDy

 

 

 

 

90