ВЫЧИСЛИТЕЛЬНЫЕ МЕТОДЫ И ПРОГРАММИРОВАНИЕ

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

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

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

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

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

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

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

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

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