Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4
Дано обоснование применимости к кинетической системе радиационных дефектов некоторых требований, характерных для компонент химической кинетической системы Гирера и Мейнхардта [1, 2] (например, автокатализ компонент). Проведенный сравнительный анализ этих кинетических систем позволил существенно уточ...
Saved in:
Date: | 2007 |
---|---|
Main Authors: | , , , , |
Format: | Article |
Language: | Russian |
Published: |
Одеський національний політехнічний університет України
2007
|
Series: | Вопросы атомной науки и техники |
Subjects: | |
Online Access: | http://dspace.nbuv.gov.ua/handle/123456789/110624 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
Journal Title: | Digital Library of Periodicals of National Academy of Sciences of Ukraine |
Cite this: | Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 / В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская, С.А. Чернеженко, В.Д. Русов // Вопросы атомной науки и техники. — 2007. — № 6. — С. 29-35. — Бібліогр.: 12 назв. — рос. |
Institution
Digital Library of Periodicals of National Academy of Sciences of Ukraineid |
irk-123456789-110624 |
---|---|
record_format |
dspace |
spelling |
irk-123456789-1106242017-01-06T03:02:50Z Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 Тарасов, В.А. Бориков, Т.Л. Крыжановская, Т.В. Чернеженко, С.А. Русов, В.Д. Физика радиационных повреждений и явлений в твердых телах Дано обоснование применимости к кинетической системе радиационных дефектов некоторых требований, характерных для компонент химической кинетической системы Гирера и Мейнхардта [1, 2] (например, автокатализ компонент). Проведенный сравнительный анализ этих кинетических систем позволил существенно уточнить кинетическую систему для дефектов физической нелинейной системы «металл + нагрузка + облучение» [3, 4]. Получена упрощенная кинетическая система для дефектов нелинейной физической системы «металл + нагрузка + облучение» сходная с кинетической системой типа Гирера-Мейнхардта. Проведены вычислительные компьютерные эксперименты по исследованию упрощенная кинетическая система для дефектов нелинейной физической системы «металл + нагрузка + облучение». Полученные результаты подтверждают образование диссипативных структур точечных дефектов. Дано обґрунтування застосовності до кінетичної системи радіаційних дефектів деяких вимог, характерних для компонентів хімічної кінетичної системи Гірера й Мейнхардта [1, 2] (наприклад, автокаталіз компонентів). Проведений порівняльний аналіз цих кінетичних систем дозволив істотно уточнити кінетичну систему для дефектів фізичної нелінійної системи «метал + навантаження + опромінення» [3, 4]. Отримано спрощену кінетичну систему для дефектів нелінійної фізичної системи «метал + навантаження + опромінення» подібну з кінетичною системою типу Гірера-Мейнхардта. Проведені обчислювальні комп'ютерні експерименти по дослідженню спрощеної кінетичної системи для дефектів нелінійної фізичної системи «метал + навантаження + опромінення». Отримані результати підтверджують утворення дисипативних структур крапкових дефектів. The applicability ground of some requirements characteristic for the components of Girer and Meinhardt chemical kinetic system [1, 2] (for example, component autocatalysis) to the radiation defect kinetic system is given. Our fulfilled comparative analysis of these kinetic systems has permitted essentially to amend the kinetic system for the defects of nonlinear physical system "meta l+ load + irradiation" [3, 4]. The simplified kinetic system for defects of nonlinear physical system "metal + load + irradiation" similar to Girer-Meinhardt kinetic system is received. Computer experiments on research the simplified kinetic system for defects of nonlinear physical system "metal + load + irradiation" are lead. The received results confirm formation of dissipative structures of point defects. 2007 Article Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 / В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская, С.А. Чернеженко, В.Д. Русов // Вопросы атомной науки и техники. — 2007. — № 6. — С. 29-35. — Бібліогр.: 12 назв. — рос. 1562-6016 http://dspace.nbuv.gov.ua/handle/123456789/110624 548.4 ru Вопросы атомной науки и техники Одеський національний політехнічний університет України |
institution |
Digital Library of Periodicals of National Academy of Sciences of Ukraine |
collection |
DSpace DC |
language |
Russian |
topic |
Физика радиационных повреждений и явлений в твердых телах Физика радиационных повреждений и явлений в твердых телах |
spellingShingle |
Физика радиационных повреждений и явлений в твердых телах Физика радиационных повреждений и явлений в твердых телах Тарасов, В.А. Бориков, Т.Л. Крыжановская, Т.В. Чернеженко, С.А. Русов, В.Д. Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 Вопросы атомной науки и техники |
description |
Дано обоснование применимости к кинетической системе радиационных дефектов некоторых требований, характерных для компонент химической кинетической системы Гирера и Мейнхардта [1, 2] (например, автокатализ компонент). Проведенный сравнительный анализ этих кинетических систем позволил существенно уточнить кинетическую систему для дефектов физической нелинейной системы «металл + нагрузка + облучение» [3, 4]. Получена упрощенная кинетическая система для дефектов нелинейной физической системы «металл + нагрузка + облучение» сходная с кинетической системой типа Гирера-Мейнхардта. Проведены вычислительные компьютерные эксперименты по исследованию упрощенная кинетическая система для дефектов нелинейной физической системы «металл + нагрузка + облучение». Полученные результаты подтверждают образование диссипативных структур точечных дефектов. |
format |
Article |
author |
Тарасов, В.А. Бориков, Т.Л. Крыжановская, Т.В. Чернеженко, С.А. Русов, В.Д. |
author_facet |
Тарасов, В.А. Бориков, Т.Л. Крыжановская, Т.В. Чернеженко, С.А. Русов, В.Д. |
author_sort |
Тарасов, В.А. |
title |
Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 |
title_short |
Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 |
title_full |
Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 |
title_fullStr |
Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 |
title_full_unstemmed |
Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 |
title_sort |
теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». часть 4 |
publisher |
Одеський національний політехнічний університет України |
publishDate |
2007 |
topic_facet |
Физика радиационных повреждений и явлений в твердых телах |
url |
http://dspace.nbuv.gov.ua/handle/123456789/110624 |
citation_txt |
Теория диссипативных структур кинетической системы для дефектов нелинейной физической системы «металл+нагрузка+облучение». Часть 4 / В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская, С.А. Чернеженко, В.Д. Русов // Вопросы атомной науки и техники. — 2007. — № 6. — С. 29-35. — Бібліогр.: 12 назв. — рос. |
series |
Вопросы атомной науки и техники |
work_keys_str_mv |
AT tarasovva teoriâdissipativnyhstrukturkinetičeskojsistemydlâdefektovnelinejnojfizičeskojsistemymetallnagruzkaoblučeniečastʹ4 AT borikovtl teoriâdissipativnyhstrukturkinetičeskojsistemydlâdefektovnelinejnojfizičeskojsistemymetallnagruzkaoblučeniečastʹ4 AT kryžanovskaâtv teoriâdissipativnyhstrukturkinetičeskojsistemydlâdefektovnelinejnojfizičeskojsistemymetallnagruzkaoblučeniečastʹ4 AT černeženkosa teoriâdissipativnyhstrukturkinetičeskojsistemydlâdefektovnelinejnojfizičeskojsistemymetallnagruzkaoblučeniečastʹ4 AT rusovvd teoriâdissipativnyhstrukturkinetičeskojsistemydlâdefektovnelinejnojfizičeskojsistemymetallnagruzkaoblučeniečastʹ4 |
first_indexed |
2025-07-08T00:53:55Z |
last_indexed |
2025-07-08T00:53:55Z |
_version_ |
1837038081346109440 |
fulltext |
УДК 548.4
ТЕОРИЯ ДИССИПАТИВНЫХ СТРУКТУР КИНЕТИЧЕСКОЙ
СИСТЕМЫ ДЛЯ ДЕФЕКТОВ НЕЛИНЕЙНОЙ ФИЗИЧЕСКОЙ
СИСТЕМЫ «МЕТАЛЛ+НАГРУЗКА+ОБЛУЧЕНИЕ».
ЧАСТЬ 4
В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская, С.А. Чернеженко, В.Д. Русов
Одесский национальный политехнический университет, г. Одесса, Украина;
E-mail: siiis@te.net.ua, тел. (0482)64-16-72
Дано обоснование применимости к кинетической системе радиационных дефектов некоторых
требований, характерных для компонент химической кинетической системы Гирера и Мейнхардта [1, 2]
(например, автокатализ компонент). Проведенный сравнительный анализ этих кинетических систем
позволил существенно уточнить кинетическую систему для дефектов физической нелинейной системы
«металл + нагрузка + облучение» [3, 4]. Получена упрощенная кинетическая система для дефектов
нелинейной физической системы «металл + нагрузка + облучение» сходная с кинетической системой типа
Гирера-Мейнхардта. Проведены вычислительные компьютерные эксперименты по исследованию
упрощенная кинетическая система для дефектов нелинейной физической системы «металл + нагрузка +
облучение». Полученные результаты подтверждают образование диссипативных структур точечных
дефектов.
УПРОЩЕННАЯ КИНЕТИЧЕСКАЯ
СИСТЕМА ДЛЯ ДЕФЕКТОВ
ТОПЛИВНЫХ МЕТАЛЛОВ
И ЕЕ МОДЕЛИРОВАНИЕ
Из-за сложности кинетической системы для
точечных дефектов ((21) из части 1 статьи [10]) для
наглядного подтверждения ее автокаталитических
особенностей (сходство с системой Гирера-
Мейнхардта), что является принципиальным для
данной статьи, сначала представляется разумным
перейти к упрощенной кинетической системе для
точечных дефектов и ее компьютерному
моделированию.
Упрощенную кинетическую систему для
точечных дефектов делящегося металла под
облучением и под действием нагрузки можно
получить из (21) [10] статьи следующим образом.
При фиксированной температуре и плотности
потока нейтронов кинетическая система (21) [10]
принимает следующий вид:
( ) ( ) ( )
( ) ( )
( ) ( ) ( )
( ) ( )
( ) ( ) ( );,,
,,,~,
,
,,
,,
,,
,
trCtrCT
tTrqtrq
trE
kT
trCTrD
trCTrD
trCTrD
t
trC
VI
Irrad
I
S
I
I
II
II
II
I
rr
rr
r
rr
rr
rr
r
α−
−Φ+−
−⎥⎦
⎤
⎢⎣
⎡ ∇∇+
+∇∇+
+Δ=
∂
∂
(1)
( ) ( ) ( )
( ) ( )
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ).,,,,,~
,,
,,
,,
,,
,
trCtrCTtTrq
trqtrE
kT
trCTrD
trCTrD
trCTrD
t
trC
VI
Irrad
V
S
VV
VV
VV
VV
V
rrr
rr
rr
rr
rr
r
α−Φ+
+−⎥⎦
⎤
⎢⎣
⎡ ∇∇+
+∇∇+
+Δ=
∂
∂
(2)
Если считать диффузию изотропной и
пренебречь дрейфовой диффузией, то для (1) и (2)
получим следующий вид:
( ) ( ) ( ) ( )
( ) ( ) ( ) ( );,,,,,~
,,,
trCtrCTtTrq
trqtrCTD
t
trC
VI
Irrad
I
S
III
I
rrr
rr
r
α−Φ+
+−Δ=
∂
∂
(3)
( ) ( ) ( ) ( )
( ) ( ) ( ) ( ).,,,,,~
,,
,
trCtrCTtTrq
trqtrCTD
t
trC
VI
Irrad
V
S
VVV
V
rrr
rr
r
α−Φ+
+−Δ=
∂
∂
(4)
_______________________________________________________________________________
ВОПРОСЫ АТОМНОЙ НАУКИ И ТЕХНИКИ. 2007. № 6. 29
Серия: Физика радиационных повреждений и радиационное материаловедение (91), с. 29-35.
Последним уравнениям можно согласно
вышеизложенному [10, 11, 12] придать следующий
вид:
( ) ( ) ( ) ( )
( ) ( ) ( );,,
,,
trCtrCTS
QtrCTDK
t
trC
VII
IIII
I
rr
r
r
α−+
++Δ+Φ=
∂
∂
(5)
( ) ( ) ( ) ( )
( ) ( ) ( ),,,
,,
trCtrCTSQ
trCTDK
t
trC
VIVV
VVV
V
rr
r
r
α−++
+Δ+Φ=
∂
∂
(6)
где и – скорости генерации
облучением с плотностью потока Φ междоузельных
атомов и вакансий соответственно; и –
автокаталитические члены, связанные с рождением
междоузельных атомов и вакансий соответственно
при движении дислокаций со ступеньками; и
- сток междоузельных атомов и вакансий
соответственно на дислокации.
( )ΦIK ( )ΦVK
IQ VQ
IS VS
Дальнейшее упрощение кинетической системы
(5), (6) связано с упрощением членов , и ,
, что возможно лишь при упрощенном описании
дислокационной системы (вклад пор и других
дефектов здесь не рассматривается), т.е. при отказе
от компьютерного моделирования [12]. Поэтому
полагаем, что все дислокации всегда
прямолинейные (жесткие дислокации), их средняя
длина
IQ VQ IS
VS
l , плотность дислокаций ρ . Упрощая, также
полагаем, что автокаталитические члены, связанные
с рождением междоузельных атомов и вакансий
соответственно при движении дислокаций со
ступеньками, могут быть заданы в следующем виде:
дисл
I
I vAQ = и , (7) дисл
V
V vAQ =
где IA и VA – константы, характеризующие
скорости рождения междоузельных атомов и
вакансий единичным участком длины дислокации
соответственно; - скорость движения
дислокации.
дислv
Для скорости движения прямолинейной
дислокации можно записать приблизительное
выражение (справедливое при >> , где и
– скорость скольжения и скорость
переползания дислокации соответственно) согласно
(7) из [12] в виде
скv перv скv
перv
пердисл vv ⋅=
λ
l
, (8)
где – среднее расстояние скольжения дислокации
до остановки на препятствиях ( l ~ расстояние
между препятствиями в плоскости скольжения в
направлении скольжения); λ – расстояние, на
которое должна переползти дислокация для
преодоления препятствия (λ обычно полагается
равным размеру препятствия).
l
Теперь выражения (7) можно записать
следующим образом:
пер
I
I vAlQ
λ
ρ l
= и пер
V
V vAlQ
λ
ρ l
= . (9)
Скорость переползания дислокации задается
выражением (33) из [12]:
,
2
2
exp
1lnexpexp
)/ln(
4
2
2
exp
1lnexpexp
)/ln(
4
0
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛ −Δ
×
×
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−⋅⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛ ⋅
×
×−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛ −Δ
×
×
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛ ⋅
×
×=−=
kT
FW
C
C
LbkT
vF
bzakT
D
kT
FW
C
C
LbkT
vF
bzakT
D
VVv
ÑÒ
V
S
O
V
Va
V
óïð
V
V
SÑÒ
I
S
I
Ia
I
óïð
I
I
SV
ïåð
I
ïåðïåð
β
π
β
π
(10)
где
b
bкр=β .
Если для упрощения выкладок обозначить
константы в (33) [12] через
)/ln(
4
1 bzakT
DC
I
I
Sπ
= ,
)/ln(
4
2 bzakT
DC
I
I
Sπ
= ;
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛ −Δ
=
kT
FWCC СТ
I
S
2
2exp13 ;
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛ −Δ
=
kT
FWCC СТ
V
S
2
2exp24 и , (11) 435 CCC +=
то выражение (10) можно записать в следующем
виде:
( )
( ) ,exp
exp
504
03
C
C
C
FC
C
C
FCv
V
V
V
I
I
Iïåð
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
=
−β
β
(12)
30
где LbkT
vF
F a
I
упр
I = и LbkT
vF
F a
V
упр
V = .
Тогда согласно (9) получим
пер
I
РФпер
I
I vAlvAlQ
λ
ρρ
λ
ρ ll )( 0 −+== (13)
и
пер
V
РФпер
V
V vAlvAlQ
λ
ρρ
λ
ρ ll )( 0 −+== , (14)
где 0ρ – начальная плотность дислокаций; РФ−ρ –
плотность дислокаций, создаваемых источниками
дислокаций типа Франка-Рида.
В первом приближении можно принять, что
пер
РФ
РФ VC ⋅= −
−ρ , (15)
где . constC РФ ≈−
Поэтому для (13) и (14) с учетом (15) получим:
2
0 пер
IРФ
пер
I
I vAlCvAlQ
λλ
ρ ll −+= ; (16)
2
0 пер
VРФ
пер
V
V vAlCvAlQ
λλ
ρ ll −+= . (17)
Преобразуем (16), используя выражение (12) для
: перv
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+=
=
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
×
×+
+
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
×
×=
−
−
−
−
ββ
ββ
ββ
λ
λ
ρ
08076
2
50403
50403
0
V
VF
I
IF
V
VF
I
IF
IРФ
V
VF
I
IF
I
I
C
CeC
C
CeCC
C
C
CeC
C
CeC
AlC
C
C
CeC
C
CeC
AlQ
VI
VI
VI
l
l
( )
,
2
0
2
110
0
10
2
0
2
9
βββ
β
−
+
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
×
×⋅−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+
V
VF
V
I
I
V
FF
I
IF
C
C
eC
C
C
C
C
eC
C
C
eC
V
VII
(18)
где
.
;2
;
;2
;2
;
2
411
4310
2
39
54408
53307
2
5506
CAlCC
CCAlCC
CAlCC
CCAlCCAlC
CCAlCCAlC
CAlCCAlC
IÐÔ
IÐÔ
IÐÔ
IÐÔI
IÐÔI
IÐÔI
λ
λ
λ
λλ
ρ
λλ
ρ
λλ
ρ
l
l
l
ll
ll
ll
−
−
−
−
−
−
=
=
=
+=
⋅+⋅=
+=
Аналогичным образом можно получить
выражение для : VQ
,
2
0
2
17
0
0
)(
16
2
0
2
15014
01312
β
ββ
ββ
β
−
+
−
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+
+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅−
−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅−
−⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
⋅+=
V
VF
V
I
I
VFF
I
IF
V
VF
I
IF
V
C
C
eC
C
C
C
C
eC
C
C
eC
C
C
eC
C
C
eCCQ
V
VI
IV
I
(19)
где
.;
;;
;;
11
17
10
16
9
15
8
14
7
13
6
12
V
I
V
I
V
I
V
I
V
I
V
I
A
A
CCA
A
CC
A
A
CCA
A
CC
A
A
CCA
A
CC
⋅=⋅=
⋅=⋅=
⋅=⋅=
Теперь рассмотрим и . Потоки точечных
дефектов на единицу длины дислокации задаются
выражениями (24) и (25) [12], которые можно
записать в виде:
IS VS
,14 0
0
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅⋅=
β
π
I
IF
III C
CeCbDI I (20)
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅⋅=
−
14 0
0
β
π
V
VF
VVV C
CeCbDI V . (21)
Пусть и ,
тогда (20) и (21) имеет вид:
0
18 4 II CbDС ⋅= π 0
19 4 VV CbDС ⋅= π
;18018 Ñ
C
CeÑI
I
IF
I I −⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅=
β
22)
31
32
19019 С
C
CeСI
V
VF
V
I −⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅=
β
. (23)
Тогда для и получаем следующие
выражения:
IS VS
[ ]
,
0
0
)(
24
2
0
2
23
02202120
180180
βββ
ββ
β
ρρ
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅+−=
=
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅×⋅⋅+==
+
−
−
I
V
V
IFF
I
IF
V
VF
I
IF
I
IF
пер
РФ
II
C
C
C
C
eC
C
C
eC
C
C
eC
C
C
eCC
С
C
C
eСlvCIlS
VII
VI
I
(24)
где
41824318
2341822
318
51818021
51818020
;
;
;
;
CÑlCCCÑlC
CCÑlCC
CÑlC
CÑlCÑlC
CÑlCÑlC
ÐÔÐÔ
ÐÔ
ÐÔ
ÐÔ
ÐÔ
⋅⋅⋅=⋅⋅⋅=
=⋅⋅⋅=
⋅⋅⋅−
−⋅⋅⋅+⋅⋅=
⋅⋅⋅+⋅⋅=
−−
−
−
−
−
ρ
ρ
и
[ ]
,0
0
)(
29
2
0
2
28
02702625
19019
0
βββ
ββ
β
ρρ
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅−−=
=
⎥
⎥
⎦
⎤
⎢
⎢
⎣
⎡
−⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
⋅×
×⋅⋅+==
+
−
−
−
−
I
V
V
IFF
V
VF
V
VF
I
IF
V
VF
пер
РФ
VV
C
C
C
C
eC
C
C
eC
C
C
eC
C
C
eCC
С
C
C
eС
lvCIlS
VIV
VI
V
(25)
где
.
;;
;
;
31929
41928419
51919027
31926
51919025
CСlCC
CСlCCCСlC
CСlCСlC
CСlCC
CСlCСlC
РФ
РФРФ
РФ
РФ
РФ
⋅⋅⋅=
⋅⋅⋅=⋅⋅⋅−
−⋅⋅⋅+⋅⋅=
⋅⋅⋅=
⋅⋅⋅+⋅⋅=
−
−−
−
−
−
ρ
ρ
Подставляя полученные выражения для ,
и , в уравнения кинетической системы (5) и
(6), после алгебраических преобразований получаем
следующий вид кинетической системы:
IQ VQ
IS VS
+⋅++
+−+
+−=
∂
∂
−
−
ββ
ββ
α
22
34
22
33
3231
)()(
)()(
),(),()(~),(
V
F
I
F
V
F
I
F
VII
I
CeCCeC
CeCCeC
trCtrCФK
t
trC
VI
VI
rr
r
(26,а)
);,()(
35 trCD
C
C
eC II
V
IFF VI
r
Δ+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+ +
β
),,(
)()(
)()(
),(),()(~),(
)(
41
22
40
22
39
3837
trCD
C
C
eC
CeCCeC
CeCCeC
trCtrCÔK
t
trC
VV
V
IFF
V
F
I
F
V
F
I
F
VIV
V
VI
VI
VI
r
rr
r
Δ+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−
−⋅++
+−+
+−=
∂
∂
+
−
−
β
ββ
ββ
α
(26,б)
где введены обозначения
;)(
;)(;))((
;))((
;))((;
0
0
102435
20
1134
20
23933
0
22832
0
2173120630
β
ββ
β
β
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
−=
=−=
+=
−=+=
−
−
I
V
VI
V
I
C
CCCC
CCCCCCC
CCCC
CCCCCCC
;)(
;))((
;;))((
;))((;
0
0
231641
20
281740
1539
0
271438
0
261337251236
β
β
β
β
⎟
⎟
⎠
⎞
⎜
⎜
⎝
⎛
+=
+=
=−=
+=+= −
I
V
V
V
I
C
CCCC
CCCC
CCCCCC
CCCCCCC
.)()(~;)()(~
3630 CФKФKCФKФK VVII +=+=
Если пренебречь внутренними локальными
напряжениями по сравнению с внешним
напряжением и ввести обозначения
,
;;;
;;;
;;;
)(
4110
2
409
2
398387
376
)(
355
2
344
2
3333232311
VI
VIV
IVIV
IVI
FF
FFF
FFFF
FFF
eCK
eCKeCKeCK
eCKeCKeCK
eCKeCKeCK
+
+
=
⋅===
==⋅=
===
то получим окончательный вид кинетической
системы:
),,(
),(
),(
)),(()),((
)),(()),((
),(),()(~),(
);,(
),(
),(
)),(()),((
)),(()),((
),(),()(~),(
10
2
9
2
8
76
5
2
4
2
3
21
trCD
trC
trCK
trCKtrCK
trCKtrCK
trCtrCÔK
t
trC
trCD
trC
trCK
trCKtrCK
trCKtrCK
trCtrCÔK
t
trC
VV
V
I
VI
VI
VIV
V
II
V
I
VI
VI
VII
I
r
rr
r
r
rr
r
Δ+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
−
−++
+−+
+−=
∂
∂
Δ+⎟⎟
⎠
⎞
⎜⎜
⎝
⎛
+
+++
+−+
+−=
∂
∂
−
−
−
−
β
ββ
ββ
β
ββ
ββ
α
α
(27)
где все – константы; 10...1, =iKi β – также
константа для заданного типа дислокаций
( 0,10,0 ≤≤ β ).
Таким образом, получили кинетическую
систему, аналогичную кинетической системе типа
Гирера-Мейнхардта [1, 2].
Для решения системы надо задать начальные и
граничные условия:
( ) constCtrC VV === 00,
r
и
( ) constCtrC II === 00,
r
;
и
( ) constCtrC VграницаV == 0,r
и ( ) constCtrC IграницаI == 0,
r
. (28)
Кинетическая система (27) с начальными и
граничными условиями (28) в настоящее время
исследуется нами с помощью компьютерного
эксперимента. В процессах исследований
используются основные математические пакеты для
численных расчётов: Mathematica 5.1.1…5.2.0,
Maple 10, MatLab 7 и другие, а также специально
написанные под эту задачу приложения. С учётом
возможности распараллеливания процессов
вычислений, а также поддержки 64-битных
вычислений в последних версиях этих пакетов
основные расчёты и анализ будут проведены на
мультипроцессорных системах. Несмотря на
сложность и нестандартный тип системы, в
результате исследований её структуры были
найдены области, расчёт в которых можно провести
с высокой точностью (порядка 0.1%) и на
достаточных для предварительного анализа
интервалах. На данный момент основные расчёты
проводились для
2-мерных ситуаций при β порядка 1 (в частности, на
приведённых ниже рисунках для β =0.9), так как
естественно предположить, что именно
максимальный порядок этого показателя сможет
выявить наиболее характерные эффекты, связанные
с нелинейностью и диссипацией. Ограничения 2-
мерностью вызвано исключительно техническими
соображениями и невозможностью на момент
проведения расчётов использовать суперЭВМ.
Следует также отметить возможность ”расширения”
областей счёта, так сказать, искусственным образом
в результате использования многомерной
экстраполяции на основе полученных точек в
считаемых областях. Таким образом, можно будет
расширить область графиков в несколько раз, а
также закрыть проблемные места в самой картинке
(иногда из-за сбоя метода в нём возникают
“дырки”). Приведенные ниже графики (рис. 1-5)
являются результатом компьютерного
моделирования упрощённой системы, входящими
параметрами и начальными условиями которой
являются величины наиболее характерные, на наш
взгляд, для описываемых выше ситуаций. Помимо
этого они подбирались из соображений
максимальной наглядности иллюстрации
диссипативности структур. А именно: температуру
Т мы принимали равной 400 К, начальные
концентрации вакансий и междоузельных атомов
10 и 10 соответственно, соответствующая
плотность потока нейтронов 10 13
10− 30−
ññìíåéòð 2. ,
нагрузка 3 2ммкГ . На представленных ниже
рисунках расчёты приведены для различных
моментов времени с шагом 0.1 с, т.е. ситуация
рассматривается как бы в динамике, и мы можем
наблюдать эволюционирование процессов во
времени. Причём можно заметить, что имеется
цикличность процессов. А учитывая размеры
просчитанных областей, можно говорить о
макроскопической пространственной корреляции с
радиусом, значительно превышающим радиус
потенциала межатомного взаимодействия, что
характерно для диссипативных, нелинейных
структур.
Рис. 1. Результаты моделирования упрощённой
кинетической системы при t=0.1 c
(x, y – линейные размеры кристалла в 10-4cм;
Сi – концентрация междоузельных атомов,
умноженная на 10-10)
33
Рис. 2. Результаты моделирования упрощённой
кинетической системы при t=0.4 c (x, y – линейные
размеры кристалла в 10-4cм; Сi – концентрация
междоузельных атомов, умноженная на 10 ) 10−
Рис. 3. Результаты моделирования упрощённой
кинетической системы при t=0.6 c (x, y – линейные
размеры кристалла в 10-4cм; Сi – концентрация
междоузельных атомов, умноженная на 10-10)
Рис. 4. Результаты моделирования упрощённой
кинетической системы при t=0.8 c (x, y – линейные
размеры кристалла в 10-4cм; Сi – концентрация
междоузельных атомов, умноженная на 10-10)
Рис. 5. Результаты моделирования упрощённой
кинетической системы при t=1c (x, y – линейные
размеры кристалла в cм; Сi – концентрация
междоузельных атомов, умноженная на 10
410−
-10)
Данные расчёты являются предварительными, в
дальнейшем планируется проведение более
детального анализа уравнений в 3D с
использованием кластерных систем, а также
выявление внешних ограничений (начальные,
граничные условия и др. параметры), определяющих
степень и вид проявления диссипативности в этой
системе. На этот момент были рассмотрены
уравнения с начальными концентрациями в
диапазоне 10 8…10-10 и 10-26…10-30, причём в этих
интервалах диссипативность в той или иной степени
проявила себя и наиболее часто встречаемым
распределением для пространственной
концентрации дефектов была поверхность с двумя
или тремя пиками (аналогично рис. 3), либо
поверхность в виде достаточно резко изогнутого
плато.
3. ВЫВОДЫ
Получена упрощенная кинетическая система для
дефектов нелинейной физической системы «металл
+ нагрузка + облучение», сходная с кинетической
системой типа Гирера-Мейнхардта [1,2] (система (1)
[10].
Проведены вычислительные компьютерные
эксперименты по исследованию упрощенной
кинетической системы для дефектов нелинейной
физической системы «металл + нагрузка +
облучение». Полученные результаты подтверждают
образование диссипативных структур точечных
дефектов.
ЛИТЕРАТУРА
1. A. Gierer, H. Meinhardt. Lectures on Mathematics
in the Life Sciences. Berlin: Springer, 1974, v. 7,
p. 163–183.
2. H. Meinhardt. Pattern formation in biology: a com-
parison of models and experiments /UK: IOP Pub-
lishing Ltd. Reports On Progress In Physics, 1992,
v. 55, p. 797–849.
3. В.Д. Русов, В.А. Тарасов и др. Моделирование
диссипативных структур и волн концентраций
точечных дефектов в открытой нелинейной
физической системе «металл + нагрузка +
облучение» //Вопросы атомной науки техники.
Серия «Физика радиационных повреждений и
радиационное материаловедение» (80). 2001,
№4, с. 3–8.
4. В.Д. Русов, В.А. Тарасов и др. Механизмы
радиационной ползучести металлического
уранового топлива и ее температурные
нелинейности //Вопросы атомной науки
техники. Серия «Физика радиационных
повреждений и радиационное
материаловедение» (84). 2003, № 6, с. 20–23.
34
35
5. Режимы с обострением. Эволюция идеи: законы
коэволюции сложных структур. М.: «Наука»,
1998, 255 с.
6. В.А. Тарасов. Математическое моделирование
радиационной ползучести реакторного топлива
на примере урана и его сплавов //Вопросы
атомной науки техники. Серия «Физика
радиационных повреждений и радиационное
материаловедение» (79). 2001, №2, с. 23–30.
7. В.А. Тарасов. Комп’ютерне моделювання
радіаційной повзучості реакторного палива
//УФЖ. 2000, №10, с. 23–35.
8. W.G. Wolfer, L.K. Mansur, J.A. Sprague //Rad. Eff.
in Breeder Reactor Structural Materials. Scottsdale,
Ariz., 1977, 841 p.
9. Ш.Ш. Ибрагимов и др. Радиационное
повреждение металлов и сплавов М.:
«Энергоатомиздат», 1985, 239 с.
10. В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская,
С.А. Чернеженко, В.Д. Руссов. Теория
диссипативных структур кинетической системы
для дефектов нелинейной физической системы
«металл + нагрузка + облучение». Ч. 1 //Вопросы
атомной науки техники. Серия «Физика
радиационных повреждений и радиационное
материаловедение» (90). 2007, №2, с. 63–71.
11. В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская,
С.А. Чернеженко, В.Д. Руссов. Теория
диссипативных структур кинетической системы
для дефектов нелинейной физической системы
«металл + нагрузка + облучение». Ч. 2 //Вопросы
атомной науки техники. Серия «Физика
радиационных повреждений и радиационное
материаловедение» (90). 2007, №2, с. 72-75.
12. Наст. журнал, с. 29-34.
ТЕОРІЯ ДИСИПАТИВНИХ СТРУКТУР КІНЕТИЧНОЇ СИСТЕМИ ДЛЯ ДЕФЕКТІВ НЕЛІНІЙНОЇ
ФІЗИЧНОЇ СИСТЕМИ «МЕТАЛ+НАВАНТАЖЕННЯ+ОПРОМІНЕННЯ».
ЧАСТИНА 4
В.О. Тарасов, Т.Л. Боріков, Т.В. Крижанівська, С.А. Чернеженко, В.Д. Русов
Дано обґрунтування застосовності до кінетичної системи радіаційних дефектів деяких вимог, характерних для
компонентів хімічної кінетичної системи Гірера й Мейнхардта [1, 2] (наприклад, автокаталіз компонентів). Проведений
порівняльний аналіз цих кінетичних систем дозволив істотно уточнити кінетичну систему для дефектів фізичної
нелінійної системи «метал + навантаження + опромінення» [3, 4]. Отримано спрощену кінетичну систему для дефектів
нелінійної фізичної системи «метал + навантаження + опромінення» подібну з кінетичною системою типу Гірера-
Мейнхардта. Проведені обчислювальні комп'ютерні експерименти по дослідженню спрощеної кінетичної системи для
дефектів нелінійної фізичної системи «метал + навантаження + опромінення». Отримані результати підтверджують
утворення дисипативних структур крапкових дефектів.
THE THEORY OF DISSIPATIVE STRUCTURES OF THE KINETIC SYSTEM FOR DEFECTS
OF NONLINEAR PHYSICAL SYSTEM “METAL + LOADING + IRRADIATION”.
PART 4
V.A. Tarasov, T.L. Borikov, T.V. Kryzhanovskaja, S.A. Chernegenko, V.D. Rusov
The applicability ground of some requirements characteristic for the components of Girer and Meinhardt chemical kinetic
system [1, 2] (for example, component autocatalysis) to the radiation defect kinetic system is given. Our fulfilled comparative
analysis of these kinetic systems has permitted essentially to amend the kinetic system for the defects of nonlinear physical sys-
tem "meta l+ load + irradiation" [3, 4]. The simplified kinetic system for defects of nonlinear physical system "metal + load +
irradiation" similar to Girer-Meinhardt kinetic system is received. Computer experiments on research the simplified kinetic sys-
tem for defects of nonlinear physical system "metal + load + irradiation" are lead. The received results confirm formation of dis-
sipative structures of point defects.
В.А. Тарасов, Т.Л. Бориков, Т.В. Крыжановская, С.А. Чернеженко, В.Д. Русов
УПРОЩЕННАЯ КИНЕТИЧЕСКАЯ СИСТЕМА ДЛЯ ДЕФЕКТОВ ТОПЛИВНЫХ МЕТАЛЛОВ И ЕЕ МОДЕЛИРОВАНИЕ
В.О. Тарасов, Т.Л. Боріков, Т.В. Крижанівська, С.А. Чернеженко, В.Д. Русов
V.A. Tarasov, T.L. Borikov, T.V. Kryzhanovskaja, S.A. Chernegenko, V.D. Rusov
|