Основы гидрогазодинамики

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

В потенциальном потоке , в плоском течении

                          функция тока ψ гармоническая

Сравнение потенциала φ и ψ позволяет записать:

 -

условие Коши-Римана.

15. Комплексный потенциал, комплексная скорость

Из теории комплексной переменной известно, что если две функции φ и ψ, зависящие от х и у, удовлетворяют условиям Коши-Римана, то комплексная величина будет не просто зависеть, а являться функцией от комплексной переменной , то есть существует некоторая функция , действительной частью которой является φ, а мнимой ψ. .

Функция  имеет большое значение при изучении плоских потенциальных течений и называется комплексным потенциалом или характеристической функцией течения.

Так как  является аналитической функцией от , то ее производная не зависит от направления дифференцирования, а зависит только от положения точки в пространстве, то есть

по условию Коши-Римана:

Если вектор U разложить в комплексной плоскости годографа U, то .

Производная от комплексного потенциала дает зеркальное изображение комплексной U относительно действительной оси. Обозначим ее как

.

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

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

16. Частные случаи плоских потенциальных течений

1. Плоско параллельный поток:

Рассмотрим комплексный потенциал - , где а – действительное число

 и

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

Линии тока - семейство прямых, параллельных оси х.  - уравнение эквипотенциальных поверхностей.

Для построения поля скоростей возьмем производные

;

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

2. Источник и сток.

Рассмотрим комплексный потенциал , а – действительное число (), тогда

Уравнение для потенциала: .  - эквипотенциальные линии, семейство окружностей с центром в точке (0,0).

 - уравнение функций тока.  - семейство прямых, проходящих через точку (0,0).

Характер (вид) течения определяет знак при а. Если a>0, то это источник, если a<0, то это – сток.

 - объемный расход;

;  

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

Если их свести вместе, то получится диполь.

3. Рассмотрим комплексный потенциал:

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

Уравнение для линий тока  - семейство окружностей, проходящих через точку (0,0) с центрами на оси у.

4. Рассмотрим комплексный потенциал вида:

Г – циркуляция вектора скорости – круговое течение потока.

- семейство прямых, проходящих через точку (0,0).

Это уравнение эквипотенциальных линий.

 - функция тока;

- линии тока – семейство окружностей с центром в (0,0).

 - радиальная скорость;

Исследованный потенциал определяет течение, которое называется потенциальным вихрем.

Окружная скорость изменяется по гиперболе.

17. Безциркуляционное обтекание круглого цилиндра

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

Если приравнять  к константе получим уравнение эквипотенциальной линии.  - линии тока,  - уравнение для нулевой линии тока. Если принять , то получим уравнение для нулевой линии тока:

Оно разделится на два: 1) у=0;

2)  - окружность с радиусом

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

 ; ;

Рассмотри составляющие скорости:

Значит: , то есть окружная составляющая скорости изменяется по синусоиде (при ,  - ). Точки А и В передняя и задняя критические точки соответственно.

Максимальные значения окружной скорости  при 90˚ и 270˚ - точки С и Д.

Нулевая линии тока проходит из (-∞) в передней критической точке А, раздваивается огибает цилиндр, соединяется в задней критической точке В и уходит в (+∞).

Для определения распределения давления по поверхности воспользуемся уравнением Бернулли:

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

На поверхности существует только окружная скорость, следовательно, для поверхности:

Из полученной формулы следует, что давление на поверхности максимально в критических точках А и В () и минимально в точках С и Д ().

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

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

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

Распределение давления описывается экспериментальными линиями, которые отличаются от теоретических. С увеличением скорости распределение давления стремится как бы к теоретическому, и

18. Обобщенный закон Ньютона

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

;

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

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

Касательное трение при движении потока вдоль оси х может быть записано в виде:

При движении потока вдоль оси у:

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

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

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

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

Компоненты, учитывающие вязкость связаны с диагональными компонентами тензора скоростей деформации соотношениями:

складываем

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

обобщенный закон Ньютона для нормальных напряжений

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

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

Рассмотрим изотермическое движение вязкой несжимаемой жидкости:

, .

В этом случае нормальные напряжения примут вид:

Уравнения движения получим из уравнения движения среды в напряжениях:

- система уравнений движения вязкой несжимаемой жидкости

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

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

1)                геометрические условия – задается геометрия изучаемой системы (канала и т.д.)

2)                физические условия однозначности – задается вид движения жидкости и значения ее основных параметров .

3)                граничные условия – определяют условия течения на границе рассматриваемой системы. Часто в качестве граничных условий используют условия прилипания потока, т.е. скорость потока на поверхности равна 0. Скорость набегающего потока задается  или среднемассовая скорость  или скорость на границе пограничного слоя.

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

20. Подобие гидродинамических явлений

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

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

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

1) ; 2) по физическим условиям однозначности должны быть определены ; 3) по граничным условиям должна быть определена ; 4) по начальным условиям задано характерное время , например период определяющий темп внешних воздействий

Таким образом в уравнениях зависимые переменные определяются как функции независимых переменных x, y, z, t и параметров задающих условие однозначности

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

; ;

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

; ; ;

Выразим размерные величины через их масштабы в Уравнении Навье-Стокса:

Аналогично могут быть получены составляющие системы уравнений вдоль оси у и z.

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

 - критерий динамической гомохронности;

 - критерий Фруда;

 - число Эйлера;  - критерий Рейнольдса.

После приведения уравнений к безразмерному виду изменился их физический смысл, т.к. один и тот же вид уравнений с подобными условиями будут соответствовать не единственному условию, а целой группе подобных явлений. В соответствии с теоремой подобия Кирпичева-Гухмана гидродинамические явления будут подобными если они: 1) описываются одной системой дифуравнений; 2) имеют подобные условия однозначности; 3) имеют численно равные критерии подобия

21. Критериальные уравнения. Критерии и числа подобия

После приведения уравнения Навье-Стокса к следующему виду они стали содержать следующие типы переменных: 1) безразмерные независимые переменные ; 2) безразмерные зависимые переменные ; 3) безразмерные критерии – комплексы, состоящие из величин заданных по условиям однозначности .

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

- система обобщенных или критериальных уравнений

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

.

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

1) отношение переменной к одноименной величине, заданной по условию однозначности:

2) если по условию однозначности нельзя задать одноименную величину, то строится комплекс приводящий величину к безразмерному виду – число подобия:

В числа подобия входят определяемая величина. Критерий подобия состоит из заранее известных величин, заданных по условиям однозначности..

1) Критерий Рейнольдса  - определяет соотношение сил инерции и вязкости в однородном потоке. Это важнейший гидродинамический критерий для вынужденного движения. При движении потока в нем возникают возмущения, которые исходят от стенок канала или вносятся в поток извне. Влияние возмущений зависит от соотношения сил. Если преобладают силы вязкости возмущения гаснут и поток не меняет своей структуры. Если преобладают силы инерции возмущения развиваются дальше, поток меняет течение, изменяется его структура. Граница соотношения сил определяется по значению Reкр. Если Re<Reкр преобладают вязкие силы, Re>Reкр – силы инерции. Re характеризует движение при соизмеримости инерции и вязкости. Если в потоке преобладает какой-то один вид сил характер перестает зависеть от величины Re. В этом случае говорят, что течение автомодельно относительно критерия Re.

2) Критерий гидродинамической гомохронности  - определяет соотношение между периодом темпа внешних воздействий на поток и периодом перестройки скоростного поля. Используют только для нестационарных задач. - время, за которое проходит частица, движущаяся со скоростью V0, путь l0. Если в задаче время подлежит определению, то рассматривается не критерий, а число Струхала:

3) Критерий Фруда  - определяет соотношение между силами инерции и тяжести в потоке. Используется только в задачах, в которых гравитационные эффекты имеют важное значение. Однако в таких задачах часто сложно задать характерную скорость (при естественной конвекции), поэтому строится критерий, в котором исключается скорость:

 - критерий Галилея.

При гравитационном движении важное значение имеет параметрический критерий: .

Причем ρ и ρ0 – плотности не только в разных точках, но и в различных фазах.  - критерий Архимеда.

При гравитационном течении однофазной жидкости движение возникает в результате расширения:

 - коэффициент объемного расширения.

 - критерий Гросгофа.

4) Число Эйлера  - определяет соотношение сил давления и сил инерции; определяемая величина; т.к. часто давление в потоке неизвестно, то больший интерес представляет определение перепада давления на рассматриваемом участке .

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

22. Моделирование ГГД явлений

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

Свяжем одноименные величины с помощью констант подобия:

 - условие, определяющее выбор констант подобия

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

23. Ламинарное и турбулентное движение

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

 - общее время наблюдения за какой-то точкой

 - время существования режима в этой точке

При  - ламинарное

При  - турбулентное

 - переход

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

24. Пограничный слой и его характерные толщины

При обтекании любого тела потоком реальной жидкости поток как бы «прилипает» к поверхности. По мере удаления от поверхности скорость возрастает и, начиная с некоторого расстояния, скорость равна скорости набегающего или невозмущенного потока. В этом состоит проявление вязкости жидкости. Прандтль определил, что толщина слоя в котором проявляется вязкость увеличивается по мере продвижения потока. Чем меньше скорость набегающего потока, тем больше толщина слоя, в котором проявляется вязкость. Прандтль назвал эту часть жидкости у поверхности гидродинамическим пограничным слоем. Все течение он разбил на 3 части: 1 – гидродинамический пограничный слой – область, где сосредоточено влияние вязкости, 2 – след, 3 – невозмущенный поток. Во внешнем потоке вязкость можно не учитывать и считать жидкость идеальной, то есть без трения. Таким образом общая задача обтекания разбивается на 2 части: 1) течение жидкости в пограничном слое; 2) течение идеальной жидкости. Результаты решений должны совпадать на внешней границе пограничного слоя. Поскольку скорость в пограничном нарастает от 0 до скорости внешнего течения постепенно - асимптотически, поэтому δ определяется достаточно условно. Принято считать за δ такое значение у поперечной координаты, при котором скорость u отличается от U не более чем на 1-2%. u/U=0,98..0,99. Теория погранслоя использует и другие более точно определяемые толщины. Рассмотрим как влияет вязкость на кинематику (положение линий тока) и динамику (потерю количества движения).

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

Если бы жидкость была идеальной, то та же линия тока располагалась бы ближе к поверхности на величину δ:

Максимальное значение δ будет равняться: δ= δmax= δ*

 - толщина вытеснения.

Это более точно вычисляемая величина чем δ.

δ* может быть представлена, как линейный отрезок через который проходит секундный объемный расход идеальной жидкости, равный потере расхода через сечение погранслоя из-за торможения в реальном течении. Из-за торможения снижается и количество движения жидкости в погранслое. С количеством движения связана толщина потери импульса:

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

25. Переход ламинарного ПС в турбулентный

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

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

Если перемещать микротрубку вдоль обтекаемой поверхности на одинаковом расстоянии от нее, то по мере погружения микротрубки в ламинарный погранслой значение скорости, которое она будет показывать с увеличением х будет убывать. Вначале перехода ламинарного погранслоя в турбулентный (хкр1) профиль скорости начинает перестраиваться – становится более заполненным – прижимается к пластине, поэтому скорость возрастает. При x>хкр2 микротрубка погружается уже в турбулентный погранслой и поэтому скорость будет уменьшаться. За точку перехода принимают значение х, в которой скорость достигает минимума. Для определения границы перехода используют:

На положение хкр значительное влияние может оказать скорость изменения внешнего потока .

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

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

-                     возмущения, которые вносятся в погранслой из внешнего потока (внешний поток турбулизирован)

-                     возмущения, исходящие от поверхности (шероховатости, стыки поверхностей)

Рассмотрим влияние турбулентности внешнего потока. Если степень турбулентности внешнего потока меньше 0,1%, то возмущение не влияет:

Если , то с увеличением ε хкр1 и хкр2 уменьшаются. Влияние шероховатости приводит к более ранней смене режимов.

Рассмотрим влияние величины U. Чем больше U, тем меньше хкр.

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

 (1)

В турбулентном ядре распределение скорости может быть описано универсальным логарифмическим профилем распределения:

 (2)

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

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

При Re>Reкр на цилиндре сначала образуется ламинарный погранслой, который преобразуется в турбулентный, который срывается при угле 140°.

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

Страницы: 1, 2