§ 1. Вязкость
В предыдущей главе мы говорили о поведении воды, пренебрегая при этом эффектами вязкости. Теперь же мне хотелось бы обсудить, как вязкость влияет на течение жидкости. Рассмотрим реальное поведение жидкости. Я опишу качественно, как ведет себя жидкость в самых разных условиях, так чтобы вы получше прочувствовали эту науку. И хотя вы увидите сложные уравнения и услышите о трудных вещах, наша цель совсем не в том, чтобы изучить все тонкости. Цель этой главы скорее «общеобразовательная», просто я хочу дать вам некоторое понятие о том, как устроен мир. Однако здесь все же есть один пункт, который стоит того, чтобы его выучить: полезно знать простое определение вязкости. С него мы и начнем. Все же остальное предназначено для вашего удовольствия.
В предыдущей главе мы нашли, что законы движения жидкости содержатся в уравнении
(41.1)
В нашем приближении «сухой» воды мы отбрасывали последнее слагаемое, так что всеми эффектами вязкости мы пренебрегали. Кроме того, мы иногда делали еще дополнительное приближение, считая жидкость несжимаемой, и при этом получали дополнительное уравнение:
Это приближение часто оказывается вполне приличным, особенно когда скорость потока много меньше скорости звука. Но в реальных жидкостях мы почти никогда не можем пренебречь внутренним трением, называемым нами вязкостью; большинство интересных вещей в поведении жидкости так или иначе связано именно с этим свойством. Так, мы узнали, что циркуляция «сухой» воды никогда не изменяется: если ее не было вначале, то она никогда и не появится. Но в то же время мы повседневно сталкиваемся с циркуляцией в жидкости. Так что нашу теорию надо подправить.
Начнем с важного экспериментального факта. Когда мы занимались потоком «сухой» воды, обтекающей какой-то предмет или текущей мимо него, т. е. так называемым «потенциальным потоком», у нас не было причин запретить воде иметь составляющую скорости, тангенциальную к поверхности предмета; только нормальная компонента должна была быть равна нулю. Мы не принимали во внимание возможность возникновения сил сдвига между жидкостью и твердым телом. А вот оказывается, хотя это далеко и не очевидно, что во всех случаях, где это было проверено экспериментально, скорость жидкости на поверхности твердого тела в точности равна нулю. Вы замечали, конечно, что лопасти вентилятора собирают на себя тонкий слой пыли, и это несмотря на то, что они вращаются в воздухе. Тот же эффект можно наблюдать даже в больших аэродинамических трубах. Почему же пыль не сдувается воздухом? Несмотря на то что лопасти вентилятора быстро вращаются в воздухе, скорость воздуха относительно них, измеренная непосредственно на их поверхности, равна нулю, так что поток воздуха не возмущает даже мельчайших пылинок[60]. Мы должны модифицировать теорию так, чтобы она согласовалась с тем экспериментальным фактом, что во всех обычных жидкостях молекулы, находящиеся рядом с поверхностью, имеют нулевую скорость (относительно поверхности[61]).
Сначала мы характеризовали жидкость так, что если приложить к ней напряжение сдвига, то, сколь бы мало оно ни было, жидкость «поддается» и течет. В статическом случае никаких напряжений сдвига нет. Однако, когда равновесия еще нет, в момент, когда вы давите на жидкость, силы сдвига вполне могут быть. Вязкость как раз и описывает эти силы, возникающие в движущейся жидкости. Чтобы измерить силы сдвига в процессе движения жидкости, рассмотрим такой эксперимент. Предположим, что имеются две плоские твердые пластины, между которыми находится вода (фиг. 41.1), причем одна из пластин неподвижна, тогда как другая движется параллельно ей с малой скоростью v0.
Фиг. 41.1. Увлечение жидкости между двумя параллельными пластинками.
Если вы будете измерять силу, требуемую для поддержания движения верхней пластины, то найдете, что она пропорциональна площади пластины и отношению v0/d, где d — расстояние между пластинами. Таким образом, напряжение сдвига F/A пропорционально v0/d:
Коэффициент пропорциональности η называется коэффициентом вязкости.
Если перед нами более сложный случай, то мы всегда можем рассмотреть в воде небольшой плоский прямоугольный объем, грани которого параллельны потоку (фиг. 41.2).
Фиг. 41.2. Напряжения сдвига в вязкой жидкости.
Силы в этом объеме определяются выражением
(41.2)
Далее, ∂vx/∂y представляет скорость изменения деформаций сдвига, определенных нами в гл. 38, так что силы в жидкости пропорциональны скорости изменения деформаций сдвига.
В общем случае мы пишем
(41.3)
При равномерном вращении жидкости производная ∂vх/∂у равна ∂vy/∂x с обратным знаком, а Sxy будет равна нулю, как это и требуется, ибо в равномерно вращающейся жидкости напряжения отсутствуют. (Подобную же вещь мы проделывали в гл. 39 при определении еxy.) Разумеется, для Syz и Szx тоже есть соответствующие выражения.
В качестве примера применения этих идей рассмотрим движение жидкости между двумя коаксиальными цилиндрами. Пусть радиус внутреннего цилиндра равен а, его скорость будет vа, а радиус внешнего цилиндра пусть будет b, а скорость равна vb (фиг. 41.3).
Фиг. 41.3. Поток жидкости между двумя концентрическими цилиндрами, вращающимися с разными угловыми скоростями.
Возникает вопрос, каково распределение скоростей между цилиндрами? Чтобы ответить на него, начнем с получения формулы для вязкого сдвига в жидкости на расстоянии r от оси. Из симметрии задачи можно предположить, что поток всегда тангенциален и что его величина зависит только от r; v=v(r). Если мы понаблюдаем за соринкой в воде, расположенной на расстоянии r от оси, то ее координаты как функции времени будут
где ω=v/r. При этом х- и y-компоненты скорости равны
(41.4)
Из формулы (41.3) получаем
(41.5)
Для точек с у=0 имеем ∂ω/∂у=0, а х(∂ω/∂х) будет равно r(dω/dr). Так что в этих точках
(41.6)
(Разумно думать, что величина S должна зависеть от ∂ω/∂r, когда ω не изменяется с r, жидкость находится в состоянии равномерного вращения и напряжения в ней не возникают.)
Вычисленное нами напряжение представляет собой тангенциальный сдвиг, одинаковый повсюду вокруг цилиндра. Мы можем получить момент сил, действующий на цилиндрической поверхности радиусом r, путем умножения напряжения сдвига на плечо импульса r и площадь 2πrl:
(41.7)
Поскольку движение воды стационарно и угловое ускорение отсутствует, то полный момент, действующий на цилиндрическую поверхность воды между радиусами r и r+dr, должен быть нулем; иначе говоря, момент сил на расстоянии r должен уравновешиваться равным ему и противоположно направленным моментом сил на расстоянии r+dr, так что τ не должно зависеть от r. Другими словами, r3(dω/dr) равно некоторой постоянной, скажем А, и
(41.8)
Интегрируя, находим как ω изменяется с r:
(41.9)
Постоянные А и В должны определяться из условия, что ω=ωa в точке r=a, а ω=ωb в точке r=b. Тогда находим
(41.9)
Таким образом, ω как функция r нам известна, а стало быть, известно и v=ωr.
Если же нам нужно определить момент сил, то его можно получить из выражений (41.7) и (41.8):
или
(41.11)
Он пропорционален относительной угловой скорости двух цилиндров. Имеется стандартный прибор для измерения коэффициентов вязкости, который устроен следующим образом: один из цилиндров (скажем, внешний) посажен на ось, но удерживается в неподвижном состоянии пружинным динамометром, который измеряет действующий на него момент сил, а внутренний цилиндр вращается с постоянной угловой скоростью. Коэффициент вязкости определяется при этом из формулы (41.11).
Из определения коэффициента вязкости вы видите, что η измеряется в ньютон·сек/м2. Для воды при 20° С
Часто удобнее бывает пользоваться удельной вязкостью, которая равна η, деленной на плотность ρ. При этом величины удельных вязкостей воды и воздуха сравнимы:
(41.12)
Обычно вязкость очень сильно зависит от температуры. Например, для воды непосредственно над точкой замерзания отношение η/ρ в 1,8 больше, чем при 20° С.
§ 2. Вязкий поток
Перейдем теперь к общей теории вязкого потока, по крайней мере настолько общей, насколько это и известно человеку. Вы уже понимаете, что компоненты сдвиговых напряжений сдвига пропорциональны пространственным производным от различных компонент скорости, таких, как ∂vx/∂y или ∂vy/∂х. Однако в общем случае сжимаемой жидкости в напряжениях есть и другой член, который зависит от других производных скорости. Общее выражение имеет вид
(41.13)
где хi — какая-либо из координат х, у или z; vi — какая-либо з прямоугольных составляющих скорости. (Значок δij обозначает символ Кронекера, который равен единице при i=j и нулю при i≠j.) Ко всем диагональным элементам Sij тензора напряжений прибавляется дополнительный член η'∇·v. Если жидкость несжимаема, то ∇·v=0 и дополнительного члена не появляется, так что он действительно имеет отношение к внутренним силам при сжатии. Для описания жидкости, точно так же как и для описания однородного упругого тела, требуются две постоянные. Коэффициент η представляет «обычный» коэффициент вязкости, который мы уже учитывали. Он называется также первым коэффициентом вязкости, а новый коэффициент η' называется вторым коэффициентом вязкости.
Теперь нам предстоит найти вязкую силу fвязк, действующую на единицу объема, после чего мы сможем подставить ее в уравнение (41.1) и получить уравнение движения реальной жидкости. Сила, действующая на маленький кубический объем жидкости, представляет собой равнодействующую всех сил, действующих на все шесть граней. Взяв их по две сразу, мы получим разность, которая зависит от производных напряжений, и, следовательно, от вторых производных скоростей. Это приятный результат, ибо он приведет нас опять к векторному уравнению. Компонента вязкой силы, действующей на единицу объема в направлении оси хi, равна
(41.14)
Обычно зависимость коэффициентов вязкости от координат положения несущественна и ею можно пренебречь. Тогда вязкая сила на единицу объема содержит только вторые производные скорости. Мы видели в гл. 39, что наиболее общей формой вторых производных в векторном уравнении будет сумма Лапласиана (∇·∇)v=∇2v и градиента дивергенции (∇ (∇·v)). Выражение (41.14) представляет как раз такую сумму с коэффициентами η и (η+η'). Мы получаем
(41.15)
В случае несжимаемой жидкости ∇·v=0 и вязкая сила в единице объема будет просто равна η∇2v. Это и все, чем обычно пользуются; однако если вам понадобится вычислить поглощение звука в жидкости, то вам потребуется и второй член. Теперь мы можем закончить вывод уравнения движения реальной жидкости. Подставляя (41.15) в уравнение (41.1), получаем
Уравнение получилось, конечно, сложное, но ничего не поделаешь, такова природа.
Если мы введем Ω=∇×v, как делали это раньше, то наше уравнение можно записать в виде
(41.16)
Мы снова предполагаем, что единственными объемными силами являются консервативные силы типа сил тяжести. Чтобы понять смысл нового члена, давайте рассмотрим случай несжимаемой жидкости. Если мы возьмем ротор уравнения (41.16), то получим
(41.17)
Это напоминает (40.9) с той только разницей, что в правой части имеется еще одно слагаемое. Когда правая часть была равна нулю, то имелась теорема Гельмгольца о том, что вихри всегда движутся вместе с жидкостью. Теперь же в правой части появилось довольно сложное выражение, из которого, однако, не сразу же следуют физические выводы. Если бы мы пренебрегли членом ∇×(Ω×v), то получили бы диффузионное уравнение. Новый член означает, что вихри диффундируют в жидкости. При большом градиенте вихри расползаются в соседние области жидкости.
Именно поэтому утолщаются кольца табачного дыма. С этим же связано красивое явление, возникающее при прохождении кольца «чистого» вихря (т. е. «бездымного» кольца, созданного с помощью описанной в предыдущей главе аппаратуры) через облако дыма. Когда оно выходит из облака, к нему «прилипает» некое количество дыма и мы видим полую оболочку из дыма. Какое-то количество завихренности Ω диффундирует в окружающий дым, продолжая свое движение вперед вместе с вихрем.
§ 3. Число Рейнольдса
Посмотрим теперь, как изменяется течение жидкости из-за нового члена с вязкостью. Рассмотрим несколько подробнее две задачи. Первая — обтекание жидкостью цилиндра; эту задачу мы пытались решить в предыдущей главе, используя теорию невязкой жидкости. Оказывается, что сегодня возможно найти решение вязких уравнений только для некоторых специальных случаев. Так что кое-что из того, что я расскажу вам, основано на экспериментальных измерениях, считая, конечно, что экспериментальная модель удовлетворяла уравнению (41.17).
Математически задача состоит в следующем: мы хотим найти решение для потока несжимаемой вязкой жидкости вблизи длинного цилиндра диаметром D. Поток должен определяться уравнением (41.17) и
(41.18)
с условием, что скорость на больших расстояниях равна некоторой постоянной V (параллельной оси х), а на поверхности цилиндра равна нулю. Так что
(41.19)
при
Это полностью определяет математическую задачу.
Если вы вглядитесь в эти выражения, то увидите, что в задаче есть четыре различных параметра: η, ρ, D и V. Можно подумать, что нам придется иметь дело с целой серией решений для разных V, разных D и т. д. Вовсе нет. Все возможные различные решения соответствуют разным значениям одного параметра. Такова наиболее важная общая вещь, которую мы можем сказать о вязком потоке. А чтобы понять, почему это так, заметьте сначала, что вязкость и плотность появляются в виде отношения η/ρ, т. е. удельной вязкости. Это уменьшает число независимых параметров до трех. Предположим теперь, что все расстояния мы измеряем в единицах той единственной длины, которая появляется в задаче: диаметра цилиндра D, т. е. вместо х, у, z мы вводим новые переменные х', у', z', причем
При этом параметр D из (41.19) исчезает. Точно так же если будем измерять все скорости в единицах V, т. е. если мы положим v=v'V, то избавимся от V, а v' на больших расстояниях будет просто равно единице. Поскольку мы фиксировали наши единицы длины и скорости, то единицей времени теперь должно быть D/V, так что мы должны сделать подстановку:
(41.20)
В наших новых переменных производные в уравнении (41.18) тоже изменятся: так, ∂/∂х перейдет в (1/D)(∂/∂х') и т. д., так что уравнение (41.18) превратится в
(41.21)
А наше основное уравнение (41.17) перейдет в
Все постоянные при этом собираются в один множитель, который мы, следуя традиции, обозначим через 1/ℛ:
(41.22)
Если теперь мы просто запомним, что все наши уравнения должны выписываться для величин, измеряемых в новых единицах, то все штрихи можно опустить. Тогда уравнения для потока примут вид
(41.23)
и
с условиями,
для
(41.24)
и
для
Что все это значит? Если, например, мы решили задачу для потока с одной скоростью V1 и некоторого цилиндра диаметром D1, а затем интересуемся обтеканием цилиндра другого диаметра D2 другой жидкостью, то поток будет одним и тем же при такой скорости V2, которая отвечает тому же самому числу Рейнольдса, т. е. когда
(41.25)
В любых случаях, когда числа Рейнольдса одинаковы, поток при выборе надлежащего масштаба х', у', z' и t' будет «выглядеть» одинаково. Это очень важное утверждение, ибо оно означает, что мы можем определить поведение потока воздуха при обтекании крыла самолета, не строя самого самолета и не испытывая его. Вместо этого мы можем сделать модель и провести измерения, используя скорость, которая дает то же самое число Рейнольдса. Именно этот принцип позволяет нам применять результаты измерений над маленькой моделью самолета в аэродинамической трубе или результаты, полученные с моделью корабля, к настоящим объектам. Напомню, однако, что это можно делать только при условии, что сжимаемостью жидкости можно пренебречь. В противном случае войдет новая величина — скорость звука. При этом различные модели будут действительно соответствовать друг другу только тогда, когда отношение V к скорости звука тоже приблизительно одинаково. Отношение скорости V к скорости звука называется числом Маха. Таким образом, для скоростей, близких к скорости звука или больших, поток в двух задачах будет выглядеть одинаково, если и число Маха и число Рейнольдса в обеих ситуациях одинаковы.
§ 4. Обтекание кругового цилиндра
Вернемся теперь обратно к задаче об обтекании цилиндра медленным (почти несжимаемым) потоком. Я дам вам качественное описание потока реальной жидкости. О таком потоке нам необходимо знать множество вещей. Например, какая увлекающая сила действует на цилиндр? Сила, увлекающая цилиндр, показана на фиг. 41.4 как функция величины ℛ, которая пропорциональна скорости V, если все остальное фиксировано.
Фиг. 41.4. Коэффициент увлечения Сdкругового цилиндра как функция числа Рейнольдса.
Фактически на рисунке отложен коэффициент увлечения Сd — безразмерное число, равное отношению силы к 1/2ρV2Dl (d — диаметр, l —длина цилиндра, а ρ —плотность жидкости):
Коэффициент увлечения изменяется довольно сложным образом, как бы намекая нам на то, что в потоке происходит нечто интересное и сложное. Свойства потока полезно описывать для различных областей изменения числа Рейнольдса. Прежде всего, когда число Рейнольдса очень мало, поток вполне стационарен, скорость в любой точке потока постоянна и он плавно обтекает цилиндр. Однако распределение линий потока не похоже на их распределение в потенциальном потоке. Они описывают решение несколько другого уравнения. Когда скорость очень мала или, что эквивалентно, вязкость очень велика, так что вещество по своей консистенции напоминает мед, можно отбросить инерционные члены и описать поток уравнением
Это уравнение впервые было решено Стоксом. Он также решил задачу для сферы. Когда маленькая сфера движется при малых числах Рейнольдса, то к ней приложена сила, равная 6πηaV, где а — радиус сферы, а V — его скорость.
Это очень полезная формула: она говорит нам о скорости, с которой мельчайшие частички, которые приближенно можно считать шариками, движутся в жидкости под действием данной силы, как, например, в центрифуге, или при осаждении, или, наконец, в процессе диффузии. В области малых чисел Рейнольдса, т. е. при ℛ<1, линии v вокруг цилиндра имеют такой вид, как на фиг. 41.5.
Фиг. 41.5. Вязкий поток вблизи цилиндра (малая вязкость).
Если теперь мы увеличим скорость потока, так что число Рейнольдса станет несколько больше единицы, то увидим, что поток изменится.
Фиг. 41.6. Поток, обтекающий цилиндр, при различных числах Рейнольдса.
Как показано на фиг. 41.6, б, за сферой возникнут вихри. До сих пор неясно, существовали ли вихри и при малых числах Рейнольдса или же они возникли неожиданно при некотором определенном числе? Обычно считали, что циркуляция нарастает постепенно. Однако теперь думают, что скорее она проявляется неожиданно и возрастает с увеличением ℛ. Во всяком случае, поток в районе от ℛ=10 до ℛ=30 меняет свой характер. За цилиндром образуется пара вихрей.
Когда число Рейнольдса проходит через значения в районе 40, поток снова меняется. Характер движения претерпевает неожиданное и резкое изменение. Один из вихрей за цилиндром становится настолько длинным, что он отрывается и плывет вниз по течению вместе с жидкостью. При этом жидкость за цилиндром снова закручивается и возникает новый вихрь. Эти вихри поочередно отслаиваются то с одной, то с другой стороны, так что в какой-то момент поток выглядит приблизительно так, как показано на фиг. 41.6, в. Такой поток вихрей называется вихревой цепочкой Кармана. Она всегда появляется для чисел Рейнольдса ℛ>40. Фотография такого потока показана на фиг. 41.7.
Фиг. 41.7. Фотография цепочки вихрей в потоке за цилиндром.
Разница в режиме между двумя потоками, изображенными на фиг. 41.6, а, б или в, очень велика. На фиг. 41.6, а и б скорость постоянна, тогда как на фиг. 41.6, в скорость в любой точке изменяется со временем. Выше ℛ=40 стационарное решение отсутствует; граница перехода отмечена на фиг. 41.4 пунктирной линией. Для таких более высоких чисел поток изменяется со временем некоторым регулярным периодическим образом. Создаются вихри.
Можно представить себе физическую причину возникновения этих вихрей. Мы знаем, что на поверхности цилиндра скорость жидкости должна быть равна нулю, но при удалении от поверхности скорость быстро возрастает. Это большое местное изменение скорости жидкости и создает вихри. Когда скорость основного потока достаточно мала, у вихрей хватает времени, чтобы продиффундировать из тонкого слоя вблизи поверхности твердого тела, где они создаются, и «расплыться» на большую область. Эта физическая картина должна подготовить нас к следующему изменению природы потока, когда скорость основного потока или число ℛ увеличивается еще больше.
По мере возрастания скорости у вихря остается все меньше и меньше времени, чтобы «расплываться» на большую область жидкости. К тому моменту, когда число Рейнольдса достигнет нескольких тысяч, вихри начинают заполнять тонкую ленту (фиг. 41.6, г). В таком слое поток хаотичен и нерегулярен. Такая область называется пограничным слоем, и этот нерегулярный поток с увеличением ℛ пробивает себе путь все дальше и дальше вниз по течению. В области турбулентности скорости очень нерегулярны и «беспорядочны», вдобавок поток больше не двумерный — он крутится во всех трех измерениях. Кроме того, на турбулентное движение налагается еще регулярное переменное движение.
При дальнейшем увеличении числа Рейнольдса область турбулентности пробирается вперед, пока при потоке с ℛ, превышающим 105, не достигнет места, где линии тока огибают цилиндр. При этом поток будет похож на то, что показано на фиг. 41.6, д, и мы получаем так называемый «турбулентный след». Кроме того, происходят еще коренные изменения в силе увлечения — она, как видно из фиг. 41.4, сильно падает. При таких скоростях увлекающая сила с возрастанием скорости действительно уменьшается. По-видимому, здесь проявляется некоторое стремление к периодичности.
А что происходит при еще больших числах Рейнольдса? С дальнейшим увеличением скорости размер области турбулентности снова увеличивается и сила сопротивления возрастает. Последние эксперименты, которые дошли до области ℛ=107 или несколько больше, показывают, что в турбулентной области появляется новая периодичность, быть может, потому, что вся область колеблется вперед и назад в общем движении, а может быть, из-за нового сорта вихрей, которые появляются вместе с нерегулярным «шумовым» движением. Детали его полностью еще не ясны, и они до сих пор изучаются экспериментально.
§ 5. Предел пулевой вязкости
Мне бы хотелось подчеркнуть, что ни один из описанных нами потоков ни в каком отношении не похож на решение уравнения потенциального потока, о котором говорилось в предыдущей главе. На первый взгляд это очень удивительно. Ведь ℛ в конце концов пропорционально 1/η. Так что предел η→0 эквивалентен пределу ℛ→∞. И если мы перейдем к пределу больших ℛ в (41.23), то избавимся от правой части и получим как раз уравнения из предыдущей главы. Но все же трудно поверить, что сильно турбулентный поток с ℛ=107 хоть в какой-то степени приближается к гладкому потоку, вычисленному из уравнений «сухой» воды. Как может случиться, что при ℛ=∞ поток, описываемый уравнением (41.23), дает решение, полностью отличное от решения, полученного при η=0, с которого мы начали? Ответ очень интересен. Обратите внимание, что в правой части (41.23) стоит произведение 1/ ℛ на вторую производную. Это наиболее высокая степень производной в уравнении: слева только первые производные. Получается так, что, хотя коэффициент 1/ ℛ становится малым, Ω в пространстве вблизи поверхности претерпевает очень быстрые изменения. Эти резкие изменения компенсируют малость коэффициента, и произведение с увеличением R не стремится к нулю. Поэтому, хотя коэффициент при ∇2Ω стремится к нулю, решения не приближаются к предельному случаю.
Вас может удивить: «Что же такое мелкомасштабная турбулентность и как она может поддерживать сама себя? Как завихренность, которая создается где-то на краях цилиндра, приводит к такому шуму позади него?». Ответ снова очень интересен. Завихренность имеет тенденцию к самоусилению. Если мы на минуту забудем о диффузии завихренности, которая обусловливает потери, то законы потока говорят (как мы уже видели), что линии вихря переносятся вместе с жидкостью со скоростью v. Представьте себе некоторое количество линий Ω, которые возмущаются и скручиваются очень сложной картиной скоростей потока v. Прежде простые линии спутаются и сожмутся. Величина завихренности будет возрастать, равно как и ее нерегулярности (положительные и отрицательные), которые, вообще говоря, тоже будут увеличиваться. Таким образом, завихренность в трех измерениях по мере перемешивания жидкости будет возрастать.
Вы можете также спросить: «Когда же в конце концов справедлива теория потенциального потока?» Прежде всего она удовлетворительна вне турбулентной области, куда проникновение завихренности из-за диффузии незначительно. Изготовляя специальные обтекаемые тела, мы стараемся сделать область турбулентности как можно меньше. Поток, обтекающий крылья самолета, которые имеют специальную рассчитанную форму, — почти настоящий потенциальный поток.
§ 6. Поток Куеттэ
Можно показать, что сложный и изменчивый характер потока мимо цилиндра не исключение и что такое разнообразие возможностей получается и в общем случае. В § 1 мы нашли решение для вязкой жидкости между двумя цилиндрами и можем сравнить эти результаты с тем, что получается на самом деле. Если мы возьмем два концентрических цилиндра и заполним пространство между ними маслом с добавленной в него мелкой алюминиевой пудрой, то поток можно легко наблюдать. Если начнем медленно вращать внешний цилиндр, то ничего неожиданного не произойдет (фиг. 41.8, а).
Фиг. 41.8. Виды потока жидкости между двумя прозрачными вращающимися цилиндрами.
Можно медленно вращать и внутренний цилиндр, все равно ничего потрясающего не будет. А вот если мы начнем очень быстро вращать внутренний цилиндр — случится нечто удивительное. Жидкость разобьется на горизонтальные полосы (фиг. 41.8, б). Если с подобной же скоростью мы будем вращать внешний цилиндр, а внутренний оставим в покое, то никакого похожего эффекта не возникает. Как же получается, что не все равно, какой цилиндр вращать — внутренний или внешний. Ведь в конце концов вид потока, который мы нашли в § 1, зависел только от ωb-ωа. Ответ можно получить, взглянув на сечение цилиндра, изображенного на фиг. 41.9.
Фиг. 41.9. Вот почему поток разбивается на полосы.
Когда внутренние слои жидкости движутся быстрее, чем внешние, они стремятся двигаться наружу: центробежная сила становится больше удерживающего давления. Но весь слой целиком не может двигаться равномерно, так как на его пути стоят внешние слои. Поэтому они разбиваются на клетки и циркулируют, как показано на фиг. 41.9, б. Это напоминает конвекционные токи в комнате, где на уровне пола имеется слой теплого воздуха. Когда внутренний цилиндр находится в покое, а внешний цилиндр вращается с большой скоростью, центробежные силы создают градиент давления, который удерживает все в равновесии (фиг. 41.9, в), как теплый воздух, находящийся у потолка.
Теперь ускорим внутренний цилиндр. Сначала число полос увеличится. Затем неожиданно полосы станут волнистыми (см. фиг. 41.8,в), и волны эти начнут обтекать цилиндр. Скорость этих волн легко измерить. При больших скоростях вращения она приближается к 1/3 от скорости внутреннего цилиндра, а почему, никто не знает. Здесь есть над чем подумать. Простое число 1/3 и полное отсутствие объяснения! Вообще говоря, весь механизм образования волн тоже далеко не ясен, хотя мы имеем дело со стационарным ламинарным потоком.
Если теперь мы еще начнем вращать и внешний цилиндр, но в противоположную сторону, то картина потока начнет разбиваться. Волновые области начнут чередоваться со спокойными на вид областями, образуя спиральную картину (см. фиг. 41.8, г). Однако в этих «спокойных» областях, как можно заметить, поток на самом деле совсем не регулярен; он полностью турбулентен. Кроме того, в волновых областях начинает еще появляться нерегулярный турбулентный поток. Если цилиндры вращаются еще быстрее, то весь поток становится хаотическим турбулентным.
Этот простой эксперимент показал нам много интересных режимов потока, совершенно отличных один от другого и все же содержащихся в нашем простом уравнении при различных величинах одного-единственного параметра ℛ. С помощью наших вращающихся цилиндров мы можем наблюдать многие эффекты, проявляющиеся в потоке, проходящем мимо цилиндра: во-первых, это стационарный поток, во-вторых, целый набор потоков, которые изменяются со временем, но регулярным гладким образом, и, наконец, поток становится полностью нерегулярным. Те же самые эффекты каждый из вас видел в столбике табачного дыма, струящегося от сигареты, когда воздух спокоен. Сначала этот столбик гладкий, затем он как-то скручивается, поток дыма начинает разрушаться, и, наконец, все заканчивается беспорядочными клубами.
Основное, что вам следует вынести из всего сказанного, заключается в том, что в одном простом наборе уравнений (41.23) скрывается огромное разнообразие поведений. Все это решения одного и того же уравнения при различных значениях ℛ. У нас нет причин думать, что в этом уравнении мы потеряли какие-то слагаемые. Единственная трудность заключается в том, что нам сегодня не хватает математических знаний, чтобы проанализировать уравнение, за исключением очень малых чисел Рейнольдса, т. е. в случае очень вязкой жидкости. Написав уравнение, мы не отняли у потока жидкости ни его чарующей прелести, ни его таинственности, ни его поразительности.
Что ожидает нас в более сложных уравнениях, если даже в таком простом уравнении с одним-единственным параметром мы видим такое разнообразие возможностей! Вполне возможно, что основное уравнение, которое описывает завихрение туманностей, или образование вращений, или взрыв звезд и галактик, будет всего-навсего простым уравнением гидродинамики почти чистого водорода. Часто люди в каком-то неоправданном страхе перед физикой говорят, что невозможно написать уравнение жизни. А может быть, и можно. Очень возможно, что на самом деле мы уже располагаем достаточно хорошим приближением, когда пишем уравнение квантовой механики
Только что мы видели, как явления во всей их сложности легко и поразительно получаются из простых уравнений, которые описывают их. Не подозревая о возможностях простых уравнений, люди часто заключают, что для объяснения всей сложности мира требуется нечто данное от бога, а не просто уравнения.
Мы написали уравнения для течения воды. Но из нашего опыта у нас сложились какие-то понятия и приближения, пользуясь которыми, мы можем обсуждать разные решения — цепочку вихрей, турбулентный след, пограничный слой. Когда подобные уравнения встречаются нам в менее знакомой ситуации, где мы еще не можем экспериментировать, то мы пытаемся решать такие уравнения примитивным, извилистым и запутанным путем, стремясь определить, какие же качественные явления можно получить из него или какие новые качественные формы являются следствием этого уравнения. Наши уравнения для Солнца, например, представляющие его как водородный шар, описывают Солнце без солнечных пятен, без зернистой структуры его поверхности, без неровностей и короны. Тем не менее все это действительно находится в уравнениях, только у нас нет еще способа вытащить их оттуда.
Есть такие люди, которые будут очень расстроены, если на других планетах не будет найдено жизни. Я не принадлежу к их числу. И я никогда не смогу перестать удивляться и радоваться результатам межпланетных исследований, обнаруживающих бесконечное разнообразие и новизну явлений, порожденных одними и теми же самыми простыми принципами. Критерий науки — ее способность предсказывать. Могли бы вы предсказать бури, вулканы, океанские волны, зори и красочные закаты, если бы вы никогда не были на Земле?
Драгоценным сокровищем для нас будет все, что мы узнаем о происходящем на каждой из мертвых планет, каждого из десятка шаров, образовавшихся из того же самого облака пыли и подчиняющихся тем же самым законам физики, что и наша планета.
Грядущая великая эра пробуждения человеческого разума принесет с собой метод понимания качественного содержания уравнений. Сегодня еще мы не способны на это. Сегодня мы не можем увидеть в уравнениях потока воды такие вещи, как спиральное строение турбулентности, которую мы видим между вращающимися цилиндрами. Сегодня мы не можем сказать с уверенностью, содержит ли уравнение Шредингера и лягушек, и композиторов, и даже мораль или там ничего похожего и быть не может. Мы не можем сказать, требуется ли что-либо сверх уравнения, вроде каких-то богов, или нет. Поэтому каждый из нас может иметь на этот счет свое особое мнение.