Презентация на тему: Турбулентное движение

* Численные методы Автор курса лекций: Породнов Борис Трифонович, д. ф.-м. н., профессор кафедры молекулярной физики УГТУ-УПИ Екатеринбург 2008 решения задач механики сплошных сред 4. Турбулентность. Пограничный слой 900igr.net

* Модуль 5 ТУРБУЛЕНТНОСТЬ 2008. Численные методы…Лекция 19

ТЕОРИЯ ПРАНДТЛЯ. ТУРБУЛЕНТНОЕ ДВИЖЕНИЕ В ТРУБАХ Лекция 19 2008. Численные методы…Лекция 19

* Цели изучения: Теория Прандтля турбулентного движения: – определение профиля скорости, – понятие вязкого подслоя. Эксперименты по изучению свойств турбулентного движения: – формулы Блазиуса для коэффициента сопротивления гладких труб, – опыты Никурадзе с шероховатыми трубами. Развитая турбулентность: – масштабирование пульсаций, – энергия крупно- мелкомасштабных вихрей, – закон Колмогорова-Обухова. 2008. Численные методы…Лекция 19

* Содержание 10.5. Теория Прандтля: 10.5.1. Постановка задачи, 10.5.2. Средняя длина перемешивания. Профиль скорости, 10.5.3. Вязкий подслой. 10.6. Турбулентное движение в трубах: 10.6.1. Движение жидкости в трубах с гладкими стенками, 10.6.2. Зависимость коэффициента гидравлического сопротивления от числа Рейнольдса, 10.6.3. Эмпирический закон Никурадзе для коэффициента сопротивления. 10.7. Развитая турбулентность: 10.7.1. Масштабирование турбулентных пульсаций, 10.7.2. Диссипируемая энергия при крупномасштабном движении, 10.7.3. Закон Колмогорова-Обухова 2008. Численные методы…Лекция 19

* 10.5. Теория Прандтля 10.5.1. Постановка задачи Теория Прандтля турбулентного движения основывается на аналогии с молекулярным движением в газе. Рассмотрим плоское турбулентное движение несжимаемой среды вдоль бесконечной стенки (рис. 10.6). Выберем два слоя, параллельных оси x. Среднюю скорость движения частиц среды в них обозначим через u и u+du. Предположим, что скорость u определяется средней пульсационной скоростью , которая является функцией y, а средняя пульсационная у-составляющая скорости равна 0. Прандтль по аналогии с кинетической теорией газов предполагает, что расстояние между слоями - это некоторый аналог средней длины свободного пробега молекул в газе. На этом расстоянии частицы среды могут переходить из слоя в слой лишь благодаря своим пульсационным скоростям . Причем, импульсы таких частиц сохраняются при прохождении длины . Частица среды, попадающая из нижнего слоя в верхний, имеет среднюю скорость , отличающуюся от средней скорости верхнего слоя на du. Очевидно, "сталкиваясь" в верхнем слое с его частицами, скорость последних так же может изменяться на величину порядка du. Это изменение создает пульсацию порядка du: (10.5.1)

* 10.5.2. Средняя длина перемешивания. Профиль скорости Предположим, что турбулентность изотропна, по крайней мере, вдали от стенки: Тогда тензор турбулентных напряжений можно записать в виде: (10.5.2) Величину l называют длиной пути перемешивания, которая может зависеть от y. Далее Прандтль предполагает, что длина пути перемешивания прямо пропорциональна расстоянию от стенки: (10.5.3) Здесь к - некоторый безразмерный коэффициент пропорциональности. Предположение (10.5.2) можно подкрепить рассуждением о том, что стенка препятствует развитию поперечных пульсационных скоростей и, очевидно, на самой стенке длина пути перемешивания или пропорциональная ей "длина свободного пробега" должна быть равна нулю. Как и при ламинарном режиме течения, также предполагается, что тензор турбулентных напряжений не изменяется по нормали к стенке и равен его значению на стенке τ0, т.е Для определения средней скорости u, используя (10.5.2, 3), получаем уравнение Решение этого уравнения имеет вид: (10.5.4)

* 10.5.3. Вязкий подслой Из (10.5.4) следует, что средняя скорость турбулентного движения вязкой несжимаемой среды вдоль бесконечной стенки распределена по логарифмическому закону. Величина в (10.5.4) имеет размерность скорости и её обычно обозначают через . Формулу (10.5.4) удобно записать в безразмерной виде. Для этого составим безразмерную комбинацию из величин (y, ν, ). Такой комбинацией является выражение . Тогда формулу (10.5.4) можно записать для безразмерной скорости в виде: (10.5.5) Такой логарифмический профиль скорости в пограничном слое называют профилем Блазиуса Очевидно, что и при турбулентном движении среды на поверхности обтекаемого тела вследствие условий прилипания всегда существует тонкий слой, в котором среды движется ламинарно. Этот слой называют вязким подслоем. Ясно, что формула (10.5.5) справедлива лишь вне вязкого подслоя. Постоянная с должна быть определена экспериментально.

* 10.6. Турбулентное движение в трубах. 10.6.1. Движение жидкости в трубах с гладкими стенками Движение среды в трубах, как очень важное для практики, интенсивно изучалось экспериментально многими исследователями, начиная с Пуазейля (1840-42 гг) и Рейнольдса (1876-83 гг). Более ранние опыты создавали, как в настоящее время сказали бы, экспериментальную базу данных для развития науки. Исследователи зачастую представляли только опытные зависимости от некоторых отдельных параметров, характеризующих движение (диаметр трубы, плотность жидкости, скорость движения и т.п.), а не имея представления о подобии движений, не могли провести их обобщение. Возможность правильной постановки опытов и обобщения их результатов появилась лишь после введения в теорию критерия динамического подобия потоков - числа Рейнольдса. Шведский учёный Блазиус в 1911 г. в основу обобщения экспериментальных результатов ввел число Рейнольдса, что позволило ему эмпирически найти формулу для коэффициента сопротивления труб при турбулентном движении, пригодную для любых сред в широком диапазоне скорости движения среды и диаметров труб. В теме 9 было показано, что для расчета гидравлического сопротивления ΔР гладких, круглых цилиндрических труб используется формула вида (10.6.1) Здесь - средняя по сечению трубы скорость движения среды.

* Зависимость коэффициента гидравлического сопротивления от числа Рейнольдса Для ламинарного движения коэффициент λ гидравлического сопротивления равен: (10.6.2) Для турбулентного движения Блазиус получил формулу (10.6.3) Из этих формул видно, что сопротивление труб ΔР при ламинарном движении пропорционально первой степени средней скорости, тогда как при турбулентном движении пропорционально или примерно квадрату скорости. Предположим, что при турбулентном движении среды в трубе обеспечивается кинематическое подобие при изменении скорости и радиуса трубы. Это означает, что, если в какой-то точке трубы скорость изменилась в некоторое число раз, то во всех других точках она изменится во столько же раз. Если за характерную выбрать скорость на оси трубы , то это условие кинематического подобия можно записать следующим образом: (10.6.4) Из закона сопротивления трубы при турбулентном течении при некоторых предположениях можно вывести закон распределения скорости по сечению трубы. Можно показать, используя формулу (10.6.3) , что этот закон имеет вид: (10.6.5) Как видно из (10.6.5), при турбулентном движении средняя скорость движения среды по трубе мало отличается от максимальной в отличие от ламинарного движения, при котором .

* Профили скорости при ламинарном и турбулентном режимах течения Нa рис. 10.7 представлены примерные эпюры скоростей при ламинарном и турбулентном движении жидкости в трубе. Как было показано в теме 8, для ламинарного движения профиль скорости определяется формулой Как видно из рисунка 10.7, основное изменение скорости турбулентного движения наблюдается вблизи поверхности, вблизи же оси это изменение мало. Рис. 10.7 В последствии турбулентное движение жидкости в трубах изучалось многими исследователями. Для коэффициента сопротивления трубы в настоящее время существует более десятка различных эмпирических формул. Так, в 1932 г. Никурадзе, обработав свои экспериментальные результаты согласно теории турбулентности Прандтля, получил для распределения скоростей эмпирическую формулу: (10.6.6) Здесь у - расстояние от стенки трубы. Характерная скорость по экспериментально измеряемым величинам может быть определена следующим образом:

* Эмпирический закон Никурадзе для турбулентного движения в трубах с гладкими стенками Формула (10.6.6) хорошо описывает экспериментальные результаты во всех точках сечения трубы, кроме, конечно, тонкого ламинарного подслоя, примыкающего непосредственно к стенке. Используя формулу (10.6.6), можно получить уравнение для определения коэффициента λ гидравлического сопротивления трубы: (10.6.7) Формула Никурадзе (10.6.7) справедлива в более широком диапазоне изменения чисел Рейнольдса, но более сложна для практического использования, чем формула Блазиуса (10.6.3). Рис.10.8 График зависимости коэффициента сопротивления гладких труб от числа Рейнольдса изображен на рис.10.8. Как видно из рисунка, при переходе от ламинарного режима течения к турбулентному коэффициент сопротивления трубы возрастает скачком. Тщательные опыты Никурадзе показывают, что показатель степени в степенной формуле (10.6.5) для распределения скорости по сечению трубы непрерывно убывает от 1/7 при Re = 4·103 до 1/10 при Re = 3,2 106. Это соответствует тому, что сопротивление трубы λ с увеличением числа Рейнольдса приближается к квадратичному закону в зависимости от средней по сечению скорости движения.

* 10.6.2.Турбулентное движение в трубах с шероховатыми стенками Выше было рассмотрено турбулентное движение среды в трубах с гладкими стенками. Опыт показывает, что степень чистоты внутренней поверхности трубы существенным образом влияет на её сопротивление, причем тем сильнее, чем больше число Рейнольдса. Так, трубы после длительной эксплуатации (вследствие коррозии) увеличивают коэффициент сопротивления более, чем в два раза. Очевидно, формулы для распределения скоростей и для коэффициента сопротивления каким-то образом должны учитывать шероховатость стенок труб. Это тем более очевидно, поскольку шероховатость труб является геометрическим параметром, необходимым для динамического подобия потоков. Обозначим среднюю высоту бугорков шероховатостей на стенках через h. В опытах Никурадзе шероховатость создавалась искусственно при помощи приклеивания к стенкам трубы песчинок, имеющих некоторый средний размер зерен h. Очевидно, что трубы будут геометрически, а потоки - динамически подобными, если для них сохраняется отношение h/r0 . Следовательно, формулы для Рис. 10.9 распределения скоростей и коэффициента сопротивления шероховатых труб должны быть некоторыми функциями от относительной шероховатости . Обрабатывая многочисленные экспериментальные данные, Никурадзе получил результаты, которые качественно иллюстрируют зависимость на рис.10.9.

* Эмпирический закон Никурадзе для турбулентного движения в трубах с шероховатыми стенками Всю область чисел Рейнольдса можно разбить на три области: область I ламинарного движения, переходную область II и область III, в которой коэффициент гидравлического сопротивления λ не зависит от числа Рейнольдса. Как видно, трудно экспериментальные кривые описать некоторой единой формулой. Однако для больших чисел Рейнольдса (область III), когда коэффициент сопротивления не зависит от числа Рейнольдса, а зависит лишь от относительной шероховатости, Никурадзе приводит следующие эмпирические формулы: (10.6.8) Вышеприведенное условие применимости формул (10.6.8) на практике почти всегда выполняется. Конечно, можно и не пользоваться формулами (10.6.8), а для практических расчетов брать данные прямо с графика на рис.10.9. Сложный характер кривых на рис.10.9 можно объяснить физически. При ламинарном движении (область I), если шероховатости не загромождают сечения трубы, их влияние никак не сказывается на сопротивлении трубы. При этом, очевидно, бугорки шероховатостей погружены в неподвижную между ними среду, играющую роль смазки (рис.10.10а). Даже при турбулентном движении в основном потоке около стенки всегда существует ламинарный подслой, который прикрывает бугорки шероховатостей.

* Физические причины независимости коэффициента гидравлического сопротивления от числа Рейнольдса При дальнейшем увеличении числа Рейнольдса (область II) вершины бугорков начинают обнажаться и с них срываются вихри, что создает дополнительную турбулизацию потока, а, следовательно, увеличение коэффициента гидравлического сопротивления (рис. 10б). Рис. 10.10 Причем, это произойдет при тем больших числах Рейнольдса, чем меньше относительная шероховатость трубы h/r0 . Наконец, при больших числах Рейнольдса или скорости движения среды (область III) бугорки шероховатостей полностью обнажаются (рис.10.10в), и степень турбулизации потока полностью определяется вихрями, срывающимися с их вершин, и, следовательно, коэффициент сопротивления не зависит от числа Рейнольдса, характеризующего основное движение в трубе. Очевидно, что чем больше относительная шероховатость стенки h/r0 , тем при меньших числах Рейнольдса наблюдается отклонение коэффициента λ от рассчитанного по формуле Блазиуса для гладких труб, и тем большее его значение наблюдается в области III (см. рис. 10.9).

* 10.7. Развитая турбулентность. 10.7.1. Масштабирование турбулентных пульсаций Рассмотрим свойства турбулентного движения, пользуясь методом анализа размерностей физических величин. После потери устойчивости движения в среде сначала появляются крупные вихри, создающие крупные нерегулярные пульсации скорости: частицы среды двигаются по сравнительно плавным траекториям, хотя их траектории сильно перепутаны. По мере дальнейшего увеличения числа Рейнольдса траектории частиц становятся менее плавными, появляются более мелкие вихри. На фоне плавных траекторий появляются более мелкие их искривления, а скорость частиц приобретает более мелкие нерегулярные пульсации. При больших числах Рейнольдса турбулентное движение содержит вихри от самых крупных до самых мелких, в которых частицы среды перемещаются как целое с одинаковыми скоростями – такое движение называют развитым турбулентным движением. Очевидно, турбулентное движение можно рассматривать как результат наложения движений различных масштабов. Под масштабом движения понимают порядок тех расстояний, на которых скорость частиц среды изменяется существенно. При появляются сначала крупномасштабные пульсации, масштаб которых порядка размеров самой области движения l (для трубы порядка её диаметра d). На меньших расстояниях изменение скорости незначительно. Крупномасштабные движения обладают и наибольшей амплитудой пульсаций, т.е. разность скоростей на концах отрезка l велика и может быть сравнимой со средней скоростью движения на расстоянии l , т.е.

* Крупно- и мелкомасштабные пульсации скорости При увеличении числа Рейнольдса появляются всё более и более мелкомасштабные пульсации вплоть до пульсаций самого мелкого масштаба λ0. Области среды размером двигаются при заданном числе Рейнольдса как целое: внутри них частицы движутся с одинаковыми скоростями. На больших расстояниях порядка l изменение скорости определяется, как отмечалось выше, в основном крупномасштабными пульсациями, и по порядку величины оно равно . На малых же расстояниях изменение скорости определяется мелкомасштабными пульсациями. Наряду с числом Рейнольдса , характеризующим турбулентное движение в целом, а так же крупномасштабное движение, можно ввести понятие о числе Рейнольдса Reλ , характеризующим движение в любом масштабе λ от самого мелкого до самого крупного: (10.7.1) В (10.7.1) λ определяет масштаб движения, а υλ - порядок величины скорости пульсационного движения этого масштаба . Движение среды при больших числах Рейнольдса как крупномасштабное движение можно рассматривать как движение маловязкой, в пределе идеальной среды. описываемой уравнениями Эйлера. Отсюда следует, что в крупномасштабных движениях не должно происходить заметной диссипации энергии движения, т.к.вклад вязких сил мал (η ≈ 0). С уменьшением масштаба движения λ уменьшается и скорость υλ, а, тем более, Reλ. Поэтому при уменьшении масштаба движения вязкость играет всё большую роль вплоть до самого мелкого масштаба λ0. Именно в наиболее мелкомасштабных движениях масштаба λ0 происходит основная диссипация энергии турбулентного движения.

* 10.7.2. Диссипируемая энергия при крупномасшщтабном движении Т. о., можно представить следующую схему диссипации энергии механического движения среды при турбулентном движении. Почти вся энергия турбулентного движения содержится в крупномасштабных движениях. Эта энергия передается в движения самых мелких масштабов порядка λ0, в которых и происходит её диссипация в тепло. Т.е., существует направленный поток энергии от крупномасштабных к мелкомасштабным движениям, причем этот поток, проходя промежуточные масштабы, остается почти неизменным. Определим порядок величины диссипируемой энергии ε в единицу времени в единице массы среды. Т.к. крупномасштабное движение можно рассматривать как движение идеальной среды (η = 0), то параметрами, определяющими это движение, являются ρ, Δu, l . Из этих определяющих параметров можно образовать единственную комбинацию, имеющую размерность ε : (10.7.2) Как уже отмечалось, турбулентно движущуюся среду можно рассматривать как среду, обладающую некоторой турбулентной вязкостью ηТ . Очевидно, что турбулентная вязкость также должна определяться крупномасштабным движением, т.е. теми же определяющими параметрами ρ, Δu, l , поскольку характер осредненного движения, очевидно, в основном определяется крупномасштабными пульсациями. Из них можно образовать лишь одну комбинацию, имеющую размерность коэффициента турбулентной вязкости ηТ , порядок величины которого определяется выражением (10.7.3)

* Порядок величины диссипируемой энергии Разделив правую и левую часть (10.7.3) на коэффициент динамической вязкости получим: (10.7.4) В формулах (10.7.3,4), вообще говоря, может стоять значительный безразмерный численный коэффициент. Очевидно, что в формуле (10.7.4) этот численный коэффициент значительно ближе к единице, если поделить правую часть (10.7.4) на Reкр. Тогда при Re ≈ Reкр , что следует и из физических соображений: (10.7.5) Тогда величина диссипируемой энергии может быть определена из (10.7.2,3) следующим образом: (10.7.6) Интересно сравнить полученное выражение с аналогичным выражением для ε при ламинарном движении. Как показано в теме 8, энергия диссипации в единице массы в единицу времени ламинарно движущейся несжимаемой среды определяется по формуле: (10.7.7) Из (10.7.6,7) видно, что формулы, определяющие порядок величины диссипируемой энергии при турбулентном и ламинарном движениях совершенно идентичны при замене η на ηТ .

* 10.7.3. Закон Колмогорова-Обухова Рассмотрим движение в промежуточных масштабах λ0 (l > λ > λ0). Очевидно, в качестве параметров, определяющих движение в промежуточных масштабах можно выбрать ρ, λ и ε - поток энергии, “переходящий” от крупномасштабных движений к движению самых мелких масштабов, в которых он диссипируется в тепло. По-видимому, вязкость для движений в промежуточных масштабах еще не играет заметной роли. Оценим порядок υλ , т.е. величину изменения скорости на концах отрезка порядка λ. Из вышеуказанных параметров можно образовать лишь одну величину с размерностью скорости - . Т.о., имеем: (10.7.8) Из (10.7.8) видно, что изменение скорости на расстояниях порядка λ пропорционально кубическому корню из этого расстояния. Соотношение (10.7.8) называют законом Колмогорова-Обухова для турбулентного движения. Перейдем к рассмотрению самых мелких масштабов движения λ0 , в которых и происходит диссипация энергии турбулентного движения вследствие вязкости среды. Определим порядок величины самых мелких масштабов движения λ0 . Рассмотрим число Рейнольдса масштаба λ : (10.7.9)

* Величина мелкого масштаба Очевидно, для самых мелких масштабов движения, когда вязкость играет определяющую роль, . Тогда из (10.7.9) следует: (10.7.10) Формулу (10.7.10) можно получить, рассматривая в качестве определяющих параметров движения в масштабах λ0 параметры ρ , ε и η. Формулой (10.7.10) определяется размер областей среды, которые при заданном числе Рейнольдса двигаются как целое. Чем больше число Рейнольдса, тем меньше размер этих областей. Эксперименты подтверждают увеличение мелкомасштабных вихрей в турбулентном движении среды с увеличением числа Рейнольдса в области Re > Reкр . .

* Выводы Введены основные определения и получены основные соотношения, описывающие турбулентное движение вязкой среды: логарифмический профиль скорости, вязкий подслой вблизи обтекаемой поверхности, зависимость коэффициента гидравлического сопротивления от числа Рейнольдса, уравнения Никурадзе для шероховатых труб, крупно- и мелкомасштабные пульсации скорости, закон Колмогорова-Обухова для скорости турбулентного движения. 2008. Численные методы…Лекция 19

* Информационное обеспечение лекции 1.Бабаков А.В., Белоцерковский О.М., Северинов Л.И. Численное исследование течения вязкого теплопроводного газа у тупого тела конечных размеров.- Изв. АН СССР. Сер. Механ. жидк. и газа, 1975, №3, с. 112-123. 2. Бабаков А.В. Численное моделирование некоторых задач аэрогидродинамики.-М.: ВЦ АН СССР, 1986, с. 56. 3. Белоцерковский О.М., Северинов Л.И. Консервативный метод потоков и расчет обтекания тела конечных размеров вязким теплопроводным газом.- ЖВМ и МФ, 1973, 13, №2, с. 385-397. 4. Самарский А.А. О консервативных разностных схемах.- В кн.: Проблемы прикл. матем. и механ. М.: Наука, 1971, с.129-136. 2008. Численные методы…Лекция 19

* Справочные данные Курс лекций является частью учебно-методического комплекса «Численные методы расчета задач механики сплошных сред. 1. Теория упругости и идеальная среда». Автор: Породнов Борис Трифонович, д. ф. – м. н., профессор кафедры молекулярной физики УГТУ-УПИ. Учебно-методический комплекс подготовлен на кафедре МФ ФТФ ГОУ ВПО УГТУ-УПИ. электронный адрес: [email protected]