![]() |
ИСТИНА |
Войти в систему Регистрация |
ИСТИНА ЦЭМИ РАН |
||
Целью данного проекта является создание явного и легко масштабируемого балансно-характеристического метода высокой разрешающей способности для решения различных задач индустриальной математики, требующих постановки в смешанных эйлерово-лагранжевых переменных. В частности, метод планируется использовать для решения задач взаимодействия потока жидкости или газа с деформируемыми объектами (FSI – fluid-structure interaction). Примером одной из таких задач является задача об обтекании частично закрепленных в подвижном канале деформируемых цилиндров, являющаяся важной модельной задачей в атомной энергетике (ОКБМ Африкантов). В рамках проекта планируется выполнение следующих задач: 1) Разработка явного балансно-характеристического метода решения систем квазилинейных уравнений в частных производных гиперболического типа в лагранжевых переменных на блочно-структурированных сетках. 2) Верификация разработанного метода для лагранжевых переменных отдельно на системах уравнений газодинамики и динамической упругости, решение модельных и тестовых задач. 3) Разработка “монолитного” балансно-характеристического алгоритма решения систем квазилинейных уравнений в частных производных гиперболического типа в смешанных эйлерово-лагранжевых переменных на блочно-структурированных сетках путем “бесшовного” соединения схемы КАБАРЕ для эйлеровых переменных и разработанного ранее метода для лагранжевых переменных. 4) Верификация “монолитного” метода для эйлерово-лагранжевых переменных на модельных задачах о взаимодействии потока с деформируемым телом. 5) Масштабирование разработанных алгоритмов на системы с общей и распределенной памятью. 6) Валидация разработанного “монолитного” метода на задаче об обтекании частично закрепленных в подвижном канале деформирующихся цилиндров.
The project result will be the creation of a new explicit conservative-characteristic method for solving equations of hyperbolic type in mixed Eulerian-Lagrangian variables in two-dimensional and three-dimensional cases. The constructed method will be adapted to problems of the interaction of deformable bodies with a gas or liquid flow (FSI - fluid-structure interaction). The method will be the representative of “monolithic” methods for FSI problems, in which the equations of dynamic elasticity and gas dynamics will be solved by a single “seamless” algorithm that processes the boundary between deformable bodies and flow using Riemannian invariants that arrive at the boundary by characteristics from both sides. The resulting “monolithic” high-resolution algorithm will possess the properties of conservation and time reversibility, unique to modern methods of solving FSI problems. Also, unlike the overwhelming majority of FSI numerical methods, the conservative-characteristic algorithm will be explicit with a minimal computational template and, as a result, easily and efficiently scalable. Based on the constructed method, a parallel computer program will be developed. Using it, the method will be tested on model problems in two-dimensional and three-dimensional cases. Using the Lomonosov Moscow State University supercomputers the flow around deformable cylinders partially fixed in a moving channel will be modelled. This problem is an important model problem in nuclear energy. The results will be compared with experimental data. The developed program can subsequently be used to calculate large industrial problems.
Будет построен новый явный балансно-характеристический метод для решения уравнений гиперболического типа в смешанных эйлерово-лагранжевых переменных в двухмерном и трехмерном случаях на четырехугольных и гексаэдральных сетках соответственно. Построенный метод будет адаптирован для решения задач о взаимодействии деформируемых тел с потоком газа или жидкости (FSI). Полученный “монолитный” алгоритм высокой разрешающей способности будет обладать уникальными для современных методов решения задач FSI свойствами консервативности и временной обратимости. Также в отличие от подавляющего числа численных методов FSI балансно-характеристический алгоритм будет явным с минимальных вычислительным шаблоном и, как следствие, легко и эффективно масштабируемым. На основе построенного метода будет разработана параллельная компьютерная программа. С помощью нее метод будет протестирован на модельных задачах в двумерном и трехмерном случаях. С помощью программы на суперкомпьютерах МГУ им. М.В. Ломоносова будет рассчитана задача об обтекании частично закрепленных в подвижном канале деформируемых цилиндров, являющаяся важной модельной задачей в атомной энергетике. Полученные результаты будут сравнены с экспериментальными данными. Разработанную в рамках проекта программу впоследствии можно будет использовать для расчетов крупных индустриальных задач.
Балансно-характеристические схемы и схема КАБАРЕ в частности активно развиваются Головизниным В.М. и соавторами в течение последних двух десятилетий. С помощью них были решены многие задачи океанологии, турбулентной газодинамики, конвективного теплопереноса, атомной энергетики, аэроакустики и множество других серьезных задач индустриальной математики. В 2013 году в издательстве МГУ имени М.В. Ломоносова вышла монография, подводящая итог пятнадцатилетней работе по этому направлению: Головизнин В.М., Зайцев М.А., Карабасов С.А., Короткин И.А. Новые алгоритмы вычислительной гидродинамики для многопроцессорных вычислительных систем. По своим характеристикам новые алгоритмы, разрабатываемые научным коллективом во главе с Головизниным В.М., превосходят зарубежные аналоги, и могут быть названы CFD-алгоритмами нового поколения (Goloviznin V. M., Chetverushkin B. N. New generation algorithms for computational fluid dynamics // COMPUTATIONS MATHEMATICS AND MATHEMATICAL PHYSICS. — 2018.). Афанасьев Н.А. последние 4 года работал под руководством Головизнина В.М. над развитием балансно-характеристических схем в применении к задачам индустриальной математики. За первый год аспирантуры Афанасьевым Н.А. была успешно решена проблема обработки звуковых точек в балансно-характеристических схемах, что позволило моделировать газодинамические трансзвуковые течения с более высокой точностью. Также аспирант имеет опыт в решении крупных актуальных задач индустриальной математики, а именно задач трехмерного конвективного перемешивания и пристеночного теплообмена (в рамках НИР “Математическое моделирование отдельных физических процессов в ЗО, определяющих формирование локальных областей с высокой концентрацией водорода”) и исследования термоакустической неустойчивости камер сгорания газовых турбин (в рамках НИР “Разработка математических моделей и компьютерных программ для анализа термоакустической неустойчивости камер сгорания газовых турбин”).
МГУ имени М.В.Ломоносова | Координатор |
грант РФФИ |
# | Сроки | Название |
1 | 1 октября 2020 г.-31 декабря 2020 г. | Консервативные разностные схемы нового поколения для задач взаимодействия деформируемых тел с потоками жидкости и газа |
Результаты этапа: Написан обзор о современных методах решения задач о взаимодействии деформируемых тел с потоком жидкости или газа (FSI – fluid structrure interaction). Написан обзор о современных разностных методах решения систем квазилинейных уравнений гиперболического типа в лагранжевых переменных. Разработан явный балансно-характеристический метод решения двумерных уравнений газовой динамики в лагранжевых переменных в слабосжимаемом случае. | ||
2 | 1 января 2021 г.-31 декабря 2021 г. | Консервативные разностные схемы нового поколения для задач взаимодействия деформируемых тел с потоками жидкости и газа |
Результаты этапа: Разработаны явные балансно-характеристические методы решения двумерных уравнений газовой динамики в смешанных эйлерово-лагранжевых переменных и уравнений динамической упругости в лагранжевых переменных. Оба метода построены на основе идеи сохранения свойства временной обратимости, выполняемого для дифференциальных уравнений гиперболического типа. В отличие от разработанных ранее балансно-характеристических методов на подвижных сетках, указанные методы полностью обратимы по времени (при отключенных процедурах монотонизации) засчет обратимого по времени алгоритма переноса узлов сетки. Такой алгоритм удалось разработать путем интерполяции скоростей с ребер сетки на узлы перед каждым передвижением сетки. Разработанные методы для уравнений газовой динамики и динамической упругости объединены в "бесшовный" алгоритм расчета задач FSI. Получившийся алгоритм является монолитным (оба физических процесса моделируются по схожим схемам, явно учитывающим передвижение границы раздела), явным, легко масштабируемым, а также обладает минимальным вычислительным шаблоном в одну пространственно-временную ячейку и свойством временной обратимости. Других явных монолитных методов решения задач FSI с аналогичными свойствами не существует. Бесшовность алгоритма гарантирована характеристической фазой схемы, которая позволяет относительно просто обрабатывать границу раздела с помощью значений инвариантов Римана, переносимых на граничное ребро из газодинамической и упругой ячеек по обе стороны от ребра. Все описанные выше разработанные методы верифицированы на модельных и тестовых задачах, а также масштабированы на системы с общей (OpenMP) и распределенной (MPI) памятью. Алгоритм для решения трехмерных уравнений газовой динамики и динамической упругости на гексагональных лагранжевых сетках построен, но еще не верифицирован. | ||
3 | 1 января 2022 г.-30 сентября 2022 г. | Консервативные разностные схемы нового поколения для задач взаимодействия деформируемых тел с потоками жидкости и газа |
Результаты этапа: Разработаны явные балансно-характеристические методы решения двумерных и трехмерных уравнений газовой динамики в смешанных эйлерово-лагранжевых переменных и уравнений динамической упругости в лагранжевых переменных. Разработан новый явный алгоритм передвижения свободной границы и границы между упругим телом и газом с помощью аппроксимации уравнения движения лагранжевой границы балансно-характеристическим методом КАБАРЕ. Разработанные методы объединены в "бесшовный" явный балансно-характеристический алгоритм расчета задач взаимодействия потока с деформируемыми объектами (FSI). Произведена верификация указанных методов на модельных и тестовых задачах. |
Для прикрепления результата сначала выберете тип результата (статьи, книги, ...). После чего введите несколько символов в поле поиска прикрепляемого результата, затем выберете один из предложенных и нажмите кнопку "Добавить".