Главная > Математика > Численные методы Монте-Карло
<< Предыдущий параграф
Следующий параграф >>
<< Предыдущий параграф Следующий параграф >>
Макеты страниц

§ 3. Пример: рассеяние частиц

3.1. Основное уравнение теории рассеяния.

Методы Монте-Карло часто используются для расчета различных задач, связанных с прохождением частиц (нейтронов, гамма-квантов и др.) через вещество. Мы не будем касаться здесь специальных вопросов, а рассмотрим лишь общую схему рассеяния, предполагая, что частицы при столкновении с атомами (точнее, с ядрами атомов) среды могут либо рассеиваться, либо поглощаться. Обозначим через Р точку шестимерного фазового пространства координат и скоростей v частицы. Обозначим элемент объема этого пространства через

Пусть количество первых столкновений, количество всех столкновений в элементе объема около точки Р (за единицу времени). Функцию можно явно вычислить, если задан источник частиц. Функцию называемую плотностью столкновений, требуется найти.

Введем ядро столкновений , которое определяется следующим условием: это вероятность того, что частица, испытавшая столкновение в точке испытает следующее столкновение в элементе объема около точки Р (за единицу времени). Конкретный вид ядра столкновений в одногрупповой теории переноса нейтронов имеется на стр. 223.

Нетрудно составить интегральное уравнение, которому подчиняется плотность столкновений:

ибо столкновение в окрестности точки Р может быть либо первым столкновением, либо следует за столкновением в окрестности некоторой точки , а количество таких столкновений (за единицу времени) равно Область интегрирования в (49) — все пространство. Введем сопряженное ядро . Тогда уравнение (49) совпадает с уравнением (25), а сопряженное уравнение (46) можно записать в виде

В рассматриваемом случае итерации - функции имеют простой физический смысл: это плотности вторых, третьих и т. д. столкновений. И ряд Неймана

означает, что плотность столкновении есть сумма плотностей первых, вторых и т. д. столкновений.

Обычно требуется вычислить какие-нибудь функционалы вида . Это можно сделать любым из методов, указанных в § 2.

<< Предыдущий параграф Следующий параграф >>
Оглавление