Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Скачиваний:
61
Добавлен:
22.08.2013
Размер:
357.38 Кб
Скачать

Московский инженерно-физический институт (государственный университет) Физико-технический факультет

Лекция 10

Многогрупповой подход. Многогрупповое уравнение диффузии. Внутренние и внешние итерации. Программы нейтронно-физического расчета.

Коэффициенты чувствительности коэффициента размножения к изменению параметров реактора.

Ф6-01Н

Теория переноса излучений

Многогрупповой подход

Многогрупповой подход при решении уравнения диффузии с энергетической зависимостью основан на посылке, что в пределах каждого энергетического диапазона {Eg} функция

плотности потока нейтронов обладает свойством подобия по энергии.

Тогда исходное уравнение диффузии заменяется системой групповых уравнений.

Ф6-01Н

Теория переноса излучений

Многогрупповой подход

В групповых уравнениях диффузии величина gx

групповое сечение процесса типа x определяется:

dE x (E)S(E)

gx = Eg

dES(E)

Eg

где S(E) – спектр свертки (известная по форме функция внутри группы).

Ф6-01Н

Теория переноса излучений

Многогрупповое уравнение диффузии

Стационарное групповое уравнение диффузии нейтронов имеет вид:

Dg( r ) g ( r ) + g0(r ) Фg( r ) =

NG

g g '

NG

=g ' 1

s,0

(r ) g ' (r ) + g ( f f)g’( r) g '( r) + Qg( r ),

 

 

g ' 1

где

Ф g( r) и Qg(r ) - соответственно поток и

внешний источник нейтронов в группе g.

Ф6-01Н

Теория переноса излучений

Элементарные частицы

Групповое уравнение диффузии можно заменить системой конечно-разностных уравнений по пространственной переменной r. Запишем в определенной последовательности значения потока Ф по дискретным значениям в виде вектора .

Для описания компонент вектора используется единственный индекс: j = 1,2, … , - l) NG.

Тогда имеем операторный (матричный) вид уравнения диффузии:

L +

= S + F + Q

(1)

Обозначив A =

L

+

 

F

 

 

 

S – , рассмотрим уравнение (1) в виде:

A = Q

Ф6-01Н

Теория переноса излучений

Внутренние и внешние итерации

Запишем матрицу

A в виде суммы трех матриц:

 

A =

D U V

где

D – диагональная .матрица (отличные от нуля элементы

находятся только на основной диагонали),

U

– верхняя треугольная матрица (отличные от нуля элементы

находятся только выше основной диагонали),

V

– нижняя треугольная матрица (отличные от нуля элементы

находятся ниже основной диагонали).

Ф6-01Н

Теория переноса излучений

Внутренние и внешние итерации

Итерационный процесс можно определить следующим образом:

(i 1) = D 1 (U + V ) (i) + D 1 Q

Итерационный процесс продолжается до тех пор, пока разность между потоками (i) и (i 1) на двух последующих итерациях не будет меньше заданного критерия.

Ф6-01Н

Теория переноса излучений

Программы нейтронно-физического расчета

Программы нейтронно-физического расчета создаются для расчета некоторых характеристик рассматриваемой системы. Их можно классифицировать по ряду параметров.

Современные программы нейтронно-физического расчета представляют из себя комплексы программ, нацеленные на получение широкого круга характеристик рассматриваемой системы.

Ф6-01Н

Теория переноса излучений

Коэффициенты чувствительности коэффициента размножения к изменению параметров реактора

Основной характеристикой размножающих свойств среды является коэффициент размножения нейтронов Keff, который

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

Keff 1 Keff

Если реактор критичен, то реактивность равна нулю, если подкритичен, то реактивность меньше нуля, если надкритичен, то реактивность больше нуля.

Ф6-01Н

Теория переноса излучений

Коэффициенты чувствительности коэффициента размножения к изменению параметров реактора

Изменение реактивности, вызванное изменением физических параметров компонентов активной зоны, называют обратными связями.

Ф6-01Н

Теория переноса излучений

Соседние файлы в папке Prezentacii