CC BY
1
Math-Net.Ru
N. V. БигшаэЬеуа, Е. А. Ьаппа, Е. Уи. Ргоэукуакоу, Неоднородное течение Куэтта двухслойной жидкости, Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2023, номер 3, 530-543
001: 10.14498/уБ^и1968
Использование Общероссийского математического портала Math-Net.Ru подразумевает, что вы прочитали и согласны с пользовательским соглашением
http://www.mathnet.ru/rus/agreement
Параметры загрузки:
IP: 109.252.33.182
29 сентября 2024 г., 12:13:51
Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki
[J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 3, pp. 530-543 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1968
MSC: 35C10, 76D05, 35G20
Inhomogeneous Couette flows for a two-layer fluid
N. V. Burmasheva1'2, E. A. Larina1'2, E. Yu. Prosviryakov1'2'3
1 Ural Federal University named after the first President of Russia B. N. Yeltsin, 19, Mira st., Ekaterinburg, 620002, Russian Federation.
2 Institute of Engineering Science, Ural Branch of RAS,
34, Komsomolskaya st., Ekaterinburg, 620049, Russian Federation.
3 Udmurt Federal Research Center, Ural Branch of RAS, 34, T. Baramzina st., Izhevsk, 426067, Russian Federation.
Abstract
The paper presents a new exact solution to the Navier-Stokes equations which describes a steady shearing isothermal flow of an incompressible two-layer fluid stratified in terms of density and/or viscosity, the vertical velocity of the fluid being zero. This exact solution belongs to the class of functions linear in terms of spatial coordinates, and it is an extension of the classical Couette flow in an extended horizontal layer to the case of non-one-dimensional non-uniform flows. The solution constructed for each layer is studied for the ability to describe the appearance of stagnation points in the velocity field and the generation of counterflows. It has been found that the flow of a two-layer fluid is stratified into two zones where the fluid flows in counter directions. It is also shown that the tangential stress tensor components can change their sign.
Keywords: stratified viscous fluid, exact solution, field stratification, coun-tercurrent.
Received: 13th December, 2022 / Revised: 28th June, 2023 / Accepted: 17th July, 2023 / First online: 25th September, 2023
Mathematical Modeling, Numerical Methods and Software Complexes Research Article
© Authors, 2023
© Samara State Technical University, 2023 (Compilation, Design, and Layout)
3 ©® The content is published under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/) Please cite this article in press as:
Burmasheva N. V., Larina E. A., Prosviryakov E. Yu. Inhomogeneous Couette flows for a two-layer fluid, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 3, pp. 530-543. EDN: YOWQMI. DOI: 10.14498/vsgtu1968. Authors' Details:
Natalya V. Burmasheva https://orcid.org/0000-0003-4711-1894
Cand. Tech. Sci.; Associate Professor; Dept. of Information Technology and Automation1; Senior Researcher; Sect. of Nonlinear Vortex Hydrodynamics2; e-mail: [email protected] Ekaterina A. Larina © https://orcid.org/0009-0000-7883-0803
Assistant; Dept. of Information Technology and Automation1; Engineer; Sect. of Nonlinear Vortex Hydrodynamics2; e-mail: [email protected] Evgeniy Yu. Prosviryakov © https://orcid.org/0000-0002-2349-7801
Dr. Phys. & Math. Sci.; Professor; Dept. of Information Technology and Automation1; Head of Sector; Sect. of Nonlinear Vortex Hydrodynamics2; Leading Researcher; Lab. of Physical and Chemical Mechanics3; e-mail: [email protected]
Introduction. The Couette flow is one of the first examples of an exact solution to the Navier-Stokes equations [1-3]. It describes the isobaric flow of a viscous incompressible fluid, which is induced due to the motion of one or both boundaries of an infinite horizontal fluid layer. Recall that the steady Couette flow is described by a linear velocity profile, and this has predetermined the popularity of this solution in the theory of hydrodynamic stability for studying secondary flows generated by different disturbance classes. There is an interesting observation for the unsteady Couette flow. It is described by the simplest linear parabolic equation having a general solution to an extensive class of functions [2,3]. Thus, the first and second Stokes problems are described by the non-stationary Couette profile, and they are its particular case [2,3].
Besides studying the hydrodynamic stability of the Couette flow for a viscous incompressible fluid, there are various modifications of the exact Couette solution for regions without plane symmetry. Note that there exists the well-known exact Taylor-Couette solution on the fluid flow in the gap between coaxial cylinders [2,4-8], as well as the solution describing isobaric fluid flow on a sphere [9,10]. The Couette flow is fundamental in the study of fluids with non-Newtonian properties [11,12]. The three-dimensional Couette flow, which is potential and non-isobaric, has been recently exemplified [13].
It is difficult to study Couette flows, different from unidirectional ones, since the reduced Navier-Stokes system becomes overdetermined [14-17]. The overde-termined Navier-Stokes equation system describing two-dimensional flows of viscous fluids began to be studied in [17], a complete list of exact solutions for two-dimensional hydrodynamics being given in [18]. An example of a nontrivial exact solution to the Navier-Stokes equation system for incompressible fluids with nonstationary and steady two-dimensional velocity fields depending on three coordinates was found in [15,16]. The first exact solution describing the non-uniform Couette flow was constructed in the class of solutions for velocities linearly dependent on two coordinates (the Lin-Sidorov-Aristov family) [19-21]. The studies along this line were continued and summarized in [16].
This paper studies a boundary value problem for 2.5D Navier-Stokes equations, which describes steady flows of a two-layer fluid. The study is based on the exact solution of the Navier-Stokes equations for incompressible fluids [15,16], which was constructed by functional variable separation. It was shown in [22] that the exact solutions found in [15,16] and describing non-uniform Couette flows can be used to describe isobaric multilayer fluids. It was reported in [15,16] that exact solutions with a velocity field linear in coordinates describe equatorial countercur-rents in the World Ocean [2,15,21]. Multilayer fluids are often used to model large oceanic flows [2,21]. In [22] it was found useful to extend the study of boundary value problems of steady flows from single-layer streams [15,16] to multilayer fluids. This should be done in order to have a store of exact solutions for studying the hydrodynamic stability of flows, for comparing model representations with natural observations, and most importantly to understand the applicability of the substitution of continuous density stratification by a discrete saltus function.
1. Problem Statement. Consider a steady flow of a viscous two-layer fluid in an extended horizontal layer. We denote the lower layer by the subscript "1" and the upper one by "2". Each of the two layers of the two-layer fluid can have its own thickness (hi and h2, respectively, see Fig. 1), density (p\ and p2), and viscosity
z, h
h2 hi
hy------------------- / ' / 1 / 1
Q / / / /
/ "y
Figure 1. The scheme of a two-layer fluid
(^i and ^2). Note that the heavier phase is located below in density-stratified fluids, i.e., p1 > p2.
It is assumed here that the flow takes place at a constant temperature and in the absence of external forces, except for gravity. The sum of pressure (ratioed to density) and the gravity potential is constant through the flow; therefore, the gradient of this sum, entering the Navier-Stokes equations, is zero. The steady flow of fluids of this type is described by two systems of nonlinear equations (each layer is described by its own system) consisting of the vector Navier-Stokes equation and the scalar incompressibility equation [22-24]:
Pi(V(1), V)V(1) = ^AV(1), V ■ V(1) = 0; (1)
P2(V(2), V)V(2) = ^AV(2), V-V(2) =0. (2)
Systems (1), (2) have the following notations: V(1) = (V^1),Vy(1),Vz(1)) and
v(2) = (v(2),vy(2),vz(2)) are the vector velocity fields for the lower and upper layers, respectively; V = (d/dx,d/dy,d/dz) is the Hamilton operator; A = (d2/dx2+d2/dy2 +d2/dz2) is the Laplace operator; (V, V) = (Vxd/dx+Vyd/dy+
Vzd/dz) is a convective derivative. Note that both system (1) and system (2) are
(j) (j)
overdetermined, i.e., the required velocities Vx , Vy for each layer must satisfy three scalar equations (the third Navier-Stokes equation is fulfilled identically since a flow with zero vertical velocity is considered). Note that the equations V ■ V(1) = 0 and V ■ V(2) = 0 (velocity divergence is zero) is termed in two ways in the scientific literature. They are termed the continuity equation in the physical literature and the incompressibility equation in the hydrodynamic literature [25,26].
In what follows, the solution of systems (1), (2) is sought in the form of functions linearly dependent on one of the horizontal coordinates [22-24]:
v;(1) = u(1)(z) + u(1)(z)y, Vy(1) = V (1)(z); (3)
v;(2) = u(2)(z) + u(2)(z)y, Vy(2) = V(2) (z). (4)
For the class (3), (4) the incompressibility equation in systems (1), (2) is satisfied identically. This circumstance dismantles the problem of overdetermination
of systems (1) and (2), i.e., each system is now reduced to finding two projections of the velocity vector from two ordinary differential equations (projections of the Navier-Stokes equation onto the axes Ox and Oy):
P1V (1)u(1) = ^(u(1)" + u(1)"y), mV(1)" = 0; (5)
PlV (2)u(2) = m(u(2)" + u(2)"y), V2V(2)" = 0. (6)
Hereinafter, the double prime marks derivation with respect to the vertical coordinate z. In view of the independence of the spatial coordinates x and y of the selected Cartesian system, equations (5) and (6) can be represented as
u(1)" = 0, V(1)" = 0, u(1)" = ^V (1)u(1); (7)
1
u(2)" = 0, V(2)" = 0, u(2)" = —V (2)u(2). (8)
2
The first two equations in both system (7) and system (8) are isolated, and the solution of the third equations in these systems is the last to be found. Double integration of systems (7) and (8) results in their general solution
u(1) = C2(1)z + C1(1), V(1) + «1(1),
(1) = J1_ 12 r]1
+ ß2(1)z + A(1); (9)
U- = (C2(1W + 2(C2(1)«1(1) + C1(1)Ö2(1))^3 + 6C1(1W) +
u(2) = C2(2)z + C1(2), V(2) = + «1(2),
U(2) = (C2(2)a2(2)Z4 + 2(C2(2)«1(2) + C1(2)«2(2))^3 + 6Cl(2)«1(2)Z2) +
12 r]2
+ &(2)* + A(2). (10)
The solutions represented by equations (9) and (10) are polynomial, the highest degree of these polynomials corresponds to expressions for velocities U(1), U(2), and this is attributable to the sequence of integration of the equations in systems (7) and (8). The constants c(1), c(2), a(1), a(2), and ^(2) (i = 1,2) in the exponential solutions (9) and (10) must be found from the boundary conditions; therefore, it is necessary to formulate twelve conditions for the determination of these values.
2. Boundary Conditions. Since the vertical fluid velocity is assumed to be zero, the fluids of the different layers do not intermix in the shear flow under study. In other words, the interlayer boundary (the boundary z = h1) is here considered to be rigid (Fig. 1). For convenience, in what follows, h = h1 + h2.
Assume that a no-slip condition is set at the lower boundary z = 0 [15]:
Vi1)(0) = 0, v„(1) (0) = 0.
Taking into account the representation (3), (4), we arrive at three conditions U (1)(0) = 0, u(1)(0) = 0, V (1)(0) = 0. (11)
A velocity field (wind effect) is set at the upper boundary z = h [15]:
V(2) (h) = W cos tp + Qy, Vjj2) (h) = W sin p.
(2)
Here, Q is the horizontal gradient of the velocity Vx ) (spatial acceleration) at the upper boundary; W is the absolute value of the uniform velocity component at
the mobile boundary z = h of a two-layer fluid; <p is the angle between the uniform
(2)
velocity component Vx ) and the axis Ox. Taking into account representation (4), we obtain three more conditions
U (2)(h) = W cos <p, u(2)(h)=tt, F(2) (h) = W sin
(12)
Besides, it is required that two additional conditions be met at the interlayer boundary 2 = hi (conditions for "sewing together" the solutions for the two layers). They are
- the solution continuity condition for velocities
Fj1)(hi) = Fj2)(hi), ^(hi) = V^2)(hi ),
and taking into account the structure of classes (3) and (4), we have three equalities
U (i)(hi) = U (2)(hi), u(i)(hi)= u(2) (hi), V (i)(hi ) = V (2)(hi); (13)
- the solution continuity condition for tangential stresses
)(hi)= 42 )(hi), r«(hi)= r|)(hi).
As distinct from expressions (13) resulting (in view of the independence of the spatial coordinates) directly from the condition of equality of the velocities at the interlayer boundary, the case with the continuity condition for tangential stresses is not as apparent, the continuity condition for the velocity field proves to be insufficient. By the Newton law, the relation of the stress tensor components to velocities is known,
- P(l) + Vi (dVxW , dVy(i K V dy dx ) Vi (dVx(i I dz d x
r« = m (dvj* , dVy (• K V dy 1 dx J V(0 +2r,dVyW (dVy ( I dz + dVz (i) \ 1 dy J
\n% (dVj^ , dVz (l h V dz dy / m (dVz(i) | dVy (• K V dy dz J p(o + y
/ -p(i)
•q.i u(t)
, dVx(€> % dz
-P(i)
dVy w ^ dz
(3),(4)
n-dVx ^ 'h dz, dVy W
Vi dz -p(i)
where is hydrostatic pressure (the pressure of a moving fluid column, depending not only on the transverse coordinate, i.e. changing only with depth, similarly to the main equation of hydrostatics) in the layer under study. This representation of the stress tensor can easily yield the following equivalent of the continuity condition for tangential stresses:
du(1)
du(2)
= ^2" ,
z=h! dZ
d V(1)
d u(1)
.z=hi dZ
dV(2) V2-
du(2)
= ^2" ,
z=hi dZ
z=hi dz
z=hi
z=hi
(14)
Thus the boundary value problem (9)-(14) becomes closed.
3. An Exact Solution to the Boundary Value Problem. Conditions (11)-(14) allow us to find a particular solution to systems (7), (8) for each layer, which would meet the selected boundary conditions:
u°(Z) = h2 ri^hm2, (15)
^)=SWh-l 2, (16)
u(1)(Z) =------3 [QP1^(h2i?1 + h^h3 sin(Z3+
12 r?1(h2??1 + h^)
+12r]1ï]2(h2r]1 + h^)2 cos ( - Q sin ((4h{h2rnri^p1 + h\rp3p1+
+h2 rj3P2 + 4h1h2r?2r?2P2 + 6h1h2i?1i?2p2)] , (17)
u(2)(2) = l^Th^ ^ + ^ - (18)
V(2)(2) = + ^ - *)), (19)
u(2)(Z) =- w --3 {h47]2^(h2V1 + hm)p2Qsiny24-
12 f]2 (h2V1 + h1î]2)
-4h3h1ip(f]1 - T]2)(h2ij1 + h1ï]2)p2Qsm(pZ3+ +6h2h1(r]1 - r]2)2(h2T]1 + h1ï]2) P2Q sin (Z2--hZ[(-3h1r?3p1 + (h2(4h3 + h2)r?3 + 4h1(h3 - 3h2h2 + h2)^2+ +6h2(-2h1 + 2h1h2 + h2)^1??2 + 12h1 Jj|)p2^ siny-
-12r]1Ï]2(h2lj1 + h1^2)2 cos ( + -h1(h1 + h2)Wsin((3h1r?3p1 - hfh2r]1p2 + h1h2r?3p2 - h2r?3p2--h\r}lr}2p2 + 5h1h2r?|r?2p2(h1 - h2) + h3^2^?2p2 + 4h31rnr]'2p2--10h21h2r]1rl^p2 + 4h1h2^1i?ip2 + 6h2r?3p2(h2 - h^)
-12h1(TJ1 - ÏJ2)m(h2V1 + h^)2 cosy}. (20)
535
Solutions (15)-(20) involve the substitution Z = z/h e [0,1]. The motion of the lower layer (Z e [0; h1 /h]) is described by expressions (15)-(17). The values of Z e [h1/h;1] correspond to the other (upper) layer, where the flow velocity is determined by solutios (18)-(20).
It is of interest that, despite seemingly simple structure of boundary conditions (11), (12) and velocity field representation (11), (12), the obtained solution (15)-(20) depends on the parameters of the boundary problem under study and the physical characteristics of the fluid in an extremely non-trivial way. Therefore, in a general form, it is impossible to make any conclusions about the effect of a specific parameter on the structure of the final solution. First of all, this concerns the uniform components U(1) and U(2). The only conclusion (fairly obvious) that can be made concerns the effect of the values of W and Q on the value of velocity in the Oy direction (the velocities v (1) and V(2)) and the non-uniform velocity component along the Ox axis (the spatial gradients -u(1) and u(2)) since functions (15), (16) and (18), (19) depend in direct proportion on these parameters.
Note also that, when W = 0, the velocity field determined by solutions (15)-(17) and (18)-(20) assumes a trivial form; therefore, it is considered hereinafter that W = 0.
Besides, note that the functions -u(1) and u(2) become zero when Q = 0, i. e., we obtain an extension of the classical Couette flow to a non-one-dimensional case. The velocities V(1) and V(2) become zero only if sin <p = 0. In this case, there is a unidirectional flow along the Ox axis with a non-uniform velocity distribution. In addition, the exact solutions for the components -u(1) and V(1) are linear functions of Z, whose zero is only the point Z = 0. The stagnation (zero) points of the fluid are of particular importance due to the fact that both simply flow suppression and flow reversal (the appearance of counterflows) are possible at these points.
4. Analysis of the Exact Solution for the Lower Layer. Expression (17) for velocity U(1) can be written as
u(1) = [aZ 3 + b], (21)
hWZ 12^1 (h2^1 + h^)
where the following nomenclature is introduced:
a = Qp1T]'^(h2^1 + h1^2)h3 sin p, b = 12rp^(h2 r/1 + h1i?2 )2 cos p-
Q Q O O O Q OO ON
-Q sin tp(4h1h2i?1i?2P1 + h1%p1 + h2r]lp2 + ^2P2 + 6h1h^rni]^p2).
The case a = 0 in Eq. (21) is not discussed here due to its triviality. Assume further that a = 0, then the velocity u (1) becomes zero inside the layer under study only if the point Z0 = — tfbja falls within this layer. The latter condition is equivalent to the fulfillment of the inequality
6(ah? + bh3) < 0.
The corresponding velocity profiles U(1) are shown in Figs. 2 and 3. These figures show the value of angular velocity, corresponding to the first Coriolis parameter, for the equatorial zone of the World Ocean (the latitude is ^/12 rad).
Figure 2. The U-velocity profiles (the dashed line Ufor the lower layer and the solid line U(2) for the upper layer) for W = 10 m/s, ip = ^/12, h1 = 0.4 m, h2 = 0.6 m, O = 1.4584• 10-5 rad/s, pi = 1100 kg/m3, p2 = 1057.6 kg/m3, -qi = 1.7 • 10-5 Pa • s, ^2 = 1.2 • 10-5 Pa • s z
Figure 3. The U-velocity profiles (the dashed line U(1) for the lower layer and the solid line U(2) for the upper layer) for W = 10 m/s, ip = ^/12, h1 = 0.6 m, h2 = 0.4 m, O = 1.4584• 10-5 rad/s, pi = 1100 kg/m3, p2 = 733.1 kg/m3, = 1.7 • 10-5 Pa • s, ^2 = 1.2 • 10-5 Pa • s
In view of the above notations, the tangential stress is defined as
r£> =-W-3 [4aZ3 + 6] T m ^ .
xz 12(h2m + h^) + h1m
Note that, when W = 0, the tangential stress Tx1 assumes a constant value (different for each section y). Therefore, it is further assumed that W = 0, and
hence the stress Tx1 can be represented as
r« =-hW-3 [4 a Z3 + 6] + h7]l ^
xz 12(h2V 1+ h1?72) h2m + h1^2
h w
h W [4 a Z3 + C], (22)
12(h2^1 + h№)3
where C denotes the following expression:
12 r]1î]2(h2V1 + h1^2 )2Qy
C=b+ W .
Note that the structure of expression (22) is similar to formula (21). Hence, by analogy with the analysis of the expression (21), it follows that, when a = 0, the
stress Tx1 becomes zero only at one point inside the layer, namely at the point Z = - ^Cj(Aa). This value belongs to the interval of interest only if
C(4ah{ + Ch3) < 0.
Thus, the profile of the tangential stress represented by expression (22) is defined by a cubic parabola, which cannot have more than one point in common with the OZ axis (Figs. 4 and 5). The profiles of the tangential stress rxz in Figs. 4 and 5 are constructed for those values of the boundary conditions, layer thicknesses, and physical fluid parameters for which a calculation is made and presented in Figs. 2 and 3, respectively.
z 1.0
0.8
0.6
o./j
6.2
-0.01 0 0.01 0.02 0.03 0.04 0.05 0.06 txz, kg/(m ■ s2)
Figure 4. The profile of the stress txz (the dashed line ri^ for the lower layer and the solid line t^2 for the upper layer) for W = 10 m/s, ip = ^/12, hi = 0.4 m, h2 = 0.6 m, O = 1.4584-10-6 rad/s, pi = 1100 kg/m3, p2 = 1057.6 kg/m3, -qi = 1.7-10-5 Pa-s, ^ = 1.2-10-5 Pa-s
z 1.0
0.8
0.6
0.4
0.2
0 0.01 0.02 0.03 0.04
0.05 txz, kg/(m • s2)
(the dashed line ri^ for the lower layer and the solid
Figure 5. The profile of the stress tx line riJ' for the upper layer) for W = 10 m/s, ip = ^/12, hi = 0.6 m, h2 = 0.4 m, O = 1.4584• 10-5 rad/s, pi = 900 kg/m3, p2 = 733.1 kg/m3, ^ = 1.7• 10-5 Pa• s, = 1.2• 10-5 Pa• s
Note that, in view of the structure of solutions(16) and (19), the tangential stress ryz assumes a constant value determined by the problem parameters:
dV(1) r/idV(1) r)ir)2W sin^> (2) yz = m dz = ~h dZ = h22m+h^]2 = Tyz '
5. Analysis of the Exact Solution for the Upper Layer. The expression^) for the velocity U(2) can be represented as
u(2) = 10 W ,-^ [kZ4 + nZ3 + mZ2 + qZ + r].
I2r]2(h2m + hiT]2)3
The coefficients k, n, m, q and r can be easily written from exact solution (20).
Let us now study the properties of the tangential stress defined by the velocity Vx(2). Using the above notations, we obtain
tX2 = £ {-W-_ [4kZ 3 + 3nZ 2 + 2mZ + q] + ^ ).
h 112r]2(h2m + h^2) h2V1 + hn^ J
(2)
Note that, when W = 0, the stress tX2) assumes a constant (through the selected section) value, hence the stratification of this field does not occur. In the
(2)
assumption that W = 0, the stress tXz can be represented as W
X2 = ?-A 4k Z3 + 3nZ2 + 2mZ + q+
x h 12 V2(h2V1 + h^)3 1
+ 12hr]1i]2(h2V1 + h^) Qy\
+ w y
and this illustrates stratification on the change in the sign of the tangential stress
(2>
relative to the Z-coordinate. The examples of the tangential stress Txz corresponding to this solution are given in Figs. 4 and 5.
Conclusion. The problem describing the isothermal flow of a viscous incompressible two-layer fluid in a horizontal layer has been studied. The properties of the layers differ in thickness, density and/or viscosity. An exact solution has been obtained to describe velocities in each of the layers for a set of boundary conditions describing no-slip of the fluid at the lower boundary of the fluid flow region under study and the non-uniform effect of wind at the upper boundary of this region. At the boundary between the layers of the two-layer fluid, it was required that the velocities and stresses be equal. The analysis of the obtained solution for the layers has shown that this exact solution is able to describe the appearance of reverse flow zones and stratification of the tangential stress field.
Competing interests. Authors declare that they have no competing interests. Authors' contributions and responsibilities. Each author has participated in the development of the concept of the article and in the writing of the manuscript. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of the manuscript. Funding. The study received no funding.
References
1. Couette M. Etudes sur le frottement des liquids, Ann. de Chim. et Phys. (6), 1890, vol. 21, pp. 433-510 (In French).
2. Drazin P., Riley N. The Navier-Stokes equations: A classification of flows and exact solutions, London Mathematical Society Lecture Note Serie, vol. 334. Cambridge, Cambridge Univ. Press, 2006, x+196 pp. DOI: https://doi.org/10.1017/CB09780511526459.
3. Aristov S. N., Knyazev D. V., Polyanin A. D. Exact solutions of the Navier-Stokes equations with the linear dependence of velocity components on two space variables, Theor. Found. Chem. Eng., 2009, vol.43, no. 5, pp. 642-662. EDN: LPGJRJ. DOI: https://doi.org/10.1134/ S0040579509050066.
4. Taylor G. I. Stability of a viscous liquid contained between two rotating cylinders, J. Phil. Trans. Royal Society A, 1923, vol.102, no. 718, pp. 541-542. DOI: https://doi.org/10. 1098/rsta.1923.0013.
5. Shlikhting G. Teoriia pogranichnogo sloia [Boundary-Layer Theory]. Moscow, Nauka, 1974, 712 pp. (In Russian)
6. Ilin K. I., Morgulis A. B. Critical curves for the Couette-Taylor throughflow, Izv. Vuzov. Severo-Kavkazskii Reg. Natural Sci., 2019, vol. 1, pp. 10-16 (In Russian). EDN: ZBZEWT.
7. Devisilov V. A., Sharay E. Yu. Hydrodynamic filtration, Safety in Technosphere, 2015, vol. 4, no. 3, pp. 68-80 (In Russian). EDN: TZKMCL. DOI: https://doi.org/10.12737/11885.
8. Lebiga V. A., Zinoviev V. N., Pak A. Yu., Zharov I. R. The circular gap Couette flow modeling, Vestn. Novosib. Gos. Univ. Ser. Fizika, 2016, vol. 11, no. 4, pp. 52-60 (In Russian). EDN: XRUKWH.
9. Astaf'ev N. M. Structures formed in a rotating spherical layer under the influence of conditions simulating global heat fluxes in the atmosphere, Sovr. Probl. DZZ Kosm., 2006, vol. 3, no. 1, pp. 245-256 (In Russian). EDN: NDPOWR.
10. Belyaev Yu. N., Monakhov A. A., Yavorskaya I. M. Stability of spherical Couette flow in thick layers when the inner sphere revolves, Fluid. Dyn., 1978, vol.13, no. 2, pp. 162-168. DOI: https://doi.org/10.1007/BF01091664.
11. Puhnachev V. V., Puhnacheva T. P. The Couette problem for a Kelvin-Voigt medium, J. Math. Sci., 2012, vol.186, no. 3, pp. 495-510. EDN: RGHYXH. DOI: https://doi.org/10. 1007/s10958-012-1003-0.
12. Skul'skiy O. I., Aristov S. N. Mekhanika anomal'no viazkikh zhidkostei [Mechanics of Abnormally Viscous Fluids]. Moscow, Izhevsk, R&C Dynamics, 2003, 156 pp. (In Russian). EDN: OWGAFL.
13. Prosviryakov E. Yu. Exact solutions for three-dimensional potential and vorticity Couette flows of an incompressible viscous fluid, Vestnik natsional'nogo issledovatel'skogo yadernogo universiteta "MIFI", 2015, vol.4, no. 6, pp. 501-506. EDN: VLRBSB. DOI: https://doi.org/ 10.1134/S2304487X15060127.
14. Aristov S. N., Prosviryakov E. Yu. On laminar flows of planar free convection, Nelin. Dinam., 2013, vol.9, no. 4, pp. 651-657 (In Russian). EDN: SAHAZJ.
15. Aristov S. N., Prosviryakov E. Yu. Inhomogeneous Couette flow, Nelin. Dinam., 2014, vol. 10, no. 2, pp. 177-182 (In Russian). EDN: SHVJGX.
16. Zubarev N. M., Prosviryakov E. Yu. Exact solutions for layered three-dimensional nonstationary isobaric flows of a viscous incompressible fluid, J. Appl. Mech. Tech. Phys., 2019, vol.60, no. 6, pp. 1031-1037. EDN: VKJMMQ. DOI: https://doi.org/10.1134/ S0021894419060075.
17. Berker R. Sur quelques cas d'intégration des équations du mouvement d'un fluide visqueux incompressible. Paris-Lille, Imprimerie A. Taffin-Lefort, 1936, 161 pp. (In French)
18. Shmyglevskii Yu. D. On isobaric planar flows of a viscous incompressible liquid, US'SR Comput. Math. Math. Phys., 1985, vol. 25, no. 6, pp. 191-193. DOI: https://doi.org/10. 1016/0041-5553(85)90030-8.
19. Lin C. C. Note on a class of exact solutions in magneto-hydrodynamics, Arch. Rational Mech. Anal., 1957, vol. 1, no. 1, pp. 391-395. DOI: https://doi.org/10.1007/BF00298016.
20. Sidorov A. F. Two classes of solutions of the fluid and gas mechanics equations and their connection to traveling wave theory, J. Appl. Mech. Tech. Phys., 1989, vol. 30, no. 2, pp. 197203. EDN: LZHPWE. DOI: https://doi.org/10.1007/BF00852164.
21. Aristov S. N. Eddy Currents in Thin Liquid Layers, Dr. Sci. [Phys.-Math.] Dissertation. Vladivostok, 1990, 303 pp. (In Russian)
22. Burmasheva N. V., Prosviryakov E. Yu. Exact solutions to the Navier-Stokes equations describing stratified fluid flows, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2021, vol.25, no. 3, pp. 491-507. EDN: JKXFDQ. DOI: https://doi.org/10.14498/vsgtu1860.
23. Burmasheva N. V., Prosviryakov E. Yu. A large-scale layered stationary convection of a incompressible viscous fluid under the action of shear stresses at the upper boundary. Temperature and presure field investigation, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2017, vol. 21, no. 4, pp. 736-751 (In Russian). EDN: YUGZXQ. DOI: https://doi.org/10.14498/vsgtu1568.
24. Burmasheva N. V., Prosviryakov E. Yu. Studying the stratification of hydrodynamic fields for laminar flows of vertically swirling fluids, Diagnostics, Resource and Mechanics of Materials and Structures, 2020, no. 4, pp. 62-78. DOI: https://doi.org/10.17804/2410-9908. 2020.4.062-078.
25. Landau L. D., Lifshitz E. M. Fluid Mechanics, Course of Theoretical Physics, vol. 6. Oxford, Pergamon Press, 1963, xii+536 pp.
26. Kochin N. K., Kibel I. A., Roze N. V. Theoretical Hydromechanics. New York, John Wiley and Sons, 1964, v+577 pp.
Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2023. Т. 27, № 3. С. 530-543 ISSN: 2310-7081 (online), 1991-8615 (print) d https://doi.org/10.14498/vsgtu1968
EDN: YOWQMI
УДК 532.5
Неоднородное течение Куэтта двухслойной жидкости
Н. В. Бурмашева1'2, Е. А. Ларина12, Е. Ю. Просвиряков1'2'3
1 Уральский федеральный университет
имени первого Президента России Б. Н. Ельцина, Россия, 620002, Екатеринбург, ул. Мира, 19.
2 Институт машиноведения УрО РАН,
Россия, 620049, Екатеринбург, ул. Комсомольская, 34.
3 Удмуртский федеральный исследовательский центр УрО РАН, Россия, 426067, Ижевск, ул. Т. Барамзиной, 34.
Аннотация
Предложено новое точное решение уравнений Навье-Стокса, описывающее установившееся изобарическое изотермическое течение стратифицированной по плотности и/или вязкости несжимаемой двуслойной жидкости. Указанное точное решение принадлежит классу функций, линейных по части пространственных координат, и является обобщением классического течения Куэтта в протяженном горизонтальном слое на случай неодномерных неоднородных течений. В качестве системы краевых условий рассмотрена связка «условие прилипания + воздействие параболического ветра». На общей границе двух слоев заявлено выполнение требования гладкости и непрерывности решения. Построенное для каждого слоя решение было исследовано на предмет возможности описывать возникновение застойных точек поля скорости и генерации противотечений. Строго показано, что указанное решение при
Математическое моделирование, численные методы и комплексы программ Научная статья
© Коллектив авторов, 2023 © СамГТУ, 2023 (составление, дизайн, макет)
3 ©® Контент публикуется на условиях лицензии Creative Commons Attribution 4.0 International (https://creativecommons.org/licenses/by/4.0/deed.ru)
Образец для цитирования
Burmasheva N. V., Larina E. A., Prosviryakov E. Yu. Inhomogeneous Couette flows for a two-layer fluid, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. Math. Sci.], 2023, vol. 27, no. 3, pp. 530-543. EDN: YOWQMI. DOI: 10.14498/vsgtu1968.
Сведения об авторах
Наталья Владимировна Бурмашева Ä https://orcid.org/0000-0003-4711-1894 кандидат технических наук; доцент; департамент информационных технологий и автома-тики1; старший научный сотрудник; сектор нелинейной вихревой гидродинамики2; e-mail: [email protected]
Екатерина Александровна Ларина © https://orcid.org/0009-0000-7883-0803 ассистент; департамент информационных технологий и автоматики1; инженер; сектор нелинейной вихревой гидродинамики2; e-mail: [email protected] Евгений Юрьевич Просвиряков © https://orcid.org/0000-0002-2349-7801 доктор физико-математических наук; профессор; департамент информационных технологий и автоматики1; заведующий сектором; сектор нелинейной вихревой гидродинами-ки2; ведущий научный сотрудник; лаб. физико-химической механики3; e-mail: [email protected]
определенном граничном управлении и варьировании геометрико-физи-ческих характеристик слоя отвечает множественной стратификации как поля скорости, так и порождаемого им поля касательных напряжений.
Ключевые слова: многослойная вязкая жидкость, точное решение, стратификация полей, противотечение.
Получение: 13 декабря 2022 г. / Исправление: 28 июня 2023 г. / Принятие: 17 июля 2023 г. / Публикация онлайн: 25 сентября 2023 г.
Конкурирующие интересы. Конфликта интересов в отношении авторства и публикации этой статьи нет.
Авторский вклад и ответственность. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами.
Финансирование. Исследование выполнялось без финансирования.
