Порядок выполнения работы и обработка результатов
ОГЛАВЛЕНИЕ
ВВЕДЕНИЕ ................................................................................................. 4
ЭКСПЕРИМЕНТАЛЬНОЕ ОПРЕДЕЛЕНИЕ ДЛИНЫ ЭКСТРАПОЛЯЦИИ В ПОЛИЭТИЛЕНЕ ........................................................................... 5
ВЛИЯНИЕ ОТРАЖАТЕЛЯ НА РАСПРЕДЕЛЕНИЯ ПОТОКОВ НЕЙТРОНОВ ТЕПЛОВОЙ И НАДТЕПЛОВОЙ ЭНЕРГИЙ .................................................................................................................... 16
МЕТОД Экспоненциальной ПРИЗМЫ .......................................... 23
ИЗМЕРЕНИЕ ДОЛИ ПОГЛОЩЕНИЙ ТЕПЛОВЫХ НЕЙТРОНОВ МЕТОДОМ ЭКСПОНЕНЦИАЛЬНОЙ ПРИЗМЫ.................................................................................................................... 34
МОДЕЛИРОВАНИЕ РЕАКТОРА НА ЭЛЕКТРИЧЕСКОЙ СЕТОЧНОЙ МОДЕЛИ.................................................................................................................... 41
МОДЕЛИРОВАНИЕ ЭЛЕМЕНТАРНОЙ ЯЧЕЙКИ РЕАКТОРА............ 50
МОДЕЛИРОВАНИЕ НЕЙТРОННОГО ПОТОКА В ГРАФИТОВОЙ ПРИЗМЕ.................................................................................................................... 57
ВВДЕНИЕ
Курс «Физика и расчет ядерного реактора» является одной из профилирующих дисциплин специальности 070500 «Ядерные реакторы и энергетические установки». Для удовлетворения растущих энергетических потребностей требуется научно, экономически, экологически и политически обоснованный подход к использованию ядерной энергии, учитывающий как связанные с ней опасности, так и большие потенциальные возможности. В курсе рассматриваются основы физики ядерных реакторов, излагаются современные методы расчета ядерно-энергетических установок.
Лабораторный практикум является неотъемлемой составной частью дисциплины иимеет своей целью приобретение навыков и умений в обращении с радиоактивными источниками, тестовой и контролирующей аппаратурой, правильной и корректной обработке полученных результатов, понимании физического смысла результатов математического расчета; все исследования, проводимые с радиоактивными источниками, производятся под непосредственным контролем преподавателя и соблюдением правил техники безопасности; оформленный отчет, содержащий выводы о проделанной работе, является критерием выполнения работы и приобретенных навыков. Он включает семь лабораторных работ:
· экспериментальное определение длины экстраполяции в полиэтилене;
· влияние отражателя на распределения потоков нейтронов тепловой и надтепловой энергии;
· метод экспоненциальной призмы;
· измерение доли поглощений тепловых нейтронов методом экспоненциальной призмы;
· моделирование реактора на электрической сеточной модели;
· моделирование элементарной ячейки реактора;
· моделирование нейтронного потока в графитовой призме;
и своей тематикой охватывает основные разделы курса «Физика и расчет ядерного реактора».
ЭКСПЕРИМЕНТАЛЬНОЕ ОПРЕДЕЛЕНИЕ ДЛИНЫ ЭКСТРАПОЛЯЦИИ В ПОЛИЭТИЛЕНЕ
Цель работы: экспериментальное определение экстраполированных размеров среды в воздухе (вакууме) на плоской границе раздела «среда – воздух»
Теоретические основы
Одной из важнейших задач физики ядерного реактора на тепловых нейтронах является определение пространственного распределения тепловых нейтронов в активной зоне. Эту задачу можно решить, используя уравнение диффузии, которое для стационарного случая имеет следующий вид:
, (1)
где Ф(r) – плотность потока нейтронов; D – коэффициент диффузии; – макроскопическое сечение поглощения нейтронов в среде; S(r)-функция, описывающая распределение потока нейтронов источника.
Уравнение диффузии достаточно точно описывает поведение нейтронов в слабопоглощающих средах , вдали от источников нейтронов, границ раздела сред с различными свойствами и сильных поглотителей (на расстояниях, равных примерно (2¸3)
).
Так как уравнение диффузии – дифференциальное уравнение второго порядка и, соответственно, его общее решение содержит две произвольные постоянные, то для единственного решения любой конкретной задачи необходимо иметь два граничных условия. Кроме того, физически очевидно, что нейтронный поток должен быть конечным и неотрицательным, не исключая его обращения в нуль. Согласно диффузионного приближения граничные условия формулируются следующим образом: на границе раздела двух сред потоки и плотности результирующих токов нейтронов равны.
Особый случай возникает, если среда, где диффундируют нейтроны, граничит с вакуумом, так как вакуум обладает тем свойством, что попавший в него нейтрон никогда не возвращается обратно в среду. В этом случае для решения уравнения диффузии используют формальный математический прием: вводят в рассмотрение экстраполированную границу среды Rэ = R + d (где R – реальный размер среды; d – длина экстраполяции) и полагают, что поток нейтронов на экстраполированной границе обращается в нуль Ф(Rэ) = 0.
![]() |
Как уже отмечалось, диффузионное приближение не дает достаточно точных результатов вблизи границы раздела сред. В этой области имеет место искажение нейтронного потока. В этом случае истинное распределение потока нейтронов можно представить в виде суперпозиции двух слагаемых, одно из которых описывает распределение потока нейтронов в рамках диффузионного приближение во внутренних областях среды, другое – распределение потока в приграничных областях. Так, на рис.1 приведено распределение потока нейтронов вблизи границы слабопоглощающей среды с вакуумом в одномерном случае.
Тогда полный поток нейтронов можно представить следующим образом:
Ф(x) = Фас(x) + Фпер(x), (2)
где Фас – асимптотическая часть полного потока, описывающая нейтронный поток в области, где справедливо диффузионное приближение; Фпер – переходная часть полного потока, учитывающая отклонение от диффузионного приближение в приграничных областях.
При таком рассмотрении под длиной экстраполяции понимается расстояние, отсчитываемое от реального размера среды, на котором асимптотическая часть полного потока обращается в нуль при условии, что в вакууме Фас является линейной функцией от расстояния (рис.1).
Задача определения величины длины экстраполяции может быть решена теоретически. В соответствии с теорией диффузии для слабопоглощающей среды с плоской границей длина экстраполяции связаны с транспортной длиной следующим образом:
. (3)
Кроме того, данную задачу можно решить экспериментальным путем. Для этого необходимо измерить величины нейтронных потоков в различных областях среды, аппроксимировать полученные значения функцией, получаемой при решении уравнения диффузии для данного типа среды и источника, по результатам аппроксимации определить асимптотический участок в распределении нейтронного потока и провести его линейную экстраполяцию.
Если имеется слабопоглощающая среда (например, полиэтилен), в которой размещен лабораторный нейтронный источник (например, плутоний-бериллиевый), то на некотором расстоянии он может рассматриваться как точечный. В этом случае решение уравнения диффузии дает следующий результат:
, (4)
где – скорость генерации тепловых нейтронов источником, нейтр./с (для конкретного источника величина постоянная); D – коэффициент диффузии, см;
– параметр, определяющий свойства среды, см-1:
. (5)
Одним из наиболее простых способов определения зависимости Ф(r) является способ, основанный на использовании активационных детекторов, которые размещаются на различных расстояниях от источника и облучаются в потоке нейтронов. При этом детекторы, размещенные в разных точках приобретут разную активность. Причем активность детекторов A(r) прямо пропорциональна величине потока на данном расстоянии: A(r) = kФ(r) (k – коэффициент пропорциональности). Тогда выражение (4) примет следующий вид:
, (6)
где B – постоянная:
. (7)
Таким образом, измерив активность детекторов на различных расстояниях от источника, аппроксимировав полученные экспериментальные точки зависимостью (6) и проанализировав полученный результат, появляется возможность определения асимптотического участка в распределении A(r) (а значит, и потока нейтронов), его линейной экстраполяции и нахождения длины экстраполяции.
В этом случае основной проблемой становится определение постоянных в аппроксимационной зависимости (6) B и . Одним из наиболее эффективных методов аппроксимации экспериментальных данных зависимостями того или иного вида является метод наименьших квадратов, сущность которого состоит в том, что неизвестные параметры аппроксимирующей функции должны быть выбраны так, чтобы сумма поделенных на дисперсии квадратов отклонений от функции экспериментальных точек была минимальной:
, (8)
где и
– экспериментальные значения и теоретические значения аппроксимирующей функции в k-ой точке;
– дисперсия экспериментальной k-ой точки; N – число экспериментальных точек.
Условием минимума функции M является равенство нулю ее производных по искомым параметрам. Таким образом, применительно к исследуемому случаю исходным условием для нахождения неизвестных постоянных функции (6) является следующая система уравнений:
. (9)
Прологарифмировав функцию (6), условие (8) примет вид:
, (10)
где Aэ(rk) – экспериментально полученное среднее значение активности детектора, расположенного на расстоянии rk от источника; b – постоянная, b = lnB; – дисперсия определения величины lnAэ(rk),
.
Дифференцируя функцию (10) по искомым параметрам b и и проведя ряд преобразований, исходная система уравнений примет вид:
, (11)
где коэффициенты ;
;
;
;
;
.
Данная система уравнений может быть решена методом Крамера:
,
где D – определитель матрицы коэффициентов при неизвестных величинах:
;
и
– определители, полученные из определителя D путем замены на столбец свободных членов первого и второго столбца рассматриваемой матрицы, соответственно:
Таким образом, искомые постоянные b и определяются следующим образом:
. (12)
Полученные значения b и представляют собой наиболее вероятные значения искомых величин. Подстановка их в исходное уравнение (10) даст выполнение условия минимума функции M. Очевидно, что определение неизвестных характеризуется своей погрешностью. Согласно теории [4] среднеквадратичная погрешность определения неизвестной величины x рассчитывается следующим образом:
, (13)
где – наиболее вероятные значения среднеквадратичной погрешности величины x; px – статистический вес определения величины x.
Статистический вес неизвестного равен определителю матрицы коэффициентов при неизвестных величинах, деленному на алгебраическое дополнение того элемента главной диагонали этого определителя, который является коэффициентом при рассматриваемом неизвестном. В данном случае статистический вес величин b и равен:
. (14)
Наиболее вероятное значения среднеквадратичной погрешности неизвестной величины определяется следующим образом:
, (15)
где M0 – минимальное значение функции M, получаемое при подстановке в соотношение (10) найденных средних значений b и ; N – число экспериментальных точек; m – число неизвестных (в рассматриваемом случае m=2). Тогда дисперсии определения неизвестных постоянных b и
определяются соотношениями вида:
. (16)
Зная, что b = lnB, искомые коэффициенты аппроксимации B и в конечном виде могут быть представлены следующим образом:
, (17)
где величины b и определяются по формулам (12), а значения
и
– по соотношениям (16).
Аппроксимируя экспериментальные точки зависимостью (6) с найденными значениями постоянных B и , можно определить асимптотический участок в распределении нейтронного потока, а затем, проведя его линейную экстраполяцию, рассчитать длину линейной экстраполяции d и найти экстраполированные размеры среды в воздухе.
Материалы и оборудование
В качестве слабопоглощающей среды используется полиэтиленовая призма, в которой обеспечивается расстояние от источника нейтронов, помещенного внутрь призмы, до поверхности, равное 24 см. Источником нейтронов служит лабораторный плутоний-бериллиевый источник мощностью 106¸107 нейтр./с. На различных расстояниях от источника размещаются активационные детекторы, представляющие собой индиевые фольги (период полураспада индия 54,5 мин.). Причем на каждом расстоянии располагается по два детектора: один – заключен в кадмиевый экран, другой – находится без экрана. Активность облученных детекторов измеряется торцевым бета-счетчиком, помещенным в защитный домик. Сигналы с детектора регистрируются пересчетным прибором типа ПСО2-4.
Порядок выполнения работы и обработка результатов
1. Предварительно преподавателем производится установка источника нейтронов и размещение активационных детекторов на различных расстояниях от источника в полиэтиленовой призме.
2. Измеряется фоновая скорость счета. Измерения проводятся 5¸7 раз, затем исключаются наибольшее и наименьшее значение измеренного фона, а по оставшимся значениям определяется среднее значение фоновой скорости счета и погрешность его определения.
3. Производится измерение активности детекторов, размещенных на различных расстояниях от источника. Активность фольг определяется по скорости регистрации излучения, испускаемого детекторами. Для каждой фольги измерения проводятся аналогично п. 2. При этом необходимо учитывать тот факт, что измеренная активность фольг меньше активности насыщения из-за переносов детекторов до места регистрации, нахождения детекторов в защитном домике. Если время облучения достаточно велико, то скорость счета, соответствующая активности насыщения детектора, может быть определена по следующей формуле:
, (18)
где N – зарегистрированная скорость счета; tв – время выдержки (время между извлечением фольги из призмы и началом счета); tc – время счета. Необходимо помнить, что, так как измерения одной фольги проводятся несколько раз, время выдержки для второго измерения складывается из времени между извлечением фольги из призмы и началом счета (времени выдержки первого измерения) и времени счета первого измерения, для третьего – из времени выдержки второго и времени счета второго измерения и т.д.
После определения среднего значения скорости счета, соответствующей активности насыщения детектора, и его погрешности из полученного результата исключается фон с учетом погрешности его определения. По полученным значениям для каждого расстояния рассчитывается активность, наведенная тепловыми нейтронами.
Результаты измерений и расчетов заносятся в таблицу.
r, см | Номер измерения | ![]() | ![]() | ![]() | ![]() | ![]() |
2 | ||||||
. . . | ||||||
4 | ||||||
. . . | ||||||
. . . |
r – расстояние от источника нейтронов до места расположения детектора, см; N и Ncd – зарегистрированные скорости счета активности фольг без кадмиевого экрана и с ним, соответственно, с-1; и
– среднее значение величины N и погрешность ее определения, соответственно;
и
– среднее значение величины Ncd и погрешность ее определения, соответственно;
и
– среднее значение скорости счета, обусловленной тепловыми нейтронами, и погрешность ее определения, соответственно:
;
.
4. На графике Nт = f(r) наносятся экспериментально полученные значения Nт со своими доверительными интервалами. Производится аппроксимация экспериментальных точек зависимостью (6). С этой целью рассчитываются постоянные аппроксимации и погрешности их определения по соотношениям (12) и (16). Если аппроксимирующая кривая проходит через доверительные интервалы всех экспериментальных точек, можно утверждать, что на всем расстоянии "источник-поверхность" справедливо диффузионное приближение. В противном случае необходимо последовательно исключать из рассмотрения точки, расположенные вблизи границы и источника, каждый раз проводя аппроксимацию оставшихся экспериментальных точек зависимостью (6). Указанную процедуру следует проводить до тех пор, пока аппроксимирующая кривая удовлетворительно не опишет оставшиеся точки. Такой участок является асимптотическим в полном распределении нейтронного потока.
5. Производится линейная экстраполяция асимптотического участка путем проведения касательной к полученной кривой в последней точке асимптотического участка до пересечения ее с осью r. Рассчитывается длина линейной экстраполяции и погрешность ее определения.
6. Составляется отчет о работе, который должен включать следующее:
– кратко необходимые теоретические сведения и соотношения;
– математическую процедуру получения постоянных аппроксимации экспериментальных данных по методу наименьших квадратов;
– теоретический расчет величин ,
, d по табулированным значениям ядерно-физических констант для полиэтилена;
– результаты измерений активности фольг, размещенных на различных расстояниях от источника;
– результаты расчетов активности, обусловленной тепловыми нейтронами;
– процедуру определения асимптотического участка в распределении потока тепловых нейтронов, расчет постоянных аппроксимации B и и погрешностей их определения;
– результаты линейной экстраполяции асимптотического участка, расчет длины линейной экстраполяции и погрешности ее определения;
– сравнение полученных экспериментальных данных по величинам , d с теоретически рассчитанными;
– аргументированные выводы по работе.
Литература
1. Фейнберг С.М., Шихов С.Б., Троянский В.В. Теория ядерных реакторов: Учебник для вузов. – М.: Атомиздат, 1978.
2. Бартоломей Г.Г., Бать Г.А. и др. Основы теории и методы расчета ядерных энергетических реакторов: Учебное пособие для вузов. – М.: Энергоатомиздат, 1989.
3. Белл Д., Глесстон С. Теория ядерных реакторов. – М.: Атомиздат, 1974.
4. Щиголев Б.М. Математическая обработка наблюдений. – М.: Гос. изд-во физико-математической литературы, 1962.
Контрольные вопросы и задания
1. Общая характеристика диффузии молекул и диффузии нейтронов.
2. Основные понятия и допущения диффузионного приближения.
3. Понятие диффузионной плотности тока.
4. Что такое асимптотическое распределение плотности потока нейтронов?
5. Что такое переходное распределение плотности потока нейтронов?
6. Записать стационарное и нестационарное уравнение диффузии.
7. Дать физическую интерпретацию каждому входящему в уравнение диффузии слагаемому.
8. Граничные условия уравнения диффузии.
9. Десять последовательно установленных индиевых фольг одинакового размера и толщиной 0,2 г/см2 каждая облучаются в параллельном потоке тепловых нейтронов, падающем по нормали к поверхности фольг. Определить, во сколько раз активность десятой фольги меньше активности первой.
10. Что такое длина экстраполяции?
11. Что такое и с какой целью вводится понятие "экстраполированный размер активной зоны"?
12. Как он связан с параметрами среды?
13. Рассчитать коэффициенты диффузии, длины диффузии и длины линейной экстраполяции в воздухе для полиэтилена и воды.
14. Транспортные параметры.
15. Длина линейной экстраполяции в диффузионном и газокинетическом приближении.
16. Зависимость длины линейной экстраполяции от геометрических и материальных параметров среды.
17. Параметры нейтронных полей.
18. Получить решение стационарного уравнения диффузии для точечного источника нейтронов в однородной бесконечной среде.
19. Использование длины экстраполяции в физике ядерных реакторов.
20. Экспериментальные методы определения экстраполированных размеров среды.