CC BY
6
УДК 536.25
Modeling of Two-layer Fluid Flows with Evaporation at the Interface in the Presence of the Anomalous Thermocapillary Effect
Olga N. Goncharova* Ekaterina V. Rezanova^
Altai State University Lenina, 61, Barnaul, 656049
Russia
Received 02.10.2015, received in revised form 16.12.2015, accepted 20.01.2016 Stationary convective flows of two immiscible viscous incompressible fluids (liquid and gas) under action of the transverse gravity field and longitudinal temperature gradient along the interface are studied analytically. Mathematical model of the fluid flows with the effects of evaporation at the interface is based on exact solutions to the Navier-Stokes equations in the Oberbeck-Boussinesq approximation. The effects of the thermodiffusion and diffusive heat conductivity in the gas-vapor layer are taken into consideration. The obtained solutions are used to model the flows in the two-layer gas-liquid system in the case when a liquid exhibits the anomalous thermocapillary effect. Examples of the two-layer fluid flows are presented for various values of the gas flow rate, longitudinal temperature gradient at the interface and the gravity force acceleration.
Keywords: mathematical model, interface, evaporation, exact solution, anomalous thermocapillary effect. DOI: 10.17516/1997-1397-2016-9-1-48-59.
Introduction
Stationary two-layer fluid flows are studied on the basis of exact solutions to the Navier-Stokes equations in the Oberbeck-Boussinesq approximation in two-dimensional case [1-3]. The solutions are obtained in the case when only the longitudinal velocity is not equal to zero and it depends on the transverse coordinate. The temperature distribution and pressure in the upper and lower layers have the components which linearly depend on the longitudinal coordinate. The longitudinal gradients of temperature and vapor concentration depend linearly on the transverse coordinate. The longitudinal pressure gradient is a quadratic function with respect to the longitudinal coordinate. One of the first results devoted to construction of an exact solution in the infinite layer with evaporation at the "liquid-liquid" interface were obtained in [4]. The thermo-capillarity effect at the interface has not been taken into consideration. The solution has been used to model the binary flows with diffusion of the light impurity at interface [4].
The mentioned here solutions (see [1,2,4]) are solutions of the Birikh-Ostroumov type [5]. The solution of Birikh [5] has been constructed to describe the stationary convection in an infinite horizontal strip with solid boundaries or with a non-deformable free boundary under the action
* [email protected] [email protected] © Siberian Federal University. All rights reserved
of the gravity field and constant longitudinal temperature gradient (see also [6-8]). For the first time the similar problem has been studied in [9].
In the present paper the Soret and Dufour effects [10,11] are taken into consideration in the gas-vapor layer (see [3]). The exact solutions allow us to analyze the problem statements and peculiarities of modeling of flow effects including the effect of anomalous thermocapillary which attracts attention due to its importance for many liquids [12-16]. One of the main questions in modeling of fluid flows in channels is the formulation of a condition for vapor concentration on the upper solid wall of the channel. In the present paper the condition of the absence of the vapor flux is assumed to be fulfilled on this boundary of the flow domain. It is also shown that the constructed solutions can be used to model flows in the two-layer system when liquid exhibits the anomalous thermocapillary effect.
1. Governing equation
The coordinate system (x,y) is chosen so that the vector of the gravity acceleration g is directed opposite to the Oy axis (g = (0, -g)). In order to formulate our problem in dimensionless form we introduce characteristic values of the problem as follows: the thickness of the liquid layer l is the characteristic length, u* = vi/l is the characteristic velocity, T* is the characteristic temperature drop, p* = piv'[/l2 is the characteristic pressure; here pi and vi are the density and the coefficient of the kinematic viscosity of the liquid, respectively. The system of two viscous incompressible fluids (liquid and gas-vapor mixture) fills the infinite horizontal layers with the upper and lower boundaries y = h, y = — 1 and plane y = 0 which is the thermocapillary interface (see Fig. 1). Index 1 (superscript or subscript) is used to characterize the properties of the liquid and index 2 is used to characterize the gas properties.
; i y h ' g
P — >
/ ^ ; x
Fig. 1. Geometry of fluid flow domains
The two-dimensional stationary convective fluid flows are characterized by the velocity field (ul,vl), pressure pl, temperature T1 and vapor concentration C which satisfy the Oberbeck-Boussinesq approximation of the Navier-Stokes equations. The convective flows of the liquid in the lower layer are described by the system of equations which can be presented in the following dimensionless form:
u1uX + v1uy = —pX + Au1, u1v1X + v1 v1y = —py + Av1 + GrT1,
uX + vl = 0 ()
u1 T 1x + vX = AT1.
The stationary flows in the upper gas-vapor layer (see Fig. 1) are described by the following equations:
22, 22 1—2 , — A 2
u ux + v uy = - -px + vAu ,
! p _
u2vX + v2v2 = --pi + VAv2 + GrfiT2 + GrcV2C, y p y
u2x + vy =0, (2)
u2Ti + v2Tl = PX- (AT2 + ¿AC,
V
Sc
Here p is the modified pressure, p1 = p1 + Gr/(l1T*)y, p2 = p2 + Grp/(l1T*)y. The following dimensionless parameters arise in the problem: Gr = /31T* gl3/v2 is the Grashof number, Grc = /v2, Pr = vi/xi is the Prandtl number, Sc = v2/D is the Schmidt number, V = v2/v1 is the ratio of the coefficients of the kinematic viscosity, p = p2/p1 is the ratio of the gas and liquid densities, X = X2/X1 is the ratio of the coefficients of thermal conductivity of gas and liquid, D is the coefficient of vapor diffusion, | is the thermal expansion coefficient (i = 1,2), y is the concentration expansion coefficient, 6 = S/T*, a = aT*, coefficients 6 and a characterize the Dufour and Soret effects, respectively.
Let us assume that solution of systems (1) and (2) has the form
ul = ul(y), vl = 0,
' (3)
T = (a\ + a2y)x + di(y), C = (61 + b2y)x + $(y).
Because temperature is continuous at the interface y = 0 we have a1 = a2 = A, A = const (see section 2.). Then exact solutions (3) of system of equations (1), (2) can be written as follows:
4 3 2
u1 = y4 Gra1 + GrA + ^ C1 + yc2 + C3. (4)
7 n„/„1\2 „.6 n„.A„1 „5
u2Cx + v2Cy = — (AC + aAT2).
T1 = (A + aly)x + G-T^T + — + —T^, (Gr(A)2 + 3a2c^ +
V T -¿yj 1008 (1/Pr) 144 (1/Pr) 120 (1/Pr) V 2 J
y4 1 1 y3 1 1 y2 A
+--t—, „ , (Ac-i + 2a2c2) +---—7——-( Ac2 + a2c3) +---—7—— c3 + y c4 + c5,
24 (1/Pr) 2 2> 6 (1/Pry 2 A> 2 (1/Pr) y 5'
p1 = (It Gra2 + yGrA + c^ x + y ^7 + yk6 + y ^5 + y ^4 + y ^3 + y ^2 + y'
u2 = I41 (Grflal + GrcV2b2) + y31 (GrfiA + Gr^) + y2 c1 + yc2 + c3, (7)
T2 = (A + a2y)x + -y- B (Grflal + GrcV2b2) +
y
+ — 720
1008 V
y°
6 [ b __b __] y5 [ b _
V(Orl3a2 + GrcV2b2) + 4^(GrfiA + Gr^) + ^ V (Gr|A + Gr^) + (8)
120
+3B2c1
4 3 2
y4 y3 y2
+ h [Bc + 2B202] + \B1c2 + B2c3] + Bc + yc4 + c5, 24 6 2
7 ( ) 6
y7 E __y6 1 ( _
C = (b1 + b2y)x + ^^ (Grflal + GrcV2b2) + ^ E1(Gr|a2 + GrcV2b2)+
5
+4E2 (GrfiA + GrcV2b1)l + (GrlA + Gr^) + 3E2cJ + (9)
4 J 3 12^ V 2 J
E1c1 + E2c^ + y {e1c2 + E2c^ + y E{c3 + yc6 + c7,
P2
y2 _ _ ■ — (pGr^a2 + pGr^V2 62) + y(pGr@A + pGrcV2b1) + pvc1
x+
y8— y7— y& — y5 — y4— y3— y2— — _
+ k7 + — ke + — k5 + — kA + — k3 + — k2 + — + yk0 + C8. 8 7 0 0 4 6 2
(10)
The coefficients * = — ^^, ^2 = (V/^c)a2 — (x/Pr)Sb_2 , El
(V/Sc)x2(l — aS) ' (V/Sc)(x/Pr)(l — aS)'
b2
b1
aB
E2
— a*2
Y(v/Sc)
ki, ki (i = 0,..., 7) in expressions (6), (10) for p and in expressions
l(V/Sc)
(8), (9) for T2 and C depend on the initial parameters of the problem such as pi,vi,xi,Pi, D, y, g (i = 1, 2). The values of the integration coefficients ci,ci (i = 1,..., 8), and relations for parameters b1,b2,ajj (i,j = 1,2) (see the form of exact solution (3)) are determined by the boundary conditions.
2. Boundary conditions and algorithm for computing the integration constants
f h
Let us assume that the gas flow rate in the upper layer Q is given Q = / u2(y)dy. The
J 0
characteristic value Q* of the gas flow rate is equal to p2uj. On the solid walls y = —1 and y = h of the channel the no-slip conditions should be fulfilled: u1\y=-1 = 0, v?\y=h = 0. We assume that the vapor concentration on the upper boundary y = h satisfies the condition of the absence of the vapor flux (Cy)\ h = 0 and the temperature boundary regime is determined as follows. The lower solid boundary y = —1 is assumed to be the heat-insulating wall: (Ty)\y=-1 = 0. The constant heat flux is maintained on the upper wall y = h: (T^)\y=h = 0+.
We suppose that the interface y = 0 is the thermocapillary surface that remains to be undeformed one. The continuity of the velocity and temperature are written as follows: u1\y=o = u2\y=o, T1 \y=o = T 2\y =0. The heat transfer condition includes the diffusive vapor flux
= —XM. The mass balance equation at the interface is
y=0
of the form: M = —(Cy + aT2) . Here k = k2/k1, Ki is the heat conductivity coefficient of
_ y=o
the fluid, X = XM*1/(k1T*), X is the latent heat of evaporation, M* is the characteristic value of the mass of the liquid evaporated from unit surface area per unit time, M* = Dp2/l.
The saturated vapor concentration can be found with the help of relation C\y=0 = C*[1 + e(T2\y=0 — T0)], C* is the saturated vapor concentration at T2 = T0, e = sT*, e is the parameter that includes the molar mass of the evaporated liquid, the universal gas constant, the latent heat of evaporation and T0 [1,2].
We assume that surface tension linearly depend on temperature:
a = 1 — aT (T1 — T). (11)
Here aT = MaCa/Pr, Ma = aTT*l/(p1v1x1) is the Marangoni number, Ca = u*p1v1 /a0 is the capillary number, Pr = v1/x1 is the Prandtl number. The dimensional parameters a0,aT are constants, a0 is the value of surface tension at some reference temperature and aT is the temperature coefficient of the surface tension.
In the case of normal thermocapillary effect at the interface we have
aT = a (a > 0), for T1 < T. (12)
at the interface: (T^ — kT2 — S nCy)
If the liquid exhibits the anomalous thermocapillary effect the coefficient of the surface tension can grow with temperature [10,12,13] (see also [14,17]), so that the coefficient ïïT (see (11)) has the following value:
ïïT = —ïï (a > 0), for T1 > T. (13)
The parabolic dependence of the surface tension on temperature with a minimum at a temperature interval was presented [12-14].
The kinematic condition at the interface is fulfilled automatically because the normal component of the velocity vector is equal to zero at the interface y = 0 due to (3). Projection of the dynamic condition (see, for example, [10]) on the normal and tangential vectors to the interface y = 0 gives the following relations: p1 = p2, u,y = pvu2y — ïïTA.
We note that T ^^ = Ax + c5. Then conditions (12), (13) allow us to find the coordinate
x =(T — c5)/A. (14)
The inequality Ti < T (see (12)) is fulfilled for x < x if A > 0. If A < 0, this inequality is fulfilled for x > x. In this case the choice of the value aT is determined according to (12).
2.1. Algorithm for computing the integration constants
A1. In the case when the condition of the flux absence is prescribed for vapor concentration C at the wall y = h, the equalities b2 =0 and ft(h) = 0 should be satisfied (the prime denotes the derivative with respect to y).
A2. The continuity conditions for velocity and temperature at the interface y = 0 require
c3 = C3, c5 = c5 and a1 = a-f = A.
A3. Taking into account the thermal insulation of the wall y = —1 and the boundary regime on y = h, we have a2 = 0, 0'1 ( —1) = 0 and a2, = 0, 0'2(h) = 0+, respectively.
A4. The equation of mass balance and the heat transfer condition at the interface y = 0 result in the following relations: M = —(c6 + aC4), b2 + aa2 = 0, a2 — Ka2 — SKb2 = 0, c4 — KC4 — SKC6 = —XM (two of them are fulfilled identically).
A5. The condition for the saturated vapor concentration is fulfilled if the longitudinal temperature and concentration gradients A and b1, constants C7 and C5 are related as follows:
b1 = C*eA, C7 = C* + C*e(C5 — T0).
A6. The dynamic conditions at the interface result in relations between ci and Ci (i = 1,2): C2 = C2pv — WtA, C1 = CipV.
A7. The system of linear algebraic equations for calculation of unknown constants Ci (i = 1,2, 3) follows from the no-slip conditions and the given gas flow rate. When Ci (i = 1,2, 3) are known the constants c1, c2, c3 can be found.
A8. The integration constant C6 is determined from the equation 0'(h) = 0 with the use of C
(i = 1, 2, 3).
A9. Constant c4 is determined from the equation 0'1 ( —1) = 0 when c1, c2, c3 are known.
A10. The integration constant C4 should satisfy the relation C4 = c4/(K + Xa) — C6 (SK + X)/(K + Xa). This relation is a consequence of the heat transfer and mass balance conditions at the interface y = 0.
A11. The given heat flux at the boundary y=h leads to the relation
h4 1 _ h3 h2
0+ = —Bi = (GrpA + CrcV2b1) + —Bici + —B{c2 + hB{c3 + C4. (15)
24 v 6 2
This relation should be considered as a compatibility condition that follows from the form of exact solution (3).
A12. The mass velocity of the evaporated liquid at the interface can be found with the use of the mass balance equation when c4 and c6 are known.
Assume that c5 is equal to zero (c5 = 0). Then the coordinate x on the Ox axis is x = T/A (see 14). The presented algorithm (A1)-(A12) of calculation of the integration constants cj (j = 1,..., 5) and ci (i = 1,..., 7) should be used twice in order to get two sets of the integration constants. We denote them by (c+, c+) for oT = & and (c-, c-) for aT = —o (j = 1, 2, 3,4, i = 1,2,3,4, 6). Equation (15) is the condition for the parameter 0+. It determines the temperature boundary regime on the upper solid wall of the channel.
3. Examples of the two-layer fluid flows
We consider the two-layer system, where the lower and upper layers with equal thicknesses l are filled by n-heptan and nitrogen, respectively [12,13,18]. The values of the physico-chemical parameters of the liquid and gas are presented in Tab. 1.
The Dufour and Soret parameters are S = 10-3(K) and a = —5 • 10-3 (1/K), respectively.
Table 1. Values of physico-chemical parameters of the liquid and gas
Parameter n-heptan nitrogen
p, g/cm3 0.68 1.2 • 10~3
v, cm2/s 0.568 • 10~2 0.15
k, kal/(cm s K) 0.3346 • 10~3 0.65 • 10~4
X, cm2/s - 10~3 0.3
ct0, dyne/cm 44
aT, dyne/(cm K) 0.133
D, cm2/s - 10-1
A, kal/g - 102
£, 1/K - 10~2 - 10-1
- 10-1 - 1
13, 1/K -10-3 0.337 • 10~2
We chose the thickness of the liquid layer l = 0.5 cm as the characteristic scale of the flow domain. Then u* = 0.011 cm/s. Tab. 2 presents the values of the dimensionless parameters of the gas-liquid system. If the characteristic temperature drop is equal to T* = 10 K then Gr & 40289, Grc « —3379 in the case of normal gravity (g = 981 cm/s2) and Gr « 402.89, Grc & —33.79 in the case of low gravity (g = 9.81 cm/s2). According to [13] (see also [14]) we chose the dimensional value T equal to 35oC. It defines the behavior of the surface tension with respect to the temperature at the interface.
We model the two-layer fluid flows for the system "n-heptane-nitrogen" in the case of positive values of the longitudinal temperature gradient A: A = 0.1 K/cm, A = 0.01 K/cm. The temperature coefficient of the surface tension aT changes sign at the point x. The dimensional velocity profiles are shown in Fig. 2-6. The figure captions contain dimensional values of longitudinal temperature gradient, gas flow rate denoted with quantities and gravity acceleration g.
Pictures 2 and 3 demonstrate the velocity profiles at various values of the longitudinal temperature gradient A at the interface. The qualitative and quantitative differences of the flow structure and sharp differences in the fluid flows at the right and the left of the point x are observed with the increase of the longitudinal gradient A (dimensional analog X of the point x is equal to 3500 (cm) for flows demonstrated in Fig. 2 and X = 350 (cm) for flows in Fig. 3). The velocity value at the interface is equal to U = 0.25 cm/s if A = 0.1 K/cm (see Fig. 3, solid line) and U = 1.16 cm/s if A = 0.01 K/cm (see Fig. 2, solid line). Here the temperature coefficient of surface tension is equal to aT = 0.133 dyne/(cm K) (in the dimensionless case we have aT = 0.03). In the case of anomalous temperature dependence of the surface tension when aT = —0.133 dyne/(cm K) (aT = —0.03) the value of velocity at the interface is U = 2.85 cm/s for A = 0.1 K/cm (Fig. 3, dashed line) and U = 1.42 cm/s for A = 0.01 K/cm (Fig. 2, dashed line). Figs. 2 and 3 present the flows under conditions of normal gravity (g = 981 cm/s2) in the case that the thickness of the gas and liquid layer is equal to 0.5 cm, dimensional gas flow rate is Q = 1.8 • 10-3 g/(cm s). The dimensionless parameters are h =1, l = 1, Q = 246.1, Gr « 40289. The values of velocity at the interface U and of the mass rate of evaporated liquid M are presented in Tab. 3 for various values of the longitudinal temperature gradients. More intensive evaporation is observed for larger value of the longitudinal temperature gradient and also in the case of anomalous effect of thermocapillarity. Fig. 4 demonstrates the flow structure in the system away from the point of change of thermocapillarity character x in the case that the longitudinal temperature gradient is equal to A = 0.1 K/cm.
A decrease of intensity of the gas flow rate results in quantitative and some qualitative differences in the flow (see Figs. 2 and 5). The velocity at the interface has the value U = 1.16 (Fig. 2, solid line) and U = 0.025 cm/s (Fig. 5, solid line) for positive value of temperature coefficient of the surface tension aT = 0.133 (aT = 0.03) if the gas flow rate is equal to Q = 1.8 • 10-3 g/(cm s) and Q = 1.8 • 10-4 g/(cm s), respectively. For negative value aT = —0.133 (aT = —0.03) the values of velocity at the interface are equal to U = 1.42 (Fig. 2, dashed line) and U = 0.28 cm/s (Fig. 5, dashed line) in the case of gas flow rate Q = 1.8 • 10-3 g/(cm s) and Q = 1.8 • 10-4 g/(cm s), respectively. Here X = 3500 (cm) for flows demonstrated in Fig. 2 and Fig. 5. Tab. 4 presents the values of the interface velocity U and the mass rate of evaporated liquid M for two values of the gas flow rate Q = 1.8 • 10-3 g/(cm s) and Q = 1.8 • 10-4 g/(cm s). It should be noted that for a larger value of the specific mass gas flow rate Q the flow and
the same symbols A and Q as dimensionless
Table 2. Values of dimensionless parameters of the gas-liquid system
Parameter Value
P 1.76 • 10~3
V 26.41
K 0.19
X 397.35
ß 3.18
Pr 7.52
Sc 1.63
a 5•10-2
S 10-4
£ 0.49
X 2.78
<7 0.03
evaporation processes are more intensive.
Fig. 2. Velocity profiles for positive and negative values of aT. Here A = 0.01 K/cm, Q = 1.8 • 10~3 g/(cm s), g = 981 cm/s2
0 1 2 3 4 5
u(y), cm/s
Fig. 3. Velocity profiles for positive and negative values of aT. Here A = 0.1 K/cm, Q = 1.8 • 10~3 g/(cm s), g = 981 cm/s2
The velocity profiles are presented in Fig. 6 in the case of low gravity g = 9.81 cm/s2 (Gr « 402.89). Quantitative flow characteristics are changed slightly in comparison with the case g = 981 cm/s2 (see Fig. 5). One can observe appearance of a reverse flow near the interface. The values U and M under normal (g = 981 cm/s2) and low gravity (g = 9.81 cm/s2) are presented in Tab. 5 in the case that Q = 1.8 • 10~4 g/(cm s), A = 0.01 K/cm. The value of the
velocity at the interface U is equal to —0.003 cm/s (Fig. 6, solid line). The value of parameter X is equal to 3500 (cm) for flows demonstrated in Fig. 5 and Fig. 6. Values of the mass of evaporated liquid are of the same order in the cases of normal and low gravity.
340 345 350 355 350 365
x, cm
Fig. 4. Velocity profiles for positive and negative values of aT (right and left). Here A = 0.1 K/cm, Q = 1.8 • 10~3 g/(cm s), g = 981 cm/s2
D 0.1 0.2 0.3 0.4 0.5
u(y), cm/s
Fig. 5. Velocity profiles for positive and negative values of aT. Here A = 0.01 K/cm, Q = 1.8 • 10~4 g/(cm s), g = 981 cm/s2
-0.1 0 0.1 0.2 0.3 0.4 0.5
u(y), cm/s
Fig. 6. Velocity profiles for positive and negative values of aT. Here A = 0.01 K/cm, Q = 1.8 • 10~4 g/(cm s), g = 9.81 cm/s2
Table 3. Values of interface velocity U and mass rate of evaporated liquid M; Q = 1.8 • 10 3
g/(cm s), g = 981cm/s2
aT = 0.133 dyne/(cm K) aT = -0.133 dyne/(cm K)
A, K/cm U, cm/s M, g/(cm2 s) U, cm/s M, g/(cm2 s)
0.1 0.25 3.29 • 10~4 2.85 1.45 • 10~3
0.01 1.16 6.45 • 10~5 1.42 7.56 • 10-5
Table 4. Values of interface velocity U and mass rate of evaporated liquid M; A = 0.01 K/cm, g = 981cm/s2
aT = 0.133 dyne/(cm K) aT = -0.133 dyne/(cm K)
Q, g/(cm s) U, cm/s M, g/(cm2 s) U, cm/s M, g/(cm2 s)
1.8 • 10~3 1.16 6.45 • 10~5 1.42 7.56 • 10~5
1.8 • 10~4 0.025 3.29 • 10~6 0.28 1.45 • 10~5
Table 5. Values of interface velocity U and mass rate of evaporated liquid M; Q = 1.8 • 10 4
g/(cm s), A = 0.01 K/cm
aT = 0.133 dyne/(cm K) aT = -0.133 dyne/(cm K)
g, cm/s2 U, cm/s M, g/(cm2 s) U, cm/s M, g/(cm2 s)
981 0.025 3.29 • 10~6 0.28 1.45 • 10~5
9.81 -0.003 1.24 • 10~6 0.26 1.24 • 10~5
Conclusion
The constructed solutions to the Navier-Stokes equations in the Oberbeck-Boussinesq approximation allow us to model the two-layer fluid flows with the thermocapillary interface including the evaporation effects at the interface and the Dufour and Soret effects in the gas-vapor phase. The evaporation effect is modeled by the use of the thermal boundary condition and mass balance equation at the interface. The anomalous thermocapillarity effect at the interface is also taken into consideration. For the n-heptanol-nitrogen systems the velocity profiles that characterize the fluid flows are presented. The effects of gravity and the longitudinal temperature gradient along the interface on the fluid flows are investigated. The obtained results demonstrate the quantitative and qualitative differences in the flow of liquid in the case of normal and anomalous thermocapillary effect and also in the conditions of normal and low gravity.
Authors gratefully acknowledge the support of this work by the Russian Foundation of Basic Research, grant 14-08-00163.
References
[1] O.N.Goncharova, M.Hennenberg, E.V.Rezanova, O.A.Kabov, Modeling of the convective fuid fows with evaporation in the two-layer systems, Interfacial Phenomena and Heat Transfer (IHMT), 1(2013), 389-396.
[2] O.N.Goncharova, E.V.Rezanova, Example of an exact solution of the stationary problem of two-layer flows with evaporation at the interface, Journal of Applied Mechanics and Technical Physics, 55(2014), no. 2, 247-257.
[3] O.N.Goncharova, Yu.V.Lyulin, E.V.Rezanova, O.A.Kabov, Modeling of the stationary problem of two-layer flows of the liquid and gas by evaporation, Journal of Thermophysics and Aeromechanics, 22(2015), no. 5, 655-661.
[4] M.I.Shliomis, V.I.Yakushin, Convection in a two-layers binary system with an evaporation, Collected papers: Uchenye zapiski Permskogo Gosuniversiteta, seriya Gidrodinamika, 4(1972), 129-140 (in Russian).
[5] R.V.Birikh, Thermocapillary convection in a horizontal layer of liquid, J. Appl. Mech. Tech. Phys, 3(1969), 43-45.
[6] Yu.V.Sanochkin, Some problems of thermocapillary fluid flows, Prikladnaya Mechanika i Technicheskaya Fizika, 5(1989), 83-88 (in Russian).
[7] V.V.Popov, Mixed convection in a two-layer liquid, Theoretical principles of chemical technology, XV(3)(1981), 398-404 (in Russian).
[8] V.K.Andreev, V.B.Bekezhanova, Stability of the non-isothermal fluids, Krasnoyarsk, Siberian Federal University, 2010 (in Russian).
[9] G.A.Ostroumov, Free convection in the conditions of the internal problem, Moscow, Leningrad: Gos. izd-vo tehniko-teoreticheskoy literatury, 1952 (in Russian).
[10] V.K.Andreev, Yu.A.Gaponenko, O.N.Goncharova, V.V.Pukhnachov, Mathematical models of convection (de Gruyter Studies in Mathematical Physics), Berlin/Boston: De Gruyter, 2012.
[11] B.Gebhart, Y.Jaluria, R.L.Mahajan, B.Sammakia, Bouyancy-induced flows and transport, Berlin, Heidelberg, New York, London, Paris, Tokyo, Springer Verlag, 1988.
[12] M.Cl.Limbourg-Fountaine, G.Petre, J.-C.Legros, Texus 8 experiment: Effects of a surface tension minimum on thermocapillary convection, Physicochemical Hydrodynamics, 6(1985), no. 3, 301-310.
[13] J.-C.Legros, M.Cl.Limbourg-Fountaine, G.Petre, Surface tension induced convection in presence of a surface tension minimum, Physicochemical Hydrodynamics, Interfacial Phenomena. In M.G. Velarde Ed., NATO ASI series, Series B: Physics, New York, Plenum Press, 174(1988), 209-227.
[14] G.Petre, M.Cl.Limbourg-Fountaine, J.-C.Legros, Study of the surface tension minimum of aqueous alcohol solutions and movements at the interface air/solution, Proc. of the 4th Europ. Symp. on Material and Sciences under microgravity, Madrid, Spain, 5-8 April 1983, 199-200.
[15] A.Cloot, G.Lebon, Marangoni convection induced by a nonlinear temperature-dependent surface tension, J. Physique, 47(1986), 23-29.
[16] R.Savino, A.Cecere, R.Di Paola, Surface tension-driven flow in wickless heat pipes with self-rewetting fluids, Int. J. Heat and Fluid Flow, 30(2009), no. 2, 380-388.
[17] O.N.Goncharova, Yu.O.Kabova, O.A.Kabov, Thermocapillary convection in a free liquid layer in the presence of an adjacent gas flow, Computational Thermal Sciences, 5(2011), no. 3, 389-396.
[18] A.A.Ravdel, A.M.Ponomareva, Hand-book of physico-chemical quantities, Sankt-Peterburg, Specialnaya literatura, 1998 (in Russian).
Моделирование двухслойных течений с испарением на границе раздела при наличии аномального термокапиллярного эффекта
Ольга Н. Гончарова Екатерина В. Резанова
Стационарные конвективные течения двух несмешивающихся жидкостей (жидкости и газа), находящиеся под действием поперечного поля силы тяжести и продольного градиента температуры, созданного вдоль границы раздела, изучаются аналитически. Математическое моделирование течений жидкостей c учетом эффектов испарения на границе 'раздела проводится на основе точных решений системы уравнений Навье-Стокса в приближении Обербека-Буссинеска. Учитываются также эффекты термодиффузии и диффузионной теплопроводности. Построенные решения применяются для описания течений в двухслойной системе жидкости и газа в случае, когда жидкость характеризуется аномальным термокапиллярным эффектом на границе раздела. Приводятся примеры двухслойных течений при различных значениях удельного расхода газа, продольного градиента температуры и ускорения силы тяжести.
Ключевые слова: математическая модель, граница раздела, испарение, точное решение, аномальный термокапиллярный эффект.
