Архив статей журнала
In this paper, we consider the Fredholm integral equations of the second kind and construct a new iterative scheme associated to the Nyström method, which was elaborated by Atkinson to approximate the solution over a large interval. Primarily, we demonstrate the inability to generalize the Atkinson iterative methods. Then, we describe our modified generalization in detail and discuss its advantages such as convergence of the iterative solution to the exact solution in the sense norm of the Banach space С0[a,b]. Finally, we give a numerical examples to illustrate the accuracy and reliability of our generalization.
В статье рассматривается параллельный алгоритм решения задач глобальной оптимизации и обсуждается его реализация с использованием набора инструментов Intel oneAPI. Предполагается, что целевая функция задачи задана как “черный ящик” и удовлетворяет условию Липшица. Изложенный в статье параллельный алгоритм использует схему редукции размерности на основе кривых Пеано, которые непрерывно и однозначно отображают отрезок вещественной оси на гиперкуб. В качестве средства для реализации параллельного алгоритма использован инструментарий Intel oneAPI, который позволяет писать один код как для центрального процессора, так и для графических ускорителей. Приведены результаты вычислительных экспериментов, полученные при решении серии сложных задач многоэкстремальной оптимизации.
На сегодняшний день многопроходный метод PIV (Particle Image Velocimetry) широко используется в области экспериментальной механики жидкости и газа из-за его высокой надежности при решении практических задач. Однако он имеет известное ограничение, связанное с ошибками, возникающими при вычислении производных скорости, необходимых для деформации обрабатываемых PIV-изображений при повышении производительности метода. Поскольку количество ошибок увеличивается с применением схем более высокого порядка, на практике чаще всего ограничиваются первым порядком, что в свою очередь приводит к снижению пространственного разрешения. В данной работе предлагается метод, допускающий применение схем более чем второго порядка, что позволяет заметно повысить точность измерения скорости и ее производных и тем самым увеличить пространственное разрешение. Метод не требует восстановления ошибочных векторов скорости, позволяет избежать численного расчета производных скорости и легко применим на практике.
Исследованы закономерности процессов испарения и конденсации чистого пара в методе решеточных уравнений Больцмана. Выполнено моделирование этих процессов при постоянных во времени потоках пара на границе расчетной области. Показано, что в этом случае осуществляются квазистационарные режимы испарения и конденсации. Предложен простой численно эффективный метод задания потока пара на плоской границе расчетной области путем вычисления функций распределения на входящих характеристиках метода решеточных уравнений Больцмана. В расчетах показано, что поток массы при испарении плоской поверхности пропорционален разности плотностей насыщенного и окружающего пара при данной температуре поверхности, что хорошо согласуется с законом Герца-Кнудсена. Результаты трехмерного и одномерного моделирования методом решеточных уравнений Больцмана совпадают с высокой точностью. Показано, что отношение разности плотностей к потоку вещества на границе фаз при заданной температуре линейно зависит от времени релаксации как для испарения, так и для конденсации. Исследовано влияние температуры на интенсивность потоков испарения и конденсации чистого пара. Обнаружена зависимость процессов испарения и конденсации от времени релаксации, которое определяет кинематическую вязкость флюида.
Процесс распараллеливания программ может быть затруднён ввиду их оптимизации под последовательное выполнение. Из-за этого полученная параллельная версия может быть неэффективной, а в некоторых случаях распараллеливание оказывается невозможным. Решить указанные проблемы помогают преобразования исходного кода программ. В данной статье рассматривается реализации в системе автоматизированного распараллеливания SAPFOR (System FOR Automated Parallelization) преобразований последовательных Фортран-программ, позволяющих облегчить работу пользователя в системе и существенно снизить трудоемкость распараллеливания программ. Применение реализованных преобразований в системе SAPFOR продемонстрировано на прикладной программе, решающей систему нелинейных дифференциальных уравнений в частных производных. Также было произведено сравнение производительности полученной параллельной версией с версиями, распараллелеными вручную с использованием DVM и MPI технологий.
Работа посвящена построению параллельных алгоритмов решения прямой начально-краевой задачи и обратной задачи о восстановлении правой части для уравнения диффузии с дробной производной по времени. При использовании дополнительной информации о решении в некоторой внутренней точке обратная задача сводится к прямой задаче для вспомогательного уравнения. После применения конечно-разностных схем задачи сводятся к решению систем линейных алгебраических уравнений. Разработанные алгоритмы основаны на методе параллельной прогонки и реализованы для многоядерных процессоров с использованием технологии OpenMP. Проведены численные эксперименты для исследования производительности разработанных алгоритмов.
Высокотемпературные эффекты оказывают существенное влияние на характеристики летательных аппаратов, движущихся с гиперзвуковой скоростью. В связи со сложностью постановки физического эксперимента, методы математического моделирования играют важную роль для нахождения характеристик гиперзвуковых летательных аппаратов. Обсуждается построение и реализация математической модели, предназначенной для численного моделирования гиперзвукового обтекания тела с учетом неравновесных физико-химических процессов в высокотемпературном воздухе. Математическая модель включает в себя уравнения газовой динамики, уравнения модели турбулентности и уравнения химической кинетики. Проводится численное моделирование сверх- и гиперзвукового обтекания полусферы потоком воздуха с учетом высокотемпературных эффектов. Приводится критический обзор различных моделей, которые применяются для нахождения расстояния от фронта ударной волны до поверхности сферы. Результаты расчетов, полученные с использованием разработанного численного метода, сравниваются с данными физического эксперимента и расчетными данными, имеющимися в литературе, в широком диапазоне чисел Маха набегающего потока. Разработанная модель и результаты расчетов имеют значение для моделирования обтекания тел сложной конфигурации и проектирования высокоскоростных летательных аппаратов.
The objective of this paper is to construct a generalized quadratic spectrum approximation based on the Kantorovich projection method which llows us to deal with the spectral pollution problem. For this purpose, we prove that the property U (see Eq. 3) holds under weaker conditions than the norm and the collectively compact convergence. Numerical results illustrate the effectiveness and the convergence of our method.
В данной работе начальная форма волны цунами (ниже именуемая источником цунами) представляется как решение обратной задачи математической физики на основе инверсии удаленных записей пришедшей волны, что позволяет детально изучить факторы, влияющие на результаты восстановления. Исследуемая задача является некорректной, что приводит к ожидаемой неустойчивости численного решения, существенно уменьшить которую позволяет регуляризация, основанная на методе усеченного сингулярного разложения (SVD) (далее метод r-решения). В рамках предложенного подхода на основе анализа распространения энергии олны предлагается методика выбора наиболее информативной части имеющейся системы наблюдения для реального события цунами на Соломоновых островах 6 февраля 2013 г. Метод может быть полезен при разработке новых систем мониторинга цунами.
Работа связана с изучением нелинейных параболических систем, возникающих при моделировании и управлении нестационарными процессами фильтрации в подземной гидродинамике. Одна из постановок является системой, которая включает в себя краевую задачу второго рода для квазилинейного параболического уравнения с неизвестной функцией источника в правой части, а также уравнение изменения по времени этой функции. В другой постановке рассматривается проблема управления этой системой с управляющим воздействием граничного режима. Данные постановки существенно отличаются от обычных краевых задач и задач управления для параболических уравнений, в которых предполагается, что все входные данные должны быть заданы. Полученные в работе результаты представляют не только теоретический интерес, они имеют практическое значение для исследования различных фильтрационных процессов. Приведены некоторые примеры таких приложений, связанных с движением жидкости в трещиновато-пористых средах.
В работе рассматривается предобусловливатель блочного неполного обратного треугольного разложения первого порядка “по значению” BIIC-IC1 для решения систем линейных алгебраических уравнений с симметричной положительно определенной матрицей. Рассматривается способ применения MPI+OpenMP технологии для построения и обращения предобусловливателя BIIC-IC1, при этом в предобусловливателе число блоков кратно числам используемых процессоров и используемых потоков. Предлагается способ применения MPI+OpenMP технологии для построения и обращения предобусловливателя BIIC-IC1, в котором для применения OpenMP технологии используется специальное упорядочение узлов сетки внутри подобластей, соответствующих расчетам на процессорах. Проводится сравнение времени решения задач методом сопряженных градиентов с предобусловливателем BIIC-IC1 с использованием MPI и гибридной MPI+OpenMP технологии на примере модельной задачи и ряда задач из коллекции разреженных матриц SuiteSparse.
Разработан алгоритм высокоточного численного решения эллиптического уравнения второго порядка при наличии в области нескольких интерфейсов, в том числе пересекающихся и невыпуклых. Для аппроксимации задачи в окрестности интерфейсов используются нерегулярные ячейки (н-ячейки), отсекаемые ими от регулярных ячеек прямоугольной сетки, и законтурные части этих ячеек. Для построения приближенного решения предложено: 1) выписывать дополнительные условия согласования в н-ячейках на интерфейсах, увеличивая количество согласуемых ячеек вблизи интерфейсов; 2) уменьшать общую часть интерфейса, заключенную в соседних ячейках и используемую для записи условий. Для решения краевой задачи Дирихле реализован hp-вариант метода коллокации и наименьших квадратов (hp-МКНК) в сочетании с современными алгоритмами ускорения итерационного процесса: предобуславливание; распараллеливание с помощью OpenMP; ускорение, основанное на подпространствах Крылова; многосеточный алгоритм. При решении различных тестовых задач исследованы сходимость hp-МКНК и обусловленность возникающих переопределенных систем линейных алгебраических уравнений (СЛАУ). Проведено сравнение результатов, полученных МКНК, с результатами других авторов, использовавших метод MIB (англ. matched interface and boundary).