![](/user_photo/_userpic.png)
Деев Основы расчета судовых ЯЕУ 2012
.pdf![](/html/65386/144/html_qoZK0JrvTY.YdzI/htmlconvd-H2Kesy131x1.jpg)
После коррекции изотопного состава программа возвращается к расчету kэф, и эти расчеты повторяются до тех пор, пока kэф не станет равным kэф(0) с точностью 1, т.е.
kэф kэф(0)
kэф
1.
Уравнение (5.4), описывающее пространственно-энергетическое распределение плотности потока нейтронов в двумерном цилиндрическом реакторе, в рамках многогруппового диффузионного приближения, может быть записано следующим образом:
|
1 |
|
|
rD (r,z) |
|
|
|
(r,z) |
|
D (r,z) |
|
|
|
(r,z) |
|
(r,z) |
|||
|
|
|
|
|
|
|
|
||||||||||||
|
r r |
k |
|
r |
k |
|
z k |
|
z |
k |
|
cfd,k |
|
||||||
|
|
|
|
|
1 |
|
|
|
|
|
|
|
k 1 |
|
|
|
|
||
|
|
|
|
χk νΣm(r,z)φm(r,z) Σm k (r,z)φm(r,z). (5.7) |
|||||||||||||||
|
|
|
kэф |
||||||||||||||||
|
|
|
|
|
m |
|
|
|
|
m 1 |
|
|
|
|
Предположим, что состав активной зоны не меняется в аксиальном направлении. Тогда функция плотности потока нейтронов может быть представлена в следующем виде:
k (r,z) k (r)cos |
z |
|
|
, |
|
|
||
|
Hа.з 2 эф,z |
где Hа.з – высота активной зоны, ЭФ,,Z – эффективная добавка в аксиальном направлении.
Подставляя k (r,z) в уравнение (5.7), получим:
1r ddr rDk (r) d dkr(r) cfd,k 2zDk (r) k (r)
|
1 |
k 1 |
|
|
|
k f ,m(r) m(r) m k (r) m(r), |
(5.8) |
||
kэф |
||||
|
m 1 |
|
||
|
|
131 |
|
где ω2 |
|
|
|
π |
|
2 |
|
|
|
|
|
– баклинг (или лапласиан), описывающий |
|
H |
|
2δ |
|
|||
z |
|
а.з |
|
|
||
|
|
|
|
эф,z |
аксиальную утечку нейтронов из активной зоны.
Можно показать, что уравнение (5.7) может быть обобщено на случай произвольной одномерной геометрии (плоскость, цилиндр, сфера) и записано следующим образом:
|
1 d |
|
|
d k (r) |
2 |
|
||||||
|
|
|
|
|
|
r |
|
Dk (r) |
|
|
cfd,k (r) zDk (r) k (r) |
|
r |
|
|
|
|
dr |
|
||||||
|
|
|
dr |
|
|
|
k 1 |
|
||||
|
|
|
|
|
1 |
|
|
|
|
|
||
|
|
|
|
k f ,m(r) m(r) m k (r) m(r), |
(5.9) |
|||||||
|
|
|
kэф |
|||||||||
|
|
|
|
|
|
|
|
|
m 1 |
|
где – показатель геометрии, равный нулю для плоскости, единице
для цилиндра и двум для сферы; |
2z – баклинг, описывающий |
||||
утечку |
нейтронов. |
Для |
плоскости |
и |
цилиндра |
2z [ /(Hа.з 2 эф,z)]2 |
, а для сферы 2z = 0. В более общем случае |
значения аксиального баклинга 2z могут быть взяты из других,
например, двумерных расчетов и введены в программу как исходные данные.
Физический смысл членов, входящих в уравнение (5.8), вполне объясним. Первый член левой части описывает перенос нейтронов в радиальном направлении, второй член – поглощение, замедление и утечку нейтронов в аксиальном направлении. Два члена правой части уравнения отвечают за генерацию нейтронов деления и “приход” нейтронов в рассматриваемую энергетическую группу из всех предшествующих за счет замедления.
Решение уравнения (5.9) осуществляется методом итераций источников, в соответствии с которым член, описывающий генерацию нейтронов деления, рассчитывается на нейтронных потоках предыдущей итерации:
|
1 d |
d (n) |
k 1 |
|||||
|
|
|
|
|
r Dk(r) |
k |
cfd,k 2z Dk(r) (kn) qk(n 1) |
m k (mn) , |
r |
|
|
dr |
dr |
||||
|
|
|
|
|
m 1 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
132 |
|
где n – номер текущей итерации источников;
qk(n 1)(r) |
1 |
k f ,m (mn 1)(r) ; |
(n 1) |
||
|
kэф |
|
kэф(n 1) f ,m(r) (mn 1)(r)r dr . m
Процесс итераций продолжается до тех пор, пока значение kэф не стабилизируется в рамках заданной точности , т.е.
|
|
kэф(n ) kэф(n 1) |
|
. |
|
|
k эф(n ) |
|
|
|
|
|
|
|
Таким образом, в рамках каждой итерации источников необхо- |
||||
димо решить следующее уравнение: |
|
|
|
|
|
1 d |
|
|
|
d k |
|
k 1 |
m k |
m qk', |
(5.10) |
||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
r |
|
Dk(r) |
|
|
k(r) k qk |
|
|
||||||
|
|
r |
|
|
dr |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
dr |
|
m 1 |
|
|
|
|
|||||
где |
k |
(r) |
cfd,k |
(r) 2 |
D (r); т.е. найти |
|
k |
(r) |
при известной пра- |
|||||||||||
|
|
|
|
|
z |
|
k |
|
|
|
|
|
|
|
вой части. Правая часть известна, поскольку член источника нейтронов деления определяется на нейтронных потоках предыдущей итерации, а член прихода нейтронов в группу за счет замедления определяется на только что рассчитанных нейтронных потоках текущей итерации для предшествующих энергетических групп.
Первым шагом к решению уравнения (5.10) является получение его конечно-разностного аналога, т.е. замена непрерывной пространственной зависимости k (r) значениями этой функции в дискретных пространственных точках. Для этого в каждой зоне системы вводится равномерная (т.е. с постоянным шагом) сетка пространственных точек. На границах раздела зон размещаются две точки с одинаковой координатой: одна слева от границы, другая – справа от нее, т.е. эти точки принадлежат разным зонам системы. Затем уравнение (5.10) интегрируется в окрестности произвольной точки l, т.е. в пределах интервала [rl 0,5 ri, rl 0,5 ri], где rl –
133
координата точки, а ri шаг между точками в зоне i. Такое инте-
грирование приводит к следующим результатам:
|
1 d |
|
|
|
d |
k |
|
|
|
|
|
|
|
|
d |
k |
|
|
|
|
|
d |
k |
|
|
|
||
|
|
|
|
r |
|
D |
|
|
r |
|
dr r |
|
D |
|
|
r |
|
D |
|
|
|
|
||||||
r dr |
|
dr |
|
|
dr |
|
|
|
|
|||||||||||||||||||
|
|
k |
|
|
|
|
|
|
k |
l 0,5 |
|
|
k dr |
l 0,5 |
|
|
||||||||||||
(r 0,5 r ) D |
|
|
|
k,l 1 k,l |
(r 0,5 r ) D |
|
k,l |
k,l 1 |
; |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
l |
|
|
i |
k,i |
|
|
|
ri |
|
|
|
l |
|
|
i |
k,i |
|
|
ri |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k (r) k (r)r dr k,i k,l rl ri ;
qk (r)r dr qk,l rl ri .
Окончательно конечно-разностный аналог уравнения (5.10) можно записать в следующем 3 точечном виде (опуская индекс энергетической группы k):
|
|
|
|
|
|
Al l 1 Bl l |
Cl l 1 Ql , |
|
|
(5.11) |
||||||
где |
A |
Di |
(r 0,5 r ) ; |
C |
Di |
(r 0,5 r ) ; |
|
|
|
|
||||||
r |
|
|
|
|
||||||||||||
|
l |
r |
l |
i |
l |
|
l |
i |
|
|
|
|
||||
|
|
|
i |
|
|
|
|
|
i |
|
|
|
|
|
|
|
B |
r |
r |
|
Di |
[(r 0,5 r ) |
(r 0,5 r ) ]; |
Q |
q |
r |
r . |
||||||
r |
||||||||||||||||
l |
|
i l |
i |
|
l |
i |
|
l |
|
i |
l |
|
l l |
i |
||
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
Итак, дифференциальное уравнение (5.10) преобразовалось в систему алгебраических уравнений (5.11), число которых равно количеству точек конечно-разностной сетки. Решение этой системы находится так называемым методом “двойной прогонки”, суть которого заключается в следующем.
Предположим, что значение плотности потока нейтронов в точке l может быть выражено через ее значение в следующей точке l + 1 следующим образом:
l |
l l 1 l . |
(5.12) |
134
Подставив это выражение в уравнение (5.11), получим систему рекуррентных уравнений для нахождения коэффициентов l и l
через их значения в предыдущей точке l 1 и l 1 :
Al ( l 1 l l 1) Bl l Cl l 1 Ql; (Al l 1 Bl ) l Cl l 1 (Ql Al l 1);
l |
|
|
|
Cl |
|
|
; |
(5.13) |
|||
A |
l 1 |
B |
|
||||||||
|
|
|
|
l |
|
|
l |
|
|||
|
l |
|
Ql Al l 1 |
. |
|
|
|||||
|
|
|
|||||||||
|
|
A |
l 1 |
B |
|
|
|||||
|
|
|
l |
|
l |
|
|
Для того чтобы начать вычисления коэффициентов l, l во всех пространственных точках системы, необходимо определить их значения в первой точке (значения 1, 1), а затем использовать рекуррентные соотношения (5.13).
Для вычисления значений 1 и 1 применяется краевое условие уравнения (5.4) в центре системы. В большинстве случаев расчета реакторов это условие симметрии нейтронного потока, т.е.
r k (r 0) 0.
Интегрируя уравнение (5.10) в окрестности первой точки (она размещена на расстоянии полшага от центра системы), можно получить соотношение:
[ 1 (0,5) r12 D1] 1 q1 (0,5) r12 D1 2 ,
т.е.
1 |
|
D1 |
|
|
|
; |
||
D (0,5) r2 |
||||||||
|
|
|
1 |
1 |
|
1 |
|
|
|
|
|
q (0,5) |
r2 |
||||
|
|
|
|
1 |
|
1 |
. |
|
|
|
(0,5) r2 |
|
|||||
1 |
|
|
D |
|||||
|
|
|
1 |
|
1 |
1 |
|
|
|
|
|
|
135 |
|
|
|
|
Теперь, используя рекуррентные соотношения (5.13), можно определить значения коэффициентов l, l во всех точках первой зоны, кроме последней точки. Последняя точка первой зоны n и первая точка второй зоны n + 1 имеют одну и ту же координату rn rn 1, но расположены по разные стороны границы раздела зон.
Для определения значений l, l в последней точке первой зоны и в первой точке второй зоны используются условия непрерывности плотности потока и тока нейтронов на границе раздела зон:
(rn) (rn 1);
D1 r (rn) D2 r (rn 1).
Из условия непрерывности плотности потока нейтронов следует,
что n 1; n 0.
Условие непрерывности тока нейтронов позволяет определить значения n 1 и n 1 следующим образом. Уравнение (5.10) необходимо проинтегрировать сначала слева от границы зоны, т.е. в
интервале [rn 0,5 r1, rn], а потом справа от границы зоны, |
т.е. в |
||||||||||||||||||||||||||||||||||||
интервале [rn 1, |
rn 1 0,5 r2], и результаты сложить. |
|
|
|
|
|
|||||||||||||||||||||||||||||||
Первое интегрирование даст |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
r |
D |
|
d |
|
r |
|
|
|
|
D |
|
n n 1 |
|
0,5 |
r |
|
r |
0,5q |
r |
|
r . |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
n |
1 |
|
|
n 0,5 |
|
|
1 |
|
|
r1 |
|
|
|
1 n |
n |
|
1 |
|
|
|
n n |
|
1 |
|
||||||||||||
|
|
|
dr n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Второе интегрирование даст |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
r |
|
D |
n 2 n 1 |
r |
|
|
D |
d |
0,5 |
r |
|
|
r |
|
0,5q |
|
r |
|
r . |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
n 1,5 |
2 |
|
|
|
|
|
|
|
n 1 |
2 |
dr |
|
|
|
2 n 1 n 1 |
2 |
|
|
n 1 |
n 1 |
|
2 |
|||||||||||||||
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
После сложения результатов интегрирования получим: |
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
r |
|
|
|
|
|
|
n |
|
n 1 |
|
r |
|
|
|
|
n 1 |
|
n 2 |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
2 |
|
|
n 0,5 |
|
|
D |
|
|
|
|
|
|
2 |
|
n 1,5 |
|
D |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
r |
|
|
|
|
1 |
|
|
|
r |
|
|
|
|
|
r |
|
|
2 |
|
|
r |
|
|
|
|
|
|
||||||
|
|
|
n |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
n |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
136 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
( 1 r1 2 r2) n qn r1 qn 1 r2 .
Учитывая, |
что n 1 N, |
а т 1 т 1 т т 1 т 1 т 1 т 1 , |
|||||||||||||||||||||||||||||||||||
можно связать значения n 1 |
|
и n 2 следующим уравнением: |
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
r |
|
|
|
|
D |
|
|
|
|
|
|
|
|
r |
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
n 0,5 |
|
|
|
|
|
|
|
|
|
n 1,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
n 1 |
2 |
|
|
|
|
|
|
1 |
(1 |
n 1 |
) 2 |
|
|
|
|
|
|
2 |
r |
2 |
r |
|
||||||||||||||
r |
|
|
|
|
|
|
r |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
r |
1 1 |
|
|
|
|
2 |
|
|
||||||||||
|
|
|
|
n |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
n |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
D |
|
|
|
|
|
||||
|
|
|
|
|
|
n 1,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n 0,5 |
|
|
|
|
|
|
||||||||||
|
n 2 |
2 |
|
|
|
|
|
|
|
2 |
|
(q |
|
r q |
n 1 |
r |
|
) 2 |
|
|
|
|
2 |
|
n 1 |
. |
|||||||||||
|
r |
|
|
r |
|
|
r |
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
n |
1 |
|
|
2 |
|
|
|
|
|
r |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
n |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
2 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из последнего уравнения легко найти значения коэффициентовl , l в первой точке второй зоны:
n 1
n 1
|
|
|
|
|
|
|
|
|
|
|
|
|
rn 1,5 |
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
rn |
|
|
|
r2 |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
|||||||
|
|
rn 0,5 |
|
|
D |
|
|
|
|
|
|
|
|
rn 1,5 |
|
|
|
D |
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
2 |
|
|
|
|
|
1 |
|
(1 |
n 1 |
) 2 |
|
|
|
|
|
|
|
|
2 |
|
r |
r |
|||||||||||||||
r |
|
|
|
r |
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
r |
1 |
1 |
2 |
2 (5.14) |
|||||||||||||||||
|
|
n |
1 |
|
|
|
|
|
|
|
|
n |
2 |
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
rn 0,5 |
|
|
|
D |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
q |
r q |
|
|
|
|
r |
2 |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r1 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
n |
1 |
|
n 1 |
|
2 |
|
|
|
rn |
|
|
|
|
|
|
. |
|||||||||||||
|
|
rn 0,5 |
|
|
D |
|
|
|
|
|
|
|
|
rn 1,5 |
|
|
|
D |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
2 |
|
|
|
|
|
|
1 |
(1 |
n 1 |
) 2 |
|
|
|
|
|
|
|
|
|
|
2 |
r |
2 |
r |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
rn |
|
r1 |
|
|
|
|
|
|
|
rn |
r2 |
1 |
1 |
2 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Затем можно снова воспользоваться рекуррентными формулами (5.13) для нахождения значений коэффициентов l , l во всех точках второй зоны, кроме последней ее точки. В последней точке этой зоны коэффициенты 1, 0 , а в первой точке следующей
зоны надо воспользоваться соотношениями (5.14). В результате использования такого алгоритма в каждой зоне можно будет найти значения коэффициентов l , l для всех пространственных точек системы.
Таким образом, так называемой “прямой” прогонкой, т.е. в направлении от центра к периферии, вычисляются коэффициенты
137
l , l для уравнения (5.12). Но для того чтобы воспользоваться этим уравнением, необходимо “зацепиться” за значение плотности потока нейтронов в последней точке системы и рассчитать значения l во всех предыдущих точках уже в направлении от периферии к центру, т.е. осуществить “обратную” прогонку.
Значение плотности потока нейтронов в последней точке может быть найдено из краевого условия задачи. Например, в большинстве расчетов используется нулевое граничное условие, т.е.N 0, где N – номер последней пространственной точки, находящейся на экстраполированной границе системы. Тогда, используя уравнение (5.12), можно найти значения l во всех предшествующих точках. Однако программа TIME26 делает возможным использование и другого краевого условия, а именно условия симметрии на границе реактора, т.е. r (rN ) 0. Для нахождения зна-
чения N необходимо проинтегрировать уравнение (5.10) в окрестности последней точки. В этом случае она отстоит на полшага от границы системы, т.е. интегрирование производится в интер-
вале [rN 0,5 ri, rN 0,5 ri]. |
Интегрирование |
|
дает следующий |
|||||||||||||||||||||||
результат: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
r |
|
Di |
|
|
|
|
|
|
|
i N |
rN |
ri qNrN |
ri . |
|
|||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
dr N 0,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
d (rN 0,5) |
|
N N 1 |
|
N (1 N 1) N 1 |
|
||||||||||||||||||||
Поскольку |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, значе- |
|
|
dr |
|
|
|
|
|
ri |
|
|
|
|
|
|
|
|
ri |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
ние N может быть вычислено по следующей формуле: |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
rN 0,5 |
|
D |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
q |
N |
|
|
|
|
|
|
|
|
i |
|
N 1 |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
r |
|
|
|
|
r2 |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
N |
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
i |
|
|
. |
|
||||
|
|
rN 0,5 |
|
D |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
i |
|
(1 |
N 1 |
) |
i |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
r |
|
|
r |
2 |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
i |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
138 |
|
|
|
|
|
|
|
|
|
|
|
Итак, “прямой” прогонкой, т.е. в направлении от центра к периферии, находятся значения коэффициентов l и l во всех пространственных точках системы. Затем из краевых условий определяется значение плотности потока нейтронов в последней точке системы, и “обратной” прогонкой, т.е. в направлении от периферии к центру, вычисляются значения l во всех точках. Этот алгоритм и получил название метода “двойной” прогонки.
Ключевым моментом в решении уравнения (5.5), описывающего временное поведение концентраций изотопов (t) при работе реактора на постоянной мощности, является вычисление экспоненци-
альной матрицы exp(Aˆ t). Постоянство мощности системы в пределах временного интервала t означает то, что элементы матрицы
Aˆ также не зависят от времени.
Для большинства цепочек изотопных переходов, имеющих ме-
сто при работе ядерного реактора на мощности, матрица Aˆ может быть записана в нижне-треугольном виде, т.е. ее элементы ai, j 0 при i j. Хотя программа TIME26 способна решать урав-
нения выгорания и для произвольной матрицы, здесь будет описан алгоритм вычисления экспоненты от нижне-треугольной матри-
цы Aˆ .
Известно, что матрицу Aˆ можно представить в следующем ви-
де:
Aˆ Dˆ diag( ) Dˆ 1,
где Dˆ – матрица, столбцами которой являются собственные векто-
ра матрицы Aˆ ; diag( ) – диагональная матрица, элементами кото-
рой являются собственные значения матрицы Aˆ .
Нетрудно показать, что, если известны собственные векторы и
собственные значения матрицы Aˆ , то экспонента от нее может быть вычислена следующим образом:
139
|
|
|
|
exp(Aˆ) Dˆ diag[exp( )] Dˆ 1 . |
|
|
(5.15) |
||||||||||
Действительно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
ˆ |
|
ˆ |
ˆ |
|
Aˆ2 |
|
|
Aˆn |
|
|
ˆ ˆ ˆ 1 |
ˆ |
ˆ |
1 |
|
|
exp(A) E |
A |
|
|
... |
|
|
DED |
|
Ddiag( )D |
|
|
||||||
2! |
n! |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
2 |
|
|
|
|
|
|
|
n |
|
|
|
|
|
||
Dˆ |
diag |
|
Dˆ |
1 |
... Dˆ |
diag |
|
|
Dˆ 1 |
Dˆ diag[exp( )]Dˆ 1 . |
|||||||
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2! |
|
|
|
|
|
n! |
|
|
|
|
|
В случае нижне-треугольной матрицы определение ее собственных значений и собственных векторов существенно упрощается. Легко показать, что собственными значениями нижне-треугольной матрицы являются ее диагональные элементы ai,i , а элементы мат-
риц Dˆ и Dˆ 1 можно вычислить по следующим формулам:
|
|
|
|
|
i 1 |
|
|
|
|
|
||
(Dˆ)i, j di, j |
|
ai,l dl, j |
|
|||||||||
|
a |
|
a |
|
; di,i 1; |
|||||||
|
|
|
|
l j |
|
|
|
|
|
|||
|
|
|
|
|
|
|
j, j |
i,i |
|
|||
|
|
|
|
|
|
|
i |
|
|
|
|
|
(Dˆ |
1 |
)i, j di, j |
|
|
|
al, j |
di,l |
1. |
||||
|
a |
a |
|
|
; di,i |
|||||||
|
|
|
|
|
l j 1 |
|
|
|
|
|||
|
|
|
|
|
|
|
i,i |
|
j, j |
|
После определения экспоненты от матрицы Aˆ t новый изотопный состав может быть определен по формуле
(t0 t) exp(Aˆ t) (t0).
В программе TIME26 имеется возможность оценить чувствительность эффективного коэффициента размножения нейтронов kэф к изменению концентрации изотопа l в зоне i. Для этого рассчиты-
140