Дифференциальные уравнения
Уравнение Лапласа
1828 год: Джордж Грин - мельник без университетского образования - публикует работу об электростатике и доказывает, что любое стационарное поле полностью определяется граничными значениями. Никакой информации о внутренности. Это уравнение Лапласа: ∇²u = 0, и максимум всегда на границе, никогда внутри.
- **Электростатическое моделирование:** распределение потенциала в полупроводниковых устройствах - прямое решение ∇²V = -ρ/ε. Используется в симуляторах чипов (SPICE, COMSOL).
- **Потенциальное обтекание:** в аэродинамике при малых скоростях поток вокруг крыла описывается уравнением Лапласа. Формула Жуковского для подъёмной силы - через комплексный потенциал.
- **Метод граничных элементов (BEM):** решение задач Лапласа только на границе (не во всём объёме) через функции Грина. Используется в акустике, электромагнетизме и МЭМС.
Предварительные знания
Уравнение Лапласа и гармонические функции
**Уравнение Лапласа** ∇²u = 0 описывает стационарные состояния: установившееся распределение температуры (∂T/∂t = 0), электростатический потенциал (нет заряда), несжимаемое безвихревое течение. Решения называются **гармоническими функциями** и обладают замечательными свойствами.
**Принцип максимума:** гармоническая функция не может достигать максимума (или минимума) внутри области - только на границе. Следствие: если ∇²u = 0 и u = 0 на границе, то u ≡ 0 внутри. Это гарантирует единственность решения.
| Свойство | Формулировка |
|---|---|
| Принцип максимума | max u и min u достигаются на границе |
| Среднее значение | u(x₀) = среднее u по сфере вокруг x₀ |
| Единственность | Задача Дирихле имеет единственное решение |
| Гладкость | ∇²u = 0 → u бесконечно дифференцируема |
| Связь с аналитич. функциями | Re(f(z)) и Im(f(z)) гармоничны |
Квадратная пластина с температурами: верх = 100°, остальные стороны = 0°. Какова температура в центре (принцип среднего)?
Уравнение Пуассона
**Уравнение Пуассона** ∇²u = f(x) - уравнение Лапласа с правой частью. Физика: f = -ρ/ε₀ (электростатика: плотность заряда), f = -q (гидростатика: источник), f = -μJ (магнетизм). Решение = гармоническая часть + конкретное решение.
**Связь с комплексным анализом:** если w = u + iv - аналитическая функция, то u и v - гармоничны и связаны условиями Коши-Римана. Потенциальное течение жидкости описывается аналитической функцией: u - потенциал скорости, v - функция тока.
Уравнение Пуассона ∇²u = f. Если решение задачи Лапласа (f=0) единственно, единственно ли решение задачи Пуассона?
Функции Грина
**Функция Грина G(x, y)** - это «отклик» уравнения Лапласа на точечный источник δ(x-y): ∇²G(x,y) = δ(x-y). Зная G, решение уравнения Пуассона ∇²u = f записывается как: **u(x) = ∫ G(x,y) f(y) dy**. Это аналог передаточной функции - но для PDE.
| Измерение | G(x,y) | Поведение |
|---|---|---|
| 1D: G(x,y) | -|x-y|/2 | Кусочно-линейная |
| 2D: G(x,y) | (1/2π) ln|x-y| | Логарифмическая сингулярность |
| 3D: G(x,y) | -1/(4π|x-y|) | Кулоновский потенциал 1/r |
В 2D функция Грина G = (1/2π)ln|x-y|. Как убывает потенциал одиночного заряда в 2D при r → ∞?
Численное решение: метод Гаусса-Зейделя
На дискретной сетке Δx = Δy = h, уравнение Лапласа ∇²u = 0 аппроксимируется: **(u[i+1,j] + u[i-1,j] + u[i,j+1] + u[i,j-1] - 4u[i,j]) / h² = 0**, откуда **u[i,j] = (u[i+1,j] + u[i-1,j] + u[i,j+1] + u[i,j-1]) / 4**. Метод Гаусса-Зейделя: итерировать это правило до сходимости.
**Скорость сходимости:** метод Якоби сходится медленно (O(h⁻²) итераций). Ускорение: **SOR** (последовательная релаксация), **многосеточный метод** (multigrid, сходимость за O(N log N)), **FFT-метод** (аналитическое решение за O(N log N)).
Джордж Грин - мельник, который изобрёл функцию Грина (1828)
Грин самостоятельно учился математике, работая на мельнице отца в Ноттингеме. В 1828 году он за свой счёт опубликовал «Essay on the Application of Mathematical Analysis to the Theories of Electricity and Magnetism» - 72 страницы, заложивших теорию потенциала, формулы Грина, и сам термин «функция Грина». Тираж - 51 экземпляр. Только в 1850 году, через 9 лет после смерти Грина, Уильям Томсон (лорд Кельвин) обнаружил эту работу и сделал её знаменитой.
Без функций Грина нет ни современной квантовой теории поля (фейнмановские диаграммы), ни метода граничных элементов в инженерии, ни сетевых пропагаторов в физике конденсированного состояния.
Почему метод Гаусса-Зейделя для уравнения Лапласа сходится?
Ключевые идеи
- **∇²u = 0** - стационарные состояния. Принцип максимума: экстремум только на границе.
- **∇²u = f** - уравнение Пуассона. Решение единственно при заданных граничных условиях.
- **Функции Грина:** G(x,y) - отклик на точечный источник. В 3D: G = -1/(4πr). В 2D: G = ln(r)/2π.
- **Метод Гаусса-Зейделя:** u[i,j] = среднее четырёх соседей. Итерировать до сходимости.
Связанные темы
Уравнение Лапласа - прототип эллиптических PDE:
- Уравнение теплопроводности — При t→∞ теплопроводность приходит к решению Лапласа - стационарному распределению
- Классификация PDE — Уравнение Лапласа - эталон эллиптических PDE; для них итерационные солверы оптимальны
- Метод конечных элементов — FEM для уравнения Лапласа - слабая формулировка ∫∇u·∇v = 0
Вопросы для размышления
- Принцип среднего значения: гармоническая функция в точке = среднее по сфере. Как это связывает уравнение Лапласа с броуновским движением?
- В 2D функция Грина G = ln(r)/2π не убывает при r → ∞. Что это означает физически для двумерной электростатики?
- Метод Гаусса-Зейделя сходится медленно для больших сеток. Объясните интуитивно, почему многосеточный метод (multigrid) быстрее.