Предложена стратегия построения обучающего набора данных для подавляющей численную дисперсию нейронной сети NDM-net (numerical dispersion mitigation network), заключающаяся в расчете полного набора сейсмограмм методом конечных разностей на грубой сетке и в расчете обучающей выборки с применением более мелкой сетки. Обучающая выборка представляет собой малый набор сейсмограмм с определенным пространственным размещением источников волнового поля. После обучения сеть NDM-net позволяет аппроксимировать низкокачественные сейсмограммы, рассчитанные на грубой сетке, в сейсмограммы с меньшим шагом дискретизации. Оптимизация процесса построения репрезентативной обучающей выборки сейсмограмм основана на минимизации метрики Хаусдорфа между обучающей выборкой и полным набором сейсмограмм. Применение нейронной сети NDM-net позволяет уменьшить временные затраты при расчетах волновых полей на мелкой сетке.
В работе рассматриваются вопросы численного моделирования двухфазных течений с применением эйлеровой многоскоростной модели типа Баера-Нунциато. Представлено описание математической модели, подробно рассмотрен численный алгоритм решения задачи на основе разрывного метода Галеркина. Представлено описание разработанного программного комплекса, основное назначение которого - математическое моделирование двухфазных течений с прямым разрешением динамики границ раздела фаз. Особенностью предложенных алгоритмов является лимитирование простых и консервативных переменных, гарантирующее как отсутствие нефизичных осцилляций, так и допустимые значения физических полей с применением лимитера WENO-S. Основная цель работы - дать исчерпывающее описание предложенного комплекса алгоритмов для решения задач рассматриваемого класса и пригодного для решения задач в реалистичных постановках. На примере модельной задачи демонстрируется возможность его применения для решения задач на сетках больших размеров.
В работе рассматривается уравнение типа Бюргерса с полиномиальной нелинейностью и нулевыми краевыми условиями. Для интересующего диапазона параметров тождественно нулевое решение задачи является локально неустойчивым, и в его окрестности существует устойчивое многообразие, имеющее конечную коразмерность. Для приближенного построения указанного многообразия предложен комбинированный итерационный алгоритм, начальное условие для которого строится аналитическим методом и имеет квадратичную точность. Численно показано, насколько существенно данная модификация позволяет уменьшить для типичных значений параметров вычислительную сложность проецирования на искомое многообразие по сравнению со стандартным линейным приближением. Полученные результаты допускают обобщение на многомерные диссипативные уравнения широкого класса и могут применяться при решении задач асимптотической стабилизации по начальным данным, краевым условиям и правой части.
В статье рассматривается явно-неявная балансно-характеристическая схема CABARETI-NH (CABARET Implicit Non-Hydrostatic), основанная на схеме КАБАРЕ, для решения гиперболизированной системы уравнений Навье-Стокса. Неявность вдоль одного пространственного направления позволяет значительно увеличить шаг по времени на вычислительных сетках с большим аспектным отношением ячеек. Для разрешения введенной неявности используется метод гиперболической прогонки. Это позволяет сохранить вычислительную эффективность алгоритма на уровне явных схем. Приводятся результаты валидации модели на лабораторном эксперименте трехмерного гравитационного течения стратифицированной жидкости.
В работе изложены алгоритмы и приведены компактные программные модули на языке С для быстрого вычисления показательной функции с помощью таблиц для процессоров архитектуры x86-64. Выполнена оценка точности и проведено сравнение быстродействия для некоторых процессоров AMD и Intel. Реализовано и протестировано обобщение табличного подхода для некоторых тригонометрических функций. В среднем предложенные функции работают в 10 раз быстрее соответствующих аналогов из стандартной математической библиотеки с прототипами в math.h.
Проведены расчеты тестовой задачи, связанной с моделированием течения в идеализированном медицинском устройстве, в программном комплексе FlowVision. Расчеты проводились для ламинарного, турбулентного и переходного режимов течения. Исследована масштабируемость задачи. На основе решения тестовой задачи сделан вывод о возможности применения программного комплекса FlowVision к решению проблем гемодинамики.
A recurrent formula for estimating an impedance of the cathode catalyst layer with fast oxygen transport in a polymer electrolyte fuel cell is derived. The catalyst layer is divided into N sub-layers and application of the charge conservation law enables to obtain the nonlinear recurrent relation Zn = f(Zn-1), where Zn is the accumulated impedance of all sub-layers up to the n-th one. Numerical solution of this relation gives the total impedance of the catalyst layer with taking into account variation of the static overpotential of the oxygen reduction reaction along the depth of this layer. The model is simple, robust and two orders of magnitude faster than the standard model based on numerical solution of the differential equation.
В статье рассматривается задача планирования работ на многоядерном процессоре с учетом их замедления при совместном выполнении. Предложена постановка задачи и модель частично целочисленного линейного программирования, доказана NP-трудность задачи при числе ядер, ограниченном константой. Результаты планировщика Intel TBB и жадного алгоритма сравниваются с результатами, полученными в соответствии с предложенной моделью с помощью пакета CPLEX. Проведенный эксперимент показал преимущества предложенного подхода по времени завершения всех работ.
Проект создания цифровой платформы Algo500 направлен на решение задачи совместного анализа свойств алгоритмов и особенностей архитектур суперкомпьютеров. В статье на основе методологии онтологического анализа рассматриваются и предлагаются понятия, модели и метамодели данных, обосновываются подходы к описанию некоторых понятий из мира высокопроизводительных вычислений (HPC), устанавливаются новые требования к моделям данных, которые должны обеспечить выполнение задач, поставленных при создании платформы Algo500.
Представлен алгоритм численного моделирования квазистатической нагрузки пористого флюидонасыщенного образца для решения задачи апскейлинга трещиновато-пористой флюидонасыщенной среды. В результате апскейлинга восстанавливается анизотропная однородная вязкоупругая эквивалентная среда, определяемая комплекснозначным частотно-зависимым тензором жесткости. Компоненты восстановленного тензора жесткости эквивалентной среды используются для оценки частотно-зависимого сейсмического затухания и фазовой скорости сейсмических волн. Численный апскейлинг включает в себя численное решение краевой задачи для системы уравнений Био анизотропной пороупругой флюидонасыщенной среды в частотной области для набора частот и различных граничных условий. Численное решение системы уравнений Био основано на конечно-разностной аппроксимации уравнений в квазистатической постановке и проводится с использованием прямого решателя результирующей СЛАУ. Используемый прямой решатель позволяет эффективно решать СЛАУ для набора правых частей, что необходимо при численном апскейлинге. Представленная реализация алгоритма позволяет численно решать двумерную задачу апскейлинга на расчетной сетке с размерами до 2000 × 2000 узлов на персональном компьютере, что обеспечивает возможность восстанавливать эквивалентную вязкоупругую модель для детализированных репрезентативных образцов трещиновато-пористой среды. Для демонстрации применимости алгоритма выполнено несколько наборов численных экспериментов, направленных на выявление влияния связности трещин и микромасштабной анизотропии пороупругого материала внутри трещин на затухание, вызванное индуцированными волной потоками флюидов, и дисперсию сейсмической волны в сложной трещиновато-пористой флюидонасыщенной среде.