Авторефераты по всем темам  >>  Авторефераты по разным специальностям

На правах рукописи

Хребтов Михаил Юрьевич

СПЕКТРАЛЬНЫЙ ПЕРЕНОС ЭНЕРГИИ ТУРБУЛЕНТНОСТИ В КРУГЛОЙ ЗАТОПЛЕННОЙ СТРУЕ

01.02.05 - механика жидкости, газа и плазмы

Автореферат диссертации на соискание учёной степени кандидата физико-математических наук

Новосибирск - 2012

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте теплофизики им. С.С. Кутателадзе Сибирского отделения РАН

Научный консультант: доктор физико-математических наук Илюшин Борис Борисович

Официальные оппоненты:

Чекмарев Сергей Федорович - доктор физико-математических наук, профессор, ФГБУН Институт теплофизики им. С.С. Кутателадзе Сибирского отделения РАН, главный научный сотрудник.

Яковенко Сергей Николаевич - кандидат физико-математических наук, ФГБУН Институт теоретической и прикладной механики им С.А.

Христиановича СО РАН, старший научный сотрудник.

Ведущая организация:

ФГБУН Институт Вычислительных Технологий СО РАН

Защита состоится 03 октября 2012 г. в 11 часов 30 минут на заседании совета Д 003.053.01 по защите диссертаций на соискание ученой степени кандидата наук, на соискание ученой степени доктора наук, в ФГБУН Институте теплофизики им. С.С. Кутателадзе СО РАН по адресу: 630090, г.Новосибирск, пр. Лаврентьева, 1.

С диссертацией можно ознакомиться в библиотеке Института теплофизики СО РАН.

Автореферат разослан 31 августа 2012 г.

Учёный секретарь диссертационного совета д.ф.-м.н. Кузнецов Владимир Васильевич

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

В области численного моделирования наиболее теоретически обоснованным является направление прямого численного решения уравнений Навье-Стокса (DNS - Direct Numerical Simulation), однако текущий уровень производительности суперкомпьютерных систем позволяет проводить расчеты по такой методике только для относительно небольших чисел Рейнольдса (до нескольких тысяч). Другим методом, активно используемым в инженерных задачах, является метод решения осредненных уравнений Навье-Стокса с использование разных видов замыканий (RANS - Reynolds-Averaged Navier-Stokes). В случаях, когда используемые для замыкания гипотезы не выполняются, этот метод может давать существенное расхождение с экспериментом. Промежуточное положение между этими методами занимает метод крупных вихрей (LES - Large Eddy Simulation), в котором мгновенные поля скорости представляются как сумма пульсаций мелких УподсеточныхФ масштабов и отфильтрованной на масштабе сетки скорости, соответствующей крупномасштабным пульсациям. Уравнения для отфильтрованной скорости разрешаются напрямую, а влияние подсеточных масштабов учитывается с помощью введения параметризации на основе какой-либо алгебраической модели. Такой метод позволяет моделировать течения с большими числами Рейнольдса чем DNS, при этом, давая возможность наблюдать за эволюцией основных энергосодержащих масштабов потока во времени. Таким образом, метод LES является на данный момент разумным компромиссом для исследования спектральной передачи энергии турбулентности.

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

Основные задачи исследования На основе LES-моделирования исследовались: направление и амплитуда спектрального переноса энергии турбулентности в невозмущенной струе, механизм формирования зон обратного потока энергии. Также исследовалось влияние на спектральный перенос в свободной струе параметров течения (число Рейнольдса, закрутка, входной профиль, низкоамплитудные возмущения, спутный поток).

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

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

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

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

Проведены исследования изменения временных спектров в слое смешения при добавлении внешних низкоамплитудных возмущений. В спектрах на границе слоя смешения присутствует интервал с постоянным наклоном близким к k-3, наблюдающимся в двумерной турбулентности.

Проведено исследование влияния числа Рейнольдса потока на спектральный перенос в струе. Показано, что с увеличением числа Рейнольдса, амплитуда спектрального переноса (как прямого, так и обратного участков) увеличивается. При этом максимум амплитуды обратного потока смещается с увеличением числа Рейнольдса вверх по потоку. Кроме того, уменьшение числа Рейнольдса меньше 25000 приводит к появлению дополнительного участка с обратным переносом на внешней границе слоя смешения, отсутствующего для больших Re.

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

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

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

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

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

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

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

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

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

Результаты анализа влияния закрутки и спутного потока на спектральный перенос в струе.

ичный вклад соискателя Соискателем были проведены: модификация расчетного модуля, добавление и тестирование новой модели подсеточной вязкости, добавление модуля для расчета спектрального потока энергии и блока для визуализации вихревых структур, модификация расчетной сетки, а также изменение начальных и граничных условий расчета. Соискателем было проведено LESмоделирование свободных струй с закруткой, внешним возбуждением, спутным потоком. Из расчетов была получена подробная база статистических и мгновенных характеристик течений, в том числе распределение спектрального потока энергии для диапазона пространственных масштабов. Были проведены обработка и анализ результатов.

Апробация работы Основные результаты диссертации докладывались и обсуждались на всероссийских и международных конференциях. Опубликованы в трех журнальных статьях (2 в журналах из перечня ВАК), и 13 тезисах докладов.

Объем работы Диссертация состоит из введения, 5 глав, заключения, списка литературы.

Общий объем диссертации - 145 страниц текста, 89 рисунков и 1 таблица.

Список литературы насчитывает 99 наименований.

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

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

Приведены ссылки на современные работы в этом направлении.

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

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

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

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

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

Приведено выражение для вычисления спектрального потока энергии через волновое число k - K. Выражение получается путем применения операции низкочастотного фильтра к уравнению Навье-Стокса, и последующего его домножения на пульсационную скорость, отфильтрованную на том же масштабе. Для метода LES необходимо модифицировать это выражение, чтобы учесть поток от подсеточных масштабов:

K K K K K LES sgs vi vj vi vi Ti K K j ui uj Ti sgs xj xj xi Здесь vi Цпульсационная компонента скорости ui, угловые скобки означают осреднение по времени, а черта сверху - низкочастотную фильтрацию с граничным волновым числом K.K имеет смысл количества энергии переносимого в единицу времени от пульсаций с волновыми числами меньшими K, к пульсациям с волновыми числами большими K. Если K имеет отрицательный знак, это означает, что процесс идет в обратном направлении. Поток энергии был обезразмерен через входную скорость и диаметр сопла. Структурная функция третьего порядка Brrr обезразмеривалась через входной диаметр, входную скорость, и шаг сетки.

Din Din Brrr (r) * K ; B* K Uin Uin r Во второй главе приведено подробное описание объекта исследования и сравнение получаемых в расчете данных с экспериментом. Объектом исследования являлась затопленная осесиммтеричная струя. Сопло (Din=10мм) было ограничено шайбой (54мм). Приведены параметры расчетной области и сетки (цилиндр, габаритами 21x20Din, сетка до 6x1узлов). Числа Рейнольдса варьировались от 12500 до 200000. Закрутка S варьировалась от нуля до 0.84. Скорость спутного потока варьировалась от нуля до 0.7 от входной скорости струи. Частота внешнего возбуждения варьировалась от Sh=0.25 до Sh=1.0 (при амплитуде в диапазоне 0.25%-1% от входной скорости).

Рис. 1. Геометрия течения.

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

Рис.2. Радиальный профиль третьего момента продольной скорости (разные модели подсеточной вязкости и шаг сетки). x/Din=16.5, Re=25000.

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

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

Рис.3. Структурная функция третьего порядка в свободной струе из PIV эксперимента (показаны только положительные значения), слева, и полученная в расчетах (модель Cмагоринского) справа.

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

При смещении вниз по потоку от сопла видно увеличение уровня турбулентных пульсаций одновременно с увеличением амплитуды B*.

Приведена полученная в расчетах вихревая структура начального участка струи, из которой видно, что максимум амплитуды структурной функции расположен там (x/D=2.5), где в результате развития неустойчивости, уже появилось достаточно большое число разномасштабных вихрей.

Приведены временные спектры энергии турбулентных пульсаций, полученные в расчете, в различных точках струи. Видно, что в сечении, где наблюдается максимум B* (x/D=2.5) спектры являются достаточно наполненными (что согласуется с экспериментальными исследованиями других авторов), кроме того, на периферии слоя смешения, в спектре наблюдается участок с наклоном E~k-3, там же наблюдаются и положительные значения B*. В центре же слоя смешения наклон близок к колмогоровскому и B* имеет там отрицательные значения. С удалением от сопла, инерционный интервал с наклоном k-5/3 в спектрах расширяется, наклон -3 перестает наблюдаться.

Рис.4. Вихревые структуры на начальном участке струи, визуализация по 2критерию. Эллипсом обозначена область максимума B*. Re=25000.

Анализ распределения коэффициента эксцесса скорости показывает, что максимум B* расположен в радиальном сечении между максимумами коэффициента эксцесса, указывающими на зоны интенсивной перемежаемости (и выделяющими границу струи). В области максимума B* коэффициент эксцесса продольной скорости имеет незначительную амплитуду, что говорит о приемлемой статистической однородности распределения пульсаций в этой области.

Вычисление распределения потока энергии по спектру K в сечении максимума B* показывает, что в частотном распределении присутствует интервал, где поток энергии имеет отрицательный знак. Причем этот интервал расположен в области относительно высоких частот (минимум соответствует Sh=4). Пик из-за прохождения когерентных структур, наблюдающийся в спектрах энергии, соответствует участку с положительным знаком потока энергии.

Рис..5. Спектры энергии турбулентности на разных расстояниях от сопла, полученные в расчете. (использовалась модель Времана, Re=25000).

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

При этом следует отметить, что для рассматриваемого сечения (x/D=2.5) поток нигде не остается постоянным по спектру, отсутствуют протяженные горизонтальные участки. Однако ниже по течению, например в сечении x/D=8, поток в частотном распределении практически постоянен, образуя протяженный горизонтальный участок, что согласуется с идеей об инерционном интервале, и видом спектра в этой области.

Рис. 6. Поток энергии по спектру в зависимости от волнового числа фильтрации для различных точек в струе (Re=25000).

Приведено пространственное распределение потока энергии через три выделенных масштаба течения. Распределение показывает концентрацию обратного спектрального переноса на начальном участке струи, на внутренней границе слоя смешения. При этом область обратного переноса простирается даже дальше смыкания слоя смешения вплоть до сечения x/D=7.

Рис. 7. Поток энергии через масштаб L=0.25D в продольном сечении струи (Re=25000).

Раздел 3.2. представляет собой попытку дать объяснение возможных механизмов появления обратного переноса на начальном участке струи.

Анализ вихревой структуры струи в зоне максимального обратного переноса на внутренней границе слоя смешения показывает, что там происходит интенсивная деформация вихрей в плоскости, перпендикулярной направлению их вращения. Такая деформация характерна для двумерных течений. За счет расширения слоя смешения средний сдвиг направляет вихри друг к другу, что приводит к интенсификации их взаимодействия между собой. Так как вихри, в связи с близостью к границе сопла, все еще имеют выделенное направление вращения, из-за осевой симметрии геометрии течения, то их взаимодействие будет похожим на взаимодействие вихрей в двумерной турбулентности. При этом следует отметить, что при сопоставлении вихревой структуры струи с распределением спектрального переноса было выявлено, что ядра когерентных структур практически не попадают в область обратного потока. В основном, в область обратного переноса попадают образующиеся в результате вторичных неустойчивостей вихри, с поперечным масштабом в несколько раз меньше чем основные когерентные структуры (дающие пик в спектре на частотах Sh~1), что согласуется с приведенным выше анализом графика частотного распределения потока энергии (с максимумом обратного потока на частоте Sh~4).

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

Раздел 3.3 посвящен исследованию влияния числа Рейнольдса на спектральный перенос энергии в струе. С этой целью были выполнены расчеты струи для чисел Рейнольдса в диапазоне 12500-200000. На приведенных графиках пространственного распределения спектрального переноса видно, что увеличение числа Рейнольдса приводит к уменьшению области обратного переноса на внешней границе слоя смешения, но в то же время, интенсивность обратного переноса на внутренней границе слоя увеличивается. Кроме того, с увеличением Re происходит расширение области обратного переноса на внутренней границе слоя с одновременным смещением максимума амплитуды ближе к соплу.

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

В главе 4 исследуется влияние внешнего гармонического возмущения на спектральный перенос энергии турбулентности в струе.

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

Приводятся экспериментальные данные исследования свободных струй. Смещение пика в спектрах, связанного с прохождением когерентных структур заканчивается в момент смыкания слоя смешения (x/D=4-5), и фиксируется на частоте около Sh=0.3. Эта частота затем проявляется в спектрах вплоть до x/D=8-10 и называется собственной частотой струи. При этом, когерентные структуры в виде клубков вихрей на внешней границе слоя смешения проходят с данной частотой для широкого диапазона чисел Рейнольдса.

Рис. 8. Распределение спектрального потока энергии через масштаб 0.25D в радиальном сечении невозмущенной струи для различных чисел Рейнольдса (x/D=2.5) Доминирование одной частоты на достаточно протяженном участке струи является логичным основанием для изучения возмущения потока именно на данной частоте. Как следует из литературы, возмущения на собственной частоте струи оказывают наиболее сильное влияние на амплитуду турбулентных пульсаций. Однако воздействие выходит на насыщение для амплитуд возмущения ~4% от входной скорости и более.

Несмотря на такую небольшую амплитуду входных возмущений, изменение итоговой амплитуды пульсаций из-за такого воздействия может быть в пять раз больше Для других частот максимум влияния возмущений на поток меньше по амплитуде. Однако если исследовать фиксированные небольшие амплитуды возмущений, то для них воздействие больших частот оказывается более сильным, чем воздействие собственной частоты. В данной работе исследовался диапазон (Sh=0.25-1.0) для амплитуд (0.25%-1%).

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

Анализ вихревой структуры показывает, что внешнее возмущение оказывает упорядочивающее влияние на струю, приводя к доминированию одной моды пульсаций и подавлению всех остальных. При этом можно видеть, что для частоты Sh=0.5, несмотря на подавление фоновых пульсаций, вихри не имеют строго осесимметричной формы, напоминая спирали, что можно объяснить тем что эта частота попадает между двумя дискретными кольцевыми модами струи. На частоте Sh=0.75 вихри имеют кольцевую форму на всем протяжении начального участка. На частоте Sh=1.0 можно видеть, что кольцевые вихри начинают объединяться второй раз попарно (чего не наблюдается для других частот), что иллюстрирует сценарий появления половинной субгармоники.

Анализ распределения спектрального потока энергии показывает, что для частоты Sh=0.75 влияние возмущения максимально. При этом области обратного потока расширяются в струе, кроме того, максимум обратного потока для крупных масштабов (0.6Din и более) теперь выходит на внешнюю границу слоя смешения чего не наблюдается для других частот. Увеличение амплитуды максимума составляет до двух раз, по сравнению с невозмущенным случаем.

Рис. 9. Вихревые структуры на начальном участке струи для различных частот возмущения (амплитуда 0.5%, Re=25000), изоповерхности 2=-60.

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

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

Отметим, что для точки x/D=5, y/D=1.0 спектр перед изломом имеет наклон близкий к -5/3, а после излома близкий к -3. Что напоминает крейчнановский спектр двумерной турбулентности.

При увеличении амплитуды пульсаций (0.5-1%) влияние входных возмущений проявляется в виде отчетливых пиков в спектрах, при этом в пиках сосредотачивается существенная часть энергии потока, так что остальной спектр становится ниже. Как и ожидалось, пики проявляются тем отчетливее, чем ближе мы находимся к оси сопла. Различия для разных частот состоят в различных амплитудах субгармоники основной частоты.

Так, например, для частот Sh=0.75, Sh=1.0 половинная субгармоника довольно быстро начинает превосходить по амплитуде основную гармонику, в то время как для частоты Sh=0.5 такого сильного роста субгармоники не происходит, ее амплитуда остается примерно в два раза меньше основной.

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

Рис. 10. Распределение спектрального потока энергии через масштаб 0.25Din в радиальном сечении для случая невозмущенной струи (сплошная линия) и с наложением внешнего возмущения Sh=0.75, amp 0.5% невозмущенной струи для различных чисел Рейнольдса (x/D=2.5) Рис. 11. Спектры энергии турбулентности для различных частот и амплитуд возмущения.

В главе 5 в первой части описывается исследование вопроса о спектральном переносе энергии в закрученных струях с входным профилем продольной скорости, взятым из невозмущенной струи и близким к гиперболическому тангенсу. Такой профиль был выбран, чтобы выделить влияние на спектральный перенос самой закрутки, исключая влияние изменения входного профиля. Проведены расчеты с различными интенсивностями закрутки. Профиль окружной компоненты на первом этапе был близок к твердотельному вращению. Приведены распределения компонент средней скорости и пульсаций для закрутки, изменяющейся в пределах (S=00.84). Видно, что с увеличением закрутки происходит усиление присоединенного вихря на начальном участке течения. Кроме того, наблюдается тенденция к смещению всех стадий развития потока ближе к соплу и более быстрое расширение струи, по сравнению с незакрученным случаем. Амплитуда пульсаций также увеличивается с увеличением закрутки. Однако для относительно небольшой закрутки (вплоть до S=0.56) схема потока остается качественно такой же, как и в незакрученном случае.

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

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

Рис. 12. Характерный вид вихревых структур на начальном участке струи для режимов течения с закруткой S=0.56.

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

Рис. 13. Распределение спектрального потока энергии через масштаб 0.25D в радиальном сечении для различных величин закрутки (Re=25000).

Причем, чем больше закрутка, тем сильнее эта интенсификация. Если сравнивать амплитуду абсолютного минимума мелкомасштабного спектрального переноса, то для закрутки S=0.56 амплитуда увеличивается примерно в два раза по сравнению с незакрученной струей, хотя сам минимуму сдвигается вверх по течению примерно на калибр. Как было сказано выше, эта интенсификация, по-видимому, вызвана взаимодействием двух частей вихревой структуры струи между собой.

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

Это приводит к тому, что вместо одного максимума в распределении энергии турбулентности появляется два - один в центре слоя смешения, другой - вблизи оси. Изменение профиля одной из компонент существенно влияет на другие компоненты средней скорости, в частности, в несколько раз уменьшается амплитуда продольной компоненты вблизи оси сопла.

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

Рис. 14. Влияние вида профиля окружной компоненты скорости на распределение пульсаций, и поток энергии в струе. S=0.5, Re=25000.

Далее в главе исследовалось влияние профиля продольной скорости и спутного потока на спектральный перенос энергии.

Профили входной и окружной скорости были взяты из PIV эксперимента для закрученной струи. Число Рейнольдса было снижено в два раза (Re=12500). Кроме того, за границей шайбы (y/D=2.71) задавался спутный поток различной амплитуды (10%-60% от входной скорости).

Для случая закрутки S=0.43, несмотря на существенные отличия в распределении средних скоростей, картина спектрального переноса в целом похожа на то, что наблюдается для невозмущенного случая за исключением того, что области прямого и обратного переноса начинаются прямо от оси сопла. Если бы энергия для обратного переноса отбиралась вихрями от среднего течения, то изменение градиента средней скорости должно было сильно повлиять на спектральный перенос, однако расчет показывает, что амплитуда обратного переноса остается практически такой же при том, что градиент средней скорости поменялся в несколько раз. Это говорит о том, что обратный перенос является следствием динамического взаимодействия вихрей между собой.

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

При изменении амплитуды спутного потока с 10% до 20% (от входной скорости струи) этот сдвиг составляет примерно 1 калибр. При этом, амплитуда обратного переноса в этой области также увеличивается.

Далее представлены результаты расчетов для случая большей закрутки (S=0.75) и больших амплитуд спутного потока (36%, 60%). Распределение спектрального переноса показывает, что теперь все области как прямого, так и обратного переноса струи сократились в несколько раз. Взаимодействие со спутным потоком приводит к тому, что струя притягивает спутный поток, так что он уменьшает свой радиус, и как следствие, завихренность переносимая спутным потоком должна испытывать эффект концентрации, это должно приводить к появлению областей обратного переноса в спутном потоке.

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

Рис. 15. Спектральный перенос энергии через масштаб 0.25D для различных величин закрутки и спутного потока.

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

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

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

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

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

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

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

Сформулированы выводы:

1. Выполнено LES-моделирование круглой затопленной струи с использованием моделей подсеточной турбулентности Смагоринского и Времана.

2. Результаты вычислений сопоставлены с данными эксперимента и дали хорошее соответствие по средним, вторым и третьим моментам скорости.

3. Вычислены продольные структурные функции потока третьего порядка, а также спектральный поток энергии турбулентности.

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

4. С целью исследования механизмов формирования зон обратного переноса энергии турбулентности выполнено LES-моделирование струи в широком диапазоне чисел Рейнольдса, а также с внешними наложенными низкоамплитудными гармоническими возмущениями и закруткой.

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

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

Основные результаты диссертации опубликованы в следующих работах:

1. M.Yu. Hrebtov, B.B. Ilyushin, D.V. Krasinsky Inverse Energy Cascade in Turbulent Round Jet //Phys.Rev.E 81, 016315 (2010) (из перечня ВАК) 2. Б.Б. Илюшин, М.Ю. Хребтов LES-моделирование обратного каскада энергии в 3D турбулентности. //Доклады РАН (ДАН). 2011. Т.437, № 4, с.485Ц487 (из перечня ВАК) 3. M.Yu. Hrebtov, B.B. Ilyushin Formation of Inverse Energy Cascade in Free Turbulent Jet // VI-th Int. conf. Solitons Collapses and Turbulence, Novosibirsk, Russia, 4-8 June 2012, pp.79-4. B.B. Ilyushin, M.Hrebtov Inverse energy transfer mechanism in turbulent round jets // Turbulence, Heat and Mass Transfer 7: Proc. of the 6-th Int. Symp., Palermo, Italy, 23-27 September 2012, pp.431-45. B.B. Ilyushin, D.V. Krasinsky and M.Yu. Hrebtov LES study of turbulent energy flux in a free circular jet // Turbulence, Heat and Mass Transfer 6: Proc. of the 6-th Int. Symp., Rome, Italy, 14-18 September 2009, pp.355-36. М.Ю.Хребтов. Спектральный перенос энергии в закрученных струях.

//IX Всероссийская конференция молодых ученых Проблемы механики:

теория, эксперимент и новые технологии Новосибирск, 2012, с. 35Ц37.

7.М.Ю.Хребтов Механизмы спектрального переноса энергии турбулентности в свободной струе. //тезисы XI Всероссийской школыконференции молодых учёных Актуальные вопросы теплофизики и физической гидрогазодинамики Новосибирск, 2010, c.22-26.

8.М.Ю.Хребтов Сравнение критериев идентификации вихрей на основе моделирования турбулентной струи методом LES // Тез. Докл. IX Всеросс.

школы-конф. молодых ученых Актуальные вопросы теплофизики и физической гидрогазодинамики, Новосибирск, 2006, с.129-130.

9.Б.Б. Илюшин, Д.В. Красинский, М.Ю. Хребтов Исследование структуры турбулентной струи методом крупных вихрей (LES) // Тез. Докл.

Всеросс. семинара Современные проблемы теоретической и прикладной механики (под ред. проф. В.Я. Рудяка), Новосибирск,2007, с.45Ц46.

10.М.Ю.Хребтов Определение вихревой структуры круглой турбулентной струи на основе метода LES // Материалы XLV Междунар.

Науч.Студ.Конф. Студент и научно-технический прогресс, Новосибирск, 2007, с. 44-45.

11. B.B. Ilyushin, D.V. Krasinsky, and M.Yu. Hrebtov LES Study of the Vortex Structure of Turbulent Round Jet // Proc. of 13-th Int. Conf. on the Methods of Aerophys. Research (ICMAR-2007), Novosibirsk, 2007, pp.141Ц146.

12.М.Ю.Хребтов Исследование эволюции крупномасштабных вихревых структур в свободной турбулентной струе // Тез. докл. Всеросс.

школы-семинара молодых ученых Физика неравновесных процессов в энергетике и наноиндустрии, Новосибирск, 2007, с. 131-132.

13. М.Ю. Хребтов Исследование статистических характеристик свободной турбулентной струи // Тезисы ХLVI Междунар. научной студенческой конференции Студент и научно-технический прогресс, Новосибирск, 2008, с.112-113.

14. М.Ю. Хребтов Влияние внешнего возбуждения на статистические характеристики турбулентной струи // Тезисы ХLVII Междунар. научной студенческой конференции Студент и научно-технический прогресс, Новосибирск, 2009, с.55.

15. Б.Б.Илюшин, М.Ю.Хребтов Формирование зон обратного каскада энергии в осесимметричной струе. // Журнал "Современная наука:

исследования, идеи, результаты, технологии" V.7, 2011, c.183-187.

16. М.Ю.Хребтов "Численное решение уравнения Навье-Стокса в калибровочных переменных". XII Всеросс. конф. мол. уч. по мат. мод. и информационным технологиям, Новосибирск, 2011, c.55-58.

Авторефераты по всем темам  >>  Авторефераты по разным специальностям