Все выпуски

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

We consider the model of chaotic motion of a plate in a viscous fluid, described by an oscillatory system of three ordinary differential equations with a quadratic nonlinearity. In the course of the bifurcation study of singular points of the system, maps of the types of singular points are constructed and a surface equation is found in the space of dissipation and circulation parameters on which the Andronov-Hopf bifurcation of the limit cycle creation takes place. With a further change in the parameters near the Andronov-Hopf surface, cascades of the period doubling doubling of the Feigenbaum cycle and the Sharkovsky subharmonic cascades, ending with the creation of a cycle of period three, are found. Expressions are obtained for saddle numbers of the saddle-node and two saddle-foci and their plots are plotted in the parameter space. It is shown that homoclinic cascades of bifurcations are realized in the system with the destruction of homoclinic trajectories of saddle-foci. The existence of homoclinic trajectories of saddle-foci is proved by a numerical-analytical method. The graphs of the largest Lyapunov exponent and the bifurcation diagrams show that when the dissipation coefficients change, the system switches to chaos in several stages.

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

This article studies the Keldysh problem for a three-dimensional equation of mixed type with three singular coefficients in a semi-infinite parallelepiped. Based on the completeness property of eigenfunction systems of two one-dimensional spectral problems, the uniqueness theorem is proved. To prove the existence of a solution to the problem, the Fourier spectral method based on the separation of variables is used. The solution to this problem is constructed in the form of a sum of a double Fourier-Bessel series. In substantiating the uniform convergence of the constructed series, we used asymptotic estimates of the Bessel functions of the real and imaginary argument. Based on them, estimates were obtained for each member of the series, which made it possible to prove the convergence of the series and its derivatives to the second order inclusive, as well as the existence theorem in the class of regular solutions

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

In this article, we study the growth of solutions of homogeneous and non-homogeneous complex linear differential equations where the coefficients are analytic functions in the extended complex plane except a finite singular point with finite logarithmic order. We extend some previous results obtained very recently by Fettouch and Hamouda.

Математическое моделирование композиционных материалов играет важную роль в современной технике, а решение и исследование обратных граничных задач теплообмена невозможно без использования систем собственных функций задачи Штурма-Лиувилля для дифференциального уравнения с разрывными коэффициентами. Одним из важнейших свойств таких систем является их полнота в соответствующих пространствах. Это свойство систем позволяет доказать теоремы существования и единственности как для прямых задач, так и обратных граничных задач теплопроводности, а также обосновать численные методы решения таких задач. В настоящей статье доказана полнота в пространстве $L_2[r_0,r_2]$ задачи Штурма-Лиувилля для дифференциального оператора второго порядка с разрывным коэффициентом. Эта задача возникает при исследовании и решении обратной граничной задачи теплопроводности для полого шара, состоящего из двух шаров с различными коэффициентами температуропроводности. Доказана самосопряженность, инъективность, а также положительная определенность этого оператора.

Mathematical modeling of composite materials plays an important role in modern technology, and the solution and study of inverse boundary value problems of heat transfer is impossible without the use of systems of eigenfunctions of the Sturm-Liouville problem for the differential equation with discontinuous coefficients. One of the most important properties of such systems is their completeness in the corresponding spaces. This property of systems allows to prove theorems of existence and uniqueness of both direct problems and inverse boundary value problems of thermal conductivity, and also to prove numerical methods of solving such problems. In this paper, we prove the completeness of the Sturm-Liouville problem in the space $L_2[r_0,r_2]$ for a second-order differential operator with a discontinuous coefficient. This problem arises when investigating and solving the inverse boundary problem of thermal conductivity for a hollow ball consisting of two balls with different temperature conductivity coefficients. Self-conjugacy, injectivity, and positive definiteness of this operator are proved.

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

We investigate the spectral properties of the discrete Schrödinger operator for the infinite band with zero boundary conditions. We prove that the eigenvalues and resonances arise for the small decreasing potentials near singularities of the non-perturbed Green function (boundary points of the subbands) and we find their asymptotic behavior. The scattering picture is described: the diffraction (i.e. the scattering mainly in the finite number of preferential directions) transforms into probability waves in time of the reflection and propagation in the considered quasi-1D system. The simple formulas for these probabilities are obtained near boundary points of the subbands (this corresponds to small velocities of the quantum particles) for the small potentials.

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

It is proved that a general cubic form over the field of complex numbers can be transformed into a form without monomials of exactly two variables by means of a non-degenerate linear transformation of coordinates. If the coefficients of monomials in only one variable are equal to one, and the remaining coefficients belong to sufficiently small polydisc near zero, then the transformation can be approximated by iterative algorithm. Under these restrictions the same result holds over the reals. This result generalizes the Levy-Desplanques theorem on strictly diagonally dominant matrices. We discuss in detail the properties of reducible cubic forms. So we prove the existence of a reducible real cubic form that is not equivalent to any form with all monomials in only one variable and without any monomials in exactly two variables. We suggest a sufficient condition for the existence of a singular point on a projective cubic hypersurface. The computational complexity of singular points recognition is discussed.

В статье рассмотрено параболо-гиперболическое уравнение с сингулярным коэффициентом и спектральным параметром в области, состоящей из характеристического треугольника и полуполосы. Сформулирована задача с нелокальным условием, связывающим значения искомой функции в точках двух граничных характеристик и линии изменения типа уравнения с помощью двух операторов, один из которых зависит от коэффициента сингулярности, а другой — от спектрального параметра. Поставленная задача исследована сведением ее к системе уравнений относительно следа искомой функции и еe производной по $x$ на линии изменения типа уравнения. Единственность решения доказана с использованием метода интегралов энергии, при этом использованы интегральные представления гамма-функции Эйлера и функции Бесселя первого рода. Существование решения задачи доказано методом интегральных уравнений, при этом поставленная задача эквивалентно сведена к интегральному уравнению Фредгольма второго рода, разрешимость которого следует из единственности решения задачи. Выявлены достаточные условия, которые обеспечивают однозначную разрешимость поставленной задачи.

In the paper, a parabolic-hyperbolic equation with a singular coefficient and a spectral parameter in the domain which consists of a characteristic triangle and a half strip has been considered. A nonlocal problem connecting the values of the desired function at the two points of boundary characteristics and the line of equation type changing by means of two operators, the first of which depends on the coefficient of the singularity and the second one - on the spectral parameters, is formulated. The considered problem is investigated by reducing it to the system of equations in the trace of the desired function and its derivative with respect to $x$ on the line of equation type changing. The uniqueness of the solution is proved by the method of energy integrals, for this we use integral representations of Euler gamma-function and Bessel function of the first kind. The existence of the solution is proved by the method of integral equations, for this we equivalently reduce the considered problem to the Fredholm integral equation of the second kind which solvability follows from the uniqueness of the problem solution. Sufficient conditions for unique solvability of the considered problem are found.

В работе рассмотрен локально-неравновесный процесс затвердевания переохлажденного бинарного расплава. В целях простоты предполагается, что затвердевающая бинарная система находится при постоянных температуре и давлении и имеет две фазы, соответствующие твердому и жидкому состояниям. Математическое описание процесса затвердевания основано на модели фазового поля, обобщающей подход Плаппа (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.