Главная >> Фейнмановские лекции по физике >> Том 7 >> Глава 40. Течение «сухой» воды

Уравнение движения

Сначала обсудим движение жидкости с чисто абстрактной теоретической стороны, а затем рассмотрим некоторые частные примеры. Чтобы описать движение жидкости, мы должны задать в каждой точке ее некие свойства. Например, вода (будем называть жидкость просто «водой») в разных местах движется с различными скоростями. Следовательно, чтобы определить характер потока, мы должны в каждой точке и в любой момент времени задать три компоненты скорости. Если нам удастся найти уравнения, определяющие скорость, то мы будем знать, как в любой момент движется жидкость. Но скорость — не единственная характеристика жидкости, которая меняется от точки к точке. Только что мы изучали изменение давления от точки к точке. А есть еще и другие переменные. От точки к точке может меняться также плотность. Вдобавок жидкость может быть проводником и переносить электрический ток, плотность которого j изменяется от точки к точке как по величине, так и по направлению. От точки к точке может меняться температура, магнитное поле и т. д. Так что число полей, необходимых для полного описания ситуации, зависит от сложности задачи. Очень интересные явления возникают, когда доминирующую роль в определении поведения жидкости играют токи и магнетизм. Эта наука носит название магнитогидродинамика. В настоящее время ей уделяется очень большое внимание. Но мы не собираемся рассматривать эти весьма сложные случаи, ибо имеется немало менее сложных, но столь же интересных явлений, и даже этот более элементарный уровень будет достаточно труден.
 
Возьмем случай, когда нет ни магнитного поля, ни проводимости и нам, кроме того, не следует беспокоиться о температурах, ибо мы предположим, что температура в любой точке единственным образом определяется плотностью и давлением. Фактически мы уменьшим сложность нашей работы, допустив, что плотность постоянна, т. е. что жидкость существенно не сжижаема. Другими словами, мы предполагаем, что изменения давлений настолько малы, что производимыми ими изменениями плотности можно пренебречь. Если бы это было не так, то в дополнение к явлениям, рассмотренным здесь, необходимо было бы учитывать и другие явления, скажем распространение звуковых или ударных волн. Распространение звуковых и ударных волн мы уже в какой-то степени изучали, так что при нашем рассмотрении гидродинамики мы изолируемся от этих явлений, допустив, что приближенно плотность ρ постоянная. Легко определить, когда такое предположение о постоянстве ρ будет хорошим. Если скорость потока гораздо меньше скорости звуковой волны, то нам не нужно заботиться об изменениях плотности. Тот факт, что вода ускользает от нас при попытке понять ее, не связан с этим приближением постоянной плотности. Усложнения, которые все-таки позволили ей остаться непонятой, мы обсудим в следующей главе.
 
Общую теорию жидкостей мы должны начать с уравнения состояния жидкости, связывающего давление и плотность; в нашем приближении оно имеет очень простой вид:

Это и есть первое уравнение для наших переменных. Следующее соотношение выражает сохранение вещества. Когда вещество утекает из какой-то точки, то количество его в этой точке должно уменьшаться. Если скорость жидкости равна v, то масса, которая протекает за единичное время через единицу площади поверхности, равна нормальной к поверхности компоненте ρv. Подобное соотношение у нас получалось уже в теории упругости. Из знакомства с электричеством мы знаем также, что дивергенция такой величины определяется скоростью уменьшения  плотности. Также  и здесь уравнение

выражает сохранение массы жидкости: это гидродинамическое уравнение непрерывности. В нашем приближении, т. е. в приближении несжимаемой жидкости, плотность ρ постоянна и уравнение непрерывности превращается просто в

Дивергенция скорости жидкости v, как и магнитного поля В, равна нулю. (Гидродинамические уравнения очень часто оказываются аналогичными уравнениям электродинамики; вот почему мы сначала изучали электродинамику. Некоторые предпочитают другой путь, считая, что сначала следует изучать гидродинамику, чтобы потом было легче понять электричество. На самом же деле электродинамика гораздо проще, чем гидродинамика.)
 
Следующее уравнение мы получим из закона Ньютона; оно говорит нам, как происходит изменение скорости в  результате действия сил. Произведение массы элемента объема жидкости на ускорение должно быть равно силам, действующим на этот элемент. Выбирая в качестве элемента объема единичный объем и обозначая силу, действующую на единичный объем, через f, получаем

Плотность сил можно записать в виде суммы трех слагаемых. Одно из них, силу давления на единицу объема — (Vp), мы уже рассматривали. Но есть еще действующие на расстоянии «внешние» силы, подобные тяжести или электричеству. Если эти силы консервативные с потенциалом, отнесенным к единице массы, равным φ, то они приводят к плотности сил —ρ (Vφ). (Если же внешние силы не консервативные, то мы вынуждены писать внешнюю силу, приходящуюся на единицу объема, как f_внешн.) Кроме нее, на единицу объема действует еще одна «внутренняя» сила, которая возникает из-за того, что в текущей жидкости могут действовать сдвиговые силы. Они называются силами вязкости, и мы будем обозначать их через f_вязк. Тогда наше уравнение движения приобретает вид

В этой главе мы будем предполагать, что наша вода «жидкая» в том смысле, что ее вязкость несущественна, так что слагаемое f_вязк будет опускаться. Выбрасывая слагаемое с вязкостью, мы делаем приближение, которое описывает некое идеальное вещество, а не реальную воду. Об огромной разнице, возникающей в зависимости от того, оставляем ли мы слагаемое с вязкостью или нет, в свое время хорошо знал Джон фон Нейманн. Известно ему было и то, что во времена наибольшего расцвета гидродинамики, т. е. примерно до 1900 г., основные усилия были направлены на решение красивых математических задач в рамках именно этого приближения, которое ничего не имеет общего с реальными жидкостями. Поэтому теоретиков, которые занимались подобными веществами, он называл людьми, изучающими «сухую воду». Они отбрасывали важнейшее свойство жидкости. Именно потому, что в этой главе мы при наших вычислениях тоже этим свойством будем пренебрегать, я озаглавил ее «Течение «сухой» воды». А обсуждение настоящей, «мокрой» воды мы отложим до следующей главы.

Если мы отбросим f_вязк, то в уравнении (40.4) все нам известно, за исключением выражения для ускорения. Может показаться, что формула для ускорения частиц жидкости должна быть очень простой, ибо очевидно, что если v — скорость частицы в некотором месте жидкости, то ускорение ее будет просто равно ∂v/∂t. Но это совсем неверно, и по довольно хитрой причине. Производная ∂v/∂t выражает изменение скорости v (х, у, z, t) в фиксированной точке пространства. А нам нужно знать, как изменяется скорость данной капельки жидкости. Представьте, что мы пометили одну капельку воды цветной краской и можем наблюдать за ней. За маленький интервал времени Δt эта капелька продвинется в другое положение. Если капелька движется по некоторому пути, изображенному на фиг. 40.4, то за промежуток Δt она из точки Р₁ переместится в точку Р₂. Фактически в направлении оси х она передвинется на расстояние v_xΔt, в направлении оси у — на расстояние v_y Δt, а в направлении оси z — на расстояние v_zΔt. Мы видим, что если v (x, у, z, t) — скорость частицы в момент t, то скорость той же самой частицы в момент t+Δt представляет   величину   v (х + Δх, у + Δy,  z + Δz,   t + Δt), причем

Из определения частных производных [вспомните уравнения гл. 2, вып. 5] мы с точностью до членов первого порядка получаем

Ускорение же Δv/Δt будет равно

Считая V вектором, это можно записать символически:

Обратите внимание,  что, даже  когда  ∂v/∂t = 0, т. е.  когда скорость  в данной точке не изменяется,  ускорение все же останется. Примером может служить вода, текущая с постоянной скоростью по кругу: она ускоряется даже тогда, когда скорость в данной точке не изменяется. Причина, разумеется, состоит в том, что скорость данной капельки воды, которая первоначально находилась в одной точке, моментом позднее будет иметь другое направление — это центростремительное ускорение.
 
Остальная часть нашей теории — чисто математическая: нахождение решения уравнения движения, полученного подстановкой ускорения (40.5) в (40.4), т. е.

где слагаемое с вязкостью уже выброшено. Воспользовавшись известным тождеством из векторного анализа, это уравнение можно  переписать по-другому:

Если определить новое векторное поле Ω как ротор скорости v, т. е.

то векторное тождество можно записать так:

а наше уравнение движения (40.6) примет вид

Вы можете проверить эквивалентность уравнений (40.6) и (40.8), расписывая их по компонентам и сравнивая их, воспользовавшись при этом выражением (40.7).
 
Если Ω всюду равно нулю, то такой поток мы называем безвихревым (или потенциальным). В гл. 3, § 5 (вып. 5), мы уже определяли величину, называемую циркуляцией векторного поля. Циркуляция по любой замкнутой петле в жидкости равна криволинейному интегралу от скорости жидкости в данный момент времени вокруг этой петли:

Циркуляция на единицу площади для бесконечно малой петли по теореме Стокса будет тогда равна V х v. Таким образом, Ω представляет собой циркуляцию вокруг единичной площади (перпендикулярной направлению Ω). Кроме того, ясно, что если в любое место жидкости поместить маленькую соринку (именно соринку, а не бесконечно малую точку), то она будет вращаться с угловой скоростью Ω/2. Попытайтесь доказать это. Вы можете также попробовать доказать, что для ведра воды на вращающемся столике Ω равна удвоенной локальной угловой скорости воды.
 
Если нас интересует только поле скоростей, то из наших уравнений можно исключить давление. Взяв ротор обеих частей уравнения (40.8) и вспомнив, что ρ — величина постоянная, а ротор любого градиента равен нулю, а также использовав уравнение (40.3), находим

Это уравнение вместе с уравнениями

полностью описывают поле скоростей v. На языке математики — если в некоторый момент мы знаем Ω, то мы знаем ротор вектора скорости и, кроме того, знаем, что его дивергенция равна нулю, так что в этих физических условиях у нас есть все необходимое для определения скорости v повсюду. (Все это в точности напоминает нам знакомые условия в магнетизме, где V·B = 0 и V X В = j/ε₀c².) Таким образом, данная величина Ω определяет v точно так же, как j определяет В. Затем из известного значения v уравнение (40.9) даст нам скорость изменения Ω, откуда мы можем получить новую Ω в следующий момент. Используя снова уравнение (40.10), найдем новое значение v и т. д. Теперь вы видите, как в эти уравнения входит весь механизм, необходимый для вычисления потока. Заметьте, однако, что эта процедура дает только скорости, а всю информацию о давлении мы потеряли.
 
Отметим особое следствие нашего уравнения. Если в какой-то момент времени t повсеместно Ω=0, то ∂Ω/∂t тоже исчезает, так что Ω всюду останется равной нулю и в момент t +Δt. Отсюда следует, что поток все время остается безвихревым. Если вначале поток не вращался, то он так никогда и не начнет вращаться. При этом уравнения, которые мы должны решать, таковы:

Они в точности напоминают уравнения электростатики или магнитостатики в пустом пространстве. Позднее мы вернемся к ним и рассмотрим некоторые частные задачи.

СМОТРИТЕ ТАКЖЕ: