Уравнение Пуассона

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

Определение и общий вид

В трёхмерном пространстве уравнение Пуассона записывается как:

\[ \nabla^2 \varphi = f(\mathbf{r}) \]

или в развёрнутой форме:

\[ \frac{\partial^2 \varphi}{\partial x^2} + \frac{\partial^2 \varphi}{\partial y^2} + \frac{\partial^2 \varphi}{\partial z^2} = f(x,y,z) \]

где:

\(\varphi\) — искомая скалярная функция (потенциал),
\(\nabla^2\) — оператор Лапласа (лапласиан),
\(f(\mathbf{r})\) — заданная функция, описывающая плотность источников поля.

В случае, когда \(f(\mathbf{r}) = 0\), уравнение Пуассона переходит в уравнение Лапласа \(\nabla^2 \varphi = 0\), описывающее поле в области без источников.

Физический смысл

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

Электростатика

В электростатике уравнение Пуассона для электрического потенциала \(\varphi\) имеет вид:

\[ \nabla^2 \varphi = -\frac{\rho}{\varepsilon_0} \]

где \(\rho\) — объёмная плотность электрического заряда, \(\varepsilon_0\) — электрическая постоянная. Это уравнение является дифференциальной формой закона Кулона и теоремы Гаусса.

Гравитация

В ньютоновской теории гравитации гравитационный потенциал \(\Phi\) удовлетворяет уравнению:

\[ \nabla^2 \Phi = 4\pi G \rho \]

где \(G\) — гравитационная постоянная, \(\rho\) — плотность массы. Это уравнение лежит в основе теории гравитационного поля.

Термодинамика и теплопередача

В стационарной теплопроводности уравнение Пуассона описывает распределение температуры \(T\) при наличии внутренних источников тепла:

\[ \nabla^2 T = -\frac{q}{\kappa} \]

где \(q\) — плотность тепловыделения, \(\kappa\) — коэффициент теплопроводности.

Гидродинамика

В теории потенциальных течений несжимаемой жидкости уравнение Пуассона связывает потенциал скорости \(\phi\) с дивергенцией источника:

\[ \nabla^2 \phi = \theta \]

где \(\theta\) — объёмная скорость расширения жидкости.

История

Уравнение названо в честь французского математика и физика Симеона Дени Пуассона (1781–1840), который впервые вывел его в 1813 году в контексте электростатики. Пуассон обобщил уравнение Лапласа, введя в правую часть ненулевую функцию, соответствующую плотности заряда. Ранее, в 1782 году, Пьер-Симон Лаплас сформулировал уравнение для потенциала в пустоте, а Пуассон распространил его на области с источниками.

Математические свойства

Линейность

Уравнение Пуассона является линейным: если \(\varphi_1\) и \(\varphi_2\) — решения с правыми частями \(f_1\) и \(f_2\), то \(\varphi_1 + \varphi_2\) — решение с правой частью \(f_1 + f_2\). Это позволяет использовать принцип суперпозиции.

Тип уравнения

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

Краевые задачи

Для однозначного определения решения уравнения Пуассона необходимо задать граничные условия на границе области \(\Omega\):

Задача Дирихле: задано значение \(\varphi\) на границе \(\partial\Omega\).
Задача Неймана: задана нормальная производная \(\frac{\partial\varphi}{\partial n}\) на границе.
Смешанная задача: комбинация условий Дирихле и Неймана на разных участках границы.

Фундаментальное решение

Фундаментальное решение уравнения Пуассона в трёхмерном пространстве имеет вид:

\[ G(\mathbf{r}) = -\frac{1}{4\pi |\mathbf{r}|} \]

В двумерном случае фундаментальное решение выражается через логарифм:

\[ G(\mathbf{r}) = \frac{1}{2\pi} \ln|\mathbf{r}| \]

Используя фундаментальное решение, решение уравнения Пуассона в неограниченном пространстве можно записать в виде интеграла:

\[ \varphi(\mathbf{r}) = \int_{\mathbb{R}^3} G(\mathbf{r} - \mathbf{r}') f(\mathbf{r}') \, d^3\mathbf{r}' \]

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

Аналитические методы

Метод функции Грина: построение решения с помощью функции Грина для заданной области и граничных условий.
Метод разделения переменных: применим для областей простой формы (прямоугольник, круг, сфера) с использованием рядов Фурье, полиномов Лежандра или сферических функций.
Метод интегральных преобразований: преобразование Фурье или Лапласа для перехода к алгебраическому уравнению.

Численные методы

Метод конечных разностей: аппроксимация производных разностными схемами на равномерной сетке.
Метод конечных элементов: разбиение области на треугольники или тетраэдры и построение приближённого решения.
Метод граничных элементов: сведение задачи к интегральному уравнению на границе области.
Метод Монте-Карло: вероятностное решение с использованием случайных блужданий.

Применение

Электростатика и электротехника

Расчёт электрических полей в конденсаторах, линиях электропередач, электронных приборах.
Моделирование распределения потенциала в полупроводниковых структурах (p-n-переходы, транзисторы).
Проектирование электростатических линз и ускорителей частиц.

Гравитационное моделирование

Расчёт гравитационного поля Земли и других планет.
Моделирование распределения массы внутри небесных тел.
Геофизическая разведка: интерпретация гравитационных аномалий.

Теплопередача

Расчёт стационарного температурного поля в твёрдых телах с внутренними источниками тепла.
Проектирование систем охлаждения электроники.
Тепловой расчёт ядерных реакторов.

Гидродинамика и аэродинамика

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

Астрофизика

Моделирование гравитационного потенциала звёзд и галактик.
Расчёт структуры протопланетных дисков.
Теория гравитационной линзы.

Связь с другими уравнениями

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

Уравнение Лапласа — частный случай при \(f=0\).
Уравнение Гельмгольца \(\nabla^2 \varphi + k^2 \varphi = 0\) — обобщение для волновых процессов.
Уравнение теплопроводности \(\frac{\partial u}{\partial t} = a \nabla^2 u\) — нестационарный аналог.
Уравнение Пуассона — Больцмана — обобщение для электролитов с учётом распределения ионов.

Обобщения

Нелинейное уравнение Пуассона

В некоторых физических моделях (например, в теории электролитов или полупроводников) правая часть может зависеть от самого потенциала:

\[ \nabla^2 \varphi = F(\varphi) \]