Все выпуски

В данной работе решается задача поиска конструкции магнитной системы для создания магнитного поля с требуемыми характеристиками в заданной области. На основе анализа математической модели магнитной системы предлагается достаточно общий подход к решению нелинейной обратной задачи, которая описывается уравнением Фредгольма H(z) = ∫_SIJ(s)G(z, s)ds, z ∈ S _H, s ∈ S _I . Необходимо определить распределение плотности тока J(s), а также расстановку источников тока для создания поля H(z). В работе предлагается метод решения этих задачс помощью регуляризованных итерационных процессов. На примере конкретной магнитной системы проводится численное исследование влияния различных факторов на характер создаваемого магнитного поля.

Проведено математическое моделирование нестационарных режимов естественной конвекции в замкнутой пористой цилиндрической полости с теплопроводной оболочкой конечной толщины в условиях конвективного теплообмена с внешней средой. Краевая задача математической физики, сформулированная на основе модели Дарси–Буссинеска в безразмерных переменных «функция тока – температура», реализована численно методом конечных разностей. Детально проанализировано влияние проницаемости пористой среды 10^–5≤Da<∞, отношения толщины твердой оболочки к внутреннему радиусу цилиндра 0.1≤h/L≤0.3, относительного коэффициента теплопроводности 1≤λ_1,2≤20 и безразмерного времени 0≤τ≤1000 как на локальные распределения изолиний функции тока и температуры, так и на интегральные комплексы, отражающие интенсивность конвективного течения и теплопереноса.

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

В статье приводятся результаты моделирования деформации метеорита при попадании в плотные слои атмосферы разрывным методом Галёркина на неструктурированных треугольных сетках и методом сглаженных частиц. В качестве исходных данных брались материалы по Челябинскому метеориту. Проводилась серия расчётов, где варьировались характеристики материала метеорита и угол входа в плотные слои атмосферы.

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

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

Построена статическая трехуровневая теоретико-игровая модель системы управления судовыми балластными водами. Используются методы иерархического управления при одновременном учете условий поддержания системы в заданном состоянии. Проводится сравнение результатов исследования модели с точки зрения игр Гермейера $\Gamma_1$ и $\Gamma_2$. Приведены примеры численных расчетов в ряде характерных случаев.

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

Рассматривается математическая модель, описывающая конкуренцию за неоднородный ресурс двух близкородственных видов на одномерном ареале. Распространение популяций определяется диффузией и направленной миграцией, а рост подчиняется логистическому закону. Исследуются решения соответствующей начально-краевой задачи для нелинейных уравнений параболического типа с переменными коэффициентами (функция ресурса, параметры роста, диффузии и миграции). Для анализа формирования популяционных структур применяется подход на основе теории косимметричных динамических систем В. И. Юдовича. Аналитически получены условия на параметры системы, при выполнении которых у системы имеется нетривиальная косимметрия. В численном эксперименте подтверждено возникновение непрерывного семейства стационарных решений при выполнении условий существования косимметрии. Расчетная схема основана на конечно-разностной дискретизации по пространственной переменной с использованием интегро-интерполяционного метода и интегрировании по времени методом Рунге–Кутты. Далее численно исследовано влияние параметров диффузии и миграции на пространственно-временные сценарии развития популяций. В окрестности многообразия, соответствующего косимметрии задачи, рассчитаны нейтральные кривые диффузионных параметров, отвечающих границам устойчивости решений с одной популяцией. Для ряда значений параметров миграции и функций ресурса с одним и двумя максимумами построены карты областей параметров, которые соответствуют различным сценариям сосуществования и вытеснения видов. В частности, найдены области параметров, при которых выживание того или иного вида определяется условиями начального размещения. Отмечено, что реализуемая при этом динамика может быть нетривиальна: после начального снижения плотностей обоих видов наблюдается последующий рост одной популяции и убывание другой. Проведенный анализ показал, что области диффузионных параметров, отвечающих различным сценариям формирования популяционных структур, группируются вблизи линий, соответствующих косимметрии рассматриваемой математической модели. Полученные карты позволяют объяснить медленную динамику системы близостью к косимметричному случаю и дать трактовку эффекта выживания популяции за счет изменения диффузионной мобильности при исчерпании ресурса.

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

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

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