Математика

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

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

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

Представлены результаты исследования возможностей создания программноаппаратных комплексов на базе вычислительных систем, интегрирующих процессоры архитектуры х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 в пакете уделяется распараллеливанию метода, показано решение «проблемы грубых уровней». программ ЛОГОС Приводится сравнение эффективности, полученной при расчетах течений вязкой несжимаемой жидкости с применением алгебраического многосеточного метода и предобусловленного метода сопряженных градиентов, а также многосеточного метода без использования глобального уровня.

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

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

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

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

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

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

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

Актуальность и цели. Распределение Гаусса естественно возникает во многих приложениях и широко используется в различных теоретических построениях. Важную роль играет и нижняя срезка Q(x) несобственного интеграла от плотности стандартного гауссова распределения. Целью данной работы является получение оценок сверху для произвольной степени функции Q(x) через несобственный интеграл того же вида с нижней границей ax, где a – некоторый параметр. Материалы и методы. Для получения необходимых оценок изучалось по- ведение разности Qm(x)−Q(mx)на различных интервалах числовой оси, при этом широко использовались хорошо известные свойства гауссова распределения. Кроме того, были выведены точные неравенства для показательной функции специального вида и получены оценки сверху и снизу функции Q(x). Результаты. В работе показано, что для любого действительного x (при m > 2 ) выполняется неравенство Qm(x)

Актуальность и цели. Булевы и многозначные функции – основной объект изучения дискретной математики. Они представляют собой зависимости между величинами, принимающими конечный набор значений. Существует несколько способов описания таких зависимостей, и на практике часто встречается табличное задание функции и задание в виде полинома. Оба эти представления функций можно выразить в виде векторов. В случае табличного задания функции это вектор ее значений, в случае полиномиального задания – вектор коэффициентов полинома. Преобразование вектора значений функции в вектор коэффициентов ее полинома в булевом случае является преобразованием Мёбиуса. Неподвижные точки такого преобразования мы будем называть стационарными функциями. Пусть α – вектор, состоящий из n элементов поля E3 . α-преобразованием функции f будем называть такую функцию g =ν (f), что g(x,,x)=f(x+α,,x+α).Если ν  (f)= f , то такую

Актуальность и цель. В статье К. И. Бабенко «О некоторых задачах теории приближений и численного анализа»1 среди ряда важных проблем вычислительной математики были сформулированы две проблемы: 1) вычисление поперечников Колмогорова и Бабенко на классе Q(Ω,M)(класс Q(Ω,M)со-