Математика. Прикладная математика
УДК 517.957
Е.М. Веселова
АНАЛИЗ СВЯЗАННОЙ СИСТЕМЫ УРАВНЕНИЙ ДЛЯ ДЕФОРМАЦИИ СРЕДЫ
И КОНЦЕНТРАЦИИ ДЕФЕКТОВ
В статье приведено решение связанной системы уравнений для деформации среды и концентрации дефектов в стационарном случае. Проведен анализ решения на устойчивость и показано, что в зависимости от превышения контролируемого параметра над пороговым значением образуются периодические структуры точечных дефектов.
ANALYSIS OF A COUPLED SYSTEM OF EQUATIONS FOR CONTINUUM DEFORMATION AND DEFECT CONCENTRATION
The article presents a solution of coupled system of equations describing continuum deformation and defect concentration within the stationary approach. The stability analysis of obtained solution was performed to indicate the arising of point defect periodical structure, which depends on control parameter value with respect to estimated threshold.
Системы, в которых происходят диффузионные процессы, относятся к динамическим системам, в которых задание начальных условий полностью определяет их поведение в последующие моменты времени [1]. Под устойчивостью динамической системы понимают реакцию на малые возмущения ее состояния. Если сколь угодно малые изменения состояния системы начинают нарастать со временем, то система является неустойчивой. Если малые возмущения затухают со временем, то система устойчива [1]. Формирование новых структур всегда сопровождается потерей устойчивости предшествующих. Линейный анализ устойчивости позволяет предсказать, при каких значениях параметров состояния системы развивается ее неустойчивость [1, 2].
Будем рассматривать случай одномерного кристалла, когда тензор деформации среды Ы;;
имеет единственную х-компоненту [3]:
Введение
Математическое описание
где Ux - компонента вектора смещения среды, а ось x направлена вдоль одной из осей типа [100] кубического кристалла.
Выражение для плотности свободной энергии одномерного деформированного ангармонического кристалла с дефектами с учетом нелокальности взаимодействия атомов решетки с друг другом и с дефектами имеет вид:
F = 1 pcX(x)+pcfl20 Х + lapcfX(x)+1 bpcfX4(x)-0dnj(x)2 ox 2 2 (2)
-odnfd ^.
Здесь p - плотность среды; ci - продольная скорость звука; а, в - константы упругого ангармонизма; l0 - характеристическая длина взаимодействия атомов кристалла с друг другом; ld - характеристическая длина взаимодействия дефектов с атомами кристалла (id > i0); nd - концентрация точечных дефектов; 6d = WdKd - потенциал деформации (Wd - изменение объема кристалла при создании одного дефекта; Kd - модуль упругости); T- температура; S - энтропия системы, а<0, в>0.
Из (2) следует выражение для напряжения в кристалле с дефектами
s(x ) = —, v ; ox&
s(x) = pc2X(x) + rl2 0ХХ) + 3 apcX2 (x) + bpc2X3 (x )-0dnd. (3)
Уравнение для вектора смещения имеет вид: 0 U ct2AU + (cf - cf )grad divU.
д 2и да
Из уравнения для вектора смещения среды-¿^ =- и (3) имеем уравнение для Х(х) :
дх дх
дХ-с; Щ1+^ Щй+а ЩхЯ+р дМ -в дп. (4)
дг2 1 дх2 1 0 дх4 ; 1 дх2 ии& дх2 р дх2 Используя выражение для энергии взаимодействия точечного дефекта с упругим континуумом, следующее из (2)
И (х ) = -вл{х Уе/Л^х, получим выражение для х-компоненты нелокальной силы, действующей на дефект
х х х3
Порождаемый этой силой поток дефектов равен
£л„л (x ш+вЛот
kBt | Ox Ox
где Ва- коэффициент диффузии дефекта; кв - константа Больцмана.
С учетом (6) уравнение для концентрации дефектов запишется в виде
f \\ O x(x) ,2 O3x(x) nd (x)———+ i;— Ox d Ox 3
+ Gd-t, (7)
где 0} - скорость генерации точечных дефектов внешним источником; т} - время жизни дефекта.
Как видно из (4) и (7), плотность энтропии S, фигурирующая в (2), в рассматриваемом приближении выпадает из уравнений движений.
Введем пространственно усредненные величины
Х = J\\x(x)dx, ndо = J Jn0 {x)dx,
где Ь - размер кристалла, и представим деформацию и концентрацию дефектов в виде
Х(х) = Х + Х (х), по (х) = Па0 + п01 (х), (8)
где £ и пса(х) - пространственно неоднородные части деформации и концентрации дефектов. Пространственно однородная деформация и концентрация дефектов связаны соотношением:
Х = KTn
Система уравнений (4) и (7) имеет пространственно однородное стационарное решение ndо = ОёТё, X = Х0. Это решение оказывается неустойчивым при превышении определенной критической концентрации дефектов, возникает новое, пространственно неоднородное стационарное состояние [4]. В этом состоянии £ и па задаются формулой (8). Из (7) и (8) в пренебрежении рекомбинации дефектов (tdl = 0) получаем уравнение для n
d Э 2 ndi
Ddqdnd о
Э Х + /2 Э Х
Из (4) и (8) имеем уравнение для :
ЭХХ=с 2
,2 = С/
Э 2Xi (х)+„ Э4Хх(х)+3 a 2 Э2 (хХ (х))
- + С/2/
+ — ac 2
+ Ppc\\
э2 (ХХ (х)) ed э^
Эх2 р Эх2
Эt2 & Эх2 Из (10) следует, что signal = sign qd .
Покажем, что (9) и (10) описывают процессы самоорганизации дефектов, взаимодействующих через деформацию упругого континуума, приводящие к образованию стационарных периодических дефектно-деформационных структур.
Рассмотрим систему уравнений (9), (10) в стационарном случае, когда
Эt & дt 2
Интегрируя (9), имеем для пространственно неоднородной концентрации дефектов
Ddqdnd 0
э2Х + / 2ЭХ
Эх2 d Эх4
D Э 2ndi
Ddqdnd 0
Х + К
i + /?Э%
где полагаем две константы интегрирования равными нулю. Подставляя (11) в (10), получаем уравнение для пространственно неоднородной деформации среды:
э2Х(х)+/2Э%М+3aЭ2(хХ(х)), bэ2(хХ(х)) ednd0 (эХх) + ,2э4^1 (х)
- + —a
рС/ kBT
- + /2
Интегрируя его дважды, имеем следующее уравнение
^ +с + & & ^ & " &
дх 2 дх дх п
д& ,2 дЪХ ^ + с, + К-Ц- + сй
Х1 + с1 х + с11 + &
+ с2х + с22 +—а&2 +с3х + с33 + ЬХз + с4х + с44 е >2 д X
& + с5 х + с55 + ¡2 + с6 х + с66
д2Х(х)
&¿пс100 /2 &0
йпй0 _ /2 &(Л
&апс100 >2 &0
X3 +
&апс100 >2 &0
О 2п &¿&^ 0 /2
& 2п л >2 _ &¿&У0
^^0 _ /2
X3 +
&апс100 >2 &0
х +
-X --/X2 -+ с х + сП = 0,
где обозначено
Г п ^ Г& 2п 1 а 0 ■
&¿с у
й&1й00 ] 2
Г& 2п ^ 0 - & 2 ^ I
Г& 2п ^ 0 - & 2
а также введены критическая концентрация дефектов
Р квГвл
и коэффициенты
с, =¿"¿00 ¡2
^ 11 " Г12п
а а0 _/2
Далее рассмотрим решения уравнения (12), соответствующие случаю С1 = 0, Сп = 0 (т.е. с! = 0, с2 = 0), которые удовлетворяют условиям локализованности (кластерное решение) £ ^0 при и условию трансляционной инвариантности (периодическая структура) £\\= £>(х+т), где т -период трансляции.
Тогда из (12) имеем следующее нелинейное уравнение для самосогласованной деформации
- X - 2 /X - 8X3 = 0,
где коэффициенты а,/ 8 задаются в (12).
Решение уравнения в виде неявной функции выражается через интеграл
- ах + ;/XX + Х
дт - а Х +1/Х3+4 Х4+},
дх ; ; 4
ахх+/хх+4 х;+с
а£2 + /Х3 + 4 X4 +
где й - постоянная интегрирования. При й Ф 0 решение уравнения (15) выражается через эллиптические функции.
Рассмотрим класс решений, соответствующий случаю й = 0. Из (15) имеем
_Х_. (16)
аХХ + /Х3 + 4 XX
Интеграл (16) выражается через аналитические функции, тип которых определяется знаками коэффициента а и А, и через А обозначено
Д-4 -/1.
Используя (13), запишем А в виде
( 12П 1йПй0 _ /2 10
V пйс
Л + > 12п
( 12П 1йПй0 _ }2 10
V пйс
Л2 ( п л
V пйс
( п л 1 &&й 0
Используя связь между а и в, имеем:
< 1----< 1.
Можно полагать, что в реальных ситуациях имеет место неравенство
(I Л2 10
Л 2 а,
< 1----< 1.
Предполагая выполнение неравенства, разобьем область изменения пм на два интервала, в которых знаки а и А различны:
(I Л2 _0
а) при концентрациях дефектов ниже пороговой, когда
, (а < 0, Д < 0 )
интеграл (14) не выражается в действительных функциях, так что Х1 - 0, — - 0,
из (11) имеем пД - 0, т.е. процессов самоорганизации дефектов не происходит; б) при превышении порога по концентрации дефектов в области
> 1, (а < 0, Д > 0)
из (16) получаем с использованием ППП Mathcad для самосогласованной деформации следующее выражение:
лЦ/ + 79га _ / 2 а 2: Полученное выражение запишем в виде Х - вл
в+бЩ а2 х
где А1 - , & 1 =, В 42/
/2 л/9 га _ /2
График периодической функции изображен на рисунке.
п0 - 2пс :
(I л2 _0_
- 0,04, р = 2, 814пвй -1.
Зависимость периодической деформации от пространственной координаты в случае периодической структуры.
Видно, что в области —> 1 решение уравнения (20) описывает периодическую деформаци-пс
онную структуру. В соответствии с (20), согласованно с ней, образуется периодическая структура
дефектов, причем если вй > 0 , то дефекты локализуются в области X > 0 , если вй < 0 , то в области XI <0.
Заключение
Проведенный анализ устойчивости решения системы уравнений показывает иерархию процессов мелкомасштабной самоорганизации точечных дефектов, образуемых вследствие воздействия на твердое тело внешним энергетическим пучком. С ростом интенсивности облучения в кристалле усиливается пространственно однородная концентрация дефектов, являющаяся внешним контролируемым параметром. Имеются два пороговых значения контрольного параметра, при достижении которых в системе дефектов, взаимодействующих через поле упругих смещений, отмечаются качественно различные процессы самоорганизации. При малых концентрациях дефектов процессов самоорганизации не происходит и пространственное распределение дефектов остается однородным. При превышении второго порога происходит рост концентрации дефектов, образуется периодическая структура дефектов, описываемая формулой (20) и изображенная на рисунке.
Найдено решение связанной системы уравнений для деформации среды и концентрации дефектов в стационарном случае. Проведен анализ решения и показано, что в зависимости от превышения контролируемого параметра над пороговым значением образуются периодические структуры точечных дефектов.