Все выпуски

Рассмотрено применение барицентрического метода для численного решения задач Дирихле и Неймана для уравнения Гельмгольца в ограниченной односвязной области $\Omega\subset\mathbb{R}^2$. Основное допущение в решении заключается в задании границы $\Omega$ в кусочно-линейном представлении. Отличительная особенность барицентрического метода состоит в порядке формирования глобальной системы векторных базисных функций для $\Omega$ через барицентрические координаты. Установлены существование и единственность решения задач Дирихле и Неймана для уравнения Гельмгольца барицентрическим методом и определена оценка скорости сходимости. Уточнены особенности алгоритмической реализации метода.

The application of the barycentric method for the numerical solution of Dirichlet and Neumann problems for the Helmholtz equation in the bounded simply connected domain $\Omega\subset\mathbb{R}^2$ is considered. The main assumption in the solution is to set the $\Omega$ boundary in a piecewise linear representation. A distinctive feature of the barycentric method is the order of formation of a global system of vector basis functions for $\Omega$ via barycentric coordinates. The existence and uniqueness of the solution of Dirichlet and Neumann problems for the Helmholtz equation by the barycentric method are established and the convergence rate estimate is determined. The features of the algorithmic implementation of the method are clarified.

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

We study a multidimensional case of a nonlinear reaction-diffusion system modeled by a system of two parabolic equations with power nonlinearities. Such systems can be used to simulate the process of propagation in space of interacting distributed formations of robots of two types. Such equations also describe the processes of nonlinear diffusion in reacting two-component continuous media. An original version of the reduction method is proposed, which reduces the construction of the dependence of the exact solution on spatial variables to the solution of the Helmholtz equation, and the dependence on time to the solution of a linear system of ordinary differential equations. A number of examples of multiparameter families of exact solutions given by elementary functions are constructed.

Для современной геометрии важное значение имеет изучение геометрий максимальной подвижности. Максимальная подвижность для $n$-мерной геометрии, задаваемой функцией $f$ пары точек означает существование $n(n+1)/2$-мерной группы преобразований, оставляющей эту функцию инвариантной. Известно много геометрий максимальной подвижности (геометрия Евклида, симплектическая, Лобачевского и т.д.), но полной классификации таких геометрий нет. В данной статье методом вложения решается одна из таких классификационных задач. Суть этого метода состоит в следующем: по известной функции пары точек $g$ трехмерной геометрии находим все невырожденные функции $f$ пары точек четырехмерных геометрий, являющиеся инвариантами группы Ли преобразований размерности 10. В этой статье $g$ - это невырожденные функции пары точек двух гельмгольцевых трехмерных геометрий: $$g = 2\ln(x_i-x_j) + \dfrac{y_i-y_j}{x_i-x_j}+2z_i+2z_j,$$ $$\ln[(x_i-x_j)^2+(y_i-y_j)^2]+ 2\gamma\,\text{arctg}\dfrac{y_i-y_j}{x_i-x_j}+2z_i+2z_j.$$ Данные геометрии локально максимально подвижны, то есть их группы движений шестимерны. Задача, решаемая в этой работе, сводится к решению аналитическими методами специальных функциональных уравнений, решения которых ищутся в виде рядов Тейлора. Для перебора различных вариантов применяется пакет математических программ Maple 15. В результате получаются только вырожденные функции пары точек.

For modern geometry, the study of maximum mobility geometries is important. The maximum mobility for $n$-dimensional geometry given by the function $f$ of a pair of points means the existence of an $n(n+1)/2$-dimensional transformation group, which leaves this function invariant. Many geometries of maximum mobility are known (Euclidean, symplectic, Lobachevsky, etc.), but there is no complete classification of such geometries. In this article, the method of embedding solves one of these classification problems. The essence of this method is as follows: from the function of a pair of points $ g $ of three-dimensional geometry, we find all non-degenerate functions $f$ of a pair of points of four-dimensional geometries that are invariants of the Lie group of transformations of dimension 10. In this article, $g$ are non-degenerate functions of a pair of points of two Helmholtz three-dimensional geometries: $$g = 2\ln(x_i-x_j) + \dfrac{y_i-y_j}{x_i-x_j} + 2z_i + 2z_j, $$ $$\ln [(x_i-x_j)^2 + (y_i-y_j)^2] + 2\gamma\,\text{arctg}\dfrac{y_i-y_j}{x_i-x_j} + 2z_i + 2z_j. $$ These geometries are locally maximally mobile, that is, their groups of motions are six-dimensional. The problem solved in this work is reduced to solving special functional equations by analytical methods, the solutions of which are sought in the form of Taylor series. For searching various options, the math software package Maple 15 is used. As a result, only degenerate functions of a pair of points are obtained.

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

The stability of incompressible fluid plane-parallel flow over a layer of a saturated porous medium is studied. The results of a linear stability analysis are described at different porosity values. The considered system is bounded by solid wall from the porous layer bottom. Top fluid surface is free and rigid. A linear stability analysis of plane-parallel stationary flow is presented. It is realized for parameter area where the neutral stability curves are bimodal. The porosity variation effect on flow stability is considered. It is shown that there is a transition between two main instability modes: long-wave and short-wave. The long-wave instability mechanism is determined by inflection points within the velocity profile. The short-wave instability is due to the large transverse gradient of flow velocity near the interface between liquid and porous medium. Porosity decrease stabilizes the long wave perturbations without significant shift of the critical wavenumber. Simultaneously, the short-wave perturbations destabilize, and their critical wavenumber changes in wide range. When the porosity is less than 0.7, the inertial terms in filtration equation and magnitude of the viscous stress near the interface increase to such an extent that the Kelvin-Helmholtz analogue of instability becomes the dominant mechanism for instability development. The stability band realizes in narrow porosity area. It separates the two branches of the neutral curve.

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

Swirling laminar axisymmetric flows of viscous incompressible fluids in a potential field of body forces are considered. The study of flows is carried out in a cylindrical coordinate system. In the flows, the regions in which the axial derivative of the circumferential velocity cannot take on zero value in some open neighborhood (essentially swirling flows) and the regions in which this derivative is equal to zero (the region with layered swirl) are considered separately. It is shown that a well-known method (the method of viscous vortex domains) developed for non-swirling flows can be used for regions with layered swirling. For substantially swirling flows, a formula is obtained for calculating the radial-axial velocity of an imaginary fluid through the circumferential vorticity component, the circumferential circulation of a real fluid, and the partial derivatives of these functions. The particles of this imaginary fluid “transfer” vortex tubes of the radial-axial vorticity component while maintaining the intensity of these tubes, and also “transfer” the circumferential circulation and the product of the circular vorticity component by some function of the distance to the axis of symmetry. A non-integral method for reconstructing the velocity field from the vorticity field is proposed. It is reduced to solving a system of linear algebraic equations in two variables. The obtained result is proposed to be used to extend the method of viscous vortex domains to swirling axisymmetric flows.