Программа 3 информационные и вычислительные технологии в задачах поддержки принятия решений отчет
Вид материала | Программа |
- Информационные технологии анализа данных используются для обобщения учетной информации,, 178.99kb.
- Системы поддержки принятия решений общие сведения о системах поддержки принятия решений, 145.37kb.
- Рабочая программа дисциплины «Системы поддержки принятия решений» Рекомендуется для, 115.78kb.
- Системы поддержки принятия решений: эволюция концепции и некоторые перспективы, 98.37kb.
- Анализ принятия управленческих решений, 54.28kb.
- Темы курсовых проектов по дисциплине «Теория принятия решений», 35.1kb.
- Программа по дисциплине "Интеллектуальные информационные системы" для подготовки специалиста, 151.52kb.
- Темы курсовых проектов по дисциплине «Теории принятия решений» Можаева, 30.34kb.
- Ложения темпоральной (временной) логики для ветвящегося времени в плане ее использования, 138.35kb.
- Рабочая программа учебной дисциплины теория принятия решений Направление подготовки, 591.05kb.
РОССИЙСКАЯ АКАДЕМИЯ НАУК
СИБИРСКОЕ ОТДЕЛЕНИЕ
Приоритетное направление 4.3.
Системы автоматизации, GALS– технологии, математические модели и методы исследования сложных управляющих систем и процессов
Программа 4.3.1.
ИНФОРМАЦИОННЫЕ И ВЫЧИСЛИТЕЛЬНЫЕ ТЕХНОЛОГИИ В ЗАДАЧАХ ПОДДЕРЖКИ ПРИНЯТИЯ РЕШЕНИЙ
ОТЧЕТ
о работе программы в 2009 г.
Институты-исполнители:
Институт автоматики и электрометрии СО РАН
Институт вычислительного моделирования СО РАН
Институт вычислительных технологий СО РАН
Институт систем энергетики им. Л.А. Мелентьева СО РАН
Координатор программы:
академик Ю.И.Шокин
Ученый секретарь программы:
д.ф.-м.н. Л.Б.Чубаров
Новосибирск, 2009
Отчеты по проектам 3
Публикации по программе 27
Отчеты по проектам
Проект 4.3.1.1. Информационно-вычислительные технологии в задачах поддержки принятия решений (Институт вычислительных технологий СО РАН)
На этапе 2009 г. предполагалось:
- Выполнить трехмерное математическое моделирование эволюции электромагнитных полей в метаматериалах.
- Разработать нестационарные модели течений в гидротурбинах для переходных режимов, учитывающие гидравлические удары.
- Провести расчеты в аксиально-симметричной постановке ВЧ - разряда во фторсордержащих газах на базе гидродинамической модели.
- Провести математическое моделирование спутного турбулентного течения с варьируемым значением суммарного избыточного импульса в линейно стратифицированной среде с целью изучения тонкой микроструктуры гидрофизических полей.
- Численно исследовать процесс охлаждения высокотемпературных продуктов сгорания твердотопливных газогенераторов с помощью насыпных фильтров в системах быстрого наддува эластичных оболочек.
- Провести численное моделирование сверхзвуковых течений в воздухозаборниках летательных аппаратов с возникновением режимов выбитой ударной волны, приводящей к запиранию потока.
Основные (важнейшие) результаты
Разработка эффективных параллельных вычислительных алгоритмов для моделирования задач нанофотоники
Исполнители
Федорук М.П., д.ф.-м.н., зам. дир.;
Лебедев А.С., к.ф.-м.н., с.н.с.;
Штырина О.В., к.ф.-м.н., н.с.;
Чубаров Д.Л., м.н.с.;
Прокопьева Л.Ю., м.н.с.;
Витковский В.Э., асп.
Разработан алгоритм метода конечных объемов для решения трехмерных нестационарных уравнений Максвелла в средах с переменными значениями диэлектрической и магнитной проницаемостей, позволяющий эффективно моделировать распространение электромагнитного излучения в сложных наноструктурированных материалах. Реализована параллельная версия этого алгоритма. Расчеты на высокопроизводительных вычислительных многопроцессорных комплексах продемонстрировали практически линейное ускорение для двумерного и одномерного случаев (, ).

Рис. 1. Ускорение параллельной версии конечно-объемного алгоритма в одномерном случае,


Рис. 2. Ускорение параллельной версии конечно-объемного алгоритма в двухмерном случаях,

Эффективность предложенного алгоритма продемонстрирована на примере расчета цилиндрической гиперлинзы (Рис. 3 – Рис. 6).
![]() | ![]() |
Рис. 3. Декомпозицая расчетной области для проведения параллельных расчетов | Рис. 4. Источники в вакууме, амплитуда поля ![]() |
![]() | ![]() |
Рис. 5. Внутренняя линза, амплитуда магнитного поля ![]() | Рис. 6. Внешняя линза, амплитуда магнитного поля ![]() |
Разработан вычислительный алгоритм и выполнена его параллельная реализация для моделирования процесса записи объемных микро - и наноструктур путем воздействия фемтосекундных лазерных импульсов на оптические световоды. Для моделирования процесса записи используется система нелинейных уравнений, которая описывает эволюцию медленной огибающей электрического поля


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

![]() | ![]() |
Рис. 7. Ускорение параллельного алгоритма относительно для кластера (Itanium 2) и универсального графического ускорителя (Tesla C1060) | Рис. 8. Распределение плотности электронов плазмы при ![]() |
Разработка нестационарных моделей течения в гидротурбинах для переходных режимов, учитывающих гидравлические удары.
Исполнители раздела
- Лапин В Н., к.ф.-м.н., н.с.;
- Черных Г.Г., д.ф.-м.н., г.н.с.;
- Чирков Д В., к.ф.-м.н., н.с.;
- Банников Д.В., асп.;
- Есипов Д В., асп.;
- Ешкунова И Ф., асп.;
- Авдюшенко А Ю., асп.;
- Астракова А С., маг.;
- Панов Л В., маг.;
- Жамбалова Д.Б., маг.
Разработаны математическая модель и методика расчета переходных процессов в гидравлических турбинах, возникающих при переходе из одного режима работы в другой, учитывающая явление гидравлического удара. Модель основана на совместном решении осредненных по Рейнольдсу нестационарных 3D уравнений Навье-Стокса в проточном тракте гидротурбины, уравнения вращения рабочего колеса как твердого целого и 1D уравнений «упругого» гидравлического удара в напорном водоводе. В большинстве случаев переходные режимы связаны с открытием/закрытием направляющего аппарата. Поэтому для описания изменения геометрии проточной части со временем базовый алгоритм расчета течения жидкости адаптирован к использованию подвижной сетки. Разработаны новые граничные условия, позволяющие задавать на входном сечении не расход жидкости, а полную энергию потока, что позволило решить широкий класс задач, как стационарных, так и нестационарных, где расход заранее неизвестен.
Выполнено моделирование одного из переходных процессов – начального этапа выхода в разгонный режим работы при фиксированном положении лопаток направляющего аппарата. При мгновенном снятии нагрузки с генератора частота вращения рабочего колеса (РК) начинает быстро увеличиваться, подчиняясь закону вращения РК как твердого целого.
Влияние структуры высокочастотного (ВЧ) разряда на процесс плазмохимического травления кремния в CF4/O2.
Исполнители
- Ю.Н. Григорьев, д.ф.-м.н., г.н.с.;
- А.Г. Горобчук, к.ф.-м., н.с.
На основе неизотермической модели плазмохимического реактора травления исследовано влияние неоднородной структуры ВЧ-разряда на основные характеристики процесса плазменного травления. Внутренние характеристики плазмы ВЧ-разряда рассчитывались в диффузионно-дрейфовом приближении для фторсодержащих газовых смесей. Выполнено численное моделирование радиального плазмохимического реактора с неоднородным ВЧ-разрядом и многокомпонентной плазменной кинетикой. Определено влияние структуры ВЧ-разряда на массообмен, скорость и однородность травления кремниевых образцов в зависимости от концентрации кислорода. Показано, что однородность травления образцов существенно зависит от изменения электронной плотности в радиальном направлении, которым часто пренебрегают, рассчитывая разряд в одномерной постановке. На рисунке (Рис. 9) показано распределение концентрации фтора CF (Моль/см3) в радиальном плазмохимическом реакторе.

Рис. 9. Распределение концентрации фтора CF (Моль/см3) в радиальном плазмохимическом реакторе: а - однородное распределение электронной плотности, б - распределение электронной плотности, имеющее радиальную зависимость в виде функции Бесселя. Параметры: p = 0.5 тор, Q = 200 см3/мин, Tw = 300 K, 30% содержание O2 в CF4/O2. Направление подачи газа - к центру реактора.
Математическое моделирование спутного турбулентного течения в устойчиво стратифицированной среде.
Исполнители:
- Г. Г. Черных, д.ф.-м.н,. г.н.с.;
- О.Ф. Воропаева, д.ф.-м.н, с.н.с.
Построены усовершенствованные численные модели динамики турбулентного следа в устойчиво стратифицированной среде. Для детального описания тонкой структуры течения в дальнем следе с нулевым избыточным импульсом построена численная модель, основанная на дифференциальных уравнениях переноса тройных корреляций турбулентных флуктуаций поля скорости, записанных с учетом вклада кумулянтов четвертого порядка, и на модифицированных алгебраических аппроксимациях тройных совместных корреляций турбулентных флуктуаций полей скорости и плотности. Результаты расчетов хорошо согласуются с известными экспериментальными данными об анизотропном вырождении следа.
Для описания динамики турбулентного следа за буксируемым телом в линейно стратифицированной среде построена иерархия математических моделей турбулентности второго порядка. Наиболее сложная математическая модель включает в себя наряду с дифференциальными уравнениями переноса компонент тензора рейнольдсовых напряжений уравнение переноса тройной корреляции турбулентных флуктуаций вертикальной компоненты скорости. Показано, что, как и в случае однородной жидкости, суммарный избыточный импульс порядка 5-10 % от суммарного избыточного импульса за буксируемым телом слабо влияет на вырождение энергии турбулентности. Основное влияние малого суммарного избыточного импульса сказывается на вырождении дефекта осредненной продольной компоненты скорости. Турбулентные следы с малым суммарным избыточным импульсом генерируют внутренние волны, слабо отличающиеся от внутренних волн, генерируемых безымпульсным турбулентным следом.
Моделирование работы насыпного фильтра для охлаждения высокотемпературных продуктов сгорания твердотопливного газогенератора.
Исполнители:
- А.Д. Рычков, д.т.н., г н.с.
Разработана математическая модель и реализующий ее численный алгоритм для моделирования процесса охлаждения высокотемпературных продуктов сгорания твердотопливных газогенераторов в фильтрах активного охлаждения. Результаты численного моделирования показали более высокую эффективность охлаждения продуктов сгорания газогенератора при использовании активных фильтров по сравнению с пассивными.
![]() Рис. 10. Изменение давления в камере охлаждения при x = 0 для активного (1) и пассивного (2) фильтров | ![]() Рис. 11. Изменение температуры на выходе из камеры охлаждения для активного (1) и пассивного (2) фильтров |
Численное моделирование сверхзвуковых течений в воздухозаборниках летательных аппаратов с возникновением режимов выбитой ударной волны, приводящей к запиранию потока.
Исполнители:
- В.М. Ковеня, д.ф.-м.н., г.н.с;.
- А.В.Базовкин, аспирант;
- А.Ю. Слюняев, аспирант.
Проведена модификация численного алгоритма решения уравнений Навье-Стокса сжимаемого теплопроводного газа, позволившая повысить точность расчетов за счет применения неравномерных адаптивных сеток, сгущающихся в областях больших градиентов. Проведены расчеты сверхзвукового обтекания модельной конфигурации гиперзвукового летательного аппарата (ГЛА) в продольной плоскости сечения в двумерном приближении в широком диапазоне чисел Маха, Рейнольдса и углов атаки. Получены основные закономерности течений около несущей поверхности аппарата, в воздухозаборнике и в окрестности тела за ним при различных краевых условиях для температуры поверхности стенок. Численно подтверждено явление выбитой ударной волны – отхода ударной волны от носка воздухозаборника с увеличением угла атаки, имеющее место в режимах течений, близких к запиранию канала. При углах атаки до 9 градусов ударная волна присоединена к кромке, а при больших углах атаки отходит от тела. Иллюстрацией этого явления служит рисунок (Рис. 12), где приведены поля скоростей при



Рис. 12. Распределения полей скорости газа в области входа в канал воздухозаборника