Все выпуски

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

The formation process of a spherical shock impulse in gas and its interaction with a protective aqueous foam barrier, accompanied by formation of vortex flows, are numerically investigated. The problem is solved in a two-dimensional axisymmetric formulation using a two-phase model of a gas-liquid mixture, which includes the laws of conservation of mass, momentum and energy of the mixture and an equation for the dynamics of volume content of phases.The numerical implementation is carried out on the basis of the OpenFOAM package using the standard compressibleMultiphaseInterFoam solver, modified in accordance with the conditions of the problem and model representations. The discretization of the system of equations in the chosen solver is carried out by the method of finite volumes using the computational Pimple algorithm. A significant decrease in the intensity of the shock wave in its interaction with the aqueous foam barrier is shown and the causes leading to vortex formation in the gas region are revealed. The reliability of the results obtained is estimated by comparison with solutions of a similar problem by other numerical methods.

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

At the river bed curves, secondary flow normal to the main flow direction are formed. Depending on the channel geometry, there may be several secondary flows in the cross-section, and they may have different scales. Even a small secondary cross-section flow affects the parameters of the hydrodynamic flow and this influence must be taken into account when modeling riverbed processes and researching coast deformations of the channel. Three-dimensional modeling of such multi-scale processes requires large computational costs and is currently possible only for small model channels. Therefore, a reduced-dimensional model is proposed in this paper to study coastal processes. The performed reduction of the problem from a three-dimensional model of river flow motion to a two-dimensional one in the plane of the channel cross-section assumes that the hydrodynamic flow is quasi-stationary and the hypotheses on the asymptotic behavior of the flow along the flow coordinate are fulfilled for it. Taking into account these limitations, a mathematical model of the problem of a stationary turbulent calm river flow in a channel cross-section is formulated in this work. The problem is formulated in a mixed velocity–vortex–stream function formulation. Specifying of the boundary conditions on the flow free surface for the velocity field determined in the normal and tangential directions to the cross-section axis is required as additional conditions for the problem reduction. It is assumed that the values of this velocity field should be determined from the solution of auxiliary problems or obtained from data of natural or experimental measurements.

The finite element method in the Petrov–Galerkin formulation is used for the numerical solution of the formulated problem. A discrete analog of the problem is obtained and an algorithm for its solution is proposed. The performed numerical studies showed generally good agreement between the obtained solutions and the known experimental data. The authors associate the errors in the numerical results with the need for a more accurate determination of the radial component of the velocity field in the cross-section by selecting and calibrating a more suitable model for turbulent viscosity calculating and a more accurate determination of the boundary conditions on the cross-section free boundary.

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

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.

Данная работа посвящена исследованию динамики следующих систем большого числа точечных вихрей на плоскости:
- вихревые кольца с внешним радиусом r=1 и переменным внутренним радиусом r₀,
- вихревые эллипсы с полуосями a, b.
Основное внимание уделено изучению асимптотического поведения (t→∞) систем и проверке критериев устойчивости для непрерывных распределений завихренности с помощью компьютерного эксперимента.

With the help of mathematical modelling, we study the dynamics of many point vortices system on the plane. For this system, we consider the following cases:
- vortex rings with outer radius r=1 and variable inner radius r₀,
- vortex ellipses with semiaxes a, b.
The emphasis is on the analysis of the asymptotic (t→∞) behavior of the system and on the verification of the stability criteria for vorticity continuous distributions.