Все выпуски

Статья посвящена моделированию гидродинамических процессов мелководных водоемов на примере Азовского моря. В статье приведена математическая модель гидродинамики мелководного водоема, позволяющая вычислить трехмерные поля вектора скорости движения водной среды. Применение регуляризаторов по Б.Н. Четверушкину в уравнении неразрывности привело к изменению способа расчета поля давления, базирующегося на решении волнового уравнения. Построена дискретная конечно-разностная схема для расчета давления в области, линейные размеры которой по вертикали существенно меньше размеров по горизонтальным координатным направлениям, что является характерным для геометрии мелководных водоемов. Описаны метод и алгоритм решения сеточных уравнений с предобуславливателем трехдиагонального вида. Предложенный метод применен для решения сеточных уравнений, возникающих при расчете давления для трехмерной задачи гидродинамики Азовского моря. Показано, что предложенный метод сходится быстрее модифицированного попеременно-треугольного метода. Представлена параллельная реализация предложенного метода решения сеточных уравнений и проведены теоретические и практические оценки ускорения алгоритма с учетом времени латентности вычислительной системы. Приведены результаты вычислительных экспериментов для решения задач гидродинамики Азовского моря с использованием гибридной технологии MPI + OpenMP. Разработанные модели и алгоритмы применялись для реконструкции произошедшей в 2001 году в Азовском море экологической катастрофы и решения задачи движения водной среды в устьевых районах. Численные эксперименты проводились на гибридном вычислительном кластере К-60 ИПМ им. М.В. Келдыша РАН.

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

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

Для квантового описания поведения двуядерных систем на начальной стадии околобарьерного слияния тяжелых ядер использованы численные методы нахождения коллективных и одночастичных состояний. Коллективные возбужденные состояния в таких системах представляют собой согласованные колебания поверхностей сферических ядер. Одночастичные состояния внешних нейтронов аналогичны состояниям валентных электронов двухатомных молекул.

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

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

Представлена программа NINE (Newtonian Iteration for Nonlinear Equation) численного решения граничных задач для нелинейных дифференциальных уравнений второго порядка на основе непрерывного аналога метода Ньютона (НАМН) с использованием нумеровской конечно-разностной аппроксимации четвертого порядка относительно шага дискретизации по пространственной переменной. Обсуждаются алгоритмы вычисления ньютоновского итерационного параметра. Выполнены методические расчеты, демонстрирующие влияние выбора итерационного параметра на сходимость итерационного процесса. Представлены результаты проведенного с помощью программы NINE численного исследования положительных частицеподобных решений уравнения скалярного поля.

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

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

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