Научная литература
booksshare.net -> Добавить материал -> Физика -> Белл Д. -> "Теория ядерных реакторов" -> 52

Теория ядерных реакторов - Белл Д.

Белл Д. Теория ядерных реакторов — Москва, 1974. — 494 c.
Скачать (прямая ссылка): teoriyayadernihreaktivov1974.djvu
Предыдущая << 1 .. 46 47 48 49 50 51 < 52 > 53 54 55 56 57 58 .. 264 >> Следующая


(1) a„,n> 0; an,m<0, если тфп- я„. m = am.

(2) аПш,1+1Ф 0;

(3) |an, n|>2 la,,.», I-

тфп

Из свойства (2) следует, что матрица является неприводимой; физически э го означает, что нейтрон из произвольной данной точки может попасть в любую другую точку пространственной сетки. Свойство (3) устанавливает пренмуще-

109
ство диагональных членов, т. е. то, что сечение о0, определенное в соответствии с уравнением (3.6) и в предположении, отмеченном после уравнения (3.1), является положительным. Тогда говорят, что А является неприводимой диагонально преобладающей матрицей. Это свойство гарантирует, что матрица невырожденна и что существует обратная матрица; в этом случае и решение

для ф, определяемое уравнением (3.26), существует [71.

До сих пор ничего не говорилось об условиях на поверхностях между областями, указаннымн на рис. 3.2. В действительности эти условия автоматически удовлетворяются разностными уравнениями. Рассмотрим очень мелкую пространственную сетку, такую, что все значения А стремятся к нулю вблизи поверхности при хи. Тогда из уравнений (3.16) — (3.18) следует, что и ф, и J являются непрерывными функциями при пересечении поверхности.

При выводе конечно-разностных уравнений для Р^приближения использовались конкретные аппроксимации интегралов, входящих в уравнения (3.14) и (3.15). Другие простые аппроксимации привели бы к разностным уравнениям, аналогичным уравнению (3.20), за исключением того, что коэффициенты аПуТ11 в них были бы несколько отличными от рассмотренных, однако эти коэффициенты вновь имели бы отмеченные выше свойства. Диффузионное приближение приводит к таким же разностным уравнениям, хотя и с другими коэффициентами. Так как диффузионное приближение используется очень широко, то интересно получить соответствующие конечно-разностные уравнения.

3.2.4. КОНЕЧНО-РАЗНОСТНЫЕ УРАВНЕНИЯ

В ДИФФУЗИОННОМ ПРИБЛИЖЕНИИ

Диффузионное приближение можно рассматривать эквивалентным первому уравнению Рх-приближения, т. е. уравнению (3.7):

+O0(X) ф W = Q0W (3.30)

ах

змее те с законом Фи ка

J(x)= -D(X) (3.31)

dx

Конечно-разностное уравнение, соответствующее уравнению (3.30), имеет такой же вид, как и для уравнения (3.7), а именно (3.16). Для токов в точках

хк+і/2 и *fc-i/2 уравнение (3.31) можно аппроксимировать выражениями

/*+1/2 &-Dk+и2 фь+1~фк ; (3.32)

Aft+ 1/2

Jk-i,2~-Dk-l/2 . (3.33)

1/2

Подставляя соотношения (3.32) и (3.33) в уравнение (3.16), получаем ко-нечно-разностное уравнение, имеющее вид (3.20), с коэффициентами

п °Ь+1/2

ак,к+1— *------— ак+1, кг

afc+l/2

ак, к = ^ok ак, h-І — ак, Ь+Г»

Dk_ ,/о

к-1 — 7 — 1, к»

А/г- 1/2

Sk — Afc Qok ¦

Вновь можно показать, что элементы матрицы А имеют отмеченные выше свой ст ва.

HO
При выводе конечно-разностных уравненийрассматривалась пространственная область, расположенная между точками хк^/г и хк+1/г. Различные члены в уравнении имеют вполне определенный смысл при нахождении баланса нейтронов в этой области. Более подробно это показано ниже при рассмотрении конечно-разностных уравнений в сферической геометрии. Таким образом, конечно-разностное уравнение можно рассматривать как уравнение баланса нейтронов для небольшой области в системе. В конечно-разностных уравнениях особенно важно обеспечить это свойство сохранения числа нейтронов, чтобы можно было прослеживать судьбу всех нейтронов деления при численных расчетах. В расчетах критичности баланс между производством и потерей нейтронов носит, конечно, решающий характер, следовательно, существенно, чтобы нейтроны искусственно не возникали или не исчезали.

3.2.5. РЕШЕНИЕ УРАВНЕНИЙ Р^-ПРИБЛИЖЕНИЯ

В предыдущих разделах описывался численный метод решения уравнений P1-и диффузионного приближений в плоской геометрии. Уравнения Рдг-приближения для более высоких значений N можно привести к системе разностных уравнений подобным же образом. Существует несколько методов решения этих уравнений [8], и один из наиболее гибких методов описан в гл. 5. Кроме того, в разд. 3.5.2 показано, что в плоской геометрии «двойное P.v-приближение» оказывается лучшим, чем простое Рд>-приблнжение. Некоторые результаты, полученные этими двумя методами, представлены в гл. 5.

3.3. РАЗЛОЖЕНИЕ ПОТОКА В СФЕРИЧЕСКОЙ И ПРОИЗВОЛЬНОЙ геометриях

3.3.1. РАЗЛОЖЕНИЕ В СФЕРИЧЕСКОЙ ГЕОМЕТРИИ

До сих пор обсуждение метода сферических гармоник касалось плоской геометрии. Здесь же рассмотрено применение этого метода и к другим геометриям. Для системы, симметричной относительно некоторой точки, можно использовать сферические координаты. Ниже показано, что уравнения метода сферических гармоник в таких координатах очень похожи на те же уравнения в плоской геометрии. Такие системы рассмотрены в настоящем разделе, а более общие геометрии, для которых разложение потока нейтронов в ряды по полиномам Лежандра неприменимо, описаны в разд. 3.3.3 для Р^приближения. Использование метода сферических гармоник в цилиндрической геометрии рассмотрено в разд. 3.6.2.
Предыдущая << 1 .. 46 47 48 49 50 51 < 52 > 53 54 55 56 57 58 .. 264 >> Следующая
Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

Есть, чем поделиться? Отправьте
материал
нам
Авторские права © 2009 BooksShare.
Все права защищены.
Rambler's Top100

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed