Функция Грина для уравнения Лапласа: смысл и построение

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

Что задаёт функция Грина

Функция Грина $G(\mathbf{r}, \mathbf{r}')$ для оператора Лапласа в области $\Omega$ определяется как решение уравнения с дельта-источником в точке $\mathbf{r}'$:

$$\Delta_{\mathbf{r}} G(\mathbf{r}, \mathbf{r}') = -\delta(\mathbf{r} - \mathbf{r}')$$

с обязательным граничным условием. Для задачи Дирихле требуют, чтобы на границе функция обращалась в ноль:

$$G(\mathbf{r}, \mathbf{r}') = 0, \quad \mathbf{r} \in \partial\Omega.$$

Физически $G$ - это потенциал в точке $\mathbf{r}$, созданный единичным точечным зарядом в $\mathbf{r}'$, при условии что граница области заземлена. Минус перед дельтой - вопрос соглашения и согласован с тем, как пишут само уравнение Пуассона: источник плотности $\rho$ даёт $\Delta u = -\rho/\varepsilon_0$.

Рассчитать для своей задачи

Ключевое отличие от свободного поля: в свободном пространстве потенциал точечного заряда задаётся только фундаментальным решением. В ограниченной области к нему добавляется гладкая поправка, которая «гасит» потенциал на границе. Эта поправка и делает $G$ привязанной к конкретной геометрии.

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

Сердцевина любой функции Грина - фундаментальное решение оператора Лапласа, то есть потенциал точечного источника в неограниченном пространстве. В трёх измерениях это

$$G_0(\mathbf{r}, \mathbf{r}') = \frac{1}{4\pi |\mathbf{r} - \mathbf{r}'|},$$

а в двумерном случае - логарифм:

$$G_0(\mathbf{r}, \mathbf{r}') = -\frac{1}{2\pi} \ln |\mathbf{r} - \mathbf{r}'|.$$

Эти выражения удовлетворяют уравнению с дельта-источником всюду, кроме самой точки $\mathbf{r}'$, где есть особенность. Проверка - стандартная: интеграл от $\Delta G_0$ по малому шару вокруг источника сводится через теорему Гаусса к потоку градиента, и он равен единице независимо от радиуса.

Полная функция Грина области строится как сумма

$$G(\mathbf{r}, \mathbf{r}') = G_0(\mathbf{r}, \mathbf{r}') + H(\mathbf{r}, \mathbf{r}'),$$

где $H$ - гармоническая в $\Omega$ функция (то есть $\Delta H = 0$), подобранная так, чтобы на границе сумма обнулилась. Особенность целиком сидит в $G_0$, а $H$ гладкая - это разделение и упрощает построение.

Метод отражений для простых областей

Когда граница - плоскость или сфера, гладкую поправку $H$ удаётся угадать геометрически: это потенциал фиктивного заряда-отражения, размещённого вне области. Приём называют методом отражений (методом изображений), и он напрямую связан с электростатикой заземлённого проводника.

Рассчитать для своей задачи

Для полупространства $z > 0$ с границей-плоскостью $z = 0$ зеркальный заряд противоположного знака ставят в точку, симметричную $\mathbf{r}'$ относительно плоскости. Тогда

$$G(\mathbf{r}, \mathbf{r}') = \frac{1}{4\pi}\left( \frac{1}{|\mathbf{r} - \mathbf{r}'|} - \frac{1}{|\mathbf{r} - \mathbf{r}'^*|} \right),$$

где $\mathbf{r}'^*$ - отражённая точка. На самой плоскости расстояния до реального и мнимого зарядов совпадают, и разность зануляется автоматически.

Для шара радиуса $R$ образ сложнее: заряд-изображение ставится в инверсную точку $\mathbf{r}'^* = R^2 \mathbf{r}' / |\mathbf{r}'|^2$ и берётся с весом $R/|\mathbf{r}'|$. Этот подбор - геометрическая инверсия, тот же приём, что стоит за решением многих краевых задач для круга и шара. Логику инверсии полезно сравнить с тем, как устроен оператор Лапласа-Бельтрами на искривлённых поверхностях, где роль плоского лапласиана меняется.

Симметрия и физический смысл

Важное свойство функции Грина оператора Лапласа - симметрия по аргументам:

$$G(\mathbf{r}, \mathbf{r}') = G(\mathbf{r}', \mathbf{r}).$$

Это принцип взаимности: потенциал в точке $\mathbf{r}$ от источника в $\mathbf{r}'$ равен потенциалу в $\mathbf{r}'$ от источника в $\mathbf{r}$. Доказывается через вторую формулу Грина, применённую к двум функциям $G(\cdot, \mathbf{r}_1)$ и $G(\cdot, \mathbf{r}_2)$: граничные интегралы исчезают благодаря условию Дирихле, и остаётся равенство значений.

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

Формула решения через функцию Грина

Главная польза функции Грина - готовая формула для решения краевой задачи Дирихле. Если внутри области задан источник $f$ и на границе - значения $g$, то

$$u(\mathbf{r}) = \int_{\Omega} G(\mathbf{r}, \mathbf{r}')\, f(\mathbf{r}')\, dV' - \int_{\partial\Omega} \frac{\partial G}{\partial n'}(\mathbf{r}, \mathbf{r}')\, g(\mathbf{r}')\, dS'.$$

Первый интеграл - вклад объёмных источников, второй - вклад граничных значений через нормальную производную функции Грина. Производную $\partial G / \partial n'$ называют ядром Пуассона.

Рассчитать для своей задачи

Для шара эта схема даёт классическую формулу Пуассона, выражающую гармоническую функцию внутри через её значения на сфере. Так функция Грина превращает дифференциальную задачу в интегральную - это тот же сдвиг точки зрения, что и при переходе к интегральному уравнению Фредгольма второго рода, где неизвестное тоже ищут под знаком интеграла.

Двумерный случай и конформные отображения

В плоскости задачу часто удаётся свести к уже решённой геометрии конформным отображением. Если найдено аналитическое отображение $w = \varphi(z)$, переводящее область $\Omega$ в единичный круг, то функция Грина области выражается через расстояние образов источника и точки наблюдения внутри круга. Для верхней полуплоскости, например, отображение на круг даёт явную формулу с логарифмом отношения модулей.

Такая техника работает только в двумерии, потому что лапласиан инвариантен относительно конформных замен лишь на плоскости: гармоническая функция остаётся гармонической после подстановки $w = \varphi(z)$. Это и делает теорию потенциала в 2D особенно богатой - круг, полуплоскость, полоса и угол связаны простыми отображениями, а значит, имеют функцию Грина в замкнутой форме. В трёх измерениях такой свободы нет, и для несимметричных областей приходится переходить к численным методам или разложениям по собственным функциям.

Связь с уравнением Пуассона

Уравнение Лапласа $\Delta u = 0$ - частный случай уравнения Пуассона $\Delta u = -\rho/\varepsilon_0$ при нулевой правой части. Функция Грина строится именно для оператора, а не для конкретной правой части, поэтому одна и та же $G$ обслуживает оба уравнения. Когда источников внутри нет, объёмный интеграл исчезает и решение полностью определяется граничными данными - так гармоническая функция восстанавливается по краю.

Это объясняет, почему функция Грина так ценна в теории потенциала: она отделяет геометрию области (зашита в $G$) от конкретного распределения зарядов и граничных значений (входят как множители под интегралом). Сменился источник или граничный профиль - решение пересобирается без повторного решения дифференциального уравнения. Удобно держать рядом смысл градиента потенциала электрического поля: именно его восстанавливают, дифференцируя найденный потенциал.

Частые ошибки

• Путают функцию Грина с фундаментальным решением. $G_0$ удовлетворяет уравнению с дельтой, но не обнуляется на границе. Полная $G$ - это $G_0$ плюс гармоническая поправка под конкретную область.
• Теряют знак при определении. Соглашение $\Delta G = -\delta$ должно быть согласовано со знаком в уравнении Пуассона и с направлением нормали в граничном интеграле; иначе решение получит неверный знак.
• Применяют метод отражений к произвольной границе. Зеркальный заряд работает только для плоскости и сферы (и их комбинаций). Для эллипса или прямоугольника простого образа нет - нужны ряды или конформные отображения.
• Берут наружную нормаль не там. В формуле Пуассона производная $\partial G/\partial n'$ берётся по точке источника на границе и по внешней нормали; перепутанная нормаль даёт знаковую ошибку в поверхностном интеграле.
• Забывают про размерность. В 2D фундаментальное решение логарифмическое, а не $1/r$; формулы из трёхмерного случая нельзя переносить буквально.

FAQ

Чем функция Грина отличается от потенциала точечного заряда? Потенциал точечного заряда в свободном пространстве - это фундаментальное решение $G_0$. Функция Грина области добавляет к нему гладкую поправку, которая обнуляет потенциал на границе. То есть $G$ учитывает геометрию области, а $G_0$ - нет.

Зачем нужно граничное условие в определении? Без граничного условия задача определена неоднозначно: к решению можно добавить любую гармоническую функцию. Условие $G = 0$ на границе (для Дирихле) фиксирует единственную функцию Грина и делает формулу решения корректной.

Можно ли построить функцию Грина для любой области? Теоретически да - она существует для широкого класса областей. Но явное выражение получается только для симметричных случаев (полупространство, шар, полуплоскость). Для остальных её ищут численно или через конформные отображения в 2D.

Коротко

Функция Грина для уравнения Лапласа - это потенциал единичного точечного источника в области с заземлённой границей. Она равна фундаментальному решению плюс гармоническая поправка, обнуляющая потенциал на границе; для плоскости и сферы поправку даёт метод отражений. Функция Грина симметрична по аргументам, не зависит от конкретного источника и через объёмный и поверхностный интегралы сразу выдаёт решение краевой задачи Дирихле, превращая дифференциальную постановку в интегральную.