Гармонические функции и задача Дирихле

В 1791 году Шарль-Огюстен де Кулон записал закон взаимодействия зарядов; через несколько десятилетий Лаплас и Пуассон поняли, что потенциал точечных зарядов в свободном от зарядов пространстве удовлетворяет $\Delta u = 0$. Сегодня каждая симуляция распределения температуры в COMSOL, каждый расчёт электрического поля микросхемы в Cadence, каждая параметризация UV-карт текстур в Maya - это численное решение задачи Дирихле. От гравитационного потенциала Солнца до изгиба обтекающего профиль воздушного потока - один и тот же оператор $\Delta$, одни и те же граничные условия.

• **Электростатика:** Cadence Virtuoso и COMSOL Electrostatics решают $\Delta V = 0$ в зазорах между проводниками микросхемы
• **Стационарная теплопроводность:** ANSYS Steady-State Thermal моделирует распределение $T(x,y,z)$ в радиаторах и теплоотводах
• **3D-графика (UV-mapping):** алгоритмы Floater и LSCM строят гармонические отображения mesh на плоскость для текстурирования
• **Метод Брауна (1942):** решение задачи Дирихле через броуновское движение - $u(x_0) = \mathbb{E}[f(B_\tau)]$, где $\tau$ - время выхода

Предварительные знания

• Аналитическое продолжение

Гармонические функции

Гармоническая функция - вещественная $u(x, y)$, удовлетворяющая уравнению Лапласа $\Delta u = 0$. Уравнение появилось в 1782 году у Пьера-Симона Лапласа при анализе гравитационного потенциала Солнечной системы: в вакууме, между притягивающими массами, потенциал гармоничен. С тех пор это уравнение собрало под одной крышей электростатику, стационарную теплопроводность, потенциальные течения и - через КР - всю аналитическую функцию.

**Уравнение Лапласа:** $$\Delta u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0$$ **Связь с КА:** если $f = u + iv$ аналитична, то $u, v$ - гармонические и связаны уравнениями Коши-Римана: $$\frac{\partial u}{\partial x} = \frac{\partial v}{\partial y}, \quad \frac{\partial u}{\partial y} = -\frac{\partial v}{\partial x}$$ **Сопряжённая функция:** если $u$ гармонична в односвязной области, то существует $v$ такая, что $f = u + iv$ аналитична.

Теорема о среднем и принцип максимума

Гармонические функции обладают двумя поразительными свойствами: значение в центре круга равно среднему по граничной окружности, а максимум всегда достигается на границе. Оба факта - проявление сглаживающего характера оператора Лапласа: $\Delta u = 0$ означает, что $u$ не может локально 'выгибаться'.

**Теорема о среднем значении:** $$u(z_0) = \frac{1}{2\pi} \int_0^{2\pi} u(z_0 + r e^{i\theta})\, d\theta$$ **Принцип максимума:** гармоническая функция в области $D$ не может достигать локального максимума (или минимума) во внутренней точке, если она не константа. **Следствие (единственность):** разность двух решений задачи Дирихле гармонична и обращается в нуль на границе, значит тождественно равна нулю внутри.

Интегральная формула Пуассона

Симеон Дени Пуассон в 1820 году решил задачу Дирихле для круга в явном виде. Формула воспроизводит гармоническую функцию внутри круга по её значениям на границе - это правый ответ на вопрос 'если знаю температуру по периметру диска, что внутри?'. Ядро формулы - $P(r, \theta) = (R^2 - r^2)/(R^2 - 2Rr\cos\theta + r^2)$ - неотрицательно и нормировано: его легко перенести в численные методы и в теорию вероятностей (там же он появляется как переходная плотность броуновского движения).

**Формула Пуассона для круга $|z| < R$:** $$u(r, \theta) = \frac{1}{2\pi} \int_0^{2\pi} \frac{R^2 - r^2}{R^2 - 2Rr\cos(\theta - \varphi) + r^2}\, f(\varphi)\, d\varphi$$ где $f(\varphi) = u(R, \varphi)$ - граничные значения. **Ядро Пуассона:** $P(r, \theta - \varphi) = \dfrac{R^2 - r^2}{R^2 - 2Rr\cos(\theta - \varphi) + r^2} \geq 0$, и $\int_0^{2\pi} P\, d\varphi = 2\pi$. **Свойства:** при $r \to R$ ядро стягивается в дельта-функцию на границе - значения $u$ приближаются к $f$.

Задача Дирихле

Квадратная металлическая пластина 1×1 м: три стороны держим при $0\,^\circ$C, одну при $100\,^\circ$C. После выхода на стационарный режим температура $u(x, y)$ в каждой внутренней точке удовлетворяет $\Delta u = 0$ и совпадает с граничными значениями. Найти $u$ во всех внутренних точках - это задача Дирихле. То же уравнение описывает потенциал между проводниками в электростатике, потенциал скорости в потенциальном течении и решение броуновской задачи о выходе из области.

**Постановка задачи Дирихле:** $$\Delta u = 0 \text{ в } D, \quad u|_{\partial D} = f$$ **Единственность:** по принципу максимума разность двух решений тождественно нулевая. **Существование (метод Перрона, 1923):** строим супергармонические барьеры, берём infimum - получаем решение в любой 'регулярной' области. **Альтернатива:** формула Пуассона для круга плюс конформное отображение области на круг - получаем явное решение в широком классе областей.

Ключевые идеи

• **Гармоническая функция:** $\Delta u = 0$; вещественная и мнимая часть аналитической $f = u + iv$ - сопряжённые гармонические
• **Теорема о среднем:** $u(z_0)$ равно среднему по окружности вокруг $z_0$ - объяснение принципа максимума
• **Формула Пуассона:** явное решение задачи Дирихле для круга через интеграл с ядром $P(r, \theta - \varphi)$
• **Задача Дирихле:** $\Delta u = 0$ плюс граничные значения дают единственное решение в регулярной области

Связанные темы

Гармонические функции - мост между комплексным анализом и математической физикой:

• Комплексный анализ в физике — Гармонические функции - это потенциалы электростатики, теплопроводности, потенциальных течений
• Конформные отображения — Конформные отображения сохраняют гармоничность, позволяя переносить задачу Дирихле из сложной области в круг
• Аналитическое продолжение — Гармонические функции продолжаются по принципу симметрии Шварца - прямой аналог аналитического продолжения

Вопросы для размышления

• Почему теорема о среднем мгновенно даёт принцип максимума? Где именно используется неравенство среднего?
• Как конформное отображение области $D$ на единичный круг превращает задачу Дирихле в готовый интеграл Пуассона? Что сохраняется при отображении, а что меняется?
• Гармонические функции и мартингалы: значение $u$ в точке - это математическое ожидание $u$ от броуновского движения, выходящего из этой точки. Как это объясняет принцип максимума на вероятностном языке?

Связанные уроки

• calc-17-multivariable