Все выпуски

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

We study the influence of noise on the Goldbeter model of the enzymatic reaction, which describes the mechanism of oscillatory synthesis of cyclic adenosine monophosphate in a cell. It is shown that the model is highly sensitive to variations of parameters and initial conditions. The phenomenon of stochastic excitability in a stable equilibrium zone is demonstrated and studied. We show that the noise results in a sharp transition from low-amplitude stochastic oscillations to large-amplitude spike oscillations. For the parametric analysis of this phenomenon, the technique of stochastic sensitivity functions and the method of confidence ellipses are used. We study how the critical value of the noise intensity corresponding to the generation of large-amplitude oscillations depends on the proximity of a control parameter to a bifurcation point. For a detailed analysis of the frequency properties of stochastic oscillations, a statistical analysis of interspike intervals is carried out, and a phenomenon of coherent resonance is found.

В задаче о движении волчка Ковалевской в двойном поле (случай интегрируемости А.Г. Реймана-М.А. Семенова-Тян-Шанского) вычислен тип всех критических точек отображения момента.

In the problem of motion of the Kowalevski top on two constant fields (the A.G. Reyman-M.A. Semenov-Tian-Shansky case) the type of all critical points of the momentum map is calculated.

Хорошо известно, что преобразование ренормгруппы $\mathcal{R}$ имеет единственную неподвижную точку $f_ {cr}$ в пространстве критических $C^{3}$-гомеоморфизмов окружности с одной кубической критической точкой $x_{cr}$ и числом вращения равным золотому сечению $\overline{\rho}: =\frac{\sqrt{5} -1}{2}.$ Обозначим через $Cr(\overline{\rho})$ множество всех критических отображений окружности $C^ {1}$-сопряженных к $f_{cr}.$ Пусть $f\in Cr(\overline{\rho})$ и $\mu:=\mu_{f}$ --- единственная вероятностная инвариантная мера для $f.$ Зафиксируем $\theta \in (0,1).$ Для каждого $n\geq 1$ определим $c_{n}:=c_{n}(\theta)$ такое, что $\mu([x_{cr}, c_{n}]) = \theta\cdot\mu([x_{cr}, f^{q_{n}} (x_{cr})]),$ где $q_{n}$ --- время первого возврата линейного вращения $f_{\overline{\rho}}.$ Мы исследуем закон сходимости перемасштабированного точечного процесса времени попадания. Мы показываем, что предельное распределение сингулярно относительно меры Лебега.

It is well known that the renormalization group transformation $\mathcal{R}$ has a unique fixed point $f_{cr}$ in the space of critical $C^{3}$-circle homeomorphisms with one cubic critical point $x_{cr}$ and the golden mean rotation number $\overline{\rho}:=\frac{\sqrt{5}-1}{2}.$ Denote by $Cr(\overline{\rho})$ the set of all critical circle maps $C^{1}$-conjugated to $f_{cr}.$ Let $f\in Cr(\overline{\rho})$ and let $\mu:=\mu_{f}$ be the unique probability invariant measure of $f.$ Fix $\theta \in(0,1).$ For each $n\geq1$ define $c_{n}:=c_{n}(\theta)$ such that $\mu([x_{cr},c_{n}])=\theta\cdot\mu([x_{cr},f^{q_{n}}(x_{cr})]),$ where $q_{n}$ is the first return time of the linear rotation $f_{\overline{\rho}}.$ We study convergence in law of rescaled point process of time hitting. We show that the limit distribution is singular w.r.t. the Lebesgue measure.

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

We consider a boundary-value problem for the nonlinear evolution partial differential equation, given in renormalized form. This problem appears in rotary system mechanics and describes the transverse vibrations of the rotating rotor of a constant cross-section from a viscoelastic material whose ends are pivotally fixed. The question of the stability of the zero equilibrium state is studied, the critical value of the rotor speed is found, above which continuous oscillations occur. Exact solutions of the boundary-value problem are found in the form of single-mode functions with respect to the spatial variable and functions periodic in time. The stability conditions for such solutions are derived, and in some cases an analysis of the stability conditions is given. The paper shows the absence of multimode time-periodic solutions. The basic and important (from an applied point of view) particular cases of this nonlinear boundary-value problem are analyzed. All the results of the analysis of a nonlinear boundary-value problem are analytical. Their conclusion is based on the qualitative theory of infinite-dimensional dynamical systems.

В настоящей работе проведено исследование модели деформаций системы из $n$ стилтьесовских струн, расположенных вдоль геометрического графа-звезды, с нелинейным условием в узле. Соответствующая граничная задача имеет вид $$ \left\{\begin{array}{lll} -\left(p_iu_i^\prime\right)(x)+\displaystyle{\int_{0}^{x}}u_i\,dQ_i=F_i(x)-F_i(+0)-(p_iu_i')(+0),\quad i=1,2, \ldots, n,\\ \sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0),\\ u_1(0)=u_2(0)=\ldots=u_n(0)=u(0),\\ (p_iu_i')(l_i-0)+u_i(l_i)\Delta Q_i(l_i)=\Delta F_i(l_i),\quad i=1,2,\ldots, n. \end{array} \right. $$ Здесь функции $u_i(x)$ определяют деформации каждой из струн; $F_i(x)$ описывают распределение внешней нагрузки; $p_i(x)$ характеризуют упругость струн; $Q_i(x)$ описывают упругую реакцию внешней среды. Скачок $\Delta F_i(l_i)$ равняется сосредоточенной в точке $l_i$ внешней силе; скачок $\Delta Q_i(l_i)$ совпадает с жесткостью упругой опоры (пружины), прикрепленной к точке $l_i$. Условие $\sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0)$ возникает за счет наличия в узле ограничителя, представленного отрезком $[-m,m]$, на перемещение струн под воздействием внешней нагрузки, то есть предполагается, что $|u(0)|\leq m$. Здесь через $N_{[-m,m]}u(0)$ обозначен нормальный конус к $[-m,m]$ в точке $u(0)$. В работе проведен вариационный вывод модели; доказаны теоремы существования и единственности решения; проанализированы критические нагрузки, при которых происходит соприкосновение струн с ограничителем; приведена явная формула представления решения.

In the present paper we study a model of deformations for a system of $n$ Stieltjes strings located along a geometric graph-star with a nonlinear condition at the node. The corresponding boundary value problem has the form $$ \left\{\begin{array}{lll} -\left(p_iu_i^\prime\right)(x)+\displaystyle{\int_0^x}u_idQ_i=F_i(x)-F_i(+0)-(p_iu_i')(+0), \quad i=1,2, \ldots, n,\\ \sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0),\\ u_1(0)=u_2(0)=\ldots=u_n(0)=u(0),\\ (p_iu_i')(l_i-0)+u_i(l_i)\Delta Q_i(l_i)=\Delta F_i(l_i), \quad i=1,2,\ldots, n. \end{array} \right. $$ Here the functions $u_i(x)$ determine the deformations of each of the strings; $F_i(x)$ describe the distribution of the external load; $p_i(x)$ characterize the elasticity of strings; $Q_i(x)$ describe the elastic response of the environment. The jump $\Delta F_i(l_i)$ is equal to the external force concentrated at the point $l_i$; the jump $\Delta Q_i(l_i)$ coincides with the stiffness of the elastic support (spring) attached to the point $l_i$. The condition $\sum\limits_{i=1}^np_i(+0)u_i'(+0)\in N_{[-m,m]}u(0)$ arises due to the presence of a limiter in the node represented by the segment $ [-m,m]$, on the movement of strings under the influence of an external load, thus it is assumed that $|u(0)|\leq m$. Here $N_{[-m,m]}u(0)$ denotes the normal cone to $[-m,m]$ at the point $u(0)$. In the present paper a variational derivation of the model is carried out; existence and uniqueness theorems for solutions are proved; the critical loads at which the strings come into contact with the limiter are analyzed; an explicit formula for the representation of the solution is presented.

Представлена полная аналитическая классификация атомов гиростата Ковалевской–Яхья, возникающих в критических точках ранга 1. Найдены все разделяющие значения гиростатического момента при классификации диаграмм Смейла–Фоменко. Разработан "конструктор" графов Фоменко, применение которого дало полное описание грубой топологии этого интегрируемого случая. Доказано, что имеется девять групп эквивалентных молекул (без меток), содержащих 22 устойчивых графа и 6 неустойчивых по отношению к количеству критических окружностей на критических уровнях.

We present the complete analytical classification of the atoms arising at the critical points of rank 1 of the Kowalevski–Yehia gyrostat. To classify the Smale–Fomenko diagrams, all separating values of the gyrostatic momentum are found. We present a kind of constructor of the Fomenko graphs; its application gives the complete description of the rough topology of this integrable case. It is proved that there exists exactly nine groups of identical molecules (not considering the marks). These groups contain 22 stable types of graphs and 6 unstable ones with respect to the number of critical circles on the critical levels.

В работе рассмотрена интегрируемая гамильтонова система на алгебре Ли $so(4)$ с дополнительным интегралом четвертой степени - интегрируемый случай Адлера-ван Мёрбеке. Рассмотрены классические работы, посвященные, с одной стороны, динамике твердого тела, содержащего полости, полностью заполненные идеальной жидкостью, совершающей однородное вихревое движение, а с другой стороны, изучению геодезических потоков левоинвариантных метрик на группах Ли. Приведены уравнения движения, функция Гамильтона, скобки Ли-Пуассона, функции Казимира и фазовое пространство рассматриваемого случая. В предыдущих работах начато исследование фазовой топологии интегрируемого случая Адлера-ван Мёрбеке: приводятся в явном виде спектральная кривая, дискриминантное множество, бифуркационная диаграмма отображения момента, предъявлены характеристические показатели для определения типа критических точек ранга 0 и 1 отображения момента. В данной работе излагается алгоритм построения торов Лиувилля. Рассмотрены примеры перестроек лиувиллиевых торов при пересечении бифуркационных кривых для перестроек одного тора в два и двух торов в два.

In this paper we consider an integrable Hamiltonian system on the Lie algebra $so(4)$ with an additional integral of the fourth degree - the Adler-van Moerbeke integrable case. We discuss classical works which explore, on the one hand, the dynamics of a rigid body with cavities completely filled with an ideal fluid performing a homogeneous vortex motion and, on the other hand, are devoted to the study of geodesic flows of left-invariant metrics on Lie groups. The equations of motion, the Hamiltonian function, Lie-Poisson brackets, Casimir functions and the phase space of the case under consideration are given. In previous papers, the investigation of the phase topology of the integrable Adler-van Moerbeke case was started: a spectral curve, a discriminant set and a bifurcation diagram of the moment map are explicitly shown, and characteristic exponents for determining the type of critical points of rank 0 and 1 of the moment map are presented. In this paper we present an algorithm for constructing Liouville tori. Examples are given of bifurcations of Liouville tori at the intersection of bifurcation curves for reconstructions of one torus into two tori and of two tori into two tori.

Рассмотрены трехмерные задачи узлов для простой кубической решетки и твердых сфер, находящихся в хаотическом движении. Установлены дополнительные (к двухпоказательному скейлингу) соотношения между индексами: 2-α-γ=ν (или νd-γ=ν) и β=-2α. Определены численные значения трехмерных критических индексов: α=-2/11, η=0, β=4/11, ν=8/11, γ=16/11 и δ=5.

Three-dimensional lattice points problems for simple cubic lattice and solid spheres in chaotic motion are considered. Additional (to two-exponential scaling) relations between indices are indicated: 2-α-γ=ν (or νd-γ=ν) and β=-2α. Numerical values of three-dimensional critical indices are defined: α=-2/11, η=0, β=4/11, ν=8/11, γ=16/11 and δ=5.

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

We study the long-time evolution of axisymmetric free-surface two-layered creeping flow subject to the initial topography of its boundaries and bottom velocities. Each layer has uniform density and viscosity. The upper layer is assumed to have a smaller density than the lower layer. Based on lubrication approximation (the Reynolds equations) the nonlinear system of diffusion-type equations with respect to the surface and interface between the layers is obtained to describe this flow. Taking the dimensionless density contrast between the layers as a small parameter, we apply the method of asymptotic expansions to extract leading-term approximation for the slowly varying large-time evolution of the governing equations. An asymptotic equation relating both surface and interface displacement to the bottom velocities is derived. Based on this equation, we develop the algorithm to calculate velocity fields within layers for large time. Streamlines are used to visualize the flow. Numerical results reveal stability of the streamlines in the upper layer under variation of the bottom velocity. As geophysical applications, the developed algorithm is used to evaluate the velocity field in the crust (the upper layer) beneath the large-scale lunar multi-ring basins influenced by deep movements in the underlying mantle (the lower layer). To validate the results of modeling, we compare the calculated velocity fields with basin ridge systems obtained by experimental observations. The model comparison has shown proximity of radii of basin rings and critical points of the surface velocity.

Классическая система реакции-диффузии — система Шнакенберга — рассматривается в ограниченной области $m$-мерного пространства, на границе которой предполагаются выполненными краевые условия Неймана. Изучается диффузионная неустойчивость стационарного пространственно-однородного решения этой системы, называемая также неустойчивостью Тьюринга, возникающая при изменении коэффициента диффузии $d.$ Путем анализа линеаризованной системы в бездиффузионном и диффузионном приближениях получено аналитическое описание области необходимых и достаточных условий неустойчивости Тьюринга на плоскости параметров системы. Показано, что одна из границ области необходимых условий является огибающей семейства кривых, ограничивающих область достаточных условий. При этом точки пересечения двух соседних кривых лежат на прямой, угловой коэффициент которой зависит от собственных значений оператора Лапласа в рассматриваемой области и не зависит от коэффициента диффузии. Найдено аналитическое выражение критического коэффициента диффузии, при котором происходит потеря устойчивости положения равновесия системы. Указаны условия, в зависимости от которых множество волновых чисел, соответствующих нейтральным модам устойчивости, счетно, конечно или пусто. Показано, что полуось $d>1$ можно представить в виде счетного объединения полуинтервалов, каждому из которых соответствует минимальное волновое число, при котором происходит потеря устойчивости, причем точки разбиения полуоси выражаются через собственные значения оператора Лапласа в рассматриваемой области.

A classical reaction-diffusion system, the Schnakenberg system, is under consideration in a bounded domain $\Omega\subset\mathbb{R}^m$ with Neumann boundary conditions. We study diffusion-driven instability of a stationary spatially homogeneous solution of this system, also called the Turing instability, which arises when the diffusion coefficient $d$ changes. An analytical description of the region of necessary and sufficient conditions for the Turing instability in the parameter plane is obtained by analyzing the linearized system in diffusionless and diffusion approximations. It is shown that one of the boundaries of the region of necessary conditions is an envelope of the family of curves that bound the region of sufficient conditions. Moreover, the intersection points of two consecutive curves of this family lie on a straight line whose slope depends on the eigenvalues of the Laplace operator and does not depend on the diffusion coefficient. We find an analytical expression for the critical diffusion coefficient at which the stability of the equilibrium position of the system is lost. We derive conditions under which the set of wavenumbers corresponding to neutral stability modes is countable, finite, or empty. It is shown that the semiaxis $d>1$ can be represented as a countable union of half-intervals with split points expressed in terms of the eigenvalues of the Laplace operator; each half-interval is characterized by the minimum wavenumber of loss of stability.