Все выпуски

В работе изучаются неустойчивости Рихтмайера–Мешкова и Рэлея–Тейлора вязкопластических жидкостей (или, в частности, бингамовских жидкостей, обладающих предельным напряжением сдвига) в трехмерной постановке задачи. Анализируется развитие неустойчивостей Рихтмайера–Мешкова и Рэлея–Тейлора бингамовских жидкостей при одномодовом возмущении скорости контактной границы. Анализ проводится на основе численного моделирования с использованием метода Мак-Кормака и метода объема жидкости (метода VOF — Volume of Fluid) для отслеживания контактной границы в различные моменты времени. Представлены результаты численного моделирования неустойчивостей Рихтмайера–Мешкова и Рэлея–Тейлора бингамовской жидкости и их сравнение как с теорией, так и с результатами моделирования ньютоновской жидкости. В результате проведенных численных расчетов показано, что предел текучести вязкопластической жидкости существенно влияет на характер неустойчивости как Рэлея–Тейлора, так и Рихтмайера–Мешкова: существует критическая амплитуда начального возмущения поля скорости контактной границы, при превышении которой начинается развитие неустойчивостей. Если амплитуда начального возмущения поля скорости меньше критического значения, то это возмущение относительно быстро затухает и развития неустойчивостей не происходит. При превышении начальным возмущением критической амплитуды характер развития неустойчивостей напоминает таковой у ньютоновской жидкости. При рассмотрении неустойчивости Рихтмайера–Мешкова оцениваются критические амплитуды начального возмущения поля скорости контактной границы при различных значениях предельного напряжения сдвига бингамовской жидкости. Кроме того, наблюдается отличие поведения неньютоновской жидкости при развитии неустойчивости от плоского случая: при одном и том же зна- чении предельного напряжения сдвига в трехмерной геометрии интервал значений амплитуды начального возмущения, при котором происходит переход от покоя к движению, несколько уже. Помимо этого показано, что критическая амплитуда начального возмущения контактной границы для неустойчивости Рэлея–Тейлора ниже, чем для неустойчивости Рихтмайера–Мешкова. Это объясняется действием силы тяжести, «помогающей» развитию неустойчивости и противодействующей силам вязкого трения.

The Richtmyer–Meshkov and the Rayleigh–Taylor instabilities of viscoplastic (or the Bingham) fluids are studied in the three–dimensional formulation of the problem. A numerical modeling of the intermixing of two fluids with different rheology, whose densities differ twice, as a result of instabilities development process has been carried out. The development of the Richtmyer–Meshkov and the Rayleigh–Taylor instabilities of the Bingham fluids is analyzed utilizing the MacCormack and the Volume of Fluid (VOF) methods to reconstruct the interface during the process. Both the results of numerical simulation of the named instabilities of the Bingham liquids and their comparison with theory and the results of the Newtonian fluid simulation are presented. Critical amplitude of the initial perturbation of the contact boundary velocity field at which the development of instabilities begins was estimated. This critical amplitude presents because of the yield stress exists in the Bingham fluids. Results of numerical calculations show that the yield stress of viscoplastic fluids essentially affects the nature of the development of both Rayleigh–Taylor and Richtmyer–Meshkov instabilities. If the amplitude of the initial perturbation is less than the critical value, then the perturbation decays relatively quickly, and no instability develops.When the initial perturbation exceeds the critical amplitude, the nature of the instability development resembles that of the Newtonian fluid. In a case of the Richtmyer–Meshkov instability, the critical amplitudes of the initial perturbation of the contact boundary at different values of the yield stress are estimated. There is a distinction in behavior of the non-Newtonian fluid in a plane case: with the same value of the yield stress in three-dimensional geometry, the range of the amplitude values of the initial perturbation, when fluid starts to transit from rest to motion, is significantly narrower. In addition, it is shown that the critical amplitude of the initial perturbation of the contact boundary for the Rayleigh–Taylor instability is lower than for the Richtmyer–Meshkov instability. This is due to the action of gravity, which helps the instability to develop and counteracts the forces of viscous friction.

Рассматривается применение консервативного численного метода потоков для изучения вихревых структур в массивных, быстровращающихся компактных астрофизических объектах, находящихся в условиях самогравитации. Моделирование осуществляется для объектов с различной массой и скоростью вращения. Визуализируются картины вихревой структуры объектов. В расчетах используется газодинамическая модель, в которой газ принимается совершенным и невязким. Численная методика основана на конечно-разностной аппроксимации законов сохранения аддитивных характеристик среды для конечного объема. При этом используются upwind-аппроксимации плотностей распределения массы, компонент импульса и полной энергии. Для моделирования объектов, обладающих быстрым вращением, при эволюционном расчете осуществляется контроль сохранения компонент момента импульса, законы сохранения для которых не входят в систему основных уравнений. Эволюционный расчет осуществляется на основе параллельных алгоритмов, реализованных на вычислительном комплексе кластерной архитектуры. Алгоритмы основаны на стандартизованной системе передачи сообщений Message Passing Interface (MPI). При этом используются как блокирующие, так и неблокирующие процедуры обмена с контролем завершения операций. Осуществляется распараллеливание по пространству по двум или трем направле- ниям в зависимости от размера области интегрирования и параметров вычислительной сетки. Одновременно с распараллеливанием по пространству для каждой подобласти осуществляется распараллеливание по физическим факторам: расчет конвективного переноса и гравитационных сил реализуется параллельно на разных процессорах, что позволяет повысить эффективность алгоритмов. Показывается реальная возможность прямого вычисления гравитационных сил посредством суммирования взаимодействия между всеми конечными объемами в области интегрирования. Для методов конечного объема такой подход кажется более последовательным, чем решение уравнения Пуассона для гравитационного потенциала. Численные расчеты осуществлялись на вычислительном комплексе кластерной архитектуры с пиковой производительностью 523 TFlops. В расчетах использовалось до тысячи процессоров.

The application of a conservative numerical method of fluxes is examined for studying the vortex structures in the massive, fast-turned compact astrophysical objects, which are in self-gravity conditions. The simulation is accomplished for the objects with different mass and rotational speed. The pictures of the vortex structure of objects are visualized. In the calculations the gas-dynamic model is used, in which gas is accepted perfected and nonviscous. Numerical procedure is based on the finite-difference approximation of the conservation laws of the additive characteristics of medium for the finite volume. The “upwind” approximations of the densities of distribution of mass, components of momentum and total energy are applied. For the simulation of the objects, which possess fast-spin motion, the control of conservation for the component of moment of momentun is carried out during calculation. Evolutionary calculation is carried out on the basis of the parallel algorithms, realized on the computer complex of cluster architecture. Algorithms are based on the standardized system of message transfer Message Passing Interface (MPI). The blocking procedures of exchange and non-blocking procedures of exchange with control of the completion of operation are used. The parallelization on the space in two or three directions is carried out depending on the size of integration area and parameters of computational grid. For each subarea the parallelization based on the physical factors is carried out also: the calculations of gas dynamics part and gravitational forces are realized on the different processors, that allows to raise the efficiency of algorithms. The real possibility of the direct calculation of gravitational forces by means of the summation of interaction between all finite volumes in the integration area is shown. For the finite volume methods this approach seems to more consecutive than the solution of Poisson’s equation for the gravitational potential. Numerical calculations were carried out on the computer complex of cluster architecture with the peak productivity 523 TFlops. In the calculations up to thousand processors was used.

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

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

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

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

The paper describes an analytical study of hydrodynamic processes taking place in the course of coolant flow through a fuel assembly of the core of a fast neutron sodium-cooled reactor. Within the framework of the study, a procedure and an analytical model were developed based on program complex FlowVision of computational fluid dynamics, which, using proved simplifications, permits to obtain a coefficient of rod lifting margin of a fuel assembly and to study hydrodynamic characteristics of processes taking place in the course of simulation of different initial events influencing motion of a reactor core fuel assembly.

For analytical justification a fuel assembly model was developed, which is equivalent by hydraulic resistance values and permits not to simulate explicitly a complicated full-scale fuel assembly design, thus, decreasing a number of computational cells in the model and, as a result, reducing computational and time resources.

Hydraulic parameters of the equivalent fuel assembly model in program complex FlowVision were analyzed in two stages. At the first stage, to determine the minimum rod lifting margin coefficient of a fuel assembly, steady-state analyses were performed, where various flowrate values were assigned at the model inlet and forces acting upon the assembly were analyzed. A series of dynamic mode analyses was performed at the second stage. Jump-like pressure increase being the initial event which could occur hypothetically in the fast neutron sodium cooled reactor plant was assigned in these modes. Hydrodynamic parameters and forces acting upon the fuel assembly were determined.

The results of the first stage of the analytical study proved the minimum coefficient of rod lifting margin of a fuel assembly of the fast neutron reactor justified in reactor plant design documentation. As a result of the second stage of the study, conclusions were made on impossibility for the fuel assembly to move at the initial event associated with jump-like pressure increase in the reactor pressure chamber.

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

The article focuses on the problem of mechanisms’ co-design for robotic systems to perform adaptive physical interaction with an unstructured environment, including physical human robot interaction. The co-design means simultaneous optimization of mechanics and control system, ensuring optimal behavior and performance of the system. Mechanics optimization refers to the search for optimal structure, geometric parameters, mass distribution among the links and their compliance; control refers to the search for motion trajectories for mechanism’s joints. The paper presents a generalized method of structural-parametric synthesis of underactuated mechanisms with closed kinematics for robotic systems for various purposes, e. g., it was previously used for the co-design of fingers’ mechanisms for anthropomorphic gripper and legs’ mechanisms for galloping robots. The method implements the concept of morphological computation of control laws due to the features of mechanical design, minimizing the control effort from the algorithmic component of the control system, which reduces the requirements for the level of technical equipment and reduces energy consumption. In this paper, the proposed method is used to optimize the structure and geometric parameters of the passive mechanism of the back support module of an industrial exosuit. Human movements are diverse and non-deterministic when compared with the movements of autonomous robots, which complicates the design of wearable robotic devices. To reduce injuries, fatigue and increase the productivity of workers, the synthesized industrial exosuit should not only compensate for loads, but also not interfere with the natural human motions. To test the developed exosuit, kinematic datasets from motion capture of an entire human body during industrial operations were used. The proposed method of structural-parametric synthesis was used to improve the ergonomics of a wearable robotic device. Verification of the synthesized mechanism was carried out using simulation: the passive module of the back is attached to two geometric primitives that move the chest and pelvis of the exosuit operator in accordance with the motion capture data. The ergonomics of the back module is quantified by the distance between the joints connecting the upper and bottom parts of the exosuit; minimizing deviation from the average value corresponds to a lesser limitation of the operator’s movement, i. e. greater ergonomics. The article provides a detailed description of the method of structural-parametric synthesis, an example of synthesis of an exosuit module and the results of simulation.

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

Based on the modeling method using a mesh system, an algorithm is implemented for solving a unsteady problem with moving bodies The algorithm takes into account the movement and rotation of bodies according to a given law of motion. The algorithm is applied to analysis the flow around an infinite composed of cylinders with an elliptical cross-section, which either move across the flow or rotate with a change in the angle of attack. To simulate the flow of bodies with a sharp edge, characteristic of the profiles of gas turbine machines, an algorithm for constructing a mesh of type C with the inclusion of a certain area behind the profile is implemented. The program for modeling the flow near the profile is implemented within the framework of models of Euler equations, Navier – Stokes equations in the approximation of a thin layer with laminar viscosity and turbulent viscosity in the framework of an algebraic viscosity model. The program has also been adapted to solve the problems of internal gas dynamics of turbomachines. For this purpose, the method of setting the boundary conditions at the entrance and exit from the calculated area from the velocity to the pressure drop, as well as at the lateral boundaries from the free flow to the periodicity, was changed. This made it possible to simulate the flow of gas in the inter-blade channels of compressors and turbines of gas turbine engines. To refine the algorithm, a series of calculations of the aerodynamic parameters of several turbine cascades in various subsonic and supersonic modes and their comparison with the experiment were carried out. Calculations of turbine grating parameters were carried out within the framework of the inviscid and viscous gas model. The calculation and experiment were compared by the distribution of gas parameters near the profile, as well as by the energy losses of the flow in the cascade. Calculations have shown the applicability and correctness of the program to solve this class of problems. To test the program on the problems of external subsonic aerodynamics, calculations of the aerodynamic characteristics of an isolated airfoil in an undisturbed flow were performed. The results obtained allow us to assert the applicability of the hybrid grid method to various classes of problems of applied gas dynamics.

Предложен новый набор ключевых баллистических параметров: b₀ = tgθ₀, θ₀ — угол вылета, R_a — вершинный радиус кривизны траектории и β₀ — безразмерный квадрат разворотной скорости, и на его основе разработан новый прием приближенного интегрирования уравнений динамики материальной точки в среде с квадратичным сопротивлением (α = R/mg = 0,5…1,5) при tgθ₀ < 0,5. Способ базируется на преобразованиях Лежандра, и он дает формулы с автоматически подстраиваемой точностью как для текущих координат x(b), y(b) и времени t(b), b = tgθ — текущий наклон траектории, так и для основных параметров (время T, дальность L, положение вершины L_a) траектории в диапазоне, далеко выходящем за малоугловую область прицельной стрельбы. Точность формул выверялась при помощи продукта Maple.

New key parameters, namely b₀ = tgθ₀, θ₀ — angle of throwing, R_a — top curvature radius and β₀ — dimensionless speed square on the top of low angular trajectory were suggested in classic problem of integrating nonlinear equations of point mass projectile motion with quadratic air drag. Very precise formulae were obtained in a new way for coordinates x(b), y(b) and fly time t(b), b = tgθ where θ is inclination angle. This method is based on Legendre transformation and its precision is automatically improved in wide range of the θ₀ values and drag force parameters α. The precision was monitored by Maple computing product.

В различных областях науки, техники, защиты окружающей среды, в строительстве актуальными являются вопросы изучения процессов взаимодействия пористых материалов с веществами, находящимися в различных агрегатных состояниях. Особенно актуальными с точки зрения экологии и защиты окружающей среды являются исследования процессов взаимодействия пористых материалов с водой в жидкой и газообразной фазе. Поскольку в одном моле воды содержится 6.022140857 · 10²³ молекул H₂O, для описания свойств, например, водяного пара в поре в основном используются макроскопические подходы, рассматривающие водяной пар как сплошную среду в рамках аэродинамики. В данной работе построена и использовалась для моделирования макроскопическая двумерная диффузионная модель [Bitsadze, Kalinichenko, 1980] поведения водяного пара внутри изолированной поры. Наряду с макроскопической моделью в работе предложена микроскопическая модель поведения водяного пара внутри изолированной поры, построенная в рамках молекулярно-динамического подхода [Gould et al., 2005]. В данной модели на основе классической механики Ньютона описывается движение каждой молекулы воды, взаимодействующей как с другими молекулами воды, так и со стенками поры. Рассматривается эволюция системы «водяной пар – пора» с течением времени. В зависимости от внешних по отношению к поре условий система эволюционирует к различным состояниям равновесия, которые характеризуются различными значениями макроскопических характеристик, таких как температура, плотность, давление. Сравнение результатов молекулярно-динамического моделирования с результатами вычислений на основе макроскопической диффузионной модели и экспериментальными данными позволяет сделать вывод о необходимости сочетания макроскопического и микроскопического подхода для адекватного и более точного описания процессов взаимодействия водяного пара с пористыми материалами.

In various areas of science, technology, environment protection, construction, it is very important to study processes of porous materials interaction with different substances in different aggregation states. From the point of view of ecology and environmental protection it is particularly actual to investigate processes of porous materials interaction with water in liquid and gaseous phases. Since one mole of water contains 6.022140857 · 10²³ molecules of H₂O, macroscopic approaches considering the water vapor as continuum media in the framework of classical aerodynamics are mainly used to describe properties, for example properties of water vapor in the pore. In this paper we construct and use for simulation the macroscopic two-dimensional diffusion model [Bitsadze, Kalinichenko, 1980] describing the behavior of water vapor inside the isolated pore. Together with the macroscopic model it is proposed microscopic model of the behavior of water vapor inside the isolated pores. This microscopic model is built within the molecular dynamics approach [Gould et al., 2005]. In the microscopic model a description of each water molecule motion is based on Newton classical mechanics considering interactions with other molecules and pore walls. Time evolution of “water vapor – pore” system is explored. Depending on the external to the pore conditions the system evolves to various states of equilibrium, characterized by different values of the macroscopic characteristics such as temperature, density, pressure. Comparisons of results of molecular dynamic simulations with the results of calculations based on the macroscopic diffusion model and experimental data allow to conclude that the combination of macroscopic and microscopic approach could produce more adequate and more accurate description of processes of water vapor interaction with porous materials.

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

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

В качестве тела рассмотрен однородный эллиптический цилиндр с килем, расположенным вдоль большей диагонали. Внутри цилиндра расположены две материальные точечные массы, перемещающиеся по окружностям. Центры окружностей лежат на наименьшей диагонали эллипса на равном удалении от центра масс.

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

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

This article is devoted to the study of self-propulsion of bodies in a fluid by the action of internal mechanisms, without changing the external shape of the body. The paper presents an overview of theoretical papers that justify the possibility of this displacement in ideal and viscous liquids.

A special case of self-propulsion of a rigid body along the surface of a liquid is considered due to the motion of two internal masses along the circles. The paper presents a mathematical model of the motion of a solid body with moving internal masses in a three-dimensional formulation. This model takes into account the three-dimensional vibrations of the body during motion, which arise under the action of external forces-gravity force, Archimedes force and forces acting on the body, from the side of a viscous fluid.

The body is a homogeneous elliptical cylinder with a keel located along the larger diagonal. Inside the cylinder there are two material point masses moving along the circles. The centers of the circles lie on the smallest diagonal of the ellipse at an equal distance from the center of mass.

Equations of motion of the system (a body with two material points, placed in a fluid) are represented as Kirchhoff equations with the addition of external forces and moments acting on the body. The phenomenological model of viscous friction is quadratic in velocity used to describe the forces of resistance to motion in a fluid. The coefficients of resistance to movement were determined experimentally. The forces acting on the keel were determined by numerical modeling of the keel oscillations in a viscous liquid using the Navier – Stokes equations.

In this paper, an experimental verification of the proposed mathematical model was carried out. Several series of experiments on self-propulsion of a body in a liquid by means of rotation of internal masses with different speeds of rotation are presented. The dependence of the average propagation velocity, the amplitude of the transverse oscillations as a function of the rotational speed of internal masses is investigated. The obtained experimental data are compared with the results obtained within the framework of the proposed mathematical model.

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

На основе базовых соотношений распада разрыва в равновесной гетерогенной среде получено асимптотически точное автомодельное решение задачи взаимодействия ударной волны со слоем газовзвеси, к которому сходится численное решение двухскоростной двухтемпературной динамики газовзвеси при уменьшении размеровди сперсных частиц.

Изучены динамика движения скачка уплотнения в газе и его взаимодействия с ограниченным слоем газовзвеси для различных размеров дисперсных частиц: 0.1, 2 и 20 мкм. Задача характеризуется двумя распадами разрывов: отраженной и преломленной ударными волнами на левой границе слоя, отраженной волной разрежения и прошедшим скачком уплотнения на правой контактной границе. Обсуждено влияние релаксационных процессов (безразмерных времен релаксации фаз) на характер течения газовзвеси. Для мелких частиц времена выравнивания скоростей и температур фаз малы, а зоны релаксации являются подсеточными. Численное решение в характерных точках с относительной точностью $O\, (10^{−4})$  сходится к автомодельным решениям.

For a non-homogeneous model transport equation with source terms, the stability analysis of a linear hybrid scheme (a combination of upwind and central approximations) is performed. Stability conditions are obtained that depend on the hybridity parameter, the source intensity factor (the product of intensity per time step), and the weight coefficient of the linear combination of source power on the lower- and upper-time layer. In a nonlinear case for the non-equilibrium by velocities and temperatures equations of gas suspension motion, the linear stability analysis was confirmed by calculation. It is established that the maximum permissible Courant number of the hybrid large-particle method of the second order of accuracy in space and time with an implicit account of friction and heat exchange between gas and particles does not depend on the intensity factor of interface interactions, the grid spacing and the relaxation times of phases (K-stability). In the traditional case of an explicit method for calculating the source terms, when a dimensionless intensity factor greater than 10, there is a catastrophic (by several orders of magnitude) decrease in the maximum permissible Courant number, in which the calculated time step becomes unacceptably small.

On the basic ratios of Riemann’s problem in the equilibrium heterogeneous medium, we obtained an asymptotically exact self-similar solution of the problem of interaction of a shock wave with a layer of gas-suspension to which converge the numerical solution of two-velocity two-temperature dynamics of gassuspension when reducing the size of dispersed particles.

The dynamics of the shock wave in gas and its interaction with a limited gas suspension layer for different sizes of dispersed particles: 0.1, 2, and 20 ìm were studied. The problem is characterized by two discontinuities decay: reflected and refracted shock waves at the left boundary of the layer, reflected rarefaction wave, and a past shock wave at the right contact edge. The influence of relaxation processes (dimensionless phase relaxation times) to the flow of a gas suspension is discussed. For small particles, the times of equalization of the velocities and temperatures of the phases are small, and the relaxation zones are sub-grid. The numerical solution at characteristic points converges with relative accuracy $O \, (10^{-4})$ to self-similar solutions.

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

The paper deals with the mathematical model formulation for studying the nonlinear hydro-elastic response of the narrow channel wall supported by a spring with cubic nonlinearity and interacting with a pulsating viscous liquid filling the channel. In contrast to the known approaches, within the framework of the proposed mathematical model, the inertial and dissipative properties of the viscous incompressible liquid and the restoring force nonlinearity of the supporting spring were simultaneously taken into account. The mathematical model was an equations system for the coupled plane hydroelasticity problem, including the motion equations of a viscous incompressible liquid, with the corresponding boundary conditions, and the channel wall motion equation as a single-degree-of-freedom model with a cubic nonlinear restoring force. Initially, the viscous liquid dynamics was investigated within the framework of the hydrodynamic lubrication theory, i. e. without taking into account the liquid motion inertia. At the next stage, the iteration method was used to take into account the motion inertia of the viscous liquid. The distribution laws of the hydrodynamic parameters for the viscous liquid in the channel were found which made it possible to determine its reaction acting on the channel wall. As a result, it was shown that the original hydroelasticity problem is reduced to a single nonlinear equation that coincides with the Duffing equation. In this equation, the damping coefficient is determined by the liquid physical properties and the channel geometric dimensions, and taking into account the liquid motion inertia lead to the appearance of an added mass. The nonlinear equation study for hydroelastic oscillations was carried out by the harmonic balance method for the main frequency of viscous liquid pulsations. As a result, the primary steady-state hydroelastic response for the channel wall supported by a spring with softening or hardening cubic nonlinearity was found. Numerical modeling of the channel wall hydroelastic response showed the possibility of a jumping change in the amplitudes of channel wall oscillations, and also made it possible to assess the effect of the liquid motion inertia on the frequency range in which these amplitude jumps are observed.