Численное моделирование неоднородного течения Куэтта - Пуазейля с учетом трения Рэлея

K V Gubareva; E Yu Prosviryakov

doi:10.22213/2410-9304-2026-1-78-88

Authors

K. V. Gubareva Samara State Technical University
E. Y. Prosviryakov Institute of Engineering Science UB RAS

DOI:

https://doi.org/10.22213/2410-9304-2026-1-78-88

Keywords:

Rayleigh friction, numerical modeling, reverse flows, finite difference method, viscous fluid, plane channel, Couette-Poiseuille flow

Abstract

The results of numerical modeling of steady viscous incompressible fluid flow in a flat channel, taking into account linear Rayleigh friction, are presented. The study focuses on a generalized model of hybrid Couette-Poiseuille flow, in which the longitudinal velocity linearly depends on the transverse coordinate, while the coefficients of this expansion vary along the channel height. The mathematical model is based on modified Navier-Stokes equations with an additional dissipative term proportional to velocity, which allows accounting for distributed drag forces arising during flow through porous media or in the presence of dispersed phases. A second-order finite difference method implemented in the MATLAB 2023b environment was used to solve the boundary value problem. Verification of the numerical method was carried out by comparing the results with analytical solutions for limiting cases. The influence of the Rayleigh friction coefficient on the flow structure was investigated for two technical fluids with different rheological properties: polyalphaolefin PAO-40 and Tyfocor LS heat transfer fluid at 40°C. The results demonstrate that an increase in the friction coefficient leads to a significant reduction in flow velocity and the formation of complex spatial structures. The emergence of reverse flow zones (counterflows) at certain parameter values was discovered. Different sensitivities of the fluids to changes in dissipative characteristics were revealed: the less viscous fluid exhibits a stronger dependence of velocity on the friction parameter. It is shown that qualitative changes in the flow structure are determined by the type of differential operator of the system. The obtained results have practical significance for the design of heat exchange systems, microfluidic devices, and technological equipment where additional dissipative mechanisms must be taken into account. The proposed calculation methodology allows optimizing the parameters of such systems considering volumetric drag forces.

Author Biographies

K. V. Gubareva, Samara State Technical University

Candidate of Engineering Sciences, Associate Professor

E. Y. Prosviryakov, Institute of Engineering Science UB RAS

Dr. Sci. (Phys.-Math.), Associate Professor, Professor, Ural Federal University; Head of Sector, Chief Researcher

References

Галкин В. А., Смородинов А. Д., Моргун Д. А. Решение уравнения Навье - Стокса для сталкивающихся потоков // Успехи кибернетики. 2023. Т. 4, № 2. С. 8-15. DOI: 10.51790/2712-9942-2023-4-2-01.

Ekman V.W. On the Influence of the Earth’s Rotation on Ocean-Currents. Arkiv för Matematik, Astronomi och Fysik, 1905, vol. 2, no. 11, pp. 1-53.

Губарева К. В., Просвиряков Е. Ю., Еремин А. В. Точное квадратичное полиномиальное решение для описания неоднородного течения Куэтта - Пуазейля в бесконечном горизонтальном слое с проницаемыми границами // Вестник Южно-Уральского государственного университета. Серия: Математика. Механика. Физика. 2026. Т. 18, № 1. С. 47-62. DOI: 10.14529/mmph260106.

Ландау Л. Д., Лифшиц Е. М. Гидродинамика. М. : Физматлит, 2001. 731 с.

Burmasheva N., Prosviryakov E. Exact Solutions to the Oberbeck-Boussinesq Equations for Shear Flows of a Viscous Binary Fluid with Allowance Made for the Soret Effect. The Bulletin of Irkutsk State University. Series: Mathematics, 2021, vol. 37, pp. 17-30. DOI 10.26516/1997-7670.2021.37.17.

Ferro C., Cafaro M., Maggiore P. Optimizing Solid Rocket Missile Trajectories: A Hybrid Approach Using an Evolutionary Algorithm and Machine Learning. Aerospace, 2024, vol. 11, p. 912. DOI 10.3390/aerospace11100912.

Аристов С. Н., Фрик П. Г. Крупномасштабная турбулентность в конвекции Рэлея - Бенара // Известия Академии наук СССР. Механика жидкости и газа. 1989. № 5. С. 43-48.

Бурмашева Н. В., Просвиряков Е. Ю. Крупномасштабная слоистая стационарная конвекция вязкой несжимаемой жидкости под действием касательных напряжений на верхней границе. Исследование полей температуры и давления // Вестник Самарского государственного технического университета. Серия: Физико-математические науки. - 2017. Т. 21, № 4. С. 736-751. DOI: 10.14498/vsgtu1568.

Полянин А. Д., Аристов С. Н. Системы уравнений гидродинамического типа: точные решения, преобразования, нелинейная устойчивость // Доклады Академии наук. 2009. Т. 428, № 2. С. 180-185.

Pedlosky J. Geophysical Fluid Dynamics. 2nd ed. New York, Springer, 1987, 710 p.

Полянин А. Д., Зайцев В. Ф. Справочник по нелинейным уравнениям математической физики. М. : Физматлит, 2002. 432 с.

Batchelor G.K. An Introduction to Fluid Dynamics. Cambridge, Cambridge University Press, 2000, 615 p.

Boyd J.P. Chebyshev and Fourier Spectral Methods. 2nd ed. New York, Dover Publications, 2001, 688 p.

Drazin P.G., Riley N. The Navier-Stokes Equations: A Classification of Flows and Exact Solutions. Cambridge, Cambridge University Press, 2006, 226 p.

Rayleigh L. On the Dynamics of Revolving Fluids. Proceedings of the Royal Society of London. Series A, 1916, vol. 93, no. 648, pp. 148-154.

Burmasheva N., Ershkov S., Prosviryakov E., Leshchenko D. Exact Solutions of Navier-Stokes Equations for Quasi-Two-Dimensional Flows with Rayleigh Friction. Fluids, 2023, vol. 8, p. 123. DOI 10.3390/fluids8040123