- •Введение в диффузию моноэнергетических нейтронов
- •Обзор предыдущего материала
- •Сведение уравнения Больцмана к диффузионному приближению
- •Условие моноэнергетичности
- •Условие изотропности источника и рассеяния
- •Закон Фика и диффузионное приближение
- •Основное приближение диффузии нейтронов
- •Физический смысл закона Фика
- •Условия применимости диффузионного приближения
- •Основные условия
- •Большая среда
- •Слабопоглощающая среда
- •Гомогенность среды
- •Далеко от локальных неоднородностей
- •Транспортная поправка
- •Проблема анизотропии рассеяния
- •Введение транспортной поправки
- •Уравнение диффузии
- •Стационарное уравнение диффузии в гомогенной среде
- •Физический смысл длины диффузии
- •Решение уравнения диффузии
- •Условия однозначного выбора решения уравнения диффузии
- •Основные условия
- •Условия на границе раздела двух сред
- •Экстраполированная граница
- •Условия для локализованного источника
- •Методы решения уравнения диффузии
- •Аналитические решения
- •Численные методы
- •Метод конечных разностей
- •Синтетические методы
5.2Физический смысл длины диффузии
(38:55) Длина диффузии L имеет четкий физический смысл - это среднее квадратичное смещение нейтрона от точки рождения до точки поглощения.
(39:37) Если нейтрон родился в одной точке, совершил несколько столкновений и был поглощен в другой точке, то средний квадрат расстояния между этими точками связан с
L2.
5.3Решение уравнения диффузии
(40:29) Уравнение диффузии является дифференциальным уравнением второго порядка. Для его решения необходимо задать граничные и начальные условия.
(41:09) В общем случае уравнение диффузии может быть неоднородным и в частных производных. Такие уравнения решаются известными методами, изучаемыми в курсах математики.
(42:35) После формулировки уравнения диффузии физическая часть задачи завершается, и начинается математическая часть - решение дифференциального уравнения.
(43:10) При решении дифференциального уравнения второго порядка в решении будут две неизвестные константы. Для их определения необходимы граничные и начальные условия.
6Условия однозначного выбора решения уравнения диффузии
6.1Основные условия
(51:20) Для однозначного выбора решения уравнения диффузии необходимо выполнение следующих условий:
1.Функция потока является неотрицательной и ограниченной: 0 ≤ Φ(r, t) < ∞
2.Для нестационарной задачи функция потока в начальный момент времени известна:
Φ(r, t = 0) = f(r)
3.Условия на границе раздела двух сред
4.Условия на невогнутой границе среда-вакуум
5.Условия для локализованного источника
6.2Условия на границе раздела двух сред
(53:29) На границе раздела двух сред должны выполняться условия равенства односторонних токов:
(
+ |
+ |
(rs) |
|
|
in(1) |
(rs) = in(2) |
(9.6.1) |
||
i− |
(rs) = i− (rs) |
|||
|
||||
n(1) |
n(2) |
|
|
|
6
где rs - точка на границе раздела.
(54:45) Из этих условий следует, что на границе раздела двух сред поток нейтронов непрерывен:
Φ1(rs) = Φ2(rs) |
(9.6.2) |
(54:55) А производная потока терпит разрыв на величину, зависящую от коэффициентов диффузии:
D1 nΦ1(rs) = D2 nΦ2(rs) |
(9.6.3) |
(57:19) В реальности функция потока не имеет изломов на границе, а ведет себя плавно. Излом в диффузионном приближении связан с тем, что на границе среды диффузия работает плохо.
6.3Условия на невогнутой границе среда-вакуум
(57:50) На невогнутой границе среда-вакуум должно выполняться условие отсутствия возвратного потока нейтронов:
in−(rs) = 0 |
(9.6.4) |
(1:00:15) С учетом выражения для одностороннего тока это условие можно переписать в виде:
|
|
Φ(rs) |
D |
nΦ(rs) = 0 |
|
|
(9.6.5) |
|||||||||||
|
|
|
|
|
|
+ |
|
|
|
|
||||||||
|
|
|
|
4 |
|
2 |
|
|
|
|||||||||
(1:00:43) Преобразуем это выражение: |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
nΦ(rs) |
= |
− |
1 |
= |
− |
3Σtr |
= |
− |
1 |
= |
− |
1 |
(9.6.6) |
||||
|
|
|
|
|
|
|||||||||||||
|
Φ(rs) |
|
|
2D |
|
2 |
|
32 λtr |
α |
|||||||||
где α - длина линейной экстраполяции потока в вакуум.
6.3.1Экстраполированная граница
(1:08:30) Для упрощения решения дифференциального уравнения вводится понятие экстраполированной границы. Рассмотрим линейную аппроксимацию потока вблизи границы:
|
|
|
Φ(x) = C1x + C2 |
|
|
|
(9.6.7) |
||||
(1:09:05) На экстраполированной границе x = x0 поток обращается в ноль: |
|
||||||||||
|
C1x0 + C2 = 0 C2 = −C1x0 |
|
|
|
(9.6.8) |
||||||
(1:09:45) Подставим эти выражения в формулу для граничного условия: |
|
||||||||||
|
Φ′(x = xs) |
|
C1 |
1 |
|
1 |
|
(9.6.9) |
|||
|
|
|
= |
|
= − |
|
|
≡ − |
|
|
|
|
Φ(x = xs) |
C1xs + C2 |
x0 − xs |
α |
|||||||
7
