Дифференциальные уравнения
Метод характеристик для PDE первого порядка
Ударные волны в аэродинамике сверхзвуковых самолётов, цунами в океане, пробки на дороге - всё это описывается уравнением Бюргерса и методом характеристик. Понимание формирования разрывов - не академическое упражнение: от этого зависит проектирование носового обтекателя и корректность численных схем.
- Аэродинамика: ударные волны при сверхзвуковом полёте, конус Маха
- Океанология: распространение цунами как мелководной волны
- Транспорт: модель ЛайтхиллаУизема-Ричардса для автомобильных пробок
- Газовая динамика: сжатие газа за поршнем, детонационные волны
- Численные методы: схемы годунова и Roe для гиперболических законов сохранения
Цели урока
- Применять метод характеристик для сведения PDE первого порядка к семейству ОДУ
- Определять время и место формирования разрыва из начальных данных для уравнения Бюргерса
- Использовать условие Ренкина-Гюгонио для нахождения скорости ударной волны
Предварительные знания
- Уравнения в частных производных первого порядка
- Обыкновенные дифференциальные уравнения: задача Коши
- Законы сохранения в интегральной форме
Метод характеристик
Для квазилинейного уравнения $u_t + c(u)u_x = 0$ характеристики - кривые $x = x(t)$ в плоскости $(x,t)$, вдоль которых $u$ постоянна: $\frac{dx}{dt} = c(u)$. Это превращает PDE в систему ОДУ: $\frac{du}{dt} = 0$ вдоль $\frac{dx}{dt} = c(u)$. Для линейного уравнения переноса $u_t + au_x = 0$: $u(x,t) = u_0(x-at)$ - бегущая волна.
Уравнение Бюргерса и формирование ударной волны
Уравнение Бюргерса: $u_t + uu_x = 0$. Характеристики: прямые $x = x_0 + u_0(x_0)t$ с наклоном $1/u_0(x_0)$. При убывающем начальном профиле $u_0'(x_0) < 0$ характеристики сходятся. Время формирования разрыва: $t_{break} = -1/\min_x u_0'(x)$. В этот момент производная $u_x \to -\infty$.
Характеристики: кривые, вдоль которых PDE = ODE
Unreal Engine 5 и Blender Cycles решают уравнение переноса излучения на 3D-сцене через метод характеристик: 1 луч = 1 характеристика. PDE в 3D сводится к ODE вдоль луча. Путь трассировка на RTX 4090 - 4096 лучей на пиксель при 60 fps - математически это 60*4096 решений ODE в секунду на пиксель.
Геометрический смысл: поверхность u(x,t) над плоскостью (x,t). PDE ограничивает наклон поверхности. Характеристика - кривая на плоскости (x,t), двигаясь вдоль которой скользишь по поверхности u без нарушения уравнения.
В чём ключевая идея метода характеристик для PDE первого порядка?
Метод характеристик находит кривые в пространстве (x,t), вдоль которых PDE первого порядка распадается на систему ОДУ. Это позволяет решать PDE аналитически, двигаясь вдоль этих кривых.
Линейное транспортное уравнение: бегущая волна
u_t + c*u_x = 0 - простейший случай. Физически: u(x,t) - концентрация загрязнения в реке, переносимого потоком со скоростью c. Начальный профиль u_0(x) переносится без изменения формы.
В ray tracing: каждый фотон движется вдоль характеристики du/ds = -sigma(x)*u, где sigma - коэффициент поглощения. Решение: u(s) = u_0 * exp(-integral sigma ds) - закон Бугера-Ламберта-Бера.
Чему равно u(5, 2) для задачи u_t + 3*u_x = 0, u(x,0) = cos(x)?
Транспортное уравнение u_t + c*u_x = 0 имеет решение u(x,t) = u_0(x - c*t): волна движется вправо со скоростью c=3. При x=5, t=2: аргумент 5-6=-1, cos(-1)=cos(1)~0.540.
Уравнение Бюргерса: самофокусировка волны
u_t + u*u_x = 0 - скорость переноса c = u зависит от самого решения. Более высокие значения u движутся быстрее. Если профиль убывает, «горб» волны нагоняет «впадину». Характеристики из разных точек имеют разные наклоны и могут пересечься. Пересечение = ударная волна.
Почему характеристики уравнения Бюргерса могут пересекаться, а линейного транспортного - нет?
В линейном транспортном уравнении скорость c фиксирована - все характеристики параллельны. В уравнении Бюргерса каждая характеристика несёт скорость u_0(x_0): если u_0 убывает, быстрые левые характеристики догоняют медленные правые и формируется ударная волна.
Ударные волны и условие Ренкина-Гюгонио
После t_break классическое решение перестаёт существовать - u многозначна. Физически: ударная волна (shock wave). В трафике - пробка. В газодинамике - сверхзвуковой скачок. Скорость разрыва определяется условием Ренкина-Гюгонио.
Модель трафика ПДЕ: d_rho/dt + d(rho*v(rho))/dx = 0. Скорость v(rho) = v_max*(1 - rho/rho_max). Пробка = ударная волна: плотность машин резко скачет. Скорость пробки по Ренкину-Гюгонио: s = (f(rho_L) - f(rho_R))/(rho_L - rho_R) - часто отрицательная (пробка движется навстречу потоку машин).
Для уравнения Бюргерса с u_L=2, u_R=0. Какова скорость ударной волны по Ренкину-Гюгонио?
Ренкин-Гюгонио для u_t + (u^2/2)_x = 0: s = (u_L^2/2 - u_R^2/2)/(u_L - u_R) = (u_L + u_R)/2 = 1. Условие Лакса: u_L > s > u_R - выполнено, шок физически корректен.
Условие Ренкина-Гюгонио
После формирования разрыва (ударной волны) скорость $s$ его распространения: $s = \frac{f(u_R) - f(u_L)}{u_R - u_L}$, где $f(u)$ - поток, $u_L$ и $u_R$ - значения слева и справа от разрыва. Для Бюргерса ($f = u^2/2$): $s = (u_L + u_R)/2$. Это следствие закона сохранения в интегральной форме.
Итоги
- Метод характеристик: вдоль $\frac{dx}{dt} = c(u)$ решение $u$ постоянно; PDE сводится к семейству ОДУ
- Уравнение Бюргерса: разрыв формируется в момент $t_{break} = -1/\min u_0'$, когда характеристики пересекаются
- Условие Ренкина-Гюгонио $s = \frac{[f]}{[u]}$ определяет скорость ударной волны из закона сохранения
Связь с другими темами
Метод характеристик обобщается на системы законов сохранения (уравнения Эйлера для газа), где характеристики - волны трёх типов. Энтропийные условия (Лакс, Олейник) выбирают физически правильное слабое решение среди нескольких кандидатов после формирования разрыва.
- Связанные темы — развивает
Вопросы для размышления
- Линейное уравнение переноса $u_t + au_x = 0$ имеет гладкое решение для любых начальных данных. Почему нелинейный случай принципиально отличается?
- Условие Ренкина-Гюгонио определяет скорость разрыва, но не единственность слабого решения. Что дополнительно нужно для единственности?
- Уравнение Бюргерса с вязкостью $u_t + uu_x = \varepsilon u_{xx}$: при малом $\varepsilon$ ударная волна 'размазывается'. Каков масштаб ширины переходной зоны?