CC BY
5Современные инновации, системы и технологии // Modern Innovations, Systems and Technologies
2023; 3(3) eISSN: 2782-2818 https://www.oajmist.com
УДК: 532.5 EDN: NHKFXR
DOI: https://doi.org/10.47813/2782-2818-2023-3-3-0501-0509
Study of pulsating flow of viscous two-phase media in a flat porous pipe considering particle gravity force
K. Z. Abidov, B. T. Ergashev, О. K. Salieva
Bukhara Engineering-Technological Institute, Bukhara, Uzbekistan.
Annotation. This paper examines the issues of unsteady pulsating motion of a viscous two-phase medium in a flat porous pipe taking into account the gravity of a solid particle in a two-dimensional formulation. Changes in the longitudinal velocity profiles of mixtures and the distribution of the volumetric content of the carrier and solid phases in a flat pipe, as well as the rate of liquid infiltration, the distribution of pressure and the volumetric content of the carrier and solid phases in pipes with porous walls have been studied.
Keywords: pulsating flow, viscous two-phase media, particle gravity force.
For citation: Abidov, K. Z., Ergashev, B. T., & Salieva О. K. (2023). Study of pulsating flow of viscous two-phase media in a flat porous pipe considering particle gravity force. Современные инновации, системы и технологии - Modern Innovations, Systems and Technologies, 3(3), 0501-0509. https://doi.org/10.47813/2782-2818-2023-3-3-0501-0509
INTRODUCTION
Recently, the interest of researchers in the issues of pulsating flow of single- and multiphase media has increased significantly. Scientific interest in this problem is due to the fact that pulsating flows of mixtures are becoming increasingly widely used in many branches of human activity. Similar processes are often observed in the food and chemical industries, in the processes of extraction, dissolution, drying, crystallization, as well as when cleaning the internal surfaces of equipment pipes from deposits and sediments.
In addition, in pulsating flows, at appropriate frequencies and amplitudes of fluid flow oscillations, the area of turbulence can be significantly moved away, hydraulic resistance can be reduced, the flow rate can be changed, the filtration of liquid through the permeable wall of the pipeline can be accelerated, and the processes of extraction of essential oils, medicines, synthetic materials, etc.
© K. Z. Abidov, B. T. Ergashev, О. K. Salieva, 2023
0501
Although a large number of scientific papers have been published on the hydrodynamics of pulsating flows, many issues have not yet been sufficiently studied. These include, for example, the unsteady movement of pulsating flows of multiphase media in pipes with porous walls. Individual works and theories for studying these processes are usually based on the use of simplified models.
THEORETICAL BASIS
Let us consider the unsteady two-dimensional pulsating motion of a viscous two-phase medium in a flat pipe with porous walls. Let us choose a Cartesian coordinate system so that the axis is directed along the bottom wall of the pipe. Let's place the origin of coordinates at the intersection of the bottom wall and the inlet section.
In the hydrodynamic formulation, we will study the flow process on the basis of a system of dimensionless two-dimensional unsteady equations of motion of two-phase media for a plane flow in an interpenetrating model taking into account gravity, the continuity equation for phases and the mixture concentration ratio [1-4]:
дщ 1 дщ дщ M 1 дР fi^T + f^iT^T + = ~ — fi7— +
д
L дх7
дх1
ai'i L дх2
1
+ —
a.j
4 д д ui f д дх1 д vi f 2 д д vi f 1 д L д х1 д vi f
3 L2 дх2 [fi дх2] Г дх1] 3 L д х2 Г дх1] [fi дх2]
+
м
+ — K(uZ-i-ui) + Fr2f1-,
(1)
д vi 1 д vi д vi M 1 д Р
fi^T + fuij+ fvi=--fi T^— +
дЬ L дх2 дх1 ai L дх2
1
+ —
a
1 д
L2 д х
fi
д vi
д х
2
4 д 3 дх1
д vi
fiд.7Ь
1 д + т-
L д х
д ui
fio71.
2 д 3 L д х1
fi
д u
д х
+
м
+ ^TiK(V3-i-Vi)+Fr2fi;
(2)
f + dfivi I 1 dfiui _ 0-
dt dx1 L dx2 '
(3)
f1 + f2 = 1; i = 1,2.
(4)
To solve the given system of equations, it is necessary to formulate the initial and boundary conditions. A flat pipe with porous walls may have either a top or bottom permeable wall. A schematic picture of such flat-porous pipes is shown in Figure 1.
Figure 1. Schematic picture of flow in a flat-permeable pipe
Let a flat pipe have a permeable upper wall. At the initial moment of time, the two-
phase mixture is at rest. Starting from t>0, pulsating motion occurs due to the periodicity of the
pressure gradient at the inlet of the pipe.
Then the initial and boundary conditions will take the following form:
at t = 0: ui = 0, vt = 0, ft = fi0 for 0 < x1 < 1, 0 < x1 < 1; at t > 0: x2 = 0: P = P0 + A cos ut, f1= fi0, vt = 0, ui = u01(A(ai) cos Mt + B(ai) sin Mt), u2 = Au1 for 0 < x1 < 1; at x1 = 0: du1/u2 = 0, vt = 0 for 0 < x1 < 1; at x2 = 0: U[ = 0, Vi = 0 for 0 < x2 < 1; at x1 = 1: Ui = 0, v1 = a(1 — f2\x=n)(.P — Pc) for 0 < x2 < 1;
If a flat pipe has a bottom permeable wall, then the condition for the permeability of the
liquid phase is set at x1=0 for 0<x2<1. Taking into account the process of blocking small holes
in the permeable wall, this condition will be written as
Vi = a(1 — f2\x=H)(P — Pc). For flat pipes having top and bottom porous walls, the permeability condition will exist
for both the top and bottom walls. Under such initial and boundary conditions, we solve the
system of equations (1) - (4) using the iterative method of alternating directions at ^=0.
It is necessary to obtain a qualitative picture of the flow, i.e. changes in the velocity field, phase concentration and pressure depending on ai, Fr2, p/[j,2, K, a.
RESULTS AND DISCUSSION
Let us consider the change in the longitudinal velocities of a two-phase flow in the section X2 = 0.5 of a flat pipe having upper and lower porous walls at a1 = 2.2, a2 = 0.84, Re1 = 459, Re2 = 66, F^ = \,^1/\i2 = 0.06, a = 10-3,K = 60.
The longitudinal components of the velocities of the first (solid line) and second (dotted) phases for rat=0 (1) and rat=rc/3 (2) are shown in Figure 2. In the initial section, the phase velocity profiles coincide, because in the calculations, X=1 was assumed. With distance from the inlet section, the more viscous phase begins to lag behind the less viscous one, and the profiles tend to be parabolic. A similar pattern was established for axisymmetric pipes, where at low pulsation frequencies the profiles tend to a parabolic shape. Numerical experiments have shown that at low pulsation frequencies, particles of the second phase begin to settle in the lower part of the pipe due to gravity. This process to some extent affects the symmetry of the phase velocity distribution. In the upper part of the flat pipe, the velocities p will be greater than at the corresponding points in the lower part of the pipe.
The longitudinal components of the velocities of the first (solid line) and second (dotted) phases for rat=0 (1) u rat=rc/3 (2) are shown in Figure 2. In the initial section, the phase velocity profiles coincide, because in the calculations, X=1 was assumed. With distance from the inlet section, the more viscous phase begins to lag behind the less viscous one, and the profiles tend to be parabolic. A similar pattern was established for axisymmetric pipes, where at low pulsation frequencies the profiles tend to be parabolic [5,6]. Numerical experiments have shown that at low pulsation frequencies, particles of the second phase begin to settle in the lower part of the pipe due to gravity. This process to some extent affects the symmetry of the phase velocity distribution. At the top of a flat pipe, particle velocities will be greater than at the corresponding points at the bottom of the pipe.
Figure 2. Longitudinal velocity profiles in a flat permeable pipe.
Figure 3. Change in longitudinal velocity profiles in a flat pipe with a permeable bottom wall.
Figure 3 shows the manifestation of the Richardson effect in a flat pipe with a permeable bottom wall. The solution was obtained with ai=6, a2=2.24, Fr2=l; a=10"3; for X2=0,5; rat=rc/12 (1), rat=0 (2). Comparing the curves, we see that with the formation of the Richardson effect, the maximum value of the relative difference u1-u2 moves from the axis closer to the wall [7,8].
The permeability of the wall reduces the total resistance of the pulsating flow. With an increase in the permeability coefficient, the pulsating flow can turn into non-pulsating laminar movement.
Therefore, if it is necessary to weaken the flow pulsation, then it is necessary to create an additional permeable inner wall of the pipe.
Numerical experiments have shown that to reduce the degree of siltation of a permeable pipe, the frequency of pulsation of the mixture along the flow should be increased.
In the near-wall region of the upper and lower porous walls, a decrease in the volumetric content of the second phase is observed. In this case, the rates of liquid penetration through the upper and lower pores of the wall are almost equal.
Figure 4 shows the change in the rate of liquid infiltration along the length of a flat pipe with a lower porous wall. Curve 1 was obtained at ai=8, a2=3.05, and curve 2 at ai=2,2; a2=0.84; with the same values of the other parameters, when Fr2=l; M1/M2 =0.062, K=50, a=10-3.
It can be seen that the rate of seepage along the length of the pipe decreases. At small values of a, particles of the solid phase begin to settle on the walls of the pipe and thereby intensify the process of clogging the pores of the porous wall. This pattern confirms the advisability of taking into account the dependence of the permeability coefficient on f2|x=0.
As the permeability coefficient increases, the internal pressure approaches the ambient pressure, and the seepage length decreases. Figure 5 shows the pressure distribution across the layers of a flat pipe with a permeable bottom wall. The solution was obtained with the same initial data for ai=2,2; a2=0,84; the numbers correspond to the pipe layer. In a permeable pipe, the pressure inside the pipe is less than in an impermeable pipe. This is due to the fact that in a permeable pipe the pressure loses its strength to overcome the resulting additional pressure gradient caused by the permeability of the wall. This effect is confirmed in other works [9, 10].
Let us consider the redistribution of concentration in a flat pipe with porous walls taking into account gravity. In Fig. Figure 6 shows graphs of the distribution of the volume concentration of the carrier phase in the section X2=0,6 of a flat pipe with a lower permeable wall. Curve 1 was obtained at ai=2,2; a2=0,84; and curve 2 at ai=8; a2=3,05; for Fr2=l; K=60; M1/M2 = 0.062; a=10"3; X2=0.6; rat=0. Permeability affects the flow concentration; the symmetry of the distribution is disrupted.
Figure 4. Variation in the rate of liquid infiltration along the length in a flat pipe with a
permeable bottom wall.
The movement of the mixture near the pipe wall is unstable with respect to the flow parameters. With increasing pulsation frequency (hence, with increasing a;), the volumetric content of the carrier phase at the permeable wall increases. From curve 2 it is clear that at low
pulsation frequencies, solid particles settle on the permeable wall and, thereby, displace liquid particles. This process leads to a decrease in the rate of liquid penetration through the pores of the wall and reduces the pressure gradient at the permeable pipe wall.
Figure 7 illustrates the change in the volume concentration of the more viscous phase in the initial sections of a flat pipe with upper and lower permeable walls. The numbers correspond to the pipe cross-section. The picture of the £2 distribution in sections of a flat porous pipe was obtained when the Richardson effect took place. The calculation was made for a1 = 8, a2 = 3.05, Re1 = 411, Re2 = 3.05, gnn Fr2 = 1, K = 60,^1/^2 = 0.062, a = 10-3, rat=0.. From the calculation results it is clear that from the initial cross sections f2 undergoes large fluctuations. As we move away from the inlet section, the two f2 peaks formed in the flow core gradually smooth out. In this case, the oscillatory nature of the concentration distribution is preserved.
Calculations have shown that at low pulsation values, the rate of liquid leakage through the upper porous wall is greater than through the lower one. This is due to the fact that near the upper permeable wall the concentration of solid particles is less than near the lower one.
0 0.66 0,70 0,74 f.
Figure 6. Distribution of the volumetric content of the carrier phase in a flat pipe with a
permeable bottom wall.
О 0,24 0,26 0.28 0,30 0,32 0.34 0,34 0,36 f2
Figure 7. Distribution of volume concentration of the solid phase in the initial sections of a
flat permeable pipe.
CONCLUSION
With a pulsating flow of two-phase media in a flat porous pipe, increasing the Womersley frequency parameter, the speed of particle transfer near the wall increases, and in the core of the flow there is practically no transfer of particles. A significant increase in the transfer of inclusion particles is observed during the flow of a less viscous phase. With low-frequency pulsations, a zone of highly concentrated mixture is formed in the near-pipe area. An increase in the pulsation frequency leads to a displacement of this zone in the direction of the axis and, thereby, forms a slightly concentrated mixture in the near-wall region.
The results obtained confirm the conclusion that it is necessary to take into account the influence of gravity of solid particles and the variability of the permeability coefficient on the nature of the distribution of phase concentrations, speed and pressure when the mixture moves in a flat porous pipe.
REFERENCES
[1] Abidov K.Z., Ergashev B.T. Hydrodynamic model for controlling the movement of a two-phase mixture during transportation in a pipe. Universum: technical sciences. 2021; 4(85): 72-77. https://doi.org/10.32743/UniTech.2021.85.4-1.72-76
[2] Mazitov A.A., Bobreneva Yu.O., Gubaidullin I.M. Mathematical modeling of multiphase flows using modern technologies. In: Supercomputer technologies of mathematical modeling: abstracts of the V international conference. Edited by V.I. Vasilyeva. Yakutsk: NEFU Publishing House, 2022. 26. https://doi .org/10.20948/prepr-2022-62
[3] Smirnov D.D., Bobreneva Yu.O., Mazitov A.A., Marchenko M.A., Gubaiddulin I.M., Chernykh I.G. Parallel algorithm for a numerical method for modeling mass transfer in a
fractured-pore reservoir on a supercomputer. Parallel Computational Technologies (PAVT 2020). 2020; 255-264.
[4] Umarov A.I., Akhmedov Sh.H. Two-dimensional problems of hydrodynamics of multiphase media. Tashkent: Fan. 1989. 94.
[5] Bobreneva Yu.O., Rahimly P.I., Poveshchenko Yu.A., Podryga V.O., Enikeeva L.V. On one method of numerical modeling of piezoconductive processes of a two-phase fluid system in a fractured-porous reservoir. Journal of Physics: Conference Series. 2021; 2131: 022001. https://doi.org/10.1088/1742-6596/2131/27022001
[6] Abidov K.Z., Ergashev B.T. Hydrodynamic model of motion control of pulsiting flow of two-phase media in a pipe. IOP Conference Series: Earth and Environmental Science. 2022; 2388: 012178. https://doi.org/10.1088/1742-6596/2388/1/012178
[7] Rihardson E.G., Tyler E. The transverse velocity gradient near the mouths of pipes in which an altemating or continuous flow of air is established. Proc.Phys. Soc. 1929; 42(1): 115. https://doi.org/10.1088/0959-5309/42/1/302
[8] Womersley I.R. Method for the calculation ob velocity rate ob flow and viscous drag In arteries when the pressure gradient is khowh. J. Physiol. 1955; 127(3): 533-563. https://doi.org/10.1113/jphysiol.1955.sp005276
[9] Mazitov A.A., et al. Mathematical modeling of multiphase flow in a single-pore reservoir. Preprints of IPM n.a. M.V.Keldysh. 2022; 62: 14. https://doi.org/10.20948/prepr-2022-62
[10] Bedarev I.A., Temerbekov V.M., Numerical modeling of sublimation of solid porous objects in an axisymmetric formulation. International seminar on dynamics of multiphase media. Berdsk, 2023. 5-6.
[11]
INFORMATION ABOUT THE AUTHORS
K. Z. Abidov, PhD, Associate Professor, Bakhtiyor Ergashev, Senior Lecturer at the
Bukhara Engineering-Technological Institute, Department of TZHBACT, Bukhara Bukhara, Uzbekistan Engineering-Technological Institute, Bukhara,
Uzbekistan
О. K. Salieva, Bukhara Engineering-Technological Institute, Bukhara, Uzbekistan
Статья поступила в редакцию 20.08.2023; одобрена после рецензирования 28.09.2023; принята
к публикации 29.09.2023.
The article was submitted 20.08.2023; approved after reviewing 28.09.2023; accepted for publication
29.09.2023.
