Все выпуски

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

We consider a problem with data on the characteristics for a loaded system of hyperbolic equations of the second order on a rectangular domain. The questions of the existence and uniqueness of the classical solution of the considered problem, as well as the continuity dependence of the solution on the initial data, are investigated. We propose a new approach to solving the problem with data on the characteristics for the loaded system of hyperbolic equations second order based on the introduction new functions. By introducing new unknown functions the problem is reduced to an equivalent family of Cauchy problems for a loaded system of differential with a parameters and integral relations. An algorithm for finding an approximate solution to the equivalent problem is proposed and its convergence is proved. Conditions for the unique solvability of the problem with data on the characteristics for the loaded system of hyperbolic equations of the second order are established in the terms of coefficient's system.

Рассматривается система реакции-диффузии с кубической нелинейностью, которая является бесконечномерным аналогом классической системы Рэлея и частным случаем системы Фитцью-Нагумо. Предполагается, что пространственная переменная изменяется на отрезке, на концах которого заданы однородные краевые условия Неймана. Известно, что в данном случае в системе Рэлея с диффузией существует пространственно-однородный автоколебательный режим, совпадающий с предельным циклом классической системы Рэлея. В настоящей работе показано существование счетного множества критических значений управляющего параметра, при которых возникают пространственно-неоднородные автоколебательные и стационарные режимы. Данные режимы устойчивы относительно возмущений, принадлежащих некоторым бесконечномерным инвариантным подпространствам системы, но неустойчивы во всем фазовом пространстве. Это свойство объясняет, почему в результате численных экспериментов при некоторых значениях параметра различным начальным условиям соответствуют нулевое, периодическое по времени или стационарное решение. Асимптотика вторичных решений построена методом Ляпунова-Шмидта. Явно найдены первые члены разложения, проанализированы формулы для общего члена асимптотики. Показано, что на инвариантных подпространствах происходит мягкая потеря устойчивости нулевого равновесия. Эволюция вторичных режимов при увеличении значений надкритичности исследована численно. Установлено, что с ростом значений надкритичности вторичные автоколебательные режимы постепенно сменяются стационарными. Амплитуда стационарных решений растет по мере увеличения надкритичности, а профиль асимптотически стремится к профилю меандра.

We consider a reaction-diffusion system with a cubic nonlinear term, which is a special case of the Fitzhugh-Nagumo system and an infinite-dimensional version of the classical Rayleigh system. We assume that the spatial variable belongs to an interval, supplemented with Neumann boundary conditions. It is well-known that in that specific case there exists a spatially-homogeneous oscillatory regime, which coincides with the time-periodic solution of the classical Rayleigh system. We show that there exists a countable set of critical values of the control parameter, where each critical value corresponds to the branching of new spatially-inhomogeneous auto-oscillatory or stationary regimes. These regimes are stable with respect to small perturbations from some infinite-dimensional invariant subspaces of the system under study. This, in particular, explains the convergence of numerical solution to zero, periodic or stationary solution, which is observed for some specific initial conditions and control parameter values. We construct the asymptotics for branching solutions by using Lyapunov-Schmidt reduction. We find explicitly the first terms of asymptotic expansions and study the formulas for general terms of asymptotics. It is shown that a soft loss of stability occurs in invariant subspaces. We study numerically the evolution of secondary regimes due to the increase of control parameter values and observe that the secondary periodic solutions are transformed into stationary ones as the control parameter value increases. Next, the amplitude of stationary solutions continues to grow and the solution asymptotically converges to the square wave regime.

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

We study an initial-boundary value problem for a multidimensional pseudoparabolic equation with variable coefficients and boundary conditions of the third kind. The multidimensional pseudoparabolic equation is reduced to an integro-differential equation with a small parameter. It is shown that as the small parameter tends to zero, the solution of the resulting modified problem converges to the solution of the original problem. For an approximate solution of the obtained problem, a locally one-dimensional difference scheme by A. A. Samarsky is constructed. An a priori estimate is obtained by the method of energy inequalities, from which the uniqueness, stability, and convergence of the solution of the locally one-dimensional difference scheme to the solution of the original differential problem follow. For a two-dimensional problem, an algorithm for the numerical solution of the initial-boundary value problem for a pseudoparabolic equation with conditions of the third kind is developed.

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

We study the equilibrium problem of a transversely isotropic plate with rigid inclusions. It is assumed that the plate deforms under hypotheses of classical elasticity. The problems are formulated as the minimization of the plate energy functional on the convex and closed subset of the Sobolev space. It is established that, as the geometric parameter (the size) of the volume inclusion tends to zero, the solutions converge to the solution of an equilibrium problem of a plate with a thin rigid inclusion. Also the case of the delamination of an inclusion is investigated when a crack in the plate is located along one of the inclusion edges. In the problem of a plate with a delaminated inclusion the nonlinear condition of nonpenetration is given. This condition takes the form of a Signorini-type inequality and describes the mutual nonpenetration of the crack edges. For the problem with a delaminated inclusion, the equivalence of variational and differential statements is proved provided a sufficiently smooth solution. For each considered variation problem, unique solvability is established.

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

The work is devoted to the study of the second initial-boundary value problem for a general-form third-order differential equation of pseudoparabolic type with variable coefficients in a multidimensional domain with an arbitrary boundary. In this paper, a multidimensional pseudoparabolic equation is reduced to an integro-differential equation with a small parameter, and a locally one-dimensional difference scheme by A.A. Samarskii is used. Using the maximum principle, an a priori estimate is obtained for the solution of a locally one-dimensional difference scheme in the uniform metric in the $C$ norm. The stability and convergence of the locally one-dimensional difference scheme are proved.

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

The paper is devoted to the analysis of dynamical bimatrix games with integral indices discounted on an infinite time interval. The system dynamics is described by differential equations in which players' behavior changes according to incoming control signals. For this game, a problem of construction of equilibrium trajectories is considered in the framework of minimax approach proposed by N.N. Krasovskii and A.I. Subbotin in the differential games theory. The game solution is based on the structure of dynamical Nash equilibrium developed in papers by A.F. Kleimenov. The maximum principle of L.S. Pontryagin in combination with the method of characteristics for Hamilton-Jacobi equations are applied for the synthesis of optimal control strategies. These methods provide analytical formulas for switching curves of optimal control strategies. The sensitivity analysis for equilibrium solutions is implemented with respect to the discount parameter in the integral payoff functional. It is shown that equilibrium trajectories in the problem with the discounted payoff functional asymptotically converge to the solution of a dynamical bimatrix game with average integral payoff functionals examined in papers by V.I. Arnold. Obtained results are applied to a dynamical model of investments on financial markets.

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

In the article a numerical technique based on the finite difference method is proposed for the approximate solution of a second order nonlocal boundary value problem for ordinary differential equations. It is clear that a bridge designed with two support points at each end point leads to a standard two-point local boundary value condition, and a bridge contrived with multi-point supports corresponds to a multi-point boundary value condition. At the same time if non-local boundary conditions can be set up near each endpoint of a multi-point support bridge, a two-point nonlocal boundary condition arises. The computational results for the nonlinear model problem are presented to validate the proposed idea. The effect of parameters variation on the convergence of the proposed method is analyzed.

Хорошо известно, что методы сопряженного градиента полезны при решении масштабных задач нелинейной оптимизации без ограничений. В данной работе мы рассматриваем объединение лучших свойств двух методов сопряженного градиента. В частности, мы даем новый метод сопряженного градиента, основанный на гибридизации полезных методов DY (Dai-Yuan) и HZ (Hager-Zhang). Параметры гибрида выбираются таким образом, чтобы предложенный метод удовлетворял условиям сопряженности и достаточного спуска. Показано, что новый метод сохраняет свойство глобальной сходимости двух вышеупомянутых методов. Описаны численные результаты для набора стандартных тестовых задач. Показано, что в большинстве случаев эффективность предложенного метода выше, чем у DY и HZ.

It is well known that conjugate gradient methods are useful for solving large-scale unconstrained nonlinear optimization problems. In this paper, we consider combining the best features of two conjugate gradient methods. In particular, we give a new conjugate gradient method, based on the hybridization of the useful DY (Dai-Yuan), and HZ (Hager-Zhang) methods. The hybrid parameters are chosen such that the proposed method satisfies the conjugacy and sufficient descent conditions. It is shown that the new method maintains the global convergence property of the above two methods. The numerical results are described for a set of standard test problems. It is shown that the performance of the proposed method is better than that of the DY and HZ methods in most cases.