ИСТИНА |
Войти в систему Регистрация |
|
ФНКЦ РР |
||
Проект направлен на экспериментальное, численное и теоретическое исследование развития возмущений и перехода к турбулентности в ламинарных затопленных струях, поиск новых методов управления струйными течениями. В настоящее время достигнут значительный прогресс в понимании механизмов развития возмущений для ряда важных течений, таких как пограничные слои, течения в трубах и каналах, течения между вращающимися цилиндрами и некоторых других. Вместе с тем, представления о характере развития возмущений на этапе ламинарно-турбулентного перехода в струйных течениях находятся на гораздо более низком уровне. Это связано с тем, что в большинстве струй переход наступает уже на расстоянии в единицы калибров от среза сопла, что чрезвычайно затрудняет экспериментальное изучение обстоятельств, сопровождающих переход. Участниками проекта ранее был разработан новый метод формирования затопленных струй диаметром 120 мм с длинным ламинарным участком (600 – 700 мм) при больших числах Рейнольдса (2000 – 13000). Такие струи будут объектом исследования в данном Проекте. Помимо экспериментальных исследований, теоретически будут исследованы линейные механизмы роста возмущений (как модальный, так и алгебраический). Параллельно будет проводиться прямое численное моделирование (DNS) изучаемых течений, что по нашему опыту существенно повышает достоверность и надежность получаемых результатов. По итогам выполнения проекта будут созданы научные основы для разработки способов управления ламинарными струйными течениями, как с целью затягивания перехода к турбулентности, так и для его ускорения. Эти способы могут быть использованы в ряде отраслей промышленности для повышения эффективности ряда технологических процессов.
The project is aimed at experimental, numerical and theoretical studies of the development of disturbances and the transition to turbulence in laminar submerged jets, the search for new methods of controlling jet flows. Currently, significant progress has been made in understanding the mechanisms of perturbation development for a number of important flows, such as boundary layers, flows in pipes and channels, flows between rotating cylinders, and some others. At the same time, ideas about the nature of the development of perturbations at the stage of the laminar-turbulent transition in jet flows are at a much lower level. This is due to the fact that in most jets the transition occurs already at a distance of units of calibers from the nozzle exit, which greatly complicates the experimental study of the circumstances that accompany the transition. Project participants previously developed a new method for forming submerged jets with a diameter of 120 mm with a long laminar section (600 - 700 mm) with large Reynolds numbers (2000 - 13000). Such jets will be the object of study in this Project. In addition to experimental studies, the linear mechanisms of perturbation growth (both modal and algebraic) will be theoretically investigated. In parallel, direct numerical simulation (DNS) of the studied flows will be carried out, which in our experience significantly increases the reliability and reliability of the results. Based on the results of the project, scientific foundations will be created to develop methods for controlling laminar jet flows, both with the aim of delaying the transition to turbulence, and to accelerate it. These methods can be used in a number of industries to increase the efficiency of a number of technological processes.
Впервые будут экспериментально получены скорости роста контролируемых возмущений в ламинарной струе и проведено сравнение результатов с теоретическими расчётами модальной теории устойчивости. С помощью прямого численного моделирования будет изучено развитие таких возмущений на нелинейной стадии. Будет исследована применимость e^n-метода оценки расстояния до зоны перехода, широко используемого в пограничных слоях, но применимость которого к струйным течениям неизвестна. Будут экспериментально и численно изучены особенности перехода к турбулентности посредством немодального роста возмущений. Будет теоретически и экспериментально исследована абсолютная неустойчивость струи без возвратного течения вокруг нее, чего ранее не проводилось, ввиду того, что струи круглого поперечного сечения с автомодельным профилем скорости (и близкими к нему) неустойчивы конвективно. Численно и экспериментально будет установлено, к какому типу течения приведет развитие абсолютной неустойчивости – ко вторичному течению или к турбулентности. Понимание механизмов перехода к турбулентности имеет не только общенаучное, но и важное прикладное значение в контексте поиска методов управления длиной ламинарного участка струи. Струи с длинными ламинарными участками могут использоваться для создания локальных чистых зон: струя очищенного воздуха, не перемешивающаяся с окружающим воздухом, помещается над операционным полем (местом сборки или обработки в электронике, оптике, фармацевтических, химических и микробиологических лабораториях) и обеспечивает больший уровень чистоты, нежели в окружающем помещении. Подобные не отгороженные от оператора твердыми стенками чистые зоны в ряде случаев сильно облегчают манипуляции с защищаемым объектом по сравнению с ламинарными шкафами и боксами, при этом обеспечивают защиту объекта от загрязнений. Турбулизация же в непосредственной близости от начального сечения струи имеет приложение в технологиях, где необходимо эффективное перемешивание (например, в форсунках камер сгорания авиационных двигателей). Таким образом, результаты Проекта составят научную основу для создания новых инновационных устройств и технологий.
1. Разработана методика затягивания перехода к турбулентности в круглых затопленных струях. Проведены детальные измерения профилей скорости и турбулентных пульсаций получаемой струи в диапазоне чисел Рейнольдса 2000…13000 и визуализация получаемых струй. Создано формирующее струю диаметром 120 мм устройство. Формируемая струя имеет ламинарный участок 660 мм при числе Рейнольдса ~ 10000 (Зайко и др. Изв. РАН. МЖГ. 2018). 2. Проведены первые теоретические расчёты в рамках линейного анализа устойчивости струи с профилями, полученными в экспериментах. Показано, что частоты наиболее быстро растущих возмущений, найденные при линейном анализе, совпадают с частотами из экспериментального спектра для соответствующего скоростного режима (Zayko, et al. Physics of Fluids. 2018). 3. Экспериментально показано, что при частотах колебаний полосы фольги, расположенной под струей, из диапазона наиболее быстро растущих возмущений ламинарный участок струи значительно сокращается, тогда как на меньших и больших частотах колебаний фольги влияние на струю незначительное (http://www.fiv2018.com/ocs/index.php/fiv2018/fivcan/paper/view/109/37). 4. Проведено предварительное исследование характера неустойчивости и показано, что в струйных течениях с профилями скорости специального вида возможна абсолютная неустойчивость течения (V. Vedeneev, J. Zayko. On absolute instability of free jets// Journal of Physics: Conference Series. 2018. Vol. 1129. 012037.). 5. У коллектива имеется весь необходимый методический багаж для проведения численных исследований заявленной программы. Численное интегрирование уравнений Навье-Стокса будет проводиться с использованием методик расчета турбулентных течений, разработанных участниками проекта. Большинство имеющихся программ для ЭВМ написаны в двух версиях - для применения на единичных процессорах и на многопроцессорных системах. У коллектива имеется также опыт разработки алгоритмов, ориентированных на применение в многопроцессорном режиме.
грант РНФ |
# | Сроки | Название |
1 | 1 апреля 2020 г.-31 декабря 2020 г. | Исследование развития возмущений в ламинарных струйных течениях и разработка научных основ новых способов управления такими течениями. |
Результаты этапа: Проведено экспериментальное изучение развития волн в затопленной струе круглого поперечного сечения с длинным ламинарным участком. Указанная струя была получена с помощью формирующего устройства, созданного в лаборатории экспериментальной гидродинамики НИИ механики МГУ. Результаты экспериментов сравнены с предсказаниями линейной невязкой теории устойчивости, применённой в квазипараллельном приближении к профилям скорости струи, полученным в экспериментах и аппроксимированным аналитическими функциями. Показано, что линейная невязкая теория устойчивости способна предсказать свойства волн в ламинарной затопленной струе круглого поперечного сечения при числе Рейнольдса ~5000. В экспериментах усиливались волны, соответствующие полученным теоретически собственным модам обеих ветвей растущих возмущений. Длины указанных волн, установленные экспериментально, близки к длинам волн, предсказанным модальной теорией устойчивости. Таким образом, обосновано применение невязкой модальной теории устойчивости для прогнозирования поведения собственных мод затопленной струи. Написана и отлажена программа для решения линейной задачи устойчивости заданного профиля скорости (в параллельном приближении, задача Орра-Зоммерфельда). Исследована устойчивость начального профиля скорости для Re=5700 по отношению к осесимметричным возмущениям (n=0). Построены зависимости скорости нарастания возмущений от частоты. Как и в невязком случае, обнаружено существование двух ветвей. Первая ветвь характеризуется относительно небольшой скоростью роста возмущений, которая имеет непрерывную зависимость от частоты. На второй ветви характерные скорости роста значительно большие, а зависимость от частоты имеет разрыв. Результаты расчетов задачи Орра-Зоммерфельда для n=0 хорошо согласуется с результатами, полученными ранее в рамках уравнения Рэлея. Проведено исследование развития малых возмущений в струях с различными профилями на входе и при разных значениях числа Рейнольдса. Типичное поведение растущих возмущений выглядит как рост в начальных сечениях с последующим замедлением скорости роста и затуханием на больших расстояниях. Исследовано поведение возмущений в течениях с входным профилем Пуазейля и со ступенчатым профилем при Re=2000. Последнее оказывается существенно более неустойчивым, чем первое. Наибольший рост возмущений происходит при n=1. Также исследованы течения с экспериментальными входными профилями при Re=2000, 5700 и 9200. | ||
2 | 1 января 2021 г.-31 декабря 2021 г. | Исследование развития возмущений в ламинарных струйных течениях и разработка научных основ новых способов управления такими течениями. |
Результаты этапа: Используя алгоритм построения оптимальных возмущений, был теоретически исследован немодальный рост возмущений для профиля скорости струи, соответствующего экспериментальному. Для значений азимутального волнового числа m=0,1,…,5 построен график зависимости энергии оптимального возмущения от безразмерного расстояния вниз по потоку z/D. Визуализированы оптимальные возмущения для различных значений z. Подтверждены основные свойства lift-up механизма немодального роста, ранее изученного для пристенных течений, в случае круглой струи. Обнаружено, что наибольшая энергия роста возмущения для заданного интервала вниз по потоку отвечает азимутальному волновому числу 3, для волновых чисел m ≥5 энергия оптимальных возмущений значительно ниже, чем для промежутка m∈[1,4]. Модифицирована система внесения возмущений в струю, изготовлена серия дефлекторов, соответствующих возмущениям на различных модах неустойчивости течения. Проведены серии экспериментов: с визуализацией течения лазерным ножом; с определением величин продольных скоростей термоанемометром; с измерением распределения поперечной скорости методом PIV. Полученные экспериментальные данные позволяют положительно судить о немодальном росте вносимых возмущений и присутствию lift-up механизма роста. Обнаружена структура вихревого поперечного движения. Показано, что каждой волне дефлектора соответствует пара («основных») вихрей с противоположными завихренностями, причем, как правило, к каждому из них примыкает вихрь с меньшей интенсивностью того же знака, расположенный ближе к оси струи. Проведено численное исследование устойчивости затопленной несжимаемой струи, вытекающей из круглого сопла в полупространство. Найденное двумерное стационарное поле скорости (основное течение) исследуется на устойчивость по отношению к осесимметричным (n=0) и трёхмерным (n≠0) возмущениям. Исследование устойчивости проведено в двух постановках: квазипараллельной и пространственной. В первом случае изменение струи вниз по потоку учитывается параметрически: профиль скорости в каждом из сечений считается независящим от продольной координаты x. Во второй, пространственной постановке, линеаризованные относительно возмущений уравнения Навье-Стокса решаются во всей двумерной области. Таким образом определяется эволюция возмущений вниз по потоку при каждой заданной частоте. В расчётах выявлены две ветви неустойчивости, связанные с особенностями экспериментального профиля скорости на срезе сопла. Первая ветвь характеризуется большим коэффициентом нарастания возмущений и возникает из-за наличия точек перегиба в профиле скорости на периферии струи. Эта неустойчивость наблюдается только на начальном участке развития струи до x/D≈3. Вторая ветвь характеризуется меньшим коэффициентом нарастания возмущений, но сохраняется (при каждом фиксированном значении частоты ω) до больших удалений от сопла. Получено качественное и количественное согласие с имеющимися экспериментальными данными. Определены возможности каждого из использованных подходов. Разработан и отлажен алгоритм оптимизации геометрии параметрически заданного канала диффузора. В качестве минимизируемого функционала могут выступать близость к заранее заданному профилю скорости, или же характеристики неустойчивости профиля скорости, получаемого на выходе из формирующего устройства. Характеристики неустойчивости течения с данным профилем скорости рассчитываются в рамках невязкого или вязкого приближений, для чего на каждой итерации решаются задачи на собственные значения для уравнения Рэлея или Орра-Зоммерфельда, соответственно. Проведено исследование характера неустойчивости плоских струй вязкой жидкости, задаваемых профилем осевой скорости U(y). Численным решением задачи на собственные значения для уравнения Рэлея (в невязком приближении) или Орра-Зоммерфельда (с учётом вязкости) анализируется положение седловой точки функции ω(α) в комплексной плоскости, где ω – комплексная частота, α – комплексное волновое число. Продемонстрировано, как в невязком приближении происходит движение седловой точки при гладкой деформации профиля скорости и, соответственно, смена конвективного характера неустойчивости на абсолютный. Для профилей скорости, абсолютно неустойчивых в невязком приближении, показано, что с учётом вязкости неустойчивость снова становится конвективной, если число Рейнольдса ниже некоторого критического значения. Найдена зависимость этого критического числа Рейнольдса от профиля скорости. Проведены прямые численные исследования характера неустойчивости в струйных течениях. А именно, для заданного оптимизированного профиля решалась задача Коши для линеаризованных уравнений Навье-Стокса с заданным при t=0 локализованным возмущением. В результате развития этого возмущения обнаружена возможность распространения области, занятой растущим возмущением, вверх по потоку, то есть возникновение абсолютной неустойчивости в однонаправленной струе (т.е. в струе без противотока). В течение второго года выполнения проекта: опубликованы 3 статьи в научных журналах, индексируемых WoS/Scopus: Journal of Fluid Mechanics (Q1), Доклады РАН, Теплофизика высоких температур; опубликованы 3 статьи в сборниках, индексируемых Scopus: 2 статьи в Proceedings of the 8th International Conference on Fluid Flow, Heat and Mass Transfer (FFHMT’21), 1 статья в Proceedings of the 72nd International Astronautical Congress (IAC); Готовятся к публикации 5 статей, запланированные к подаче в журналы в течение ближайших 3 месяцев (из них три – в журналы Q1), по научным результатам, полученным в течение отчётного периода; сделаны 13 докладов на 8 российских и международных конференциях, в том числе 3 приглашённых доклада. Члены коллектива – аспирант Л.Р. Гареев и студент 6 курса К.Э. Абдульманов – заняли первые места среди докладов аспирантов и студентов, соответственно, на конференции-конкурса молодых учёных НИИ механики МГУ. | ||
3 | 1 января 2022 г.-31 декабря 2022 г. | Исследование развития возмущений в ламинарных струйных течениях и разработка научных основ новых способов управления такими течениями. |
Результаты этапа: I. Теоретическая часть. В двух постановках проведено исследование границы абсолютной/конвективной неустойчивости затопленных струй вязкой несжимаемой жидкости. Рассматривалось двухпараметрическое семейство профилей скорости, в которых задаётся скорость в точке перегиба и её градиент. В первой постановке с помощью критерия Бриггса (связанного с седловой точкой функции omega(alpha)) найдена граница абсолютной неустойчивости в пределе числа Рейнольдса, стремящегося к бесконечности. Для струй, абсолютно неустойчивых в этом пределе, найдены числа Рейнольдса, ниже которых неустойчивость снова становится конвективной. Таким образом, впервые теоретически показана возможность существования абсолютно неустойчивых плоских струй вязкой несжимаемой жидкости. Во второй постановке проводилось прямое численное моделирование (DNS) линеаризованных уравнений Навье-Стокса для плоской струй, в которую вносилось произвольное возмущений. Численно определялось поведение возмущения в фиксированной точке пространства. Показано, что при превышении критического числа Рейнольдса, полученного аналитически, возмущение в точке растёт, и ниже него – затухает. Таким образом, в численном моделировании подтверждена граница смены характера неустойчивости с абсолютной на конвективную. Факторами, способствующими абсолютной неустойчивости, являются понижение скорости и повышение градиента скорости в точке перегиба. В рамках этого же семейства профилей была проверена возможность возникновения абсолютной неустойчивости в круглых струях, и ни в одном из рассмотренных случаев с положением точки перегиба u>0.1 Umax абсолютную неустойчивость получить не удалось. Однако, более низкие скорости в точке перегиба делают практически невозможной реализацию абсолютной неустойчивости на практике, т.к. точка перегиба, находясь вблизи границы струи, быстро окажется “поднятой” из-за вязкого размывания струи. Таким образом, получено, что в плоских струях имеется гораздо большая предрасположенность струи к абсолютной неустойчивости, чем в круглых струях. II. Экспериментальная часть. Были проведены расчёты по проектированию внутреннего канала плоской струйной установки с целью понижения скорости в точки перегиба профиля струи. Форма канала параметризована сплайном; в результате решения стационарных уравнений Навье-Стокса рассчитывался выходной профиль скорости струи. Выбранный профиль канала был сконструирован в виде сменных вставок установки, напечатан на 3D-принтере и установлен в устройство, формирующее струю. Результаты измерения профиля скорости термоанемометром на выходе из струи показали удовлетворительное соответствие расчёту. Проведены эксперименты с визуализацией развития полученной струи. В отличие от наблюдавшегося ранее “классического” разрушения ламинарной струи через рост и снос вниз валиков Кельвина-Гельмгольца, в данной струе зафиксировано движение вверх по потоку возмущений, представляющих собой замкнутые вихревые зоны. Таким образом, можно утверждать, что впервые экспериментально сформирована абсолютно неустойчивая плоская затопленная струя. III. Расчётная часть Проведено численное исследование устойчивости и нелинейной эволюции возмущений в круглой струе. Воспроизведены условия лабораторного эксперимента при Re = 2850, проведенного ранее в НИИ механики МГУ им. М.В. Ломоносова (Зайко и др., ДАН, 2021; Gareev et al., JFM, 2022). Характерная особенность изучаемой струи – наличие трех точек перегиба во входном профиле скорости, что существенно определяет свойства течения. Линейная устойчивость исследуется с использованием двух подходов: квазипараллельного и пространственного. В первом изменение основного течения вдоль потока учитывается параметрически: задача устойчивости решается для серии профилей основного течения, соответствующих разным расстояниям от входа. В пространственном подходе линеаризованные уравнения Навье–Стокса для возмущений интегрируются с учетом изменения основного течения. Получено хорошее количественное согласие результатов двух подходов, а также согласие с результатами невязкой теории. Основное качественное свойство устойчивости течения состоит в существовании двух мод растущих возмущений: внутренней и внешней, отвечающих внутренней (r=0.5R) и паре внешних (r=0.85R) точек перегиба. Скорость роста возмущений внешней моды на порядок превышает аналогичную скорость роста возмущений внутренней моды. Однако эта неустойчивость пропадает уже на ранних стадиях развития струи вместе с исчезновением внешних точек перегиба в профиле скорости основного течения. Неустойчивость к внутренним модам сохраняется на всем протяжении развития струи. Проведены численные расчеты с целью объяснения и интерпретации результатов лабораторного эксперимента, в котором обнаружено изменение протяженности зоны перехода к турбулентности в струе под действием периодических во времени осесимметричных возмущений. В расчетах показано, что осесимметричные возмущения даже значительной начальной амплитуды не приводят к переходу к турбулентности. Колебания, вызываемые развитием таких возмущений, сохраняют регулярный, близкий к гармоническому, характер. Наблюдающийся в эксперименте переход к турбулентности может объясняться наличием неконтролируемых трехмерных возмущений, усиливающихся на фоне достаточно интенсивных искусственных возмущений. Проведенные в данной работе исследования подтверждают выдвинутую гипотезу. Входные возмущения в виде комбинации осесимметричного и малого трехмерного возмущения приводят к развитию турбулентности на расстоянии около 5 диаметров струи. Так же, как и в эксперименте, протяженность области перехода сокращается в тех случаях, когда частота входного возмущения близка к частоте наиболее быстро растущих возмущений по линейной теории. В описанной картине инициатором перехода служат трехмерные возмущения, а осесимметричные возмущения при достаточно большой амплитуде лишь ускоряют их рост. Отсюда следует, что более интенсивные начальные трехмерные возмущения могут быть способными обеспечивать еще более быстрый переход даже при отсутствии осесимметричной составляющей. Последнее предположение подтверждается расчетами с трехмерными начальными возмущениями конечной амплитуды. Длина зоны перехода сокращается в этом случае до 1–2 диаметров струи, что объясняет причину быстрого перехода, полученного в экспериментах. В течение третьего года выполнения проекта: • опубликованы 6 статей (9 статей с учётом квартиля Q1) в научных журналах, индексируемых WoS/Scopus: Journal of Fluid Mechanics (Q1, две статьи), Physics of Fluids (Q1), Доклады РАН, Известия РАН. МЖГ (2 статьи); • Сданы в печать 2 статьи в Journal of Fluid Mechanics (Q1), находящиеся в данный момент на втором круге рецензирования. • Сделаны 10 докладов на 7 российских и международных конференциях. |
Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".