О распределении нейтронов в слабо размножающих средах
Макроскопическое сечение утечки.
Возьмем для примера сферический реактор.
Если известен коэффициент k, то шар может быть реактором только в случае k>1. При выполнении этого условия найдем положительное значение В2 и экстраполированный размер шара RЭ = π/В, необходимый для осуществления в шаре стационарной цепной ядерной реакции. Реальный размер критического реактора будет несколько меньше, чем RЭ. Значит, если из того же материала сделать шар радиусом RЭ, то получится реактор с некоторым превышением Кэф над единицей (с положительной реактивностью).
Геометрический параметр:
B2 = (π/Rэ)2 (10.9)
и К0 - значение коэффициента k, необходимое для осуществления стационарной цепной ядерной реакции. Значение К0 найдем по формуле
Для вероятности избежать утечки тепловых нейтронов можно записать:
Анализируя эту формулу можно отметить, что величина DB2 играет роль «макроскопического сечения утечки».
Из тепловой группы нейтроны выбывают только двумя путями: или поглощаются, или «утекают». Следовательно, отношение Σа/(Σа+DB2) есть доля поглощаемых нейтронов по отношению ко всем выбывающим из группы нейтронам или, говоря иначе, вероятность избежания утечки во время пребывания в тепловой группе.
При этом отношение Σf/Σа есть вероятность для всех поглощаемых нейтронов быть захваченными в делящемся нуклиде с последующим делением ядра. Вероятность поглощения с делением ядра для всех нейтронов группы будет равна произведению этих двух вероятностей, т.е. отношению Σf/(Σа+DB2).
В критическом реакторе В2 = χ2 или Вг2 = Вм2
Постоянная величина В2 получила в английской литературе название «баклинг» (buckling), что наиболее близко по смыслу слову «выгиб» (функции
Ψ(r), если изобразить ее графиком).
В случае, когда баклинг - положительное число (B2 может быть и отрицательным, т.е. величина В - мнимой), ограниченным в центре решением уравнения распределения Ф в сферическом реакторе будет функция
Ф(r)=С r-1 Sin(Br),
которая подходит для описания плотности потока нейтронов в шаре. График этой функции «выгнут» вверх, т.е. будучи положительной в точке r=0, она убывает при отклонении r от нуля и при r=π/В обращается в нуль.
Слабо размножающей назовем среду, в которой k - 1 < 0. В такой среде стационарная плотность потока нейтронов, функция Ф(r), может существовать только с отрицательным баклингом (т.е. на графике с «выгибом» вниз).
В одномерной сферической геометрии примером неразмножающей или слабо размножающей среды может быть отражатель шарообразного реактора или вставка (для какой-либо специальной цели) с коэффициентом k < 1 в виде сплошного шара в центре реактора.
В неразмножающем отражателе (т.е. где k = 0) будем иметь
(10.12)
Вспомним, что величина L2 =D/Σa называется квадратом длины диффузии. Ее неформальный смысл раскрывается равенством
где r2 - средний квадрат смещения (по прямой) теплового нейтрона за время диффузии. Соответственно величину L называют диффузионной длиной. Для сравнения диффузионные длины в реакторных материалах-замедлителях приведены в табл. 10.1.
Таблица 10.1 Диффузионные характеристики замедлителей
Замедлитель | Плотность, г/см3 | L, см | Σa, см-1 | D, см |
Вода (Н2О) | 1.00 | 2.88 | 0.017 | 0.142 |
Тяжелая вода (D2O) | 1.1 | 0.00008 | 0.80 | |
Графит | 1.62 | 50.2 | 0.00036 | 0.903 |
Если считать для простоты, что толщина отражателя в несколько раз больше чем L, то плотность потока нейтронов в сферическом отражателе будет функцией такого вида:
(10.14)
Коэффициент С' (точнее, отношение С'/C) определяется условием сопряжения функций в активной зоне и в отражателе на границе Ra3 между ними. Характерный общий вид функции Ф(r) в реакторе с отражателем показан на рисунке (10.1).
Рис. 10.1. Профиль плотности потока нейтронов в реакторе с отражателем (односкоростная модель):
Rаст - радиус слабо размножающей вставки;
Rа.з. - радиус активной зоны (граница с отражателем);
Rэ - экстраполированный радиус «голого» реактора;
Rотр- внешний радиус отражателя; Ф(r) - активная зона однородна;
Ф '(r) - в активной зоне вставка.
Величина DB2 как «сечение утечки» в отражателе отрицательна. А это значит, что в балансе элементарного слоя отражателя формальный член утечки нейтронов, по существу, будет прибылью (за счет диффузии частиц со стороны активной зоны реактора).
Это в общем верно и для слабо размножающей вставки в центре реактора. Там токнейтронов направлен от границы к центру, а в самом центре обращается в нуль (т.е. во вставке I'(r)≤0). Если в центральной вставке K≥0 (но k<1), то плотность потока нейтронов во вставке будет иметь вид функции
(10.16)