'

Проектирование и эффективная реализация средств численного моделирования в задачах вихревой гидроаэродинамики

Понравилась презентация – покажи это...





Слайд 0

Проектирование и эффективная реализация средств численного моделирования в задачах вихревой гидроаэродинамики С.В. Гувернюк, Д.А. Григоренко, А.И. Гирча, П.Р. Андронов, Г.Я. Дынникова НИИ механики МГУ


Слайд 1

Плоские нестационарные течения вязкой (или идеальной) несжимаемой жидкости Re = 2000 Связанная задача: совместное решение уравнений гидродинамики и динамики при движении тела в жидкости


Слайд 2

Используемые математические модели В качестве базовых моделей среды используется две: модель идеальной несжимаемой жидкости (двумерные уравнения Эйлера) модель вязкой несжимаемой жидкости (двумерные уравнения Навье-Стокса) Во всех случаях начальное состояние среды – покой во всем пространстве Движущиеся тела – двумерные, произвольной формы конечной толщины Граничные условия: условие прилипания (в вязкой жидкости) условие непротекания с априорным заданием точек отрыва (в идеальной жидкости)


Слайд 3

«Вихревая» формулировка плоской сопряжённой задачи динамики и аэрогидродинамики R В начальный момент среда и тело покоятся в неограниченном пространстве. Тело начинает двигаться под действием внешних сил , моментов и наложенных связей, испытывая также дополнительное влияние аэродинамических сил и моментов со стороны среды, возмущённой движением тела. - гидродинамический импульс среды (Г. Ламб, 1947) - вращательный импульс среды (Г. Ламб, 1947) - P.G.Saffman. Vortex Dynamics (1992) Для безвихревого движения теория развита Томсоном, Тэтом, Кирхгофом (1869).


Слайд 4

Лагранжевы методы Метод дискретных вихрей (МДВ) С.М. Белоцерковский и др. Модификации Кратные цепочки С.Я. Герценштейн, С.В. Гувернюк, Х. Исванд Метод вязких вихревых доменов (ВВД) П.Р. Андронов, С.В. Гувернюк, Г.Я. Дынникова Уравнения Навье-Стокса Генерация вихрей Диффузионная скорость Y. Ogami, T. Akamatsu


Слайд 5


Слайд 6

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


Слайд 7

Требования Необходимо обобщение Различные типы тел Различные типы взаимодействия Эффективность Переносимость Расширяемость Отчуждаемость


Слайд 8

Структура комплекса Средства подготовки заданий Скрипты Средства расчета Средства обработки результатов Скрипты Дополнительные средства (GUI) Сторонние программы


Слайд 9

Обобщенный алгоритм Единые обозначения величин Подробное описание, вспомогательные формулы Замкнутые двумерные тела произвольной формы Случаи идеальной и вязкой среды Число степеней свободы 0, 1, 2, 3 Возможны внешние силы Различные сценарии схода вихрей Несколько не сталкивающихся тел Множество рассчитываемых параметров (Cx,Cy,F,V,P,t…)


Слайд 10

Архитектура Счетная программа: ООП С++ Графический интерфейс: Python Qt Концепция: Open Source GPL


Слайд 11

Структура программы


Слайд 12

Эффективность Факторы: Прозрачность кода Переносимость Основная проблема: расчет поля скоростей Пути повышения: Оптимизация Профилирование Использование возможностей компилятора Распараллеливание Изменение алгоритма


Слайд 13

Алгоритмы Точные Приближенные (до ) Отличия от классической задачи n тел: Аннигиляция вихрей Нет симметрии


Слайд 14

Использование Т-системы Автоматическая балансировка загрузки Гетерогенность (сериализация) Новые возможности (в сравнении с MPI) Недостатки: Лишние пересылки данных Тонкие ограничения на C++ Модификатор static


Слайд 15

Время работы MPI


Слайд 16

T-система vs MPI Вычислительные машины: AMD Athlon(tm) 64 Processor 3000+ Intel(R) Celeron(R) CPU 2.53GHz Результаты T-система (Celeron+Athlon): 70 секунд MPI (Celeron+Athlon): 105 секунд Athlon: 73 секунды Celeron: 171 секунда Идеально: 51 секунда


Слайд 17

Переносимость Windows (MSVC, Cygwin + g++) Linux (g++) AIX (xlC) MPI T-система Трудности при переносе возможны


Слайд 18

Достоверность результатов Теоретическое обоснование Проверка инвариантов Сравнение с расчетами Библиотека GMP Другие методы Сравнение с экспериментом


Слайд 19

Нестационарное обтекание цилиндра квадратного сечения при Re = 57. Мгновенные линии тока. Слева – расчёт В.П. Шкадовой и др. (2004) сеточным методом при Re=56,6; справа – расчёт методом ВВД. Интервалы безразмерного времени между последовательными картинками слева и справа одинаковые.


Слайд 20

Результаты тестовых расчётов на основе метода ВВД при тех же начальных данных для случая Re=1000 при разной степени дискретизации по времени и по пространству: чёрная линия – при разбиении контура пластины на 164 отрезка и при шаге по времени 0,025, красная линия – при разбиении контура пластины на 328 отрезков и при шаге по времени 0,0125. Максимальное наблюдавшееся отличие угловой скорости при разной степени дискретизации составляет около 5%.


Слайд 21

Сравнение расчётной и экспериментальной картин обтекания колеблющегося крылового профиля. Re=44 000 Чёрные точки – расчётная визуализация положения частиц пассивной примеси. Красные точки – положительные вихри, синие точки – отрицательные вихри. Поток направлен слева. Светлые линии – дымовая визуализация мгновенного поля течения в эксперименте. Степень дискретизации профиля в расчёте - 247 отрезков.


Слайд 22

Расчётная (слева) и экспериментальная (справа) картины течения при авторотации вертушки с двумя лопастями для различных фаз поворота (поток слева направо, вращение против часовой стрелки) Sh=0,064 Sh=0,089 Re = 128 000 НИИ механики МГУ, аэродинамическая труба А-8; 2003-2005 гг.


Слайд 23

Сравнение с экспериментом (3 лопасти) Кратные цепочки вихрей Автоколебания Скорость потока U=30 м/с, длина лопасти h= 1,4 R


Слайд 24

Сравнение с экспериментом (3 лопасти) Кратные цепочки вихрей Авторотация Скорость потока U=30 м/с, длина лопасти h= 1,4 R


Слайд 25

Практическое применение Расчет мощности ветроэнергетических установок Ветровые нагрузки на высокие строения Нестационарные аэродинамические нагрузки на крыло Нестационарная теплопередача при обтекании нагретых цилиндров


Слайд 26


Слайд 27


Слайд 28


Слайд 29

Расчёт методом ВВД эволюции поля завихренности при авторотации пластины для случая: Re = 1000; безразмерный момент инерции J = 20; начальный угол поворота ? = 10 градусов. Шкала от синего до красного цвета соответствует изменению завихренности от –0,5 до +0,5. а) б) в) г) д) е)


Слайд 30


Слайд 31


Слайд 32

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


Слайд 33

Спасибо за внимание Авторы выражают благодарность Васенину В.А., Водомерову А.Н.


×

HTML:





Ссылка: