Дифференциальные уравнения
Классификация PDE
Написать в коде неправильный метод для своего типа PDE - как использовать молоток для завинчивания шурупов. Формально работает, но катастрофически медленно или неустойчиво. Классификация PDE - это карта, указывающая нужный инструмент.
- **Симуляция физики в играх:** волновое уравнение (звук, взрывы) - upwind; теплопроводность (огонь) - неявная схема; статика (деформации) - эллиптический солвер. Все три типа в одном движке.
- **Числовой прогноз погоды:** уравнения Навье-Стокса - смесь гиперболических (адвекция) и эллиптических (давление) частей. Оператор расщепления разделяет их и решает разными методами.
- **Обратные задачи (ML):** восстановление источника тепла из наблюдений - некорректная обратная задача. Регуляризация Тихонова = добавление штрафного члена для стабилизации.
Предварительные знания
Классификация PDE второго порядка
Линейное PDE второго порядка: **A·u_xx + 2B·u_xy + C·u_yy + D·u_x + E·u_y + F·u = G**. Тип определяется знаком дискриминанта **Δ = B² - AC**: если Δ < 0 - **эллиптическое** (Лаплас, Пуассон), Δ = 0 - **параболическое** (теплопроводность), Δ > 0 - **гиперболическое** (волновое).
| Тип | Δ = B²-AC | Прототип | Физика |
|---|---|---|---|
| Эллиптическое | < 0 | ∇²u = 0 (Лаплас) | Стационарные поля |
| Параболическое | = 0 | ∂u/∂t = α∇²u | Диффузия, тепло |
| Гиперболическое | > 0 | ∂²u/∂t² = c²∇²u | Волны, распространение |
**Переменный тип:** уравнение Трикоми (u_xx + y·u_yy = 0) эллиптическое при y > 0 и гиперболическое при y < 0. Возникает в трансзвуковой аэродинамике при переходе через скорость звука. Моделирование требует разных методов в разных областях.
Жак Адамар и корректность задач (1902-1923)
В 1902 году Адамар в лекциях в Йельском университете сформулировал три аксиомы 'разумной' физической задачи: существование решения, единственность, непрерывная зависимость от данных. В 1923 в «Lectures on Cauchy's Problem» он показал, что задача Коши для уравнения Лапласа этим условиям не удовлетворяет: малое возмущение начальных данных порядка e^(-n) приводит к решению порядка e^n. Это убило надежды на обратное решение многих задач 'наивным' способом и породило целую отрасль - регуляризацию некорректных задач (Тихонов, 1963).
Современная регуляризация в ML (L2, dropout, weight decay) - наследник тихоновской регуляризации, выросшей из адамаровского анализа.
Уравнение u_xx + 4u_xy + 4u_yy = 0. Какого оно типа?
Характеристические кривые
**Характеристики** - это специальные кривые в плоскости (x,y), вдоль которых информация распространяется. Для гиперболического уравнения - 2 семейства характеристик. Для параболического - 1. Для эллиптического - характеристик нет (информация приходит «со всех сторон»).
Для волнового уравнения u_tt = c²u_xx характеристики: **x - ct = const** и **x + ct = const** - это именно «конусы влияния». Точка (x,t) получает информацию только из отрезка [x-ct, x+ct] начальных данных - это принцип конечной скорости распространения.
**Парадокс теплопроводности:** уравнение ∂u/∂t = α∂²u/∂x² даёт мгновенное распространение тепла - малое нагревание в одной точке сразу (при t=0+) влияет на всё пространство. Физически нереалистично, но математически удобно. Гиперболические тепловые уравнения это исправляют.
Для задачи Коши волнового уравнения с начальными данными на отрезке [-2, 2]. При t=1, c=1: какой отрезок влияет на точку x=3?
Корректность по Адамару
Задача PDE называется **корректной по Адамару**, если: 1. решение существует 2. единственно 3. непрерывно зависит от данных (устойчивость). Нарушение любого из этих условий делает задачу некорректной - «физически бессмысленной» с точки зрения моделирования.
| Тип PDE | Корректная постановка | Некорректная |
|---|---|---|
| Эллиптическое | Граничные условия везде | Задача Коши (неустойчива) |
| Параболическое | НУ при t=0, ГУ по x, прямое по t | Обратная задача (t→-∞) |
| Гиперболическое | НУ при t=0 (u и ∂u/∂t), прямое по t | ГУ на двух поверхностях |
Задача Коши для уравнения Лапласа (начальные данные на линии, а не на границе). Корректна ли она?
Выбор численного метода по типу PDE
Тип PDE определяет, какой численный метод эффективен. **Эллиптические** - нет «направления времени», решение глобальное → итерационные и прямые солверы (Гаусс-Зейдель, MG, FFT). **Параболические** - маршевые по времени, явные или неявные схемы. **Гиперболические** - маршевые с характеристиками, upwind-схемы.
| Тип | Метод | Особенности |
|---|---|---|
| Эллиптический | Итерационные (CG, MG), прямые (LU) | Нет явного направления маршировки |
| Параболический | Явный (FTCS), неявный (Кранк-Николсон) | Устойчивость: r≤0.5 для явных |
| Гиперболический | Upwind, Lax-Wendroff, WENO | CFL-условие: c·Δt/Δx ≤ 1 |
| Смешанный | Разные методы в разных областях | Трансзвуковая аэродинамика |
Задача на установившееся состояние (∂u/∂t = 0) электростатики в 3D. Какой метод предпочтителен?
Ключевые идеи
- **Дискриминант Δ = B² - AC:** < 0 эллиптическое, = 0 параболическое, > 0 гиперболическое.
- **Характеристики:** гиперболическое - 2 семейства (конус влияния), параболическое - 1, эллиптическое - нет.
- **Корректность по Адамару:** существование + единственность + устойчивость. Обратные задачи часто некорректны.
- **Численный метод:** эллиптическое → итерационные солверы; параболическое → маршировка по t; гиперболическое → upwind + CFL.
Связанные темы
Классификация объединяет все PDE уроки в единую систему:
- Уравнение теплопроводности — Параболический тип: одно семейство характеристик, маршировка по времени
- Волновое уравнение — Гиперболический тип: два семейства характеристик, конечная скорость
- Уравнение Лапласа — Эллиптический тип: нет характеристик, глобальные солверы
Вопросы для размышления
- Уравнение Трикоми меняет тип при y=0. Как подойти к численному решению в области, содержащей обе части?
- Обратная задача теплопроводности некорректна: маленький шум → огромные ошибки. Как ML-методы (нейросети) справляются с такими задачами?
- CFL-условие для гиперболических уравнений: c·Δt/Δx ≤ 1. Почему его нарушение приводит к нестабильности? Объясните через характеристики.