Аннотация результатов, полученных в 2019 году
Разработаны программные реализации алгоритмов управления модельными течениями основанные на оптимальном управлении и управлении с использованием машинного обучения с подкреплением. Модифицирован пакет расчета гидродинамики TFlowsS для возможности моделирования течений в режиме с активным управлением с обратной связью. Пакет протестирован на расчетной задаче снижения силы аэродинамического сопротивления при обтекании цилиндра (2D LES, Re=100). Получены данные для коэффициента лобового сопротивления и подъемной силы для синусоидального во времени возмущения – вращательных колебаний стенки цилиндра вокруг своей оси. Показана высокая эффективность данного метода управления для высоких чисел Рейнольдса. По данным прямого численного моделирования распространения затопленной струи газ, вытекающей в газ другой плотности, был произведен анализ координаты распространения фронта струи в сравнении с литературными данными. Для целей разработки трёхжидкостной континуальной модели для концентрации, скорости и напряжений Рейнольдса инерционных частиц в пристенной зоне турбулентного течения, предложен простой вариант учёта неравновесности функции плотности вероятности скорости частицы с помощью подхода использующего двухжидкостную континуальную модель переноса напряжений Рейнольдса Л.И.Зайчика и известный в вычислительной гидродинамике (CFD) метод пристеночных функций для концентрации и второго момента скорости частиц, учитывающих их разделение на две группы: диффузионных и баллистических. Разработана математическая модель для описания теплообмена в нестационарных газокапельных струях при их взаимодействии с импактной поверхностью преграды. Проведено моделирование турбулентных импульсных импактных газокапельных струй, получены зависимости коэффициента теплообмена от расхода жидкой фазы при вариации продолжительности импульса, диаметра капель и их массовой доли. Показано согласие между данными моделирования и эксперимента в режимах с малым временем импульса и отличие при большом времени, установлены причины такого рассогласования. Показано, что поперечная пульсирующая струя является эффективным смесителем, поскольку она захватывает больше жидкости, чем обычная струя. Эксперименты по визуализации течения с круглыми турбулентными струями в поперечном потоке показали, что в структуре поперечной струи преобладает образование криволинейного сдвигового слоя, состоящего из отдельных вихревых петель вокруг струи в ближнем поле. В работе апробирован подход машинного обучения для решения задачи прогнозирования предотказных и аварийных состояний авиадвигателей. На основе баз данных эксперимента NASA проведено исследование математических моделей многоклассовой классификации режимов работы авиадвигателей и регрессии числа рабочих циклов до отказа. Протестированы различные семейства алгоритмов: от линейных до нейросетевых алгоритмов. Показана эффективность предлагаемого подхода, разработаны оптимальные прогнозные модели, доступные в виде фреймворка. Получены распределения турбулентных характеристик изотермического и реагирующего потока в модельной камере ГТУ. Проведена декомпозиция мгновенных полей скорости POD для различных режимов течения, которая указала на наличие глобальной моды, соответствующей прецессирующему вихревому ядру с сильными пульсациями давления. Исследована связь этой структуры и распределения CO и NOx в продуктах горения при значительном избытке воздуха и атмосферном давлении. В рамках проекта разработаны и усовершенствованы подходы для качественной и количественной диагностики распыла топлива форсунками ГТД. На образце топливной форсунки ГТД проведены предварительные исследования, получено описание структуры факела распыла и характеристик спрея. Кроме того, предложены способы модуляции подачи рабочей жидкости и воздействия на факел распыла для реализации методов управления на последующих этапах проекта. Разработана математическая модель и методика расчёта распыла, на основе комбинированного эйлерово-лагранжева подхода, основанная на моделировании потока жидкости внутри сопла и вблизи него на базе Эйлерового VOF метода, а образующуюся на удалении от сопла дисперсную среду – при помощи Лагранжевого метода дискретных частиц. Методика расчета распыла валидирована по данным прямого численного эксперимента. Показано, что использование LES модели турбулентности позволяет более корректно описывать распыл жидкой фазы по сравнению с гибридной RANS/LES моделью турбулентности. Наблюдается хорошее совпадение по характеристикам потока и дисперсной фазы, а также по границе раздела газа и жидкости. Проведено экспериментальное исследование восприимчивости обдуваемых газом пленок жидкости к наложенным гармоническим возмущениям. Исследован диапазон среднерасходных скоростей газа от 0 до 100 м/с и чисел Рейнольдса жидкости от 30 до 400, частота возбуждения варьировалась от 20 до 100 Гц. Показано, что наложенные возмущения существуют на протяжении некоторого расстояния от входа в канал и разрушаются вследствие слияний волн. Это расстояние слабо снижается с увеличением расхода жидкости и сильно - с увеличением скорости обдува. В экспериментах на немодифицированной модели направляющей лопатки с тупой задней кромкой (без устройства управления) была проведена высокоскоростная визуализация паровых каверн и измерены двумерные распределения скорости течения методом стробоскопической трассерной визуализации для различных условий течения по углу атаки и числу кавитации. В результате была составлена карта режимов кавитационного обтекания, проведен анализ зарождения и развития паровых каверн, их пространственной структуры и эволюции, исследовано влияние режима кавитационного обтекания на гидродинамические характеристики течения. На основе углубленного анализа литературы была определена оптимальная конструкция управляющего устройства в методе управления течением на основе периодического воздействия с помощью наложения внешнего электромагнитного поля (количество электродов, их геометрические размеры и взаимное расположение). Проведено численное моделирование в задаче обтекания гидропрофиля при высоких числах Рейнольдса для до- и кавитационных режимов. Была исследована динамика срыва кавитационного облака, частотные характеристики которого хорошо согласуются с данными визуализации из эксперимента. Предложена модификация сопла для завесного охлаждения плоской пластины, с каналом подвода охладителя содержащим изгиб и уступ. Наличие зоны рециркуляции в подобном сопле приводит к формированию когерентных структур, с возможностью управления частотой их схода. Расчеты показывают, что данная конфигурация может приводить к улучшениям тепловой эффективности охлаждения, и улучшает равномерность охлаждения в горизонтальной плоскости.

Аннотация результатов, полученных в 2020 году
С использованием современных методов экспериментальной диагностики, физико-математического и численного моделирования, в 2020 году в рамках проекта получены важные научные результаты о возможностях применения методов анализа больших массивов данных и машинного обучения для реализации эффективных методов управления тепломассопереносом в потоках, реализующихся в энергетическом оборудовании. В частности, акцент сделан на возможности управления развитием гидродинамической неустойчивости, интенсивностью тепломассопереноса и теплообмена, диспергирования и смесеобразования топлива, значительного влияния управляющих параметров на локальную скорость фазовых превращений (включая кавитацию) и горение. - Развита предложенная ранее схема управления течением в следе за цилиндром, основанная на наложении вращательных осцилляций цилиндра. Инструментом являлась совместная работа кода вычислительной гидродинамики T-Flows и нейронной сети OpenAI: baselines. Апробация такого метода управления с обратной связью была проведена на ламинарном течении при малых числах Рейнольдса (Re = 100). Определены критерии эффективного управления сопротивлением при помощи различных методов наложения внешних вращательных осцилляций. - Проведено прямое численное моделирование (DNS) импульсной струи воздуха, истекающей в пространство, заполненное гелием. Проанализировано влияние отношения плотностей смешивающихся газов на скорость распространения фронта струи. Изучено влияние начального входного профиля скорости на эволюцию течения. С использованием современных оптико-информационных методов диагностики выполнена серия экспериментов по исследованию влияния многомодового возбуждения на динамику вихревых структур на начальном участке затопленной струи жидкости. Частота периодического воздействия на двух гармониках выбрана так, чтобы основная гармоника соответствовала собственной частоте неустойчивости струи - значению числа Струхаля 0,5, а вторая частота обеспечивала низкочастотную модуляцию амплитуды сигнала. Кроме того, исследована эффективность многомодового возбуждения расхода струи топлива на возможность интенсификации смешения струи с окружающим воздухом и его влияние на процесс горения. - Экспериментально исследована динамика прямоточной газовой струи в поперечном потоке воздуха в условиях периодического возмущения расхода. Показано, что вертикальная пульсирующая струя в поперечном потоке является эффективным смесителем, поскольку захватывает больше окружающего газа, чем струя в стационарных условиях. Была изучена эффективность работы струйного осциллятора как инструмента управления потоком. Установлено, что в сносящем потоке он является эффективным устройством для внесения периодических возмущений в поток без использования механических подвижных деталей или модуляции расхода. В результате сравнения с данными исследования прямоточной струи в сносящем потоке, даже в условиях модуляции её расхода, струйный осциллятор демонстрирует значительно большую эффективность по интенсивности тепломассобмена между струей и основным потоком. Проведено моделирование воздействия устройства на базе барьерного разряда (ДБР-актуатора) на пристенную струю газовой завесы. Было обнаружено, что воздействие электрическими импульсами на выделенной частоте, близкой к одной из собственных частот струи, приводит к существенному уменьшению зоны рециркуляции и подавлению отрыва, при мощности воздействия незначительной по сравнению с потоком кинетической энергии струи. - Разработана гибридная модель переноса концентрации частиц, основанная на совместном применении диффузионно-инерционной модели в ядре потока и трёхжидкостной континуальной модели для концентрации и напряжений Рейнольдса частиц в пристенной зоне, учитывающей степенную зависимость для концентрации частиц у стенки, связанную с явлением кластеризации инерционных частиц. Предложено новое граничное условие типа пристеночной функции для рейнольдсовых напряжений частиц. Разработанная модель дает более реалистичные результаты для высокоинерционных частиц, адекватно воспроизводя пристенный пик концентрации частиц в вязком подслое у стенки. Выполнено численное моделирование влияния формы импульсов (прямоугольная, треугольная и гармоническая) на теплообмен в точке торможения импактной пульсирующей газокапельной струи (спрея). Для описания динамики и теплопереноса импульсного осесимметричного течения применялась система нестационарных осредненных по Рейнольдсу уравнений Навье–Стокса (URANS). Турбулентность несущей фазы описывалась с помощью модели переноса компонент рейнольдсовых напряжений. Показано, что в нестационарной импактной струе возможно как увеличение, так и подавление теплообмена по сравнению со стационарным течением для всех исследованных форм следования импульсов. Увеличение числа Рейнольдса вызывает снижение параметра интенсификации теплообмена, и данные для всех частот приближаются к стационарному режиму течения. Интенсификация теплообмена обусловлена вкладом испарительного охлаждения, который в свою очередь зависит от термического сопротивления жидкой пленки. Взаимодействие этих факторов, противоположно влияющих на пристенный тепломассоперенос, вызывает локальный экстремум в зависимости параметра эффективности от длительности импульса подачи спрея. Показано, что длительностях импульса около 10 мс наблюдается максимум коэффициента тепловой эффективности спрея. - Методами PIV/PLIF проведено детальное исследование смесеобразования и горения в модельной камере сгорания ГТУ (при значительном избытке воздуха) в условиях резонансных пульсаций давления, вызывающих модуляцию расхода топлива. Для идентификации когерентных структур в потоке данные обработаны методом POD. Проведена скоростная визуализация интенсивности хемилюминесценции (OH*) в модельной камере сгорания ГТУ с одновременной регистрацией пульсаций давления и измерением концентрации CO и NOx в продуктах горения на выходе из камеры. Установлено, что в камере сгорания без выходного сопла имеют место пульсации давления, вызванные прецессиией вихревого ядра потока, сопряженные с вращением крупномасштабных винтовых вихрей вокруг оси камеры сгорания. Данная мода неустойчивости соответствовала первой азимутальной моде и числу Струхаля 0,75. Аналогичные результаты были получены в результате математического моделирования методом LES: в отсутствие дежурного факела в потоке характерно наличие присоединенной к завихрителю рециркуляционной зоны и наличие прецессирующего вихревого ядра. В условиях резонанса пульсаций давления и расхода топлива, при наличии поджатия на выходе камеры сгорания с узким выходным соплом, в потоке были зафиксированы интенсивные пульсации интенсивности излучения OH*, соответствующие нулевой азимутальной моде. Они были вызваны модуляцией расхода топлива через топливные каналы и, соответственно, периодическому изменению соотношения топливо/окислитель, поступающей в камеру частично перемешанной смеси. Обратная связь в таком резонансом режиме горения обусловлена пульсациями тепловыделения и, соответственно пульсациями плотности, вызывающими пульсации давления. - В рамках проекта разработан стенд для исследования методов управления распылом жидкости за счёт акустических возмущений. Исследованы подходы для количественной панорамной диагностики дисперсного состава спрея на основе лазерно-индуцированной флуоресценции и упругого рассеяния излучения на каплях. На образце модельной топливной форсунки проведены предварительные исследования структуры факела распыла при воздействии акустических возмущений. Разработана методика нестационарного расчёта распыла, на основе комбинированного эйлерово-лагранжева подхода. Поток жидкости в сопле и начальной области струи описывается при помощи Эйлерова подхода (VOF), а образовавшиеся капли далее моделируются с использованием метода Лагранжа. Методика нестационарного расчета распыла с режимом модуляции расхода была валидирована по данным эксперимента распыла кольцевой струи воды с сопутствующими воздушными потоками. На основе метода теневой фотографии, высокоскоростной визуализации с использованием длиннофокусного макроскопического объектива изучены характеристики капель различных видов топлива, распыляемого струей перегретого водяного пара. Впервые проведена визуализация процесса распада разных типов жидкого углеводородного топлива при взаимодействии со сверхзвуковой газовой струей. Также впервые показана возможность эффективного диспергирования керосина паровой струей, такой метод может быть использован при проектировании перспективных газотурбинных двигателей. - Для дисперсно-кольцевого течения проведено исследование изменения записей толщины пленки вблизи режимной границы начала срыва капель. Показано, что в этой зоне резко повышается крутизна фронтов волн, что приводит к появлению мнимых пиков высокой интенсивности в записях толщины пленки, полученных методом лазерно-индуцированной флюоресценции. Таким образом, близость режимной границы можно оценить по поведению ЛИФ-записей до появления высоких волн. поведению ЛИФ-записей до появления высоких волн. - Выполнена серия расчетов течения жидкости при различных режимах обтекания гидропрофиля. Получены детальные согласующиеся с экспериментом результаты расчетов докавитационного режима течения при помощи метода RANS, а также для кавитационного режима методом LES с применением условного осреднения. Подробно изучена динамика потока при помощи временного спектрального анализа, а также POD разложения, позволяющего идентифицировать наиболее энергонесущие вихревые структуры. Идентифицированы низкочастотные пульсации, соответствующие срыву кавитационного облака с гидрокрыла. Оптимизирована конструкция модифицированного гидрокрыла со встроенным элементом управления, представляющим собой "сэндвич" из электродов и постоянных магнитов, созданы его трехмерная модель и лабораторный образец. С использованием ансамбля измеренных полей мгновенной скорости течения получены средние и пульсационные характеристики течения для различных кавитационных режимов. Анализ первичной информации о распределениях пульсаций скорости показал, что турбулентная структура течения существенно изменяется при появлении кавитации. Рассчитаны спектры мощности пульсаций давления для нестационарных режимов кавитационного обтекания и определены характерные частоты и амплитуды пульсаций давления в следе гидропрофиля, связанные с осцилляциями длины присоединенной каверны и отрывами крупномасштабных облачных каверн. Исследована эффективность управления динамикой кавитирующего потока на основе инжектирования щелевой струи на поверхности обтекаемого профиля. Показана возможность эффективного управления размером зоны рециркуляции, положением точки отрыва потока и размером кавитационного облака. Результаты исследований 2020 года опубликованы в ряде научных статей, в том числе в высокорейтинговых журналах первого квартиля.

Аннотация результатов, полученных в 2021 году
Проведен DNS расчет струи воздуха при числе Рейнольдса 5300, истекающей в пространство, заполненное углекислым газом. Скорость распространения фронта, а также мгновенные поля продольной скорости во фронте струи, сравнены для всех трех пар газов, включая “воздух-воздух” и “воздух-гелий”, рассмотренные за два предыдущих года выполнения проекта. Разработана и протестирована измерительная PIV система реального времени с регистрацией пар изображений VGA размера и одновременной GPU обработкой данных на частоте 25Гц. Разрешение получаемого поля скорости соответствует размеру входных изображений. Система может применяться для количественной визуализации течения в реальном времени, для целей обучения и как датчик скорости в системах активного управления потоками с характерными временами процессов не менее 0,1с. Созданная система протестирована на ряде течений в жидкости, в том числе в реальном эксперименте с пульсирующей затопленной струей, натекающей на плоскую горизонтальную нагреваемую стенку. Разработана континуальная модель для концентрации, скорости и напряжений Рейнольдса инерционных частиц в пристенной зоне для условий неупругого отскока и осаждения частиц. Модель учитывает степенную особенность концентрации частиц у стенки, связанную как с явлением кластеризации инерционных частиц, так и с явлением неупругого коллапса траекторий инерционных частиц. Разработанная модель имеет хорошее совпадение с данными стохастического лагранжева моделирования для турбулентного течения с взвешенными инерционными частицами в канале для условий неупругого отскока и осаждения частиц, адекватно воспроизводя пристенный пик концентрации частиц в вязком подслое у стенки. Проникновение импульсной струи в поперечный поток оказалось выше, чем в случае стационарной струи во всех случаях. Показано, что пульсирующая поперечная струя является эффективным смесителем. Управление частотой пульсаций и расходом импульсной струи показало возможность изменения (регулирования) подъема струи и эффективного перемешивания ее с поперечным потоком. В результате исследования влияния внешнего периодического возмущения на динамику поля скорости потока и положение фронта факела богатой смеси получена база данных PIV изображений и изображений спектрозональной визуализации пламени (изображения OH и CH). Применение метода главных компонент позволило получить количественную информацию об энергии когерентных пульсаций, вызванных периодическим движением тороидальных вихревых структур в потоке без горения и с горением. Впервые установлено, что для струйного факела с горением за фронтом пламени происходит резкое падение частоты следования вихрей практически в два раза. При этом скорость следования вихрей возрастает. Данный эффект ранее обнаружен не был и может быть объяснением того, что для повышения эффективности сжигания топлива в струйном факеле необходимо накладывать возмущения с частотой в два раза больше частоты основной гармоники без горения. С использованием изотропной низкорейнольдсовой k-эпсилон модели турбулентности Hwang and Lin (AIAA J., 1998), алгебраической модели турбулентности Girimaji (Theoret. Comput. Fluid Dynamics, 1996) и модели переноса рейнольдсовых напряжений Craft and Launder (AIAA J., 1992) выполнено исследование двухфазных нестационарных импактных струй. Было показано, что только модель Craft and Launder (AIAA J., 1992) дает неплохое согласие с данными измерений для одно- и двухфазного режимов. Были выполнены исследования импактных импульсных газокапельных струй с использованием эйлерова двухжидкостного и лагранжева траекторного подходов. Эйлерово описание движения и теплообмена в дисперсной фазе является более экономичным в вычислительном плане. В целом, оба подхода расчета дают близкие результаты. Основное отличие между эйлеровым и лагранжевым методами проявляется в профилях концентрации капель, и оно достигает 15%. Показано, что в центральной области теплообменника при натекании спрея на поверхность образуются застойные зоны со слабым разнонаправленным течением и области с интенсивным волновым течением, которое направлено от центра к периферии. Анализ данных о течении в центральной и периферийных областях поверхности теплообменника показывает, что в центральной области преобладает хаотичное сдвиговое волнообразование с разнонаправленным радиальным движением волн. С ростом средней толщины пленки возрастает амплитуда волн в горизонтальном движении к периферии поверхности, но скорость волн при этом остается низкой. В результате работ по анализу нестационарных явлений в модели камеры сгорания ГТУ создан и испытан оптический отсек диагностики поля скорости и формы фронта пламени с использованием методом PIV и OH PLIF. Стенд также обеспечивает контроль уровня пульсаций давления и содержание CO и NOx в продуктах сгорания. Проведена серия экспериментов при бедном сжигании метана с использованием модельного фронтового устройства для давлений до 4 атм. Акцент в исследовании сделан на анализе режимов нестационарного горения при сжигании газового топлива при значительном избытке воздуха. Отлажен метод численного моделирования методом LES. Установлено, что для нестационарных режимов горения с интенсивными пульсациями давления в камере сгорания характерны следующие процессы: периодический срыв (но не унос) фронта пламени и смещение вниз по потоку, накопление топлива в центральной зоне рециркуляции, резкое возрастание скорости распространения пламени, его продвижение вниз по потоку и быстрое сжигания топлива, накопившегося в центральной зоне рециркуляции. Последний процесс сопровождается быстрым ростом тепловыделения в камере сгорания, возрастанием давления и запиранием топливоподводящих каналов, что и обуславливает обратную связь системы. Аналогичный механизм термоакустических пульсаций имел место и для пламени при повышенном давлении. Влияние повышенного давления приводило к повышению сажеобразования в камере сгорания и повышению уровня NOx в продуктах горения. При помощи современных численных подходов в рамках метода крупных вихрей и моделей турбулентного горения была проведена серия параметрических расчетов модельной камеры сгорания при числе Рейнольдса 15000. Была проведена верификация результатов численного моделирования на данных эксперимента. Для изотермического и реагирующего случая наблюдается хорошее согласие с экспериментальными данными. Проанализировано распределение концентраций CO, CO2, NO и т.д. Показано, например, что CO2 образуется в зоне рециркуляции, в то время как высокая концентрация NO наблюдается уже в следе за рециркуляционной областью. Выполнена диагностика структуры факела распыла и характеристик спрея промышленной центробежной форсунки при воздействии акустического возмущения, вносимого посредством четырёх симметрично расположенных динамиков высокой мощности. Эксперименты проведены при суммарной мощности акустического воздействия до 1200 Вт в диапазоне частот возмущений от 100 Гц до 500 Гц. Показано, что акустическое возмущение в исследованных условиях не оказывает существенного влияния на осреднённые по времени характеристики факела. Идентифицированы собственные моды пульсаций в спрее, акустическое воздействие на частотах выделенных собственных мод способно приводить к их интенсификации. На основе разработанного ранее нестационарного эйлерово-лагранжева подхода для моделирования распыла, для промышленной форсунки получены результаты по структуре факела и дисперсности капель в режиме модуляции расхода в диапазоне частот возмущений от 20 Гц до 1000 Гц. Показано влияние модуляции расхода на динамику первичного распада струи. Проведены экспериментальные исследования характеристик горения керосина (марка ТС-1) при распылении высокоскоростной струёй перегретого водяного пара на лабораторной горелке мощностью 5-20 кВт. Впервые показана возможность сжигания керосина в струе водяного пара, достигнуто устойчивое воспламенение, найдены режимы, обеспечивающие низкие концентрации CO и NOx в конечных продуктах сгорания. Проведено сопоставление результатов для режимов с распылением керосина струей нагретого воздуха, показаны преимущества использования пара по сравнению с воздухом. Предложенный способ сжигания топлива в струе водяного пара может быть использован при проектировании малоэмиссионных газотурбинных двигателей. Проведено исследование эффекта наложенных периодических возмущений на волновую структуру пленки жидкости в дисперсно-кольцевом потоке. Управление дисперсно-кольцевым потоком при помощи наложения пульсаций расхода жидкости возможно в определенном диапазоне частот. Разработаны автоматические методы определения восприимчивости на основе спектрального подхода. Ширина диапазона восприимчивости составляет 15-20 Гц. Положение диапазона примерно соответствует частоте следования волн возмущения вдали от входа в канал. В этом диапазоне управление потоком заключается в значительной регуляризации структуры потока, повышении предсказуемости волновых характеристик, возможности изменять частоту, скорость и амплитуду волн возмущения, а также в повышении расхода дисперсной фазы. Проведено комплексное экспериментальное исследование влияния внешних возмущений на эволюцию кавитационных каверн. На основе измеренных пульсаций давления в следе за крылом были получены спектры мощности. Проведен первичный анализ спектров мощности пульсаций давления для нестационарных режимов кавитационного обтекания в случае внешнего возмущения. Продемонстрирована возможность задержки развития нестационарностей на облачных режимах обтекания модельного гидрокрыла. Разработан и апробирован новый метод извлечения информации о среднем распределении пара в кавитирующем потоке на основе обработанных PIV данных для жидкой фазы (Pervunin et al., Exp. Fluids, 2021). С использованием вихреразрешающих методов проведено численное моделирование кавитационного течения около гидропрофиля при больших числах Рейнольдса. Проведена верификация полученных данных с применением предложенной процедуры условного усреднения. Анализ безразмерного коэффициента подъемной силы, действующей на гидропрофиль, совместно с усреднением по характерным моментам времени показал, что коэффициент подъемной силы минимален при наличии наибольшего количества паровой фазы на поверхности гидропрофиля. В ходе реализации данного проекта в 2021 г. произведена модернизация разработанной ранее трехмерной численной модели расчета в эйлеровом приближении для расчета взаимодействия основного однофазного газового течения и двухфазного турбулентного пристенного потока, также выполнена реализация численного кода. Также было проведено моделирование тепловой эффективности газокапельной завесы при ее вдуве через наклонные цилиндрические отверстия в поперечную траншею. Получено значительное увеличение эффективности тепловой завесы при добавлении капель в поток охладителя (до 2-х раз в сравнении с подачей однофазного охладителя через отверстия в траншею). Основное увеличение тепловой эффективности пристенной двухфазной завесы наблюдается на расстоянии x/b = 15–20. Далее за счет испарения мелких капель происходит резкое снижение величины эффективности двухфазной завесы. Проведена оценка эффективности использования электрического поля для улучшения эффективности теплозащиты. В пакете моделирования сплошных сред OpenFoam была проведена модификация решателя для двухфазных течений по методу жидких границ для введения учета эффекта электростатического поля. Распределения электрического потенциала находились на каждой временной итерации солвера путем решения уравнения Пуассона (с учетом диэлектрической проницаемости жидкого охладителя). По найденным распределениям электрического потенциала строилось поле силы, действующей на охладитель (через поляризацию диэлектрика). Максимальное воздействие на диэлектрик (жидкий охладитель) будет оказано в области сильной неоднородности электростатического поля, при этом диэлектрик будет втягиваться в область максимальной концентрации силовых линий. С точки зрения увеличения эффективности теплозащиты данный эффект может быть использован для предотвращения отрыва пленки охладителя от охлаждаемой поверхности. Рассмотрена задача о пленочном охлаждении плоской пластины при подаче воды через цилиндрический канал под углом 35 градусов. Предложена модификация данной конфигурации, путем введения продольного выступа, идущего по нижней поверхности канала вдоль потока, с острой кромкой концентратором электростатического поля. При подаче постоянной разности потенциалов между верхней и нижней стенками канала электрическое поле, концентрирующееся на выступе стремится удержать струю охладителя на поверхности стенки, при этом не мешая ее продвижению в продольном направлении. Разработана установка по созданию газокапельной взвеси с целью управления потоком щелевой тепловой завесы с использованием электрического поля. Распределение дисперсной фазы по размеру исследовано с использованием лазерного допплеровского анализатора частиц. Проведены первые PIV-PLIF измерения по анализу изменения скорости движения частиц под действием постоянного электрического поля с целью оценки возможного влияния электрического поля на эффективность теплозащиты с учетом эффектов электрической поляризации капель.

Аннотация результатов, полученных в 2022 году
Проведена серия экспериментальных исследований по активному управлению течением за цилиндром со щелевым вдувом в корму цилиндра с использованием пульсирующей струи. Показано, что интегральное сопротивление изменяется при вариации давления, создающего пульсирующую струю и слабо зависит от частоты пульсаций. При постоянном вдуве наблюдались снижение сопротивления по сравнению с пульсирующими режимами, вследствие большего расхода воздуха через щель при равных давлениях. Для исследования конфигурации течения применялся скоростной PIV и дальнейший корреляционный анализ изображений частиц для оценки мгновенных полей скорости в отрывной зоне за цилиндром. Hеализовано активное управление потоком при обтекании цилиндра на основе обратной связи по датчику термоанемометра, установленному в следе за цилиндром напротив щели. Поиск в пространстве функций выполнялся методом генетического программирования из класса эволюционных алгоритмов, где из поколения в поколение переходят наиболее приспособленные законы управления. Реализованный метод управления позволил снизить сопротивление на более чем 20%. В ходе работы по проекту разработана модель обучаемой нейронной сети на базе архитектуры PWC-Net для оценки полей скорости по изображениям частиц. Представлены результаты оценки мгновенного поля скорости нейросетевым алгоритмом и алгоритмами оптического потока и точным, но медленным алгоритмом Dual TV L1. На основе вихреразрешающего численного моделирования, верифицированного на основе экспериментальных данных проведен анализ влияния входных условий, закрутки потока и соотношения топлива-окислитель на распределение концентрации CO и NOx в ближней области модельной камеры сгорания. Показана эффективность снижения вредных выбросов и уменьшении пульсаций давления при подавлении естественных гидродинамических мод в модельной камере сгорания ГТУ при горении метана со значительным избытком воздуха и повышенном давлении на основе многомодового воздействия и применения машинного обучения, в том числе с использованием искусственных нейронных сетей. Для случаев наиболее эффективного управления проанализированы динамика потока и гидродинамические моды на основе панорамных оптических измерений и вихреразрешающего численного моделирования с использованием методов понижения размерности стохастических динамических систем. В результате получен научно-технический задел по интеллектуальному управлению смешением и горением в камерах сгорания реальных газовых турбин. Реализована методика внесения азимутально-неравномерного акустического возмущения при помощи системы из четырёх динамиков с фазовым сдвигом сигнала для каждого динамика. Выполнена диагностика структуры факела распыла образца промышленно производимой центробежной форсунки при воздействии указанного возмущения в различных конфигурациях. Получены данные о распределении скорости капель в факеле распыла, геометрии факела и распределении относительной концентрации жидкости в факеле при воздействии азимутально неравномерного возмущения. Отмечается, что акустическое возмущение в исследованных условиях не оказывало существенного влияния на структуру факела и распределение скорости капель в нём. Отмечается определённое влияние возмущения на распределение концентрации жидкости и траекторию движения мелкой фракции распыла. В целом можно заключить, что в исследовавшихся конфигурациях влияние на спрей акустического возмущения является слабым, и проявляется только в осреднённых характеристиках. Проведены экспериментальные исследования влияния параметров перегретого водяного пара, при сжигании газообразных (пропан-бутановая смесь) и жидких углеводородов (керосин) в условиях паровой газификации, на концентрацию токсичных продуктов сгорания (CO, NOx). Получено, что при подаче пара обеспечивается значительное снижение содержания монооксида углерода (в 1.7 раз для жидкого топлива и 1.6 раза для газообразного топлива) и оксидов азота в продуктах сгорания (в 2 раза для жидкого топлива и 1.8 раза для газообразного топлива) по сравнению с подачей вместо пара нагретого воздуха. Полученные результаты демонстрируют, что метод сжигания газообразных и жидких углеводородов в условиях паровой газификации оказывается предпочтительным для повышения экологических показателей. Предложенный способ сжигания топлива в струе водяного пара может быть использован при проектировании малоэмиссионных газотурбинных двигателей. В данном разделе проекта проведено численное моделирование импульсной импактной газокапельной струи для различных форм импульсов (прямоугольная, синусоидальная и треугольная). Для описания динамики и теплопереноса нестационарного осесимметричного течения была использована система нестационарных осесимметричных осредненных по Рейнольдсу уравнений Навье–Стокса. Для описания турбулентности газа будут использованы модель переноса компонент рейнольдсовых напряжений, записанная с учетом двухфазности потока. Динамика и теплообмен в газокапельном потоке моделировались в эйлеровом приближении. Были проведены работы по созданию и обобщению базы данных численных исследований для нестационарных одно- и двухфазных турбулентных импактных струй. Проведен комплекс экспериментальных исследований нестационарного теплообмена при натекании импульсного двухфазного водно-воздушного спрея на вертикальную теплообменную поверхность. Проведена вариация параметров спрея в широких пределах: длительности импульса подачи жидкой фазы, частоты следования импульсов и скорости спутного потока воздуха, а также перепадов давления на жидкостных и газовых форсунках. Все опыты проведены в режиме постоянной температуры поверхности в 70 градусов Цельсия. Установлено, что спутный поток воздуха оказывает существенное влияние на интенсификацию теплообмена за счет турбулизации пристенного течения, разрушения жидкой пленки и возврата капель жидкости, отраженных от поверхности. Индивидуальный вклад импактный струи воздуха в суммарный теплообмен не велик и на балансе тепловой энергии практически не сказывается. Экспериментально продемонстрирована возможность расширения области существования волн возмущения и срыва капель в дисперсно-кольцевых потоках при помощи наложения низкочастотных пульсаций расхода жидкости. Осуществимо как уменьшение критических расходов жидкости при высоких скоростях газа, так и снижение критической скорости газа практически до нуля при высоких расходах жидкости. Проведено исследование истечения струи, генерируемой струйным осциллятором, в щелевой канал. Результаты исследования показали, что при равных значениях высоты канала и ширины горловины выходного сопла истечение струи в щелевой канал сопровождается формированием крупномасштабных квазидвумерных вихревых структур на основной частоте работы струйного осциллятора. Впервые при помощи 3D PIV измерений обнаружены крупномасштабные продольные вихревые структуры. Показано, что динамика этих структур связана с меандрирующим характером течения струи и развитием крупномасштабных квазидвумерных вихревых структур. Получена экспериментальная база средних и турбулентных характеристик течения для кавитационного обтекания гидрокрыла при наложенных внешних согласованных возмущениях. Проведен анализ результатов всего комплекса экспериментальных данных, полученных для модифицированного и оригинального гидрокрыла. Показана эффективность подавления нестационарности и кавитации реализованным методом. Для экспериментальных данных и данных численного моделирования кавитации проведен анализ POD и спектральных характеристик кавитационного течения около гидропрофиля. Для задачи завесного охлаждения пластины проведена многопараметрическая оптимизация формы поверхности на основе метода сопряженных полей и вихреразрешающего численного моделирования для получения максимальной эффективности теплозащиты. Получены результаты моделирования влияния нестационарности входной концентрации частиц на пристенную аккумуляцию частиц в турбулентном течении в плоском канале.