Дифференциальные уравнения движения идеальной жидкости

Размещено на http://www.allbest.ru/

Лекция 1. Дифференциальные уравнения движения идеальной жидкости

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

Рассмотрим движение элементарного объема жидкости в виде параллелепипеда со сторонами dx, dy, dz, который мысленно выделим в движущейся жидкости (рис. 1).

Рис. 1

В самом общем случае этот элементарный объем жидкости будет двигаться с каким-то ускорением, а его центр тяжести будет описывать какую-то криволинейную траекторию. Но согласно первому закону Ньютона, если материальная точка (представляемая элементарным объемом) совершает движение, отличное от прямолинейного и равномерного, то на нее действует некоторая сила dP, величина которой согласно второму закону Ньютона, равна произведению массы dm на векторную величину ее ускорения

d = dm .

Будем полагать в дальнейшем, что рассматривается движение идеальной жидкости. Силы, действующие на элементарный объем идеальной жидкости, будут складываться из массовых и поверхностных сил. Если плотность жидкости , то масса ее в элементарном объеме dW = dxdydz будет равна dm=dW= dxdydz. Пусть на каждую единицу массы жидкости в элементарном объеме действует массовая сила Ф. Тогда ко всей элементарной массе будет приложена массовая сила величиной Фdm=Фdxdydz. Обозначим проекцию единичной массовой силы Ф на оси х, у, z соответственно через X, У, Z.. Тогда проекции массовой силы на те же оси будут равны:

dxdydz; Ydxdydz и Zdxdydz.

Из поверхностных сил при движении идеальной жидкости на рассматриваемый объем жидкости будут действовать лишь силы гидродинамического давления, которые будут прикладываться к шести граням параллелепипеда. Найдем величину сил давления жидкости на грани, параллельные координатной плоскости zоу. Если гидродинамическое давление на левой грани будет равно р, то при помощи разложения в ряд Тейлора, отбрасывая члены высшего порядка малости, получим на правой грани гидродинамическое давление, равное р + dx.

Тогда давление жидкости на рассматриваемые грани площадью dydz будет равно соответственно pdydz и

dydz.

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



,

Полученные дифференциальные уравнения являются дифференциальными уравнениями движения идеальной жидкости. Эти уравнения были получены выдающимся математиком Л.Эйлером в 1755 году.

В случае движения капельной жидкости () в эти уравнения входят четыре неизвестных: P, нx, нy, нz.

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

Уравнение неразрывности (сплошности) потока

Рис. 2. К выводу дифференциального уравнения неразрывности потока.

Установим общую зависимость между скоростями в потоке жидкости, для которого соблюдается условие сплошности, или неразрывности движения, т. е. не образуется пустот, не заполненных жидкостью.

Выделим внутри потока элементарный параллелепипед объемом dV = dxdydz, ребра которого ориентированы параллельно осям координат (рис. *****).

Пусть составляющая скорости потока вдоль оси х в точках, лежащих на левой грани параллелепипеда площадью dS = dydz, равна Тогда через; эту грань в параллелепипед войдет вдоль оси х за единицу времени масса жидкости , а за промежуток времени -- масса жидкости

,

где -- плотность жидкости на левой грани параллелепипеда.

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

и .

Тогда через правую грань параллелепипеда за то же время выйдет масса жидкости

,

Приращение массы жидкости в параллелепипеде вдоль оси х:

= dx

Если составляющие скорости вдоль осей у и z равны и соответственно, то приращения массы в элементарном объеме вдоль этих осей по аналогии составят:

dx

Общее накопление массы жидкости в параллелепипеде за время равно сумме ее приращений вдоль всех осей координат:

dx

Вместе с тем изменение массы в полностью заполненном жидкостью объеме параллелепипеда возможно только вследствие изменения плотности жидкости в этом объеме. Поэтому

dx

Приравнивая оба выражения , сокращая на (--dx) и перенося - в левую часть уравнения, окончательно получим

Уравнение (5) представляет собой дифференциальное уравнение неразрывности потока для неустановившегося движения сжимаемой жидкости. Уравнение (5) может быть записано и в несколько иной форме. Проводя дифференцирование произведений , получим



,

Или

+

где -- субстанциональная производная плотности.

В установившемся потоке плотность не изменяется во времени, т. е. , и уравнение (5) принимает вид:

Для капельных жидкостей, которые практически несжимаемы, а также для газов в условиях изотермического потока при скоростях, значительно меньших скорости звука, = const и, следовательно

Уравнение (7) является дифференциальным уравнением неразрывности потока несжимаемой жидкости.

Сумма изменений скорости вдоль осей координат в левой части уравнения (7) называется дивергенцией вектора скорости и обозначается через div . Поэтому данное уравнение можно представить как

div = 0

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

Если бы площадь сечения трубопровода не изменялась, то для установившегося однонаправленного движения (в направлении оси х) интегрирование уравнения (6) дало бы зависимость

,

где -- средняя скорость жидкости.

Если же площадь сечения S трубопровода переменна, то, интегрируя также по площади, получим

Для трех различных сечений (1-1, 2-2 и 3-3) трубопровода, изображенного на рис. 2, имеем

Или

,

где М = -- массовый расход жидкости, кг/сек.

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

Рис. 2. К выводу уравнения постоянства расхода.

Для капельных жидкостей , и уравнение (8) принимает вид:

,

Следовательно

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

Согласно уравнению (8), массовый расход жидкости через начальное сечение трубопровода равен ее расходу через конечное сечение трубопровода. Таким образом, уравнение постоянства расхода является частным случаем закона сохранения массы и выражает материальный баланс потока.

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

Дифференциальные уравнения движения Навье--Стокса

При движении реальной (вязкой) жидкости в потоке жидкости помимо сил давления и тяжести действуют также силы трения.

Действие сил трения Т на выделенный в потоке вязкой жидкости элементарный параллелепипед (рис. 3) проявляется в возникновении на его поверхности касательных напряжений .

Рис. 3. К выводу уравнений Навье-Стокса.

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

В этих условиях касательные напряжения возникают лишь на поверхностях dF верхней и нижней граней элементарного параллелепипеда, причем dF = dxdy. Если касательное напряжение на нижней грани параллелепипеда равно , то на верхней оно составляет ( + . Производная выражает изменение касательного напряжения вдоль оси z в точках, лежащих на нижней грани параллелепипеда, а представляет собой изменение этого напряжения вдоль всей длины dz ребра параллелепипеда.

Указанные на рис. *** стрелками направления сил трения, приложенных к параллелепипеду на его нижней и верхней гранях, обусловлены, например, тем, что более медленные вышележащие слои жидкости затормаживают слой, в котором находится параллелепипед, а более быстрые нижележащие слои «разгоняют» его.

Тогда проекция равнодействующей сил трения на ось х

dxdy = -

Подставив в это выражение значение касательного напряжения по [ = - где -- вязкость жидкости], получим

,

В более общем случае трехмерного потока составляющая скорости х будет изменяться не только в направлении z, но и в направлениях всех трех осей координат. Тогда проекция равнодействующей сил трения на ось х примет вид

,

Сумму вторых производных по осям координат называют оператором Лапласа:

,

Следовательно, проекция равнодействующей сил трения на ось x; может быть представлена как

,

Соответственно проекции равнодействующей сил трения:

на ось у

,

на ось z

,

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

на ось х

,

на ось у

,

на ось z

,



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

,

Уравнения (***) представляют собой уравнения Навье -- Стокса, описывающие движение вязкой капельной жидкости.

При движении сжимаемой жидкости (газа) в ней дополнительно возникают вызванные трением силы сжатия и растяжения, поэтому уравнения Навье--Стокса принимают вид:

,

где частные производные выражают изменения скорости по осям х, у, z, связанные с действием сил сжатия и растяжения, причем

,

Левые части уравнений (10) выражают произведение массы единицы объема на проекцию ее ускорения, т. е. представляют собой проекции равнодействующей сил инерции, возникающих в движущейся жидкости.

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

Каждый член уравнений (10) имеет размерность соответствующей силы (тяжести, давления, трения или инерции), отнесенной к единице объема жидкости.

При движении идеальной жидкости, когда силы трения отсутствуют, уравнения движения Эйлера можно получить как частный случай уравнений Навье--Стокса.

Полное описание движения вязкой жидкости в его наиболее общей форме возможно путем решения уравнений Навье--Стокса совместно с уравнением неразрывности потока. Однако уравнения Навье--Стокса не могут быть решены в общем виде. Получены решения этой сложной системы уравнений только для некоторых частных случаев. Так, для установившегося ламинарного движения жидкости решение уравнений Навье-- Стокса позволяет вывести уравнение Пуазейля, полученное другим способом. ,

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

Лекция 2. Уравнение Д. Бернули для элементарной струйки идеальной жидкости

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

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

из внешних массовых сил действует лишь сила тяжести;

гидродинамическое давление является функцией координат и не зависит от времени;

жидкость является несжимаемой ()

Умножим уравнения Д. Эйлера соответственно на dx,, dy и dz.

,

,

,

Проекции ускорения массовой силы (в данном случае силы тяжести) примут следующие значения при выбранном; направлении осей координат:

X=0; Y=0; Z=-g.

После преобразования получим:

,

или:

,

Деля на g, получим:

,

Интегрируя это дифференциальное уравнение в полных дифференциалах, придем к следующему результату:

,

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

Для двух произвольных сечений элементарной струйки:

,

Это и есть уравнение Д. Бернули.

Геометрический и энергетический смысл уравнения

Д. Бернулли

Все члены, входящие в уравнение Д. Бернулли, имеют линейную размерность, поэтому их принято называть высотами. Соответственно общеприняты следующие названия для этих членов:

z - геометрическая или геодезическая высота;

- пьезометрическая высота или высота давления;

- динамический или скоростной напор;

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

Можно трактовать смысл отдельных членов уравнения

Бернулли иначе. Выше было показано, что сумма

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

z - есть удельная энергия положения;

- энергия давления;

- есть удельная кинетическая энергия.

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

Рис. 1

Полная удельная энергия (т. е. потенциальная + кинетическая) называется гидродинамическим напором и обозначается . Таким образом, уравнение Бернулли показывает, что при установившемся движении идеальной жидкости для данной струйки гидродинамический напор есть величина постоянная. На графике линия гидродинамического напора изображается горизонтальной линией.

Уравнение Д. Бернулли для элементарной струйки реальной жидкости. Пьезометрический и гидравлический уклоны.

При движении реальной жидкости между соседними струйками возникают силы трения, на преодоление которых затрачивается часть энергии жидкости. Поэтому удельная энергия жидкости в сечении элементарной струйки 2--2 будет менее удельной энергии жидкости в сечении 1--1 на некоторую величину , которую называют потерянной высотой или потерянной удельной энергией, затрачиваемой на преодоление гидравлических сопротивлений. Аналитически это положение запишется таким образом:

,

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

Легко изобразить уравнение Бернулли для рассматриваемого случая графически. Для этого следует, выбрав произвольную горизонтальную плоскость сравнения, отложить на ней в каждом сечении высоты ; ; ; и . Концы отрезков z, соединенные плавной кривой, покажут положение оси струйки. Соединяя концы отрезков плавной кривой, получим так называемую пьезометрическую линию. Отложив в каждом сечении вверх от пьезометрической линии отрезки, равные скоростным напорам , и соединив их концы плавной кривой, получим линию гидродинамического напора или, как ее часто называют, гидравлическую линию (рис. 2). Отрезки, равные расстояниям по вертикали от гидравлической линии до горизонтальной плоскости, проходящей над плоскостью сравнения на высоте, равной начальной удельной энергии, представляют собой потери энергии на гидравлические сопротивления на участке от начального до рассматриваемого сечения.



Рис. 2

Падение гидравлической линии на единицу длины элементарной струйки назовем гидравлическим уклоном I:

,

Гидравлический уклон (рис.3) есть всегда величина положительная, так как полная удельная энергия движущейся части жидкости постепенно уменьшается по мере ее продвижения вдоль элементарной струйки, затрачиваясь на преодоление сил трения, превращаясь в тепловую энергию и рассеиваясь.

Рис. 3

Понятие о плавно изменяющемся (медленно изменяющемся) движении жидкости

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

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

,

Уравнение Д. Бернули для потока реальной жидкости.

Условия применимости уравнения Д. Бернули.

Распространим уравнение Бернулли на установившийся поток реальной жидкости. Для этого выберем на слабо деформированном участке потока живое сечение, вблизи которого движение можно считать плавноизменяющимся.

Через это сечение каждой элементарной струйкой за время dt вносится энергия, которая в соответствии с вышеизложенным оказывается равной:

,

Вынося за скобки вес жидкости, прошедшей через поперечное сечение элементарной струйки за ;время dt, равный , перепишем это выражение в следующем виде:

,

Рис. 4

Найдем полную энергию, проносимую потоком жидкости через живое сечение 1-1. Для этого необходимо, очевидно, просуммировать полученное выражение по всем струйкам данного живого сечения. Тогда получим:

,

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

Запишем второй интеграл в следующем виде:

,

Этот интеграл представляет собой, как уже указывалось, кинетическую энергию, проносимую потоком через сечение 1-1 за время dt. Для его вычисления необходимо знать, каким образом распределяются скорости движения частиц жидкости по живому сечению. Если вычислить кинетическую энергию потока в предположении о постоянстве этих скоростей (другими словами, по средней скорости потока в данном живом, сечении ), то получим:

,

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

,,

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

,

Тогда:

,

После интегрирования и подстановки получим:

,

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

Полученное уравнение и есть уравнение Д. Бернули для установившегося потока реальной жидкости.

Лекция 3. Гидравлические сопротивления

Классификация гидравлических сопротивлений и потерь напора.

При движении реальной жидкости часть энергии потока затрачивается на преодоление гидравлических сопротивлений, которые подразделяются на два вида:

1) сопротивления по длине потока;

2) местные сопротивления.

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

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

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

1) потери напора по длине потока, которые вызываются гидравлическими сопротивлениями по длине потока (hf);

2) местные потери напора, которые вызываются местными гидравлическими сопротивлениями (hi). Общие потери напора:



hщ = .

Потери напора существенно зависят от режима движения жидкости.

Ламинарный и турбулентный режимы движения жидкости.

Различают два режима движения жидкости: ламинарный и турбулентный.

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

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

Слово «ламинарный» произошло от латинского слова lamina -- пластинка, полоска, слой; слово «турбулентный» произошло от латинского слова turbulentus - беспорядочный.

На существование в природе двух режимов движения жидкости впервые указал выдающийся русский ученый профессор Д. И. Менделеев в 1880 г. в работе «О сопротивлении жидкости и воздухоплавании».

Экспериментальное изучение режимов движения было проведено английским ученым О. Рейнольдсом в 1883 г.

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

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

Рис. 1

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

На основании экспериментального изучения режимов движения О.Рейнольдc дал критерий для установления того или иного режима движения.

Критерием для определения режима движения жидкости является так называемое число Рейнольдса, которое обозначается через Re и находится по формуле:



,

где V - средняя скорость движения потока;

l - характерный геометрический размер живого сечения потока;

- кинематический коэффициент вязкости.

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

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

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

Red(нк) = 20002320;

Red(вк) = 1000013000.

Лекция 4. Ламинарный режим движения жидкости в трубах

Распределение скоростей по живому сечению.

Рассмотрим равномерное ламинарное движение жидкости в горизонтальном напорном трубопроводе с радиусом r0 и диаметром d (рис. 1).

Рис. 1



Выделим двумя сечениями 1-1 и 2-2 отсек жидкости длиной l. В этом отсеке около оси трубы выделим цилиндрик жидкости с радиусом r и применим к нему основное уравнение равномерного движения:

.

В соответствии с законом внутреннего трения Ньютона:

.

Приравняв правые части двух последних равенств и разделив переменные, после интегрирования получим:

н = .

По этой формуле можно определить скорость в любой точке живого сечения потока при ламинарном движении жидкости в круглой трубе.

Из формулы (1) видно, что при ламинарном движении жидкости в круглой трубе распределение скоростей по живому сечению потока происходит по параболическому закону.

При r=0 н = нmax= (на оси трубы).

При r=r0 н=0 (у стенок трубы).

Эпюра скоростей показана на рис.1.

Распределение напряжения силы трения по живому сечению.

Напряжение силы трения:



.

При r=0 (на оси трубы).

При r=r0 max = 0 = 0 (на стенках трубы).

Разделим одно выражение на другое:

= ,

Откуда

.

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

Эпюра показана на рис.2.

Определение расхода и средней скорости.

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

н = .

Площадь поперечного сечения слоя dщ = 2rdr.

Количество жидкости, проходящее через данный слой в единицу времени, т.е. элементарный расход:



dQ = dщ н = 2rdr.

Полный расход:

Q = .

Интегрируя это выражение, окончательно получаем:

Q = = 4.

Средняя скорость:

V = ;

щ = .

Окончательно получаем:

V = 2.

Из сравнения V с нmax видно, что V = нmax, т.е. при ламинарном движении жидкости в круглой трубе средняя скорость равна половине максимальной.

Определение потерь напора на гидравлические сопротивления по длине.

Из формулы Пуазейля можно получить следующее выражение для потерь напора:



hпот = V.

Обозначим через

.

Тогда

hпот = ,

где коэффициент гидравлического трения.

Лекция 5. Турбулентный режим движения жидкости в трубах

I. Механизм и структура турбулентного потока

На рис. 1, а показан продольный разрез трубопровода, в котором жидкость движется турбулентно.

Экспериментальными исследованиями установлено, что при турбулентном режиме движения около стенок имеется тонкий ламинарный слой, толщина которого измеряется долями миллиметра. За ламинарным слоем, как показали исследования Г. А. Гуржиенко, располагается очень тонкий переходный слой, являющийся переходной зоной от ламинарного движения к турбулентному. Оба слоя вместе - ламинарный и переходный - составляют так называемый пограничный слой. За переходным слоем располагается ядро течения, в котором движение жидкости является уже турбулентным. Такова структура турбулентного потока.



Рис. 1

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

Исследования М. А. Великанова показали, что мгновенные местные скорости изменяются, колеблясь около некоторой постоянной величины, называемой осредненной местной скоростью.

Осредненная местная скорость равна:

,

Рис. 2



Осредненные местные скорости, вычисленные за достаточно продолжительный промежуток времени, являются почти постоянными и направленными параллельно оси потока. Следовательно, учет осредненных местных скоростей позволяет свести турбулентное движение к движению фиктивно-параллельному, к которому уже можно применять уравнение Бернулли.

Не следует смешивать понятий «средняя скорость потока» и «осредненная местная скорость»: средняя скорость потока, равная отношению расхода к площади живого сечения, постоянна во времени; осредненная же местная скорость в данной точке - это средняя скорость течения за определенный промежуток времени.

II. Толщина ламинарного слоя и скорость касательного напряжения

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

Напряжение сил трения на стенке:

,

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

,

где - толщина ламинарного слоя.

Решая это уравнение, после ряда преобразований окончательно получаем:

,

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

III. Напряжение сил трения в турбулентном потоке

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

,

Где - основное вязкостное напряжение, которое действует между отдельными слоями жидкости и находится по формуле Ньютона,

- дополнительное напряжение, возникающее вследствие перемешивания частиц жидкости.

Из теоретических работ Прандтля и Кармана следует:

,

Величину l Прандтль назвал длиной пути перемешивания и принял ее пропорциональной расстоянию п данной точки от стенки трубы, то есть:

l=Kn ,,

где К - коэффициент, постоянный для разных жидкостей и названный Карманом универсальной постоянной. По Карману, К=0,360,435.

Окончательно получаем:

,

,

Влияние каждого из членов, входящих в это выражение, в различных потоках и в разных точках живого сечения сказывается по-разному, а именно:

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

При ламинарном режиме вследствие отсутствия перемешивания жидкости величина l = 0 и формула (4) превращается в формулу (5). В этом случае пропорционально градиенту скорости, а следовательно, и средней скорости потока.

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

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

IV. Распределение скоростей по живому сечению потока

Для установления закона распределения скоростей по живому сечению потока при турбулентном режиме движения в круглой трубе воспользуемся уравнением (3). Подставим в него значение l = K n, разделим переменные и проинтегрируем:

,

Рис. 3

Окончательно получим:

,

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

Размещено на Allbest.ru