Решение уравнений кинетики с учетом всех групп запаздывающих нейтронов является довольно громоздким. Некоторое представление об относительной значимости установившегося и переходных членов уравнения (12.38) можно получить, рассматривая простой случай, когда имеется только одна группа запаздывающих нейтронов, характеризуемых постоянной распада λ. Величину λ определим как средневзвешенное значение по всем группам запаздывающих нейтронов. Для случая одной группы запаздывающих нейтронов уравнение (12.36) принимает вид
(12.43)
где - суммарная доля запаздывающих нейтронов; - среднее значение постоянной распада осколков деления, излучающих запаздывающие нейтроны.
Для реактора, работающего на тепловых нейтронах, в котором в качестве делящегося изотопа используется U235, суммарная доля запаздывающих нейтронов =0,0064, а постоянная распада λ=0,077 с-1 . Для упрощения дальнейших преобразований положим, что реактивность ρ настолько мала, что ≈1. Тогда уравнение (12.43) запишется так:
(12.44)
Полученное уравнение является квадратным уравнением относительно и может быть записано следующим образом:
Отсюда
(12.45)
Корень квадратный вида при можно разложить в ряд Тейлора
Если принять l=10-3 с, ρ=0,003, то 4lλ р =4 10-3 0,077 0,003 = 0,925 10-6, откуда
Тогда, ограничившись в разложении двумя первыми членами, выражение (12.45) можно представить в следующем виде:
Отсюда
При определении пренебрегли вторым членом в квадратных скобках по сравнению с единицей. Если учесть теперь, что =0,077 10-3 , а
β—ρ=0,0064-0,0030=0,034 10-1 , то произведением по сравнению с (β-ρ) также можно пренебречь, и решение относительно со окончательно будет иметь вид
(12.46)
(12.47)
При положительном скачке реактивности и (β-ρ)>0 значение положительно, а отрицательно.
Выражение (12.38), описывающее изменение плотности нейтронного потока во времени, с учетом одной средневзвешенной группы запаздывающих нейтронов примет вид
(12.48)
а концентрация предшественников запаздывающих нейтронов в соответствии с (12.31) будет выражаться соотношением
(12.49)
Коэффициенты А и В в (12.48) и (12.49) являются постоянными интегрирования и определяются из начальных условий.
С учетом одной группы запаздывающих нейтронов уравнение (12.13).принимает вид
(12.50)
в котором связаны величины Ф и с.
Подставляя в (12.50) вместо Ф и с их выражения, согласно (12.48) и (12.49), получаем
Дифференцируя теперь выражение (12.49), будем иметь
Из сравнения последних двух уравнений видно, что правые части их будут равны, если приравнять соответствующие коэффициенты при одинаковых экспонентах:
отсюда коэффициенты B 0 и B1будут равны:
Подставляя в эти выражения вместо и их значения в соответствии с (12.46) и (12.47) и пренебрегая величиной по сравнению с величиной (β-ρ),
получаем:
(12.51)
(12.52)
Как видно из формул (12.51) и (12.52), связь между соответствующими постоянными А и В выражается посредством величин, определяющихся свойствами данного реактора, поэтому остается определить только две произвольные постоянные. Для этой цели требуются два начальных условия.
Условия в момент времени t =0 в соответствии с (12.48) и (12.49):
(12.53)
(12.54)
где Ф0 и c0 — соответственно плотность потока нейтронов и концентрация предшественников запаздывающих нейтронов до внесения возмущения.
Для нахождения связи между Ф0 и c 0 , учитывая, что при условии стационарной работы реактора концентрация предшественников не должна меняться со временем, потребуем, чтобы в момент времени t =0 = 0 или, согласно (12.50),
(12.55)
Из уравнения (12.55) следует, что
(12.56)
Подставив теперь в (12.54) вместо c0 его выражение из (12.56), а вместо B0 и B1 - их выражения в соответствии с (12.51) и (12.52), будем иметь
Используя (12.53), заменим в последнем выражении А1 на (Ф0 — A 0 ). Тогда
Второй член в скобках правой и левой частей, , имеющий порядок
10-3 , мал по сравнению с первыми, имеющими порядок 10-1 , и им можно пренебречь. С учетом этого:
(12.57)
. (12.58)
Подставив этот результат, а также (12.46) и (12.47) в выражение (12.48), получим формулу для плотности нейтронного потока в зависимости от времени
(12.59)
При положительном значении ρ и β-ρ>0 плотность нейтронного потока представляет собой разность двух членов, из которых положительный член растет со временем, а отрицательный быстро падает.
Для иллюстрации подставим в формулу (12.59) численные значения величин. Предположим, как и ранее, что l=10-3 с, ρ=0,003, β=0,0064, λ=0,077 с-1 , тогда
(12.60)
На рис. 12.3 сплошной линией показано изменение отношения Ф/Ф0 по формуле (12.60). На этом же рисунке пунктирными линиями показано изменение каждого члена правой части формулы (12.60) в отдельности.
Кривая 1 соответствует изменению первого члена, а кривая 2 — изменению второго. Из рисунка видно, что второй член очень быстро убывает и в данном случае уже к исходу первой секунды вклад его составляет около 0,2% и им можно пренебречь. Поэтому по истечении очень малого промежутка времени (порядка секунды) изменение плотности нейтронного потока практически определяется первым членом, и можно записать, что
Из последнего выражения следует, что установившийся период реактора равен
(12.62)
Если подставить численные значения β, λ и ρ, которые использовались ранее, то получим
Установившийся период реактора с учетом запаздывающих нейтронов составляет 14,7 с. Для этих же условий, но без учета запаздывающих нейтронов период равен 0,33 с. Таким образом, переходный процесс будет идти значительно медленней, чем без учета запаздывающих нейтронов.
При очень малых значениях реактивности (порядка сотых долей процента) в формуле (12.62) величиной р по сравнению с величиной ρ можно пренебречь и установившийся период реактора будет обратно пропорционален реактивности
Величина представляет собой средневзвешенное время жизни запаздывающих нейтронов τзн , т.е.
Это выражение аналогично формуле (12.12), только вместо времени жизни мгновенных нейтронов здесь фигурирует время жизни запаздывающих нейтронов. Поэтому установившийся период реактора полностью определяется свойствами запаздывающих нейтронов.
Период реактора является очень важной эксплуатационной характеристикой и всегда контролируется в процессе работы аппарата, в особенности во время пуска реактора. Величина , обратная периоду реактора, представляет собой скорость разгона аппарата. Чем больше период, тем меньше скорость разгона, и наоборот, с уменьшением периода реактора скорость разгона возрастает. При больших значениях реактивности период реактора мал и скорость разгона может оказаться настолько большой, что управление реактором будет невозможно. Если значение , то запаздывающие нейтроны практически не оказывают влияния на переходный процесс.
Вводится понятие так называемой мгновенной критичности [2], которая соответствует состоянию . Это значит, что реактор критичен уже по отношению только к мгновенным нейтронам. А так как наряду с генерацией мгновенных нейтронов идет генерация и запаздывающих нейтронов, то реактор в целом надкритичен с положительной реактивностью , и разгон его идет очень быстро с периодом, определяемым мгновенными нейтронами. Таким образом, значение ρ всегда должно оставаться меньше β.
При положительном скачкообразном изменении реактивности и плотность нейтронного потока в первый момент времени после внесения возмущения быстро возрастает, а затем скорость изменения плотности потока переходит постепенно к скорости, соответствующей установившемуся периоду реактора. Такой характер изменения плотности потока определяется тем, что мгновенные нейтроны реагируют на скачкообразное изменение реактивности за очень короткий промежуток времени. Однако при дальнейшее быстрое возрастание плотности потока прекращается, так как реактор по отношению к мгновенным нейтронам подкритичен и последующее более плавное изменение обусловлено действием запаздывающих нейтронов.
Легко показать, что изменение плотности потока в первый момент времени после внесения возмущения обусловлено действием мгновенных нейтронов. Для этого выражение (12.59) продифференцируем по времени и возьмем производную в момент t =0:
В полученном выражении левая часть представляет собой обратную величину периода. В этом легко убедиться, если продифференцировать обе части формулы (12.11). Тогда
При небольших значениях ρ первый член в правой части мал по сравнению
со вторым и им можно пренебречь. В этом случае
Данный результат совпадает с формулой (12.12), когда все нейтроны только мгновенные.
На рис. 12.4 показано изменение плотности потока нейтронов при скачкообразном изменении реактивности для различных значений р.
С увеличением реактивности возрастает первоначальный скачок плотности потока нейтронов и последующее установившееся изменение потока идет быстрее. Случай ρ=0 соответствует стационарному режиму, плотность потока нейтронов не изменяется. При приближении ρ к β плотность потока возрастает катастрофически быстро, и реактор становится практически неуправляем.
Дата: 2019-04-23, просмотров: 381.