Укрощение бесконечности. История математики от первых чисел до теории хаоса - Йен Стюарт

В 1807 г. Жозеф Фурье представил свой доклад о путях распространения тепла Французской академии наук, но его не утвердили из-за серьезных недоработок. Чтобы подтолкнуть ученого к дальнейшим исследованиям, Академия даже учредила к 1812 г. Большую премию за изучение теплопроводности. Поскольку о награде стало известно заранее, Фурье уже в 1811 г. успел оформить свои идеи в виде доклада на соискание премии и выиграл ее. Но его труд всё равно жестоко критиковали за недостаток логической строгости, и Академия не разрешила его публиковать в виде научной статьи. Фурье, в ярости из-за такого отношения, написал труд «Аналитическая теория тепла», изданный в 1822 г. Туда почти полностью и без изменений вошел доклад 1811 г., но было и много новых материалов. Наконец в 1824 г. ученого оценили по заслугам: он был избран секретарем Академии и уже без помех опубликовал в виде научной статьи свой доклад от 1811 г.

Первым шагом Фурье был вывод ДУЧП для описания теплопроводности. Там имелось множество упрощений и допущений: тело должно быть однородным (с одинаковыми свойствами по всему объему) и изотропным (его свойства не зависят от направления) и т. д. В итоге он получил выражение, которое теперь известно как уравнение теплопроводности. Оно описывает распределение температуры в любой точке трехмерного тела и в зависимости от времени. Уравнение теплопроводности очень похоже с виду на уравнение Лапласа и на волновое уравнение, но с частной производной по времени в первой степени, а не во второй. Это небольшое отличие очень важно для математики ДУЧП в целом.

Были выведены такие же уравнения для одномерных и двумерных тел (стержень и плоскость), полученные удалением переменной z (для двумерного тела) и y (для одномерного). Фурье решил уравнение теплопроводности для стержня (чью длину мы принимаем равной π), на концах которого сохраняется неизменная температура, с условием, что в момент времени t = 0 (начальное состояние) температура в точке x стержня принимает вид:

(выражение, полученное с помощью предварительных вычислений), и сделал вывод, что температуру должно описывать схожее, но более сложное выражение, где каждый член умножается на соответствующую экспоненциальную функцию. Аналогия с гармониками в волновом уравнении поразительна. Но там каждая мода задана обычной синусоидой, колеблется бесконечно с одинаковой амплитудой, а здесь каждая синусоидальная мода распределения температуры убывает экспоненциально по времени, и более высокие моды убывают быстрее.

Представление с помощью ряда Фурье прямоугольной волны: вверху ее компоненты, синусоиды, внизу – их сумма

Причина такого различия в том, что в волновом уравнении энергия сохраняется, и поэтому колебание не затухает. А в уравнении теплопроводности тепло распространяется по всему стержню и теряется на его концах, потому что они охлаждаются.

Результатом работы Фурье стало то, что мы можем разложить начальное распределение температуры в ряд Фурье – сумму синусов и косинусов, похожую на приведенную выше формулу, а значит, способны немедленно описать, как тепло распространяется по телу со временем. Фурье считал очевидным, что такое выражение можно составить для любого начального распределения температуры, – и здесь-то начинались его неприятности. Мало кому из современников ученого было интересно, какое отношение теплопроводность имеет к волнам. Ее изучение казалось гораздо более сложным занятием.

Доводы Фурье в пользу возможности разложить функцию на синусы и косинусы были слишком сложными, запутанными и недостаточно строгими. Ему пришлось воспользоваться всеми разделами математики, чтобы в конце концов получить простые выражения для коэффициентов b1, b2, b3 и т. д. Обозначив начальное распределение температуры как f(x), он получил:

В 1777 г. Эйлер уже вывел эту формулу во время работы над волновым уравнением для звука и доказал ее с помощью мудрого наблюдения, заметив, что разные моды, sin mπx и sin nπx, являются ортогональными, т. е.

равен 0, если m и n – разные целые числа, не равные 0, т. е. на самом деле равен π/2, если m = n. Если предположить, что f(x) можно разложить в ряд Фурье, то, умножив обе стороны выражения на sin nx и проинтегрировав, мы избавимся от всех слагаемых, кроме одного, и в остатке получим формулу Фурье для bn.

Ни одно обсуждение ДУЧП в математической физике не будет полным без упоминания гидродинамики. И правда, эта область очень важна для практического применения, поскольку уравнения описывают, как вода обтекает подводные лодки или воздух – воздушные суда, и даже показывают сопротивление воздуха во время гонок «Формулы-1».

Эйлер сделал первые шаги в этой области в 1757 г., выведя ДУЧП для движения жидкости с нулевой вязкостью («липкостью»). Это уравнение остается в силе для некоторых жидкостей, но из-за излишней упрощенности не очень практично. Уравнения для вязких жидкостей вывел в 1821 г. Клод Навье, а потом их получил в 1829 г. Пуассон. Уравнения включают различные частные производные скорости движения жидкости. В 1845 г. Джордж Стокс вывел те же уравнения исходя из базовых физических принципов, и в итоге они получили название уравнения Навье – Стокса.