Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Авдеев Е.Ф., Смирнова В.О. Конспект лекций по курсу Механика жидкости и газа

.pdf
Скачиваний:
177
Добавлен:
24.03.2021
Размер:
716.34 Кб
Скачать

М И Н И С Т Е Р С Т В О О Б Р А З О В А Н И Я И Н А У К И Р О С С И Й С К О Й Ф Е Д Е Р А Ц И И

ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ

«Национальный исследовательский ядерный университет «МИФИ»

Обнинский институт атомной энергетики –

филиал федерального государственного автономного образовательного учреждения высшего образования «Национальный исследовательский ядерный университет «МИФИ»

(ИАТЭ НИЯУ МИФИ)

Отделение Ядерной Физики и Технологий

Авдеев Е.Ф., Смирнова В.О.

Конспект

лекций по курсу «Механика жидкости и газа»

Обнинск, 2018г

СОДЕРЖАНИЕ

 

Предисловие ................................................................................................................................................

4

1. Введение. Аксиоматика механики жидкости и газа ........................................................................

5

2. Основные понятия кинематики и образы кинематики жидкости...................................................

5

2.1 Задание движения. Выражение для полного ускорения................................................................

5

2.2 Линии тока и трубки тока. Объемные и массовые расходы. Средняя и массовая скорости .....

6

2.3 Понятие циркуляции скорости.........................................................................................................

6

2.4 Вихревые линии и вихревые трубки. Распределение скоростей и давлений в сечении вихрей7

3. Математическое представление закона сохранения массы и его следствия .................................

8

3.1. Уравнение неразрывности...............................................................................................................

8

3.2. Потенциальные течения. Характеристическая функция течения................................................

9

3.3. Особенности кинематики турбулентных течений. Критическое число Рейнольдса .................

9

4. Математическое представление закона сохранение импульса. Интеграл Д. Бернулли и

 

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

9

4.1. Классификация сил, действующих в жидкости ..........................................................................

10

4.2. Уравнения динамики в напряжениях ...........................................................................................

11

4.3. Понятие давления...........................................................................................................................

11

4.4. уравнения движения вязкой жидкости (Навье-Стокса)..............................................................

12

4.5. уравнения движения идеальной жидкости Л. Эйлера ................................................................

12

4.6. Интеграл Д. Бернулли и область его действительности.............................................................

12

4.7. Частные случаи интеграла Бернулли ...........................................................................................

13

4.7.1. Изотермический процесс в газе .............................................................................................

13

4.7.2. Адиабатический процесс в газе .............................................................................................

13

4.8. Соотношения для покоящейся жидкости ....................................................................................

14

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

14

4.9.1. Плавно изменяющееся течение и его основные свойства ...................................................

14

4.9.2. Обобщение интеграла Бернулли на поток конечного сечения ...........................................

15

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

15

5.1. местные гидравлические сопротивления .....................................................................................

16

5.2. Сопротивление трения в каналах..................................................................................................

17

5.2.1. Ламинарное течение................................................................................................................

17

6. Тепловая формула интеграла Бернулли. Закон сохранения и превращения. Одномерный поток

 

газа. .............................................................................................................................................................

20

6.1 дифференциальное уравнение для закона сохранения и превращения энергии .......................

20

6.2. Тепловая формула интеграла Бернулли. Изоэнтропические формулы ....................................

21

6.3. Понятие скорости звука и числа M. Изоэнтропический формулы............................................

21

6.3.1 Постановка общей задачи........................................................................................................

21

6.3.2. Понятие скорости звука и числа M. ......................................................................................

21

6.3.3 Изоэнтропические формулы ...................................................................................................

22

6.4. Решения общей задачи...................................................................................................................

24

6.4.1. Истечение из конфузорного сопла ........................................................................................

25

6.4.2. Истечение через сопло Лаваля ...............................................................................................

26

6.4.3. Ударные волны и скачки уплотнения ...................................................................................

27

7. Особенности гидравлического удара во вскипающем теплоносителе. Кавитация. .......................

28

7.1. Прямые и непрямые гидроудары ..................................................................................................

29

7.2. Вскипающий теплоноситель .........................................................................................................

29

7.3. Кавитация........................................................................................................................................

29

8. Пограничный слой ................................................................................................................................

29

8.1. Физические представления о пограничном слое. Уравнения пограничного слоя Прандтля..

30

8.2. Интегральные соотношения Т.Кармана .......................................................................................

31

Предисловие

Изложение материала представлено по дедуктивному методу – от общего к частно-

му.

Чтобы не затенять физику процесса, промежуточные математические выкладки либо приведены частично, либо опущены.

Выдержана практическая направленность курса для специальности «Атомные электростанции: эксплуатация, проектирование и инжиниринг» или смежных направлений подготовки выпускников.

Неоценимую помощь в подготовке конспекта оказала студентка 5-го курса, Смирнова В.О., которую я заслуженно поставил соавтором.

Авдеев Е.Ф.

1. Введение. Аксиоматика механики жидкости и газа

Невозможно без установления основных положений МЖГ приступать к изучению следующих дисциплин: «Тепломасообмен», «Парогенераторы», «Паротурбинные установки» и т.д.

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

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

Введение модели идеальной жидкости (отсутствие трения) является вынужденной необходимостью, так как в общем случае для реальной среды будет 11 неизвестных (в числе которых 6 напряжений), а количество уравнений для их определения только 6. Введение модели идеальной жидкости сводит число напряжений к одному понятию давления, число неизвестных (6) приводится к числу уравнений. Если модель идеальной жидкости для внутренних задач течения в каналах практически ничего не дает (за исключением течения газа с малой вязкостью), то для внешних задач обтекания она дает практически 100% совпадение с опытом. Влияние сжимаемости за счет скорости начинается примерно с 2/5 скорости звука в среде.

Сжимаемость за счет скорости оказывает влияние, начиная примерно с 2/5 скорости звука в среде, до этих значений скоростей закономерности в капельной жидкости и газе одинаковы (если нет сильного влияния давления и температуры), поэтому под термином «жидкость» можно понимать и капельные жидкости и газ: объяснение многих реальных процессов можно получить только с привлечением особенностей акустики (скорости звука) в среде. Поэтому, приступающим к изучению дисциплины, кроме математического анализа, уравнений математической физики, термодинамики, необходимы первоначальные знания раздела физики «колебания и волны».

2. Основные понятия кинематики и образы кинематики жидкости

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

2.1 Задание движения. Выражение для полного ускорения

Вотличие от метода Лагранжа в МЖГ принят метод задания движения по Эйлеру

заданием поля скоростей в фиксированных точках пространства в любой момент вре-

мени t:

 

 

̅ ̅(

)

(2.1)

где x,y,z – координаты фиктивных точках, мимо которых движутся жидкие части-

цы.

Отсюда полное ускорение частицы, например вдоль оси X будет:

(2.2)

Где – называют локальным ускорением, а другие 3 слагаемые – конвективным

ускорением. Далее в уравнениях движения будет использоваться запись полного ускорения в виде:

̅

 

̅

(

 

)

̅

̅

(2.3)

 

 

 

 

 

 

 

 

2.2 Линии тока и трубки тока. Объемные и массовые расходы. Средняя и массовая скорости

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

.

(2.4)

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

(

 

)

(2.5)

 

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

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

Чтобы ввести величину, которая вдоль трубки тока всегда будет постоянной, объемный расход умножают на плотность и приходят к пониманию массового расхода G:

(

 

)

(2.6)

 

Его приравнивают к произведению массовой скорости

и площади сечения.

Массовая скорость играет важную роль в инженерных расчѐтов, так как остается постоянной при постоянстве сечения.

2.3 Понятие циркуляции скорости

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

Циркуляция скорости Г на участке произвольной линии тока АВ называют сумму произведений проекции скорости на направление касательной в данной точке на элемент дуги ds:

 

 

(2.7)

Она может быть представлена скалярным произведением вектора скорости и эле-

ментарного вектора касательной

̅, так как

.

 

∫ ̅ ̅

если

| ̅|. В случае потенциальных (безвихревых) течений, существует по-

тенциал скорости , и циркуляция скорости будет равна разности потенциалов в точке В и в точке A.

,то есть совпадает с потенциалом скорости. Наложение циркуляции скорости на симметричные профили искажает симметричность распределения скоростей

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

2.4 Вихревые линии и вихревые трубки. Распределение скоростей и давлений в сечении вихрей

Понятие вихревой линии вводится аналогично линиям тока, только теперь вдоль этой линии по касательной направлены угловые скорости ̅. Последнее означает, что поступательного движения двух вихревой линии нет, только вращение вокруг нее элементарных объемов. Конечно, может иметь место, когда ̅ параллельно ̅ , но тогда речь пойдет о вихревом движении. Аналогично поверхностям тока водится понятия вихревой поверхности. Тогда жидкость находится внутри вихревой поверхности, называемой вихревой трубкой. Интенсивностью вихревой трубки i (вихря) называют сумму произведений нормальных составляющих угловых скоростей на элементарную площадку сечения вихря , то есть поток вектора угловой скорости:

(2.8)

Так как поток вектора угловой скорости на боковые поверхности вихревой трубки равен нулю ( ), интенсивность вихревой трубки сохраняется постоянной вдоль неѐ.

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

-либо замкнутые торообразные;

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

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

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

 

Обозначив через радиус вихря, и вычислив циркуляцию скорости по контору ра-

диуса

(вихревая трубка одна) получим:

, откуда следует гиперболическое

уменьшение скорости от поверхности вихря Va , до нуля на бесконечности:

(2.9)

А внутри вихря будет линейное уменьшение скорости от Va до нуля в центре вихря. Для целостности изложения, позаимствуем связь между давлением и скоростью из

раздела 4.

Вне вихря движение безвихревое и для любых точек:

Определив постоянную из условия на бесконечности ( ), получим гиперболическое (в квадрате) уменьшение давления при приближении к вихрю:

(2.10)

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

(2.11)

Внутри вихря связь давления и окружной скорости иная:

(2.12)

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

(2.13)

В центре вихря (V=0) давление будет ниже давления окружающей среды на вели-

чину

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

Если вихря два, то, в зависимости от направления их закрутки, они будут двигаться как парный вихрь или параллельно, вызывая «пляску» расходов по ТВС.

3. Математическое представление закона сохранения массы и его следствия

3.1. Уравнение неразрывности

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

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

 

 

 

 

 

 

̅

(3.1)

 

 

 

 

или

 

 

 

 

 

 

 

̅

(3.1’)

 

 

 

 

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

жидкости, так как согласно (3.1), при

 

 

 

 

 

 

̅

(3.2)

Для плоских течений (3.2) допускает существование функции (

), скорость че-

рез которую выражается в виде:

 

 

 

 

 

 

 

 

 

 

(3.3)

 

 

 

 

 

 

 

 

 

Подставляя (3.3) в дифференциальное уравнение линии тока, после интегрирование

получаем уравнение семейства линий тока через функцию ( ):

 

.

 

(3.4)

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

 

 

(3.5)

Если

, нулевая линия тока, совпадающая со стенкой, тогда значение

бу-

дет объемным расходом.

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

3.2. Потенциальные течения. Характеристическая функция течения

Течения с потенциалом скорости вихревое). Из последнего следует, что осям:

называют потенциальным, когда

̅

(без-

̅

, значит составляющие скорости по

 

(3.6)

 

 

Связь потенциала скорости и функции тока (3.6) совпадает с условиями Коши-

Римана, являющимися необходимыми для существования функции

 

( )

(3.7)

зависящий от комплексного аргумента z=x+iy. Функцию (

) называют характери-

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

 

 

 

 

 

̅

(3.8)

 

 

 

 

 

где ̅ – сопряженная скорость.

 

 

Мнимая часть интеграла

̅

, действительная часть интегра-

ла

̅

 

 

 

 

3.3. Особенности кинематики турбулентных течений. Критическое число Рейнольд-

са

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

В течении столетия не удавалось определить критерий, по которому нужно устанавливать режим течения – ламинарное оно или турбулентное. Это не удалось сделать даже Д.И. Менделееву, который целенаправленно, по поручению правительства, занимался этой проблемой.

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

Для труб – средняя скорость, – кинематическая вязкость.

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

4. Математическое представление закона сохранение импульса. Интеграл нулли и уравнение Бернулли для потока конечного сечения

Теорему об изменении количества движения можно представить в виде:

( ̅) ̅

Д. Бер-

(4.1)

где ̅ – вектор внешних сил, а произведение вектора силы на элементарное время ( ̅ ) называют импульсу силы.

4.1. Классификация сил, действующих в жидкости

Влегко деформируемой жидкой среде (воде) силы классифицируются на массовые

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

 

̅

̅

(4.2)

 

 

 

 

 

 

где

– объем элементарной массы

 

 

. Согласно (4.2) вектор, с одной сто-

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

Для дальнейшего определим потенциалы массовых сил, в случае действия сил тя-

жести и центробежных сил инерции: ̅

̅

̅.

Исходя из понимания потенциала векторного поля:

̅

 

(4.3)

Проекции силы ̅ на оси будут связаны с потенциалом U зависимостями:

А полный дифференциал dU может быть записан через проекции силы:

(4.4)

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

Потенциал массовых сил после интегрирования (4.4):

 

 

.

(4.5)

При отсутствии вращения:

 

 

 

 

.

 

 

 

 

(4.6)

Поверхностные силы задаются вектором напряжения на площадках, ориентация ко-

торых определяется направлением нормали к ним:

 

̅

 

̅̅̅̅

.

(4.7)

 

 

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

̅

̅

̅

̅

(4.8)

Здесь индексы у векторов

напряжений

означают направления

к площадкам:

В проекциях на оси (4.8) дает три скалярных равенства:

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