🌐 Матфизика

Метод характеристик

англ. Method of characteristics

Приём решения УЧП первого порядка: кривые-характеристики превращают уравнение в частных производных в обычное ОДУ вдоль кривой; данные едут по характеристикам, как багаж по рельсам.

Дым из трубы несёт ветром: каждая порция едет вправо со скоростью ветра, не меняя ни формы, ни концентрации. Это уравнение переноса $u_{t} + a u_{x} = 0$ , и метод характеристик решает его одним наблюдением: вдоль прямой $x - a t = const$ полная производная $\frac{d u}{d t} = u_{t} + a u_{x}$ обращается в ноль — решение на такой прямой постоянно. Кривые, вдоль которых УЧП вырождается в обыкновенное дифференциальное уравнение, называют характеристиками; отыскать их и прокатить по ним начальные данные — вся работа метода. Вместо борьбы с двумя частными производными решаем пару ОДУ: одно рисует траекторию, второе показывает, как меняется решение на ней.

Пример с числами: $u_{t} + 2 u_{x} = 0$ , $u (x, 0) = e^{- x^{2}}$ . Характеристики — прямые $x - 2 t = C$ , значения по ним не меняются, поэтому начальное значение из точки $ξ$ доезжает до $x = ξ + 2 t$ : ответ $u (x, t) = e^{- (x - 2 t)^{2}}$ — профиль уехал вправо на $2 t$ , не исказившись. Для общего уравнения $a (x, t, u) u_{x} + b (x, t, u) u_{t} = c (x, t, u)$ метод даёт характеристическую систему из трёх ОДУ — на этом языке переписывается и формула Даламбера для струны. О кривых-характеристиках уравнений второго порядка рассказано отдельно, в материале характеристики гиперболического уравнения; а когда характеристики изгибаются и начинают пересекаться, начинается территория квазилинейного уравнения. Полный разбор метода — в уроке метод характеристик, а вводная база по УЧП — в уроке уравнения в частных производных.

u_{t} + a u_{x} = 0 : \frac{d x}{d t} = a, \frac{d u}{d t} = u_{t} + a u_{x} = 0 \Rightarrow u (x, t) = u_{0} (x - a t)

уравнение переноса: вдоль характеристики x − at = const производная по времени полная нулевая, и начальный профиль u₀ просто сдвигается со скоростью a без изменения формы

a (x, t, u) u_{x} + b (x, t, u) u_{t} = c (x, t, u) : \frac{d x}{d s} = a, \frac{d t}{d s} = b, \frac{d u}{d s} = c

характеристическая система общего уравнения первого порядка: два уравнения рисуют кривую-траекторию, третье показывает, как меняется решение вдоль неё

Частые вопросы

Чем характеристики первого порядка связаны с характеристиками волнового уравнения?

Это одно понятие на двух уровнях. Для уравнения первого порядка характеристики — кривые, вдоль которых УЧП становится ОДУ, и метод ими решает задачу целиком. Для волнового уравнения второго порядка характеристики $x \pm a t = const$ — прямые, вдоль которых бегут возмущения; замена на характеристические координаты превращает волновое уравнение в $u_{η ξ} = 0$ и выводит формулу Даламбера. Механика одна: есть направления преимущественного переноса информации.

Когда метод характеристик перестаёт работать?

Когда характеристики пересекаются: тогда в одну точку приезжают два разных значения, и гладкого решения не существует. У линейного уравнения с постоянным коэффициентом $a$ все характеристики параллельны — пересечений нет. У квазилинейного скорость зависит от решения, характеристики изгибаются и могут столкнуться: возникает градиентная катастрофа и разрыв. После этого нужен аппарат слабых решений и законов сохранения — сам метод честно указывает, где гладкость кончилась.