ISSN 2658–5782
DOI 10.21662
Electronic Scientific Journal


© Институт механики
им. Р.Р. Мавлютова
УФИЦ РАН


Galeeva D.R., Kireev V.N. Modeling the dynamics of a drop moving in a plane channel with a viscous barrier. Multiphase Systems. 2025;20(3):129–134 (in Russian).
2025. Vol. 20. Issue 3, Pp. 129–134
URL: http://mfs.uimech.org/mfs2025.3.017
DOI: 10.21662/mfs2025.3.017
Modeling the dynamics of a drop moving in a plane channel with a viscous barrier
D.R. Galeeva🖂, V.N. Kireev
Ufa University of Science and Technology, Ufa, Russia

Abstract

This article examines the influence of a highly viscous region on the deformation of a single droplet moving in a thermoviscous fluid flow with an anomalous viscosity-temperature dependence. Numerical modeling is performed using the phase field method, consisting of the Cahn-Hilliard-Navier-Stokes system of equations, supplemented by an equation for the temperature field. The results of numerical modeling for laminar flow with a steady-state temperature field in a flat channel are presented. The influence of a viscous barrier located in the middle of the channel on the droplet velocity and its deformation is analyzed for three capillary numbers. Graphs of the droplet velocity dependence on the capillary number are plotted. It is shown that for all cases of capillary numbers, the viscous barrier formed in the channel imparts additional acceleration to the droplet. For all calculations, the total normalized mass of the system was calculated over time, and the maximum error is shown to not exceed 5%. The droplet shape at the channel exit is shown for different cases of a viscous barrier. It is shown that if the capillary number is high, inclusions of the surrounding liquid can form in the droplet. However, if the droplet moves through a channel with a viscous barrier, the high-viscosity region prevents the formation of inclusions of other liquids due to more intense droplet deformation. It is established that the relative velocity of a droplet in a flow of an anomalously viscous liquid is higher than in a flow of liquid with a constant minimum viscosity. Graphs of the droplet’s relative velocity versus the liquid anomaly parameter are presented for three capillary numbers. It is shown that the higher the liquid anomaly parameter, the higher the relative velocity of the droplet. It is established that the viscous barrier formed in the channel accelerates the droplet and increases its relative velocity.

Keywords

phase field model,
Cahn–Hilliard–Navier–Stokes equation,
anomalously viscous liquid,
surface tension,
droplet motion,
laminar flow

Article outline

In recent years, the study of the properties of anomalously viscous liquids has attracted the interest of many scientists from various fields of science, from oil refining to medicine. Unlike ordinary liquids, the viscosity of which decreases with increasing temperature, the viscosity of anomalously viscous media changes non-monotonically: it can increase sharply and then decrease in a certain temperature range. Under certain temperature conditions, an interesting phenomenon occurs in such liquids - a viscous barrier. This is a localized region with high viscosity, which can affect the movement of the main liquid flow, up to blockage. In this paper, we study the effect of a viscous barrier on the velocity of a droplet moving in a flow of an anomalously viscous liquid in a flat channel. Numerical modeling is performed using the phase field method, consisting of a system of Cahn–Hilliard–Navier–Stokes equations, supplemented by an equation for the temperature field. The purpose of the modeling is to study the effect of a viscous barrier located in the center of the channel on the velocity and deformation of a droplet. Three values of the capillary number Ca = 0.1, Ca = 1, Ca = 10 and three parameters of the viscosity anomaly A = 0 (without a viscous barrier), A = 5 (a small viscous barrier), A = 10 (a pronounced viscous barrier) are considered.

Graphs of the droplet velocity dependence on the capillary number are compiled. It is shown that for all cases of capillary numbers, the viscous barrier formed in the channel imparts additional acceleration to the droplet. For all calculations, the total normalized mass of the system was calculated over time, it is shown that the maximum error does not exceed 5%. Images of the droplet shape at the channel outlet for different cases of a viscous barrier are given. It is shown that if the capillary number is large (Ca > 1), then inclusions of the surrounding liquid can form in the droplet, however, if the droplet moves in a channel with a viscous barrier, then the highly viscous region does not allow inclusions of another liquid to form in the droplet due to more intense deformation of the droplet. It is established that the relative velocity of a drop in a flow of an abnormally viscous liquid is higher than in a flow of liquid with a constant minimum viscosity. Graphs of the dependence of the relative velocity of a drop on the liquid anomaly parameter for three capillary numbers are presented. It is shown that in a channel without a viscous barrier, the relative velocity of a drop is the lowest for all capillary numbers. In a channel with the largest viscous barrier (A = 10), the relative velocity of a drop is the highest. It is evident that the larger the viscous barrier, the stronger the acceleration it imparts to a drop: a viscous barrier acts as a kind of ”viscous gun“, imparting acceleration to a drop and increasing its velocity. It is established that a viscous barrier formed in a channel imparts acceleration to a drop and increases the relative velocity of a drop.

References

  1. Laurano R, Boffito M. Thermosensitive micellar hydrogels as vehicles to deliver drugs with different wettability. Frontiers in Bioengineering and Biotechnology. 2020;8:00708. DOI: 10.3389/fbioe.2020.00708
  2. Домнина ЮМ, Жаворонок ЕС, Суслов ВВ, Решетняк ДВ, Кедик СА. Термообратимое повышение вязкости при нагревании водных растворов блок-сополимеров этилен- и пропиленоксида, не сопровождающееся образованием геля. Высокомолекулярные соединения (серия А). 2022;64(3):163–170.
    Domnina YM, Zhavoronok ES, Suslov VV, Reshetnyak DV, Kedik SA. Thermoreversible Growth Of Viscosity Upon Heating Of The Aqueous Solutions Of Ethylene Oxide-Propylene Oxide Block Copolymers Unaccompanied By Gelation. Polymer Science, Series A. 2022;64(3):163–170. (in Russian) DOI: 10.31857/S2308112022200019
  3. Киреев ВН, Урманчеев СФ. Течение жидкостей с температурной аномалией вязкости. Труды Института механики Уфимского научного центра РАН. 2003;3:232–245.
    Kireev VN, Urmancheev SF. Flow of liquids with temperature anomaly of viscosity. Proceedings of the Institute of Mechanics of Ufa Branch of RAS. 2003;3:232–245. (in Russian) DOI: 10.21662/uim2003.1.017
  4. Мухутдинова АА, Низамова АД, Киреев ВН, Урманчеев СФ. Экспериментальная установка для исследования устойчивости течения жидкости. Многофазные системы. 2024;19(1):35–39.
    Mukhutdinova AA, Nizamova AD, Kireev VN, Urmancheev SF. Experimental setup for studying the stability of fluid flow. Multiphase systems. 2024;19(1):35–39. (in Russian) DOI: 10.21662/mfs2024.1.005
  5. Nizamova AD, Kireev VN, Urmancheev SF. The Influence Of The Annular Gap Thickness On The Critical Reynolds Number During The Flow Of Thermoviscous Liquids. Lobachevskii Journal of Mathematics. 2024;45(5):2119–2127. DOI: 10.1134/S1995080224602315
  6. Badalassi VE, Ceniceros HD, Banerjee S. Computation of multiphase systems with phase field models. Journal of computational physics. 2003;190(2):371–397. DOI: 10.1016/S0021-9991(03)00280-8
  7. Галеева ДР, Киреев ВН, Галеева ГЯ, УУНиТ. Программа для расчета скорости капли в плоском канале методом фазового поля. Свидетельство о государственной регистрации программы для ЭВМ № 2025614539. 2025.
    Galeeva DR, Kireev VN, Galeeva GYa, UUS\&T. Program for calculating the velocity of a drop in a flat channel using the phase field method. Certificate of state registration of the computer program No. 2025614539. 2025. (in Russian) EDN: UPVEMH
  8. Галеева ДР, Киреев ВН, Ковалева ЛА, Мусин АА. Консервативная численная схема для решения уравнения Кана–Хилларда. Прикладная математика и механика. 2025;89(1):136–148.
    Galeeva DR, Kireev VN, Kovaleva LA, Musin AA. A Conservative Numerical Method For Solving The Cahn-Hilliard Equation. Applied Mathematics And Mechanics. 2025. 2025;89(1):136–148. (in Russian) DOI: 10.31857/S0032823525010101