Все выпуски

Рассматривается задача о назначении спектра показателей Ляпунова линейной управляемой системы с дискретным временем $$x(m+1)=A(m)x(m)+B(m)u(m),\quad m\in\mathbb N,\ x\in\mathbb R^{n},\ u\in\mathbb R^{k}, \qquad (1)$$ посредством линейной по фазовым переменным обратной связи $u(m)=U(m)x(m)$ в малой окрестности спектра показателей свободной системы $$x(m+1)=A(m)x(m),\quad m\in\mathbb N,\ x\in\mathbb R^{n}. \qquad (2)$$ Дополнительно требуется, чтобы норма матрицы обратной связи $U(\cdot)$ удовлетворяла липшицевой оценке по отношению к требуемому смещению показателей. Это свойство называется пропорциональной локальной управляемостью полного спектра показателей Ляпунова замкнутой системы $$x(m+1)=\bigl(A(m)+B(m)U(m)\bigr)x(m),\quad m\in\mathbb N,\ x\in\mathbb R^{n}. \qquad (3)$$ Построен пример, показывающий, что найденные ранее достаточные условия пропорциональной локальной управляемости полного спектра показателей Ляпунова системы (3) (равномерная полная управляемость системы (1) и устойчивость показателей Ляпунова свободной системы (2)) не являются необходимыми.

We consider a problem of assigning the Lyapunov spectrum for a linear control discrete-time system $$x(m+1)=A(m)x(m)+B(m)u(m),\quad m\in\mathbb N,\ x\in\mathbb R^{n},\ u\in\mathbb R^{k}, \qquad (1)$$ in a small neighborhood of the Lyapunov spectrum of the free system $$x(m+1)=A(m)x(m),\quad m\in\mathbb N,\ x\in\mathbb R^{n},\qquad (2) $$ by means of linear feedback $u(m)=U(m)x(m)$. We assume that the norm of the feedback matrix $U(\cdot)$ satisfies the Lipschitz estimate with respect to the required shift of the Lyapunov spectrum. This property is called proportional local assignability of the Lyapunov spectrum of the closed-loop system $$x(m+1)=\bigl(A(m)+B(m)U(m)\bigr)x(m),\quad m\in\mathbb N,\ x\in\mathbb R^{n}. \qquad (3)$$ We previously proved that uniform complete controllability of system (1) and stability of the Lyapunov spectrum of free system (2) are sufficient conditions for proportional local assignability of the Lyapunov spectrum of closed-loop system (3). In this paper we give an example demonstrating that these conditions are not necessary.

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

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.

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

The issues of uniqueness and stability of solutions to the control problems for the model of electron-induced charging of an inhomogeneous polar dielectric are studied. Sufficient conditions for the uniqueness and stability of optimal solutions to the considered extremum problems are established, and the local estimates of their stability with respect to small perturbations of the cost functionals are derived.

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

В работе рассмотрен локально-неравновесный процесс затвердевания переохлажденного бинарного расплава. В целях простоты предполагается, что затвердевающая бинарная система находится при постоянных температуре и давлении и имеет две фазы, соответствующие твердому и жидкому состояниям. Математическое описание процесса затвердевания основано на модели фазового поля, обобщающей подход Плаппа (M. Plapp, Phys. Rev. E 84, 031601 (2011)) на случай локально-неравновесных процессов. Для вывода термодинамически согласованных уравнений модели использован метод расширенной необратимой термодинамики в отличие от феноменологического подхода Плаппа. Другое различие с моделью Плаппа состоит в использовании в качестве динамической переменной концентрации, а не химпотенциала примеси. В рамках полученной модели показана эквивалентность описания процесса затвердевания через концентрационное поле и через химпотенциал системы. В силу малости времен релаксации представленная модель сводится к сингулярно-возмущенной системе уравнений в частных производных параболического типа, описывающих динамику фазового и концентрационного полей. В работе предполагается известным описание термодинамических равновесных состояний на основе экспериментально полученных потенциалов Гиббса.

Для проверки полученной модели проведено численное моделирование одномерной задачи затвердевания в приближении разбавленного расплава Si-As, ранее неоднократно исследовавшегося экспериментально. Чтобы численно решить систему сингулярно-возмущенных уравнений, в работе предложен градиентно-устойчивый явный метод интегрирования уравнений второго порядка точности по времени. Для сведения бесконечного пространственного интервала к конечному использован метод «периодического сдвига». Оценка устойчивости получена из численных экспериментов.

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

We consider a locally nonequilibrium process of solidification for a supercooled binary melt. For sake of simplicity, it is assumed, that the solidifying binary system is at constant temperature and pressure. Also there are two phases corresponding to the solid and the liquid states. The mathematical description of the solidification process is based on the phase-field model that generalizes the approach of Plapp (M. Plapp, Phys. Rev. E 84, 031601 (2011)) to the case of locally nonequilibrium processes. We use the method of extended irreversible thermodynamics to derive thermodynamically consistent equations of the model, in contrast to the phenomenological approach of Plapp. A concentration as a dynamic variable (and not the chemical potential of the impurity) is another difference from Plapp's model. The equivalence of describing the process of solidification through the concentration field and through the chemical potential of the system is shown in the framework of the resulting model. In view of the smallness of the relaxation times, the present model is reduced to the singular-perturbed system of partial differential parabolic equations describing the dynamics of concentration and phase fields. In the paper, it is assumed that the description of the thermodynamic equilibrium states on the basis of the experimentally obtained Gibbs potentials is given.

To verify the model, the numerical simulation of the one-dimensional problem of solidification of the melt was performed in the approximation of the diluted melt Si-As, which had been repeatedly investigated experimentally. In this paper, we propose a gradient-stable explicit method of integrating equations of the second order of accuracy in time in order to solve the system of singularly-perturbed equations numerically. We reduced an infinite space interval to a finite interval by the method of «periodic translation». The estimation of stability was performed using numerical experiments.

The concentration profile, the phase-field profile and the distribution coefficient of the impurity at the front of solidification depending upon the value of supercooling were obtained from the numerical simulation of the solidification process for diluted melt Si-As. An analytical expression for the distribution coefficient as a function of supercooling that follows from the locally nonequilibrium model with a sharp interface was used to test the adequacy of the results of numerical experiments. The effect of the model parameters on the solidification process and behavior of the numerical solutions near the diffuse boundary were investigated.