Тема 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. Нет ограничений, связанных с нелинейностью2. Описание сложных физических процессов3. Описание эволюции течения 1. Погрешность округления2. Проблемы задания ГУ3. Стоимость ЭВМ
Основные этапы математичкеского моделирования
Рис.
Структурные элементы математического моделирования вместе со связями показаны на рисунке.
Математическая постановка задачи базируется на физической модели рассматриваемых течений, которая строится на основе имеющихся данных об объекте исследования.
Характеризующие математическую модель исходные уравнения и граничные условия с помощью конечно-разностных методов преобразуются в дискретную модель.
В результате реализации дискретной модели на одном из... программирования программу для ЭВМ. Решение тестовых задач и анализ результатов позволяет убедиться в работоспособности разработанных алгоритмов и программ.
Решение конкретных задач и анализ полученных результатов позволяет судить об эффективности и применимости разработанных алгоритмов.
Если обнаружится несоответствие расчётных и экспериментальных данных - это значит, что физическая модель, математическая модель или дискретная модель не адекватны изучаемому объекту. В этом случае проводятся дополнительные исследования. Процесс исследования продолжается до момента устранения.
4. Три этапа развития вычислительной аэрогидродинамики
Этап(урав.) Результ.расчёт. Год расчёта ЭВМ Времярас-чёта
профиль реальн.ком.
IУр.потен. 1. Распр. давл.2. Индук. сопр. 1930 1968 IBM 360CDC 6600 5 м
IIУр.Эйлера 1. Трансзвук2. Гиперзвук 1971 1976 3707600 5
IIIУр. Н.-С. 1. Отр. потока2. Турб. 1975 1985 CRAV 5
Рис. Рис.
Рис. Методы расчета параметров течения
Тема 2
Физические свойства жидкостей и газов
1. Различные состояния вещества. Твёрдые тела, жидкости и газы. Силы, действующие на жидкости.
2. Основные свойства реальных жидкостей.
3. Поверхностное натяжение.
4. Уравнение состояния. Адиабата Тэйда.
1. Различные состояния вещества. Твёрдые тела, жидкости и газы
В природе различают четыре агрегатных состояния вещества: твёрдое, жидкое, газообразное и плазменное. Жидкость занимает промежуточное положение между твёрдыми телами и газами. Свойства жидкостей при низкой температуре и высоком давлении близки к свойствам твёрдых тел, а при высокой температуре и низком давлении - к свойствам газов.
Жидкость, как всякое тело, имеет молекулярное строение, т.е. состоит из молекул, объём пустот между атомами намного превосходит объём самих молекул. Причём в жидкостях и твёрдых телах объём пустот между молекулами меньше, а межмолекулярные силы больше, чем в газах. Виду бесконечной малости молекул и пустот между ними по сравнению с рассматриваемыми объёмами жидкости можно рассматривать жидкости и газы в виде... сплошной среды, придавая ей свойства непрерывности.
Жидкость - это физическое тело, обладающее лёгкой подвижностью частиц, текучестью и способное изменять свою форму под воздействием внешних сил.
Жидкости разделяют на сжимаемые (газообразные) и несжимаемые или весьма малосжимаемые (капельные).
Для облегчения изучения законов движения жидкости вводят понятия идеальной и реальной жидкости.
Идеальные - невязкие жидкости, обладающие абсолютной подвижностью, т.е. отсутствием сил трения и касательных напряжений и абсолютной неизменностью а объёме под воздействием внешних сил.
Реальные - вязкие жидкости, обладающие сжимаемостью, сопротивлением растягивающим и сдвигающим усилиям и достаточной подвижностью, т.е. наличием сил трения и касательных напряжений.
Реальные жидкости могут быть ньютоновскими и неньютоновскими (бингамовскими). В ньютоновских жидкостях при движении одного слоя жидкости относительно другого величина касательного напряжения пропорциональна скорости сдвига. При относительном покое эти напряжения равны нулю. Такая закономерность была установлена Ньютоном в 1686 году, поэтому эти жидкости (вода, масло, бензин, керосин, глицерин и др.) называют ньютоновскими жидкостями.
Неньютоновские жидкости не обладают большой подвижностью и отличаются от ньютоновских жидкостей наличием касательных напряжений (внутреннего трения) в состоянии покоя. Эта особенность была подмечена Ф.Н.Шведовым (1889), а затем Бингемом (1916), поэтому такие жидкости (битум, гидросмеси, глинистый раствор, коллоиды, нефтепродукты при температуре, близкой к температуре застывания) получили и другое название - бингемовские.
Силы, действующие в жидкости, принято делить на внутренние и внешние.
Внутренние силы представляют собой силы взаимодействия частиц жидкости, внешние силы делятся на силы поверхностные и объёмные.
Поверхностные силы (сжатие, давление, растяжение, силы трения) приложены к поверхностям, ограничивающим объём жидкости.
Объёмные силы (например, сила тяжести, сила инерции, электромагнитная сила) распределяются по всему объёму жидкости.
С термодинамической точки зрения состояние жидкости или пара характеризуется тремя параметрами: давлением..., плотностью... и температурой Т, связанными между собой уравнением состояния.
Исходной единицей давления в Международной системе единиц СИ является паскаль:
...
На практике используют более крупные единицы - гектапаскаль (1гПа =... Па), килопаскаль (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. Виды давления ( барометрическое, абсолютное, избыточное, манометрическое )
Различают следующие виды давления: барометрическое, абсолютное, манометрическое и вакуумметрическое.
Барометрическое (или атмосферное) давление... зависит от места над уровнем моря и от погоды. За нормальное барометрическое давление принимают давление, равное 760 мм рт. ст., что соответствует 101325... С высотой барометрическое давление убывает. В глубоких шахтах барометрическое давление значительно больше, чем на уровне моря.
Давление, вычисляемое по соотношению
...
называется абсолютным.
Аналогичное давление в точке равно сумме внешнего поверхностного и весового давления.
Если к свободной поверхности приложено барометрическое давление..., то есть... и основное уравнение гидростатики перепишем так
...
Давление
...
носит название манометрического или избыточного. Таким образом, манометрическим давлением называется разность между абсолютным давлением... и барометрическим..., если...
Если в данной точке жидкости абсолютное давление меньше барометрического, то разность между барометрическим и абсолютным давлениями называется вакуумметрическим давлением...
Итак, если..., то
...
Абсолютное давление... отрицательным быть не может, поэтому вакуумметрическое давление не может быть больше барометрического.
5. Приборы для измерения давления
Приборами для измерения барометрического давления служат барометры различных конструкций.
Для измерения манометрического давления служит манометр. Манометрическое давление можно измерить высотой столба жидкости. Сосуд наполнен жидкостью с плотностью.... Давление на свободной поверхности...
Пусть необходимо измерить давление на уровне 1-1. Если на этом уровне сделать отверстие и присоединить к нему стеклянную трубку П, то жидкость в этой трубе поднимется под действием давления на некоторую высоту... По основному уравнению гидростатики
...
откуда
...
Этой высотой... поднятия жидкости в трубке П можно измерять манометрическое давление. Трубка П называется пьезометром.
Hайдем соотношение между 1..., 1 м вод. ст. и 1 мм рт. ст.
При высоте вод. столба... = 1 м давление
...
При высоте ртутного столба... = 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) причину повышения подъемной силы крыла при заданной картине линий тока
...
Уменьшение давления в точках, где скорость потока больше, положено в основу водоструйного насоса. Струя воды подается в трубку, открывающуюся в атмосферу, так что на выходе их трубки давление равно атмосферному. В трубке имеется сужение, по которому вода идет с большой скоростью, вследствие чего давление в этом месте оказывается меньше атмосферного. Такое же давление устанавливается и в охватывающей трубку камере насоса, которая сообщается с трубкой через разрыв, имеющийся в узкой части трубки. Подсоединив к камере насоса откачиваемый объект, из него можно откачать воздух (или какой-либо другой газ) до давления порядка 100 мм рт. ст. Откачиваемый воздух захватывается струей воды и уносится в атмосферу.
...
Тема 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) сила тяжести жидкости...
Первые две силы горизонтальны, третья вертикальна.
Сума проекций импульсов этих сил на направление движения, т.е. на горизонтальную ось равна
...
Согласно теореме об изменении количества движения получаем
...
Сокращая на..., имеем
...
Откуда
...
Обозначив повышение давления...-... буквой Р, находим
...
Формула называется формулой Н.Е.Жуковского, который первый дал теорию гидравлического удара.
Разделим последнее соотношение на..., получим
... или ...
Из формулы видно, что при гидравлическом ударе повышение напора в трубопроводе равно....
Численное значение величины С также выведено Н.Е.Жуковским и определяется по следующей формуле
...
где... - плотность жидкости,.. - модуль упругости жидкости,... - модуль упругости стенок трубы,... - внутренний диаметр трубы,... - толщина стенки трубы.
Рассмотрим пример.
Пример. Определить повышение напора при гидравлическом ударе в чугунной трубе диаметром..., если толщина стенки трубы... = 0.0105 мм, модуль упругости воды... = 2..., модуль упругости чугуна....., а скорость течения....
Решение.
По формуле находим скорость распространения ударной волны
...
Найдём повышение напора
...
Гидравлический удар может повредить трубы. Для предотвращения разрушения труб применяются следующие меры.
1. Из формулы... видно, что увеличение давления пропор-
ционально скорости течения..., поэтому в трубопроводах не следует
допускать больших скоростей без принятия соответствующих предохранительных мер.
2. Причиной гидравлического удара является быстрое закрытие крана. При продолжительности закрытия.......... повышение давления равно... (так называемый прямой гидравлический удар). При продолжительности......... повышение давления меньше...... (непрямой гидравлический удар).
Продолжительность закрытия... (в секундах) может быть подсчитана по формуле Н.Е.Жуковского
..... или...................,
где... - плотность жидкости,... - скорость течения,... - длина трубопровода,... - допустимое повышение напора столба жидкости (в метрах).
Время закрытия трубопровода... прямо пропорционально длине трубопровода.... Т.е. чем длинее трубопровод, тем длительнее должно быть закрытие кранов и задвижек.
3. Для уменьшения вредного действия давления при гидравлическом ударе ставят предохранительные клапаны, которые, открываясь при определённом давлении, предохраняют провод от разрушения.
4. Кроме предохранительных клапанов, для уменьшения давления применяют воздушные колпаки. В момент повышения давления жидкость входит в колпак и сжимает находящийся в нём воздух, что уменьшает повышение давления.
Пример. Определить продолжительность закрытия задвижки на трубопроводе, если длина трубопровода... = 800 м,... = 3..., допускаемое давление в трубопроводе 1 000 000..., а гидростатическое давление Р = 200 000....
Решение.
Допускаемое повышение давление от гидростатического удара
... = 1 000 000 - 200 000 = 800 000...
Продолжительность закрытия задвижки
...
2. Вытекание жидкости при переменном уровне
Рассмотрим случай истечения жидкости из открытого сосуда в атмосферу через отверстие площадью....
Струя при вытекании через отверстие постепенно сжимается. Ближайшее к отверстию наименьшее живое сечение С-С, в котором движение можно рассматривать плавно изменяющимся, называется сжатым сечением. Обозначим площадь сжатого сечения С-С...
...
...
Отношение
...
(... = 0.64 для круглого отверстия)
называется коэффициентом сжатия.
Обозначим через... высоту уровня жидкости над центром тяжести отверстия,... - скорость в сжатом сечении.
Запишем уравнение Бернулли для сечений О-О и сжатого сечения С-С.
...
где... - скорость свободной поверхности,
... - потери напора при вытекании через отверстие, они
определяются из соотношения
...
Пренебрегая величиной... (ввиду её малости по сравнению с Н), получаем
...
отсюда скорость истечения
...
где... - коэффициент скорости (.....0.97).
Для определения расхода надо скорость умножить на площадь сжатого сечения:
...
по формуле... , откуда
...
тогда расход, выраженный через... равен
...
где... - коэффициент расхода (... = 0.62).
Рассмотрим вытекание жидкости из ёмкости при переменном уровне. Движение в данном случае является неустановившимся. С достаточной для практики точностью можно считать, что в каждый момент времени скорость вытекания определяется соответствующим этому моменту напором Н так же, как и при установившемся движении.
...
...
Определим время, в течение которого жидкость опустится на...-...
Рассмотрим промежуточное положение уровня с напором Н. За время... вытечет объём жидкости, равный
...
За это время... напор изменится на (-...Н). Объём жидкости, вытекшей из сосуда, равен
...
где... - площадь свободной поверхности в сосуде.
Приравнивая выражения, получаем
...
откуда
...
Интегрируя, находим
...
При постоянной площади свободной поверхности
...
Пример. Вычислить продолжительность опорожнения цистерны при её диаметре... = 2 м и длине... = 5 м, если диаметр сливного отверстия... = 0.1 м, а коэффициент расхода... = 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 кг газа в этом процессе, выразим через приращение температуры газа... :
...
Множитель С, представляющий собой количество теплоты, необходимое для подогрева 1 кг газа на 1 град в данном процессе, называется удельной теплоёмкостью.
Удельная теплоёмкость существенно зависит от характера процесса.
Рассмотрим теплоёмкости, соответствующие процессам, происходящим при постоянном объёме... и давлении.... Зависимость между удельными теплоёмкостями идеального газа... и... определяется следующим соотношением.
...
В термодинамике и газодинамике важное значение имеет отношение теплоёмкостей...... Величина... зависит от структуры молекулы газа. Так, для идеальных одноатомных газов... = 1.66, для двухатомных газов, в том числе и для воздуха,... = 1.4.
3. Первый закон термодинамики.
Пусть некоторое количество газа находится в равновесии. Обозначим через... количество подведённой к газу извне теплоты. В общем случае подвод теплоты приводит к изменению внутренней энергии газа... и объёма. ПРи изменении объёма газ совершает внешнюю работу, равную.... Поэтому
...
или, относя все величины к 1 кг массы газа, получаем
...
где... - суммарная теплота, подведенная к 1 кг массы газа извне,... - изменение внутренней энергии 1 кг массы газа,...... - работа, затрачиваемая на расширение (... - объём, занимаемый 1 кг массы газа).
При постоянном объёме... = 0,... = 0 или......., т.е. вся теплота, подводимая к газу,..... тратится на увеличение его внутренней энергии. Поэтому
...
Пренебрегая зависимостью... от температуры и имея в виду, что при.......0... = 0, имеем
...
Внутрення энергия является одной из функций состояния газа. Используя формулы
...
Уравнение является математическим выражением первого закона термодинамики.
Энтальпия. Введём ещё одну функцию состояния..., определяемую соотношением
...
Или, пренебрегая изменением...,
...
Эта функция называется энтальпией. Из определения энтальпии следует, что её приращение... представляет собой приращение теплоты... в процессе... =... Имея это в виду, из первого закона термодинамики (...........................), интегрируя его в предположении..........., получим
...
Используя уравнение состояния (......) и соотношение......., имеем
...
Энтропия. При изучении течения газа часто используют понятие энтропии. Эта функция определяется дифференциальным соотношением
...
Найдём связь между энтропией и энтальпией
...
из первого закона термодинамики
...
следует
...
...
...
...
... - тензор плоскости импульса.
...
...
Течение в трубе.
...
Оператор Лапласа
...