§5.2. Основные свойства вынужденного деления

Рассмотрим основные свойства вынужденного деления ядер нейтронами, когда энергия возбуждения составного ядра, образованного при захвате нейтрона, превышает энергетический барьер деления, т.е.

,

(5.2.1)

а согласно (4.5.32)

.

(5.2.2)

Из (5.2.1) и (5.2.2) следует, что эффективное протекание реакции деления составного ядра возможно тогда, когда кинетическая энергия нейтрона

.

(5.2.3)

Выполнение неравенства (5.2.3) возможно в двух случаях.

1). Если S_n(C) > W_f, то из (5.2.3) следует, что реакция деления составного ядра не имеет энергетического порога и деление возможно при любой энергии налетающего нейтрона. Нуклиды, захват нейтрона которыми ведет к образованию составного ядра, обладающего такими свойствами, называютсяделящимися или топливными нуклидами. Топливные нуклиды используются для производства энергии в ядерных реакторах и в ядерном оружии.

2).W_f>S_n(C). И хотя и в этом случае деление составного ядра является экзоэнергетическим процессом и формально не имеет энергетического порога (ядро может разделиться подбарьерным путем, но вероятность такого процесса крайне мала), однако эффективно они могут делиться только тогда, когда кинетическая энергия нейтронов

.

(5.2.4)

Нуклиды, образующие при захвате нейтрона составные ядра, для которых выполняется это условие, называются сырьевыминуклидами, так как из них возможно получение топливных нуклидов.

В таблице 5.2.1 приведены характеристики некоторых тяжелых нуклидов, способных делиться под действием нейтронов.

Таблица 5.2.1

Нуклид	232Th	233U	235U	238U	239Pu
Составное ядро	233Th	234U	236U	239U	240Pu
Энергетический барьер Wf,МэВ	5,9	5,5	5,75	5,85	5,5
Энергия связи нейтрона Sn,МэВ	5,07	6,77	6,4	4,76	6,38

Из этой таблицы видно, что нуклиды ²³³U,²³⁵Uи²³⁹Puявляются топливными нуклидами, а²³²Thи²³⁸U– сырьевыми. Переработка сырьевых нуклидов в топливные основана на реакции радиационного захвата нейтрона (4.9.11):

n + 232Th → γ + 233Th ,	(5.2.5)
n + 238U → γ + .	(5.2.6)

В реакции (5.2.6) были получены первые трансурановые элементы.

На рис.5.2.1 приведены зависимости сечения деления от кинетической энергии нейтронов, показывающие существенное различие топливных и сырьевых нуклидов относительно вынужденного деления. В то время как топливные нуклиды в тепловой области имеют сечения ~ 1000барн, а сырьевые вообще не делятся, то в области быстрых нейтронов сечение деления сырьевых нуклидов быстро возрастает, образуя своеобразный порог в районе ~1,4МэВ, и сравнивается по порядку величины с сечением деления топливных нуклидов.Выясним, чем вызвано такое различие на примере нуклидов²³⁵Uи²³⁸U. Во-первых, составное ядро²³⁶Uимеет большее значение параметра делимости, нежели составное ядро²³⁹Uи, следовательно, меньшую высоту потенциального барьераW_f (см. таблицу 5.2.1). Во-вторых, энергия связи нейтрона в четно-четном ядре²³⁶Uбольше, чем в четно-нечетном ядре²³⁹U,согласно пятому члену формулы Вейцзеккера (2.1.1). Другими словами, энергия связи парного нейтрона всегда больше, чем непарного.

Измерения кинетической энергии осколков показывают, что деление носит асимметричный характер, а образующиеся осколки имеют различные величины кинетической энергии и, следовательно, различные массы. Распределение осколков по энергиям для случая деления²³⁵Uтепловыми нейтронами представлено на рис. 5.2.2, из которого следует распределение осколков по массам (кривая «а» на рис. 5.2.3). Действительно, если принять начальный импульс системы ядро + нейтрон, по величине,Р₁=Р₂, откуда следует, что равным нулю, то импульсы осколков должны быть равны друг другу

.

(5.2.7)

При делении образуется несколько десятков пар осколков преимущественно неравной массы. Наиболее вероятным (~ 6 ÷ 7%, площадь под кривой на рис. 5.2.3. нормирована на 200%) оказывается выход Yосколков с массовыми числами 95 и 141, т.е. массы осколков относятся как 2 : 3. Вероятность симметричного деления в 600 раз меньше. С ростом энергии нейтронов асимметрия в распределении масс осколков уменьшается (кривая «б» на рис. 5.2.3).

Объяснение асимметрии деления при помощи капельной модели предполагает деформацию делящегося ядра в виде груши, но деление капли на две равные части оказывается наиболее энергетически выгодным. Одно из возможных объяснений асимметрии деления может быть получено с привлечением модели ядерных оболочек, как результат преимущественного образования ядер-осколков с близкими к магическим (50 и 82) числами заполнения протонных и нейтронных оболочек.

В процессе деления выделяется энергия примерно равная 200 МэВ. Подавляющая часть этой энергии приходится на кинетическую энергиюQ_fкосколков, приобретаемую ими в результате кулоновского отталкивания. Энергия кулоновского взаимодействия осколков в момент их образования (позиция 4 на рис. 5.1.3) определяется кулоновским барьером (1.9.2) и составляет

,

(5.2.8)

где: Z₁иZ₂–заряды осколков, аR₁иR₂- их радиусы. Подсчет энергии по этой формуле для пары наиболее вероятных осколков дает величину ~ 170МэВ.

Образующиеся осколки должны быть радиоактивными и могут испускать нейтроны. В различных областях периодической системы элементов наблюдаются следующие соотношения между числом нейтронов и протонов для β-устойчивых ядер, лежащих на дорожке стабильности (рис. 1.1.2):

Ядро
N/Z	1,0	1,3	1,45	1,55

Из приведенных данных следует, что осколки перегружены нейтронами, так как они образуются из ядер ²³⁶U, у которыхN/Z≈ 1,55. Поэтому они располагаются ниже дорожки стабильности и должны быть β^--активными. Из-за большого избытка нейтронов выход на дорожку стабильности должен осуществляться путем нескольких последовательных β^--распадов, т.е. осколки деления должны давать начало длинным цепочкам радиоактивных ядер. Пример распада наиболее вероятных осколков приведен на рис. 3.5.4.

Таким образом, часть энергии деления освобождается в виде энергии Q_fβпри выходе на дорожку стабильности путем β^--распадов.

Кроме того, осколки будут испускать нейтроны непосредственно после деления, когда они находятся на расстоянии ~ 10^-8см. что соответствует ~ 10^-17спосле захвата нейтрона ядром. Этивторичные нейтроны, возникающие в результате деления, в отличие отпервичных нейтронов, вызывающих деление, называютсямгновенными нейтронами деления, и также уносят часть энергииQ_nf. Впервые вторичные нейтроны деления наблюдались Ж.Кюри в 1939 г. при делении²³⁵Uтепловыми нейтронами, а в 1941 г. Э.Ферми были зарегистрированы нейтроны, испускаемые при спонтанном делении ядер²³⁸U. Измерения показали, что на один акт деления возникает от одного до пяти вторичных нейтронов.

Некоторая часть вторичных нейтронов выделяется из осколков деления спустя время от 0,1 до 50 с. Эти нейтроны были названызапаздывающими. Физическая причина появления запаздывающих нейтронов обсуждается в §3.5.

Полное число вторичных нейтронов есть сумма числа (p–prompt) мгновенных нейтронов и числа(d–delay) запаздывающих нейтронов:

,

(5.2.9)

Доля запаздывающих нейтронов в полном числе вторичных нейтронов деления

,

(5.2.10)

Несмотря на незначительную долю, запаздывающие нейтроны играют определяющую роль в управлении цепной реакции деления в ядерном реакторе.

Многочисленные экспериментальные исследования дают основания полагать, что величина почти линейно увеличивается с ростом энергии первичных нейтронов

,

(5.2.11)

где - средний выход мгновенных нейтронов при тепловой энергии первичных нейтронов, а коэффициента≈ (0,1 ÷ 0,15)МэВ^-1для большинства ядер. Данные по числу вторичных нейтронов деления для ряда нуклидов приведены в таблице5.2.2.

П

Таблица 5.2.2.

Данные для топливных нуклидов – деление тепловыми нейтронами,

для сырьевых – деление нейтронами с энергией 2 МэВ.

Нуклид	233U	235U	239Pu	232Th	238U
	2,480	2,407	2,874	2,37	2,75
	0,0066	0, 016	0,0061	0,051	0,041
	2,47	2,42	2,88	2,42	2,79
	0,0026	0,0065	0,0021	0,021	0,015
a, МэВ-1	0,123	0,133	0,135	-	-

оскольку в разных актах деления случайным образом возникают различные по массам пары осколков и различное число нейтронов, то кинетическая энергия вторичных нейтронов будет так же случайной величиной, но при наблюдении за большим количеством делящихся ядер будет наблюдаться уже вполне закономерныйэнергетический спектр вторичных нейтронов. Измерения распределения вторичных нейтронов по энергиям позволило установить, что форма распределения очень слабо зависит от кинетической энергии первичных нейтронов и сорта делящихся ядер и может быть удовлетворительно представлена максвелловским распределением. В случае деления²³⁵U тепловыми нейтронами спектр вторичных нейтронов имеет вид (см. рис 5.2.4):

.

(5.2.12)

В настоящее время измерение спектра доведены до 17 МэВ. Средняя энергия вторичных нейтронов при этом близка к 2МэВ, а наиболее вероятная (в максимуме распределения) составляет около 0,7МэВ. Так как нейтроны уносят энергию, то при деления²³⁵Uтепловыми нейтронамиМэВ.

При изучении процесса деления были обнаружены мгновенные γ-кванты, испускаемые возбужденными осколками за время ~ 10^-14спосле вылета из них мгновенных нейтронов. Энергетический спектр γ-излучения непрерывный и убывает с ростом энергии, максимальная энергия γ-квантов составляет около 7МэВ. В процессе деления²³⁵Uтепловыми нейтронами на один акт деления возникает в среднем примерно 7 γ-квантов со средней энергией около 1МэВ, а среднее количество энергии, уносимой γ‑квантами составляет 7МэВ.

Примерное распределение высвобождаемой энергии для случая делении ядер ²³⁵Uтепловыми нейтронами приведено в таблице 5.2.3.

§

Таблица 5.2.3

Мгновенное энерговыделение	Кинетическая энергия осколков		168	МэВ
	Энергия вторичных нейтронов деления		5	“
	Энергия мгновенных γ‑квантов		7	“
Запаздывающее энерговыделение	Энергия β-частиц осколков деления		8	“
	Энергия γ‑излучения осколков деления		7	“
	Энергия антинейтрино при β-распадах		10	“
		Всего	205 МэВ

5.3. Цепная реакция деления

Возникновение вторичных нейтронов в процессе деления тяжелых ядер нейтронами дает возможность осуществить цепной процесс деления ядер. Цепной процесс характерен тем, что в его основе лежит экзоэнергетическая реакция, возбуждаемая нейтроном, которая порождает вторичные нейтроны. В этом случае появление нейтрона в делящейся среде вызывает цепь следующих друг за другом реакций деления, которая продолжается до обрыва вследствие потери нейтрона – носителя процесса. Основных причин потерь нейтрона две: поглощение первичного нейтрона ядром без испускания вторичных (например, радиационный захват) или уход нейтрона за пределы объема среды (называемыйактивной зоной), в которой протекает цепной процесс деления. Если в результате реакции возникает более одного нейтрона, которые в свою очередь вызывают деления, то такая реакция являетсяразветвленной реакцией. Осуществление цепного процесса оказывается возможным только с помощью разветвленной реакции. Средняя длина пробега нейтрона от точки рождения до точки, в которой нейтрон производит деление, является макроскопической величиной. Поэтому цепная реакция деления является макроскопическим процессом. Каждый нейтрон, участвующий в цепном процессе, проходит цикл обращения: рождается в реакции деления, некоторое время существует в свободном состоянии, затем либо теряется, либо порождает новый акт деления и дает нейтроны следующего поколения. Нейтрону необходимо, хотя и малое, но конечное время для прохождения через цикл обращения. Среднее время τ, полученное усреднением по большому числу циклов деления нейтронами, называетсявременем цикла обращения нейтронаилисреднимвременем жизни нейтронов. Таким образом, цепной процесс деления можно представит как последовательность следующих друг за другом лавин или поколений, разделенных промежутком времени τ:

N0→N1→N2→ …→Ni→ Ni+1→… ,

(5.3.1)

где N– число нейтронов в данном поколении. Отношение числа нейтронов последующего поколения к их числу в предшествующем поколении во всем объеме активной зоны называется коэффициентом размножения нейтронов:

k = Ni+1/ Ni .

(5.3.2)

Величины τ и kполностью определяют развитие цепного процесса во времени. Пустьk - постоянная величина, т.е. не зависит от времени. Тогда число нейтронов в следующем поколенииN_i+1=k N_i, затем, через промежуток времени τ количество нейтроновN_i+2=k N_i+1=k² N_i, через время 2 τ количество нейтронов составитN_i+3=k N_i+2=k² N_i+1=k³ N_iи т.д. Количество нейтронов в поколении под номеромm(число нейтронных циклов) составит

Nm=N0km,

(5.3.3)

если число нейтронов в начальный момент времени наблюдения равно N₀. Время наблюдения при этом составитt = mτ, что позволяет записать зависимость (5.3.3) в явном виде от времени:

N(t) =N0kt/τ.

(5.3.4)

Однако, выражения (5.3.3) и (5.3.4) верны только приблизительно, поскольку случаи рождения и исчезновения нейтронов происходят случайным образом, и в любой момент времени в активной зоне присутствуют нейтроны из разных поколений, т.е. процесс изменения числа нейтронов в активной зоне происходит непрерывно.

Приращение числа нейтронов в цепном процессе за время τ цикла обращения нейтрона составит:

ΔN=Nk – N= N(k – 1),

(5.3.5)

а скорость изменения числа нейтронов будет равна

.

(5.3.6)

Уравнение (5.3.6) называется точечным уравнением кинетикибез запаздывающих нейтронов. Считаяk постоянной величиной и разделяя в (5.3.6) переменные, получаем решение этого уравнения:

,

(5.3.7)

где N₀=N(t= 0) - число нейтронов в начальный момент наблюдения.

Проанализируем полученное выражение (5.3.7).

Если k > 1, то число нейтронов в активной зоне будет непрерывно увеличиваться и процесс цепной реакции, раз возникнув, будет сам собой развиваться во времени. Процесс сk > 1 носит названиенадкритического режима.

При k = 1 количество нейтронов в активной зоне и число происходящих в единицу времени делений не изменяются со временем и остаются постоянными. Такой режим носит названиекритического режима.

Наконец, если k < 1, то процесс размножения нейтронов затухает и называется соответственноподкритическим режимом.

Таким образом, для протекания самоподдерживаемой цепной реакции деления необходимо, чтобы k ≥ 1. Возможность осуществления цепной реакции обычно характеризует коэффициентk_∞размножения в среде с бесконечным объемом, когда можно пренебречь утечкой нейтронов через поверхность активной зоны. Тогда для активной зоны конечных размеров

k =κk∞,

(5.3.8)

где κ – вероятность нейтрону избежать утечки из активной зоны конечного объема. Если существует некоторый конечный объем, то конфигурация, нуклидный состав и масса активной зоны, при которых выполняется условие

k = κk∞ ≥ 1,

(5.3.9)

называются критическими параметрами. Величинаkзависит от многих параметров: нуклидного состава активной зоны, ее формы и размера, от энергетического спектра нейтронов, вызывающих деление. Расчет величиныkявляется сложной инженерно-физической задачей и требует знания огромного числа констант, определяющих протекание цепного процесса в реальных активных зонах.

Рассмотрим цепной процесс, у которого время нейтронного цикла τ ≈ 10^-3с. Такая величина τ характерна для реакторов на тепловых нейтронах. Если предположить, что коэффициент размноженияk = 1,005, то за одну секунду число нейтронов, согласно (5.3.7) увеличится в

(5.3.10)

В такое же число раз возрастет количество делений в единицу времени и, следовательно, мощность установки. Ясно, что управлять таким процессом невозможно и в контролируемой установке превышение kнад единицей всего на 0,5% недопустимо. Однако, приведенная оценка не учитывает запаздывающих нейтронов и поэтому является завышенной. Действительно, число нейтронов в активной зоне в данный момент времени может быть представлено следующим образом:

Ni+1 = Np + Nd,

(5.3.11)

где: N_p– количество мгновенных нейтронов, возникших непосредственно в момент деления ядер;N_d– количество запаздывающих нейтронов, возникших из осколков деления в результате запаздывающего энерговыделения. Разделив равенство (5.3.11) слева и справа на количество нейтроновN_i в предыдущем цикле и учитывая определение (5.3.2), получим, что коэффициент размножения может быть представлен в виде суммы двух слагаемых:

k = kp + kd,

(5.3.12)

из которых первое является коэффициентом размножения на мгновенных нейтронах, а второе –коэффициентом размножения на запаздывающих. Вспомним, что доля запаздывающих нейтронов в полном числе вторичных нейтронов деления обозначаетсяβ (см. определение (5.2.9) и таблицу (5.2.2.)). Тогда в цепном процессе, идущем в²³⁵Uпод действием тепловых нейтронов

kp = k(1 – β) = 1,005(1 – 0,0065) = 0,9985,	(5.3.12)
kd = k β=1,005·0,0065=0,0653.	(5.3.12)

Таким образом, цепной процесс на одних только мгновенных нейтрона является подкритическим, и управление процессом может осуществляться с помощью изменения количества запаздывающих нейтронов. Если k_pстановится равным или больше единицы, что соответствуетk≥ (1 +β), то цепной процесс становится неконтролируемым.

Найдем среднее время τ₀нейтронного цикла с учетом запаздывающих нейтронов. По правилу нахождения среднего

τ0= (1 – β)τp+ βτd,

(5.3.14)

где τ_p– среднее время жизни мгновенных нейтронов, а τ_d– запаздывающих. Например, при делении²³⁵Uтепловыми нейтронами среднее время жизни запаздывающих нейтронов около 13с, а τ_p ≈ 10^-3си

τ0= 10-3+ 0,085 ≈ 0,085с.

(5.3.15)

Из приведенного примера следует важный вывод о том, что среднее время цикла обращения нейтронов в управляемом цепном процессе определяется средним временем жизни запаздывающих нейтронов и не зависит от времени жизни быстрых, но при условии k< (1 +β). Используя в примере (5.3.10) время τ = τ₀= 0,085сполучим, что за одну секунду мощность цепного процесса увеличится всего на 6 %, что не представляет проблем для регулирования. В теории регулирования цепного процесса обычно используется величинаT,называемая периодом реактора, которая есть время, в течение которого количество нейтронов в активной зоне увеличивается в «е» раз. Из (5.3.7) получаем, что

.

(5.3.16)

Если опять же принять k= 1,005, а τ = 0,085с, то период реактораТ= 17с. Приk → 1T→ ∞, что следует непосредственно из (5.3.16).

Рассмотрим кинетику цепного процесса с учетом запаздывающих нейтронов. Скорость приращения мгновенных нейтронов, по аналогии с (5.3.6), будет равна

,

(5.3.17)

а скорость приращения запаздывающих в соответствии с законом радиоактивного распада осколков относительно испускания запаздывающих нейтронов (λ– постоянная такого распада) равна

λс,

(5.3.18)

где с– количество накопившихся в предыдущих поколениях ядер‑источников запаздывающих нейтронов. Полная скорость изменения числа нейтронов

.

(5.3.19)

Это уравнение необходимо дополнить уравнением для скорости образования ядер предшественников:

.

(5.3.20)

Уравнения (5.3.19) и (5.3.20) образуют систему связанных линейных дифференциальных уравнений точечной кинетики с учетом запаздывающих нейтронов. При более точном рассмотрении учитывают шесть временных групп запаздывающих нейтронов и получают систему из семи уравнений.

Покажем, что на тепловых нейтронах можно организовать цепной процесс деления в металлическом уране природного состава, в котором относительное атомное содержание ²³⁸Uсоставляет 99,28 %,²³⁵U– 0,714 % и²³⁴U– 0,006 %. Тепловыми нейтронами делятся только нуклид²³⁵U. Ввиду ничтожного содержания²³⁴Uего участие в цепном процессе учитывать не будем. Среднее число η вторичных нейтронов на один акт поглощения теплового нейтрона ураном природного состава будет равно (дробь определяет вероятность нейтрону произвести деление):

,

(5.3.21)

где: - среднее число вторичных нейтронов на один акт деления;n –концентрация ядер нуклида²³⁵Uили²³⁸U(с соответствующими верхними индексами); σ_а–сечение захвата нейтронов ядрами²³⁵Uили²³⁸U;⁵σ_f – сечение деления ядер²³⁵Uнейтронами. Для тепловых нейтронов эти величины равны:= 2,41;⁵σ_f= 582,6 барн (рис. 5.2.1);⁵σ_а= 680,1барндля ядер²³⁵U;⁸σ_а= 2,71барндля ядер²³⁸U. Для природного урана⁸n/⁵n = 99,28/0,714 = 139. Подставив эти значения в формулу (5.3.21), получим η = 1,33. Следовательно, цепной процесс на ядрах²³⁵Uв составе металлического урана природного изотопного состава возможно осуществить, если в результате утечки и замедления вторичных нейтронов деления до тепловых энергий потерять в среднем не более 0,3 нейтрона.

Однако самопроизвольный цепной процесс деления в металлическом уране природного состава невозможен и вот почему. При делении ядер средняя энергия вторичных нейтронов составляет ~ 2 МэВ, т.е. превышает фактический порог деления ядер²³⁸U, равный 1,35МэВ (см. рис. 5.2.1). Поэтому целесообразно энергетический спектр нейтронов деления разделить на две энергетические группы с энергиями 0 <Т_n< 1,35МэВи 1,35 <Т_n< 6МэВ. Относительная доля нейтронов в первой группе составляетf₁ = 0,42 иf₂ = 0,55 во второй. Для нейтронов первой энергетической группы средние значения величин, входящих в формулу (5.3.21), приведены в таблице 5.3.1. Тогда среднее число вторичных нейтронов при делении ядер²³⁵U на один захваченный нейтрон из первой энергетической группы составит

Таблица 5.3.1 Нуклид ,б ,б 2,50 - 1,61 - 1,70 0,042

(5.3.22)

г

Таблица 5.3.2

Нуклид		,б	,б	,б
	2,68 2,9	1,32 0,57	1,37 0,6	1,8 2,3

де. Величины, необходимые для оценки среднего числавторичных нейтронов деления ядер²³⁵U при захвате нейтрона из второй группы, приведены в таблице 5.3.1. Следует учесть, что практически каждый из актов неупругого рассеяния приводит к тому, что энергии нейтрона становится меньше порога деления и переводят их в первую группу. Поэтому, а. Подставив необходимые величины в (5.3.22) получим, что0.0055. Теперь необходимо учесть деление ядер²³⁸Uнейтронами с энергиейТ_n> 1,35МэВ. В этом случае

.

(5.3.23)

Доля нейтронов из второй группы, равная

.

(5.3.24)

в результате неупругих взаимодействий с ядрами ²³⁵U и²³⁸U перейдет в группу с энергией 0 <Т_n< 1,35МэВи вызовет деление ядер²³⁵U. Так как, то на один захваченный нейтрон из этой группы будет в среднем рождатьсявторичных нейтронов.

Окончательно, полное число вторичных нейтронов на один захваченный составит

(5.3.25)

Таким образом, даже пренебрегая утечкой нейтронов, осуществить цепной процесс деления ядер в металлическом уране природного изотопного состава нельзя.

Осуществить самоподдерживаемую цепную реакцию деления, возбуждаемую быстрыми нейтронами, оказывается возможным, если увеличить относительную концентрацию ядер ²³⁵Uв металлическом уране до~4%.

Если обогащение нуклидом ²³⁵Uсоставляет 30% то, используя выше приведенные данные для быстрых нейтронов,.

Таким образом, если один из нейтронов необходим для поддержания цепного процесса, то еще один из нейтронов можно использовать для воспроизводства ядерного горючегов реакции (5.2.6), так как имеется запас 0,3 нейтронов для компенсации неизбежных потерь нейтронов в результате утечки и захвата ядрами конструкционных материалов.

1Исключение составляет случай I = 0, когда все три проекции на осиX, Y, Z имеют определенное значение, равное нулю.

1Изомерами называются метастабильные ядра, имеющие период полураспада, отличный от периода полураспада ядра в основном состоянии.

1 Плазмой называют газ, в котором 5 ÷ 10 % атомов или молекул находятся в ионизованном состоянии. С ростом температуры до ~ 10⁷К газ становится полностью ионизованным и состоит из ядер и электронов.

222