Что такое уравнения в частных производных

Когда неизвестная функция зависит не от одной переменной, а от двух, трёх или четырёх сразу, обыкновенные дифференциальные уравнения заканчиваются — и начинается мир уравнений в частных производных.

Уравнение в частных производных (ДУЧП) — это уравнение, связывающее неизвестную функцию нескольких переменных и её частные производные по этим переменным.

Представьте металлический стержень, который вы нагрели в одном месте. Что описывает его состояние? Не одно число и не одна функция времени, а целое поле: для каждой точки x вдоль стержня и для каждого момента t существует своя температура. Математически это функция двух переменных u(x, t) — температура в точке x в момент t. Именно такие объекты — функции нескольких переменных, подчиняющиеся дифференциальным законам, — изучает теория ДУЧП. И почти вся физика сплошных сред, от тепла и звука до электромагнетизма и квантовой механики, записана на этом языке.

Частные производные: производная по одной переменной из многих

Если функция зависит от нескольких аргументов, то у неё есть не одна производная, а несколько — по каждой переменной отдельно. Такая производная называется частной. Для нашей температуры u(x, t) можно спросить два разных вопроса. Первый: как быстро меняется температура в фиксированной точке стержня с течением времени? Ответ даёт частная производная по времени, которую обозначают du/dt (правильнее писать ∂u/∂t, но мы будем пользоваться буквой d). Второй вопрос: как температура меняется вдоль стержня в фиксированный момент? Это частная производная по координате du/dx.

Ключевая идея частной производной проста: дифференцируем по одной переменной, считая все остальные временно замороженными константами. Когда мы берём du/dt, мы стоим в точке x и смотрим только на изменение во времени. Когда берём du/dx, мы фиксируем момент t и смотрим на профиль вдоль стержня — как фотографию.

Особую роль играет вторая частная производная по пространству, d^2u/dx^2. Она измеряет кривизну профиля — насколько он выгнут. Если профиль температуры в какой-то точке имеет форму выпуклого вверх горба, его вторая производная отрицательна; если форму впадины — положительна; если профиль линеен, вторая производная равна нулю. Эта величина окажется главной движущей силой во всех уравнениях этого раздела: природа стремится сглаживать кривизну.

Три классических типа ДУЧП

Линейные ДУЧП второго порядка делятся на три семейства, и это деление — не формальность, а отражение трёх принципиально разных типов поведения природы.

Стоит сразу подчеркнуть, зачем вообще нужна эта классификация. Тип уравнения (параболическое, гиперболическое или эллиптическое) определяется не словами, а структурой старших производных — и именно он диктует, какой численный метод вообще имеет смысл применять. Параболическое уравнение разумно решать шагами по времени с явной или неявной схемой, потому что у него есть выделенное направление времени и состояние эволюционирует вперёд. Гиперболическое требует осторожности с распространением волн и часто специальных схем, уважающих конечную скорость сигнала. Эллиптическое же вообще не шагают по времени — там сразу ищут установившееся поле, решая большую систему уравнений итерациями. Применить «тепловой» алгоритм к волновому уравнению или наоборот — значит получить либо взрыв, либо физически бессмысленный результат. Поэтому первое, что делает вычислитель, получив ДУЧП, — определяет его тип, и только потом выбирает инструмент.

Параболические — диффузия и сглаживание

Главный представитель — уравнение теплопроводности du/dt = α·d^2u/dx^2. Слева скорость изменения во времени, справа кривизна профиля. Смысл: где профиль выгнут впадиной (положительная кривизна), там температура растёт; где горбом — падает. Результат — неумолимое сглаживание: острые пики расплываются, перепады размываются, система движется к равновесию. Так распространяется тепло, диффундируют молекулы примеси в растворе, расплывается капля чернил.

Гиперболические — волны и распространение

Главный представитель — волновое уравнение d^2u/dt^2 = c^2·d^2u/dx^2. Здесь по времени стоит вторая производная (ускорение), и поведение совсем другое: возмущение не размывается, а бежит. Толкните струну — и горб поедет вдоль неё с конечной скоростью c, отразится от концов, вернётся. Так распространяются звук, свет, рябь на воде, сейсмические волны. Характерная черта — конечная скорость распространения сигнала.

Эллиптические — равновесие и стационарность

Главный представитель — уравнение Лапласа d^2u/dx^2 + d^2u/dy^2 = 0 (и его версия с источником — уравнение Пуассона). Здесь нет времени вообще: уравнение описывает не процесс, а установившееся состояние. Решение Лапласа — это профиль, у которого в каждой внутренней точке значение равно среднему по соседям; ни горбов, ни впадин внутри быть не может. Так выглядит распределение электрического потенциала в области без зарядов, форма провисшей мыльной плёнки, стационарная температура тела, которое уже перестало остывать.

Как работает под капотом

Почему ДУЧП радикально труднее обыкновенных уравнений? Дело в размерности задачи. Когда мы решали ОДУ, состояние системы в каждый момент описывалось конечным набором чисел — например, координатой и скоростью. Мы делали шаги по времени, и на каждом шаге обновляли эти несколько чисел. Конечномерная задача.

В ДУЧП состояние в каждый момент — это целая функция от пространства, то есть бесконечное число значений (по одному в каждой точке континуума). Формально это бесконечномерный объект. Чтобы загнать его в компьютер, мы вынуждены дискретизировать пространство: разбить стержень на конечное число узлов сетки и хранить температуру только в них. Тогда состояние снова становится конечным вектором — но теперь длинным, по числу узлов. А если задача двумерная (пластина, мембрана), узлов будет уже сетка N×N, и их число растёт как квадрат разрешения.

Чтобы почувствовать масштаб, прикинем числа на бытовых примерах каждого из трёх типов. Возьмём остывание стальной детали (параболическое уравнение теплопроводности): чтобы разрешить тонкий прогретый слой, по каждой из трёх пространственных осей берут хотя бы сотню узлов — это уже миллион неизвестных, и все их надо обновлять на каждом из тысяч шагов по времени. Прогрев двигателя, теплозащита спускаемого аппарата, остывание отливки — всё это та же история про миллионы клеток. Распространение звука в зале или сейсмической волны в породе (гиперболическое волновое уравнение) и расчёт антенны или электромагнитного поля живут по тем же законам размерности: трёхмерная область плюс время дают десятки и сотни миллионов неизвестных. А обтекание крыла потоком или распределение электрического потенциала в детали без зарядов (эллиптическое уравнение Лапласа) хоть и не зависят от времени, зато требуют разом решить гигантскую систему уравнений, связывающую все узлы сразу. Вот почему ДУЧП — это не «ОДУ, только посложнее», а отдельная вычислительная вселенная, ради которой строят суперкомпьютеры: даже скромная на вид физическая задача в трёх измерениях мгновенно упирается в миллионы и миллиарды чисел.

Дальше всё зависит от типа. В параболических и гиперболических уравнениях есть время, и мы шагаем по нему, на каждом шаге пересчитывая значения во всех узлах пространственной сетки сразу. То есть у нас две дискретизации: по пространству И по времени, и они взаимодействуют. В эллиптических уравнениях времени нет — там нужно сразу найти такое поле, которое всюду удовлетворяет уравнению равновесия; это приводит к большим системам уравнений, которые решают итерациями. Каждый из этих трёх миров требует своего численного подхода, и именно поэтому деление на типы так важно: оно подсказывает, какой алгоритм вообще имеет шанс работать.

Частые ошибки

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

Вторая ошибка — игнорировать краевые и начальные условия. У ОДУ хватало начального значения. У ДУЧП этого мало: нужно задать и начальный профиль (распределение в момент t=0), и условия на границах области на всё время (что происходит на концах стержня). Без них задача не определена.

Третья ошибка — считать, что раз мы аккуратно дискретизировали уравнение, ответ обязательно получится. Как мы увидим в последнем уроке, явная схема может взорваться, если шаг по времени слишком велик относительно шага по пространству. Дискретизация — это не безобидная замена, а тонкая операция со своими подводными камнями.

Итоги

• ДУЧП связывает функцию нескольких переменных и её частные производные; пример — температура u(x, t).
• Частная производная — это производная по одной переменной при замороженных остальных; вторая производная d^2u/dx^2 измеряет кривизну профиля.
• Три типа: параболические (диффузия, сглаживание), гиперболические (бегущие волны), эллиптические (равновесие).
• ДУЧП труднее ОДУ из-за бесконечномерности: приходится строить сетку по пространству и (для процессов) ещё и по времени.
• Обязательны и начальные, и краевые условия — иначе задача не определена.