А.Б. Мазо, К.А. Поташев

ГИДРОДИНАМИКА

Казань – 2008

-2А.Б. Мазо, К.А. Поташев

ГИДРОДИНАМИКА

Учебное пособие

для студентов нематематических факультетов

Казань – 2008

-3УДК 534

Печатается по решению учебно-методической комиссии

механико-математического факультета Казанского государственного университета Протокол № 1 от 2 октября 2008 года Рецензент – член-корреспондент РАН, доктор физико-математических наук Д. А. Губайдуллин.

Мазо А.Б., Поташев К.А. Гидродинамика. Учебное пособие. – Казань:

КГУ, 2008. – 126 с.

Учебное пособие основано на специальном курсе лекций и лабораторных занятий для студентов III курса географического факультета Казанского университета. Основной задачей авторов при написании данного учебного пособия было представление теоретических основ гидродинамики в форме, доступной студентам, для которых математика и механика не являются профильными дисциплинами. Это обусловило несколько упрощенный стиль изложения материала; в частности, в пособии практически не используется аппарат тензорного исчисления, почти все уравнения выводятся только в декартовой системе координат. Тем не менее, все традиционные модели гидродинамики и теплообмена – от простейших потенциальных течений до осредненных по Рейнольдсу уравнений Навье-Стокса – нашли свое место в учебном пособии.

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

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

© Казанский государственный университет. -4Оглавление Лекция 1. Введение

ВВЕДЕНИЕ

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

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

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

Жидкость характеризуется несколькими физическими параметрами.

,, c, Это плотность динамическая вязкость теплоемкость теплопроводность. Дадим краткое неформальное описание этих свойств и укажем их размерности.

Плотность жидкости определяется как предел отношения массы жидкости к ее объему:

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

-6При этом скорость деформации (скорость сдвига) равна e = du / dy, с 1, а для реализации такого течения необходимо приложить напряжение, Па=H/м2.

Рис. 1.1. К определению сдвиговой вязкости жидкости Закон трения Стокса выражает линейную связь между напряжением и скоростью деформации, причем коэффициент пропорциональности – вязкость очевидно, различен для различных жидкостей. По опыту мы знаем, что вязкость масла больше, чем вязкость воды. Вязкость жидкости обычно зависит от ее температуры T ; обычно убывает с ростом T.

изменении его температуры на 1 градус.

Теплопроводность, Вт/(м K) характеризует способность жидкости dT / dx = lim (T2 T1 ) /( x2 x1 ) в линейном законе Фурье Знак минус указывает на то, что тепловой поток направлен от горячего к холодному: если температура возрастает с ростом x, то dT / dx > 0, а q < 0, т.е. тепло течет в противоположном направлении от направления роста температуры.

При теоретическом изучении течений жидкости в каждой точке следующие величины: скорость v = ( v x, v y, v z ) = ( v1, v2, v3 ) = (u, v, w), температуру T, давление p и внутренние напряжения. По этим характеристикам определяют гидродинамическое взаимодействие потока жидкости с границами потока и обтекаемыми телами. Например, сопротивление труб, подъемная сила крыла самолета, плавучесть корабля и его сопротивление при движении. В применении к расчетам природных процессов данные характеристики позволяют определить, например, разрушительное действие потоков на поверхности (подвижки ледников, сели, прибой морских волн, цунами, ветровая эрозия и пр.) В гидродинамике рассматриваются математические модели течений жидкости и газа в различных условиях. Эти модели, как правило, представляют собой системы дифференциальных уравнений в частных производных и краевые (начальные и граничные) условия. Чтобы с помощью модели предсказать, как будет проходить (или как проходило ранее) течение жидкости, надо решить эти уравнения, а точнее, краевые задачи. В некоторых простейших частных случаях решение этих задач получается в аналитической форме (т.е. в виде формул), однако большинство реальных течений моделируется на основе численного решения задач на компьютерах.

';

-8Задачи к лекции Задача 1.1. Над твердой поверхностью протекает слой воды толщиной h = 10см. Определите величину сдвигового напряжения в жидкости, если известно, что через поперечное сечение шириной a = 10см жидкость протекает с объемным расходом Q = 200мл/c. Вязкость воды равна = 1.5 103 Па с. Профиль продольной скорости считать линейным по высоте (рис. 1).

Задача 1.2. Зависимость вязкости воды от температуры может быть приближенно описана уравнением ( Па с ) e в температурном диапазоне от 5 C до 50 C. Определить, как изменится величина напряжения в условиях задачи 1.1 при температуре воды T = 20 C.

Задача 1.3. Записать уравнение для определения полного теплового потока вдоль отрезка, соединяющего точки с температурой T0 и T1.

КРАТКИЕ СВЕДЕНИЯ ИЗ МАТЕМАТИКИ И МЕХАНИКИ

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

Скаляр – это величина, значением которой является число.

Функция f ( x ) – отображает числовой аргумент x в числовое значение функции y = f ( x ). Функцию представляют в виде графиков и таблиц (в одной колонке аргумент, в другой – значение функции). Наряду с функцией одной переменной, рассматривают функции нескольких переменных y = f ( x1, x1,... xn ). Функция многих переменных отображает набор числовых аргументов x1, x1,... xn в одно числовое значение функции y = f ( x1, x1,... xn ). Например, уравнение состояния газа задает давление как функцию плотности и температуры: p = F (, T ). Функция двух переменных обычно изображается в виде карты изолиний (линий уровня) функции. Например, функция T ( x, y ), задающая температуру на поверхности в декартовых координатах x, y обычно представляется картой изотерм T ( x, y ) = const, см. рис. 2.1.

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

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

- 10 y Вектор v - это объект, который характеризуется абсолютной величиной v и ориентацией в пространстве (направлением). Сам вектор не зависит от системы координат, однако если какая-либо система координат выбрана, то можно определить проекции вектора на три координатные оси, и это будут компоненты вектора в данной системе координат. В прямоугольной декартовой системе координат x, y, z с ортами i, j, k компоненты вектора ( v1, v2, v3 ) = v x, v y, v z - это скалярные функции, которые однозначно задают вектор а его модуль (длина) v вычисляется по формуле Поле векторной величины – распределение вектора над областью изменения аргумента. Для определения векторного поля v достаточно задать систему координат и скалярные поля компонент вектора. В случае плоских течений это два скалярных поля ( v x, v y ), а в случае объемных течений - три скалярных поля ( v x, v y, v z ) его компонентов.

1. Умножение на скаляр: w = a v, где a – число или скалярная величина. Каждая компонента w есть каждая компонента v, умноженная на a : wi = avi, i = 1, 2,3.

2. Сложение двух векторов. u = v + w означает ui = vi + wi.

3. Скалярное произведение векторов определяется формулами Очевидно, скалярное произведение ортогональных векторов равно нулю, а скалярное произведение колинеарных (параллельных) векторов равно произведению их длин.

4. Квадрат вектора есть скалярное произведение его на себя:

5. Векторное произведение u v = w есть вектор, определяемый как определитель Вектор w направлен ортогонально к плоскости, натянутой на вектора u, v, а его длина равна произведение колинеарных векторов равно нулю.

Векторное поле изображают разными способами: 1) стрелками; 2) линиями тока; 3) эпюрами в проекциях.

Оператор действует на скалярную функцию или вектор, в результате - 12 d дифференцирования, действуя на функцию f ( x ) дает новую функцию – ее производную g ( x ) = f ( x ).

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

Например, оператор набла ( ) в декартовой системе трех координат можно представить как вектор, компонентами которого являются операторы дифференцирования по координатам, Если функция f(x,y,z) - скалярная, то применяя этот оператор, получаем вектор Градиентом скалярной функции f называется вектор наискорейшего роста функции.

Если оператор набла применяется к вектору, то результат равен скалярному произведению оператора-вектора на вектор v = v x, v y, v z :

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

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

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

Уравнения движения в форме Громеки - Ламба Запишем уравнение Навье-Стокса для x -компоненты скорости.

Левая часть представляет собой субстанциональную производную.

Формулировка Громеки-Ламба основана на том, чтобы выразить субстанциональную производную скорости через кинетическую энергию V 2 / 2 и завихренность.

поэтому Последние два слагаемых можно записать в виде - 44 и окончательно получаем уравнение движения в форме Громеки – Ламба для х - составляющей вектора скорости Остальные два уравнения для v и w получаются аналогично и имеют вид Удобно записать все три уравнения движения в индексной форме:

Замечание. В научной литературе можно встретить формулировку, при которой используется 2 ( v w ) x. Это связано с тем, что раньше завихренность находили как от того, что используем мы.

Уравнения в форме Громеки – Ламба удобны при рассмотрении невязких течений, когда член ui исчезает. Как увидим позднее, из этой формы получается интеграл Бернулли.

Уравнение движения в форме переноса завихренности Эта форма особенно удобна и применяется для описания плоских течений несжимаемой жидкости (w=0, =const). Основная идея введения новых преобразованных переменных состоит в том, чтобы вместо трех уравнений в естественных переменных (двух уравнений Навье - Стокса и уравнения неразрывности) получить 2 уравнения для завихренности Из этого определения следует, что функция может быть задана с точностью до произвольной постоянной. В самом деле, Такие функции, непосредственные значения которых не важны (именно поэтому они определяются с точностью до произвольной постоянной), а важны лишь их производные, называются потенциалами.

Физически задает расход жидкости. Изолинии функции тока =const являются линиями тока, т.е. такими линиям, касательные к которым в каждой точке совпадают по направлению с вектором скорости.

Чем чаще линии тока, тем больше расход в данном сечении.

Рис. 5.2. Пример линий тока в термоконвективных течениях.

Функция тока замечательна тем, что определенные по формуле (5.7) неразрывности. В самом деле, - 46 Для функции тока получается простое уравнение. Продифференцируем u = / y по y, а v = / x по x и вычтем второе из первого.

Получим уравнение Пуассона с вихрем в правой части, или Осталось вывести уравнение для завихренности. Используем для этого уравнения Навье-Стокса Продифференцируем (5.10) по y, а (5.11) – по x и вычтем из второго первое. Получим Итак, вместо трех уравнений Навье-Стокса (2 уравнения для u, v и уравнение неразрывности для преобразованных переменных – завихренности (5.12) и функции тока (5.9):

При этом давление не входит в постановку задачи. Если требуется найти u, v, p, то в рамках модели переноса завихренности поступают так:

заодно по формулам (5.8) и скорости u, v. После того, как задача решена, для отдельных моментов времени решается самостоятельная задача для давления. Уравнение для p получают из уравнений импульсов (5.10), (5.11). Первое дифференцируют по x, второе по y и складывают. В результате получается уравнение Пуассона для давления Краевые задачи для уравнений в преобразованных переменных В параграфе 4.3 мы формулировали начальные и граничные условия для уравнений Навье-Стокса в естественных переменных для задачи обтекания тела в канале. Ниже будет дана постановка краевых условий для той же задачи в преобразованных переменных.

Начальные условия. Поскольку уравнение для стационарно, начальное условие нужно лишь для. Если в начальный момент жидкость покоится, то и вихрей нет, Граничные условия на твердых стенках. Это верхняя и нижняя стенки канала, y = ± H, а также поверхность тела. Задается условие прилипания, т.е нормальная и касательная компоненты скорости равны нулю. В терминах функции тока, согласно определению (5.7), эти компоненты запишутся так:

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

Здесь Q - полный расход жидкости в канале.

Граничное условие (5.18) справедливо лишь для случая, когда тело стационарен. В общем случае определение значения функции тока на контуре представляет собой отдельную проблему. Вообще говоря, математическое решение сформулированной задачи существует при любом заданном значении, поэтому для однозначного выбора этой константы необходимо дополнительное условие. Данное условие сводится к требованию непрерывности давления преобразованных переменных, воспользуемся исходными уравнениями Навье-Стокса (5.10), (5.11), записав их в непосредственной близости к контуру и выбрав в качестве пространственных координат длину дуги s и нормаль n вместо x и y. Поскольку на контуре из-за прилипания скорость обращается в ноль, получим Здесь для простоты мы отбросили массовые силы. В выбранной локальной системе координат s, n с началом в точке контура оператор Лапласа от компонент скорости будет содержать лишь производные по нормали:

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

Условия (5.16) – (5.18), когда на границе задана искомая функция, называются граничными условиями 1-го рода, или главными граничными условиями, или условиями Дирихле. Все эти граничные условия относятся к уравнению для функции тока, в то время как для вихря нет физически содержательных граничных условий на стенках. Вместо них имеются «лишние» граничные условия Неймана для, выраженные вторым равенством (5.14), а именно - 50 Граничные условия на входе и выходе канала.

На входе x = L1 задан профиль продольной скорости u0 ( y ), причем такой, что расход равен Q. Это значит, что Согласно определениям функции тока и завихренности на входе имеем поэтому граничные условия для входного сечения канала можно записать так:

В частности, если поток однороден, то u0 = Q /(2 H ), и условия (5.20) принимают вид На выходе обычно задают «мягкие» граничные условия вида Либо, если заранее известно, что течение там установилось, используют аналитические решения для и.

Задачи к лекции Задача 5.1. Дано плоское течение несжимаемой жидкости:

неразрывности и является безвихревым.

Задача 5.2. Записать уравнения движения в форме Громеки-Ламба для случая плоского течения жидкости u, v 0, w = 0 и удостовериться в том, что все слагаемые имеют одинаковую размерность Задача 5.3. Доказать, что при установившемся движении среды ui t = 0 линии тока и траектории совпадают.

Задача 5.4. Определить и отобразить на графике линии тока плоского течения и описать движение среды по полю скоростей u = K x, Задача 5.5. Имеется плоский поток несжимаемой жидкости, в котором u = A y / r 2, где r 2 = x 2 + y 2. Найти во всем потоке компоненту v, если v ( x = 0, y ) = 0, при всех значениях y. Показать, что движение безвихревое, а линии тока – окружности.

- 52 Лекция 6.

ВАЖНЫЕ ЧАСТНЫЕ ФОРМУЛИРОВКИ УРАВНЕНИЙ

ТЕРМОГИДРОДИНАМИКИ

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

1. Уравнения Навье – Стокса (3 уравнения для компонент вектора скорости u, v, w).

2. Уравнение неразрывности (для давления p ).

3. Уравнение конвективной теплопроводности (для расчета температуры T ).

4. Уравнение состояния, которое связывает p,, T (для идеальных газов это уравнение имеет вид p = RT ).

Имеем 6 уравнений и 6 неизвестных (u, v, w,, T, p).

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

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

Установившееся (квазистационарное) движение Движение называется установившимся, если в каждой точке пространства (x,y,z) его параметры (u, v, w,, T, p) не меняются со временем.

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

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

Для трехмерных задач Навье – Стокса в естественных переменных это приводит к уравнениям Важное свойство установившихся течений состоит в том, что линии тока неизменны и совпадают с траекториями движения частиц. Частным случаем установившегося течения является равновесие жидкости. Это происходит, когда течения нет, т.е. v 0.

Гидростатика Раздел гидродинамики, который изучает равновесие жидкости называется гидростатикой. Уравнения гидростатики получаются из общих уравнений (6.1), если в них формально положить v 0. При этом из (6.1) получаем 4 уравнения или Здесь gx, gy, gz – компоненты вектора ускорения массовой силы.

Закон Паскаля В случае, когда массовые силы отсутствуют, т.е. g=0, из этих уравнений получается, что откуда следует, что p = const. Это решение носит название закона Паскаля, который гласит, что в покоящейся жидкости (газе) при отсутствии массовых сил давление постоянно.

Гидростатическое давление у поверхности Земли Пусть массовая сила представляет собой гравитационное поле Земли, тогда ускорение равно Ось z направлена вверх, а ускорение вниз, поэтому знак минус. Из (6.2) получаем Первые два уравнения констатируют, что давление не зависит от плановых координат x, y, а изменяется только с высотой z. Если проинтегрировать третье уравнение (6.4) от 0 до h, где ноль совпадает с поверхностью земли, а через p0 обозначить давление на этой поверхности, то получим интеграл в правой части выражения (6.5) представляет собой вес столба жидкости высотой h на основании 1 м2. Из (6.5) мы видим, что давление жидкости в сосудах не зависит от их формы, а определяется только высотой. Данный факт иллюстрируют следующие примеры.

2) На чашки весов, выполненные в виде свободно скользящих по стенкам поршней, давят одинаковые силы gh. Весы будут в равновесии.

- 56 h Равновесие идеального газа в поле силы тяжести.

Барометрическая формула Уравнение состояния идеального газа плотность газа в зависимости от температуры Подставим (6.6) в (6.4) и получим Это обыкновенное дифференциальное уравнение первого порядка для давления p( z ). Оно решается методом разделения переменных:

Отсюда получаем барометрическую формулу Она даёт закон изменения давления с высотой, если известно распределение температуры по высоте. В частном случае, если считать атмосферу изотермической, когда Т=const, из формулы (6.8) получаем убывает с высотой по экспоненциальному закону.

Стандартная атмосфера Если атмосфера неизотермична, источников тепла нет ( f = 0 ), а коэффициент температуропроводности постоянен ( a =const), то уравнение для температуры (6.2) будет иметь следующий вид:

Общее решение этого уравнения где T0 - температура на поверхности Земли, k – коэффициент, который определяется на основе измерений. Для атмосферы k = 0.650/100 м.

Эту формулу можно применять до высоты 11 км (тропосфера), а выше нее полагают T=const=-560C (стратосфера). Описанную модель называют стандартной атмосферой, ее используют для имитации высотных полетов. Но процессы в атмосфере всегда нестационарные, на распределение давления и температуры влияют излучение Солнца, геотермический поток, изменение состава воздуха при ионизации.

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

Несжимаемая атмосфера Пусть плотность атмосферы постоянна, =const. Интегрируем уравнение (6.4) по z и получаем линейное убывание давления с высотой:

С помощью формулы (6.11) легко вычислить точку, в которой давление равно нулю:

Это грубая оценка для толщины атмосферы.

Закон Архимеда Рассмотрим сосуд, заполненный жидкостью. Поместим в него тело объемом V, с поверхностью S, внешнюю нормаль к S обозначим через n. Считаем, что жидкость покоится, поэтому сила A, действующая на тело, сводится к давлению:

Здесь применена формула Остроградского-Гаусса. Учтем, что у поверхности Земли при гидростатике p = g = ( 0,0, g ). Получим выражения для компонент вектора силы Архимеда объеме тела и действует вверх. Точнее говоря, сила Архимеда направлена в сторону, противоположную вектору массовых сил. В случае, когда в качестве такой силы выступает инерция (например, при стационарном вращении жидкости в центрифуге), сила Архимеда будет направлена к оси вращения. На этом основаны сепараторы, предназначенные для быстрого извлечения примесей из жидкости.

Эта формула справедлива только для гидростатики. Кроме того, для применения формулы (6.14) необходимо, чтобы поверхность тела S была замкнута.

Пример 1. Подводная лодка.

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

- 60 Пример 2. «Вечный двигатель».

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

Кажется, что за счет силы легкого (плавучего) материала, будет вращаться, т.к. его левая часть «захочет» всплыть. Но на самом деле он не вращается, так давления направлена через центр момент равен нулю. Это явление называется парадоксом Жуковского.

Задачи к лекции давлений в точках, находящихся на осях труб А и В (рис 1), заполненных водой, если разность уровней ртути h = 30 см. Разность рт = 13600 кг / м3.

Ответ: 46844 Па.

по дифференциальному манометру (рис. 2), равна h = 50 см. Определить разность давлений в точках А и В, если H = 30 см. Трубы заполнены водой.

Ответ: 58800 Па.

Задача 6.3. Определить глубину погружения h плавающего в воде деревянного ( Д = 800 кг / м3 ) конуса, если его высота равна 3 м.

Ответ: 2.785 м (по закону Архимеда вес вытесненной воды равен весу конуса).

Задача 6.4. Определить равнодействующую величину давления и положение центра давления воды на вертикальную пластину высотой h и шириной b, верхняя граница которой совпадает с поверхностью воды.

Ответ: x=0, y=2h/3.

Задача 6.5. Известно, что плотность среды зависит от температуры.

При этом коэффициент теплового расширения определяется следующим образом откуда 1 = 0 e t T или 1 = 0 (1 t T ) мало при малых изменениях плотности.

Определить плотность воды в открытой емкости при температуре окружающей среды t = 30°C, если t воды = 6 104 K 1.

- 62 Задача 6.6. Известно, что плотность жидкости зависит от давления.

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

p воды = 0.45 109 Па 1.

Ответ: ( / 0 = 1.046, 1бар = 105 Па ).

Задача 6.7. Какова будет плотность воды в паровом котле при температуре t1 = 152°C и давлении p1 = 5 бар ?

Ответ: = 909 916 кг / м3.

УСТАНОВИВШИЕСЯ ТЕЧЕНИЯ ИДЕАЛЬНОЙ ЖИДКОСТИ.

ИНТЕГРАЛ БЕРНУЛЛИ

характеристики не меняются со временем. Математически это выражается в том, что частные производные по времени равны нулю, /t=0, где – разные гидродинамические величины (скорость, давление и т.д.).

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

Запишем уравнения движения в форме Громеки – Ламба Для установившегося течения v = 0, т.к. вязкие напряжения в идеальной жидкости отсутствуют.

Выделим в жидкости линию тока. Вектор скорости направлен по касательной к этой линии.

Рис. 7.1. Ориентация векторов скорости v, завихренности - 64 Запишем проекцию уравнения движения (7.1) на выделенную линию тока. Оператор градиента перейдет в производную по направлению:

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

допущениях.

1) Массовые силы потенциальны, т.е. существует такая функция G, что F = grad G, следовательно 2) Можно ввести функцию давления Р, такую, что Подставив (7.3) и (7.4) в (7.2), получим уравнение Раз производная вдоль линии тока равна нулю, значит, выражение в фигурных скобках постоянно вдоль этой линии, интегрирования зависит от выбранной линии. Равенство (7.5) называют первым интегралом уравнения движения или интегралом Бернулли.

Примеры функции давления Р.

1. Несжимаемая жидкость, = const. Очевидно, что P = p /, и 2. Совершенный (идеальный) газ, p = RT, = p / RT. Подставив это выражение в определение функции давления, получим В изотермическом случае, когда T=const (не зависит от давления), имеем Пример потенциала массовых сил G.

Пусть массовая сила F – это сила тяжести, g - ее ускорение (ось z направлена вертикально вверх). Тогда Использование интеграла Бернулли при решении задач гидравлики Интеграл Бернулли справедлив для любых течений, включая сильные разрывы в виде скачков уплотнения. Поэтому первый интеграл широко используется при решении самых разнообразных задач гидродинамики и гидравлики.

Если значение постоянной Бернулли известно в какой-либо точке потока, то уравнение (7.5) позволяет вычислить скорость и давление в любой другой точке на линии тока. Ниже представлены несколько примеров.

- 66 Пример 1. Истечение несжимаемой жидкости из бака.

Считаем, что =const, следовательно, функция давления и потенциал массовых сил равны В точке 1 на свободной поверхности воды давление равно атмосферному, а скорость жидкости близка к нулю: p1 = p атм, v1 = 0. Вычисляем константу Бернулли (7.5) для произвольной линии тока, проходящей с поверхности до выходного отверстия (точка 2 на рисунке).

Но в точке 2 давление также равно атмосферному, поэтому следовательно вытекать из нее со скоростью Пример 2. Водослив через плотину.

Требуется вычислить скорость воды в точке 2. Свободную поверхность потока принимаем за контрольную линию тока в интеграле потенциальность массовых сил. Вычислим постоянную Бернулли i для точек 1 и 2.

Если точка 1 удалена от плотины достаточно далеко вверх по течению, то v2, так что можно считать, что v1 = 0. Получаем - 68 Пример 3. Обтекание крылового профиля.

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

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

Возьмем две точки – 1 и 2 на верхней и нижней поверхности крыла, пусть толщина профиля 1 м. Экспериментальные исследования показали, что v1 > v, а v 2 < v. Предположим, что v1 =105 м/с, v 2 = 95 м/с.

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

Перепад давления, обеспечивающий подъемную силу крыла, равен Видно, что p > 0, т.е. p 2 > p1.

Перепад давления состоит из двух частей: гидростатический перепад гидростатический перепад p гст =1·10·1=10 Па. При v1 =105 м/с, v 2 = м/с динамический перепад давления равен В задачах обтекания тел динамический перепад давления обычно намного выше, чем гидростатический.

Чтобы вычислить подъемную силу F, надо умножить p на площадь крыла S. При S = 10 м 2 получаем Задачи к лекции Задача 7.1. Записать уравнение Бернулли и убедиться в том, что размерности всех его слагаемых совпадают.

значение u1=1 м/с, датчик в сечении (2) показывает значение скорости u2=3 м/с.

1. оценить среднюю величину ускорения на участке dx, если dx=10 см;

2. предполагая жидкость несжимаемой, определить расстояние между поверхностями в сечении (1), если в сечении (2) это расстояние составляет 5 см.

Задача 7.3. Вода бьет струей из фонтана на высоту 8 м. Фонтан имеет форму усеченного конуса (малое выходное сечение диаметром - 70 мм, большее (входное) сечение диаметром 50 мм. Определить расход воды и давление у нижнего основания конуса.

p = pa + 0.802 105 Н / м Задача 7.4. По горизонтальной трубе переменного сечения протекает d1 = d 2 = 100 мм, d 3 = 25мм, p1 = 0.3МПа. Трением пренебречь.

Ответ: h1 = h3 = 32.2 м, h2 = 11.12 м.

Задача 7.5*. Определить форму осесимметричного сосуда, уровень жидкости в котором при ее истечении из дна опускается равномерно (водяные часы – клепсиндры).

Задача 7.6. Жидкость с плотностью вращается как твердое тело в вертикальном цилиндрическом круговом сосуде с постоянной угловой скоростью вокруг оси цилиндра. В состоянии покоя жидкость имеет уровень h от дна сосуда и на свободной поверхности давление равно нулю. Радиус основания равен R.

жидкости б) Определить расстояние z0 вдоль центральной оси сосуда от его дна до поверхности жидкости.

в) С какой угловой скоростью надо вращать жидкость, чтобы обнажился центр дна сосуда?

г) Найти давление на дно сосуда.

- 72 Лекция 8.

ПОТЕНЦИАЛЬНЫЕ ТЕЧЕНИЯ ИДЕАЛЬНОЙ ЖИДКОСТИ.

ИНТЕГРАЛ КОШИ-ЛАГРАНЖА

Рассмотрим уравнение Громеки – Ламба при следующих предположениях.

1. Движение потенциально, т.е. существует такая функция – потенциал течения, что вектор скорости v = ( u, v, w ) определяется как v = grad =, или в декартовой системе координат Если движение потенциально, то в нем нет вихрей. В самом деле, вычислим одну компоненту вектора завихренности:

Аналогично получается, что x = 0 и y = 0, т.о. = 0.

2. Для всего потока можно ввести единую функцию давления Р, такую Заметим, что ранее при рассмотрении интеграла Бернулли функция давления вводилась для каждой линии тока. Для несжимаемой жидкости P = p/.

В общем случае ввести единую функцию давления можно для т.н.

баротропных течений, когда давление зависит только от плотности, p=p().

При сделанных предположениях из уравнений (8.1) получим Если массовые силы также потенциальны, т.е. F = G, то из уравнения (8.4) получается равенство Оно означает, что выражение в фигурных скобках не зависит от координат, но может зависеть от времени. Поэтому можно записать Здесь f – произвольная функция времени. Уравнение (8.5) называют интегралом Коши – Лагранжа.

Легко заметить, что для установившихся течений, когда /t=0, а f не зависит от времени, интеграл Коши – Лагранжа совпадает с интегралом Бернулли. Использование интеграла Коши-Лагранжа преследует те же цели, что и интеграл Бернулли, а именно: если потенциал известен, то можно найти давление. Обычно правую часть f определяют по известным значениям скорости и давления на границе потока. Так как функция f одинакова для всех точек потока, достаточно знать ее значение хотя бы в одной точке. Таким образом, в задаче является неизвестным лишь одна функция - потенциал скорости. Чтобы его найти, используют уравнение неразрывности и определение скорости через потенциал (8.2):

Подставив (8.7) в (8.6), получаем классическое уравнение Лапласа для потенциала - 74 В декартовой системе координат оно имеет вид Это уравнение стационарно, т.е. не зависит от времени. Потенциальное течение представляет важный класс физических содержательных потоков.

Например: волновые движения воды; движения воздуха при распространении в нем акустических волн; непрерывное движение жидкости при обтекании твердых тел; струйные течения. Потенциальное течение является идеализацией реальных течений и в реальной жизни встречается крайне редко. Потенциальность нарушается из-за вязкости, которая провоцирует появление вихрей.

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

Рассмотрим пример потенциального течения.

Движение шара в идеальной жидкости Будем изучать течение жидкости, вызванное движением шара радиуса a со скоростью V в декартовой системе координат x, y, z, причем шар движется в направлении оси x.

Рис. 8.1. К постановке задачи о движении шара в вязкой жидкости Выберем произвольную точку M на поверхности шара. Её скорость V можно разложить на две компоненты: нормальную vr = V cos и касательную v = V sin. Здесь - угол между радиус-вектором или нормалью в точке M и направлением движения (или осью x ). Очевидно, Требуется определить поле скоростей и давления в окружающей жидкости, вызванной этим движением. Воспользуемся уравнениями потенциального течения несжимаемой жидкости Это дифференциальное уравнение 2-го порядка имеет бесконечно много решений, выбрать одно из них можно с помощью граничных условий. Вдали от тела возмущение, вносимое шаром в неподвижную жидкость, затухает. То есть в бесконечности можно поставить граничное условие r = : | v |= 0. Бесконечность удобно представить как шар большого радиуса, и вместо градиента использовать производную по нормали:

Для того, чтобы сформулировать граничные условия на поверхности (в произвольной ее точке M ), разложим вектор скорости v = (V,0,0) на три vr = V cos, v z = 0, v = V sin. Нас интересует лишь радиальная (она же нормальная) компонента скорости Т.о. для потенциала мы получили классическую задачу Неймана. Для того, чтобы найти решение задачи, изучим вначале некоторые свойства гармонических функций и частные решения уравнения Лапласа удовлетворяющие уравнению Лапласа называют гармоническими.

гармоническая функция. В самом деле, координатам сами являются гармоническими.

2) Функция 2 = C1 + B, где и В – постоянные, удовлетворяет уравнению Лапласа, т.е. является гармонической функцией.

3) Функция уравнения Лапласа (8.10). В самом деле, Аналогично, Вычисляем вторые производные:

Сложив их, получим требуемое равенство Итак, мы показали, что функция является решением уравнения Лапласа. Такое решение описывает течение жидкости от точечного источника. Согласно первому свойству гармонических функций производная также будет гармонической. Такое решение уравнения Лапласа описывает течение от диполя с осью в направлении х. По второму свойству функция Cx / r 3 также гармоническая. С ее помощью мы построим краевой задачи (8.10), (8.11), (9.12), удовлетворив граничные условия за счет выбора константы C. Итак, ищем решение задачи в виде - 78 C x C cos Условию на бесконечности (9.11) функция удовлетворяет при любой константе C, поскольку Сможем ли мы с помощью выбора постоянной величины С удовлетворить условию (8.12) на поверхности движущегося шара? Для этого надо вычислить производную (8.16) при r = a и подставить её в граничное условие (8.12). Имеем откуда находим константу Подставив ее в формулу (8.15), окончательно получим потенциал течения координатам.

В плоскости z = 0 формулы (8.18) дают следующую картину течения Такие линии тока течение имеет с точки зрения Эйлера, когда наблюдатель (кинокамера) фиксирован в пространстве, а мимо него проплывает шар. Если же наблюдатель (кинокамера) связан с движущейся сферой, то он будет видеть совершенно другую картину (точка зрения Лагранжа), которая получается с помощью суперпозиции (наложения) двух течений:

1) течения (8.17), (8.18) – с точки зрения Эйлера, и 2) поступательного движения со скоростью ( V,0,0) и потенциалом Это значит, что обтекание сферы в системе координат, связанной с ее центром, будет описываться потенциалом - 80 Поле скоростей этого течения показано на следующем рисунке Рис. 8.3. Поле скорости при движении шара в идеальной жидкости.

Чтобы найти давление используем интеграл Коши – Лагранжа (8.5), который для несжимаемой жидкости без массовых сил примет вид Из условий на бесконечности ( v = 0, p = p ) находим f = p /. Для простоты в дальнейшем положим = 1.

Уравнение (8.20) справедливо в системе координат, когда отсчитываются от центра шара. Чтобы верно найти давление в жидкости, следует перейти к фиксированным в координатам. При постоянной скорости движения это достигается заменой Тогда интеграл Коши-Лагранжа (8.20) при установившемся движении шара в несжимаемой жидкости принимает вид Отсюда находим поле давления которое показано на рисунке 8.4.

Рис. 8.4. Поле давления при движении шара в невязкой жидкости Из формулы (8.21) и рис. 9.4 видно, что функция p симметрична относительно центра шара, а значит интеграл от давления по поверхности шара равен нулю. Это значит, что на шар, плывущий в идеальной жидкости не действуют силы сопротивления. Такая не соответствующая жизненному опыту и научным измерениям ситуация является следствием модели идеальной жидкости и называется парадоксом Даламбера. Этот - 82 парадокс относится к движению любого тела в идеальной жидкости, а не только шара.

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

Рассмотрим применение этого приёма на задаче о движении шара в бесконечности на него набегает поток со скоростью u = V, v = 0, w = 0.

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

Решение этой задачи в силу ее линейности можно найти, используя суперпозицию. Суперпозиция – это когда решение одной задачи представлено в виде суммы других, более легких задач. Справедлива только для линейных задач. Это значит, что на решение предыдущей задачи равномерное движение потока со скоростью u = V. Получим решение задачи обтекания в следующем виде:

Проверим, удовлетворяет ли этот потенциал граничным условиям непротекания на поверхности шара. Подставим r = a в (8.23) и получим Таким образом, граничные условия оказались верными, поверхность шара непроницаема.

Задачи к лекции Задача 8.1. Задано распределение скоростей для течения жидкости:

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

Задача 8.2. Движение идеальной жидкости задано уравнениями:

Является ли движение безвихревым? Если да, то найти потенциал скорости. Является ли данная среда несжимаемой?

Задача 8.3. Дано поле скоростей плоского установившегося потока жидкости Найти: 1) потенциал скорости, 2) уравнение эквипотенциалей, 3) divV.

линии постоянных.

ПОЛЗУЧИЕ ВЯЗКИЕ ТЕЧЕНИЯ

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

Приведем несколько примеров вязких ползучих течений.

1. Движение мелких частиц суспензии.

2. Туманы (мелкие капли жидкости в газе).

многокилометровые массы льда движутся очень медленно).

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

которые выводятся из общих уравнений Навье – Стокса Для ползучих течений ускорения настолько малы, что ими пренебрегают. Это означает, что силы давления и массовые силы уравновешиваются силами вязкости. В уравнениях (9.1) остаются лишь соответствующие члены, Необходимо добавить уравнение неразрывности (считаем плотность постоянной) - 86 Эта система уравнений называется уравнениями Стокса.

Давление в ползущих течениях является гармонической функцией. В самом деле, применим оператор дивергенции к векторному уравнению (9.2). Получим Так как для несжимаемой жидкости в силу (9.3) div v =0, то (div v ) = 0.

Если к тому же div g = 0 (массовые силы соленоидальны), то Простым примером массовых сил, удовлетворяющих условию div g = 0, является сила тяжести. Eё ускорение – это вектор g =(0,0,-g). Очевидно, что Рассмотрим два примера ползучих течений: первый из них относится внешним, второй – к внутренним течениям.

Пример 1. Течение Стокса около шара (1815 г).

Схема обтекания и обозначения показаны на рисунке 9.1.

Рис. 9.1. Схема обтекания шара и основные обозначения тяжести не учитывается. Радиус шара a, cos = x / a, r = x 2 + y 2 + z 2.

Поле скоростей описывается следующими формулами.

При r = a получаем скорости на поверхности шара. Легко видеть, что u=0, v=0, w=0, это граничное условие прилипания.

Течение в плоскости z = 0 показано на рис. 9.2 с помощью стрелок.

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

Рис. 9.2. Течение вязкой жидкости около шара. Поле скорости Давление в жидкости равно - 88 Карта изобар показана на рисунке.

Рис. 9.3. Течение вязкой жидкости около шара. Поле давления Перед лобовой частью шара давление резко возрастает, затем оно падает и в т.н. донной точке, лежащий напротив лобовой точки, достигает минимума. Вдали от тела как вверх, так и вниз по потоку давление стремится к p. Следовательно, парадокса Даламбера нет, и сила сопротивления будет равна интегралу давления (9.6) плюс интеграл касательных напряжений по поверхности шара. Сложные вычисления дают величину силы сопротивления Здесь d = 2a – диаметр шара, безразмерный комплекс Re d = u d / – число Рейнольдса, С x – коэффициент сопротивления (безразмерная величина). Формула (9.7) говорит о том, что сила сопротивления пропорциональна площади сечения шара S = а 2 и скоростному напору обратно пропорционален числу Рейнольдса.

Число Рейнольдса Re d играет фундаментальную роль в механике вязкой жидкости, все характеристики течения рассматриваются как функции числа Red.

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

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

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

скалярной функции f ( r, z, ) определяются как Интегрирование функции H равен Оператор Лапласа в цилиндрической системе координат имеет вид Компоненты вектора ускорения в координатах r, z, записывается следующим образом.

Обратим внимание на наличие двух дополнительных (по сравнению с прямоугольной декартовой системой координат) слагаемых для компонент в ar и a z. Первое из них, u / r - это центростремительное объясняется, в частности, смещение русел рек вправо по направлению течения.

Итак, уравнения Навье – Стокса и уравнение неразрывности в цилиндрической системе координат имеют следующий вид:

При расчете течения в трубе мы ограничимся рассмотрением простейшего симметричного по стационарного течения. В этом случае наша система существенно упрощается, поскольку В результате получаем два уравнения импульса Кроме того, так же как и при течении в плоском канале, отсутствует поперечная компонента течения (в направлении радиуса), следовательно выпадает всё первое уравнение, а второе уравнение импульса и уравнение неразрывности приобретают вид Последнее равенство означает, что скорость в любой точке сечения трубы (z=const) одинакова, то есть продольная скорость uz зависит только от r.

Окончательно получаем единственное уравнение или Это уравнение необходимо дополнить граничными условиями симметрии на оси трубы и прилипания на стенке:

Математическую задачу на основе соотношений (9.8), (9.9) можно сформулировать двумя способами.

1-й способ. Считаем заданным перепад давления, тогда краевая задача (9.8), (9.9) однозначно определяет все характеристики течения, а именно:

- среднюю скорость течения в трубе - напряжение трения на стенке 2-й способ. Задан расход через трубу Q = U R 2. В этом случае градиент подлежит определению.

Решим вначале первую задачу. Для этого удобно перейти к безразмерным переменным Уравнение (9.9) в переменных (9.12) имеет вид Поделив это уравнение на, получаем формулировку задачи в безразмерных переменных Безразмерный комплекс Re -число Рейнольдса - является единственным параметром задачи. При этом исходная задача определяется не четырьмя независимыми размерными параметрами (,, R, U), а единственным безразмерным параметром Re. Правая часть уравнения (9.13) в данной постановке считается заданной константой; обозначим ее через Заметим, что при течении в направлении z давление падает, dp / dz < 0, поэтому константа C отрицательна.

проинтегрируем по y от нуля до y. С учетом первого из граничных условий (9.13) получим Сократим на y и интегрируем ещё раз по y от 1 до y. Будем иметь Используя граничное условие прилипания на стенке, окончательно получим параболический профиль скорости в трубе - 94 C Главным вопросом при изучении течения по трубам является напряжения на стенке (9.11) к кинетической энергии потока, При записи формулы (9.16) использовалось решение (9.15) для эпюры скорости и формулы (9.12), связывающие размерные и безразмерные величины. Итак, коэффициент сопротивления трубы падает с ростом числа Рейнольдса по закону Первая задача о течении в трубе с заданным градиентом давления решена.

Перейдем к решению второй задачи, в которой задан расход Q = US, или, что равнозначно, средняя скорость U = Q / S. Для того чтобы определить градиент давления, необходимый для прокачки жидкости с такой скоростью, воспользуемся формулой (9.10), которую запишем в безразмерной форме Подставив сюда решение для скорости (9.15), найдем константу C и коэффициент сопротивления а затем и градиент давления Полученное решение называется решением Хагена – Пуазейля для гладкой трубы.

- 96 Лекция 10.

УРАВНЕНИЯ ПОГРАНИЧНОГО СЛОЯ

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

Это используется для расчета обтекания тел по следующим причинам.

1) Потенциальное течение описывается простыми уравнениями типа интегралов Коши – Лагранжа или Бернулли.

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

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

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

Рис. 10.1. Ламинарный пограничный слой при продольном обтекании Выведем уравнения пограничного слоя из двумерных уравнений Навье-Стокса Введем масштабы длины в продольном направлении l, длины в поперечном направлении ; продольной скорости u0 ; поперечной скорости v0 ; давления p0 ; времени t0. Среди этих масштабов u0 выбираем как среднюю скорость внешнего течения u0 = ue, l ассоциируем с продольной координатой x, а значения, v0, p0, t0 пока не известны.

Перепишем систему (10.1) в безразмерных переменных Из уравнения неразрывности получаем - 98 откуда находим масштаб поперечной скорости С учетом этого соотношения преобразуем уравнение импульса в направлении x.

Разделим это уравнение на u0 / l и обозначим через Re безразмерный комплекс u0l /. Будем иметь Выбор масштабов времени t0 = l / u0 и давления p0 = u0 теперь очевиден, причем последний вполне согласуется с давлением вне пограничного слоя, где течение происходит по уравнению Бернулли Сущность пограничного слоя заключается в том, что в нем инерция потока в продольном направлении x компенсируется сдвиговыми вязкими напряжениями, пропорциональными градиенту скорости в поперечном направлении y. Это означает, что в уравнении (10.4) необходимо положить коэффициент при 2U / Y 2 равным единице. Тогда Соотношение (10.6) показывает, что при больших числах Рейнольдса пограничный слой является тонким, 1, причем его толщина растет как квадратный корень от расстояния в продольном направлении (масштаб l можно заменить на x ), и тем больше, чем выше вязкость жидкости и ниже ее скорость.

движения (10.4) в пограничном слое принимает вид Преобразуем теперь второе уравнение импульсов.

Поделим это уравнение на ue /. Будем иметь Таким образом, из второго уравнения движения следует, что с точностью до величин порядка O ( 2 ) получаем Это значит, что давление поперек пограничного слоя не меняется и равно давлению во внешнем потоке, где p подчиняется уравнению Бернулли (10.5). Это позволяет исключить давление из рассмотрения и записать уравнения плоского пограничного слоя в которых только две неизвестных: компоненты безразмерной скорости U и V, а скорость внешнего потока U e ( x ) 1 задана.

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

Вместо толщины динамического пограничного слоя мы определяем толщину теплового погранслоя T. Будем изучать то же обтекание пластины, но добавим условие, что внешний поток имеет температуру Te, координаты вводятся так же, как в (10.2), а безразмерная температура определяется следующим образом:

Уравнение (10.11) преобразуется к виду откуда по аналогии с (10.6) получаем Здесь T = T / l, число Пекле Pe = u0l / a. Сравнивая формулы (10.13) и (10.6), легко сопоставить толщины динамического и теплового пограничного слоев:

Число Прандтля Pr является теплофизическим свойством жидкости; в общем случае Pr зависит от температуры. Понятно, что при Pr = 1 и отсутствии градиента давления безразмерные уравнения для продольной скорости и температуры совпадают, а значит совпадают и их графики.

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

Мы рассмотрим три варианта.

1) Задана температура стенки. y = 0, T = Tw. В безразмерном виде это эквивалентно условию Y = 0, = 0.

2) Задана плотность теплового потока. y = 0, T / y = q. В безразмерном виде это условие имеет вид где число Стантона определяется формулой 3) Задан теплообмен между стенкой и потоком по закону Ньютона.

y = 0, T / y = (T Tw ). Здесь, Вт/м 2 К – коэффициент теплообмена, который может зависеть от температуры, а также чисел Рейнольдса и Прандтля. В безразмерном виде это условие имеет вид где число Нуссельта определяется формулой концентрирует в себе информацию о теплообменных процессах в погранслое; при использовании граничного условия третьего рода безразмерное решение задачи не обязано равняться нулю на стенке (для размерной температуры это означает T (0) Tw ).

- 102 Задачи к лекции Задача 10.1. Записать уравнения плоского пограничного слоя для установившегося потенциального течения, выразив компоненты скорости через функцию тока.

Задача 10.2. Записать уравнения пограничного слоя, полученные в задаче 10.1, и граничные условия для случая пристенного пограничного слоя на продольно обтекаемой пластинке, считая профиль скоростей набегающего на пластинку потока невозмущенным. Предполагать, что ось OX направлена вдоль поверхности пластинки, ось OY по перпендикуляру к поверхности, а начало координат выбрать в передней кромке пластинки (задача Блазиуса).

МАТЕМАТИЧЕСКОЕ ОПИСАНИЕ ТУРБУЛЕНТНЫХ ТЕЧЕНИЙ

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

рис. 11.1 и 11.2). Переход от ламинарного течения к турбулентному определяется критическим числом Рейнольдса наблюдается ламинарное, а при Re > Re* - турбулентное течение. Для течения в трубе Re* 2300.

Рис. 11.1. Безотрывное стационарное ламинарное обтекание шара при Re=8.15.

- 104 Рис. 11.2. Обтекание шара при Re=15 000. Ламинарный пограничный слой отрывается перед экватором, но остается ламинарным на длине, почти равной радиусу. Затем слой становится неустойчивым и превращается в турбулентный.

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

Уравнения Рейнольдса Замеры реального поведения характеристик турбулентного потока дают следующую картину, рис. 11.3, которая позволила О. Рейнольдсу представить мгновенную характеристику как сумму среднего u и пульсации u.

Представление вида u = u + u тем более продуктивно, что в инженерных приложениях интерес представляют именно средние поля.

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

Проведем осреднение по Рейнольдсу уравнений Навье-Стокса для несжимаемой жидкости Подставим в него разложение ui = ui + ui, p = p + p. Будем иметь Применим к (11.3) оператор осреднения. Из-за первого и третьего свойств осреднения (11.2) все линейные вхождения пульсаций после осреднения пропадут. Останется Используя последнее уравнение (11.4), можем записать Произведение пульсаций называют их корреляцией второго порядка, а величины которые образуют тензором напряжений Рейнольдса. Уравнения (11.4), переписанные в виде – уравнениями Рейнольдса (черта над средними опущена). Сравнивая (11.5) с исходной системой уравнений Навье-Стокса, легко заметить, что в осредненных уравнениях появились дополнительные градиенты напряжений. Очевидно, модель (11.5) теперь не замкнута, и требуется дополнительно 6 уравнений для компонент симметричного тензора турбулентных напряжений. Простейший способ замыкания уравнений Рейнольдса состоит в том, чтобы выразить ij через характеристики осредненного течения, например Здесь коэффициент T называется турбулентной вязкостью, который в отличие от молекулярной вязкости не является свойством жидкости, в турбулентной вязкости через средние характеристики потока составляет предмет т.н. полуэмпирических моделей турбулентности.

Тем не менее, выражение (11.6) можно подставить в уравнения (3.5).

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

где l – длина пути смешения, e - тензор скоростей деформаций, 0.4 – постоянная Прандтля, y – расстояние до стенки.

Осредненное уравнение конвективной теплопроводности Температурное поле в турбулентном потоке тоже испытывает пульсации:

T = T + T. Если подставить это разложение в линейное уравнение теплопроводности, получим Раскроем скобки в конвективном члене.

При осреднении (11.2) уравнения (11.9) слагаемые с линейными вхождениями пульсаций, как и при выводе уравнений Рейнольдса, обратятся в ноль. В результате получим Используя уравнение неразрывности для пульсаций скорости, можем представить нелинейный член справа в виде дивергенции турбулентного потока тепла ( uiT ) / xi. В отличие от молекулярного потока Фурье qi = T / xi этот поток имеет конвективную природу: пульсация температуры T переносится пульсацией скорости ui.

Если умножить уравнение (11.10) на c и вернуться от оператора Лапласа к дивергенции теплового потока q, получим осредненное уравнение конвективной теплопроводности при турбулентном течении По аналогии с тем, как был введен коэффициент турбулентной вязкости в уравнениях Рейнольдса, можно ввести коэффициент турбулентной теплопроводности T выразив пульсационный поток qiT = c uiT через градиент среднего температурного поля:

Подставив (11.12) в (11.10), опуская черту над средними и вновь деля на c, окончательно получим уравнение для средней температуры Здесь aT – турбулентная температуропроводность. Обычно используют турбулентное число Прандля, чтобы выразить aT через вихревую вязкость В простейших сдвиговых течениях значение PrT близко к единице (равно 0.5 0.9).

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

Модель Прандтля применяется обычно применяется к простым потокам, когда средняя скорость имеет одну компоненту (пограничные слои, трубы, каналы). Пусть это скорость u в направлении x, а поперечная координата y, отсчитываемая от стенки.

Формула Прандтля для турбулентной вязкости имеет вид в которой l – длина пути смешения.

Длина пути смешения l в разных течениях задается по-разному.

Например, в свободном сдвиговом слое шириной d полагают l = Cd, C 0.1. Для течения жидкости в канале или трубе l = Cy 1 y / d ( d – глубина открытого канала, полуширина закрытого канала или радиус трубы).

- 110 Для турбулентного пограничного слоя (ТПС) экспериментально найдена зависимость l = y, 0.41 – постоянная Прандтля. Более коэффициента Ван-Дриста l = y [1 exp( y / A)].

Универсальный профиль скорости в ТПС получается из уравнения с нулевым градиентом давления. Мы учитываем, что пограничный слой имеет составную структуру: тонкий ламинарный подслой 0 < y < l и собственно турбулентный слой. Уравнение импульса в ТПС имеет вид Интегрируя это уравнение по y от l до y и обозначая через w = ( l ) напряжение сдвига на границе вязкого подслоя, имеем Если подставить сюда выражение (11.15) при l = y, получим Величина u = w / называется динамической скоростью. Обозначим переменные u + = u / u, y + = y / y0 и перепишем (11.18) в этих переменных.

Будем иметь du + / dy + = 1/ ( y + ), откуда следует т.н. логарифмический профиль скорости в ТПС Константа интегрирования u0, согласно экспериментальным данным, равна 5.0 5.2. Чтобы придать безразмерной координате y + реальный смысл, надо определить масштаб y0. Это можно сделать, рассмотрев течение в вязком ламинарном подслое и приняв, что касательное напряжение в нем постоянно и равно w, а скорость на границе имеет молекулярная кинематическая вязкость, находим Из определения следует, что y0 имеет порядок толщины ламинарного подслоя l. Понятно, что в действительности переход от вязкого подслоя к ТПС происходит плавно, и полученные результаты верны на некотором удалении от переходной зоны. Экспериментально установлено, что 50 < y + < 0.2, где – толщина ТПС.

Для внешней области (т.н. слой следа) 0.2 < y + < характерно гораздо более медленное изменение гидродинамических параметров. В качестве масштаба скорости в этой области принято использовать скорость на внешней границе пограничного слоя ue, а в качестве линейного масштаба – толщину (радиус трубы). В ряде моделей турбулентная вязкость во внешней области ТПС предполагается постоянной величиной. Клаузер предложил для течений в следе формулу T = ue, где = 0.0168 – эмпирический коэффициент.

На рис. 11.4 показан типичный профиль скорости в ТПС, развивающемся на плоской пластине без градиента давления.

Рис. 11.4. Три участка разбиения профиля: вязкий подслой, логарифмический слой и внешний слой (слой следа) внутренней части слоя (ТПС + ламинарный подслой) и в предположении (подтверждаемом практикой), что суммарное касательное напряжение постоянно поперек слоя и равно w = u2. Получим уравнение с граничным условием прилипания на стенке Перейдем к безразмерным величинам y + = u y /, u + = u / u и получим квадратное уравнение относительно производной z = du + / dy +, а именно z + ( y + ) z 2 = 1. Положительное решение этого уравнения имеет вид z = 1 + 1 + 4( y + ) 2 / ( 2( y + ) 2 ). Умножив числитель и знаменатель на 1 + 1 + 4( y + ) 2, получим задачу Коши для скорости Эту задача решается численно. Решение задачи показано на рисунке 11.5 в логарифмических координатах. Представлены два случая: решение уравнения (11.21), в котором длина пути смешения определена согласно демпфирующим коэффициентом Ван-Дриста Заметно существенное отличие профилей скорости.

Рис. 11.5. Универсальный профиль средней скорости во внутренней части турбулентного пограничного слоя при определении длины пути смешения по Фон Карману (внизу) и Ван-Дристу (вверху) В заключение отметим, что в настоящее время для расчета турбулентных течений используются более сложные дифференциальные модели турбулентной вязкости.

- 114 Задачи к лекции Задача 11.1. Какой должна быть скорость и кинематическая вязкость жидкости при моделировании движения судна при условии, что длина модели в 100 раз меньше длины судна. В качестве критериев подобия использовать постоянство чисел Рейнольдса Re = и Фруда Fr =.

Задача 11.2. Конденсатор паровой турбины состоит из 250 трубок для протока охлаждающей воды. Определить максимально допустимый диаметр трубок, обеспечивающий турбулентный режим течения при общем объемном расходе воды через все трубки Коэффициент кинематической вязкости принять равным Предполагается, что устойчивый турбулентный режим течения воды в трубках удерживается при Re > 104, то есть при числах Рейнольдса, превышающих критическое значение.

Ответ: d 14 мм.

Решение задачи 1.1. Напряжение в жидкости может быть рассчитано из уравнения =. В случае линейного профиля продольной скорости по высоте зависимость u ( y ) можно представить в виде u = u0 y h, поскольку на твердой поверхности скорость равна нулю. Здесь u0 – величина продольной скорости на высоте h от твердой поверхности.

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

Решение задачи 1.3.

Решение задачи 2.2. Расчет произведений векторов и тензоров основывается на механизме перемножения матриц.

Пусть a D = v. Тогда - 116 Пусть D b = w. Тогда Решение задачи 2.3. Изменение содержания вещества сплошной среды за единицу времени в объеме V есть разность между втекшим и вытекшим через ограничивающую объем поверхность количеством использованием формулы Остроградского-Гаусса:

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

Таким образом, для заданного поля скоростей в любой части полосы, ограниченной плоскостями x = 1 и x = 2, изменения количества вещества не происходит.

Решение задачи 3.2. Компоненты напряжения в точке P принимают значения Единичный вектор нормали в точке P определяется вектором Следовательно, в точке P Тогда единичный вектор нормали в точке P есть n = j 2 + 3 k 2. Вектор напряжения на площадке, перпендикулярной к n в точке P, равен Решение задачи 4.1. Для определения функции давления следует записать уравнения Навье-Стокса для соответствующей модели сплошной среды Последнее уравнение, являющееся выражением закона неразрывности для несжимаемой жидкости выполняется автоматически для заданного поля скоростей. Два первых уравнения перепишутся в виде откуда можно записать выражение полного дифференциала давления Из последнего выражения найдем функцию давления Решение задачи 5.1. Уравнение неразрывности для несжимаемой жидкости имеет вид div u = 0 или ui xi = 0. В данном случае Складывая, получаем:

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

Для безвихревого течения необходимо выполнение условия rot u = 0.

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

Решение задачи 5.4. Из определений функции тока и линий тока следует, что линии тока в данный момент времени t являются решениями дифференциальных уравнений:

В случае плоских течений остается одно дифференциальное уравнение:

которое для данной задачи приводится к виду:

Решением данного уравнения является зависимость x y = C – уравнение гиперболы. Различные значения постоянной C будут соответствовать различным линиям тока (рис. 2). Направление течения можно определить по заданным скоростям. Предполагая K положительным параметром, получим, что при положительных значениях y компонента v = K y положительной. Направления течения в других четвертях координатной плоскости определяются аналогично (см. рис).

Рис. к задаче 5.4. Изображение линий тока и направления течения Решение задачи 6.3. Вес конуса в Вес воды в объеме, равном объему погруженной в нее части конуса Fg B = r 2 h B g. Из рисунка видно, что После подстановки rK в формулу для Fg K и учитывая, что по закону Архимеда Fg K = Fg B, получаем - 120 Д Решение задачи 7.4. Пьезометрические высоты в сечениях (1) и (3) будут Для определения пьезометрической высоты в сечении (2) запишем уравнение Бернулли для первого и второго сечений:

Здесь z1 = z2, так как труба горизонтальна. Скорости Тогда искомая пьезометрическая высота во втором сечении будет:

К решению задачи 7.6. Центробежная сила: F = m R 2. Уравнение установившегося движения: P x i = Fi. Уравнение поверхности найдем интегрирования – из условия p ( 0,0, z0 = 0 ), z0 - из объема жидкости z0 = h ( 2 R 2 / ( 4 g ) ). Ответ: P = g h R 2 - совпадает с весом жидкости.

Решение задачи 10.1. Уравнения плоского пограничного слоя для установившегося течения имеют вид V выражаются через функцию тока следующим образом:

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

Решение задачи 10.2. Заметим, что в данном случае внешний по отношению к пограничному слою поток будет однородным, с одинаковой U e = const, dU e dx = 0 и уравнения пограничного слоя примут вид:

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

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

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

Записанные уравнения пограничного слоя и граничные условия являются формулировкой задачи Блазиуса.

- 122 Контрольные задания сферическая, а – прямоугольная декартова СК. Найти якобиан преобразования xi x j.

цилиндрическая, а { xi } – прямоугольная декартова СК. Найти якобиан преобразования xi x j.

где r 2 = x 2 + y 2 + z 2 и A – произвольная константа, удовлетворяет уравнению неразрывности несжимаемой жидкости.

4. Определить линии тока и описать движение среды по полю скоростей u = K y, v = K x.

5. Доказать, что для течения ui = xi (1 + t ) линии тока и траектории совпадают.

линии тока будут окружностями.

7. Доказать, что течение с полем скоростей безвихревым?

8. Найти разность высот h (рис 3), если H = 40 см, Ответ: 100 см образовавшегося круглого отверстия равен 0.8 м. Определить силу избыточного давления морской воды ( = 1030 кг / м3 ) на «пластырь», заведенный на пробоину. Ответ: 15.2 кН.

10. На переднюю часть погруженного в реку тела, находящегося на глубине 5 м, действует максимальное избыточное давление, равное 0.815*10^5 Па. Определить скорость течения реки.

11. Определить число Рейнольдса для движения воды в трубе диаметром d = 10 см с объемным расходом Q = 10 л / с. Коэффициент кинематической вязкости воды = 1.5 мм 2 / с. Ответ: 84666.

12. Определить число Рейнольдса для течения моторного топлива по трубопроводу диаметром 100 мм при расходе 10 м3 / ч. Кинематическая вязкость топлива = 35 106 м 2 / с. Ответ: Re=1011.

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

14. Бочка заполнена смазочным маслом и плотно закрыта при температуре t1 = 10°C и давлении p1 = 105 Па. Определить давление в заполненной бочке = const ).

15. При гидравлических испытаниях трубопровода диаметром и длиной l = 4 м в системе поднимают давление. Определить D = 0.5 м массу воды, которую требуется дополнительно закачать в трубопровод, чтобы давление повысилось на величину p = 2.5 бар.

Ответ: 0.0883 кг.

- 124 Чему равно относительное изменение плотности морской воды (в процентах) при вертикальном погружении на глубину h = 300 м ?

Ответ: 0.13 %.

17. Записать уравнение неразрывности для несжимаемой жидкости в сферической системе координат.

18. Записать уравнение неразрывности для сжимаемой жидкости в цилиндрической системе координат.

1. Седов Л.И. Механика сплошной среды. Т.1,2. М.: Наука. 1978.

2. Шлихтинг Г. Теория пограничного слоя. М.: Изд-во иностранной литературы. 1956. 528 с.

3. Себиси Т., Брэдшоу П. Конвективный теплообмен. М.: Мир. 1987.

4. Беннет К.О., Майерс Дж. Е. Гидродинамика, теплообмен и массообмен. М.: Недра. 1966. 725 с.

5. Ван-Дайк М. Альбом течений жидкости и газа. М.: Мир. 1986. 182 с.

6. Фрик П.Г. Турбулентность: подходы и модели. Москва-Ижевск:

Институт компьютерных исследований. 2003. 292 с.

7. Мазо А.Б. Моделирование турбулентных течений несжимаемой университета. 2007. 109 с.

8. Лойцянский Л.Г. Механика жидкости и газа. М.: «Наука», 1987, – 9. Мейз Дж. Теория и задачи механики сплошных сред. М.: «Мир», гидромеханики. М.: «РКонсульт», 2003 151 с.

- 126 -