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

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

Белл Д. Теория ядерных реакторов — Москва, 1974. — 494 c.
Скачать (прямая ссылка): teoriyayadernihreaktivov1974.djvu
Предыдущая << 1 .. 69 70 71 72 73 74 < 75 > 76 77 78 79 80 81 .. 264 >> Следующая


(Bg' -*¦ g )т,т — °s0 , g' -* g (хт)’

а т-я компонента вектора vFg есть AmvFg (хт). В общем виде vFg можно записать следующим образом:

vFg = 2 Cg'_>g ф (4.59)

g'

где Cg'-g — диагональная матрица с неотрицательными компонентами, которые для плоской геометрии определяются выражением

(?g' -*g)m.m = v07, g' C*m)-

Если уравнение (4.58) подставить в (4.57), то получим, что

АеФі= 2 вг'-*g Фі' + vFg/&. (4.60)

g'?=g

Система конечно-разностных уравнений (4.59) и (4.60) есть многогрупповая задача на собственное значение, которую следует решать методом внешних итераций, описанным в разд. 4.4.4. Решение дает эффективный коэффициент размножения вместе с соответствующей собственной функцией ф g для каждой группы, т. е. g — 1,2, ...,G. В рассматриваемом случае схема внешних итераций, которая была представлена в Р^приближении уравнениями (4.45) —

(4.48), выражается уравнением

А ,?*“’= (4.61)

gVg

соответствующим уравнениям (4.45) и (4.46), где

vFgl-!) — У ф{3п~1)- (4.62)

г'

Кроме того, из уравнения (4.47) можно рассчитать приближенное значение k, т. е. k{n), в виде

kin) =-----f (4.63)

Vf(«- *»/*<«—¦)

где в принятой системе обозначений vF{n) — полный источник нейтронов деления — можно представить как

Vf^=SvFy0-I. (4.64)

g

Здесь I — вектор, имеющий единичные компоненты, такой, что vFg-І представляет собой сумму по элементам объема или объемный интеграл от vFg.

Как обсуждалось в разд. 4.4.4, метод решения уравнений (4.61)—(4.64)

состоит в выборе пробной функции x'Fg0>/k<°\ расчете <pg1} из уравнения (4.61) и последующем получении величины VFglyIkP) из уравнений (4.62), (4.63) и (4.64). Итерационный процесс продолжается до тех пор, пока не будет достигнута его сходимость.

.Бьыа исследована система уравнений (4.59) и (4.60) и показано [20], что наибольшее собственное значение k является положительным и простым, а также, что соответствующий ему единственный собственный вектор может быть выбран таким образом, чтобы иметь неотрицательные компоненты. Кроме того, было доказано, что метод итераций по источникам деления сходится к этому собственному вектору. Эти выводы аналогичны описанным в разд. 4.4.3 для многогрупповых уравнений с непрерывной пространственной зависимостью потока нейтронов. К тому же они обеспечивают прочную основу для использования метода внешних итераций. Как и в случае внутренних итераций, имеются различные методы для ускорения сходимости внешних итераций [21].

153
Были рассмотрены также дискретные нестационарные многогрупповые уравнения, полученные добавлением к левой части уравнения (4.54) члена (alv„)d<p Jdt при k = 1 [22]. Решение этой краевой задачи имеет экспоненциальную временную зависимость, пропорциональную ехр (а^) при t-*¦ оо. Следовательно, критическое состояние системы можно определить, основываясь на знаке а. Результаты, приведенные в разд. 1.5 для общей теории переноса нейтронов и разд. 4.4.3 для многогруппового диффузионного приближении с непрерывной пространственной зависимостью потока нейтронов, распространяются и на многогрупповое диффузионное приближение с дискретным пространственным представлением потока нейтронов. Кроме того, коэффициент перед экспоненциальным решением дается в виде произведения вектора начального потока нейтронов и нормированного положительного собственного вектора сопряженных уравнений (см. гл. 6). Когда в уравнении присутствует источник, то ограниченное нестационарное решение при t—*-oо можно получить только для подкритической системы, что находится в соответствии с физическими соображениями, изложенными в разд. 1.5.4.

4.4.7. ВНЕШНИЕ ИТЕРАЦИИ В МНОГОГРУППОВОМ

P1-ПРИБЛИЖЕНИИ

Проведенный в предыдущем разделе анализ" был основан на многогрупповом диффузионном приближении. Для большинства других приближений, включая и Рх-приближение, также проведен соответствующий математический анализ. Часто нельзя сделать никаких выводов, поскольку система конечно-разностных уравнений не будет соответствовать положительному оператору [23]. Тем не менее общая стратегия внешних итераций успешно применялась в большинстве многогрупповых задач, включая, например, основанные на разложении в ряд по сферическим гармоникам или на методе дискретных ординат (см. разд. 5.4.3), в которых метод внешних итераций не имеет прочной математической основы. При таких условиях метод внешних итераций не всегда должен приводить к устойчивому численному решению; тем не менее на практике он оказывается очень плодотворным.

В связи с соотношением между P1- и диффузионным приближениями интересно выяснить, почему результаты анализа диффузионного приближения, проведенного в разд. 4.4.6, неприменимы к Р^приближению. Объяснение этому можно найти в выражении (3.24), которое описывает компоненту вектора источника для односкоростной задачи в Р^приближении. Соответствующие многогрупповые значения Q0 и Q1 определяются уравнениями (4.52) и (4.53) соответственно. Комбинируя эти результаты, находим, что вектор источника в Р^приближении невозможно выразить в простом виде уравнения (4.58), применимого в диффузионном приближении.
Предыдущая << 1 .. 69 70 71 72 73 74 < 75 > 76 77 78 79 80 81 .. 264 >> Следующая
Реклама

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed

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

c1c0fc952cf0704ad12d6af2ad3bf47e03017fed