Все выпуски

Предложен бессеточно-спектральный метод расчета динамики плоских вихревых течений невязкой несжимаемой жидкости в геофизических приближениях с учетом планетарного вращения. Математически задача описывается системой двух уравнений в частных производных относительно функций тока и завихренности с различными граничными условиями (замкнутая область течения и периодические условия). В основе метода лежат следующие положения: поле завихренности задано значениями на множестве частиц; функция завихренности приближается с помощью кусочно-непрерывной аппроксимации кубическими полиномами от двух пространственных переменных; коэффициенты полиномов находятся методом наименьших квадратов; функция тока на каждом временном шаге находится методом Бубнова–Галёркина; динамика жидких частиц рассчитывается псевдосимплектическим методом Рунге–Кутты. В статье впервые подробно описан вариант метода для периодических граничных условий. Адекватность численной схемы проверена на тестовых примерах.

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

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

Данная работа лежит в области неравновесной статистической механики. Она посвящена изучению пространственно-временного поведения изолированных многочастичных систем. Сформулирована феноменологическая модель, сочетающая законы классической динамики с присутствием случайных переходов между различными микроскопическими состояниями системы. Все процессы, происходящие с такой системой, имеют место на замкнутой гиперповерхности в фазовом пространстве координат-импульсов. Эта поверхность соответствует фиксированному значению энергии системы; она обозначена как ${\mathrm Surf}_{\mathrm E}^{}$. Кинетическое уравнение для плотности вероятностей распределения системы в фазовом пространстве $f(\vec{x},\,t)$ приведено к форме, удобной для нахождения решений этого уравнения ($\vec{x}$ — совокупность фазовых переменных). Дано доказательство того, что экстремум энтропии системы, который находится методом вариационного исчисления, является максимумом. Вблизи своего максимума энтропия пропорциональна интегралу от квадрата вариации $f(\vec{x},\,t)$ по всему объему гиперповерхности ${\mathrm Surf}_{\mathrm E}^{}$. Найдена зависимость энтропии от вариации плотности распределения $f(\vec{x},\,t)$, когда эта вариация имеет постоянные значения в двух подобластях гиперповерхности, причем эти значения могут быть весьма велики. Сформулировано кинетическое уравнение для $D(\vec{x},\,t)$ — отклонения плотности распределения частиц $f(\vec{x},\,t)$ от ее равновесного значения. Найдены аналитические решения этого уравнения для нескольких вариантов зависимости скорости случайных переходов от расстояния между начальным и конечным микроскопическими состояниями перехода. Если эта скорость не зависит от вышеназванного расстояния, то исходная форма $D(\vec{x},\,t)$ сохраняется в течение всего процесса экспоненциального затухания величины $D(\vec{x},\,t)$. Рассмотрен переходный процесс в более общем случае, когда случайные переходы возможны только в ограниченном диапазоне расстояний между микросостояниями. Тогда величина $D(\vec{x},\,t)$ равна сумме пространственных мод, каждая из которых при стремлении плотности распределения к равновесному значению затухает со своей скоростью. Чем выше порядок моды, тем быстрее она затухает. Поэтому в данном случае переходный процесс приводит не только к общему уменьшению $D(\vec{x},\,t)$, но и к сглаживанию первоначальной формы $D(\vec{x},\,t)$. Найдена общая зависимость разности между равновесным и неравновесным значениями энтропии от величины $D(\vec{x},\,t)$.