автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.16, диссертация на тему:Пакет прикладных программ GAS DYNAMICS TOOL и его применение в задачах численного моделирования газодинамических процессов

Введение.

Глава 1. Цель и задачи работы, анализ состояния вопроса.

1.1 Особенности исследования нелинейных процессов. газовой динамики. Актуальность работы.

1.2 Цель работы и задачи исследования.

1.3 Выбор базового метода численного интегрирования.

Р> Анализ состояния вопроса.

Глава 2. Пакет прикладных программ GasDynamicsTool

2.1. Базовая физико-математическая модель.

2.2. Метод решения и разностная схема.

2.3. Общие функциональные возможности пакета.

2.4. Система ввода и вывода информации.

2.5. Быстродействие и работа с памятью.

2.6. Выводы и результаты.

Глава 3. Тестирование кода

3.1. Общие замечания.

3.2. Взаимодействие ударных волн с поверхностями.

3.3 Нестационарное истечение сверхзвуковой струи в затопленное пространство.

3.4. Распределение газодинамических параметров при взаимодействии потоков с поверхностями.

3.5. Распад разрыва.

3.6. Выводы и результаты.

Глава 4. Численное моделирование многомерных газодинамических ^ процессов при функционировании ракетных и артиллерийских систем.

4.1. Общие замечания.

4.2. Моделирование процессов в модельном дульном тормозе.

4.3. Период последействия пороховых газов на снаряд при выстреле.

4.4. Дульный глушитель.

4.5. Газодинамика старта реактивных систем.

4.6. Моделирование распространения струйных течений в помещениях ф» ограниченного объема.

4.7. Стационарные газодинамические процессы для элементов ракет и снарядов.

4.8. Выводы и результаты.

Глава 5. Моделирование детонационных и ударно - волновых процессов при взрывах гомогенных составов конденсированных взрывчатых веществ, низкоплотных взрывчатых составов и топливо - воздушных смесей.

5.1. Фугасное поле и его воздействие на различные объекты.

0 5.2. Физическая и математическая модели процессов. детонации ТВС и КВВ в воздухе.

5.2.1. Основные принципы построения физических и математических моделей для детонации ТВС и КВВ.

5.2.2. Особенности физических моделей процесса детонации КВВ и ТВС.

5.2.3. Математическая постановка задачи, система уравнений и основные допущения.

5.3. Базовое тестирование и результаты численных экспериментов.

5.3.1. Взрывы ТВС в замкнутых объемах.

5.3.2. Детонационный взрыв КВВ в замкнутых объемах. 5.3.4. Сравнение с решением автомодельных задач.

5.3.5. Взрыв заряда КВВ около жестких поверхностей.

5.3.6. Взрыв ТВ С в типовой квартире.

5.3.7. Взрыв противотанковой мины.

5.4 Выводы и результаты.

Глава 6. Анализ процессов дефлаграции и детонации многокомпонентных гомогенных и гетерогенных систем.

6.1. Анализ механизма и методов решения задач дефлаграции.

6.1.1. Основные представления о горении. многокомпонентных систем.

6.1.2. Режимы химических реакций в многокомпонентных смесях.

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

6.1.4. Анализ методов решения задачи дефлаграции многокомпонентных систем.

6.2. Разработка и предварительное тестирование солверов и системы представления данных для задач дефлаграции.

6.2.1. Математическая формулировка задачи, обоснование метода решения и соответствующей разностной схемы.

6.2.2. Результаты предварительных тестовых расчетов, сравнение с известными данными.

6.2.3. Проведение тестов по сравнению с аналитическими решениями.

6.3. Анализ возможности многомерного моделирования физических процессов в многофазных средах. 6.3.1. Математические модели для описания течения гетерогенных многокомпонентных сред.

6.3.2. Приближенная модель горения алюминиевых частиц.

6.4. Физическая и математическая формулировка односкоростной многофазной модели многокомпонентных энергоносителей, содержащих твердую фазу (алюминиевые частицы).

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

6.4.2. Правые части уравнений.

6.4.3. Кинетика воспламенения и горения алюминия.

6.4.4. Режимы химических реакций продуктов детонации углеводородных топлив в присутствии алюминиевых частиц.

6.4.5. Выбор численных коэффициентов для определяющих с соотношений модели горения алюминиевых частиц.

6.5. Анализ уравнений состояния продуктов взрыва высокоплотных составов СаНаОсН*

6.5.1 Общие положения.

6.5.2 Обобщенное уравнение состояния для расчета температуры продуктов детонации мультифазных составов произвольной плотности.

6.5.3. Модель детонации многофазных систем произвольной плотности.

6.5.4 Модификация алгоритма расчета параметров детонационной волны в топливовоздушной смеси (ТВС).

6.6. Тестирование компьютерной версии односкоростной многофазной модели.

6.6.1. Выбор и обоснование тестовых задач.

6.6.3. Анализ результатов решения тестовых задач.

6.7. Выводы и результаты.

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

Эти заключительные слова из книги «Численное решение многомерных задач газовой динамики» под редакцией С.К.Годунова [118] написаны более четверти века назад. По сравнению со временем развития численных методов интервал значительный. В то время еще не было персональных компьютеров и программы носили не на компакт дисках а в портфеле в пачках перфокарт. А писали эти программы не на С++, а в машинных кодах.

Но несмотря на столь разные условия в создании средств интегрирования нелинейных систем газовой динамики тогда и сейчас, эти слова не потеряли актуальность.

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

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

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

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

Данная диссертационная работа посвящена разработке универсального пакета прикладных программ GasDynamicsTool и завершает 12-летний период по его созданию.

Проект разрабатывается с 1989 года. Изначально он был предназначен для моделирования нелинейных процессов при истечении газов в затопленное пространство, при наличии в области течения тел со сложной (произвольной) конфигурацией. Пакет был реализован на языке FORTRAN и использовался на машинах типа ЕС-1060. В последствии, проект значительно расширился и появилась реальная возможность его более интенсивного применения в различных приложениях. К 1992г. была создана первая версия программы (версия 2.36, [54]) предназначенная для использования на персональном компьютере с операционной средой DOS. Эта программа требовала минимальной подготовки пользователя и была способна решать широкий круг задач для нестационарных течений системы невязких сжимаемых газов на основе интегрирования уравнений Эйлера методом крупных частиц [8, 9], который является развитием метода частиц в ячейках Ф. Харлоу. Максимальная сетка для программы составляла 16000 ячеек. Данное ограничение было связано с организацией работы с памятью в операционной системе DOS и аппаратной базой компьютеров того времени. Программа GDT 2.36 имела хорошо развитый графический интерфейс пользователя, написанный на языке MODULA 2, некоторые характерные черты которого сохранились в последующих версиях. Версия 2.36 была завершена и протестирована в декабре 1992 г.

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

В 1996г. была создана первая 32-битная версия GDT (версия 3.0, язык BORLAND С++, компилятор 5.01, [66, 67]), в которой были сняты ограничения на количество расчетных ячеек, реализован алгоритм решения уравнений Навье-Стокса [63] с учетом теплопроводности, и диффузии, в схему введен механизм учета химических реакций в газах, что открыло возможность изучения таких практически важных явлений как горение и детонация, введена возможность учета воздействия гравитационных или массовых сил. В феврале 2000 г. была завершена разработка 4-й версии пакета [56]. Она включила серию принципиально новых разработок. Была реализована система трехмерной визуализации на основе применения технологии воксельной графики, разработана развитая технология импорта геометрических объектов из стандартных графических форматов (STEP), разработан солвер для решения многофазных задач в односкоростном приближении [55, 57, 59, 62], существенно улучшен блок граничных условий и сделан ряд других существенных улучшений кода и интерфейса.

В настоящее время программа GDT реализована для нескольких операционных систем - DOS, Windows 95 / 98 / 2000 и Windows NT, и имеет возможность рассчитывать области в десятки миллионов ячеек при весьма умеренных требованиях к характеристикам компьютеров, что очень важно для использования на отечественных предприятиях, как правило, имеющих в качестве компьютерной базы персональные машины.

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

ЦЕЛЬ И ЗАДАЧИ РАБОТЫ, АНАЛИЗ СОСТОЯНИЯ ВОПРОСА.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ ПО РАБОТЕ

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

2. Проанализирован класс околодульных газодинамических процессов при наличии в области течения подвижных и неподвижных тел со сложной трехмерной геометрией.

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

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

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