Все выпуски

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

The influence of additive and parametrical noise on attractors of the one-dimensional system governed by the stochastic differential Ito equation is investigated. It is shown that unlike additive, parametrical disturbances lead to the shift of extrema of probability density function. For the value of this shift, a decomposition on small parameter of noise intensity is obtained. It is shown that the influence of the parametrical noise can change not only the arrangement, but also the quantity of extrema of probability density function. The corresponding noise-induced phenomena are studied for three dynamical models in detail. An analysis of the error for the different order estimations of the shift of extrema for the probability density function is presented by the example of a linear model. Two scenarios of the transition between unimodal and bimodal forms of the stochastic attractor are investigated for systems with different types of cubic nonlinearity.

Рассмотрено взаимное расположение вершин единичного многомерного куба, аффинного подпространства и его ортогональных проекций на координатные подпространства. Даны верхние и нижние ограничения размерности подпространства, при которых некоторая ортогональная проекция всегда сохраняет отношение инцидентности подпространства и вершин куба. Также рассмотрены некоторые косоугольные проекции. Кроме того, дан краткий обзор истории развития многомерной начертательной геометрии. Аналитические и синтетические методы в геометрии обособились с XVII века. Хотя анализ и синтез тесно переплетаются, с этого времени многие геометры и инженеры делают тонкое различие. Указания на идею о многомерном пространстве можно найти в работах XVIII века, но настоящее развитие началось с середины XIX века. Вскоре такие работы появились и на русском языке. Далее многие математики обобщали свои теории на многомерный случай. Наши новые результаты получены аналитическими и синтетическими методами. Они иллюстрируют сложность задач псевдобулева программирования, поскольку снижение размерности задачи методом ортогонального проектирования встречает препятствие в худшем случае.

We consider the arrangement of vertices of a unit multidimensional cube, an affine subspace, and its orthogonal projections onto coordinate subspaces. Upper and lower bounds on the subspace dimension are given under which some orthogonal projection always preserves the incidence relation between the subspace and cube vertices. Some oblique projections are also considered. Moreover, a brief review of the history of the development of multidimensional descriptive geometry is given. Analytic and synthetic methods in geometry diverged since the 17th century. Although both synthesis and analysis are tangled, from this time forth many geometers as well as engineers keep up a nice distinction. One can find references to the idea of higher-dimensional spaces in the 18th-century works, but proper development has been since the middle of the 19th century. Soon such works have appeared in Russian. Next, mathematicians generalized their theories to many dimensions. Our new results are obtained by both analytic and synthetic methods. They illustrate the complexity of pseudo-Boolean programming problems because reducing the problem dimension by orthogonal projection meets obstacles in the worst case.

Изучаются свойства дискретной вариационной задачи динамической аппроксимации в комплексном евклидовом (L + 1)-мерном пространстве E. Она обобщает известные задачи среднеквадратической полиномиальной аппроксимации функций, заданных своими отсчетами в конечном интервале. В рассматриваемой задаче аппроксимация последовательности y = {y_i}_L⁰ отсчетов функции y(t) ∈ L²[0, T], T = Lh на сетке I_h осуществляется решениями однородных линейных дифференциальных или разностных уравнений заданного порядка n с постоянными, но, возможно, неизвестными коэффициентами. Тем самым показано, что в последнем случае задача аппроксимации включает в себя и задачу идентификации. Анализ ее особенностей - основная тема статьи. Ставится задача нахождения вектора коэффициентов  разностного уравнения Σ_n⁰ ŷ_i+k α_i = 0, где k = 0,L − n. Оптимизируются коэффициенты  и начальные условия переходного процесса y этого уравнения. Цель оптимизации - наилучшая аппроксимация исследуемого динамического процесса y ∈ E. Критерий аппроксимации  минимум величины ||y − ŷ||²_E. Показано, что изучаемая вариационная задача сводится к задачам проектирования в E вектора y на ядра разностных операторов с неизвестными коэффициентами α ∈ ω ⊂ S ⊂ Eⁿ⁺¹. Здесь α - направление, S - сфера или гиперплоскость. Показана связь изучаемой задачи с задачами дискретизации и идентифицируемости. Тогда координаты вектора y ∈ E есть точное решение дифференциального уравнения на сетке I_h и y = ŷ. Дано сравнение изучаемой задачи вариационной идентификации с алгебраическими методами идентификации. Показано, что ортогональные дополнения к ядрам разностных операторов всегда имеют теплицев базис. Это приводит к быстрым проекционным алгоритмам вычислений. Показано, что задача нахождения оптимального вектора α сводится к задаче безусловной минимизации функционала идентификации, зависящего от направления  в Eⁿ⁺¹. Предложена итерационная процедура его минимизации на сфере с широкой областью и высокой скоростью сходимости. Изучаемую вариационную задачу можно применять при математическом моделировании в управлении и научных исследованиях. При этом на конечных интервалах может использоваться, в частности, возможность кусочно-линейной динамической аппроксимации сложных динамических процессов разностными и дифференциальными уравнениями указанного типа.

Some properties of the discrete variational problem of the dynamic approximation in the complex Euclidean (L + 1)-dimensional space are studied here. It generalizes familiar problems of the mean square polynomial approximation of the functions given on the finite interval in accordance with their references. In the problem under consideration sequence approximation y = {y_i}_L⁰ of the references of the function y(t) ∈ L²[0, T], T = Lh on the lattice I_h is achieved by solving homogeneous linear differential equations or difference equations of the given order n with constant but possibly unknown coefficients. Thus, it is shown that in the latter case the approximation problem also includes the identification problem. The analysis of its properties is the main subject of the article. The problem is set to find vector of coefficients  of difference equation Σ_n⁰ ŷ_i+k α_i = 0, where k = 0,L − n. Coefficients  and initial conditions of the transient process by of this equation are optimized. The optimization purpose is to achieve the best approximation of the dynamic process y ∈ E being considered here. The approximation criterion is a minimum of the quantity ||y − ŷ||²_E. The variational problem under study is shown to be reduced to the problem of projecting vector y in E on the kernels of the difference operators with unknown coefficients  α ∈ ω ⊂ S ⊂ Eⁿ⁺¹, where  is a direction, S is a sphere or a hyperplane. The problem under study is shown to be related to the problems of the discretization and identifiability. In this case vector coordinates y ∈ E is an exact solution of differential equation on the lattice I_h and y = ŷ. The problem of the variational identification is compared with algebraic methods of identification. The orthogonal complement to the kernels of the difference operators are shown to always have Toeplitz basis. This results in fast projecting algorithms of computation. The problem of finding optimal vector α is shown to be reduced to the problem of the absolute minimization of the identification functional depending on the direction  in Eⁿ⁺¹. The iterative procedure of its minimization on a sphere with wide domain and high speed of convergence is presented here. The variational problem considered here can be applied in mathematical modeling for control problem and research purposes. On the finite intervals, for example, it is possible to use piecewise-linear dynamic approximations of the complex dynamic processes with difference and differential equations of the specified type.

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

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.

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

One of the versions of the generalized variational Ginzburg-Landau equation is considered, supplemented by periodic boundary conditions. For such a boundary value problem, the question of existence, stability, and local bifurcations of single-mode equilibrium states is studied. It is shown that in the case of a nearly critical threefold zero eigenvalue, in the problem of stability of single-mode spatially inhomogeneous equilibrium states, subcritical bifurcations of two-dimensional invariant tori filled with spatially inhomogeneous equilibrium states are realized. The analysis of the stated problem is based on such methods of the theory of infinite-dimensional dynamical systems as the theory of invariant manifolds and the apparatus of normal forms. Asymptotic formulas are obtained for the solutions that form invariant tori.

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

The paper considers a model of chemical kinetics for which the derivation of equations does not rely on the law of mass action, but is rather based on such principles as geometric probability and joint probability. For this model a generalization is constructed for the case of reaction-diffusion systems in heterogeneous medium, with respect to the convective and diffusive transfer of heat. The construction of this generalization is carried out by an alternative methodology, which is based fully on systems of ordinary differential equations, without a transition to partial derivatives. The description of this new method is a bit similar to the finite volume method, except that it uses statistical simplifying positions and geometric probability to describe diffusion processes. Such approach allows us to greatly simplify the numerical implementation of the resulting model, as well as to simplify its quantitative analysis by dynamical systems theory methods. Moreover, the efficiency of parallel implementation of the numerical method is increased for the resulting model. In addition, the author considers an application of this model for the description of some example reaction with quasi-periodic regime, as well as an algorithm for the transition from standard models with dimensional kinetic constants to its formalism.

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

The problem of stationary flame propagation is studied by the variational method based on the non-equilibrium thermodynamic approach. The analysis has been carried out for the case of one-dimensional propagation of premixed gaseous flame. The primary aim of the analysis is formulation of the variational principle for minimization of the functional at the stationary state of the system, which is identified as a steady-state regime of flame propagation. The possibility to obtain a relationship for the prediction of flame propagation velocity as a dependent variable has been investigated through the different representations of potential.

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

Mathematical simulation of convective-conductive-radiative heat transfer in a cubical cavity filled with diathermanous medium has been carried out. The domain of interest is a closed volume having heat-conducting solid walls of finite thickness with diffuse grey inner surfaces. The outer surfaces of two vertical walls are isothermal while the other walls are adiabatic. The boundary-value problem has been formulated in dimensionless variables such as “vector potential-vorticity vector-temperature’’ in the Boussinesq approximation and taking into account the diathermancy of the continuous medium. An analysis of surface thermal radiation has been conducted on the basis of the net-radiation method in the form of Poljak. The formulated transient boundary-value problem has been solved by finite difference method in a wide range of the Rayleigh number, thermal conductivity ratio and surface emissivity. Correlations for the average convective and radiative Nusselt numbers at the characteristic internal solid-fluid interface have been obtained. The comparison between the obtained three-dimensional results and the two-dimensional data has been conducted. It has been found, that on the basis of a three-dimensional model it is possible to analyze the formation of intensive transverse flows from two vertical surfaces which are absent in a two-dimensional model. It has been also shown, that the solution of convective-radiative heat transfer problems in the conjugate statement leads to essential changes in distributions of local and integral parameters in comparison with the non-conjugate model, which first of all is related to a more correct description of the thermal radiation in diathermanous media due to taking into account the thermal conduction of the solid walls.

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

The paper studies the motions of a dynamically symmetric satellite (rigid body) relative to the center of mass in the central Newtonian gravitational field on a weakly elliptical orbit in the neighborhood of its stationary rotation (cylindrical precession). We consider the values of the parameters for which, in the limiting case of a circular orbit, one of the frequencies of small linear oscillations is equal to unity and the other is equal to zero, and the rank of the coefficient matrix of the linearized equations of the perturbed motion is equal to two, as well as a small neighborhood of this resonant point in the three-dimensional space of parameters. The resonant periodic motions of the satellite, analytical in fractional powers of a small parameter (the eccentricity of the orbit of the satellite's center of mass), are constructed. A rigorous nonlinear analysis of their stability is carried out. The methods of KAM theory are used to describe two- and three-frequency conditionally periodic motions of a satellite, with frequencies of different orders in a small parameter. A number of general theoretical issues concerning the considered multiple parametric resonance in Hamiltonian systems with two degrees of freedom that are close to autonomous and periodic in time are discussed. Several qualitatively different variants of parametric resonance regions are constructed. It is shown that in the general case the nature of nonlinear resonant oscillations of the system is determined by the first approximation system in a small parameter.

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

The concept of $\alpha$-set in a finite-dimensional Euclidean space, which is one of generalizations of the notion of a convex set, is introduced. The emergence of this concept is connected with the study of properties of attainability sets of nonlinear controlled systems. The numerical characteristic of nonconvexity degree of a set on the basis of which a classification of sets is carried out is defined in the paper. Analogs of basic concepts from the convex analysis are introduced into consideration and their properties are studied. Statements in the spirit of such theorems from the convex analysis as the theorem of existence of basic hyperplane to a convex set and theorems of separability of convex sets in Euclidean space are formulated and proved. The concept of magoriums of nonconvex sets is studied. Property of a magoriums is a sufficient condition for representation of a closed nonconvex set in the form of crossing of half-spaces in the sense of definitions entered in this work. The obtained results of the theory of separability of nonconvex sets can be extended on a case of hypograph and epigraph of the scalar functions with Lipschitz condition.