Математика

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

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

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

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

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

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

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

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

Полностью дискретным (по времени и пространству) вариационным методом получена явная консервативная разностная схема с определением всех сеточных величин на одни и те же моменты времени с целочисленными индексами. Для определенности рассматриваются трехмерные течения

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

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

Предлагается неравновесная модель процесса кристаллизации двухкомпонентных расплавов, в которой численно решается задача Стефана с двумя границами: твердая фаза ! двухфазная переходная зона и двухфазная зона ! жидкая фаза. Двухфазная зона представлена пористой средой с переменной пористостью. Учитывается дополнительная сила сопротивления течению расплава в пористой среде по аналогии с законом Дарси. Сравниваются решения с линейной и квадратичной зависимостями объемной доли жидкой фазы от температуры. Приведены результаты расчетов эксперимента по кристаллизации расплава Sn с 20%-ным содержанием Pb методом направленной вертикальной кристаллизации, которая ведется сверху вниз по направлению вектора гравитации

В течение последних десяти лет в РФЯЦ-ВНИИТФ разрабатывается программа "Молох", используемая для проведения параллельных расчетов на суперкомпьютерах методом классической молекулярной динамики с высоким коэффициентом распараллеливания. Был разработан набор тестов для анализа эффективности реализованных алгоритмов. Он может быть использован для сравнения различных молекулярнодинамических комплексов

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

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

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

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

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

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

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

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

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

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

Рассматривается сервисная часть пакета программ ЛЭГАК-3D, предназначенная для подготовки задач, сопровождения расчетов и постобработки результатов счета. Возможности представленной системы, помимо визуализации в различных видах и редактирования всех компонентов задачи, включают высокоуровневые операции над этими компонентами (конструирование) ! перенос компонентов из одной задачи в другую с добавлением или заменой имеющихся. Анализ результатов счета задач обеспечивается как операциями контекстного меню в процессе отображения данных, так и отдельными инструментальными средствами (в том числе с использованием внешних приложений)

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

Разработан алгоритм расчета сегрегации примеси при кристаллизации двухкомпонентных расплавов, характеризующихся большими размерами двухфазной области. При расчете конвекции в пористой среде учитывается добавочное сопротивление течению по аналогии с законом Дарси. Температура ликвидус линейно зависит от концентрации примеси в жидкой фазе Предложенная модель реализована в программном комплексе GIGAN-2D. Для верификации модели просчитаны два численных бенчмарка: по кристаллизации расплава Sn c 10%-ным содержанием Pb и расплава Pb с 18%-ным содержанием Sn.

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

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

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

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

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

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

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

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

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

Создан код OptModel, предназначенный для решения стандартных оптико-модельных задач, связанных с упругим рассеянием n, p, d, t, 3,4,6He, 6Li на сферических ядрах (или близких к ним). В коде используется феноменологический оптический потенциал на основе формы Вудса – Саксона. Для решения радиальных уравнений Шредингера впервые применен метод 12-го порядка точности (разработан в 2004 г. в Шанхайском университете, Китай), значительно дополненный и усовершенствованный. Расчет кулоновских функций ведется по программе RCFWN с относительной точностью 10–10, язык программирования C++ В коде используется современный графический интерфейс.

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

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

Предлагается модификация алгоритма Боуэра – Уотсона построения треугольной неструктурированной топологически двумерной сетки в областях произвольной формы и связности. Алгоритм состоит из четырех частей: 1) триангуляция области, ограниченной исходным контуром; 2) минимизация геометрической меры нерегулярности; 3) минимизация топологической меры нерегулярности; 4) измельчение полученной сетки (необязательное)

Полностью дискретным (по времени и пространству) вариационным методом получена явная консервативная разностная схема с определением всех сеточных величин на одни и те же моменты времени с целочисленными индексами. Для определенности рассматриваются трехмерные течения

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

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

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

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

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

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

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

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

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

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