Все выпуски

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

The article is devoted to the development of Lyapunov methods for analyzing the instability of the equilibrium of a dynamical system in the space of probability measures, given by the nonlocal continuity equation. We consider the case of non-smooth Lyapunov function, but barycentrically subdifferentiable only. Sufficient instability conditions are obtained, which are an analogue of the Chetaev theorem and are based on an analysis of the behavior of the non-smooth Lyapunov function in the neighbourhood of the equilibrium. Also we give an example of a dynamical system, the instability of whose equilibrium position is proved using the obtained theorem.

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

In this paper, we consider a boundary value problem for differential equations of Langevin type with the Caputo fractional derivative in a Banach space. It is assumed that the nonlinear part of the equation is a Caratheodory type map. Equations of this type generalize equations of motion in various kinds of media, for example, viscoelastic media or in media where a drag force is expressed using a fractional derivative. We will use the theory of fractional mathematical analysis, the properties of the Mittag-Leffler function, as well as the theory of measures of non-compactness and condensing operators to solve the problem. The initial problem is reduced to the problem of the existence of fixed points of the corresponding resolving integral operator in the space of continuous functions. We will use Sadovskii type fixed point theorem to prove the existence of fixed points of the resolving operator. We will show that the resolving integral operator is condensing with respect to the vector measure of non-compactness in the space of continuous functions and transforms a closed ball in this space into itself.

Пусть n,m, ℓ, s ∈ N – заданные числа, П ⊂ Rⁿ – измеримое по Лебегу множество, X, Z – банаховы идеальные пространства измеримых на П функций. Рассматривается нелинейное операторное уравнение:

x = θ + AF[x], x ∈ X^ℓ, (1)

где A : Z^m → X^ℓ – линейный ограниченный оператор, F : X^ℓ → Z^m – некоторый оператор. Уравнение (1) является естественной формой описания широкого класса сосредоточенных и распределенных систем. Ранее В.П. Политюковым был предложен метод монотонизации для обоснования разрешимости уравнения вида (1) и получения поточечных оценок решения. Суть его состояла в том, что разрешимость уравнения (1) доказывалась (помимо прочих условий) для случая, когда I) оператор F допускал поправку вида G = λI до монотонного оператора F[x] = F[θ+x]+G[x] такую, что II) (I +AG)^−1A > 0 (λ > 0, I  тождественный оператор). Как видно из примеров, приведенных в данной статье, условия I) и II) могут противоречить друг другу, что сужает сферу применения метода. Основной результат статьи в том, что в случае оператора A, обладающего свойством вольтерровости, естественным для эволюционных уравнений, требование монотонизируемости I) можно заменить требованием оценки оператора F на некотором конусном отрезке сверху и снизу через линейный оператор G плюс фиксированный элемент. Доказывается, что для глобальной разрешимости начально-краевой задачи, связанной с полулинейным эволюционным уравнением, достаточно, чтобы аналогичная начально-краевая задача, связанная с линейным уравнением, полученным путем оценки правой части исходного полулинейного уравнения на некотором конусном отрезке, имела положительное решение. В качестве иллюстрации рассматривается применение указанных результатов к системе Гурса–Дарбу, задаче Коши для волнового уравнения и первой краевой задаче для уравнения диффузии.

Let  n,m, ℓ, s ∈ N be given numbers, П ⊂ Rⁿ be a set measurable by Lebesgue and  X, Z  be some Banach ideal spaces of functions measurable on . We consider a nonlinear operator equation of the form as follows:

x = θ + AF[x], x ∈ X^ℓ, (1)

where A : Z^m → X^ℓ is bounded linear operator, F : F : X^ℓ → Z^m is some operator. Equation (1) is a natural form of lumped and distributed parameter systems from a wide enough class. Formerly, by V.P. Polityukov it was suggested monotonization method for justification of solvability of equation (1) and obtaining pointwise estimations for solutions. The matter of this method consisted in that solvability of equation (1) was proved (besides other conditions) under following: I) operator F allows some correction of the form G = λI to monotone operator F[x] = F[θ+x]+G[x] such that II) (I +AG)^−1A > 0 (λ > 0, I is identity operator). As our examples show, conditions I) and II) may be contradictory to each other, that narrows a sphere of application of the method. The main result of the paper is that for the case of operator A, possessing the Volterra property, which is natural for evolutionary equations, the requirement I) of ability to be monotonized can be replaced by the requirement of some upper and lower estimates for operator F on some cone segment through linear operator G and additional fixed element. We prove that for global solvability of a boundary value problem associated with a semilinear evolutionary equation it is sufficient that analogous boundary value problem associated with linear equation, derived from the original equation by estimating of a right-hand side on some cone segment, have a positive solution. The application of results obtained is illustrated by Goursat–Darboux system, Cauchy problem associated with wave equation and first boundary value problem associated with diffusion equation.

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

The axisymmetric model based on the Stokes equations is proposed to investigate the multi-ring pattern formation in two-layer creeping flow with variable thickness of layers. Each layer has uniform density and viscosity. The upper layer is lighter than the lower layer. The flow is generated by both surface and interface geometry. The effect of surface tension is supposed to be negligible. We apply the method of multiple scales to obtain the governing equations describing instability in the form of wave in the flow. Using the Fourier-Laplace method, we analyze the small-amplitude leading behavior of the instability. The asymptotic study reveals that this kind of instability manifests itself as axisymmetric wave which length is comparable with layer thickness; moreover, layer thicknesses play a major role in spatial distribution of wave extrema. The other model parameters alter mostly the wave amplitude. The equation relating extrema distribution to layer thicknesses is derived. We apply the obtained results to study a ring spacing rule observed for most multi-ring basins on the Moon. Using parameters of some lunar multi-ring basins we calculate the consecutive crest radii of the unstable wave and compare the results of simulation with the measured ring radii.

Рассматривается построение и исследование неявных численных схем интегрирования задач динамического контактного взаимодействия двух контактирующих трехмерных тел без трения в рамках альтернирующего метода Шварца. Приводятся результаты тестирования контактного алгоритма декомпозиции Шварца с использованием схемы HTT-$\alpha$ в комбинации с методом перераспределения массы на границе области контакта.

Implicit integration scheme for Schwarz alternating method for dynamic contact interaction problems of two interacting volumetric bodies without friction is considered. The paper presents the results of testing a contact algorithm of Schwarz domain decomposition using HTT-$\alpha$ scheme in consistent method redistribution of mass on the boundary of contact. To prevent artificial oscillations on the contact boundary together with common dissipative properties of the $\alpha$-scheme, the consistent mass redistribution method was used. The main advantage of this approach is the option to use multigrid methods to speed up the algorithm on subdomains, also there is no need for contact elements, contact parameters, Lagrange multipliers or regularization. Numerical examples including various contact zones, different materials of contact bodies and comparisons with measurements of other methods show the wide applicability of the derived algorithm.