Аэрогазодинамика
| Сдавался/использовался | 1998г. |
| Загрузить архив: | |
| Файл: vdv-0251.zip (53kb [zip], Скачиваний: 178) скачать |
Тема 1
Введение в аэрогазодинамику
1. Предмет, задачи и методы аэрогидромеханики. Задачи аэрогидро-
динамического расчёта.
2. Классификация видов и режимов движения жидкости.
3. Сравнение экспериментального, теоретического и вычислительного
подходов.
4. Вычислительная аэродинамика и этапы её развития.
1. Предмет, задачи и методы аэрогидромеханики
Однимиз основныхразделовсовременной физикиявляется
учение обаэрогидромеханике.Аэрогидромеханикаимеет дело с
жидкимиигазообразными средами.Жидкостиещё часто называют
капельными илинесжимаемымижидкостями, а вторые - газами или
сжимаемыми жидкостями.
Гидроаэромеханикаисследуетвопросы, связанныеспокоем
жидкости (гидростатика) и с её движением (гидродинамика).
Главноевнимание уделяется решению двух основных связанных
между собой задач: определения распределения скоростей и давлений
внутри жидкостииопределения силовоговзаимодействиямежду
жидкостью и окружающими её твёрдыми телами.
Теорияи экспериментявляются двумя основными подходами к
решению задач гидроаэродинамики.
Теоретическая гидроаэродинамикабазируетсяв основном на
невязкой(или так называемой идеальной) жидкости, внутри которой
отсутствует внутреннее трение.
Экспериментальная гидромеханика поставила своей целью
установить закономерности течения вязкой (реальной) жидкости.
Возникновение двух ветвей гидромеханики объяснялось
отсутствием достаточных представлений о механизметечения
жидкости и трудностями решенияуравненийдвижения вязкой
жидкости.
Всвязи с влиянием ... эффектов поток вязкой жидкости делят
надвеобласти: пограничный слой, где преобладают силы трения и
используются уравнения движения вязкой жидкости, и внешний поток,
ккоторомуможно применятьзакономерностидинамики невязкой
жидкости.
Наоснове решениязадачгидродинамики удаётсяполучить
теоретические зависимости, раскрывающие закономерности
сопротивлений, возникающийприобтекании тел (крыла и фюзеляжа
самолёта, лопаститурбины,кораблей различныхформи т.д.)
жидкостью.
Задачи аэродинамического расчёта
Процесспроектированияи конструированияЛА начинается с
проведения аэродинамического расчёта, в основу которого положены
две взаимозависимые задачи :
1) выбор аэродинамическойкомпоновки ЛА,
2) расчёт аэродинамических характеристик ЛА.
Привыборе аэродинамическойкомпоновки ЛА решаются задачи
отбора формы, размеров и взаимного расположения элементов ЛА.
В задачу расчёта АДХ ЛА входит:
1) расчёт распределения давления на поверхности ЛА,
2) расчёт составляющих аэродинамических сил и моментов,
3) расчёт аэродинамических характеристик органов управления,
4) расчёт температуры и тепловых потоков на поверхности ЛА.
Аэродинамическийрасчётобеспечивает исходныеданные для
проведения других работ в процессе проектирования ЛА.
1) расчёт тепловых режимов элементов конструкций,
2) расчёт траектории полёта,
3) расчёт динамических нагрузок,
4) расчёт управляемости и устойчивости.
2. Классификация видов движения жидкости
Проведём классификацию видов движения жидкости.
1. Классификация по признаку зависимости движения жидкости
от времени.
1.1. Установившееся (стационарное).
1.2. Неустановившееся (стационарное).
2. Классификация по признаку учёта сил трения, вязкости и
теплопроводности.
2.1. Идеальная невязкая жидкость.
2.2. Вязкая жидкость.
3. Классификация по виду движения жидкости (поступательное
или вращательного движение).
3.1. Безвихревое (потенциальное) (движение, когда вращение
отсутствует).
3.2. Вихревое движение.
4. Классификация по характеру изменения плотности в потоке.
4.1. Несжимаемая (жидкость),
4.2. Сжимаемая (газ),
5. Классификация по скорости и её отношению к скорости расши-
ряющихся возмущений (скорости звука).
5.1. Дозвуковое ()
а - скорость звука
5.2. Трансзвуковое ()
5.3. Сверхзвуковое ()
6. Классификация по режиму течения.
6.1. Ламинарный режим, ()
6.2. Турбулентный режим ()
7. Вид течения.
7.1. Свободное.
7.2. Вынужденное.
3. Сравнение экспериментального, теоретического и
вычислительного методов
┌───────────────┬──────────────────────┬────────────────────────┐
│ Метод │ Преимущества │ Недостатки │
├───────────────┼──────────────────────┼────────────────────────┤
│ │ │ │
│Эксперимент │1. Получение наиболее │1. Сложное оборудование │
│ │ близких к реальным │2. Проблемы моделирова- │
│ │ результатов │ ния │
│ │ │3. Обработка полученной │
│ │ │ информации. Кор. │
│ │ │ измер. значений │
│ │ │4. Сложность измерений │
│ │ │5. Стоимость │
├───────────────┼──────────────────────┼────────────────────────┤
│Теоретический │1. Получение информа- │1. Ограниченность про-│
│ │ ции в виде формул│ стейшими конфигураци-│
│ │ │ ями │
│ │ │2. Обычно применим толь-│
│ │ │ ко к линейным задачам│
├───────────────┼──────────────────────┼────────────────────────┤
│Численный │1. Нет ограничений, │1. Погрешность округле- │
│ │ связанных с нели-│ ния │
│ │ нейностью │2. Проблемы задания ГУ│
│ │2. Описание сложных │3. Стоимость ЭВМ │
│ │ физических процес- │ │
│ │ сов │ │
│ │3. Описание эволюции │ │
│ │ течения│ │
└───────────────┴──────────────────────┴────────────────────────┘
Основные этапы математичкеского моделирования
Рис.
Структурные элементы математического моделирования вместе со
связями показаны на рисунке.
Математическаяпостановказадачи базируется на физической
модели рассматриваемых течений, которая строится на основе имею-
щихся данных об объекте исследования.
Характеризующие математическуюмодель исходные уравнения и
граничные условия с помощью конечно-разностных методов преобразу-
ются в дискретную модель.
Врезультате реализациидискретной модели на одном из ...
программирования программу для ЭВМ. Решение тестовых задач и ана-
лизрезультатов позволяет убедиться в работоспособности разрабо-
танных алгоритмов и программ.
Решениеконкретных задачианализ полученных результатов
позволяет судитьобэффективности и применимости разработанных
алгоритмов.
Еслиобнаружится несоответствие расчётных и эксперименталь-
ныхданных-это значит, что физическая модель, математическая
модельилидискретная модель не адекватны изучаемому объекту. В
этомслучаепроводятся дополнительные исследования. Процесс ис-
следования продолжается до момента устранения.
4. Три этапа развития вычислительной аэрогидродинамики
┌─────────┬────────────────┬───────────────────┬─────────┬──────┐
│ Этап│ │ Год расчёта │ │Время │
│ (урав.) │ Результ.расчёт.├───────┬───────────┤ ЭВМ │рас-│
│ │ │профиль│реальн.ком.│ │чёта│
├─────────┼────────────────┼───────┼───────────┼─────────┼──────│
│ I │1. Распр. давл. │ │ │ IBM 360 │ │
│Ур.потен.│2. Индук. сопр. │ 1930│
1968 │ CDC 6600│
├─────────┼────────────────┼───────┼───────────┼─────────┼──────┤
│ II │1. Трансзвук │ │ │ 370 │ │
│Ур.Эйлера│2. Гиперзвук │ 1971│ 1976 │ 7600│ 5 │
├─────────┼────────────────┼───────┼───────────┼─────────┼──────┤
│ III │1. Отр. потока│ │ │ │ │
│Ур. Н.-С.│2. Турб. │ 1975│ 1985 │ CRAV │ 5 │
└─────────┴────────────────┴───────┴───────────┴─────────┴──────┘
Тема 2
Физические свойства жидкостей и газов
1. Различные состояния вещества. Твёрдые тела, жидкости и
газы. Силы, действующие на жидкости.
2. Основные свойства реальных жидкостей.
3. Поверхностное натяжение.
4. Уравнение состояния. Адиабата Тэйда.
1. Различные состояния вещества. Твёрдые тела, жидкости
и газы
В природе различают четыре агрегатных состояния вещества:
твёрдое, жидкое, газообразное и плазменное. Жидкость занимает про-
межуточное положение между твёрдыми телами и газами. Свойства жид-
костей при низкой температуре и высоком давлении близки к свойст-
вам твёрдых тел, а при высокой температуре и низком давлении - к
свойствам газов.
Жидкость, как всякое тело, имеет молекулярное строение, т.е.
состоит из молекул, объём пустот между атомами намного превосходит
объём самих молекул. Причём в жидкостях и твёрдых телах объём пус-
тот между молекулами меньше, а межмолекулярные силы больше, чем в
газах. Виду бесконечной малости молекул и пустот между ними по
сравнению с рассматриваемыми объёмами жидкости можно рассматривать
жидкости и газы в виде ... сплошной среды, придавая ей свойства
непрерывности.
Жидкость - это физическое тело, обладающее лёгкой подвижнос-
тью частиц, текучестью и способное изменять свою форму под воздей-
ствием внешних сил.
Жидкости разделяют на сжимаемые (газообразные) и несжимаемые
или весьма малосжимаемые (капельные).
Для облегчения изучения законов движения жидкости вводят по-
нятия идеальной и реальной жидкости.
Идеальные - невязкие жидкости, обладающие абсолютной подвиж-
ностью, т.е. отсутствием сил трения и касательных напряжений и аб-
солютной неизменностью а объёме под воздействием внешних сил.
Реальные - вязкие жидкости, обладающие сжимаемостью, сопро-
тивлением растягивающим и сдвигающим усилиям и достаточной подвиж-
ностью, т.е. наличием сил трения и касательных напряжений.
Реальные жидкости могут быть ньютоновскими и неньютоновскими
(бингамовскими). В ньютоновских жидкостях при движении одного слоя
жидкости относительно другого величина касательного напряжения
пропорциональна скорости сдвига. При относительном покое эти на-
пряжения равны нулю. Такая закономерность была установлена Ньюто-
ном в 1686 году, поэтому эти жидкости (вода, масло, бензин, керо-
син, глицерин и др.) называют ньютоновскими жидкостями.
Неньютоновские жидкости не обладают большой подвижностью и
отличаются от ньютоновских жидкостей наличием касательных напря-
жений (внутреннего трения) в состоянии покоя. Эта особенность бы-
ла подмечена Ф.Н.Шведовым (1889), а затем Бингемом (1916), поэтому
такие жидкости (битум, гидросмеси, глинистый раствор, коллоиды,
нефтепродукты при температуре, близкой к температуре застывания)
получили и другое название - бингемовские.
Силы, действующие в жидкости, принято делить на внутренние и
внешние.
Внутренние силы представляют собой силы взаимодействия частиц
жидкости, внешние силы делятся на силы поверхностные и объёмные.
Поверхностные силы (сжатие, давление, растяжение, силы тре-
ния) приложены к поверхностям, ограничивающим объём жидкости.
Объёмные силы (например, сила тяжести, сила инерции, электро-
магнитная сила) распределяются по всему объёму жидкости.
С термодинамической точки зрения состояние жидкости или пара
характеризуется тремя параметрами: давлением ..., плотностью ...
и температурой Т, связанными между собой уравнением состояния.
Исходной единицей давления в Международной системе единиц СИ
является паскаль:
...
На практике используют более крупные единицы - гектапаскаль
(1гПа = ... Па), килопаскаль (1кПа = ... Па), бар (1бар = ... Па)
и метапаскаль (1МПа = ... Па).
В технической литературе часто встречается другая единица из-
мерения давления - техническая атмосфера (ат).
...
Плотность выражается в единицах массы, приходящихся на едини-
цу объёма.
Исходной единицей массы в СИ служит
Размерность плотности
...
Основные гипотезы и понятия сплошной среды
Классическая гидромеханика основана на трёх утверждениях:
1) справедлива классическая механика - механика Ньютона
2) справедлива классическая термодинамика
3) справедлива гипотеза сплошности.
Первое утверждение предполагает, что изучаются движения со
скоростями, малыми по сравнению со скоростями света и рассматри-
ваются объекты, размеры которых существенно превосходят размеры
микромира.
Второе утверждение предполагает, что в окрестности каждой
точки жидкость находится в состоянии термодинамического равнове-
сия, вследствие чего можно пользоваться термодинамическими зако-
нами.
Третье утверждение предполагает замену реальной жидкости с
её дискретным молекулярным строением моделью сплошного распределе-
ния вещества по рассматриваемому объёму. Возможность такой замены
и носит название гипотезы сплошности.
Плотность жидкости в данной точке определяется как предел:
...
В системе СИ единица плотности...
В технических приложениях часто используется такая единица
СИ - вес единицы объёма или удельный вес:
...
Объёмные и поверхностные силы
Поверхностные силы (сжатие, давление, растяжение, силы тре-
ния) приложены к поверхностям, ограничивающим объём жидкости.
Объёмные силы (например, силы тяжести, сила инерции, электро-
магнитная сила) распределяются по всему объёму жидкости.
Пусть ... - главный вектор объёмных сил, действующих в
объёме ... . Тогда вводится понятие плотности распределения
объёмных сил в виде предела
...
Рассмотрим поверхностные силы. Пусть... - главный вектор
силы, приложенной с одной стороны, к площадке... . Индекс "..."
означает не проекцию силы, а указание на то, что сила действует на
площадке ... , произвольно ориентированной в пространстве. Введём
в рассмотрение вместо силы напряжение
...
Рассмотрим тетраэдр, три грани которого параллельны коорди-
натным плоскостям, а четвёртая ориентирована произвольным образом.
...
...
Обозначим площади граней ...
Ориентация площади определяется единичной нормалью ... с направля-
ющими косинусами ... . Тогда справедливы соотношения
...
Пусть высота тетраэдра равна ... . Тогда его объём равен
... . Воспользуемся вторым законом Ньютона и со-
ставим уравнение движения тетраэдра:
...
...
где ... - ускорение центра масс тетраэдра.
Переходя к пределу (устремляя ...), получим
...
Получим формулу Коши, утверждающую, что напряжения на гранях
образуют систему взаимно уравновешенных напряжений.
Проектируя векторное уравнение на оси координат, получим три
скалярных уравнения:
...
...
...
Напряжённое состояние в произвольной точке сплошной среды ха-
рактеризуется девятью компонентами, образующими тензор второго
ранга или диаду:
...
Тензор напряжений в произвольной точке пространства обладает
свойством симметрии (теорема Коши о взаимности касательных напряже-
ний)
...
Он содержит лишь шесть независимых компонент.
Рассмотрим равенство Коши для случая отсутствия касательных
напряжений, т.е. полагая ...= 0. Поскольку вязкость по гипотезе
Ньютона проявляется только при наличии неоднородного поля скорос-
тей, сделанное предположение будет соответствовать либо покою жид-
кости, либо её движению как твёрдого тела6 либо равенству нулю
вязкости (... = 0).
Итак
...
С другой стороны,
...
Сравнивая равенства, находим
...
Введём понятие давления Р согласно равенствам
...
Таким образом, в случае отсутствия касательных напряжений
давление в точке является скалярной величиной, т.е. оно не зави-
сит от ориентации площадки, проходящей через рассматриваемую точ-
ку. Знак минус означает, что давление рассматривается как сжимаю-
щее напряжение.
Температура жидкости выражается в единицах градусов абсолют-
ной шкалы
...
2. Основные свойства реальных жидкостей
Сжимаемость. При сжатии реальные жидкости незначительно умень-
шаются в объёме. Свойство жидкостей изменять объём при изменении
давления характеризуется коэффициентом объёмного сжатия ... ,
представляющим собой относительное изменение объёма жидкости ...
при изменении давления Р на единицу
...
где... - первоначальный объём жидкости, ...
... - изменение объёма ... при увеличении давления на
величину ...
Модулем объёмной упругости жидкости ... называется величина,
обратная коэффициенту объёмного сжатия... .
Для воды при атмосферном давлении он составляет около 2000 МПа.
При повышении давления на 0.1 МПа объём воды уменьшается
всего лишь на ... первоначального объёма.
Коэффициент объёмного сжатия для других капельных жидкостей
такого же порядка, поэтому в большинстве случаев сжимаемостью
капельных жидкостей можно пренебречь.
Температурное расширение
Это свойство жидкостей изменять свой объём. Характеризуется ко-
эффициентом температурного расширения ..., представляющим собой
относительное изменение объёма жидкости ... при изменении темпера-
туры ... на 1 С и постоянном давлении
...
Коэффициент температурного расширения ...при .. = 20 С и
давлении ...Па:
для воды 0.00015 С
для спирта0.00110 С
для нефти 0.00060 С
Вязкость - это способность жидкости оказывать сопротивление
скольжению одного слоя относительно другого. Силы, возниающие при
скольжении слоёв, называют силами внутреннего трения или силами
вязкости. Появление их обусловлено наличием межмолекулярных связей
между движущимися слоями. Вязкость характеризует степень подвиж-
ности частиц жидкости или текучести.
Согласно гипотезе, высказанной впервые Ньютоном в 1686 году,
а затем экспериментально обоснованной профессором Н.И.Петровым в
1863 году, силы внутреннего трения, возникающие между соседними дви-
жущимися слоями жидкости, прямо пропорциональны градиенту скорос-
ти, площади трущихся слоёв и зависит от свойств жидкости, т.е.
...
или
...
гдеТ - сила трения
... - площадь поверхности трущихся слоёв
... - динамический коэффициент вязкости
... - касательное напряжение
... - градиент скорости
Из соотношения для силы трения можно определить динамическую
вязкость
...
В гидравлических расчётах часто используется кинематическая
вязкость, равная отношению динамической вязкости ... к плотности
... жидкости:
...
Вязкость жидкостей зависит от температуры. С увеличением тем-
пературы вязкость капельной жидкости уменьшается, а вязкость га-
зов, наоборот, возрастает.
Кинематическая вязкость воды
при... = 20 имеет значение101 ...
при... = 40 имеет значение 66 ...
при... = 60 имеет значение 48 ...
Вязкость жидкостей измеряют с помощью приборов - вискозимет-
ров.
Для неньютоновских (бингемовских) жидкостей соотношение между
касательными наряжениями ... и градиентом скорости .... имеет вид
...
... - касательное напряжение в состоянии покоя.
Движение вязкопластических жидкостей начинается лишь после то-
го, как внешней силой преодолено сопротивление сдвига ... .
3. Поверхностное натяжение
Молекулы жидкости, находящиеся на свободной поверхности (тре-
ние, раздела жидкость - газ или жидкость - пар), испытывают одно-
стороннее воздействие со стороны соседних молекул. Поэтому на кри-
волниейной поверхности должны возникать растягивающие усилия. Для
количественного описания этого явления ещё в 1805 году Юнгом была
проведена классическая аналогия с упругой плёнкой. Натяжение этой
плёнки, т.е. усилие, приходящееся на единицу длины поперечного
разреза плёнки, характеризуется коэффициентом поверхностного натя-
жения
...
Сила поверхностного натяжения стремится сократить площадь
свободной поверхности. Их действие впервые обнаружено в капилярах,
поэтому эти силы до сих пор часто называют капилярными.
Величина ... зависит прежде всего от природы контактирующих
сред. Числовые значения его для некоторых пар приведены в таблице.
Таблица
┌──────────┬──────────────┬────────────┬────────────────────────┐
│ Вещество │Контактирующая│Температура,│Коэф. пов. натяжения│
│ │ среда │ К │... │
├──────────┼──────────────┼────────────┼────────────────────────┤
│ Вода │ Воздух │ 293 │ 72,8 │
│ ││ │
│ Жидкий │Пар. вещест-│ 373 │ 58,8 │
│ водород │ ва ││ │
│ │ ││ │
│ Жидкий │ │ 21 │ 2,0 │
│ кислород │ то же │ 91 │ 13,0 │
└──────────┴──────────────┴────────────┴────────────────────────┘
Коэффициент поверхностного натяжения ... падает с ростом тем-
пературы и практически не зависит от давления. Поверхностное натя-
жение может быть существенно снижено с помощью поверхностно-актив-
ных веществ, к числу которых относятся моющие средства.
Величина ... может служить мерой свободной энергии, которой
обладает граница раздела:
...
где ... - площадь свободной поверхности.
В этом случае
...
что согласуется с ранее указанной размерностью.
Существование поверхностного натяжения должно приводить к
возникновению на криволинейной поверхности перепада давлений, ко-
торые будут зависеть от конкретной геометрии поверхности.
Для объяснения этого факта рассмотрим равновесие элемента не-
плоской поверхности с линейными размерами ... и ... и главными ра-
диусами кривизны ... и ... соответственно.
...
...
Равнодействующие сил поверхностного натяжения, действующих на
границе выделенного контура, равны ... и ..., а возникающая вслед-
ствие этого сила, действующая по нормали к выделенной площадке, в
первом приближении равна
...
С учётом того, что
...
имеем выражение для силы
...
Эта величина, очевидно, и есть скачок давления на поверхности
раздела двух сред, обусловленный поверхностным натяжением.
Обозначив теперь через ... и ... давление в средах на границе
раздела из условия равновесия элементарной площадки, запишем соот-
ношение
...
которое называется формулой Лапласа.
Для цилиндрических поверхностей с круговым поперечным сечени-
ем радиуса ... имеем ... = ..., ... = ... и формула Лапласа прини-
мает вид:
...
В случае сферических поверхностей ... = ... = ... и тогда
получаем:
...
Если радиус сферической полости мал, то давления, развиваемое
поверхностным натяжением, могут стать значительными.
...
...
Весьма характерной является система газ - жидкость - твёрдая
стенка. В этом случае вводят значение краевого угла (угла ...
или угла смачивания).
Характерные значения краевых углов приведены в таблице
Таблица
┌─────────────┬───────────────────────────┬──────────────────────┐
│ Тв. вещества│ Жидкость │ Кр. угол, град │ │ │ │
├─────────────┼───────────────────────────┼──────────────────────┤
│ Сталь │ Вода │ 70 - 90 │
│ │ Жидкий водород │ 0│
│ │ Жидкий кислород │ 0│
│ Сталь │ Ртуть │ 128...148 │
└─────────────┴───────────────────────────┴──────────────────────┘
Если ... , жидкость называется смачивающей, если ...
- несмачивающей.
Высота подъёма или опускания жидкости в капиляре определяется
с помощью соотношения
...
где ... - диаметр капиляра, а ... - угол смачивания.
Уравнение состояния воды. Адиабата Тэйда
Опыт показывает, что между основными параметрами, характери-
зующими состояние газа (давление, плотность, температура) сущес-
твует определённая зависимость.
Уравнение
...
устанавливающее связь между этими параметрами, называется
уравнением состояния.
Поэтому состояние любого газа определяется двумя параметрами
(например, плотностью и температурой), так как третий параметр
(давление) можно найти из уравнения состояния.
Для идеального газа уравнение состояния можно представить в виде
...
где ... - газовая постоянная, зависящая от
относительной молекулярной массы ... . Для воздуха ... = 29,
... = 287 ... .
Существенное отклонение свойств воздуха от свойств идеального
газа наблюдается при высоких давления и низких температурах. На
состояние газа влияют такие процессы, как диссоциация и ... .
Уравнение состояния воды
Пусть в равновесном состоянии справедливо уравнение
. Тогда при малых отклонениях параметров Р и Т от
... и ... уравнение состояния воды в линейном приближении можно
записать в форме, предложенной Буссинеском:
...
где - коэффициент изотермической сжимае-
мости
- коэффициент теплового расширения
При температуре 293 К
...
Зависимость ... от давления весьма стойкая.
Адиабатические процессы, характеризующиеся отсутствием внеш-
него подвода или отвода тепла, протекают в воде практически при
постоянной температуре. Это объясняется особенностью молекулярно-
го строения жидкости. Ввиду большой плотности упаковки молекулы
жидкости помимо обмена импульсами в ... движении испы-
тывают дополнительные силы отталкивания. При сжатии жидкости даже
без нагревания развивается большое внутреннее давление нетеплово-
го происхождения. Изменение давления происходит только в результа-
те давления происходит только в результате изменения его механи-
ческой компоненты.
В случае значительных изменений давления связь между плот-
ностью и давлением становится существенно нелинейной. Наиболее
широкое распространение получило эмпирическое уравнение ...
, которое носит название уравнения Тэйда:
...
где С и ... - константы ( С ... 3200 ... Па, ... = 7.15).
Уравнение Тэйда устанавливает зависимость плотности только
от давления. Это означает, что оно описывает баротропный процесс.
Тема 3
Кинематика течений жидкости
1. Два подхода к описанию движения сплошной среды.
Переменные Эйлера и Лагранжа.
2. Траектория. Линия (поверхность) тока.
3. Кинематика вихрей. Циркуляция скорости.
Кинематикой называется раздел механики, изучающий движение
материальных тел в пространстве с геометрической точки зрения без
выяснения причин его возникновения. Все кинематические величины,
характеризующие движение твёрдого тела и движение отдельных точек
(расстояния, скорости, ускорения и т.д.), рассматриваются как фун-
кции времени.
1. Два подхода к описанию движения сплошной среды.
Переменные Эйлера и Лагранжа
Для описания движения сплошной среды возможны два подхода.
Один из них называется лагранжевым, другой - эйлеровым.
Лагранжев метод описания движения относится к типу отсчётных.
В некоторый (начальный) момент времени ... каждая из жидких частиц
маркируется путём присвоения ей значения координат в данный момент
времени.
В трёхмерном пространстве введём обозначения
...
В дальнейшем прослеживается движение каждой частицы индивиду-
ально. При таком подходе положение частицы в каждый момент времени
... будет зависеть от параметров а,б,с и ..., которые назы-
ваются переменными Лагранжа. Можно записать, что вектор положения
жидкой частицы равен
...
Скорость жидкой частицы выразится через производную ради-
ус-вектора
...
а ускорение через производную скорости
...
В последних двух формулах при дифференцировании параметры
а,б,с являются постоянными, ... и ... являются только функционала-
ми времени и в этом случае энергии дифференцирования ... и ...
тождественны.
Эйлеров метод описания движения относится к типу простран-
ственных. В каждой точке пространства с координатами ...
изучаются параметры движения в различные моменты времени ... .
Таким образом, скорость жидкости в различных точках пространства
должна быть функцией четырёх переменных ... , называемых
переменными Эйлера,
...
а её дифференциал
...
В движущейся среде приращения ... не ...
независимыми, а соответственно равны
...
Поэтому справедливо равенство
...
где
...
Это означает, что полное ускорение ... индивидуальной жид-
кой частицы, находящейся в момент времени ... в точке пространства
с координатами ... , состоит из двух частей: локального ускоре-
ния ..., обусловленного изменением скорости во времени в данной
точке, и конвективного ускорения ..., обусловленного неоднород-
ностью поля скоростей в окрестности данной точки и связанного с
этим обстоятельством конвективного переноса.
Производная ...носит название индивидуальной или субстан-
циональной производной.
Если ... , поле скоростей стационарно, однако это ещё
не означает, что в жидкости отсутствуют ускорения. Стационарность
или нестационарность поля скоростей зависит от выбора системы ко-
ординат.
Если ... = 0, поле скоростей однородно.
2. Траектория. Линия (поверхность) тока
Траекторией жидкой частицы называется геометрическое место
точек пространства, через которое частица последовательно проходит
во времени.
В переменных Лагранжа траекторию определяет уравнение
...
Если задача решена в переменных Эйлера, то известно поле ско-
ростей ... и траекторию следует находить путём решения
дифференциального уравнения
...
с начальным условием: при ... .
Линией тока называется линия, в каждой точке которой в каждый
момент времени скорость направлена по касательной к этой линии.
В векторной форме условие тангенциальности можно записать в
виде
...
В проекциях на оси координат получим систему уравнений
...
которую можно переписать также в виде
...
Время здесь является фиксированным параметром.
В стационарном случае траектория и линия тока совпадают. В
нестационарных течениях траектории отличаются от линий тока.
Поверхность тока определяется как поверхность, в каждой точке
которой в фиксированный момент времени вектор скорости лежит в ка-
сательной плоскости. Такую поверхность можно образовать, например,
путём проведения через замкнутую кривую непрерывной совокупности
линий тока. В этом случае говорят о трубке тока.
2. Кинематика вихрей
Рассмотрим вектор вихря скорости, который определяется соот-
ношением
...
называемый иногда вектором завихренности.
Линии в потоке жидкости, в каждой точке которой вектор вихря
скорости является касательным к данной линии, называются вихревыми
линиями.
...
...
Обобщение данного понятия на поверхность (вектор вихря в каж-
дой точке поверхности должен лежать в касательной плоскости) даёт
понятие вихревой поверхности или вихревого слоя.
Совокупность вихревых линий,проведенных через замкнутый кон-
тур, образует вихревую поверхность, а жидкость, заключённая внутри
вихревой поверхности, - вихревую трубку.
Интенсивность вихревой трубки удобнее выразить через циркуля-
цию вектора скорости Г.
В общем случае Г определяется как
...
где ... - вектор перемещения вдоль произвольного контура, со-
единяющего точки А и Б.
Если контур замкнут, то
...
Тема 4
Система уравнений гидростатики.
Динамика течений невязкой (идеальной) жидкости
1. Уравнение неразрывности.
2. Уравнение Эйлера.
3. Уравнение адиабатического движения жидкости.
4. Уравнения Эйлера в форме Громеки.
5. Гидростатика.
6. Уравнение Бернулли.
Система уравнений, описывающих течение жидкостей и газов, ос-
новывается на фундаментальных законах сохранения. К ним относятся
законы сохранения массы, количества движения, энергии.
Уравнения записываются в интегральной или дифференциальной
форме в зависимости от типа решаемой задачи.
Рассмотрим систему уравнений, которая описывает динамику те-
чений невязкой (идеальной ) жидкости.
Идеальной называется жидкость, у которой нет трения, т.е.
жидкие элементы, могут свободно перемещаться в касательном направ-
лении один относительно другого. В такой жидкости отсутствует теп-
лообмен между различными её участками, а тангенциальные и нормаль-
ные силы внутреннего трения не возникают.
В идеальной жидкости существуют силы только нормального да-
вления, однозначно определяемые её плотностью и температурой. Иде-
альная жидкость - абстракция, которой можно пользоваться на прак-
тике, если скорости изменения деформации в жидкости малы. Посколь-
ку касательные напряжения связаны с понятием вязкости, можно ут-
верждать, что идеальная жидкость - это невязкая жидкость.
Движение идеальной жидкости будем рассматривать в поле сил,
характеризуемых объёмной плотностью на единицу объёма жидкости.
1. Уравнение неразрывности
Вывод основных гидродинамических уравнений начнём с вывода
уравнения неразрывности, выражающего закон сохранения в гидродина-
мике.
Математическое описание состояния движущейся жидкости осуще-
ствляется с помощью функций, определяющих распределение скоростей
... и каких-либо двух термодинамических вели-
чин, например, ... - давления и ... -
плотности.
Скорость, давление и плотность жидкости будем относить к дан-
ным точкам пространства, а не к определённым частицам жидкости,
передвигающимся во времени и в пространстве. То есть будем пользо-
ваться переменными Эйлера.
...
...
Рассмотрим некоторый объём ... пространства. Количество (мас-
са) жидкости в этом объёме есть
...
Через элемент поверхности ..., ограничивающей рассматриваемый
объём, в единицу времени протекает количество ........ жидкости.
Вектор ... по абсолютной величине равен площади элемента по-
верхности и направлен по внешней нормали к ней. Тогда ...
положительно, если жидкость вытекает из объёма, и отрицательно, ес-
ли жидкость втекает в него.
Полное количество жидкости, вытекающей в единицу времени из
объёма ...
...
где ... - поверхность, ограничивающая выделенный объём ... .
С другой стороны, уменьшение количества жидкости в объёме ...
можно записать в виде
...
Приравнивая оба выражения, получаем:
...
Интеграл по поверхности преобразуем в интеграл по объёму
...
Таким образом,
...
Поскольку это равенство должно иметь место для любого выде-
ленного объёма, то должно быть равным нулю подынтегральное выраже-
ние, т.е.
...
Получили уравнение неразрывности.
... выражение ... можно записать
...
В декартовых координатах
...
Вектор
...
называют плотностью потока жидкости.
Его направление совпадает с направлением движения жидкости,
а абсолютная величина определяет количество жидкости, протекающей
в единице времени через единицу площади, расположенной перпендику-
лярно к скорости.
2. Уравнения Эйлера
Выделим в жидкости конечный объём. Полная сила, действующая
на выделенный объём жидкости, равна интегралу
...
взятому по поверхности рассматриваемого объёма. Преобразуем
его в интеграл по объёму, имеем
...
Отсюда видно, что на каждый элемент объёма ... жидкости дей-
ствует со стороны окружающей его жидкости сила - ....
Тогда на единицу объёма жидкости действует сила ... .
Мы можем теперь написать уравнение движения элемента объёма
жидкости, приравняв силу ... произведению массы ... еди-
ницы объёма жидкости на её ускорение
... (1)
Стоящая здесь производная ... определяет не изменение скорос-
ти жидкости в данной неподвижной точке пространства, а изменение
скорости определённой передвигающейся в пространстве частицы жид-
кости. Эту величину необходимо выразить через величины, относящи-
еся к неподвижным в пространстве точкам.
Изменение скорости ... данной жидкой частицы в течение време-
ни ... складывается из двух частей:
- из изменения скорости в данной точке пространства в течение
времени ...
- и из разности скоростей (в один и тот же момент времени) в
двух точках, разделённых расстоянием ..., пройденным рас-
сматриваемой частицей в течение времени ... .
Первая из этих частей равна
...
где производная берётся ...при постоянных ...,
т.е. в заданной точке пространства.
Вторая часть изменения скорости равна
...
Таким образом,
...
или, разделив обе скорости равенства на ...
...
Подставив полученное соотношение в (1), получим
...
Полученное уравнение движения жидкости - уравнение Эйлера
(1755), и является одним из основных в гидродинамике.
Если жидкость находится в поле тяжести, то на каждую единицу
её объёма действует ещё сила ..., где ... есть ускорение силы
тяжести. Эта сила должна быть прибавлена к правой стороне уравне-
ния и уравнение принимает вид:
...
При выводе уравнений движения мы совершенно не учитывали про-
цессов диссоциации энергии, которые могут иметь место в текущей
жидкости вследствие внутреннего трения (вязкости) в жидкости и
теплообмена между различными её участками.
Отсутствие теплообмена между отдельными участками жидкости
означает, что движение происходит адиабатически. Таким образом,
движение идеальной жидкости следует рассматривать как адиабатичес-
кое.
При адиабатическом движении энтропия каждого участка жидкости
остаётся постоянной при перемещении последнего в пространстве.
Обозначая ... энтропию, отнесённую к единице массы жидкости, мы
можем выразить адиабатичность движения уравнением
...
полная производная по времени означает изменение энтропии
заданного перемещающегося участка жидкости. Эту производную можно
записать в виде
...
Это есть общее уравнение, выражающее собой адиабатичность
движения идеальной жидкости. С помощью уравнения неразрывности
его можно написать в виде уравнения неразрывности для энтропии.
...
где ... - плотность потока энтропии.
Иногда это условие используют в более простой форме. Если
в некоторый момент времени энтропия одинакова во всех точках
объёма жидкости, то она остаётся везде одинаковой и неизменной
со временем и при дальнейшем движении жидкости.
В этих случаях уравнение адиабатичности записывается в виде
...
Изэнтропичностью движения можно воспользоваться и предста-
вить уравнения Эйлера в другом виде. Из термодинамических соотно-
шений известно
...
... - тепловая функция единицы массы жидкости,
... - удельный объём, Т - температура.
Поскольку...., имеем просто
...
и поэтому
...
Уравнения Эйлера можно записать в виде
...
Воспользуемся известной формулой векторного анализа
...
уравнение Эйлера можно записать в другом виде
...
К уравнениям движения необходимо добавить граничные условия,
которые должны выполняться на ограничивающих жидкость границах.
Для идеальной жидкости это условие должно выражать собой просто
тот факт, что жидкость не может проникнуть за твёрдую поверхность.
На неподвижных стенках это означает, что должна обращаться в
нуль нормальная к стенке компонента вектора скорости:
...
3. Гидростатика
Для покоящейся жидкости, находящейся в однородном поле тяжес-
ти, уравнение Эйлера принимает вид
...
Это уравнение описывает механическое равновесие жидкости.
Если внешние силы вообще отсутствуют, то уравнения равновесия
дают
...
т.е. ... .
- давление одинаково во всех точках жидкости.
Притом плоскость жидкости постоянна во всём объёме. Направим
ось ... вертикально вверх, имеем
...
Откуда
...
Если покоящаяся жидкость имеет свободную поверхность (на
высоте ...), к которой приложено одинаковое во всех точках внешнее
давление ..., то эта поверхность должна быть горизонтальной плос-
костью ....
...
...
Из условия ...при ...имеем
...
так что
...
4. Уравнение Бернулли
Уравнения гидродинамики заметно упрощаются в случае стацио-
нарного течения жидкости. Под стационарным (или установившимся)
подразумевают такое течение, при котором в каждой точке простран-
ства, занятого жидкостью, скорость течения остаётся постоянной во
времени. Скорость ... остаётся функцией только координат
...
...
Рассмотрим некоторые сведения о линиях тока. Линии тока это
линии, касательные к которым указывают направление вектора скорос-
ти в точке касания в данный момент времени. Уравнения линий тока
определяются системой дифференциальных уравнений
...
При стационарном движении жидкости линии тока остаются неиз-
менными во времени и совпадают с траекториями частиц жидкости.
При нестационарном течении такое совпадение не имеет места:
- касательные к линии тока дают направление скорости разных
частиц жидкости в последовательных точках пространства в
определённый момент времени
- касательные к траектории дают направление скорости опреде-
лённых частиц в последовательные моменты времени.
Умножим уравнение Эйлера для стационарного потока жидкости на
единичный вектор касательной к линии тока в каждой её точке ... .
Проекция градиента на некоторое направление равна производ-
ной, взятой по этому направлению. Поэтому
...
Вектор ... перпендикулярен вектору скорости, и поэтому
его проекция на направление ... равна нулю
...
Таким образом получаем
...
Откуда следует, что величина ...постоянна вдоль линии тока
...
Значение ... , вообще говоря, различно для разных линий то-
ка. Это уравнение называют уравнением Бернулли.
Если течение жидкости происходит в поле сил тяжести, то в
правой части уравнений Эйлера есть ускорение силы тяжести ... .
Выберем направление силы тяжести в качестве направления оси
..., причём положительные значения ... отсчитываются вверх. Тогда
проекция ... на ... есть
...
Соответственно этому будем иметь
...
Таким образом, уравнение Бернулли гласит, что вдоль линий
тока остаётся постоянной длина
...
Тема 5
Потенциальные и несжимаемые течения
1. Сохранение циркуляции.
2. Потенциальное движение.
3. Несжимаемая жидкость.
1. Сохранение циркуляции скорости
Интеграл
...
взятый вдоль замкнутого контура, называют циркуляцией ско-
рости вдоль этого контура.
Рассмотрим некоторый замкнутый контур, проведенный в жидкос-
ти в некоторый момент времени. Будем рассматривать его как "жид-
кий", составленный из находящихся на нём частиц жидкости. С тече-
нием времени контур перемещается.
Вычислим производную по времени от циркуляции скорости с учё-
том подвижности контура. Временное дифференцирование по координа-
там обозначим знаком ..., знак ... - дифференцирование по времени.
Будем учитывать, что меняются скорость и сам контур.
...
По определению скорость ... это производная радиус-вектора
...
...
Интеграл по замкнутому контуру от полного дифференциала ра-
вен нулю и остаётся
...
Из уравнений Эйлера имеем
...
Применим формулу Стокса, получаем тогда (поскольку ........)
...
Таким образом, переходя к прежним обозначениям, находим око-
нчательно:
... или ...
Мы приходим к результату, что в идеальной жидкости циркуляция
скорости вдоль замкнутого контура остаётся неизменной со временем.
Это утверждение называется теоремой Томсона или законом со-
хранения циркуляции скорости. Соотношение получено путём использо-
вания уравнений Эйлера с использованием предположения об изэнтро-
пичности движения жидкости.
Применим теорему Томсона к бесконечно малому замкнутому кон-
туру ... и преобразовав интеграл по теореме Стокса, получим:
...
где ... - элемент поверхности, опирающейся на контур ...
Вектор ....... часто называется завихренностью течения жидкости в
данной её точке. Постоянство произведения ......................
... можно использовать, сказав, что завихренность переносит-
ся вместе с движущейся жидкостью.
2. Потенциальное движение
Движение жидкости, при котором во всём пространстве
...
называется потенциальным (или безвихревым) в противополож-
ность вихревому движению, при котором ротор скорости отличен от
нуля.
Таким образом, мыпришли бы к выводу, что стационарное обте-
кание взятого тела натекающим из бесконечности однородным потоком
должно быть потенциальным. Поскольку на бесконечности натекающий
поток однороден, его скорость ... , так что ... = 0 на всех
линиях тока.
Однако, ввиду наличия стенки нельзя провести в жидкости за-
мкнутый контур, который охватывал бы такую линию тока.
...
...
В результате возникает картина течения, характеризующаяся
наличием отходящей от тела "поверхности тангенциального разрыва",
на которой скорость жидкости терпит разрыв непрерывности.
Как и всякое векторное поле с равным нулю ротором, скорость
потенциально движущейся жидкости может быть выражена в виде гра-
диента от некоторого скаляра, называемого потенциалом скорости
...
...
Напишем уравнения Эйлера в виде
...
и подставив в него ........, получаем
...
Откуда находим следующе равенство
...
где ... произвольная функция времени. Это равенство представ-
ляет собой первый интеграл уравнений потенциального движения.
При стационарном движении имеем ... = 0, ......... и интеграл
переходит в уравнение Бернулли
...
Отметим существенные отличия между уравнениями Бернулли в
случае потенциального и непотенциального движения. ..... в правой
части этого уравнения есть величина, постоянная вдоль каждой линии
тока, но вообще говоря, различная для разных линий тока.
При потенциальном же движении ... в уравнении Бернулли есть
величина, постоянная во всём объёме жидкости.
3. Несжимаемые жидкости
Для плоских течений жидкостей их плотность можно считать по-
стоянной вдоль всего объёма жидкости в течение всего времени движе-
ния. Такое движение называется движением несжимаемой жидкости.
Общие уравнения гидродинамики для несжимаемой жидкости упро-
щаются. Уравнение неразрывности при ......... принимает простой
вид
...
уравнения Эйлера не меняют своего вида, запишем их в виде
...
Для несжимаемой жидкости тепловая функция записывается следу-
ющим образом
...
Тогда уравнение Бернулли для несжимаемой жидкости имеет вид
...
Особенно упрощается уравнение для потенциального течения не-
сжимаемой жидкости.
При подстановке ........ в уравнение неразрывности
..... = 0, получим
...
то есть уравнение Лапласа для потенциала ... .
Граничные условия. К этому уравнению должны быть добавлены
граничные условия на поверхности соприкосновения жидкости с твёр-
дыми телами:
- на неподвижных твёрдых поверхностях нормальная к поверхнос-
ти компонента ... скорости жидкости должна быть равна нулю, для
движущихся тел ... должна быть равна проекции скорости движения те-
ла на направление той же нормали.
С другой стороны, скорость ... равна производной от потенциала
... по направлению нормали
...
Таким образом, граничные условия гласят в общем случае, что
.... является на границах заданной функцией координат и времени.
При потенциальном движении скорость связана с давлением для
несжимаемой жидкости соотношением
...
Если движение жидкости является потенциальным и вызвано дви-
жением некоторого тела то уравнение Лапласа не содержит явно вре-
мени, время входит в решение через граничные условия.
Из уравнения Бернулли ..................... видно, что при
стационарном движении несжимаемой жидкости вне поля тяжести на-
ибольшее значение давления достигается в точках, где скорость
обращается в нуль. Такая точка обычно имеется на поверхности обте-
каемого жидкостью тела (точка О) и называется критической точкой.
Если ... - скорость набегающего на тело потока жидкости (скорость
на бесконечности), а ... - давление в критической точке равно
...
Если распределение скоростей в движущейся жидкости зависит
только от двух координат, то о таком течении говорят как о двумер-
ном или плоском. Для решения задач о двумерном течении несжимаемой
жидкости иногда удобнее использовать функцию тока. Из уравнения не-
разрывности
...
видно, что компоненты скорости могут быть записаны в виде про-
изводных
...
от некоторой функции ..., называемой функцией тока. Урав-
нение неразрывности при этом удовлетворяется автоматически.
...
Зная функцию тока, можно непосредственно определить форму ли-
ний тока для стационарного движения жидкости. Дифференциальное
уравнение линий тока
...
или
...
оно выражает условие параллельности касательной к линии тока и
направления вектора скорости.
Подставляя сюда выражение для скоростей через функцию тока
...
откуда ........ Таким образом, линии тока представ-
ляют собой семейство кривых, получающихся приравниванием функции
тока ... постоянной.
Если между точками 1 и 2 в плоскости .... провести кривую,
то поток жидкости ... черезэту кривую определится разностью зна-
чений функции тока в этих точках независимо от формы кривой.
Действительно, если ... - проекция скорости на нормаль к кри-
вой в данной точке, то
...
или
...
Мощные методы решения задач о простом потенциальном обтекании
несжимаемой жидкостью различных профилей связаны с применением к
ним теории функций комплексной переменной.
Тема 6
ГИДРОСТАТИКА
1. Силы, действующие на жидкость. Давление.
Единицы измерения давления.
2. Закон Паскаля.
3. Дифференциальное уравнение равновесия жидкости.
4. Виды давления ( барометрическое, абсолютное, избыточное,
манометрическое ).
5. Приборы для измерения давления.
6. Сила давления жидкости на плоскую стенку.
7. Простейшие гидравлические машины.
8. Закон Архимеда.
9. Равновесие и остойчивость тел, полностью
погруженных в жидкость.
1. Силы, действующие на жидкость. Давление.
Единицы измерения давления
Рассечем жидкость, находящуюся в объеме ... (например,
сосуде) некоторой поверхностью на две части I и II.
Рассмотрим жидкость в объеме I. Все, что окружает этот
объем, отбросим (дно, боковые стенки и т.д.) и действие
отброшенного заменим соответствующими силами. Эти силы называются
поверхностными.
Кроме них на жидкость действуют еще массовые силы (силы
тяжести и инерции), которые пропорциональны массе тела.
Выделим из жидкости некоторый объем. Возьмем на поверхности
этого объема бесконечно малую площадку ... . Hа эту площадку
действует поверхностная сила ... . Разложим эту силу на
нормальную ... и касательную ...
Hормальная сила, приходящаяся на единицу площади,
называется давлением и обозначается буквой ... , т.е.
...
Измеряется давление в ....
Сила трения (касательная сила), приходящаяся на единицу
площади, обозначается буквой ..., т.е.
...
Сила трения обычно пропорциональна градиенту скорости ... .
Для жидкости, находящейся в равновесии (в покое), сила трения
равна нулю, так как в этом случае ... .
2. Закон Паскаля
Если в жидкости взять любую точку, то на основании основного
уравнения гидростатики
...
давление в этой точке равно давлению, приложенному к
свободной поверхности, плюс ..., где ... - глубина точки.
Таким образом
Закон Паскаля
Давление, приложенное к свободной поверхности, передается во
все точки жидкости без изменения.
3. Дифференциальное уравнение равновесия жидкости
Запишем уравнение Эйлера
...
Если жидкость покоится
...
Дифференциальные уравнения равновесия жидкости в проекции на
оси декартовой системы координат могут быть записаны так
...
Здесь ... - проекции на оси ... сил, действующих на единицу
массы рассматриваемой жидкости.
Умножая давления соответственно на ... и складывая их,
получаем
...
Левая часть уравнения представляет полный дифференциал
...
следовательно, и правая часть должна быть также полным
дифференциалом, для этого необходимо и достаточно, при постоянном
..., чтобы существовала функция ... такая что
...
Имеем
...
Проинтегрировав, получим
...
где С - постоянная интегрирования.
Если в какой-либо точке известно давление ... и
постоянная функция ..., то ...
из интеграла имеем
...
В частности, когда жидкость находится в поле сил тяжести
...
Следовательно,
...
Уравнение для давления принимает вид
...
Свободная поверхность жидкости плоская ... . При равновесии
жидкости в поле земного тяготения поверхности уровня представляют
собой горизонтальные плоскости.
Рассмотрим примеры.
Пример 1.
Определить уравнение свободной поверхности жидкости в
сосуде, движущемся горизонтально с ускорением а.
Решение. Из нескольких сил на жидкость действуют сила
тяжести и сила инерции, т.е.
...
Имеем
...
откуда
...
- уравнение прямой.
Следовательно, свободная поверхность представляет собой
плоскость, наклоненную к горизонту под углом ..., который
определяется из равенства
...
Пример 2. Определить уравнение свободной поверхности
жидкости в сосуде, вращающемся вокруг вертикальной оси с угловой
скоростью ...
Решение.
Вследствие трения о стенки сосуда жидкость будет вращаться с
такой же угловой скоростью. Жидкость будет находиться в
относительном покое. Поэтому при решении задачи применимы
уравнения равновесия.
Из массовых сил на жидкость действует центробежная сила и
сила тяжести. Центробежная сила, действующая на массу ...,
находится на расстоянии ... от оси вращения
...
Проекции силы на оси, отнесенные к единице массы, будут
...
Тогда
...
Откуда
...
т.е. свободная поверхность - параболоид вращения.
4. Виды давления ( барометрическое, абсолютное, избыточное,
манометрическое )
Различают следующие виды давления: барометрическое,
абсолютное, манометрическое и вакуумметрическое.
Барометрическое (или атмосферное) давление ... зависит от
места над уровнем моря и от погоды. За нормальное барометрическое
давление принимают давление, равное
соответствует 101325 ... С высотой барометрическое давление
убывает. В глубоких шахтах барометрическое давление значительно
больше, чем на уровне моря.
Давление, вычисляемое по соотношению
...
называется абсолютным.
Аналогичное давление в точке равно сумме внешнего
поверхностного и весового давления.
Если к свободной поверхности приложено барометрическое
давление ..., то есть ... и основное уравнение гидростатики
перепишем так
...
Давление
...
носит название манометрического или избыточного. Таким
образом, манометрическим давлением называется разность между
абсолютным давлением ... и барометрическим ..., если ...
Если в данной точке жидкости абсолютное давление меньше
барометрического, то разность между барометрическим и абсолютным
давлениями называется вакуумметрическим давлением ...
Итак, если ... , то
...
Абсолютное давление ... отрицательным быть не может, поэтому
вакуумметрическое давление не может быть больше барометрического.
5. Приборы для измерения давления
Приборами для измерения барометрического давления служат
барометры различных конструкций.
Для измерения манометрического давления служит манометр.
Манометрическое давление можно измерить высотой столба жидкости.
Сосуд наполнен жидкостью с плотностью ... . Давление на свободной
поверхности ...
Пусть необходимо измерить давление на уровне 1-1. Если на
этом уровне сделать отверстие и присоединить к нему стеклянную
трубку П, то жидкость в этой трубе поднимется под действием
давления на некоторую высоту ... По основному уравнению
гидростатики
...
откуда
...
Этой высотой ... поднятия жидкости в трубке П можно измерять
манометрическое давление. Трубка П называется пьезометром.
Hайдем
соотношение между 1 ...,
При высоте
вод. столба ... =
...
При высоте
ртутного столба ... =
...
Для измерения вакуумметрического давления применяется
вакуумметр. Допустим, что требуется измерить вакуумметрическое
давление воздуха в сосуде ..., т.е. величину ..., где ... -
абсолютное давление в сосуде.
Присоединим к сосуду изогнутую трубку, опущенную в жидкость.
Применяя основное уравнение гидростатики для точки, расположенной
в трубке на уровне свободной поверхности жидкости в резервуаре,
получим
...
Так как
..., то
...
Вакуумметрическому давлению будет соответствовать высота
подъема ... жидкости в изогнутой трубке над уровнем в резервуаре.
6. Сила давления жидкости на плоскую стенку
Гидростатическое давление представляет собой систему
параллельных сил, действующих в одну сторону и перпендикулярных к
плоскости стенок.
Такая система приводится к одной силе - равнодействующей,
равной арифметической сумме всех сил и приложенной в центре
параллельных сил. Для определения равнодействующей давлений,
приложенных к площадке ..., плоскость которой ... наклонена к
горизонту под углом ..., возьмем начало координат в плоскости
приведенного уровня на линии пересечения с плоскостью площадки,
приняв линию пересечения за ось ... на направив ось ...
вертикально вниз, кроме того в плоскости площадки возьмем
вспомогательные оси ... и ..., совместив ... и ... .
...
Последний интеграл равен площади площадки ..., умноженной на
координату центра тяжести ... .
...
Произведение ... выражает объем цилиндрического столба с
основанием ... и высотой ... и мы приходим к выводу, что давление
тяжелой жидкости на плоскую площадку измеряется весом
цилиндрического столба этой жидкости, который был бы расположен
над площадкой, если бы она лежала горизонтально на глубине своего
цента ... .
Сосуды различной формы, но с одинаковой площадью дна,
наполненные жидкостью на одну и ту же высоту H, имеют одинаковую
силу давления на дно.
7. Простейшие гидравлические машины
Жидкости практически несжимаемы и равномерно передают
давление по всему объему. Это свойство широко используется в
различных отраслях техники (гидроприводы, гидроавтоматика,
гидравлические тормоза, усилители и т.д.).
Принцип их работы основан на следующем: имеются два
соотносящиеся между собой цилиндра разного диаметра.
Приложим к поршню меньшего из цилиндров какую-то внешнюю
силу ..., мы тем самым создаем на поверхности жидкости давление
...
Это давление равномерно передается во все точки
пространства, заполненного жидкостью. Тогда на поршень большего
цилиндра будет действовать сила
...
Таким образом, чем больше разняться между собой площади
поперечного сечения цилиндров, тем большую силу мы будем получать
в таких гидравлических устройствах.
8. Закон Архимеда
Определим силу давления жидкости на погруженное тело А
объемом ...
Представим, что в жидкости выделен объем, точно такой же,
как и тело А. Этот объем жидкости находится в равновесии под
действием двух сил:
1) силы давления жидкости ... на поверхность выделенного
объема,
2) силы тяжести жидкости, равной ... и направленной
вертикально вниз.
Следовательно, сила ... равна силе тяжести выделенного
объема жидкости, направленная в обратную сторону, то есть
вертикально вверх, и приложена в центре объема, т.е. в той же
точке, в которой приложена сила тяжести выделенного объема
жидкости.
Точка ... называется центром водоизмещения.
Закон Архимеда.
Сила давления жидкости на погруженное в нее тело приложена в
центре водоизмещения, направлена вертикально вверх и равна силе
тяжести жидкости, вытесненной телом.
...
Сила ... называется архимедовой силой, ... - объемным
водоизмещением, а ... - водоизмещением.
9. Равновесие и остойчивость тел, полностью
погруженных в жидкость
Если сила тяжести ... тела А больше архимедовой силы ..., то
равнодействующая этих сил (... и ... ) направлена вниз и
заставляет тело опускаться на дно. Таким образом, если ..., тело
тонет.
Если сила тяжести ... тела меньше архимедовой силы ..., то
равнодествующая этих сил (... и...) направлена вертикально вверх
и заставляет тело подняться на поверхность. При выходе части тела
из жидкости сила давления на оставшуюся погруженную часть тела
соответственно уменьшается, благодаря чему уменьшается и величина
направленной вверх равнодействующей, заставляющей тело всплывать,
в результате при некотором частичном погружении тела
устанавливается равновесие и тело оказывается плавающим на
поверхности жидкости.
Таким образом при ... тело всплывает на поверхность
жидкости.
Для того, чтобы тело не опускалось на дно и не всплывало,
необходимо, чтобы ... .
Остойчивостью плавающего тела называется его способность
возвращаться в первоначальное положение равновесия после
приращения силы, вызвавшей крен.
Возможны три случая
1) центр тяжести С лежит ниже центра водоизмещения ...,
2) центр тяжести С находится выше центра водоизмещения ...,
3) центр тяжести С совпадает с центром водоизмещения ....
В первом случае равновесие остойчивое, так как при ...
возникает пара сил, стремящаяся вернуть тело в первоначальное
положение.
Во втором случае равновесие неустойчивое, в третьем -
Тема 7
Анализ и применение уравнения Бернулли
1. Уравнение неразрывности в гидравлике. Расход.
2. Анализ уравнения Бернулли.
3. Энергетический смысл уравнения Бернулли.
4. Предел применимости уравнения Бернулии.
5. Примеры применения уравнения Бернулли.
5.1. Расходомер Вентури.
5.2. Измерение скорости (Трубка Пито).
5.3. Кавитация.
5.4. Формула Торичелли.
6. Уравнение неразрывности в гидравлике. Расход.
1. Расход. Уравнение неразрывности в гидравлике
Рассмотрим установившийся поток между живыми сечениями 1,2.
Живым сечением называется поверхность в пределах потока,
проведенная параллельно к направлению струек. За единицу времени
через живое сечение 1 в рассматриваемый объем жидкости
...
где ... - площадь живого сечения, ... - средняя скорость в
сечении.
Через живое сечение 2 за это время вытекает объем жидкости
...
где ... - площадь живого сечения 2, ... - средняя скорость в
сечении 2.
Поскольку форма объема 1-2 с течением времени не изменяется,
жидкость несжимаемая, объем жидкости ... должен равняться объему
вытекающему ...
Поэтому можно записать
...
Это уравнение называется уравнением неразрывности
Из уравнения неразрывности следует, что
...
Средние скорости обратно пропорциональны площадям
соответствующих сечений.
2. Анализ уравнения Бернулли
Запишем уравнение Бернулли для установившегося движения
идеальной сжимаемой жидкости при условии ее баротропности (...)
в поле массовых сил
...
проинтегрировав
...
Для потенциального течения константа уравнения Бернулли постоянна
для всей области течения. При вихревом движении идеальной
жидкости константа С в интеграле Бернулли сохраняет последнее
значение только для данной вихревой линии, а не для всего
пространства, как при безвихревом течении.
Уравнение Бернулли является одним из основных в
гидрогазодинамике, так как определяет изменение основных
параметров течения - давления, скорости и высоты положения
жидкости.
Проинтегрируем дифференциальное уравнение Бернулли для
конечного участка струйки 1-2
...
Интеграл ... выражает работу ... т.е. работу сил давления по
перемещению килограмма жидкости из области 1 с давлением ... в
область 2 с давлением ... .
а) ... б)
изохорный процесс
В зависимости от типа процесса (термодинамического) который
совершает жидкость, то есть от вида зависимости ...
Рассмотрим изобарный ....
...
изохорный процесс ...
...
Для несжимаемой жидкости при течении без обмена механической
работой с внешней средой, получим, при ... из уравнения Бернулли
...
или умножив на ...
...
или разделив на ...
...
где константы имеют следующий физический смысл:
С - полная механическая энергия массы жидкости объемом в
кубический метр или полный напор,
... или Па,
... - полная механическая энергия ..., ... жидкости или
полный напор в метрах столба данной жидкости.
Все три величины имеют одинаковый физический смысл, поэтому в
учебной и технической литературе можно встретиться с тем, что
любой из них присваивают название полного напора.
Составляющие полной механической энергии жидкости наиболее
наглядно изображаются и измеряются в метрах столба жидкости,
... - потенциальная энергия положения жидкости,
отсчитываемая от произвольно выбранной ... плоскости, или
геометрический напор, ..., Па, м;
... - потенциальная энергия или гидростатический напор, ...,
Па, м;
... - кинетическая энергия жидкости или скоростной (для
жидкостей) напор, ..., Па, м.
Пьезометрический напор ... может измеряться от полного
вакуума ... или, например, от давления окружающей среды. В обеих
частях равенств должно подставляться абсолютное или избыточное
давление.
Hачало отсчета энергии произвольно, но должно быть
одинаково для обеих частей равенств.
...
3. Энергетический смысл уравнения Бернулли
Заключается в утверждении закона сохранения полной
механической энергии единицы массы несжимаемой жидкости.
а) при потенциальном течении для любой точки пространства,
б) при вихревом - только вдоль вихревой линии тока и
элементарной струйки.
Этот закон иногда формулируется в виде теоремы трех высот -
в приведенных условиях сумма трех высот - геометрической,
пьезометрической и динамической сохраняет неизменное значение.
При этом составляющие полной энергии могут
взаимопревращаться.
Следует иметь в виду, что изменение кинетической энергии
несжимаемой жидкости вдоль элементарной струйки ... не может
задаваться произвольно: в соответствии с уравнением неразрывности
это изменение однозначно определяется изменением площади
поперечного сечения канала
...
Течение в горизонтальной струйке имеет большое практическое
значение, оно реализуется в ... двигателей, уравнение Бернулли
при ...
...
Итак, увеличение скорости несжимаемой жидкости в
горизонтальной элементарной струйке всегда сопровождается
уменьшением давления, а уменьшение скорости - увеличением
давления вплоть до ... при ... Поэтому скоростной напор широко
используется, например, для подачи воды в систему охлаждения,
разрушения горных пород и т.д.
В связи с тем, что скорость несжимаемой жидкости может
уменьшаться только вследствие изменения площади сечения, приходим
к важному выводу о том, что картина линий тока при течении
несжимаемой жидкости однозначно определяет не только изменение
скорости, но и статического давления: при сгущении линий тока
давление уменьшается, при расширении - увеличивается. Это правило
широко используется при анализе движения жидкости и ее
взаимодействии с телами.
4. Предел применимости уравнений неразрывности и
Бернулли.
При течении жидкости по каналу при постоянстве ..., и при
произвольно изменяемой площади 2. Казалось бы, что
...
Однако по уравнению Бернулли при ...
...
давление ... должно было бы принять значение минус
бесконечность, что лишено смысла: абсолютное давление не может
быть меньше нуля.
Таким образом уравнения неразрывности и Бернулли справедливы
лишь до тех пор, пока минимальное давление в потоке остается
большим нуля.
5. Примеры применения уравнения Бернулли.
Рассмотрим примеры применения уравнения Бернулли.
1. Расходомер Вентури.
Для определения скорости и расхода жидкости часто
используется расходомер Вентури. Измерим статическое давление ...
и ... в поперечных сечениях с различными площадями.
Интеграл Бернулли для сечений 1 и 2 принимает вид
...
Из уравнения равенства расходов для двух сечений 1 и 2 имеем
...
Для вычисления показания дифференциального манометра запишем
условие равновесия
...
Собирая все результаты, получаем
...
Формула используется для определения скорости в трубе. Hа
практике для повышения точности иногда вводят эмпирический
коэффициент, учитывающий гидравлические ... в трубке Вентури.
2. Измерение скорости.
Для измерения кинетической энергии используется трубка
полного давления, которая устанавливается в точке измерения
открытым концом против потока жидкости.
Струйка жидкости, подтекающая к открытому концу трубки,
полностью замораживается (...=0) и весь скоростной напор
превращается в давление, которое в сумме со статическим достигает
давления торможения ... в данной точке, которое называется полным
...
откуда
...
Таким образом измерение скорости жидкости или "несжимаемого"
газа (...) основано на сопоставлении давления торможения с
давлением в невозмущенном потоке. Последнее еще называется
статистическим давлением ... Приемником давления служит
Г-образная трубка, или трубка Пито. Давление обычно измеряют с
помощью ...-образной трубки, куда залита жидкость манометрическая
(спирт, вода, ртуть).
Приемное отверстие статического давления должно находится
не слишком далеко от входа в трубку Пито, чтобы не случилось
рассеивание механической энергии за счет вязкости, и не слишком
близко, чтобы присутствие трубки Пито не искажало статическое
давление.
3. Кавитация.
Hа практике оказывается, что в жидкости давление, равное
нулю, недостижимо. Если давление ..., снижаясь, достигает
давления паров этой жидкости, насыщающих пространство при данной
температуре ..., то начинается процесс образования пузырьков пара
(кипение), и неразрывность течения капельной жидкости нарушится.
...
...
Далее смесь капельной жидкости и пузырьков пара попадает в
расширяющийся канал, давление возрастает и пузырьки пара начинают
конденсироваться.
Кавитацией называется совокупность процессов образования
пузырьков пара и их конденсация.
Кавитация может возникать не только в трубопроводах, но и
при внешнем обтекании тел в областях, где возрастают местные
скорости и уменьшается давление. Кавитации подвержены быстроходные
колеса насосов и турбин, гребные винты.
Конденсация пузырьков пара происходит на твердых
поверхностях очень быстро и завершается гидравлическим ударом,
при котором развивается местное ударное давление на твердых
поверхностях, достигающее сотен и даже тысяч атмосфер. Поэтому
кавитация сопровождается тряской, шумом, снижением КПД насосов и
турбин, эрозией твердых поверхностей, а иногда и выходом из строя
агрегатов.
Обычно работа гидравлических систем в условиях кавитации не
достигаются. Для предотвращения кавитации минимальное давление
жидкости в системе должно быть больше давления паров, насыщающих
пространство.
Одним из способов предотвращения кавитации является снижение
температуры жидкости. Это приводит к снижению давления паров,
насыщающих пространство.
Hапример, вода при 373 К кипит при давлении ... Па, а при
193 К - ... Па. При кавитации многокомпонентных жидкостей
(керосин, бензин и т.д.) вначале вскипают легкие фракции, а затем
тяжелые. Конденсация происходит в обратном порядке.
Для оценки возможности возникновения кавитации используется
безразмерный критерий - число кавитации
...
Значение, числа кавитации при котором она возникает,
называется критическим ... .
Явление используется в кавитационных регуляторах расхода.
4. Формула Торичелли
Применим интеграл Бернулли для определения скорости
истечения несжимаемой тяжелой жидкости из большого открытого
сосуда через малое отверстие.
Здесь ... - площадь свободной поверхности, ... - площадь
отверстия, ... и ... - скорости на поверхности и в отверстии.
Уравнение неразрывности принимает вид
...
Считая движение жидкости установившимся и безвихревым
применим интеграл Бернулли
...
Откуда
...
Из уравнения неразрывности
...или ...
Если отношение ... мало, то пренебрегая членом ..., получаем
для скорости истечения приближенную формулу Торичелли.
Пример.Определить форму сосуда вращения, употребляемого для
водяных часов.
...
Используя уравнение Бернулли можно объяснить принцип
действия
1) работы струйного насоса, в котором высоконапорный поток .
.. используется для подачи жидкости ... из резервуара.
2) принцип наддува топливного самолетного бака для
предотвращения кавитации в топливной системе при полетах на
большой высоте.
...
3) причину повышения подъемной силы крыла при заданной
картине линий тока
...
Уменьшение давления в точках, где скорость потока больше,
положено в основу водоструйного насоса. Струя воды подается в
трубку, открывающуюся в атмосферу, так что на выходе их трубки
давление равно атмосферному. В трубке имеется сужение, по
которому вода идет с большой скоростью, вследствие чего давление
в этом месте оказывается меньше атмосферного. Такое же давление
устанавливается и в охватывающей трубку камере насоса, которая
сообщается с трубкой через разрыв, имеющийся в узкой части
трубки. Подсоединив к камере насоса откачиваемый объект, из него
можно откачать воздух (или какой-либо другой газ) до давления
порядка
воды и уносится в атмосферу.
...
Тема 8
Потери напора
1. Классификация потерь напора. Задачи гидродинамического
расчета.
2. Потери напора по длине.
2.1. Основное уравнение равномерного движения.
2.2. Два режима течения жидкости.
2.3. Профиль скорости при ламинарном и турбулентном режимах
течения.
2.4. Критерии режима течения жидкости.
2.5. Определение потерь напора на трение.
3. Местные гидравлические сопротивления. Формула Вейсбаха.
3.1. Внезапное расширение трубопровода.
4. Гидравлический расчёт напорных трубопроводов.
4.1. Классификация трубопроводов. Задачи гидравлического
расчёта трубопроводов.
4.2. Расчёт коротких трубопроводов.
4.3. Расчёт длинных трубопроводов при последовательном
соединении труб.
4.4. Расчёт трубопровода при параллельном соединении труб.
1. Классификация потерь напора и задач
гидродинамического расчёта
Потери напора делятся на два вида: потери по длине и местные
потери.
Потерями напора по длине называются потери удельной энергии
потока на преодоление сопротивления движения напора на участке
рассматриваемой длины без учёта влияния местных сопротивлений.
Местными потерями напора называют потери удельной энергии
потока на преодоление сопротивлений движению потока, вызванных
каким-либо местным препятствием (расширение, сужение потока,
задвижка, шейка, клапан, колено и т.д.).
Потери напора обозначаются буквой ... с индексом,
определяющим их вид.
Задачи гидродинамического расчёта:
1. Определение потерь напора.
2. Определение расхода.
2. Потери напора по длине
2.1. Основное уравнение равномерного движения
Рассмотрим прямолинейное равномерное движение жидкости.
Живые сечения в этом случае могут быть произвольной формы, но
не должны изменяться по всей длине рассматриваемого участка.
В таком потоке потери напора определяются лишь потерями по длине.
Выделим из потока участок жидкости длиной ... и запишем
уравнение Бернулли для сечений 1 и 2
...
...- ординаты центра тяжести сечений 1,2
...- давление в центрах тяжести этих сечений
...- средние скорости в этих сечениях
...- потери напора по длине.
Так как давление равномерное, то ... и уравнение можно
переписать так:
...
в случае равномерного движения разность удельных потенциальных
энергий равна потере напора по длине.
Для вычисления этой разности напишем сумму проекций на
ось А-А всех сил, действующих на участке 1-2. Эти силы
следующие:
1) сила тяжести жидкости
...
2) силы давления на плоские сечения
...
3) сила трения
...
где... - сила трения на единицу площади смачиваемой поверхности
русла,
... - смоченный периметр,
4) силы давления стенок русла на жидкость,эти силы не подсчи-
тываем, так как они параллельны оси А-А и, следовательно,
их проекции на ось А-А равны нулю.
Спроектируем все эти силы на ось А-А:
...
Из рисунка
...
Подставим выражение для сил в уравнение
...
Разделим обе части этого равенства на ..., имеем
...
Сравнивая выражения (1) и (2) , находим
...
откуда
...
Отношение площади живого сечения ... к смоченному периметру
... называется гидравлическим радиусом
...
Величина ... обозначается через ...
Получаем
...
Это уравнение называется основным уравнением равномерного
движения.
Величина ... имеет размерность квадрата скорости
...
Выражение ... - называется динамической скоростью,
обозначается ...
...
2.2. Два режима течения жидкости
Величина коэффициента трения зависит от режима течения
жидкости.
Опытами было установлено, что при течении жидкости возможны
два режима: ламинарный и турбулентный.
При ламинарном режиме жидкость течёт слоями, не перемеши-
ваясь.
При турбулентном частицы жидкости интенсивно перемешиваются.
Ламинарное и турбулентное течение жидкости можно наблюдать
в стеклянной трубе В.
Питание трубы производится из бака, а скорость течения
регулируется краном С. Для наблюдения за характером движения
жидкости по толстой трубе ... в трубу В подводится подкрашен-
ная жидкость такой же плотности, как и движущаяся жидкость
(например, чернило).
При малых скоростях в трубе В струйка продолжает двигаться,
не перемешиваясь с остальной жидкостью, что указывает на лами-
нарный режим течения.
При больших скоростях в трубе струйка очень сильно
перемешивается со всей жидкостью, что указывает на турбулентный
режим.
2.3. Критерии режима течения жидкости
В 1883 году английским учёным Осборном Рейнольдсом
(1842-1912 гг.) было установлено, что критерием режима течения
жидкости является безразмерная величина, представляющая собой
отношение произведения средней скорости потока и линейного
размера, характерного для живого сечения, к кинематической
вязкости жидкости ...
...
Критерий режим течения жидкости называется числом
Рейнольдса.
При ... течении жидкости в круглых трубах за характерный
размер ... объёма принимается внутренний диаметр трубы ...,
тогда
...
Пример.
Установить, какой режим будет в трубе диаметра ...=20 см,
если средняя скорость ... , а кинематическая вязкость ...
...
Опытные данные Рейнольдса показывают наличие трёх
областей:
АК - ламинарной,
ВК - переходной или
ВС - турбулентной.
Точки К и В называются критическими точками, точками, в
которых происходит смена режима течения.
Ниже точки К режим всегда ламинарный, выше точки В -
турбулентный.
В зависимости от изменения скорости от малых значений к
большим и от больших к малым ламинарный режим удерживается
до точки В при увеличении скорости, или при уменьшении до
точки К.
Значение числа Рейнольдса, соответствующее нижней крити-
ческой точке К, называется нижним критическим числом Рейнольдса,
число ... соответствует верхней критической точке - верхним
критическим числом Рейнольдса.
Нижнее число Рейнольдса ...= 956.
Переход к турбулентному режиму зависит (помимо скорости
течения, вязкости и характерного размера) от ряда факторов -
источников питания трубопровода, шероховатости труб, местных
сопротивлений и т.д.). Верхнее число Рейнольдса обычно принимают
равным ...= 5000.
На практике ламинарный режим встречается
1) при движении очень вязких жидкостей
2) при движении жидкости в ... трубах
3) при движении воды в грунтах.
Турбулентный режим наблюдается значительно чаще: при
движении в каналах, трубах и т.д.
Профиль скорости при ламинарном и турбулентном режиме
течения
При ламинарном режиме движения жидкостидвижение жидкости
как бы разделяется на бесконечно большое число тонких ... ...
относительно оси трубопровода слоёв.
Распределение скоростей по сечению имеет вид параболы.
Скорость у стены равна нулю. При удалении от стенки скорости
возрастают и достигают максимума на оси трубы.
Определим закон распределения скорости. Выделим объём жид-
кости в виде цилиндра радиуса ... и длиной ... и составим
уравнение равновесия
...
Движение установившееся, скорости на одном радиусе
одинаковы.
...
С учётом гидравлического закона
...
имеем
...
Проинтегрируем по сечению трубы, учитывая, что при ... ...
=0, получим закон распределения скоростей в сечении
...
Минимум скорости при ...=0
...
Определим расход жидкости через трубу
...
Средняя скорость
...
Соотношение между ...и средней скоростью
...
Турбулентный режим движения жидкости характеризуется ...
движением частиц. При этом режиме частицы жидкости движутся
по произвольным траекториям и с различной скоростью. Скорость
изменяется по величине и направлению около среднего значения.
Такое изменение скорости называется пульсацией скорости.
Среднюю по времени скорость называют осреднённой скоростью.
Связь между осреднённой и мгновенной скоростью может быть
выражена зависимостью
...
где Т - период наблюдения.
Распределение скоростей течения в этом случае выглядит
иначе, чем при ламинарном режиме.
В ламинарной пленке и переходном слое скорости течения
изменяются так же, как при ламинарном режиме течения.
В переходной зоне зарождаются вихри, обусловленные
увеличением скорости движения, влиянием выступов шероховатости.
Если выступы шероховатости ламинарной пленки, стенка будет
гидравлически гладкой. При величине выступов выше толщины
ламинарной пленки, неровности стенок будут увеличивать ...
движения и стенка будет гидравлически шероховатой.
Возникающие в пограничном слое вихри проникают в
центральную часть потока и образуют ядро турбулентного течения.
В ядре потока происходит интенсивное и непрерывное перемешивание
частиц жидкости.
Для описания профиля скорости в ядре течения турбулентного
состояния используется логарифмический закон распределения
скоростей
...
2.4. Определение потерь напора на трение по длине.
Формула Дарси-Вейсбаха
Теоретические и экспериментальные исследования показывают,
что потери напора на трение по длине вычисляются по формуле
Дарси-Вейсбаха
...
где ... - расстояние между рассматриваемыми сечениями, т.е.
длина трубы,
... - скорость течения,
... - внутренний диаметр трубы,
... - коэффициент гидравлических потерь на трение по
длине,
... - относительная шероховатость.
Для ламинарного режима движения жидкости
...
Основные два вопроса, которые интересуют инженера при
рассмотрении турбулентного движения жидкости в трубах:
1) определение потерь напора,
2) распределение скоростей по поперечному сечению трубы.
Потери напора и распределение скоростей могут сильно
меняться в зависимости от диаметра трубы, скорости движения,
вязкости жидкости и шероховатости стенок труб.
Для учета шероховатости используют понятие относительной
шероховатости
...
Cистематические опыты для выяснения характера зависимости
коэффициента гидравлического трения от числа ... Рейнольдса и
шероховатости ... были проведены H.Hикурадзе в 1933 в гладких
трубах с искусственной равномерно-зернистой шероховатостью из
кварцевого песка
от ... = 0.00197 до 0.066.
При различных расходах измерялась потеря напора и вычислялся
коэффициент ... по формуле Дарси-Вейсбаха.
Результаты опытов Hикурадзе представлены в виде графика
зависимости величины ... от числа ...
...
При ламинарном режиме ... 2000, или ... 3.3 ... точки,
независимо от шероховатости стенок, ложатся на прямую линию ...
При ламинарном режиме движения шероховатость не оказывает влияния
на сопротивление.
При турбулентном режиме (... 2000, ... 3.6) ... данные
ложатся на линию ..., полученную при испытании гладких труб без
искусственной шероховатости.Малые шероховатости не оказывают
влияния на сопротивление трубы при турбулентном движении.
При больших числах Рейнольдса коэффициент гидродинамического
трения перестает зависеть от числа Рейнолдса (то есть от вязкости
жидкости) и для данного значения ... сохраняет постоянную
величину.
Полученные результаты могут иметь следующее физическое
истолкование. При малых числах Рейнольдса жидкость обтекает
выступы шероховатости без образования и отрыва вихрей благодаря
значительному влиянию вязкости жидкости, свойства поверхности
стенок труб не оказывают при этом влияния на сопротивление и
кривые ... совпадают с прямой.
С увеличением скорости (т.е. числа Рейнольдса) от бугорков
шероховатости начинают отрываться вихри, свойства поверхности уже
оказывают влияние на сопротивление и кривые ... отклоняются от
линии ... трения.
В результате многочисленных исследований были предложены
различные эмпирические формулы для определения коэффициента
гидравлического трения.
Для гидравлически гладких труб широкое распространение
получили формулы Блазиуса
...
а для вполне шероховатых труб - формулы Шифринсона
...
3. Местные гидравлические сопротивления.
В гидросистемах часто встречаются повороты, краны, вентили,
сужения, расширения и т.д. В этих местах поток деформируется,
возникают интенсивные перемешивания жидкости, поперечные потоки,
образуются застойные зоны. Все это приводит к дополнительным
потерям напора, которые называются потерями напора на местных
сопротивлениях.
Рассмотрим гидросистему
...
1 - вход в трубу,
2 - внезапное расширение,
3 - ... сетка,
4 - внезапное сужение,
5 - диффузор,
6 - диафрагма,
7 - конфузор,
8 - поворот,
9 - тройник,
10 - колено,
11 - вентили, задвижки,
12 - поворот,
13 - вход в резервуар.
Потери напора, затраченные на преодоление местного
сопротивления, принято оценивать в долях скоростного напора,
соответствующего скорости непосредственно за рассмотренным
местным сопротивлением и определять по формуле Вейсбаха
...
... - коэффициент местного сопротивления.
Коэффициенты местных сопротивлений находят, обычно, опытным
путем. Таблицы и эмпирические формулы для них содержатся во всех
инженерных справочниках по гидравлике.
Для некоторых практически важных случаев значения
коэффициента местного сопротивления удалось получить
теоретически.
3.1. Внезапное расширение трубопровода
Рассмотрим потерю напора при внезапном расширении потока.
Пусть поток несжимаемой жидкости течет в горизонтальной трубе,
претерпевающей резкое увеличение площади поперечного сечения от
величины ... до ... .
Пусть скорость течения уменьшается при этом от ... до ... .
Массовый расход остается одинаковым в обоих сечениях
...
Секундное количество движения в сечении 1, ограничивающем
рассматриваемый элемент потока слева, равен
...
где ... - поправка к количеству движения на неравномерное
распределение скоростей в сечении.
Сечение 2, ограничивающее элемент потока справа, выбираем в
таком удалении от внезапного расширения, где возмущение течения,
вызванные в потоке расширением русла, можно полагать успокоенным.
В этом сечении секундное количество движения равно
...
Сила давления, действующая на выделенный элемент потока,
равна:
...
где ... , ... - давления в сечениях 1 и 2.
В проекции на ось трубы будет иметь следующее равенство
...
или
...
откуда
...
Уравнение Бернулли для двух сечений имеет следующий вид
...
или
...
Hа основании (...) имеем
...
Если положить ...=1, что верно для большинства турбулентных
потоков, то
...
Это положение, известное под название теоремы Борда,
формулируется так
Теорема Борда.
Потеря напора при внезапном расширении потока равна
скоростному напору, вычисленному по ... скорости.
...
Для других видов местных сопротивлений потеря напора
определяется по формуле, аналогичной внезапному расширению
...
Безразмерный коэффициент ..., входящий в формулу, называется
коэффициентом местного сопротивления.
Значение этого коэффициента зависит от конструкции местного
сопротивления, которая определяет характер отрыва потока от
обтекаемых внутренних полостей и интенсивность возникающих при
этом вихреобразований.
Часто при определении потерь напора на местные сопротивления
оказывается удобным введение так называемой эквивалентной длины
детали трубопровода.
Эквивалентной длиной данного местного сопротивления называют
такую длину прямого отрезка трубы, которая создает гидравлическое
сопротивление, равное сопротивлению детали трубопровода,
обусловившей потери напора.
Пусть ... - эквивалентная длина данного местного
сопротивления, потеря напора на прямом участке трубы длиной ...
по формуле равна
...
По условию эквивалентности должно быть ..., откуда ...,
следовательно
...
Таким образом, эквивалентная длина местного сопротивления
выражается через диаметр трубы, поэтому, например, говорят, что
сопротивление углового вентиля эквивалентно сопротивлению участка
трубы того же диаметра длиной, равной 200 диаметрам трубы.
Пусть требуется определить потерю напора в трубопроводе,
состоящем из прямых отрезков труб, соединенных между собой с
помощью всевозможных ... частей, с включением различного рода
задвижек, вентилей, клапанов и т.д. Эту задачу можно решить,
определяя по формулам и таблицам из справочников, коэффициенты
местных сопротивлений ... или вычислив предварительно
эквивалентные длины местных сопротивлений.
В первом случае потеря напора может быть определена по
формуле
...
а во втором - по формуле
...
Исследованию местных коэффициентов сопротивлений посвящается
обширная литература, проделано огромное количество опытов, однако
до сих пор задача о местных сопротивлениях остается разрешенной
еще не полностью.
Можно считать доказанным, что величина местного
сопротивления при ламинарном течении меняется в зависимости от
числа ..., при турбулентном режиме она остается почти постоянной
при любых ...
...
4. Гидравлический расчет напорных трубопроводов.
4.1. Классификация трубопроводов. Задача гидравлического
расчета трубопроводов.
Трубопроводы широко применяются для перемещения различных
жидкостей (вода, нефть, бензин, различные растворы и т.д.) и
изготавливаются из металла, бетона, дерева, пластмасс.
По степени заполнения поперечного сечения жидкостью
различают напорные и безнапорные трубопроводы. В напорных
трубопроводах жидкостью заполнено полностью все поперечное
сечение, а в безнапорных - часть поперечного сечения и имеется
свободная поверхность.
По виду потерь напора бывают короткие и длинные
трубопроводы.
Короткие трубопроводы - это такие трубопроводы, у которых
местные потери напора соизмеримы с потерями напора по длине.
Длинные трубопроводы - это трубопроводы, у которых местные
потери напора незначительны и не превышают 10% от потерь по
длине.
В свою очередь, длинные трубопроводы разделяются на простые
и сложные.
Простые трубопроводы выполняют без ответвлений, сложные
изготавливают с отверстиями, переменной длины и диаметра и могут
соединяться как последовательно, так и параллельно.
Задача гидравлического расчета трубопровода заключается в
определении для заданной длины по двум величинам третьей
неизвестной величины: расхода жидкости ..., потери напора ...,
диаметра трубопровода ...
4.2. Расчет коротких трубопроводов.
Рассмотрим короткий трубопровод с местным сопротивлением,
присоединенным к резервуару, заполненному жидкостью. Истечение
жидкости в атмосферу из трубопровода длиной ... и диаметром ...
происходит под постоянным напором H.
При заданных длине и диаметре трубопровода ... необходимо
определить скорость движения жидкости ... и расход ... .
Составим уравнение Бернулли для сечений 1 и 2. При этом
считаем, что ... и ...
...
или
...
где ... - суммарные (местные и по длине) потери напора между
сечениями 1 и 2, которые можно представить в виде зависимости
...
где
...
Формулу можно записать в следующем виде
...
Отсюда найдем скорость истечения
...
где ... - коэффициент скорости.
Расход, пропускаемый коротким трубопроводом
...
4.3. Расчет длинных трубопроводов при
последовательном соединении труб.
Рассмотрим трубопровод, состоящий из последовательных
длинных труб разного диаметра ... и длины ... при постоянном
расходе жидкости по длине трубопровода.
Расчет сводится к определению суммарных потерь напора по
длине трубопровода, так как местными потерями пренебрегают.
...
Преобразуем выражение для потери напора по длине
...
где ... - расходная характеристика.
Тогда
...
Формула показывает, что трубопровод, составленный из
последовательно соединенных труб разного диаметра и длины, можно
рассматривать как простой трубопровод, суммарные потери напора в
котором равны сумме потерь напора составляющих его труб.
Формула позволяет решить и обратную задачу, т.е. при
заданных напоре, диаметре труб вычислить расход ...
...
4.4. Расчет трубопровода при параллельном соединении труб.
Особенность гидравлической схемы работы трубопровода при
параллельном соединении труб состоит в том, что все трубы
работают под действием напора ..., который необходим для
преодоления потерь напора по длине ... При этом следует иметь в
виду, что во всех ответвлениях параллельных труб потери напора
будут одинаковыми.
...
Расчет трубопровода при параллельном соединении труб
сводится к составлению для каждого ответвления уравнения
...
и общего уравнения для расхода жидкости в трубопроводе
Тема 9
Неустановившееся движение жидкости
1. Гидравлический удар в трубах.
2. Вытекание жидкости при переменном уровне.
Неустановившимся движением жидкости называется такое движе-
ние, при котором скорости в точках пространства, занятого жидкос-
тью, изменяются со временем.
С неустановившимся движением воды часто сталкиваются при про-
ектрировании трубопроводов, расчёте каналов, водопроводных сетей,
истечении жидкости при переменном уровне, расчёте гидравлического
удара в трубах.
В рассматриваемом курсе для примера выполним исследование гид-
равлического удара в трубах и истечения жидкостей при переменном
уровне.
1. Гидравлический удар в трубах
Если при напорном движении жидкости в трубе мгновенно закрыть
кран, то движущаяся жидкость остановится, кинетическая энергия по-
тока израсходуется на сжатие жидкости и расширение стенок трубы.
Вследствие сжатия жидкости и расширения стенок трубы любое
сечение А-А, взятое в жидкости, сместится по направлению движения
в положение В-В.
...
...
Аналогичные явления произойдут и со всеми остальными сечения-
ми. Таким образом, вся жидкость в трубе по окончанию деформации
окажется сжатой, а поэтому обладающей большей энергией, чем жид-
кость в баке.
В результате этого начинается обратное движение жидкости и
сечение В-В, пройдя своё первоначальное положение А-А, займёт мес-
то С-С.
Аналогичное движение совершают и все остальные сечения,
вследствие чего в трубе создаётся пониженное давление и жидкость
двинется от сосуда к крану. Затем все явление повторяется и будет
повторяться снова, пока под влиянием сопротивления оно постепенно
не прекратится.
Частицы жидкости будут совершать затухающие колебания, одно-
временно с которыми будет изменяться и давление. Изменение давле-
ния в жидкости при напорном движении, вызываемое резким изменением
скорости течения за весьма малый промежуток времени, называется
гидравлическим ударом.
Увеличение давления при гидравлическом ударе может привести к
разрыву стенок трубы. Это увеличение давления в первый момент про-
исходит непосредственно к крана, а затем оно передаётся через со-
седние слои по всей длине l трубы до её начала с некоторой скорос-
тью С. Эта скорость носит название скорости распространения удар-
ной волны.
По истечении времени ... ударная волна дойдёт до нача-
ла трубы и вся жидкость в трубе остановится.
Определим величину повышения давления в трубе при гидравли-
ческом ударе.
Пусть давление в горизонтальной трубе в сечении 1 равно ...,
а в сечении 2 - ..., площадь поперечного сечения трубы ..., рассто-
яние между сечениями 1-1 и 2-2 - ... .
Воспользуемся теоремой об изменении количества движения, со-
гласно которой приращение количества движения системы за некоторый
промежуток времени равно сумме проекций импульсов сил на направле-
ние движения.
Применим теорему к массе жидкости, заключённой между сечения-
ми 1-1 и 2-2. В момент закрытия крана количество движения жидкости
равнялось ... , где ... - масса жидкости, равная ... , ... -
скорость. Через промежуток времени ... , т. е. когда вся
жидкость в трубе остановится и скорость будет равна нулю, коли-
чество движения также будет равно нулю.
Следовательно, за время ... приращение количества дви-
жения равно ...
В течении этого времени на жидкость действовали следующие си-
лы, не считая сил трения, которыми пренебрегаем:
1) в сечении 1-1 сила ...
2) в сечении 2-2 сила ...
3) сила тяжести жидкости ...
Первые две силы горизонтальны, третья вертикальна.
Сума проекций импульсов этих сил на направление движения, т.е.
на горизонтальную ось равна
...
Согласно теореме об изменении количества движения получаем
...
Сокращая на ... , имеем
...
Откуда
...
Обозначив повышение давления ...-... буквой Р, находим
...
Формула называется формулой Н.Е.Жуковского, который первый
дал теорию гидравлического удара.
Разделим последнее соотношение на ..., получим
... или ...
Из формулы видно, что при гидравлическом ударе повышение напо-
ра в трубопроводе равно ... .
Численное значение величины С также выведено Н.Е.Жуковским и
определяется по следующей формуле
...
где ... - плотность жидкости, .. - модуль упругости жидкости,
... - модуль упругости стенок трубы, ... - внутренний диаметр
трубы, ... - толщина стенки трубы.
Рассмотрим пример.
Пример. Определить повышение напора при гидравлическом ударе
в чугунной трубе диаметром ... , если толщина
стенки
трубы ... =
... = 2 ... , модуль упругости чугуна
..... , а скорость течения ... .
Решение.
По формуле находим скорость распространения ударной волны
...
Найдём повышение напора
...
Гидравлический удар может повредить трубы. Для предотвращения
разрушения труб применяются следующие меры.
1. Из формулы ... видно, что увеличение давления пропор-
ционально скорости течения ..., поэтому в трубопроводах не следует
допускать больших скоростей без принятия соответствующих предохра-
нительных мер.
2. Причиной гидравлического удара является быстрое закрытие
крана. При продолжительности закрытия .......... повышение давле-
ния равно ... (так называемый прямой гидравлический удар). При
продолжительности ......... повышение давления меньше ......
(непрямой гидравлический удар).
Продолжительность закрытия ... (в секундах) может быть под-
считана по формуле Н.Е.Жуковского
..... или ...................,
где ... - плотность жидкости, ... - скорость течения, ... -
длина трубопровода, ... - допустимое повышение напора столба жид-
кости (в метрах).
Время закрытия трубопровода ... прямо пропорционально длине
трубопровода ... . Т.е. чем длинее трубопровод, тем длительнее
должно быть закрытие кранов и задвижек.
3. Для уменьшения вредного действия давления при гидравличес-
ком ударе ставят предохранительные клапаны, которые, открываясь
при определённом давлении, предохраняют провод от разрушения.
4. Кроме предохранительных клапанов, для уменьшения давления
применяют воздушные колпаки. В момент повышения давления жидкость
входит в колпак и сжимает находящийся в нём воздух, что уменьшает
повышение давления.
Пример. Определить продолжительность закрытия задвижки на
трубопроводе, если длина трубопровода ... =
... = 3 ..., допускаемое давление в трубопроводе
1 000 000 ..., а гидростатическое давление
Р = 200 000 ... .
Решение.
Допускаемое повышение давление от гидростатического удара
... = 1 000 000 - 200 000 = 800 000 ...
Продолжительность закрытия задвижки
...
2. Вытекание жидкости при переменном уровне
Рассмотрим случай истечения жидкости из открытого сосуда в
атмосферу через отверстие площадью ... .
Струя при вытекании через отверстие постепенно сжимается.
Ближайшее к отверстию наименьшее живое сечение С-С, в котором дви-
жение можно рассматривать плавно изменяющимся, называется сжатым
сечением. Обозначим площадь сжатого сечения С-С ...
...
...
Отношение
...
( ... = 0.64 для круглого отверстия)
называется коэффициентом сжатия.
Обозначим через ... высоту уровня жидкости над центром тяжес-
ти отверстия, ... - скорость в сжатом сечении.
Запишем уравнение Бернулли для сечений О-О и сжатого сечения
С-С.
...
где ... - скорость свободной поверхности,
... - потери напора при вытекании через отверстие, они
определяются из соотношения
...
Пренебрегая величиной ... (ввиду её малости по сравнению
с Н), получаем
...
отсюда скорость истечения
...
где ... - коэффициент скорости (.....0.97).
Для определения расхода надо скорость умножить на площадь
сжатого сечения:
...
по формуле ... , откуда
...
тогда расход, выраженный через ... равен
...
где ... - коэффициент расхода (... = 0.62).
Рассмотрим вытекание жидкости из ёмкости при переменном уров-
не. Движение в данном случае является неустановившимся. С доста-
точной для практики точностью можно считать, что в каждый момент
времени скорость вытекания определяется соответствующим этому мо-
менту напором Н так же, как и при установившемся движении.
...
...
Определим время, в течение которого жидкость опустится на
...-...
Рассмотрим промежуточное положение уровня с напором Н. За
время ... вытечет объём жидкости, равный
...
За это время ... напор изменится на (-...Н). Объём жидкости,
вытекшей из сосуда, равен
...
где ... - площадь свободной поверхности в сосуде.
Приравнивая выражения, получаем
...
откуда
...
Интегрируя, находим
...
При постоянной площади свободной поверхности
...
Пример. Вычислить продолжительность опорожнения цистерны при
её диаметре ... =
сливного отверстия ... =
да ... = 0.62.
Решение. Продолжительность опорожнения
...
...
... - переменная по высоте горизонтальная площадь сечения ци-
стерны, причём
...
Тема 10
Кинематика плоских движений жидкости
1. Сетка течения плоского потока несжимаемой жидкости.
Функция тока.
2. Примеры плоских течений.
1. Однородный равномерный поток.
2. Источник и сток.
3. Вихрь.
4. Вихреисточник.
5. Диполь.
3. Бесциркуляционное обтекание цилиндра.
1. Сетка течения плоского потока несжимаемой жидкости.
Функция тока
В гидродинамике невязкой жидкости особенно полно разработана
теория плоских стационарных (установившихся) течений.
Пусть, например, плоский безграничный поток обтекает цилин-
дрическое (или призматическое) тело, бесконечное в направлении,
перпендикулярном к скорости течения. Характер течения (обтекания)
тела будет одинаков во всех плоскостях, перпендикулярных к обра-
зующим тела.
Следовательно, для исследования кинематики и динамики такого
потока достаточно рассмотреть плоскую задачу "обтекаемого" тела.
В этом случае скорости и давления зависят только от двух коорди-
нат, пусть, например, X и Y, также функциями этих двух координат
являются проекции ... и ... скорости течения.
Пусть определена функция ... , которая удовлет-
воряет следующим условиям
...
Такая функция называется в гидромеханике функцией тока.
Уравнение линий тока в случае плоского течения имеет вид:
...
или
...
Подставляя сюда выражения проекций скорости через частные
производные функции ..., найдём
...
При установившемся течении левая часть этого выражения пред-
ставляет собой полный дифференциал функции ..., напишем
...
Отсюда следует, что ... , таким образом, функция
тока на линии тока сохраняет постоянное значение.
Предположим, что рассматриваемый плоский поток является по-
тенциальным, т.е. что во всех точках потока имеет место условие
...
В соответствии с принятыми предположениями в этом случае
...
где ... - потенциал скорости.
Из условия ... имеем
...
Подставляя сюда выражение для функции тока, получим
...
Поскольку мы рассматриваем несжимаемую жидкость, то уравнение
неразрывности принимает вид
...
или через потенциал скорости
...
Дифференциальные уравнения второго порядка, выражающее, что
сумма вторых частных производных скалярной функции равняется нулю,
являются, как известно, уравнениями Лапласа.
Таким образом, потенциал скорости и функция тока удовлетворя-
ют уравнению Лапласа.
Это уравнение обладает следующим свойством. Если имеются
функции, например, ..., ..., ...или ..., ..., ...такие, что
каждая из них в отдельности удовлетворяет уравнению Лапласа, то
ему будут удовлетворять также их линейные комбинации
...
...
где ..., ..., ..., ..., ...- постоянные.
Отсюда следует, что при наложении одного плоского потенциаль-
ного потока на другой потенциальный поток полученное движение бу-
дет также потенциальным и его потенциал скорости и функция тока
будут определяться путём суммирования значений потенциалов и функ-
ций тока слагаемых потоков.
Если построить два семейства кривых: кривые ... = К,
представляющие собой эквипотенциальные линии (т.е. линии равного
потенциала) и кривые ... = ... линии тока (здесь К и ... -
параметры), то эти семейства кривых образуют ортогональную сетку
плоского течения.
...
...
Это можно показать следующим образом. Вектор скорости ...,
совпадающий с направлением касательной к линии тока, образует с
осью абсцисс угол ..., тангенс которого с учётом выражения для
скоростей равен
...
Из уравнения же эквипотенциальной линии следует
...
и отсюда тангенс угла ..., который образует касательная к
эквипотенциальной линии с осью абсцисс, равен
...
Показать, что векторы ... взаимно перпендикулярны,
можно так
...
...
В результате перемножения получаем
...
Этому условию отвечают условные коэффициенты взаимно перпен-
дикулярных линий.
Функция тока ... имеет физический смысл. Определим расход
жидкости через сечение потока между двумя линиями тока ... и ...
(т.е. расход струйки тока, ограниченной поверхностями, для которых
названные линии тока являются образующими), размер сечения струйки
по нормали к плоскости ...будем предполагать равным единице.
...
где ... - элемент живого сечения струйки, ... - ...,
... - единичный вектор по нормали к элементу ... ,
... и ... - границы сечения.
Обозначим через ... угол, образуемый вектором ... с осью ...,
тогда ... и ... будут проекциями этого вектора на оси координат и,
следовательно,
...
но ...
поэтому
...
...
Таким образом, разность значений функции тока на двух каких-
нибудь линиях тока равна секундному объёмному расходу сквозь сече-
ние струйки тока, ограниченной соответствующими поверхностями то-
ка.
Из сопоставления
...
следует
...
Из теории функций комплексного переменного следует, что если
выполняются условия Коши-Римана, то линейная комбинация
...
функций ... и ... является функцией комплексного переменного
..., т.е.
...
Функция ...называется комплексным потенциалом, последний
удовлетворяет уравнению Лапласа.
Найдём производную от комплексного потенциала
...
причём
...
...
где ... и ... - бесконечно малые величины высшего порядка. В
пределе
...
Из этого выражения с учётом условий Коши-Римана следует
...
- это выражение называется комплексной скоростью.
Модуль комплексной скорости даёт величину скорости
...
Вводим комплексную скорость
...
сопряжённую скорость
...
Тогда
...
...
...
...
Рассмотрим
...
Тогда
... - циркуляция
... - расход.
2. Примеры плоских течений
1. Однородный равномерный поток.
Рассмотрим плоское прямолинейное и равномерное установивше-
еся течение несжимаемой жидкости с одинаковой скоростью во всём
потоке скоростью ... , параллельной оси ... . В этом случае
...
Отсюда
...
Линии равных потенциалов ... представляют собой пря-
мые, параллельные оси ординат.
Можно положить ... = 0 и ... = 0, тогда
...
Функцию тока найдём из условия
...
Сетка такого плоского течения изображается семейством ортого-
нальных прямых, параллельных осям координат, а комплексный потен-
циал равен
...
Для прямолинейного течения сжимаемой невязкой жидкости со
скоростью ..., наклонённой к оси абсцисс под углом ..., будем
иметь
...
откуда
...
и
...
Комплексный потенциал такого течения будет иметь вид
...
...
2. Источник и сток
В качестве следующего примера рассмотрим течения, которые
носят название источника и стока.
...
...
Пусть невязкая несжимаемая жидкость непрерывно возникает в
некоторой точке Р и вытекает в неограниченное пространство с по-
стоянным расходом ... и с одинаковой интенсивностью во всех на-
правлениях.
Линии тока этого воображаемого источника будут представлять
собой прямые, расходящиеся из точки Р. Это характеризует простран-
ственный источник.
Если жидкость течёт из неограниченного пространства в точку,
где непрерывно исчезает, течение называется пространственным сто-
ком.
Рассмотрим плоский источник и проведём из него как из центра
несколько концентрических окружностей различного радиуса. Уравне-
ние неразрывности - уравнение постоянства расхода через любую кон-
центрическую цилиндрическую поверхность, имеющую высоту, равную
единице, в случае несжимаемой жидкости будем считать
...
Отсюда скорость
...
и, следовательно,
...
Откуда
...
Интегрируя
...
где С -константа интегрирования, которая может быть принята
равной нулю, если полагать, что на круге ... = 1 функция ... = 0.
Для определения функции тока воспользуемся выражением
...
откуда полный дифференциал
...
После интегрирования имеем
...
...и С = 0 при ... = 0.
Следовательно
...
Потенциал скорости источника ...(...)может быть интерпрети-
рован в виде семейства концентрических кругов различного радиуса,
а функция тока ...(...) в виде пучка прямых, исходящих из источни-
ка.
3. Вихрь
Рассмотрим комплексный потенциал
...
Пусть А - действительное число
...
...
...
Линии тока лучи ...
Изопотенциальные линии - окружности.
Найдём расход
...
...
...
...
... - комплексный потенциал источника или стока мощнос-
ти ...
Пусть А - чисто мнимое.В..., где В - действительное.
...
...
4. Вихреисточник
Рассмотрим случай комплексного коэффициента при логарифме
...
Такой комплексный потенциал можно рассматривать как результат
наложения двух потоков
...
...
... - комплексный потенциал вихреисточника.
5. Диполь
Рассмотрим комплексный потенциал ...
...
...
Найдём семейство линий тока
...
...
Линии тока - окружности с центрами на оси ...
Изопотенциальные линии - окружности с центрами на оси ...
Диполь
...
где ... - момент диполя.
3. Бесциркуляционное обтекание цилиндра.
Наложим плоский параллельный оси ...однородный поток со
скоростью .... и комплексным потенциалом
...
на скоростное поле диполя с комплексным потенциалом
...
...
...
... функции тока отделим ... часть
...
Нулевая линия тока
...
Решение распадается на две кривые
1) окружность...
2) ось ... ... = 0.
Выберем произвольную до сих пор величину момента диполя рав-
ной
...
Получим нулевую линию тока в виде совокупности окружности
радиуса а с центром в начале координат и оси ... .
Остальные линии тока
...
Движение происходит в двух областях - вне и внутри круга.
Течение вне круга можем рассматривать как обтекание круглого
цилиндра а плоскопараллельным потоком, имеющим на бесконечности ско-
рость ...
Такому потоку соответствует комплексный потенциал
...
Остановимся подробнее на внешнем течении. Найдём распределе-
ние скоростей в области ...
Найдём распределение скоростей на поверхности цилиндра
...
...
...
Найдём модуль скорости на контуре круга
...
Отсюда следует, что при плоском безвихревом обтекании круго-
вого цилиндра идеальной жидкостью скорость распределена по закону
синуса.
Максимальная скорость при ...
...
Используя уравнение Бернулли, можно найти распределение да-
вления
...
...
... - коэффициент давления
...
...
Циркуляционное обтекание цилиндра
...
Тема 11
УРАВНЕНИЯ НАВЬЕ-СТОКСА
1. Тензорная запись уравнений Эйлера.
Тензор плотности потока импульса.
2. Тензор плотности потока импульса для вязких течений.
3. Уравнения Навье-Стокса в декартовых координатах.
4. Течение в трубе.
1. Тензорная запись уравнений Эйлера.
Тензор плотности потока импульса.
Определим скорость изменения импульса единицы объема
жидкости
...
Воспользуемся тензорными обозначениями
...
Из уравнения неразрывности имеем
...
Воспользуемся уравнениями Эйлера, записанными в тензорной
форме
...
Таким образом получаем
...
Член с давлением запишем в виде
...
Уравнения количества движения принимают вид
...
где тензор ...определяется как
...
Выясним физический смысл тензора ... . Проинтегрируем
уравнение количества движения по некоторому объему
...
Преобразуем интеграл в правой части в интеграл по
поверхности
...
Слева стоит изменение в единицу времени i - той компоненты
импульса в рассматриваемом объеме. Поэтому интеграл по
поверхности в правой части есть количество импульса, вытекающего
в единицу времени через ограничивающую объем поверхность.
Следовательно, ... есть i - я компонента импульса,
протекающая через элемент ... поверхности.
Тензор ... называют тензором плотности потока импульса.
2. Тензор плотности потока импульса для вязких течений
Плотность потока импульса, определяемая соотношением
...
представляет собой обратимый процесс переноса импульса,
связанный с механическим передвижением различных участков
жидкости из одного места в другое и с действующими в жидкости
силами давления.
Вязкость ( внутреннее трение ) жидкости проявляется в
наличии еще дополнительного, необратимого, переноса импульса из
мест с большей в места с меньшей скоростью.
Поэтому уравнения движения вязкой жидкости можно получить,
прибавив к "идеальному" потоку импульса дополнительный член ... ,
определяемый необратимый, "вязкий" перенос импульса в жидкости.
Таким образом, мы будем писать тензор плотности потока
импульса в вязкой жидкости в виде
...
Тензор
...
называют тензором напряжений, а ... - вязким тензором
наряжений.
... определяет ту часть потока импульса, которая не связана
с непосредственным переносом импульса вместе с массой
передвигающейся жидкости.
Процессы внутреннего трения в жидкости возникают только в
тех случаях, когда различные участки жидкости движутся с
различной скоростью, так что имеет место движение частей жидкости
друг относительно друга.
Поэтому ... должно зависеть от производных скорости по
координатам. Если градиенты скорости по координатам не очень
велики, то можно считать, что обусловленный вязкостью перенос
импульса зависит только от первых производных скорости.
Зависимость ... от производных ... можно в том же
приближении считать линейной. Не зависящие от ... члены должны
отсутствовать в выражении для ... , поскольку ... должно
обращаться в нуль при ... = const.
... должно обращаться в нуль также в том случае, когда вся
жидкость как целое совершает равномерное вращение, поскольку при
таком движении внутреннее трение не происходит. При равномерном
вращении с угловой скоростью ... скорость ... равна векторному
произведению .... Линейными комбинациями производных ... ,
обращающимися в нуль при ... , являются суммы
...
Поэтому ... должно содержать именно эти симметричные
комбинации производных ... .
Наиболее общим видом тензора второго ранга, удовлетворяющего
этим требованиям, является
...
с не зависящими от скорости коэффициентами ... и ... .
Величины ... и ... называются коэффициентами вязкости ( причем ..
часто называют второй вязкостью ).
3. Уравнения Навье-Стокса в декартовых координатах
Уравнения движения вязкой жидкости можно
теперь получить непосредственно путем прибавления выражения ... к
правой части уравнений Эйлера
...
Получаем,
...
Величины ... и ... являются в общем случае функциями
давления и температуры. Поэтому они не постоянные в объеме и не
могут быть вынесены из-под знака производной.
При постоянных значениях коэффициентов вязкости уравнения
Навье-Стокса в векторной форме имеют вид
...
Уравнения были впервые сформулированы Навье в 1827 году,
вывод уравнений близкий к современному, был дан Стоксом в 1845
году.
Если жидкость считать несжимаемой, то ... = 0 и последний
член исчезает
...
Тензор напряжений в несжимаемой жидкости принимает более
простой вид
...
Отношение... = ... называют кинематической вязкостью, ...
- динамической вязкостью.
Граничные условия.
Между поверхностью твердого тела и вязкой жидкостью
существуют силы межмолекулярного сцепления, приводящие к тому, что
прилегающие к твердой стенке слой жидкостью как бы прилипает к
ней.
Граничные условия к уравнениям движения вязкой жидкости
состоит в требовании обращения в нуль скорости жидкости на
неподвижных твердых поверхностях
...
В общем случае движущейся поверхности скорость ... должна
быть равна скорости этой поверхности.
4. Течение в трубе
Известно несколько точных решений для уравнений
Навье-Стокса. Рассмотрим одно из них - для случая стационарного
течения жидкости в трубе произвольного сечения ( одинакового
вдоль всей длины трубы ).
Ось трубы выберем в качестве оси... . Очевидно, что
скорость ... жидкости направлена везде по оси ... и является
функцией только от ... и ... .
Уравнение неразрывности удовлетворяется тождественно, а
проекции на оси ... и ... из системы уравнений Навье-Стокса дают
...
То есть давление постоянно вдоль сечения трубы. Уравнение в
проекции на ось ... дает
...
Откуда имеем, что ... = const , градиент давления можно
записать в виде ... , где ... - разность давлений на концах
трубы, а ... - ее длина.
Распределение скоростей в потоке жидкости в трубе
определяется двумерным уравнением типа
...
Уравнение должно быть решено при граничном условии ...= 0 на
контуре сечения трубы.
Решим это уравнение для трубы кругового сечения. Выбирая
начало координат в центре трубы кругового сечения и вводя
полярные координаты, имеем в силу симметрии ... .
Воспользуемся выражением для оператора Лапласа в полярных
координатах, имеем
...
Интегрируя, находим
...
Постояннуюaнадо положить равной нулю, поскольку скорость
должна оставаться конечной во всем сечении трубы, включая ее
центр.
Постояннуюbопределим из требования... = 0, при r = R (
где R - радиус трубы ) и получаем
...
Таким образом, скорость распределена по сечению трубы по
параболическому закону.
Определим расход жидкости в трубе - количество ( массу )
жидкостиQ, протекающей в 1 секунду, через поперечное сечение
трубы.
Через кольцевой элемент ... площади сечения трубы проходит
в 1 секунду количество жидкости ... .
Поэтому
...
Количество протекающей жидкости пропорционально четвертой
Тема 12.
Дозвуковое и сверхзвуковое течения газов
(основы газодинамики)
1. Адиабатически установившееся течение.
2. Уравнение состояния.
3. Удельные теплоемкости газа.
4. Первый закон термодинамики. Энтальпия. Энтропия.
5. Характеристики заторможенного потока.
6. Сопло Лаваля.
7. Скачок уплотнения.
8. Теория Ньютона.
1. Адиабатическое установившееся течение. Истечение из
резервуара. Характеристики заторможенного газа
Изучение движения газов с высокими скоростями, достигающими
скорости звука, является предметом газовой динамики.
Одной из фундаментальных задач последней является исследова-
ние течений без учёта сопротивлений и в отсутствие теплообмена
(т.е.) адиабатических. В этих условиях уравнение баланса удельной
энергии имеет вид
...
Уравнение адиабаты идеального газа представим в виде
...
Будем отмечать в дальнейшем индексом "о" величины, характери-
зующие газ, находящийся в покое, или, как говорят в газодинамике,
в заторможенном состоянии, подставим в уравнение неразрывности
...
и после интегрирования
...
...
При установившемся течении весовой расход газа во всех сече-
ниях по длине газопровода одинаков в течение всего процесса движе-
ния.
Следовательно при установившемся течении
...
что является выражением условия неразрывности при движении
газа (и также сжимаемых жидкостей). В трубопроводе постоянного се-
чения одинаковой по длине трубопровода будет также весовая ско-
рость
...
Изменение в удельном весе (плотности) идеального газа при из-
менении давления и температуры выражаются законом Клайперона-Мен-
делеева
...
где Т - абсолютная температура газа,
... - газовая постоянная.
В технике имеют особое значение изотермическое и адиабатичес-
кое течения газа. При изотермическом (Т = ...) течении идеально-
го газа зависимость между давлением и плотностью получает вид
...
при адиабатическом
...
где ... - показатель адиабаты, ... - удельная теплоёмкость
газа при постоянном давлении, ... - удельная теплоёмкость газа
при постоянном объёме.
Имея в виду последнее соотношение, можно записать
...
получаем
...
Имея в виду, что ... = 0 при ... (состояние покоя),
найдём:
...
или
...
Уравнение Гюгонио. Сопло Лаваля
Запишем уравнение Бернулли в дифференциальной форме
...
Преобразуем уравнение Бернулли для газа так, чтобы можно было
ввести число Маха. Имеем
...
квадрат скорости звука ... , тогда
...
Поделим на ... , получим
...
или в окончательном виде
...
где ..... - число Маха.
Другим уравнением, необходимым для анализа течений газа в
трубе переменного сечения, является уравнение неразрывности,
или сохранения массы.
Будем рассматривать одномерное установившееся течение
газа вдоль трубы переменного сечения, при этом предположим, что
параметры потока газа, такие, как скорость потока, давление и
плотность, одинаковы во всех точках каждого из конечных сечений,
перпендикулярных к оси трубы.
Это предположение довольно хорошо соответствует действи-
тельности для элементарной трубки тока, но его применяют и для
труб конечных размеров, используя средние величины по сечениям
трубы.
Через каждое поперечное сечение трубы в случае одномерного
течения проходит за 1 с масса газа ..........., где ... -
площадь поперечного сечения трубы, ... - скорость течения газа,
... - плотность газа. ПРи установившемся течении через все по-
перечные сечения должна пройти одна и та же масса газа, т.е.
...
Прологарифмируем это уравнение сохранения массы. Получим
...
Считая переменными величины .............., возьмём полные
дифференциалы от обеих частей. Имеем
...
Это и есть уравнение неразрывности для установившегося одно-
мерного течения идеального газа в трубе переменного сечения.
Из уравнения неразрывности и уравнения Бернулли исключим
величину ... . Получим
...
Это уравнение носит название уравнения Гюгонио.
Используя уравнение Гюгонио, проанализируем характер воз-
можных течений газа в трубе переменного сечения.
Из уравнений следует:
1) при ... 1, что соответствует дозвуковым течениям, знаки
величин ... и ... противоположны, т.е. там, где возрастает
..., в направлении течения скорость должна убывать, и на-
оборот
2) для сверхзвуковых течений ......1, знаки ... и ... оди-
наковы, т.е. сверхзвуковой поток расширяется противополож-
но дозвуковому. Чтобы увеличить его скорость, трубу следу-
ет расширить
3) при ... = 1 имеем ... = 0, т.е. в этом случае ... достига-
ет максимума или минимума. Можно показать, что ... = 1
может быть только в самом узком сечении трубы, где .......
Выводы о характере течений газа в трубах переменного сечения
нашли применение в конструкциях сопел современных ракетных двига-
телей и аэродинамических трубах больших скоростей. Для получения
больших сверхзвуковых скоростей выходящего из сопла газа следует
сначала сопло сужать, чтобы получить звуковую скорость газа в уз-
ком сечении сопла, а затем сопло надо расширять для дальнейшего
увеличения скорости выходящего из него газа.
Наибольшая скорость, которая может быть получена на выходе из
сопла, зависит от площади выходного сечения и должна обеспечивать-
ся необходимым для ... скорости давлением на входе в сопло.
1. Уравнение состояния.
Опыт показывает, что между основными параметрами, характери-
зующими состояние газа (давлением, плотностью и температурой), су-
ществует определённая зависимость.
Уравнение ............ = 0 , устанавливающее связь между
этими параметрами, называется уравнением состояния.
Поэтому состояние любого газа определяется двумя параметрами
(например, плотностью и температурой), так как третий параметр
(давление) можно найти из уравнения состояния.
Для идеального газа уравнение состояния можно записать в виде
...
где ... - газовая постоянная, зависящая от относительной
молекулярной массы газа ... . Для воздуха ... = 29, ... = 287 ...
Под идеальным газом принято понимать газ, в котором взаимо-
действие молекул между собой осуществляется посредством упругих
столкновений, а линейный размер молекулы по сравнению со средним
молекулярным расстоянием мал.
Существенное отличие свойств воздуха от свойств идеального
газа наблюдается при высоких давлениях и низких температурах.
2. Уравнение теплоёмкости газа.
Рассмотрим некоторый произвольный термодинамический процесс.
Количество теплоты ..., подведенное к
выразим через приращение температуры газа ... :
...
Множитель С, представляющий собой количество теплоты, необхо-
димое для подогрева
ется удельной теплоёмкостью.
Удельная теплоёмкость существенно зависит от характера про-
цесса.
Рассмотрим теплоёмкости, соответствующие процессам, происхо-
дящим при постоянном объёме ... и давлении ... . Зависимость между
удельными теплоёмкостями идеального газа ... и ... определяется
следующим соотношением.
...
В термодинамике и газодинамике важное значение имеет отноше-
ние теплоёмкостей ...... Величина ... зависит от структуры молеку-
лы газа. Так, для идеальных одноатомных газов ... = 1.66, для
двухатомных газов, в том числе и для воздуха, ... = 1.4.
3. Первый закон термодинамики.
Пусть некоторое количество газа находится в равновесии.
Обозначим через ... количество подведённой к газу извне теплоты.
В общем случае подвод теплоты приводит к изменению внутренней
энергии газа ... и объёма. ПРи изменении объёма газ совершает
внешнюю работу, равную ... . Поэтому
...
или, относя
все величины к
...
где ... -
суммарная теплота, подведенная к
извне, ... - изменение внутренней энергии
...... - работа, затрачиваемая на расширение (... - объём, за-
нимаемый
При постоянном объёме ... = 0, ... = 0 или ......., т.е. вся
теплота, подводимая к газу, ..... тратится на увеличение его вну-
тренней энергии. Поэтому
...
Пренебрегая зависимостью ... от температуры и имея в виду,
что при .......0 ... = 0, имеем
...
Внутрення энергия является одной из функций состояния газа.
Используя формулы
...
Уравнение является математическим выражением первого закона
термодинамики.
Энтальпия. Введём ещё одну функцию состояния ..., определяе-
мую соотношением
...
Или, пренебрегая изменением ...,
...
Эта функция называется энтальпией. Из определения энтальпии
следует, что её приращение ... представляет собой приращение те-
плоты ... в процессе ... = ... Имея это в виду, из первого закона
термодинамики (...........................), интегрируя его в
предположении ..........., получим
...
Используя уравнение состояния (......) и соотношение .......,
имеем
...
Энтропия. При изучении течения газа часто используют понятие
энтропии. Эта функция определяется дифференциальным соотношением
...
Найдём связь между энтропией и энтальпией
...
из первого закона термодинамики
...
следует
...
...
...
...
... - тензор плоскости импульса.
...
...
Течение в трубе.
...
Оператор Лапласа
...