В рамках проекта № 20-31-70034 Российского фонда фундаментальных исследований комплексно (экспериментально, численно и аналитически) изучена пондеромоторная магнитная сила плавучести Fm, действующая на твердую немагнитную или намагничивающуюся (суперпарамагнитную) сферу в цилиндрическом контейнере со стороны магнитной жидкости (МЖ), намагниченной внешним однородным магнитным полем H, ориентированным вдоль оси цилиндра.

Read more: Плавание твердых ферро-, пара- и диамагнитных тел в немагнитном контейнере, заполненном магнитной...

Решением редколлегии журнала «Известия РАН. Механика твердого тела» от 20 декабря 2021 года статья «Термодинамика упругого материала с релаксирующим потоком тепла» (Известия РАН. МТТ. 2020. № 4. С. 152-157.) признана лучшей статьей за 2020 г. В статье дается вывод уравнений, позволяющих исследовать с помощью вычислительного моделирования на компьютере интересное явление. Это появление в материале тепловой волны. Вывод уравнений осуществлен с из первого и второго законов термодинамики и требования сохранения их вида в любой инерциальной системе отсчета. В 2022 году осуществлено применение теории к конкретному явлению.

Read more: Тепловая волна

            В специальном выпуске “Advances in Computational Electromagnetics” журнала Magnetochemistry (№7 за 2021 г.) вышла статья сотрудников ИМСС УрО РАН А.А. Рогового, О.В. Столбова и О.С. Столбовой “The Microstructural Model of the Ferromagnetic Material Behavior in an External Magnetic Field”, получившая статусы Feature Papers (тематические статьи, представляющие собой самые передовые исследования со значительным потенциалом для оказания большого влияния в этой области) и Editor’s Choice (выбранные редактором статьи, недавно опубликованные в журнале, которые, по их мнению, будут особенно интересны авторам или важны в этой области).

Read more: Микроструктурное моделирование сплава Гейслера

В рамках проекта Российского научного фонда (конкурс по лабораториям мирового уровня) № 21-79-30041 «Закономерности критичности в материалах с дефектами, разработка подходов по мониторингу и прогнозированию ресурса при широкодиапазонном силовом и энергетическом воздействии (приложения в авиационном моторостроении), 2021-2024 гг.» совместно с АО «ОДК-Авиадвигатель» (индустриальный партнер) выполняется широкий класс экспериментальных и теоретических исследований по разработке методов оценки ресурса материалов и конструкций авиационного моторостроения применительно к перспективным отечественным двигателям ПД-14, ПД-35.

Read more: Прогнозирование ресурса материалов и элементов конструкций авиационных двигателей нового поколения

Договоры в рамках работ по гранту

договор 01_06 на перенос кабтрассы

договор 02_07 на подготовку площадкидоговор 27_10_2015 на строительство объекта

договор 986 на инженерные изыскания

Заявления

заявление на раздел земельного участка

заявления на получение разрешения на строительство объекта №1

заявления на получение разрешения на строительство объекта №2

Письма входящие

№01_5075_08 О рассмотрении обращения

№007_18.1_14_АС_1847 О согласии на строительство

№742_41 Об установлении адреса

№И-22-01-05_1-374 Об отказе в выдаче разрешения на строительство

Письма исходящие

№16344_01-247 С просьбой о поддержке проекта

№16344_01-306 О согласии на строительство

№16344_01_115 Предложения по использованию образуемых земельных участков

№16344_01_308 Об ответах на замечания к проектной документации и получении разрешения на строительство

Приказы

приказ 31 от 7.08.2015

приказ 32 от 7.08.2015

Проектно-сметная документация

РНФ 14-29-00172 Инженерные изыскания (RAR)

РНФ 14-29-00172 Рабочая документация (RAR)

РНФ 14-29-00172 Проектная документация (RAR)

РНФ 14-29-00172 Сметы (RAR)

РНФ 14-29-00172 ТЗ на проектирование объекта

Разрешительная документация

Приказ №32 (лист 1)

Приказ №32 (лист 2)

Приказ №32 (лист 3)

ЕГРН 1_3

ЕГРН 2_3

ЕГРН 3_3

РНФ 14-29-00172 Градостроительный план земельного участка

РНФ 14-29-00172 Разрешение на строительство объекта

доверенность на Гусева Г.Н.

разрешение на строительство 1_3

разрешение на строительство 2_3

разрешение на строительство 3_3

свидетельство СРО на проектирование

2022 07 14 11 21 55 828

 2022 07 14 11 21 56 014

2022 07 14 11 21 56 102

2022 07 14 11 21 56 182

2022 07 14 11 21 56 238

2022 07 14 11 21 56 294

2022 07 14 11 21 56 368

2022 07 14 11 21 56 434

2022 07 14 11 21 56 491

2022 07 14 11 21 56 549

2022 07 14 11 21 56 602

2022 07 14 11 21 56 663

2022 07 14 11 21 56 718

2022 07 14 11 21 56 775

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

Этап 1

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

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

Выполнены эксперименты на сепарационной ячейке, представляющей собой вертикально расположенную плоскую прямоугольную кювету из оргстекла с рядом аккумулирующих перегородок, задерживающих сепарированную примесь в области между ними, с двумя электродами для пропускания тока вдоль кюветы. Ячейка помещена в магнитное поле. Измерения проводились при помощи ультразвукового доплеровского анемометра (УДА). Для верификации расчетов гидродинамических процессов было проведено экспериментальное исследование процесса МГД-воздействия на включения. Средой являлся раствор электролита, а второй фазой и трассерами - частицы сажи. В верификационном эксперименте была специально создана неоднородность электромагнитной силы, которая привела к генерации достаточно интенсивного вихревого течения. Результаты расчетов демонстрируют неплохое согласие с экспериментом. В экспериментах по процессу электромагнитного осаждения, пространстве между перегородками возникают паразитные вторичные течения. Интенсивность таких течений невелика, но этого оказывается достаточно для постепенного вымывания примеси из области аккумуляции и затруднения анализа изучаемого процесса. Все указанные трудности приводят к тому, что моделируемый процесс хоть и демонстрирует свою работоспособность, но имеет низкую эффективность на электролите. Для работы с химически активными расплавами жидких металлов наиболее подходит индукционный метод измерений электропроводности. В ходе выполнения проекта был разработан и изготовлен портативный индукционный датчик электропроводности металлов. Было определено, что датчик способен различать между собой образцы металлов и сплавов с разной электропроводностью, а также регистрировать изменение концентрации меди в оловянно-свинцовых цилиндрах. Таким образом, выполнено экспериментальное исследование индукционной методики определения электропроводности на твердых цилиндрических образцах металлов и их сплавов с известными электрическими проводимостями.

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

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

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

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

Этап 2

С использованием комбинированной лагранжевой-эйлеровой модели выполнены расчеты сепарации в цилиндрических каналах различного диаметра (А1.2). Построены зависимости коэффициента сепарации от частоты внешнего переменного магнитного поля для частиц различного размера при течении в каналах различного диаметра. Также определена зависимость коэффициента сепарации от расхода при фиксированной частоте внешнего магнитного поля. С использованием многоскоростной модели проведены расчеты эволюции распределения примеси во вращающемся течении с целью дальнейшей верификации. Определены зависимости коэффициента гомогенизации от времени для различных топологий электромагнитной силы.

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

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

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

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

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

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

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

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

Исследован кондукционный способ генерации электромагнитной силы (B4.1). С помощью математического моделирования получены величины интенсивности течения, формы профиля и пульсационные характеристики для различных значений интенсивности электромагнитной силы. Предложен и реализован кондукционный способ измерения зависимости эффективной проводимости ограниченного объёма металлического расплава. Полученная зависимость эффективной проводимости от концентрации примеси имеет существенно нелинейный характер в области малых концентраций.

Выполнена численная и экспериментальная проверка кондукционного способа регистрации скорости на реальной среде путем сравнения с имеющимся эталонным расходомером (В5.2). Калибровка производилась на натриевом стенде в температурном диапазоне 140-210 градусов. Среднее значение относительной погрешности составило 1.3%, максимальное – 3.2%. Расчеты и эксперименты, в целом, показали, что наличие растворенных примесей не оказывает существенного влияния на показания расхода в данном способе.

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

Для лабораторной установки по электромагнитному разделению фаз в электропроводной среде изготовлен насос, обеспечивающий циркуляцию жидкого металла (B6.3). Это индукционный электромагнитный насос с цилиндрическим коаксиальным каналом. Насос имеет 12 катушек. Расчетное давление в стопорном режиме – 0.25 МПа; расход жидкого металла – 3 куб.м/ч.

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

Этап 3

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

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

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

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

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

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

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

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

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

В ходе выполнения проекта опубликовано 8 статей и сделан 31 доклад на 14 конференциях.

 

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

Руководитель проекта:

Фрик Петр Готлобович, д.ф.-м.н., зав. лаб. ИМСС УрО РАН (г. Пермь)

Краткий отчет по результатам выполнения второго года проекта РНФ №22-61-00098:
Все задачи проекта разбиты на три блока. Первый блок направлен на исследование блокирующих явлений в умеренных широтах, которые могут быть причиной опасных метеорологических явлений, связанных с волнами тепла или холода. Второй блок задач направлен на исследование влияния крупных плавающих теплоизолирующих объектов на формирование конвективных течений и распределение теплопотоков, с приложением к таким природным объектам, как массивы льда, тектонические плиты, облачные кластеры. Третий блок задач направлен на изучение формирования вынужденных и свободноконвективных течений в области городской агломерации, в том числе городского острова тепла, и их влияния на условия комфортности жизни.
Среди основных результатов второго этапа выполнения проекта можно отметить результаты лабораторного и численного моделирования общей циркуляции атмосферы при изменении условий в полярной области. Показано, что значительная вариация мощности охлаждения и граничных условий приводит к количественным изменениям в структуре и интенсивности бароклинных волн. Установлено, что отклонения энергии основных мод для разных случаев могут достигать 25 %, а локальная вариация охлаждения оказывает сильное влияние на структуру аналога полярной ячейки. Размер и граничные условия на поверхности холодильника играют решающую роль в структуре и интенсивности аналога полярной ячейки. Переход от сильного локального охлаждения в центре и равномерного охлаждения на остальной части свободной поверхности к равномерному охлаждению по всей поверхности приводит к ослаблению зональных потоков. Главным и довольно неожиданным результатом моделирования потепления Арктики является заметная трансформация поля средней температуры, а именно: центральная область и большая часть нижнего слоя становятся теплее, в то время как большая часть верхнего слоя и периферийная (экваториальная) часть нижнего слоя становятся холоднее.
С помощью стандартного метеорологического пакета WRF-ARW версии 4.2.2 проведено моделирование предельных состояний земной атмосферы для выявления роли отдельных факторов. Предельные ситуации — глобальный океан или глобальная пустыня показали роль фазовых переходов в атмосфере Земли. В первом случае вследствие парникового эффекта циркуляция не выражена из-за отсутствия значимого межширотного градиента. Во втором случае песчаная поверхность, в отличие от океана, не может аккумулировать большое количество тепла, поэтому воздух у поверхности земли быстро охлаждается в зимнем полушарии и быстро нагревается в летнем, из-за чего пассатных течений не возникает. Комбинация океанических и пустынных зон (все однородны по меридиану) приводит к тому, что крупномасштабная циркуляция атмосферы приобретает черты, близкие к земной, что позволяет рассматривать такую конфигурацию базовой для дальнейших численных экспериментов. Проведенные эксперименты показали ведущую роль поглощенной части солнечной радиации (инсоляции) и вклада водяного пара (парниковый эффект, выделение скрытого тепла) при формировании крупномасштабных атмосферных течений.
На основе лабораторных и численных экспериментов проведен систематический анализ динамики свободноплавающего вдоль вытянутой конвективной ячейки тела и структуры возникающих конвективных течений. Показано, что поведение системы определяется числом Релея, геометрическими параметрами (отношением толщины слоя к его длине и отношением размеров тела к длине слоя), и глубиной погружения тела и демонстрирует разнообразные режимы поведения тела, от покоя и строго периодических колебаний до стохастических. Режим регулярных колебаний возникает в полостях не слишком вытянутых, при умеренных нагревах и вертикальном положении диска вблизи теплового пограничного слоя. Этот режим ограничен в пространстве параметров, но не имеет резких границ. Периодические движения возникают в определенном диапазоне отношений размера тела к длине полости. С ростом числа Релея движения диска становятся все менее регулярными. При удалении диска от теплообменника режим регулярных колебаний не регистрируется вовсе. При малых нагревах диск остается у одной из стенок, а по мере роста числа Релея начинает совершать колебательные движения вблизи одной стенки. При дальнейшем росте числа Релея диск иногда совершает экскурсы вплоть до противоположной стенки, но что интересно, асимметрия движения сохраняется на протяжении всей реализации, что позволяет предположить существование доминирующей крупномасштабной циркуляции одного знака. При еще больших нагревах движение становится хаотическим, без какой-либо асимметрии.
В лабораторных экспериментах по исследованию динамики и структуры течений в конвективной ячейке со свободной верхней поверхностью и радиационным подогревом нижней поверхности, в которой на фиксированной высоте свободно перемещается в одном направлении теплоизолирующая пластина, показано, что динамика тела существенно зависит от оптических свойств его поверхности. Диск со светоотражающей поверхностью демонстрирует регулярные квазипериодические колебания, перемещаясь от одного края полости к другому. Период колебаний существенно зависит от геометрии полости и высоты фиксации диска. Если же диск поглощает или пропускает падающий на него свет, то он перемещается к одной из боковых стенок и остается там.
Усовершенствована модель переноса излучения в городской среде в рамках программного комплекса «SigmaEco», которая теперь учитывает рассеяние солнечного и теплового излучения в дисперсной среде, а также лучистый теплообмен с зелеными насаждениями. Разработана математическая модель, описывающая влияния зеленых насаждений на тепломассообменные процессы в городской среде в приближении пористой среды и на процессы образования и переноса паров воды. Отладка модели м тестирование было выполнено на данных наблюдения в г. Красноярск.
В результате экспериментов в аэродинамической трубе с набором макетов морфотипов городской застройки с применением скоростных оптических методов, в том числе трассерной визуализации верифицированы численные модели, определены условия масштабируемости задачи, подтверждены возможности физического моделирования формирования аэрационного режима и условий ветрового комфорта городской среды в лаборатории.
Методика зонирования по уровням пешеходного ветрового комфорта, была реализована на реальных жилых микрорайонах города Красноярска. Рассматривалась городская застройка преимущественно смешанного морфотипа, где здания имеют разную геометрическую форму и высоту. В ходе проведения численных исследований было определено, что на территориях дворовых пространств преобладают благоприятные условия для пребывания человека, скорости ветра не превышают 2,5 м/с и относятся к А классу комфортности среды. Рекомендовано расположение зон отдыха, детских площадок, а также расположение общественных террас. В области расположения высотных зданий наблюдается увеличение скоростей до значений, считающихся дискомфортными и опасными для жизнедеятельности, скоростной режим характеризуется классом С. В данных областях необходимо осуществлять ветрозащитные мероприятия.

Публикации по результатам проекта, вышедшие в 2023 г.:

  1.  Фрик П., Филимонов С., Гаврилов А., Попова Е., Сухановский А., Васильев А. (Frick P., Filimonov S., Gavrilov A., Popova E., Sukhanovskii A., Vasiliev A.) Rayleigh-Benard convection with immersed floating body Journal of Fluid Mechanics (2023 г.) WOS SCOPUS Q1 иные ББД
  2.  Васильев А. Ю. , Попова Е.Н., Фрик П.Г., Сухановский А.Н. (Vasiliev A. Yu. , Popova E. N., Frick P. G., Sukhanovskii A. N.) Drift of a Free-floating Body in a Convective Layer Heated by Radiation Journal of Siberian Federal University. Mathematics and Physics. (2023 г.) WOS SCOPUS RSCI РИНЦ
  3.  Дектерев Д. А., Лобасов А.С., Мешкова В.Д., Литвинцев К.Ю., Дектерев Ар.А., Дектерев А.А. (Dekterev D.A., Lobasov A.S., Meshkova V.D., Litvintsev K.Yu., Dekterev Ar.A., Dekterev A.A.) Анализ влияния масштабного фактора на результаты моделирования обтекания зданий Теплофизика и аэромеханика (2023 г.) WOS SCOPUS RSCI РИНЦ
  4.  Литвинцев К.Ю., Дектерев А.А., Мешкова В.Д., Филимонов С.А. (Litvintsev K.Yu., Dekterev A.A., Meshkova V.D., Filimonov S.A.) Влияние излучения на формирование ветрового и температурного режима в городской среде Теплофизика и аэромеханика (2023 г.) WOS SCOPUS RSCI РИНЦ
  5.  Филимонов С.А., Гаврилов А.А., Фрик П.Г., Сухановский А.Н., Васильев А.Ю. (Filimonov S.A., Gavrilov A.A., Frick P.G., Sukhanovsky A.N., Vasiliev A.Yu.) Моделирование движения погруженной пластины в развитом свободно-конвективном слое Известия вузов. Физика. (2023 г.) WOS SCOPUS RSCI РИНЦ
  6.  Калинин Н.А., Ветров А.Л. (Kalinin N.A., Vetrov A.L.) Индексы комфортности климата в Перми и Красноярске за период 1991–2020 гг. Геосферные исследования (2023 г.) WOS SCOPUS РИНЦ
  7.  Быков А.В., Ветров А.Л., Фрик П.Г., Сухановский А.Н., Калинин Н.А., Степанов Р.А. (Bykov, A.V., Vetrov, A.L., Frick, P.G., Sukhanovsky, A.N., Kalinin, N.A., Stepanov, R.A.) Численное моделирование предельных состояний планетарной атмосферы Географический вестник (2023 г.) RSCI РИНЦ
  8.  Мешкова В., Дектерев А., Литвинцев К., Филимонов С. (Meshkova V., Dekterev A., Litvintsev K., Filimonov S.) Current approaches to studying the level of pedestrian comfort in urban development E3S Web of Conferences (2023 г.) SCOPUS иные ББД

Результаты проекта были широко представлены на российских и зарубежных конференциях в 2023 году: :
   XXIII Зимняя Школа по механике сплошных сред, Пермь, 13-17 февраля 2023:
Организована и проведена подсекция «Динамическая метеорология», на которой участниками проекта представлены устные доклады:
1) Литвинцев К.Ю., Гаврилов А.А., Филимонов С.А. Влияние излучения на формирование аэрационного режима в городской застройке
2) Дектерев Д.А., Мешкова В.Д. Применение методов трассерной визуализации для исследования аэродинамики модельных групп зданий
3) Евграфова А. В., Ветров А. Л. Анализ температуры поверхности городов Пермь и Красноярск на основе спутниковых данных
4) Васильев А.Ю., Попова Е.Н., Фрик П.Г., Сухановский А.Н. Динамика свободноплавающего тела в условии естественной конвекции при радиационном нагреве
5) Филимонов С.А., Гаврилов А.А., Литвинцев К.Ю., Васильев А.Ю., Попова Е.Н. Движение погруженной пластины в развитом свободно конвективном слое
6) Сухановский А.Н., Попова Е.Н., Васильев А.Ю. Лабораторное моделирование арктического потепления
7) Степанов Р.А. Пространственная неоднородность атмосферной турбулентности по данным реанализа NCEP-CFSR
8) Калинин Н.А., Быков А.В., Ветров А.Л. Прогноз степени суровости погоды в феврале 2019 года в городе Красноярск по данным модели WRF
9) Степанов Р.А., Соколов Д.Д., Фрик П.Г. Согласованность климатических изменений различных временных масштабов в Центральной Англии и Гренландии
Всероссийская научная конференция с международным участием «Енисейская Теплофизика», Красноярск, 28-31 марта 2023 :
10) Фрик П.Г. «Особенности турбулентной конвекции жидкого металла в замкнутых объемах» - Приглашенная лекция.
11) Васильев А., Попова Е., Сухановский А., Фрик П. Динамика и теплообмен в конвективной системе с плавающим теплоизолятором при радиационном нагреве подстилающей поверхности - Устный доклад. (Тезисы докладов I      Всероссийской научной конференции с международным участием «Енисейская теплофизика». Красноярск, 2023. С.261-263.)
12) Гаврилов А. А., Платонов Д. В., Сухановский А.Н. «Прямое численное моделирование лабораторной модели общей циркуляции атмосферы» - Устный доклад (Тезисы докладов I Всероссийской научной конференции с международным участием «Енисейская теплофизика». Красноярск, 2023. С.264-266.)
13) Мешков К. Н., Дектерев А. А., Мешкова В.Д., Шульженко П.Д., «Сбор метеоданных из открытых телекоммуникационных источников и методология их использования для оценки температурного уровня комфорта городской среды» - Устный доклад (Тезисы докладов I Всероссийской научной конференции с международным участием «Енисейская теплофизика». Красноярск, 2023. С.282-284.)
14) Литвинцев К.Ю., Дектерев А.А., Филимонов А.А. Интенсификация процессов тепломассопереноса в городском квартале под действием солнечного излучения - Устный доклад (Тезисы докладов I Всероссийской научной конференции с международным участием «Енисейская теплофизика». Красноярск, 2023. С. 118-120.)
15) Дектерев Д.А., Лобасов А.С., Мешкова В.Д., Литвинцев К.Ю., Дектерев Ар.А., Дектерев А.А. Анализ влияния масштабного фактора на результаты моделирования обтекания зданий - Устный доклад. (Тезисы докладов I Всероссийской научной конференции с международным участием «Енисейская теплофизика». Красноярск, 2023. С. 99-101.) В рамках конференции "Енисейская теплофизика" организован и проведен круглый стол «ЭКОЛОГИЧЕСКИЕ ПРОБЛЕМЫ КРУПНЫХ ГОРОДСКИХ АГЛОМЕРАЦИЙ» с двумя докладами (из 7) участников проекта:
16) Ветров А.Л. и др. «Региональная модель WRF как инструмент изучения городской атмосферы»
17) Литвинцев К.Ю. и др. «Микромасштабное моделирование атмосферы города»
VII Международный симпозиум по региональной экономике (REC-2023) «Города нового времени: система GLASS», 27-29 июня 2023 г. (Екатеринбург)
18) Мешкова В.Д., Дектерев А.А., Дектерев Д.А., Литвинцев К.Ю., Филимонов С.А.Современные подходы исследования уровня пешеходного комфорта в условиях городской застройки - устный доклад.

   XIII всероссийский съезд по теоретической и прикладной механике, 21-25 августа 2023 года, Санкт-Петербург.
19) Фрик П.Г., Степанов Р.А., Шестаков А.В. О реализуемости режима Обухова-Болджиано в среде с неоднородной плотностью и неконсервативных каскадах в развитой турбулентности - Устный доклад
20) Сухановский А.Н., Васильев А.Ю., Попова Е.Н., Гаврилов А. Влияние локального охлаждения на формирование бароклинных волн во вращающемся слое - устный доклад
 VIII Всероссийская научная конференция "Теплофизика и физическая гидродинамика"

  Научная молодежная школа "Теплофизика и физическая гидродинамика: современные вызовы" (ТФГСВ2023), 4-10 сентября (Махачкала)  http://www.itp.nsc.ru/tph/2023/index.html
21) Дектерев Д.А., Литвинцев К.Ю., Мешкова В.Д., Дектрев Ар.А. «Анализ применимости численных моделей для моделирования ветровой комфортности микрорайонов» - устный доклад IX Всероссийская конференция "Пермские гидродинамические научные чтения", 4-6 октября 2023 (Пермь)
22) Фрик П.Г., Васильев А.Ю., Попова Е.Н., Сухановский А.Н., Филимонов С.А., Гаврилов А.А. Дрейф погруженного в жидкость тела при конвекции в прямоугольных полостях - устный доклад (Аннотации докладов. С.42.)
23) Васильев А.Ю., Попова Е.Н., Сухановский А.Н., Фрик П.Г. Динамика свободноплавающих тел в слое жидкости с радиационным нагревом - устный доклад (Аннотации докладов. С.12.)
24) Быков А.В., Сухановский А.Н., Калинин Н.А., Степанов Р.А., Фрик П.Г., Ветров А.Л. Моделирование предельных состояний планетарной атмосферы с использованием пакета WRF-ARW - устный доклад (Аннотации докладов. С.12.)
25) Сухановский А.Н., Гаврилов А.А., Васильев А.Ю., Попова Е.Н. Моделирование арктического потепления в лабораторной модели общей циркуляции атмосферы - устный доклад (Аннотации докладов. С.40.) 9-th International Conference on Rayleigh-Benard turbulence, 16-20 Oct.2023. Xi'an. China
26) Peter Frick “Large Scale Dynamics of Turbulent Convection in Closed Systems” // 9-th International Conference on Rayleigh-Benard turbulence, 16-20 Oct.2023. Xi'an. China (Keynote Talk) - Приглашенная лекция
27) Andrei Sukhanovskii, Andrei Gavrilov, Elena Popova, Andrei Vasiliev “Modeling of baroclinic waves with geophysical applications” // 9-th International Conference on Rayleigh-Benard turbulence, 16-20 Oct.2023. Xi'an. China - устный доклад

   Всероссийская конференция «XXXIX Сибирский теплофизический семинар» посвящённой 90-летию академика А.К. Реброва» 28-31 августа 2023 г., Новосибирск.
28) Шульженко П.Д., В.Д. Мешкова, Дектерев А.А. Численное исследование формирования ветровых зон в результате обтекания моделей разной геометрической формы - устный доклад
29) Шульженко П.Д., Мешкова В.Д., Дектерев А.А., Мешков К.Н. Численное исследование аэрационного режима жилой застройки - устный доклад

   Международная научная конференция «Безопасность жизнедеятельности и климатические риски развития территории Енисейской Сибири», 16-19 октября 2023г., Красноярск.
30) Шульженко П.Д., Мешкова В.Д., Дектерев А.А., Мешков К.Н. Численное исследование аэрационного режима городской застройки периметрального типа - устный доклад
31) Мешкова В.Д., Шульженко П.Д., Дектерев А.А., Литвинцев К.Ю., Мешков К.Н. Численное исследование формирования аэрационного режима и условий ветрового комфорта в разных морфотипах городской застройки - устный доклад.

Краткий отчет по результатам выполнения третьего года проекта РНФ №22-61-00098:

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

Блок 1. "Циркуляция и блокирование в умеренных широтах":
По результатам лабораторного и численного моделирования предложен и обоснован один из возможных физических механизмов усиления потепления в Арктике, основанный на изменениях радиальных (меридиональных) потоков тепла. Ослабление и смещение вверх лабораторного аналога полярной ячейки приводит к значительному уменьшению отрицательного теплового потока вблизи дна, что и вызывает повышение температуры в придонном слое. Наблюдаемые для лабораторного аналога тренды в изменении температуры качественно согласуются с данными натурных наблюдений (реанализа) и ансамблевых расчетов при помощи климатических моделей для арктического региона.
На основе пакета WRF-ARW реализованы численные эксперименты по динамике общей циркуляции и бароклинных волн в атмосфере идеализированной планеты земного типа. Модель, с одной стороны, точно воспроизводит реальные характеристики Земли, а с другой - предельно упрощает конфигурацию континентов и океанов, полностью игнорируя топологию планеты. Показано, что относительно небольшой экваториальный океан на фоне плоской пустынной планеты достаточен для воспроизведения общей циркуляции атмосферы и динамики среднеширотных вихревых структур планетарного и синоптического масштабов, реалистичных для Земли.
В части анализа пространственно-временной структуры вариаций климата проведено исследование согласованности климатических изменений в разных временных масштабах по обе стороны от Атлантики. Было установлено, что коротких масштабах (5-10 лет) колебания температуры в двух регионах относительно согласованы, но вариации на масштабах 40-50 лет идут со значительным сдвигом по фазе. Показано, что, несмотря на признанное отсутствие прямой корреляции между температурой и солнечной активностью, временная зависимость вейвлет-коэффициента кросс-корреляции двух температурных рядов на временном масштабе 11 лет воспроизводит долгосрочные вариации основного солнечного цикла, указывая на существование неявных связей вариаций планетарного масштаба с солнечной активностью.

Блок 2. "Динамика конвективных течений при наличии свободноплавающих включений"
По результатам численного моделирования динамики погруженного тела в конвективном слое Релея-Бенара проведен анализ динамики плавающего тела и статистических характеристик турбулентных полей скорости и температуры, полученных в рамках двухмерной и полной трехмерной численных моделей. В частности, показано, что тепловой поток в 3D случае выше, а разница в значениях числа Нуссельта для 2D и 3D моделирования увеличивается с ростом числа Рэлея. 3D-моделирование выявило асимметрию потока, связанную с крупномасштабной циркуляцией жидкости. Асимметрия более выражена в случае прямоугольной пластины, хотя асимметрия в гистограмме скоростей пластины наблюдается и для круглой пластины. Этот результат интересен тем, что в экспериментах наблюдали режимы с асимметричными движениями плавающего тела, которые объясняли наличием остаточного крупномасштабного вихря. В этом случае асимметрия хорошо выражена, поскольку плавающее тело колеблется у одной стенки, практически не достигая противоположной стенки. На отчетном этапе исследования обнаружен качественно иной режим, в котором асимметрия возникает в режиме «конвективного маятника». Асимметрия относительно слаба и проявляется только при детальном анализе динамики плавающего тела. Получены количественные зависимости интегральных параметров, характеризующих интенсивность конвективного течения жидкости в слое (число Рейнольдса), интегральный тепловой поток (число Нуссельта) и интенсивность перемещений плавающего в слое тела (средний период полного цикла) от степени радиационного нагрева (плотности падающего радиационного потока, который определяет устанавливающееся число Релея. Для пластины, частично пропускающей и частично отражающей падающее излучение, получены зависимости характеристик конвективного течения и движений пластины от коэффициента пропускания (отношение пропущенного светового потока к падающему).

Блок 3. "Конвективный тепломассоперенос в крупных городах и агломерациях"
В результате проведенных экспериментальных работ по исследованию обтекания набегающим потоком набора макетов морфотипов городской застройки получены разносторонние данные для настройки и апробации численных подходов. Установлено удовлетворительное согласование результатов расчетов и данных экспериментов. Результаты численного моделирования продемонстрировали, что разработанная методика и реализованный инструментарий визуального отображения в программе «SigmaEco» позволяют получить количественные данные по биоклиматическим критериям и отрисовать соответствующие полученным значениям зоны пешеходной комфортности.
Результаты численных моделирования ветрового режима на территории кампуса Пермского университета подтвердили, что реализуемые численные алгоритмы адекватно воспроизводят результаты прямых измерений ветрового потока. Показано, что геометрическая детализация объекта исследования необходима для качественного прогнозирования условий локального аэрационного режима территории.
Показано, что по эффективно-эквивалентной температуре климат в Перми и городах-миллионниках Сибири более суровый, чем по эффективной температуре, особенно в холодный период года. Это объясняется тем, что при температуре ниже +7°С любой ветер является охлаждающим фактором.
По данным метеорологических станций установлены особенности городского острова тепла в Красноярске. В условиях жаркой погоды городской остров тепла ночью выражен ярче, чем днем. Наряду с городским островом тепла, в Красноярске во время жаркой погоды прослеживается формирование острова сухости, когда относительная влажность уменьшается по сравнению с загородными территориями на 20 %.
Дальнейшие усилия по улучшению точности прогноза температуры в городской застройке могут быть направлены на подбор оптимальной параметризации пограничного слоя в стандартных настройках модели WRF-ARW. Путь решения задачи для ограниченной территории исследования лежит в области интеграции региональной модели WRF-ARWс микромасштабной моделью с шагом менее 10 м, в которой процесс обтекания отдельных зданий реализуется явно.

Публикации проекта вышедшие и принятые к печати в 2024 году
1. Сухановский А.Н., Гаврилов А.А., Попова Е.Н., Васильев А.Ю. The study of the impact of polar warming on global atmospheric circulation and mid-latitude baroclinic waves using a laboratory analog Weather and Climate Dynamics (2024) wos, scopus, rsci, РИНЦ, Q1
2. Быков А.В., Ветров А.Л., Калинин Н.А. Прогноз сильной жары в Красноярске с использованием региональной модели WRF-ARW Гидрометеорологические исследования и прогнозы (2024) rsci, РИНЦ, «Белый список»
3. Степанов Р.А., Соколов Д.Д., Фрик П.Г. Consistency of climatic changes at different time scales in Central England and Greenland Journal of Atmospheric and Solar-Terrestrial Physics (2024) wos, scopus, rsci, РИНЦ, «Белый список»
4. Калинин Н.А., Ветров А.Л. Особенности изменений индексов комфортности климата в Перми и городах-миллионниках Сибири в холодный и теплый периоды года Географический вестник (2024) rsci, РИНЦ, «Белый список»
5. Мешкова В.Д., Филимонов С.А., Дектерев А.А., Ветров А.Л., Шульженко П.Д. Numerical prediction of bioclimatic comfort level of urban development in street spaces E3S Web of Conferences (2024) scopus, РИНЦ
6. Филимонов С.А., Гаврилов А.А., Фрик П.Г., Сухановский А.Н., Васильев А.Ю. 2D and 3D numerical simulations of a convective flow with a free-floating immersed body Heat Transfer Research (2025) wos, scopus, РИНЦ, «Белый список»
7. Сухановский А.Н., Васильев А.Ю., Гаврилов А.А., Евграфова А.В., Попова Е.Н. Конвективный теплоперенос в крупномасштабных природных системах Инженерно-физический журнал (2025) scopus, rsci, РИНЦ, «Белый список»
8. Филимонов С.А., Гаврилов А.А., Литвинцев К.Ю., Васильев А.Ю., Сухановский А.Н., Фрик П.Г. Numerical modeling of immersed plate motion in a convective layer under radiation heating EURASIAN JOURNAL OF MATHEMATICAL AND COMPUTER APPLICATIONS (2025) scopus, rsci, РИНЦ, «Белый список»
9. Литвинцев К.Ю., Мешкова В.Д. The influence of solar and thermal radiation on the wind regime formation in urban block EURASIAN JOURNAL OF MATHEMATICAL AND COMPUTER APPLICATIONS (2025) scopus, rsci, РИНЦ, «Белый список»

Результаты проекта были широко представлены на российских и зарубежных конференциях в 2024 году

XVII Минский международный форум по тепломассообмену, Минск (Беларусь),20-24 мая 2024 г.:

1. Сухановский А.Н., Попова Е.Н., Васильев А.Ю. (Пермь), Гаврилов А.А. (Новосибирск). Конвективный теплоперенос в крупномасштабных природных системах (проблемный доклад)
2. Филимонов С.А., Гаврилов А.А. (Новосибирск), Фрик П.Г., Сухановский А.Н., Васильев А.Ю. (Пермь). Трехмерная модель движения погруженного диска в турбулентном свободно конвективном слое (устный доклад)

XX Международная конференция студентов, аспирантов и молодых ученых "Проспект Свободный - 2024" – 19 апреля 2024 – Красноярск:

3. Шульженко П.Д., Мешкова В.Д., Дектерев Д.А., Дектерев А.А. Исследование влияния строчного морфотипа городской застройки на формирование аэродинамического режима. (Диплом 2 степени)

Международная конференция по измерениям, моделированию и информационным системам для изучения окружающей среды (ENVIROMIS-2024), Томск, 1-6 июля 2024:

4. Сухановский А.Н., Степанов Р.А., Быков А.В., Ветров А.Л., Калинин Н.А., Фрик П.Г. «Бароклинные волны средних широт в случае зонально однородных граничных условий» (устный доклад)
5. Мешкова В.Д., Филимонов С.А., Дектерев А.А., Ветров А.Л., Литвинцев К.Ю. "Численные и экспериментальные исследования условий ветрового комфорта в городе Перми" (стендовый доклад)

Семинар в университете Канпура (Индия) :

6. Степанов Р.А. «Бароклинные волны средних широт в случае зонально однородных граничных условий» (устный доклад)

СТС-2024: «XL Сибирский теплофизический семинар», посвященный 110-летию со дня рождения академика С. С. Кутателадзе и 300-летию Российской академии наук:

7. Филимонов С. «Моделирование движения погруженной пластины в конвективном слое при радиационном нагреве дна» (устный доклад).
8. Шульженко П.Д., Мешкова В.Д., Дектерев А.А., Литвинцев К.Ю. "ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ФОРМИРОВАНИЯ УСЛОВИЙ ВЕТРОВОГО КОМФОРТА В РАЗЛИЧНЫХ ГОМОГЕННЫХ МОРФОТИПАХ ГОРОДСКОЙ ЗАСТРОЙКИ" (устный доклад)
9. Васильев А., «Режимы движения погруженного тела в конвективном слое» (устный доклад).
10. Мешкова В.Д., Дектерев Д.А., Шульженко П.Д., Дектерев А. А." ПРИМЕНЕНИЕ ПОВЕРХНОСТНЫХ МЕТОДОВ ВИЗУАЛИЗАЦИИ ДЛЯ АНАЛИЗА СТРУКТУРЫ ТЕЧЕНИЯ В ПРИЗЕМНОЙ ОБЛАСТИ ГОРОДСКОЙ ЗАСТРОЙКИ" (стенд).

IX Всероссийская научная конференция с элементами школы молодых ученых «Теплофизика и физическая гидродинамика: Современные вызовы», г. Сочи, 8-15 сентября 2024:

11.Мешкова В.Д., Филимонов С.А., Дектерев А.А., Ветров А.Л., Литвинцев К.Ю. «Анализ влияния детализации геометрии 3D-модели городской застройки на точность моделирования задач пешеходного ветрового комфорта» (устный доклад)
12. Литвинцев К.Ю., Мешкова В.Д. «Влияние солнечного и теплового излучения на формирование аэрационного режима в городском квартале.» (устный доклад)

Международный симпозиум «НЕРАВНОВЕСНЫЕ ПРОЦЕССЫ В СПЛОШНЫХ СРЕДА», Пермь, 1-4 октября 2024

13. Сухановский А.Н., Гаврилов А.А., Васильев А.Ю., Попова Е.Н. «Влияние арктического потепления на характеристики бароклинных волн средних широт и распределение температуры» - устный доклад
14. Ветров А.Л., Сухановский А.Н., Степанов Р.А., Быков А.В., Калинин Н.А., Фрик П.Г. «Характеристики бароклинных волн средних широт в случае зонально однородных граничных условий» - устный доклад
15. Беляева А.В., Сухановский А.Н. «Теплообмен в условиях город-ской среды на примере города Красноярска» - устный доклад

«Турбулентность, динамика атмосферы и климата», конференция, посвященная памяти академика Александра Михайловича Обухова,19–21 ноября 2024 г. Москва

16. Фрик П.Г., Сухановский А. Н., Степанов Р.А., Быков А.В., Ветров А.Л., Калинин Н.А. «Глобальная циркуляция и бароклинные волны средних широт в случае зонально однородных граничных условий» (Пленарный доклад)
17. Сухановский А.Н., Гаврилов А.А., Васильев А.Ю., Попова Е.Н. «Механизм арктического усиления в лабораторном аналоге общей циркуляции атмосферы» (Устный доклад)
18. Степанов Р.А., Фрик П.Г., Соколов Д.Д. «Согласованность климатических изменений в разных временных масштабах в отдаленных районах северной Атлантики» ( Устный доклад)

II Всероссийская научно-практическая конференция с международным участием "Цифровая география". Пермь, 25-28 сентября 2024 г.

19. Калинин Н.А., Старцева П.Ю. Оценка качества воспроизведения биоклиматических индексов моделью WRF-ARW (устный доклад).
https://xn--80aacgackdshf3de3ag2d2bxk.xn--p1ai/cifrovaya-geografiya

Russian Science Foundation Project No. 22-21-00572 " Development of laboratory model of general atmospheric circulation ". 

Project leader:

Sukhanovskii Andrei , senior researcher, D.Sc. ICMM UB RAS (Perm, Russia).

Project objectives:

The project aims to create and validate a laboratory model of the general circulation of the atmosphere.

 

At present, the most widely used model of general atmospheric circulation is a gap between two coaxial cylinders, on the walls of which different temperature values are set. This approach quite effectively models the large-scale circulation in the middle latitudes.  However, the general atmospheric circulation has a more complex structure and consists of several cells (Hadley, Ferrell, and polar), the seasonal intensity of which can vary appreciably.  It leads to the problem of the atmospheric circulation modeling  in a more realistic framework. The project participants, varying the layer thickness, the distance between heating and cooling sources, their intensity, and the bottom slope, plan to realize the meridional circulation in the form of three cells. The implemented laboratory model will serve as the basis for solving a whole set of problems linked to the baroclinic waves in atmosphere.

Supplementary materials:

video (in a rotating frame, accelerated 10 times)


0.48 rad/s Rот=0.2

0.37 rad/s Roт=0.4

0.13 rad/s Roт=3.1

0.08 rad/s Roт=8.4

 

Summary final report on the results of the second year of the RNF project No. 22-21-00572:

A new model of the general circulation of the atmosphere has been developed, which represents a rotating liquid layer with an aspect ratio noticeably less than unity, with a localized rim heater located at the periphery of the bottom and a disk-shaped cooler placed in the central part, at the upper boundary of the liquid layer. The rim heater simulates heating in the equatorial region and the cooler, cooling in the polar region. The heater is specially displaced from the side wall to minimize its influence on the formation of flows, due to non-slip condition at the boundary.

Based on the developed principal scheme, an experimental model was realized, which is a tank of square cross-section with a side of 700 mm and a height of 200 mm. To realize the cylindrical layer, an additional cylindrical wall made of Plexiglas with a thickness of 3 mm and radius R = 345 mm was inserted into the tank. At the bottom there is a rim, electric heater with a width of 25 mm. The distance from the side wall to the outer boundary of the heater is 40 mm. Cooling of the liquid is provided by heat exchange with ambient air on the free surface (room temperature is kept constant by the air conditioning system), a cooler located in the central part of the upper boundary, and heat flow through the side wall. Aluminum powder was used to visualize the flow pattern in the upper layer. Images were captured using a Bobcat 2020 4 megapixel CCD camera.

Experimental measurements provide only partial information about the system under study, therefore, to reconstruct the three-dimensional structure of the flow, a direct numerical simulation of thermal convection in a rotating cylindrical layer was performed using the freely distributable computational fluid dynamics package OpenFOAM v2106. The computational domain is a digital copy of the experimental model in terms of its geometric dimensions, location of the heater and cooler.

The solution of the main task of the project required a large number of experiments and numerical simulations. Heating power (from 8 to 120 W), rotation speed (from 0.08 to 0.48 rad/s), layer thickness (3 cm and 6 cm), physical properties of the liquid (water and PMS-5, Prandtl numbers 5.4 and 63) were varied. A large amount of data obtained by experiments and mathematical modeling was analyzed. On the basis of these results three papers were prepared and published, including in two high-level journals specializing in geophysical topics (Q2), one more paper has been prepared and is under review.

The following can be highlighted as the main results obtained in the frame of the project:

- the principal differences in the formation of mean large-scale flows and baroclinic waves in different atmospheric circulation models are shown. An important difference of the presented model from classical annulus configurations is the absence of steady waves. All wave modes, even with a regular wave structure, are characterized by strong non-periodic fluctuations. The observed baroclinic wave structures are a combination of different baroclinic modes evolving in time. This proves that the spatial distribution of heating and cooling, their location, and the type of boundary condition are important for the stability of baroclinic waves. With increasing rotational velocity, the waves in our model become irregular and represent a set of azimuthal wave modes, the main energy of which is contained in modes from m=2 to m=8, which agrees well with the data for the real atmosphere. A map of regimes in the plane of the thermal Rossby number-Taylor number is constructed.

- the influence of processes in the Ekman layer on the formation of mean flows and baroclinic waves was analyzed. The analysis was carried out on the basis of data from numerical calculations. It is shown that the main source of meridional circulation is thermal convection, and the radial flows induced in the viscous boundary layer are not of primary importance. The distribution of velocity fluctuations shows that, except for the heating region, they are mainly concentrated in the upper layer. The velocity fluctuations are damped in the Ekman layer, and as a result the baroclinic waves in our model are localized in the upper layer. It is consistent with the results of other authors, both for the laboratory system and for the real atmosphere.

- the most important result of the project is the successful realization in the laboratory model of the flows whose structure is qualitatively similar to the circulation in the real atmosphere. The model allows us to obtain the meridional circulation consisting of three cells, analogs of Hadley, Ferrel (of baroclinic wave nature) and polar cells. The structure of baroclinic waves, the location of their formation, and the mode composition are also similar to what is observed in the atmosphere. It is shown that the atmospheric mode in the laboratory model is realized only in a limited range of parameters. Decreasing the rotation rate leads to regularization of baroclinic waves, and increasing the rotation rate on the contrary leads to destruction of long waves and geostrophic turbulence. The realized model makes it possible to carry out long-term controlled experiments and can serve as an effective tool for studying various aspects of the global atmospheric circulation.

- to study the influence of the beta effect on the structure of mean flows and baroclinic waves, experiments and numerical simulations with a sloping bottom were realized. The simulations showed that a twofold change in the aspect ratio (from 0.09 to 0.17) does not lead to a qualitative change in the structure of the mean meridional circulation. A relatively small tilt (5 degrees) leads to noticeable quantitative changes in the fluid flow, preserving the overall flow structure. However, a further increase of the inclination up to 10 degrees leads to crucial changes both in the flow structure and in the intensity of wave motions. The formation of a downward cold flow due to the bottom slope (analog of slope currents) has a strong influence on the current structure. The presence of bottom slope leads to intensification of wave motions. Transformation of the meridional circulation at large tilts leads to noticeable changes of zonal flows, in particular, to the formation of an anticyclonic circulation near the bottom. The general changes in the structure and dynamics of flows lead to noticeable changes in the mean temperature distribution. Experiments were carried out for an inclination angle of 5 degrees, a fixed heating power (120 W), two layer thicknesses (35 mm and 60 mm), three rotational velocities of 0.11 rad/s (mode without baroclinic waves), 0.23 rad/s (mode of regular unsteady baroclinic waves), and 0.37 rad/s. Qualitative analysis of the observed flows showed that the beta effect has a significant influence on the structure and dynamics of the waves.

На портале Минобрнауки вышла статья, посвящённая работе учёных из ИМСС УрО РАН

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

Read more: Цифровые двойники помогут оценить рентабельность нефтяных месторождений

В ведущем научном журнале в области механики жидкости Journal of Fluid Mechanics вышла статья сотрудников ИМСС УрО РАН Мизёва А.И., Шмырова А.В. и Шмыровой А.И. «On the shear-driven surfactant layer instability». Работа посвящена изучению течений в жидкостях, на поверхности которых имеется слой молекул поверхностно-активного вещества (или сурфактанта). Межфазная гидродинамика, в рамках которой рассматриваются такого типа задачи, до сих пор далека от полного понимания явлений, несмотря давнюю историю развития. Это связано с междисциплинарным характером таких проблем, лежащих на стыке гидродинамики и физической химии. Слои сурфактанта толщиной всего в одну молекулу представляют собой двумерные объекты, характеризующиеся своей собственной диффузией и реологией, что необходимо учитывать при описании течений жидкости.

Read more: На стыке гидродинамики и физической химии

Subcategories

Здесь планируется размещать Веб-версии сборников статей