Уравнение переноса излучения¶

Уравнение переноса излучения (УПИ) — это интегро-дифференциальное уравнение, описывающее изменение интенсивности электромагнитного излучения (или других частиц, например, нейтронов) при его распространении через среду, которая может поглощать, испускать и рассеивать это излучение. Уравнение является фундаментальным для физической кинетики, астрофизики, теории переноса нейтронов, атмосферной оптики и других областей науки и техники, где необходимо моделировать взаимодействие излучения с веществом.

¶История

Формирование теории переноса излучения началось в конце XIX — начале XX века. Первоначально задача ставилась в рамках астрофизики для описания переноса энергии из недр звезд к их поверхности. В 1906 году немецкий физик Карл Шварцшильд ввел понятие лучистого равновесия и сформулировал основные уравнения переноса в звездных атмосферах. В 1920-х годах Артур Эддингтон (Великобритания) развил приближенные методы решения, в частности, метод диффузионного приближения и ввел понятие «эддингтоновского предела» светимости.

В 1940-х годах, в связи с разработкой ядерного оружия и ядерных реакторов, возникла необходимость в описании переноса нейтронов. Это привело к созданию строгой математической теории переноса, в которой уравнение переноса излучения (в данном случае — нейтронов) стало центральным. Значительный вклад внесли советские ученые: В. С. Владимиров, А. А. Самарский, а также американские ученые — С. Чандрасекар (впоследствии нобелевский лауреат), который в 1950 году опубликовал монографию «Перенос лучистой энергии», ставшую классической.

¶Основное уравнение

Уравнение переноса излучения в его наиболее общей форме записывается для спектральной интенсивности излучения \( I_\nu(\mathbf{r}, \mathbf{\Omega}, t) \), которая зависит от пространственных координат \( \mathbf{r} \), направления распространения \( \mathbf{\Omega} \) (единичный вектор), частоты \( \nu \) и времени \( t \). В стационарном случае (без зависимости от времени) оно имеет вид:

\[ \mathbf{\Omega} \cdot \nabla I_\nu = \varepsilon_\nu - \kappa_\nu I_\nu - \sigma_\nu I_\nu + \frac{\sigma_\nu}{4\pi} \int_{4\pi} p(\mathbf{\Omega}', \mathbf{\Omega}) I_\nu(\mathbf{\Omega}') d\mathbf{\Omega}' \]

где:

• \(\mathbf{\Omega} \cdot \nabla I_\nu\) — изменение интенсивности вдоль направления \(\mathbf{\Omega}\) (конвективная производная);
• \(\varepsilon_\nu\) — коэффициент излучения (эмиссии) среды;
• \(\kappa_\nu\) — коэффициент истинного поглощения;
• \(\sigma_\nu\) — коэффициент рассеяния;
• \(p(\mathbf{\Omega}', \mathbf{\Omega})\) — индикатриса рассеяния, описывающая вероятность перехода излучения из направления \(\mathbf{\Omega}'\) в направление \(\mathbf{\Omega}\).

В более компактной форме, вводя полный коэффициент ослабления \(\chi_\nu = \kappa_\nu + \sigma_\nu\) и функцию источника \(S_\nu\), уравнение записывают как:

\[ \mathbf{\Omega} \cdot \nabla I_\nu = \chi_\nu (S_\nu - I_\nu) \]

¶Физический смысл членов

Уравнение представляет собой баланс энергии в единице объема в единицу времени. Левая часть описывает изменение интенсивности за счет геометрического переноса (поток через границы элементарного объема). Правая часть состоит из трех типов процессов:

• Поглощение: излучение переходит в другие формы энергии (например, тепловую), интенсивность уменьшается.
• Излучение: среда сама испускает фотоны (например, за счет теплового движения атомов), интенсивность увеличивается.
• Рассеяние: фотоны меняют направление, при этом часть излучения уходит из данного направления, а часть приходит из других направлений (интегральный член).

¶Классификация сред и режимов

УПИ можно классифицировать по свойствам среды и характеру взаимодействия:

¶По оптической толщине

• Оптически тонкая среда: среда прозрачна, фотоны проходят без значительного поглощения или рассеяния. Уравнение упрощается до \( I_\nu \approx \int \varepsilon_\nu dl \).
• Оптически толстая среда: среда непрозрачна, излучение распространяется путем диффузии. В этом случае применимо диффузионное приближение.

¶По типу рассеяния

• Изотропное рассеяние: индикатриса \( p = 1 \) (равновероятно во все стороны).
• Анизотропное рассеяние: фотоны преимущественно рассеиваются вперед или назад (например, рассеяние Ми на крупных частицах).

¶По частоте

• Монохроматическое УПИ: рассматривается одна частота.
• Многогрупповое УПИ: спектр разбивается на интервалы (группы), для каждой группы решается свое уравнение.

¶Методы решения

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

¶Метод сферических гармоник (P_N-приближение)

Интенсивность раскладывается в ряд по сферическим функциям, что сводит УПИ к системе дифференциальных уравнений в частных производных. Наиболее популярно P_1-приближение (диффузионное), которое дает хорошие результаты для оптически толстых сред.

¶Метод дискретных ординат (S_N-приближение)

Пространство направлений \(\mathbf{\Omega}\) дискретизируется на конечное число узлов (ординат). УПИ превращается в систему уравнений переноса вдоль фиксированных направлений. Метод широко используется в кодах для расчета ядерных реакторов (например, в кодах MCNP, OpenMC).

¶Метод Монте-Карло

Статистический метод, моделирующий траектории большого числа фотонов (или нейтронов). Позволяет решать задачи в сложной геометрии с высокой точностью, но требует больших вычислительных ресурсов. Применяется в астрофизике (например, код RADMC-3D) и в ядерной физике.

¶Метод характеристик

УПИ интегрируется вдоль характеристик (лучей) в заданном направлении. Эффективен для двумерных и трехмерных задач, особенно в сочетании с методом конечных элементов.

¶Применение

¶Астрофизика

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

¶Ядерная физика и ядерная энергетика

В расчетах ядерных реакторов УПИ описывает перенос нейтронов. На его основе проектируются активные зоны, защита от излучения, системы управления. Методы решения УПИ лежат в основе кодов, таких как MCNP (США) и MCU (Россия).

¶Атмосферная оптика и климатология

УПИ моделирует распространение солнечного излучения в атмосфере Земли, учитывая рассеяние на молекулах (рэлеевское рассеяние) и аэрозолях, а также поглощение газами (озон, водяной пар, углекислый газ). Это необходимо для прогноза погоды, расчета радиационного баланса и оценки климатических изменений.

¶Медицинская физика

В лучевой терапии и диагностике (например, компьютерная томография) УПИ описывает перенос рентгеновского излучения в тканях. Это позволяет оптимизировать дозы облучения и улучшить качество изображений.

¶Оптика и фотоника

УПИ используется для расчета светораспределения в светодиодах, оптических волокнах, фотокаталитических реакторах, а также в задачах визуализации (например, в биомедицинской оптике для томографии).

¶Интересные факты

• Уравнение переноса излучения является частным случаем кинетического уравнения Больцмана для фотонов (или нейтронов), в котором пренебрегается столкновениями между частицами (так как фотоны не взаимодействуют друг с другом).
• В астрофизике существует понятие «лучистой теплопроводности» — это диффузионное приближение УПИ для оптически толстых сред, где перенос энергии описывается как теплопроводность с эффективным коэффициентом.
• Метод Монте-Карло для решения УПИ впервые был применен в 1940-х годах в Лос-Аламосской национальной лаборатории (США) в рамках Манхэттенского проекта.
• В 2023 году российские ученые из Института прикладной математики им. М. В. Келдыша РАН разработали новый численный метод решения УПИ на основе нейросетей, позволяющий ускорить расчеты в 10-100 раз по сравнению с классическими методами.

¶Критика и ограничения

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