Гидродинамика вязкой жидкости

Введение

ГИДРОДИНАМИКА (от гидро... и динамика), раздел гидромеханики, изучает движение жидкостей и воздействие их на обтекаемые ими твердые тела. Теоретические методы гидродинамики основаны на решении точных или приближенных уравнений, описывающих физические явления в движущихся жидкости или газе. В экспериментальной гидродинамике возникающие задачи исследуются на моделях, обтекаемых жидкостью или газом, при этом должны соблюдаться условия подобия теории. Результаты гидродинамики используют при проектировании кораблей, самолетов, ракет и др.

 Гидродина­мика представляет собой раздел механики сплошных сред, в кото­ром изучается движение несжимаемых жидкостей и взаимодействие несжимаемых жидкостей с твердыми телами, — использует единый подход к изучению жидкостей и газов.

В механике с большой степенью точности жидкости и газы рассматриваются как сплошные, непрерывно распределенные в занятой ими части пространства. Плотность же газов от давления зависит существенно. Из опыта известно, что сжимаемостью жидкости и газа во многих задачах можно пренебречь и пользоваться единым поняти­ем несжимаемой жидкости — жидкости, плотность которой всюду одинакова и не изменяется со временем.

Если в покоящуюся жидкость поместить тонкую пластинку, то части жидкости, находящиеся по разные стороны от нее, будут действовать на каждый ее элемент ∆S с силами ∆F, которые независимо от того, как пластинка ориентирована, будут равны по модулю и направлены перпендикулярно площадке ∆S, так как наличие касательных сил привело бы частицы жидкости в движение.

Физическая величина, определяемая нормальной силой, действующей со стороны жидкости на единицу площади, называется давлением р жидкости:   

P = ∆F/∆S.

Единица давления — паскаль (Па): 1 Па равен давлению, создаваемому силой 1 Н, равномерно распределенной по нормальной к ней поверхности площадью 1 м² (1 Па=1 Н/м²).

Давление при равновесии жидкостей (газов) подчиняется закону Паскаля*: давле­ние в любом месте покоящейся жидкости одинаково по всем направлениям, при­чем давление одинаково передается по всему объему, занятому покоящейся жидкос­тью.

1. Коэффициент вязкости. Течение по трубе

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

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

Для выяснения закономерностей, которым подчиняются силы внутреннего трения, рассмотрим следующий опыт. В жидкость погружены две параллельные друг другу пластины, линейные размеры которых значительно превосходят расстояние между ними d. Нижняя пластина удерживается на месте, верхняя приводится в движение относительно нижней с некоторой скоро­стью v₀. Опыт дает, что для перемещения верхней пластины с постоянной скоростью v₀ необходимо действовать на нее с вполне определенной постоянной по величине силой F. Раз пластина не получает ускорения, значит, действие этой силы уравновешивается равной ей по величине противоположно направленной силой, кото­рая, очевидно, есть сила трения, действующая на пластину при ее движении в жидкости. Обозначим ее F_тр.

Варьируя скорость пластины v₀, площадь пластин S и расстоя­ние между ними d, можно получить, что

                                             (1)

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

При движении жидкости в круглой трубе скорость равна нулю у стенок трубы и максимальна на оси трубы. Полагая течение ла­минарным, найдем  закон   изменения скорости   с   расстоянием   r   от   оси трубы.

Выделим воображаемый цилинд­рический объем жидкости радиуса r  и длины l. При стацио­нарном течении в трубе постоянного сечения скорости всех частиц жидкос­ти остаются неизменными. Следовательно, сумма внешних сил, приложенных к любому объему жидкости, равна нулю. На осно­вания рассматриваемого цилиндрического объема действуют си­лы давления, сумма которых равна cила дей­ствует в направлении движения жидкости. Кроме того, на боковую поверхность цилиндра действует сила трения, равная                                        duldr на расстоянии r от оси трубы). Ус­ловие стационарности имеет вид

 (1)

Скорость убывает с расстоянием от оси трубы. Следовательно, duldr отрицательна и ldu/drl=—duldr. Учтя это, преобразуем соот­ношение следующим образом:

Разделив переменные, получим уравнение:

Интегрирование дает, что

Постоянную интегрирования нужно выбрать так, чтобы скорость обращалась в нуль на стенках трубы, т. е. при r=R (R — радиус трубы). Из этого условия

Подстановка значения С в (2) приводит к формуле

Значение скорости на оси трубы равно

(4)

С учетом этого формуле (3) можно придать вид

Таким образом, при ламинарном течении скорость изменяется с рас­стоянием  от  оси  трубы  по  параболическому  закону.

2. Формула Пуазейля.

Метод Пуазейля. Этот метод основан на ламинарном течении жидкости в тонком капилляре. Рассмотрим капилляр радиусом R и длиной /. В жидкости мысленно выделим цилиндрический слой радиусом r и толщиной dr. Сила внутреннего трения , действующая на боковую поверхность этого слоя,

 где dS — боковая поверхность цилиндрического слоя; знак минус означает, что при возрастании радиуса скорость уменьшается.

Для установившегося течения жидкости сила внутреннего трения, действующая на боковую поверхность цилиндра, уравновешивается силой давления, действующей на его основание:

После интегрирования, полагая, что у стенок имеет место прилипание жидкости, т. е. скорость на расстоянии R от оси равна нулю, получаем

Отсюда видно, что скорости частиц жидкости распределяются по параболическому закону, причем вершина параболы лежит на оси трубы. За время t из трубы вытечет жидкость, объем которой

откуда вязкость

3. Формула Стокса.

Формула Стокса. При малых Re, т. е. при небольших скоростях движения (и небольших /), сопротивление среды обусловлено практически только силами трения. Стокс установил, что сила сопротивления в этом случае пропорциональна коэффициенту динамической вязкости v движения тела относительно жидкости и характерному размеру тела I: (предполагается, что расстояние от тела до границ жидкости, например до стенок со­суда, значительно больше размеров тела). Коэффициент пропор­циональности зависит от формы тела. Для шара, если в качестве / взять радиус шара r, коэффициент пропорциональности оказывается равным 6я. Следовательно, сила сопротивления движению шарика в жидкостях при небольших скоростях в соответствии с формулой Стокса равна

                                    (1)

Метод Стокса. Этот метод определения вязкости основан на измерении скорости медленно движущихся в жидкости небольших тел сферической формы.

На шарик, падающий в жидкости вертикально вниз, действуют три силы: сила тяжести  (р — плотность шарика), сила Архимеда (р' — пло­тность жидкости) и сила сопротивления, эмпирически установленная Дж. Стоксом: , где r — радиус шарика, v — его скорость. При равномерном движении шарика

  или    

Откуда

Измерив скорость равномерного движения шарика, можно определить вязкость жид­кости (газа).

4. Закон подобия.

Геометрическое, кинематическое, динамическое подобие.

Этап изучения зависимости интересующей величины от системы выбранных определяющих факторов может выполняться двумя путями: аналитическим и экспериментальным. Первый путь применим лишь для ограниченного числа задач и при том обычно лишь для упрощенных моделей явлений.

Другой путь, экспериментальный, в принципе может учесть многие факторы, но он требует научно обоснованной постановки опытов, планирования эксперимента, ограничения его объема необходимым минимумом и систематизации результатов опытов. При этом должно быть обосновано моделирование явлений.

Эти задачи позволяет решать так называемая теория подобия, т. е. подобия потоков несжимаемой жидкости.

Гидродинамическое подобие складывается из трех составляющих: геометрического подобия, кинематического и динамического.

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

Отношение двух сходственных размеров подобных русел назовем линейным масштабом и обозначим эту величину через  .Эта величина одинакова для подобных русел I и II.

Кинематическое подобие означает пропорциональность местных скоростей в сходственных точках и равенство углов, характеризующих направление этих скоростей:

 Где – масштаб скоростей, одинаковый при кинематическом подобии.

Так как T – время, – масштаб времени).

Из кинематического подобия вытекает геометрическое подобие линий тока. Очевидно, что для кинематического подобия требуется геометрическое подобие русел.

Динамическое подобие – это пропорциональность сил, действующих на сходственные объемы в кинематических подобных потоках и равенство углов, характеризующих направление этих сил.

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

5. Турбулентность.

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

Режим течения данной жидкости в данной трубе изменяется примерно при определенной средней по сечению скорости течения v и обратно пропорционально диаметру d трубы, т. е.

Входящий в эту формулу безразмерный коэффициент пропорциональности k одинаков для всех жидкостей и газов, а также для любых диаметров труб. Это означает, что изменение режима течения происходит при определенном соотношении между скоростью, диаметром и вязкостью  v:

Полученное безразмерное число называется критическим числом Рейнольдса и обозначается

.

Этот результат согласуется с изложенной ниже теорией гидродинамического подобия, и вполне закономерно, что именно число Рейнольдса является критерием, определяющим режим течения в трубах.

Как показывают опыты, для труб круглого сечения.

Таким образом, критерий подобия Рейнольдса позволяет судить о режиме течения жидкости в трубе. При  течение является ламинарным, при  турбулентным. Точнее говоря, вполне развитое турбулентное течение в трубах устанавливается лишь при  имеет место переходная, критическая область.

Смена режима течения при достижении  обусловлена тем, что одно течение теряет устойчивость, а другое – приобретает. При  ламинарное течение является вполне устойчивым: всякого рода искусственная турбулязация потока и его возмущения (сотрясения трубы, введение в поток колеблющегося тела и пр.) погашаются влиянием вязкости и ламинарное течение восстанавливается. Турбулентное течение при этом неустойчиво. При  наоборот, турбулентное течение устойчиво, а ламинарное – неустойчиво.

6. Гидродинамическая неустойчивость

Итак, переход к турбулентности связан с неустойчивостью, а неустойчивость, в свою очередь, – с возникновением и развитием возмущений. Откуда же в реальной физической системе, какой является движущая жидкость, могут зародиться возмущения? Источников возмущений очень много. Прежде всего реальная установка (канал с движущейся жидкостью) находится на лабораторном столе, которому передаются колебания от стен и пола здания – результат сотрясения из-за проехавшей по соседству машины или, может быть, даже слабого сейсмического возмущения. Далее, вход жидкости в канал практически никогда не бывает идеально гладким, на входе в жидкость вносятся входные возмущения, они движутся вдоль жидкости вместе с ней и могут при благоприятных (неблагоприятных?) условиях нарастать. Стенки канала почти никогда не бывают лишены неровностей, шероховатостей. Обтекающий эти шероховатости поток непрерывно возмущается. Этот список можно было бы продолжать долго. Но есть источник возмущений, принципиально неустранимый. Это так называемые флуктуации. Когда мы говорим, например, что в данной точке потока плотность постоянна, это лишь означает, что она постоянна в среднем. Около этого среднего значения происходят малые, но макроскопические отклонения в ту или другую сторону. Они приводят к макроскопическим (малым) отклонениям (флуктуациям) давления, температуры и скорости. Флуктуации, таким образом, являются постоянно действующим источником возмущений, в принципе неустранимым.

Поставим теперь (мысленно) эксперимент по ламинарно-турбулентному переходу в трубе конечной длины. Вход в трубу постараемся сделать, насколько это возможно, гладким и постепенным, пытаясь устранить возмущения на входе. От шероховатости стенок также попытаемся отделаться благодаря тонкой шлифовке поверхности. Тот факт, что труба имеет конечную длину, также играет важную роль: представим себе, что в потоке жидкости возникло малое возмущение, которое, во-первых, сносится потоком вниз по течению и, во-вторых, в условиях неустойчивости нарастает. Для его роста требуется некоторое характерное время. Требуется время и для сноса возмущения потоком, оно просто равно (по порядку величины) длине трубы, поделенной на скорость потока. Если характерное время нарастания возмущения больше времени сноса, то оно не успеет вырасти на рабочем участке трубы и будет вынесено за его пределы. Если поставить опыт с учетом сделанных оговорок, то получится, что такие важные источники возмущений, как вход и шероховатость стенок, почти полностью устраняются, а те возмущения, которые все-таки возникнут, будут вытеснены потоком за пределы рабочего участка. Результаты такого опыта оказываются удивительными: удается существенно отодвинуть порог возбуждения турбулентности, критическое число Рейнольдса, таким образом, удается увеличить на 2-3 порядка, происходит "затягивание порога турбулентности".

Можно поставить также опыт с регулируемой шероховатостью стенок. Уменьшить шероховатость можно лишь до определенного предела, скажем до молекулярных размеров. Но можно ее искусственно увеличить, наклеивая на стенки, допустим, мелкие кристаллики контролируемых размеров. Таким образом, удается создать целую гамму трубок с оцениваемой наперед шероховатостью. Опыт говорит, что в этих случаях порог ламинарно-турбулентного перехода также изменяется в довольно широких пределах, причем критическое число Рейнольдса возрастает с уменьшением шероховатости.

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