Все выпуски

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

The paper is devoted to a method of solving of stationary filtration problems of non-compressible fluid which follows the nonlinear multi-valued anisotropic law of filtration with limiting gradient. This problem mathematically is formulated in the form of variational inequality of the second kind in Hilbert space with inversely strongly monotone operator. The functional occurring in this variational inequality is a sum of several lower semi-continuous convex proper functionals. For solving the considered variational inequality the splitting method is offered.

Представлена детализированная математическая модель радиально-поршневого пневмомотора для анализа динамических процессов и построения механических характеристик. Показан процесс выбора параметров с помощью программы многокритериальной оптимизации MOVI.

The detail mathematical model of radial-piston pneumomotor for dynamic process analysis and calculate mechanical characteristics are presented. The procedure of parameters determine using the multicriterion optimization programme MOVI are shown.

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

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.

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

We study chaotic dynamics in a nonholonomic model of celtic stone. We show that, for certain values of parameters characterizing geometrical and physical properties of the stone, a strange Lorenz-like attractor is observed in the model. We study certain steps of appearance and break-down of this attractor.

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

In many large areas of the economy (such as metallurgy, oil production and refining, electronics), the main competition takes place among several companies that dominate the market. The first models of such markets - oligopolies were described more than a hundred years ago in articles by Cournot, Bertrand, Hotelling. Modeling of oligopolies continues in many modern works. Moreover, in 2014 Nobel Prize in Economics “for his analysis of market power and regulation in sectors with few large companies” was received by Jean Tirole - the author of one of the best modern textbooks on the theory of imperfect competition “The Theory of Industrial Organization”. The main idea of all these publications, studying the behavior of oligopolies, is that every company is primarily concerned with its profits. This approach meets the concept of Nash equilibrium and is actively used in modeling the behavior of players in a competitive market. The exact opposite of such “selfish” equilibrium is “altruistic” concept of Berge equilibrium. In this approach, each player, without having to worry about himself, choose his actions (strategies) trying to maximize the profits of all other market participants. This concept called Berge equilibrium appeared in Russia in 1994 in reference to the France Claude Berge monograph published in 1957. The first works on the concept of Berge equilibrium belong to K.S. Vaisman and V.I. Zhukovskii. Once outside Russia, the concept of “Berge equilibrium” is slowly gaining popularity. To day, the number of publications related to this balance is already measured in tens. However, all of these items are limited to purely theoretical issues, or, in general, to psychology applications. Works devoted to the study of Berge equilibrium in economic problems, were not seen until now. It's probably a consequence of Martin Shubik's review (“… no attention is paid to the application to the economy. … the book is of little interest for economists”) of the Berge's book, it “scared” economists for a long time. However, it is not so simple. In this article, Berge equilibrium is considered in Cournot oligopoly, its relation to Nash equilibrium is studied. Cases are revealed in which players gain more profit by following the concept of Berge equilibrium, than by using strategies dictated by Nash equilibrium.

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

The article investigates properties of the value function of the optimal control problem on infinite horizon with an unlimited integrand index appearing in the quality functional with a discount factor. The estimate is derived for approximating the value function in a problem with the infinite horizon by levels of value functions in problems with lengthening finite horizons. The structure of the value function is identified basing on stationary value functions which depend only on phase variables. The description is given for the asymptotic growth of the value function generated by various types of the quality functional applied in economic and financial modeling: logarithmic, power, exponential, linear functions. The property of continuity is specified for the value function and estimates are deduced for the Hölder parameters of continuity. These estimates are needed for the development of grid algorithms designed for construction of the value function in optimal control problems with infinite horizon.

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

The paper presents the model of anisotropic growth of dendritic crystallization of chemically pure and binary liquid (solution or melt) based on forced convection of the liquid phase. The dependencies of the growth rate and the radius of the top of a dendrite from under-cooling fluid in cases of a chemically pure material and alloys are presented. A comparative analysis of the influence of forced convection on the dendrite growth kinetics is carried out. Evaluation of growth rate and morphology of dendrite by high-speed crystal growth model was done. The contribution of convective flow and the anisotropic properties of the liquid-crystal boundary were taking into account. The model is also used hyperbolic diffusion equation to describe the non-equilibrium impurity capture by crystal surface, which occurs under the rapid crystals growth.

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

A discrete neuron model proposed by Rulkov is studied. In the deterministic version, this system simulates different modes of neural activity, such as quiescence, tonic and chaotic spiking. In the presence of random disturbances, another important mode of bursting characterized by the alternation of quiescence and excitement regimes can be observed. We study the probabilistic mechanisms of noise-induced transitions from quiescence to bursting in the zone of the tangent bifurcation. It is shown that such transitions are accompanied by a transformation of the system dynamics from regular to chaotic. For the analysis of these bifurcation phenomena, the stochastic sensitivity functions technique and method of confidence intervals are used.

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

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.

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

The questions of constructing admissible controls in a problem of optimal control of a nonlinear dynamic system under constraints on its current phase state are discussed. The dynamic system under consideration describes the controlled motion of a carrier rocket from the launching point to the time when the carrier rocket enters a given elliptic earth orbit. The problem consists in designing a program control for the carrier rocket that provides the maximal value of the payload mass led to the given orbit and the fulfillment of a number of additional restrictions on the current phase state of the dynamic system. The additional restrictions are due to the need to take into account the values of the dynamic velocity pressure, the attack and slip angles when the carrier rocket moves in dense layers of the atmosphere. In addition it is required to provide the fall of detachable parts of the rocket into specified regions on the earth surface. For carrier rockets of some classes, such a problem is equivalent to a nonlinear time-optimal problem with phase constraints. Two algorithms for constructing admissible controls ensuring the fulfillment of additional phase constraints are suggested. The numerical analysis of these algorithms is performed. The methodological basis of one algorithm is the application of some predictive control, which is constructed without taking into account the constraints above. Another algorithm is based on special control modes. The results of numerical modeling are presented.