Все выпуски

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

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.

На основе кусочно-квадратичной интерполяции получены полуаналитические аппроксимации потенциала двойного слоя вблизи и на границе двумерной области. Для вычисления интегралов, образующихся после интерполяции функции плотности, используется точное интегрирование по переменной $\rho=\left(r^2-d^2\right)^{1/2}$, где $d$ и $r$ — расстояния от наблюдаемой точки до границы области и до граничной точки интегрирования соответственно. Доказана устойчивая сходимость таких аппроксимаций с кубической скоростью равномерно вблизи границы класса $C^5$, а также на самой границе. Также доказано, что использование для вычисления интегралов стандартных квадратурных формул не нарушает равномерной кубической сходимости аппроксимаций прямого значения потенциала на границе класса $C^6$. При некоторых упрощениях доказано, что использование для вычисления интегралов стандартных квадратурных формул влечет отсутствие равномерной сходимости аппроксимаций потенциала внутри области вблизи любой граничной точки. Теоретические выводы подтверждены результатами численного решения задачи Дирихле для уравнения Лапласа в круговой области.

On the basis of piecewise quadratic interpolation, semi-analytical approximations of the double layer potential near and on the boundary of a two-dimensional domain are obtained. To calculate the integrals formed after the interpolation of the density function, exact integration with respect to the variable $\rho=\left(r^2-d^2\right)^{1/2}$ is used, where $d$ and $r$ are the distances from the observed point to the boundary of the domain and to the boundary point of integration, respectively. The study proves the stable convergence of such approximations with the cubic velocity uniformly near the boundary of the class $C^5$, and also on the boundary itself. It is also proved that the use of standard quadrature formulas for calculating the integrals does not violate the uniform cubic convergence of approximations of the direct value of the potential on the boundary of the class $C^6$. With some simplifications, it is proved that the use of standard quadrature formulas for calculating the integrals entails the absence of uniform convergence of potential approximations inside the domain near any boundary point. The theoretical conclusions are confirmed by the results of the numerical solution of the Dirichlet problem for the Laplace equation in a circular domain.

На основе кусочно-квадратичной интерполяции получены полуаналитические аппроксимации нормальной производной потенциала простого слоя вблизи и на границе двумерной области. Для вычисления интегралов, образующихся после интерполяции функции плотности, используется точное интегрирование по переменной $\rho =(r^{2} -d^{2} )^{1/2} $, где $d$ и $r$ — расстояния от наблюдаемой точки до границы области и до граничной точки интегрирования соответственно. Доказана устойчивая сходимость таких аппроксимаций с кубической скоростью равномерно вблизи границы класса $C^{5}$, а также на самой границе. Также доказано, что на границе аппроксимации по аналогии с точной функцией терпят разрыв, величина которого пропорциональна значениям интерполированной функции плотности, но могут быть доопределены на границе до функций, непрерывных или на замкнутой внутренней, или на замкнутой внешней приграничной области. Теоретические выводы о равномерной сходимости подтверждены результатами вычисления нормальной производной вблизи границы единичного круга.

On the basis of piecewise quadratic interpolation, semi-analytical approximations of the normal derivative of the simple layer potential near and on the boundary of a two-dimensional domain are obtained. To calculate the integrals formed after the interpolation of the density function, exact integration over the variable $\rho=(r^{2}-d^{2})^{1/2} $ is used, where $d$ and $r$ are the distances from the observed point to the boundary of the domain and to the boundary point of integration, respectively. The study proves the stable convergence of such approximations with cubic velocity uniformly near the boundary of the class $C^{5}$, as well as on the boundary itself. It is also proved that, by analogy with the exact function, the approximations suffer a discontinuity at the boundary, the magnitude of which is proportional to the values of the interpolated density function, but they can be extended on the boundary to functions that are continuous either on a closed internal border domain or on a closed external one. Theoretical conclusions about uniform convergence are confirmed by the results of calculating the normal derivative near the boundary of a unit circle.

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

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.

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

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.

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

Mathematical models were developed and the nonlinear boundary value problem of dynamics thinwalled shells of the arbitrary form under action shock pulse is formulated. Dependence of processes of deformation on speed loading, compressibility of a material, finite deformations and large displacements of a shell middle surface, formation and kinetic of plasticity zones of a material during action of a shock wave are considered. Parameterization of a shell surface is carried out by bi-cubic splines. For the description of nonlinear, time and speed dependents of a shell material behavior with anisotropic hardening the generalized model of microplasticity is developed on the account of viscosity of deformation, hysteresis losses and Baushinger's effect. The solution of boundary value problems on the basis of difference schemes is constructed. Results of modeling of nonlinear wave processes in a assemble shell under action of explosion also are presented.