Математика

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

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

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

Кратко описаны исходные уравнения и разностная схема для моделирования вязких и невязких газодинамических течений многокомпонентной среды в эйлеровых переменных в коде ЭГАК. Приведена теоретическая оценка схемной вязкости разностной схемы методики. Проведены тестовые расчеты, которые показали корректность использования теоретической оценки схемной вязкости в расчетах конкретных задач

Приведены результаты прямого численного мод6елировани6я пристенной турбулентности в плоском канале на расчетных сетках от 4!10 до 34!10 ячеек для значений числа Рейнольдса 5600, 13750, 21900. Расчеты проводились по модификации схемы КАБА- РЕ. На основе анализа корреляционных зависимостей и спектральных характеристик потока обоснован выбор геометрических размеров расчетной области и расчетных сеток.

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

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

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

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

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

Дается описание алгоритма, реализующего метод сшивания многогранных неструктурированных сеток. Данный метод реализован в программе расчета начальных данных методики ТИМ-3D. Применимость метода продемонстрирована на примерах сшивания различных типов сеток в единую счетную сетку области. Проведено исследование быстродействия алгоритма, выполняющего сшивание трехмерных сеток

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

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

Сервисная система GeomGrid2 предназначена для подготовки, формирования и расчета начальных данных двумерных задач механики сплошной среды. GeomGrid2 является подсистемой программного комплекса БАЗИС. Подсистема GeomGrid2 позволяет в визуальном интерактивном режиме формировать начальную геометрию задачи, строить расчетные сетки (структурированные и неструктурированные) и поля точечных масс (частиц), задавать физические свойства для расчетов двумерных задач механики сплошной среды

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

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

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

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

Представлен параллельный метод построения призматической неструктурированной сетки, используемой для дискретизации сложных геологических структур при численном моделировании нефтяных и гидрогеоэкологических задач. Метод позволяет проводить адаптацию сетки к различным типам объектов (скважины сложной траектории, геологические разломы, пласты и т. п.) и обладает высокой скоростью построения Данный метод реализован в программном комплексе НИМФА.

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

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

Представлены результаты разработки алгоритма обнаружения нарушений оптического волокна, используемого для передачи информационных сигналов. Алгоритм реализован на языке ассемблер в программном обеспечении для микроконтроллера PIC16C717 и на языке С++ для микроконтроллера PIC18F1320, которые входят в состав системы мониторинга прямых потерь оптического волокна по информационным или пилот-сигналам.

Представлены результаты математического моделирования процесса тепломассопереноса при перекристаллизации интерметаллического образца TiAl(Nb,Cr,Zr) методом вертикальной бестигельной зонной плавки. Для расчетов применялась неравновесная математическая модель затвердевания многокомпонентных расплавов ! обобщение ранее предложенной модели затвердевания бинарных расплавов. Целью работы являлось установление соответствия между технологическими условиями процесса кристаллизации и микроструктурой получаемых в результате слитков

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

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

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

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

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

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

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

Статья посвящена реализации в рамках программы InterVal-3D возможности пересчета газодинамических величин с одной трехмерной регулярной (структурированной) сетки на другую с учетом восстановленной границы раздела между веществами в смешанных ячейках. Дается описание программы InterVal-3D, ее возможностей, а также принципов пересчета, положенных в основу программы. Описывается реализованный алгоритм восстановления границы раздела веществ. Приводится сравнительный анализ результатов пересчета с использованием алгоритма восстановления контактных границ и без него в модельных задачах с различными типами пространственных сеток.

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

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

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

С помощью методик ЭГАК и ТРЭК исследуется развитие турбулентного перемешивания и регулярного локального возмущения на границе раздела воздух!SF6 при взаимодействии с распространяющейся в ударной трубе волной с числами Маха в SF6 от 1,7 до 9,1. Рассмотрены случаи с локальным возмущением в виде треугольной канавки и полуцилиндра. Расчеты проводились с помощью прямого двумерного и трехмерного численного моделирования. Результаты расчетов сравниваются с известными экспериментальными данными.

Методика ТИМ-2D предназначена для решения задач механики сплошной среды на неструктурированных многоугольных лагранжевых сетках произвольного вида. В методике используется трехуровневое распараллеливание. На первом и втором уровнях, при распараллеливании счета по математическим областям и внутри математических областей по параобластям, используется модель распределенной памяти и интерфейс MPI. На третьем уровне осуществляется распараллеливание итераций счетных циклов в модели общей памяти с помощью интерфейса OpenMP. Уровни распараллеливания могут использоваться как по отдельности, так и в различных сочетаниях при решении одной задачи.

Приводятся результаты прямого численного моделирования пристенной турбулентности в плоском канале на расчетной сетке 257!129!129, выполненного на суперкомпьютере "Ломоносов". Расчеты проводились при числе Рейнольдса 5600, определенного по средней скорости потока и высоте канала. В качестве расчетной схемы выбрана схема КАБАРЕ без использования подсеточных моделей турбулентности, ранее хорошо себя зарекомендовавшая при моделировании затухания однородной изотропной турбулентности. Показано, что вычисленные характеристики пристенной турбулентности ! осредненный профиль скорости, компоненты тензора турбулентных напряжений, вклады молекулярной и турбулентной составляющих тензора вязкости и коэффициент сопротивления канала ! хорошо согласуются с расчетными данными, полученными при решении аналогичной задачи на аналогичных расчетных сетках в 1987г. группой Джона Кима из исследовательского центра в Эймсе (НАСА).

Рассматриваются особенности программной реализации многосеточного метода решения систем линейных алгебраических уравнений с агрегативным способом огрубления. Реализация метода выполнена в пакете программ ЛОГОС для модуля, ориентированного на численное решение уравнений гидродинамики вязкой несжимаемой жидкости на произвольных неструктурированных сетках. Рассматриваются варианты решения скалярных систем линейных алгебраических уравнений с использованием V -, W- и F-циклов многосеточного метода. Особое внимание уделяется распараллеливанию метода, показано решение проблемы грубых уровней. Приводится сравнение эффективности, полученной при расчетах течений вязкой несжимаемой жидкости с применением алгебраического многосеточного метода и предобусловленного метода сопряженных градиентов, а также многосеточного метода без использования глобального уровня.

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

Описывается прикладная программа EFRVIEWER. Программа предназначена для просмотра, редактирования и обработки данных, которые являются результатами численного моделирования физических процессов по прикладным программным комплексам РФЯЦ-ВНИИЭФ и представляются в виде бинарных файлов унифицированного формата ЕФР. Являясь своего рода сервисной утилитой, EFRVIEWER через графический интерфейс открывает пользователю доступ к двоичному коду задачи в удобной форме посредством гипертекста, таблиц и графиков

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

Описывается разрабатываемый в рамках контракта с Министерством образования и науки РФ отечественный программный комплекс S-MPI, включающий в себя библиотеку S-MPI и программные инструменты для анализа корректности и эффективности параллельных приложений. Приводится сравнение производительности S-MPI и двух других широко распространенных реализаций MPI-2 (Open MPI 1.5.4 и MVAPICH 1.2.7)

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

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

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

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

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

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

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

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