Все выпуски

Пусть зафиксирован некоторый класс возмущений матрицы коэффициентов $A(\cdot)$ дискретной линейной однородной системы вида $$x(m+1)=A(m)x(m),\quad m\in\mathbb Z,\quad x\in\mathbb R^n,$$ с вполне ограниченной на $\mathbb Z$ матрицей $A(\cdot)$. Спектральным множеством этой системы, отвечающим заданному классу возмущений, называем совокупность полных спектров показателей Ляпунова возмущенных систем, когда возмущения пробегают весь заданный класс. Основное внимание в работе уделено классу ${\cal R}$ возмущенных систем вида $$y(m+1)=A(m)R(m)x(m),\quad m\in\mathbb Z,\quad y\in\mathbb R^n,$$ с вполне ограниченными на $\mathbb Z$ матрицами $R(\cdot)$, и его подклассам ${\cal R}_{\delta}$ с матрицами $R(\cdot)$, удовлетворяющими оценке $\sup_{m\in\mathbb Z}\|R(m)-E\|<\delta$, где $\delta>0$. Доказано, что если показатели Ляпунова исходной системы устойчивы, то спектральное множество $\lambda({\cal R})$, отвечающее классу ${\cal R}$, совпадает с множеством всех упорядоченных по возрастанию наборов из $n$ чисел, при этом для каждого $\Delta>0$ существует такое $\ell=\ell(\Delta)>0$, что для любого $\delta<\Delta$ спектральное множество $\lambda({\cal R}_{\ell\delta})$ содержит в себе $\delta$-окрестность полного спектра показателей Ляпунова невозмущенной системы.

Рассматривается дискретная линейная однородная система

$$x(m+1)=A(m)x(m),\quad m\in\mathbb Z,\quad x\in\mathbb R^n, \qquad\qquad (1)$$

с вполне ограниченной матрицей $A(\cdot)$ и полным спектром показателей Ляпунова $\lambda_1(A)\leqslant\ldots\leqslant\lambda_n(A)$. Показатели Ляпунова системы (1) называются устойчивыми, если для любого $\varepsilon>0$ найдется такое $\delta>0$, что для всякой вполне ограниченной на $\mathbb N$ $n\times n$-матрицы $R(\cdot)$, удовлетворяющей оценке $\sup_{m\in\mathbb N}\|R(m)-E\|<\delta$, для полного спектра показателей Ляпунова $\lambda_1(AR)\leqslant\ldots\leqslant\lambda_n(AR)$ возмущенной системы

$$z(m+1)=A(m)R(m)z(m),\quad m\in\mathbb Z,\quad x\in\mathbb R^n,$$

справедливо неравенство $\max_{j=1,\ldots,n}|\lambda_j(A)-\lambda_j(AR)|<\varepsilon$. В работе построен пример системы вида (1) с неустойчивыми показателями Ляпунова.

Ряд задач в теории характеристических показателей Ляпунова линейных дифференциальных систем

ẋ=A(t)x,    x∈Rⁿ,    t≥0,

сводится к изучению влияния возмущений коэффициентов на характеристические показатели и другие асимптотические инварианты возмущенных систем

ẏ=A(t)y+Q(t)y,    y∈Rⁿ,    t≥0.

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

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

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

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

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

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

В работе рассмотрен локально-неравновесный процесс затвердевания переохлажденного бинарного расплава. В целях простоты предполагается, что затвердевающая бинарная система находится при постоянных температуре и давлении и имеет две фазы, соответствующие твердому и жидкому состояниям. Математическое описание процесса затвердевания основано на модели фазового поля, обобщающей подход Плаппа (M. Plapp, Phys. Rev. E 84, 031601 (2011)) на случай локально-неравновесных процессов. Для вывода термодинамически согласованных уравнений модели использован метод расширенной необратимой термодинамики в отличие от феноменологического подхода Плаппа. Другое различие с моделью Плаппа состоит в использовании в качестве динамической переменной концентрации, а не химпотенциала примеси. В рамках полученной модели показана эквивалентность описания процесса затвердевания через концентрационное поле и через химпотенциал системы. В силу малости времен релаксации представленная модель сводится к сингулярно-возмущенной системе уравнений в частных производных параболического типа, описывающих динамику фазового и концентрационного полей. В работе предполагается известным описание термодинамических равновесных состояний на основе экспериментально полученных потенциалов Гиббса.

Для проверки полученной модели проведено численное моделирование одномерной задачи затвердевания в приближении разбавленного расплава Si-As, ранее неоднократно исследовавшегося экспериментально. Чтобы численно решить систему сингулярно-возмущенных уравнений, в работе предложен градиентно-устойчивый явный метод интегрирования уравнений второго порядка точности по времени. Для сведения бесконечного пространственного интервала к конечному использован метод «периодического сдвига». Оценка устойчивости получена из численных экспериментов.

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

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