Вычислительная математика, численный анализ

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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