Цилиндрический ЯР с боковым отражателем в одногрупповом приближении

T

H

R

В рамках одногруппового приближения аналитически может быть решена только задача, в которой отражатель расположен только на основаниях ци­линдра. В том случае, когда активная зона окружена отражателем со всех сторон, переменные r и z не разделяются и аналитическое решение в рамках одногруппового приближения не может быть получено. Такой случай будет рассмот­рен ниже при анализе двухгруппового приближения.

Рассмотрим ЯР, представляющий собой АЗ высотой Н (экстраполированный размер) и радиусом R, окруженную боковым отражателем, толщиной Т. Начало координат находится в центре симметрии. Так как на торцах ЯР отсутствует отражатель, то его влияние скажется на радиальную составляющую потоков, т.е. аксиальная составляющая потока будет точно такой, как в случае ЯР без отражателя. Предположим, что переменные r и z разделяются и будем искать решение в виде:

Условия симметрии удовлетворяются, а граничные условия на торцевых плоскостях цилиндра:

Таким образом, решение задачи сводится к нахождению функций f ₁(r) и f ₂(r).

Для функции f (r) запишем уравнения ЯР:

(1)

(2)

где – радиальная составляющая материального параметра АЗ; – радиальная составляющая материального параметра отражателя.

Определим эти составляющие.

Активная зона:,

где – полный материальный параметр; – аксиальная составляющая материального параметра. При этом вследствие того, что нет бокового отражателя, аксиальная составляющая материального параметра равна аксиальной составляющей геометрического параметра такого же реактора без отражателя (реактор критический):. Тогда получаем:

. Отсюда видно, что - величина действительная.

По аналогии рассмотрим радиальную составляющую материального параметра отражателя, учтя при этом, что в отражателе нет делящихся материалов (k _∞=0)

. Отсюда видно, что - величина мнимая. Таким образом, исходные уравнения примут вид:

(1)

(3)

Эту систему уравнений необходимо дополнить граничными условиями:

f ₂(R ₂) = 0, где R ₂= R + T (4)

f ₁(R) = f ₂(R) (5) (равенство потоков на границе АЗ-отражатель)

(6) (равенство диффузионных токов на АЗ –

отражатель)

Решение уравнения (1) известно

(7)

Решим уравнение (3). В цилиндрических координатах оно имеет вид:

(8)

Уравнение (8) помножим на r ², затем первое слагаемое умножим и разделим на, второе слагаемое -.

(9)

Уравнение (9) – уравнение Бесселя с аргументом. Причем этот аргумент – мнимый. В этом случае решением уравнения Бесселя являются модифицированные функции Бесселя нулевого порядка:

, (10)

где – модифицированная функция Бесселя нулевого порядка от мнимого аргумента первого рода; – модифицированная функция Бесселя нулевого порядка от мнимого аргумента второго рода (функция Макдональда).

Для решения (10) воспользуемся граничным условием (4):

, отсюда:

(11)

Подставим (11) в (10) и для потоков в АЗ и отражателе окончательно получим:

(12)

(13)

Для (12) и (13) используем граничные условия (5) и (6):

(14)

(15)

Разделим (15) на (14) и получим:

, (16)

где R ₂= R + T.

По аналогии с ЯР в форме пластины докажем, что условие (16) является условием критичности, т.е. если уберем отражатель, условие (16) должно прейти в условие критичности реактора без отражателя:.

Пусть нет отражателя (Т =0). Тогда знаменатель в правой части (16) обращается в 0, а сама правая часть стремится к бесконечности. Тогда

Т.к. функция ограничена при всех х, то указанное условие будет выполняться только при или (первый корень функции Бесселя первого рода нулевого порядка), что и требовалось доказать.

Как и в случае плоского ЯР введем понятие эффективной добавки за счет отражателя δ = R ₀ – R, где R ₀ – критический радиус эквивалентного цилиндрического ЯР без отражателя. Выразим R и подставим в левую часть условия критичности (16):

Рассмотрим ЯР, имеющий большие размеры, т.е. радиус велик по сравнению с величиной эффективной добавки δ << R. Тогда разложим (17) в ряд Тейлора по малому параметру, ограничившись первым числом разложения (достаточно сложные математические выкладки опустим).

(17)

С другой стороны в большом реакторе R >> M ₂, следовательно,. В этом случае для правой части условия (16) воспользуемся асимптотическим разложением модифицированных функций Бесселя:

;.

Тогда правая часть (16) примет вид:

(18)

В итоге сопоставляя (17) и (18), для больших ЯР получим условие критичности:

(19)

Это выражение с точностью до первых членов разложений совпадает с результатом, полученным для плоского ЯР. Надо заметить, что с уменьшением R необходимо учитывать больше членов разложения вследствие увеличения кривизны поверхности и полученное для δ выражение усложнится. Далее необходимо провести анализ δ при различных толщинах отражателя.

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

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