Том 27, № 1 (2019)
- Год: 2019
- Статей: 6
- URL: https://journals.rudn.ru/miph/issue/view/1243
- DOI: https://doi.org/10.22363/2658-4670-2019-27-1
Весь выпуск
Теория массового обслуживания
Тяжелая асимптотика исходящих вызовов для очереди повторных вызовов с двусторонней связью и множественными типами исходящих вызовов
Аннотация
В этой статье рассматривается модель очередей с одним сервером M|M|1|N с двумя типами вызовов: входящие вызовы и исходящие вызовы, когда входящие вызовы поступают на сервер в соответствии с процессом Пуассона. По прибытии входящий вызов немедленно занимает сервер, если он не используется, или выходит на орбиту, если сервер занят. С орбиты входящий вызов повторяет попытку занять сервер и ведет себя так же, как свежий входящий вызов. Сервер совершает исходящие звонки после экспоненциально распределенного простоя. Это можно интерпретировать как исходящие вызовы, поступающие на сервер в соответствии с процессом Пуассона. Существует N типов исходящих вызовов, длительность которых соответствует N различным экспоненциальным распределениям. Научная новизна заключается в получении асимптотики количества входящих вызовов в очереди повторных попыток в условиях высоких скоростей исходящих звонков и низкой продолжительности обслуживания каждого типа исходящих звонков. На основе полученной асимптотики построены аппроксимации распределения вероятностей количества входящих вызовов в системе.
Численное и имитационное моделирование
Параллельный алгоритм численного решения уравнения теплопроводности в сложной цилиндрической области
Аннотация
В этой статье представлен параллельный алгоритм для моделирования процесса теплопроводности внутри, так называемой, импульсной криогенной ячейки. Это моделирование важно для разработки устройства для порционной подачи рабочих газов в ионизационную камеру источника ионов. Моделирование основано на численном решении квазилинейного уравнения теплопроводности с периодическим источником в многослойной цилиндрической области. Для численного решения используется метод неявного направления (ADI). Из-за нелинейности уравнения теплопроводности был применен метод простой итерации. Для обеспечения сходимости итерационного процесса был реализован адаптивный временной шаг. Распараллеливание вычислений было реализовано с помощью прикладного программного интерфейса с разделяемой памятью OpenMP, и производительность параллельного алгоритма согласуется с примерами из литературы.
Задачи символьных вычислений, связанные с методами Рунге-Кутты и их решение в Sage
Аннотация
Схемы Рунге-Кутты играют очень важную роль в численном решении обыкновенных дифференциальных уравнений. В работе представлен пакет rk, являющийся подпрограммой Sage для вычисления матрицы Бутчера. Было проведено несколько численных экспериментов со стандартными и симплектическими схемами и проведена верификация с аналитическими результатами вычислений. Во-вторых, в Sage есть превосходные инструменты для исследования алгебраических множеств, основанные на методе базиса Грёбнера. Как известно, выбор параметров в схеме Рунге-Кутты свободен. С помощью этих инструментов получены алгебраические свойства многообразий в аффинном пространстве, координатами которых являются коэффициенты Бутчера в схеме Рунге-Кутты. Результаты приведены как для явной схемы Рунге-Кутты, так и для неявной схемы Рунге-Кутты с использованием разработанного пакета rk. Также приведены примеры для обоснования полученных результатов. Все расчеты выполнены в системе компьютерной алгебре Sage.
Программное обеспечение для численного решения дифференциальных уравнений в частных производных первого порядка
Аннотация
Дифференциальные уравнения с частными производными первого порядка, возникающие в прикладных задачах оптики и оптоэлектроники, часто содержат коэффициенты, которые не определяются одним аналитическим выражением во всей рассматриваемой области. Например, уравнение эйконала содержит показатель преломления, который описывается различными выражениями в зависимости от оптических свойств сред, которые заполняют рассматриваемую область. Этот тип уравнений не может быть проанализирован стандартными инструментами, встроенными в современные системы компьютерной алгебры, включая Maple. В статье рассматривается адаптация классического метода Коши интегрирования уравнений в частных производных первого порядка к случаю, когда коэффициенты уравнения являются различными аналитическими выражениями в подобластях G1,. , , , Gk, на которые делится рассматриваемый домен. В этом случае предполагается, что эти субдомены задаются неравенствами. Этот метод интеграции реализован как программа на Python с использованием библиотеки SymPy. Характеристики рассчитываются численно с использованием метода Рунге-Кутты, но с учетом изменения выражений для коэффициентов уравнения при переходе от одного субдомена к другому. Описаны основные функции программы, в том числе те, которые можно использовать для иллюстрации метода Коши. Проверка проводилась путем сравнения с результатами, полученными в системе компьютерной алгебры Maple.
Математическое моделирование
Математическая модель кавитации под воздействием одиночного импульса растяжения
Аннотация
В данной статье описана созданная математическая модель, позволяющая исследовать динамику кавитационных пузырьков под воздействием одиночного импульса отрицательного давления. Временная зависимость и координаты параметров несущей фазы, температуры и давления паровой фазы, концентрации и размера пузырьков определяются численно. Сделан вывод, что созданная модель дает хорошее согласие между расчетными и экспериментальными данными.
Метод объемных интегральных уравнений в задачах магнитостатики
Аннотация
В данной статье рассматриваются вопросы применения метода объемного интегрального уравнения к расчетам магнитных систем. Основным преимуществом этого подхода является то, что в этом случае нахождение решения уравнений сводится к области, заполненной ферромагнетиком. Сложность применения метода связана с особенностью ядра интегральных уравнений. По этой причине в методе коллокации используется только кусочно-постоянная аппроксимация неизвестных переменных в рамках элементов фрагментации внутри известного пакета GFUN3D. В качестве альтернативного подхода точки наблюдения могут быть заменены интегрированием по элементу фрагментации, что позволяет использовать приближение неизвестных переменных более высокого порядка. В представленной работе на примере обсуждаются основные аспекты применения этого подхода к моделированию магнитных систем. линейной аппроксимации неизвестных переменных: дискретизация исходных уравнений, разложение области вычисления на элементы, вычисление матричных элементов дискретной системы, решение полученной системы нелинейных уравнений. В рамках метода конечных элементов область расчета делится на набор тетраэдров. В начале начальная область аппроксимируется комбинацией макроблоков с предварительно построенной двумерной сеткой на их границах. После этого для каждого макроблока отдельно выполняется процедура построения сетки тетраэдра. При вычислении матричных элементов шестикратные интегралы по двум тетраэдрам сводятся к комбинации четырехкратных интегралов по треугольникам, которые рассчитываются по кубатурным формулам. Предлагается приведение сингулярных интегралов к комбинации регулярных интегралов методами, основанными на понятии однородных функций. Простые итерационные методы используются для решения нелинейных дискретизированных систем, что позволяет избежать обращения больших матриц. Результаты моделирования сравниваются с расчетами, полученными с использованием других методов.