Фейнмановские лекции по физике. 7. Физика сплошных сред

Фейнман Ричард Филлипс

Из определения коэффициента вязкости вы видите, что h измеряется в ньютон·сек/м². Для воды при 20° С

h=10³ нъютон·сек/м².

Часто удобнее бывает пользоваться удельной вязкостью, которая равна h, деленной на плотность r. При этом величины удельных вязкостей воды и воздуха сравнимы:

Вода при 20°С h/r=10^{– 6}м²/сек

Воздух

при 20°С h/r=15·10^{– 6}м²/сек. , (41.12)

Обычно вязкость очень сильно зависит от температуры. Напри­мер, для воды непосредственно над точкой замерзания отно­шение h/r в 1,8 больше, чем при 20° С.

§ 2. Вязкий поток

Перейдем теперь к общей теории вязкого потока, по крайней мере настолько общей, насколько это и известно человеку. Вы уже понимаете, что компоненты сдвиговых напряжений сдвига пропорциональны пространственным производным от раз­личных компонент скорости, таких, как dv_x/dy или dv_y/дх. Однако в общем случае сжимаемой жидкости в напряжениях есть и другой член, который зависит от других производных скорости. Общее выражение имеет вид

где х_i — какая-либо из координат х, у или z; v_i — какая-либо з прямоугольных составляющих скорости. (Значок d_ij обозна­чает символ Кронекера, который равен единице при i=j и нулю при i№j.) Ко всем диагональным элементам S_ijтензора напряжений прибавляется дополнительный член h'С·v. Если жидкость несжимаема, то С·v=0 и дополнительного члена не появляется, так что он действительно имеет отношение к внутренним силам при сжатии. Для описания жидкости, точно так же как и для описания однородного упругого тела, требуются две постоянные. Коэффициент h представляет «обыч­ный» коэффициент вязкости, который мы уже учитывали. Он называется также первым коэффициентом вязкости, а новый коэффициент h' называется вторым коэффициентом вязкости.

Теперь нам предстоит найти вязкую силу f_вязк, действую­щую на единицу объема, после чего мы сможем подставить ее в уравнение (41.1) и получить уравнение движения реальной жидкости. Сила, действующая на маленький кубический объем жидкости, представляет собой равнодействующую всех сил, действующих на все шесть граней. Взяв их по две сразу, мы получим разность, которая зависит от производных напряжений, и, следовательно, от вторых производных скоростей. Это прият­ный результат, ибо он приведет нас опять к векторному урав­нению. Компонента вязкой силы, действующей на единицу объема в направлении оси х_i, равна

Обычно зависимость коэффициентов

вязкости от координат положения несущественна и ею можно пренебречь. Тогда вяз­кая сила на единицу объема содержит только вторые производ­ные скорости. Мы видели в гл. 39, что наиболее общей формой вторых производных в векторном уравнении будет сумма Лапласиана (С·С)v = С²v и градиента дивергенции (С (С·v)). Выражение (41.14) представляет как раз такую сумму с коэф­фициентами h и (h+h'). Мы получаем

В случае несжимаемой жидкости С·v=0 и вязкая сила в еди­нице объема будет просто равна hС²v. Это и все, чем обычно пользуются; однако если вам понадобится вычислить погло­щение звука в жидкости, то вам потребуется и второй член. Теперь мы можем закончить вывод уравнения движения реальной жидкости. Подставляя (41.15) в уравнение (41.1), получаем

Уравнение получилось, конечно, сложное, но ничего не поде­лаешь, такова природа.

Если мы введем W=СXv, как делали это раньше, то наше уравнение можно записать в виде

Мы снова предполагаем, что единственными объемными силами являются консервативные силы типа сил тяжести. Чтобы понять смысл нового члена, давайте рассмотрим случай несжимаемой жидкости. Если мы возьмем ротор уравнения (41.16), то полу­чим

Это напоминает (40.9) с той только разницей, что в правой части имеется еще одно слагаемое. Когда правая часть была равна нулю, то имелась теорема Гельмгольца о том, что вихри всегда движутся вместе с жидкостью. Теперь же в правой части появилось довольно сложное выражение, из которого, однако, не сразу же следуют физические выводы. Если бы мы пренебрегли членом СX(WXv), то получили бы диффузион­ное уравнение. Новый член означает, что вихри диффундируют в жидкости. При большом градиенте вихри расползаются в со­седние области жидкости.

Именно поэтому утолщаются кольца табачного дыма. С этим же связано красивое явление, возникающее при прохождении кольца «чистого» вихря (т. е. «бездымного» кольца, созданного с помощью описанной в предыдущей главе аппаратуры) через облако дыма. Когда оно выходит из облака, к нему «прилипает» некое количество дыма и мы видим полую оболочку из дыма. Какое-то количество завихренности W диффундирует в окру­жающий дым, продолжая свое движение вперед вместе с вихрем.

§ 3. Число Рейнольдса

Посмотрим теперь, как изменяется течение жидкости из-за нового члена с вязкостью. Рассмотрим несколько подробнее две задачи. Первая — обтекание жидкостью цилиндра; эту задачу мы пытались решить в предыдущей главе, используя теорию невязкой жидкости. Оказывается, что сегодня возможно найти решение вязких уравнений только для некоторых спе­циальных случаев. Так что кое-что из того, что я расскажу вам, основано на экспериментальных измерениях, считая, конечно, что экспериментальная модель удовлетворяла урав­нению (41.17).

Математически задача состоит в следующем: мы хотим найти решение для потока несжимаемой вязкой жидкости вблизи длинного цилиндра диаметром D. Поток должен опреде­ляться уравнением (41.17) и

W=СXv (41.18)

с условием, что скорость на больших расстояниях равна не­которой постоянной V (параллельной оси х), а на поверхности цилиндра равна нулю. Так что

v_я=v_у=v_z=0 (41.19)

при