Two-dimensional convection of an incompressible viscous fluid with the heat exchange on the free border Текст научной статьи по специальности «Физика»

Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki

[J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 3, pp. 567-577

ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1483

MSC: 76F02, 76F45, 76M45, 76R05, 76U05

TWO-DIMENSIONAL CONVECTION OF AN INCOMPRESSIBLE VISCOUS FLUID WITH THE HEAT EXCHANGE ON THE FREE BORDER

S. S. Vlasova1, E. Yu. Prosviryakov2

1 Kazan National Research Technical University named after A. N. Tupolev, 10, Karl Marx st., Kazan, 420111, Russian Federation.

2 Institute of Engineering Science, Ural Branch of RAS, 34, Komsomolskaya st., Ekaterinburg, 620049, Russian Federation.

Abstract

The exact stationary solution of the boundary-value problem that describes the convective motion of an incompressible viscous fluid in the two-dimensional layer with the square heating of a free surface in Stokes's approach is found. The linearization of the Oberbeck-Boussinesq equations allows one to describe the flow of fluid in extreme points of pressure and temperature. The condition under which the counter-current flows (two counter flows) in the fluid can be observed, is introduced. If the stagnant point in the fluid exists, six non-closed whirlwinds can be observed.

Keywords: exact solution, Newton-Rikhmann law, thermal convection, Oberbeck-Boussinesq equations, counter-current flow.

Introduction. The research of the convective flows of an incompressible viscous fluid is caused by a considerable drop of temperature in a wide range of processes connected with the deformation of dissipative environments. The convection induced by the non-uniform heating of incompressible and compressed substances is the most widespread kind of gas dynamics and fluid flows in the Universe. The convective motion of the fluid heated from below in a two-dimensional horizontal layer is one of the most popular subjects of studying. The first example of self-organization of the nonlinear phenomena is Rayleigh-Benard convection. The choice of the two-dimensional layer as the abstract mathematical object is mainly connected with the fact that this geometry can be quite easily realized in an experiment and provides certain conveniences for taking thermal and optical measurements. The two-dimensional horizontal layer is a matter of great

© 2016 Samara State Technical University. Please cite this article in press as:

Vlasova S. S., Prosviryakov E. Yu. Two-dimensional convection of an incompressible viscous fluid with the heat exchange on the free border, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 3, pp. 567-577. doi: 10.14498/vsgtu1483. Authors Details:

Svetlana S. Vlasova ([email protected]; Corresponding Author), Postgraduate Student, Dept. of Machinery Engineering & Engineering Graphics.

Evgeny Yu. Prosviryakov (Cand. Phys. & Math. Sci.; [email protected]), Leading Researcher, Laboratory of Applied Mechanics.

importance in connection with the theory of convective stability applications in meteorology, geophysics and astrophysics [1,2].

The first exact solution of the natural convection for a two-dimensional layer of fluid with a constant gradient of the temperature on the borders of the layer is described in the article [3], where two cases of boundary conditions for the velocity are considered. A brief survey of the articles and reviews, in which the possibilities of finding exact solutions to free convection equations and the research of stability of these solutions were studied [4-10], is given in the papers [1,2] and the reference lists.

In the work [5], a method for the description of heat exchange with the localized parabolic heating of the border is offered. This method is offered in the context of the exact solutions class in which velocities linearly depend on horizontal coordinates, and fields of the pressure and temperature are distributed under the square law [2, 6, 11]. In this case, finding the exact solution is reduced to solving the nonlinear system of one-dimensional evolution heat conduction equations like heat conductivity, and stationary equations of the gradient type. When finding the exact solutions, which describe the convective fluid flow [1-10], an ideal heat transfer on the borders is assumed. The aim of the present work is to investigate the motion of a viscous incompressible fluid with heat exchange on the upper boundary.

1. Mathematical model and main equations. The plane layer stationary convection of a viscous incompressible fluid (Fig. 1) can be presented by the system of Oberbeck-Boussinesq equations:

TdVx 0VX dP ( d 2VX d2Vr

dx dz dx V dx2 dz2

3Vr dVZ dP Vx-V- + Vz ^ = - — + vAVz + gfiT, dx dz dz

dVx 3VZ (1)

—x +--z = 0,

dx dz

9T vdTL- ( d^L d^T )

x dx z dz v dx2 dz2 J '

In the system of equations (1), the following designations are introduced: Vx, Vz are the velocities that are parallel to x axis and z axis, respectively; P is the pressure deviation from hydrostatic pressure divided by the constant average fluid density p; T is a deviation from the average temperature; v, x, P are the dissipative coefficients of the kinematic viscosity, heat diffusivity and thermal

z

h 9

I X

Figure 1. Schematic view of the model problem

expansion of fluid, respectively; g is the acceleration of gravity;

dx2 dz2

is the two-dimensional Laplace operator written down in the Cartesian orthogonal coordinate system.

We search the stationary solution of the system (1) in the form of [2,6,11]:

Vx = xu(z), Vz = w(z),

x2 x2 (2)

T = To(z)+ Tn(z) y, P = Po(z) + Pn(z) y. ()

Note that if the temperature T is formally substituted with the concentration function C in the system (1) the solutions (2) are also true for the concentration convection. ft is the coefficient of concentration fluid extension, in this case.

We substitute the class of solutions (2) in the system of equations (1) and receive the following nonlinear system defining unknown functions u, w, T0, Tn, P0:

d2u 2 du ^ dw

= u + w— + P11, — = -u, dz2 dz dz

dP0 d2w nrT1 dw

= ^ttt + gftTo - w—, dz dz2 dz (3)

dPii RT dTo ( + d2To \ ()

— = gftTi1, w— = X(Ti1 + -d^), 0 rp ^ dT11 d2T11 2uTn + w—— = X—Tor.

dz dz2

We reduce the system of ordinary ninth-order differential equations (3) to the dimensionless form. We introduce the following characteristic geometric scale values as basis [2,6]: h is the transverse characteristic size, l is the axial characteristic size. The units of measurements are: 0 is for the temperature, VZ and Vx are for the velocities 9l3®21 and respectively. Thus, the accounting of geometrical

anisotropy of the task results in the formation of a flow and availability of two scales for the velocities, since

gft 0h4 _ gft0h3 = ;

the pressure for an incompressible fluid divided by the constant density is gft0h; 5 = h/l.

We write down the system (3) in the dimensionless form:

dP11 T dw 2 (2 T + dTn \ 1 d2Tn

—— = T11, — = -u, Gr 5 2uTn + w—- = -—-r^-dz dz \ dz J Pr dz2

r2 dT0 1 (r2 d2T0 \ c2 dw dP0 c2 d2 w ^ Gr5 w—— = — (5 ^ —¡-4- , Gr5 w— = -—0 + 52— + T0, (4) dz Pr V dz2 J dz dz dz2

Gr 52('° + 5^ \ = -P11 + 52 g.

Both dimensionless and dimensional (3) variables and functions in the system (4) are designated by the same symbols. Only dimensionless variables are used further on. In this work

_ gPeh5

Gr = w

is the modified Grashof number.

We find the solution of the system (4) in the extremum points of the temperature. For this purpose we linearize the equations (3). The linearization of the nonlinear system (4) is possible when Gr 52 ^ 1. In this case, the received exact solutions are fair with any Grashof number they satisfy the assessed value Gr ^ 1/52 where Gr e (0; 1/52). It should be noted that with an isotropic geometry = 1), we receive the classical criterion of reducing the equations (1) to Stokes equation [11].

The linearized system (3) in the dimensionless form is written as:

d2Tn dPn d2u dw

0; —r~ = Til, -r^ = Pll, -r- = -u,

dz2 ' dz ' dz2 ' dz (5)

d2To A2 t dPo T + A2 d2w ()

= —0 T11, — = To + 0

dz2 ' dz dz2 '

The equations in the system (5) are written out in the order, the integration of equations is effected.

2. Boundary conditions and exact solutions. The system of the linear differential equations (5) has the exact polynomial solution:

z2

T11 = Ciz + C2, P11 = Ci -+ C2z + C3,

z4 z3 z2 '

u = Cl — + C2 3! + C3+ C4z + C5,

y4 r3 y2

w = —Cl - — C2^ — C33! — C4^ — C5z + C6, (6)

z3 z2

To = —C1O2 3! — C2O2 2! + Czz + Cs,

.y4 r3 J2

Po = —2C1O2 - — 2C2O2 3! + (C7 — 02C3) -2 + (Cs — 02 C4) z + C9.

We formulate boundary conditions for finding the constants of integration—the coefficients of polynomials (6). The heat source is introduced on the lower solid bound:

T = A + Bx2 = -|(1 — x2) with dimensionless variables, it is written as:

T = 1(1 — x2).

The boundary conditions for defining the constants of integration of dimensionless solutions (5) can be written as:

z = —1: w = u = 0, T11 = —1, To = 1; (7)

z = 0: w = 0, P11 =0, Po = 1, dTo = - BiTo, ddTl1 = - BiTn. (8)

dz dz

The no-slip condition is set on the lower bound (z = -1), and the heat exchange under Newton-Rikhman law is set on the free plane layer bound [12]. Bi is Biot number [12]. It is obvious that the integration constants defining the structure of the exact solution (6) by boundary conditions (7) and (8) can be as follows:

C3 = Co = 0, C9 = 1. C Bi C 1 C 2(4 Bi +15) C 3 Bi +10

Cl = ———T, C2 = ————, C4 = —, , , C5 =

Bi + 1' 2 Bi + 1' 4 5!(Bi + 1) ' 5 5!(Bi + 1)' = ¿2 Bi(Bi + 3) - 3 Bi(Bi + 1) = 3(Bi + 1) - ¿2(Bi + 3)

1 = 3!(Bi + 1)2 , Cs = 3!(Bi + 1)2 .

Let's substitute (6) with the boundary conditions (7) and (8) into the equations (2), then the expressions of hydrodynamic fields are written in the following form:

/ Bi z4 z3 2(4 Bi +15)z 3 Bi +10 \

= H4!(Bi +1) - 3!(Bi +1) + 5!(Bi +1) + 5!(Bi +1)J,

/ Bi z5 z4 (4 Bi +15) z2 (3 Bi +10) z

Vz = - ( .—TTT +

V5!(Bi +1) 4!(Bi +1) 5!(Bi +1) 5!(Bi +1)

Bi ¿2z3 d2z2 Ô2 Bi (Bi +3) - 3 Bi (Bi +1)

= - 3!(Bi +1) + 2!(Bi +1) + 3!(Bi +1)2 z+

3(Bi +1) - ¿2(Bi +3) x2(Bi z - 1) + 3!(Bi +1)2 + 2!(Bi +1) :

2 Bi ô2z4 2ô2z3 2 (¿2 Bi (Bi +3) - 3 Bi (Bi +1)) z2 = -4!(Bi +1) + 3!(Bi +1) + 4! (Bi +1)2 +

3(Bi +1) - ¿2 (Bi +3) 2(4 Bi +15) ¿2 ^

5!(Bi

+ 1n2KB^î)

H 3! (Bi +1)2 5!(Bi +1) Jz+

( Bi z2 z \ x2

Bi +1 2

3. Analysis of the plane convective motion of the fluid. Since the solutions (6) of the system of equations (5)are polynomial, the analysis of flows can always be converted to the solution of the generalized Raus-Gurvits problem. We consider the characteristic properties of the velocities depending on the value of Biot number. Considering the boundary conditions (6) and (7), we find out that the variety of the function's values can be presented in the following form:

/ ^f Bi z3 (Bi + 4)z2

u = (z +i)f (z) = (z + - IKB^T) +

(Bi +4)2 3 Bi +10

4!(Bi +1) 5!(Bi +1)/'

The function f is a cubic polynomial with the coefficients depending on number Bi under the linear-fractional law. We investigate the spectral properties of a polynom f on the domain of the definition. It is known that the function contains an odd number of zero values (the quantity of stagnant points of the flow) when an inequality is being solved:

f (-1)f (0) < 0

and an even number when the opposite inequality is being solved. Thus, the function f has the only one solution in the interval with

Bi e (-to; -10/3) U (-5/2; -1) U (-1; +to) .

When Bi = 0, the assumed cubic polynomial degenerates into a linear function. Two solutions are possible in the one case with Bi e (-10/3; -5/2). Let's consider graphs of the stream function and vorticity when Bi = -2.7 (Fig. 2) and Bi = -2 (Fig. 3). The expressions are assumed for the stream function of the velocity:

Bi 22

* = XZ(Z + 1)^ Bi 22 I (2 Bi +5)2 3 Bi +10 A

* = -x2(2 + 1) 5!(Bi +1) + 5!(Bi +1) - 5!(Bi +1)J

and vorticity:

_ / Bi 23 22 2(4 Bi +15) \

ly = H3!(Bi +1) - 2!(Bi +1) + 5!(Bi +1) ).

Investigating the localization of the polynom f roots with the values of numbers Bi = -10/3 and Bi = -5/2, we find out that the function takes the zero

-0.5 0.0

X

1.0 -1.0

-0.5 0.0

X

Figure 2. Isolines of the current function (left) and the function of vorticity (right)

with Bi = -2.7

z

z

- 0.2

- 0.4

- 0.6

- 0.8

-1.0

- 0.2

- 0.4

-1.0 - 0.5 0.0 0.5 1.0

x

- 0.6

- 0.8

-1.0

Figure 3. Isolines of the current function (left) and the function of vorticity (right) with Bi = —2

values in the intervals of definition z = 0, z = -(1 + \/2l)/10 and z = —1, z = (1 — \/6)/5 respectively.

Let's note that for the velocity there is a point different from z = —1 in which the velocity does not depend on the values of Biot number. To prove this statement, we take two Biot numbers not equal among one another and substitute them in the expression for a gradient u. Subtracting one polynomial function from another, we receive the equation for the definition of this point within the interval [—1; 0]:

5z4 + 20z3 — 22z — 7 = 0.

Calculating the roots of this equation, we receive the point in which the velocity does not depend on Biot number, and the value of velocity is:

z = —0.370734, u(—0.370734) = 0.01929.

Now we consider the characteristic properties of the velocity parallel to an axis of z-coordinates. Its value as well as the value of velocity Vx does not depend on parameter S. If we consider the boundary conditions (7) and (8), we find out that the multiplicity of the function Vz values can be presented in the following form:

w = z(z + 1)2/(z) = z(z + 1)2( - JBi+T +

(2 Bi +5)z 3 Bi +10 \ 5!(Bi +1) + 5!(Bi +1) — 5!(Bi +1)J.

The function f is a square polynomial with the coefficients depending on Bi under the fractional-linear function. The function f has the only one root within the interval [—1;0] with Bi e (—10/3; —5/2). If Bi = 0, the initial square polynomial degenerates into a linear function. When an inequality f (—1)f(0) > 0 is being solved in the interval [—1;0] no roots are found. Analyzing the arrangement of square polynomial roots f with the values of numbers Bi = —10/3 and Bi = —5/2, we find out that the function takes these values in the intervals of definition z = 0 and z = —1, respectively.

z

Let's consider the additive components of the temperature. The square additive component T11 is distributed under the linear law. The variety of the function T11 values for 2 e [—1; 0] takes negative values with Biot Bi e [-1;0). The received T11 solution takes up the constant value when Bi = 0: T11 = -1. The additive component is monotonously increasing with Bi e (-to; -1) U (0; +to), otherwise it is monotonously decreasing. Function T11 vanishes in the interval of the definition in the point 2 = 1/Bi with Bi e [-1; 0).

Now let's consider the properties of the background temperature and pressure. The background temperature can be presented in the following form:

To = Bi C2523- - C2522- - Bi Cg2 + Cs.

The existence of an odd number of solutions is equivalent to an inequality being solved:

To(-1)To(0) = Cs (Cs(Bi + 1) - C3|2(Bi + 3)) < 0.

Due to the boundary conditions (7) and (8), the temperature can be presented in another form:

To = 2 + 3!(Bi++1|)2 ( 22Bi52(Bi +1) + 2d2(Bi +1)(Bi +3)- 3 Bi(Bi +1) - 52(Bi +3)) =

= 1+ (2 +1) g (2)

= 2 + 3!(Bi +1)2g (2).

Let's investigate the function g(z). The function g(z) is a square polynomial function depending on two parameters: 5 and Bi. The function g(x) takes a zero value in the interval of the definition when an inequality is being solved:

g(-1)g(0) = (3Bi + 3Bi2 + 352 + Bi52)(3Bi + 3Bi2 + 652 + 6Bi52 + 2Bi252) < 0.

The function g(2) has two solutions when the system of inequalities is being solved:

D> 0,

g(0)g(-1) > 0.

D is the discriminant of a quadratic equation. If the minimum of the function To is not more than 0.5, the number of zeros of the function To is equal to the number of zeros.

Further on, additive components of the pressure are analyzed. We provide isolines of the temperature and pressure when Bi = -2 and 5 = 1 (Fig. 4) and when Bi = -0.4 and 5 = 0.48 (Fig. 5).

It is obvious that the square law change of P11 converts into a linear dependence P11 = -2 when Bi = 0. The function P11 reduces to zero within the interval [-1; 0) in the point 2 = 2/ Bi. Since the square curve turns into a linear dependence, the function P11 monotonously decreases.

Let's consider the background pressure function P0(2). The function P0(2) is a quartic polynomial depending on two parameters: 5 and Bi. The function

Figure 4. Isolines of the functions of temperature and pressure with Bi = —2 and S = 1

Figure 5. Isolines of the function of temperature pressure with Bi = —0.4 and S = 0.48

P0(2) takes one zero value in the interval of the definition 2 e [-1;0] when the inequality is being solved:

P(-1)P(0) < 0.

The function of the pressure comes up to the maximum with some values of parameters 5 and Bi in the considered interval. The existence of the extremum points in the interval of the definition 2 e [-1;0] is confirmed by the inequality solution:

P'(-1)P/(0) < 0.

Thus, the structure of the function for the temperature shows that some locally hyperbolic level lines exist. In other words, they can not be closed with any Biot number and 5 values. The local ellipticity isolines of the pressure can be observed in a rather wide range of the dimensionless complexes defining the topology of the fluid flow.

Conclusion. In this paper the convective motion of the two-dimensional flow of a viscous incompressible fluid under the Newton-Rikhman law on one of the borders of an infinite fluid layer in Stokes's approach has been analyzed.

The assessed value of the function that makes the linearization of the Navier-Stokes equation in the Oberbeck-Boussinesq approach possible has been obtained. It is shown that the considered system of ordinary differential equations received within the announced class of exact solutions, exactly describes a fluid flow in extreme points of the temperature and pressure. The values when counter-current flows in the fluid can be observed, are found. It is shown that depending on similarity numbers in fluid, a different number of whirlwinds and the local ellipticity of pressure isolines can be observed.

z

z

x

x

z

x

x

Declaration of Financial and Other Relationships. This work was supported by the Foundation for Assistance to Small Innovative Enterprises in Science and Technology (the UMNIK program); the agreement no. 8389 GU2/2015. Each author has participated in the article concept development and in the manuscript writing. The authors are absolutely responsible for submitting the final manuscript in print. Each author has approved the final version of manuscript. The authors have not received any fee for the article.

ORCIDs

Svetlana S. Vlasova: http://orcid.org/0000-0002-1458-8760 Evgeny Yu. Prosviryakov: http://orcid.org/0000-0002-2349-7801

REFERENCES

1. Getling A. V. Formation of spatial structures in Rayleigh—Benard convection, Sov. Phys. Usp., 1991, vol.34, no. 9, pp. 737-776. doi: 10.1070/pu1991v034n09abeh002470.

2. Aristov S. N., Prosviryakov E. Yu. On one class of analytic solutions of the stationary axisymmetric convection Benard-Maragoni viscous incompreeible fluid, Vestn. Samar. Gos. Tekh. Univ., Ser. Fiz. Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2013, no. 3(32), pp. 110-118 (In Russian). doi: 10.14498/vsgtu1205.

3. Birikh R. V. Thermocapillary convection in a horizontal layer of liquid, J. Appl. Mech. Tech. Phys., 1966, vol.7, no. 3, pp. 43-44. doi: 10.1007/bf00914697.

4. Andreev V. K., Bekezhanova V. B. Stability of non-isothermal fluids (Review), J. Appl. Mech. Tech. Phys., 2013, vol.54, no. 2, pp. 171-184. doi: 10.1134/s0021894413020016.

5. Aristov S. N., Shvarts K. G. Vikhrevye techeniya advektivnoi prirody vo vrashchayushchem-sya sloe zhidkosti [Advective Eddy Flows in a Rotating Liquid Layer]. Perm, Perm. Gos. Univ., 2006 (In Russian).

6. Aristov S. N., Shvarts K. G. Vikhrevye Techeniya v tonkikh sloyakh zhidkosti [Eddy Flows in Thin Liquid Layers]. Kirov, Vyat. Gos. Univ, 2011 (In Russian).

7. Aristov S. N., Shvarts K. G. Convective heat transfer in a locally heated plane incompressible fluid layer, Fluid Dynamics, 2013, vol.48, no. 3, pp. 330-335. doi: 10.1134/ s001546281303006x.

8. Goncharova O. N., Rezanova E. V. Modeling of two-layer fluid flows with evaporation at the interface in the presence of the anomalous thermocapillary effect, J. Sib. Fed. Univ. Math. Phys., 2016, vol.9, no. 1, pp. 48-59. doi: 10.17516/1997-1397-2016-9-1-48-59.

9. Efimova M. V. On one two-dimensional stationary flow of a binary mixture and viscous fluid in a plane layer, J. Sib. Fed. Univ. Math. Phys., 2016, vol. 9, no. 1, pp. 30-36. doi: 10. 17516/1997-1397-2016-9-1-30-36.

10. Goncharova O. N., Kabov O. A., Pukhnachov V. V. Solutions of special type describing the three dimensional thermocapillary flows with an interface, Int. J. Heat Mass Transfer., 2012, vol.55, no.4, pp. 715-725. doi: 10.1016/j.ijheatmasstransfer.2011.10.038.

11. Aristov S. N., Prosviryakov E. Yu. A New Class of Exact Solutions for Three Dimensional Thermal Diffusion Equations, Theor. Found. Chem. Eng., 2016, vol.50, no. 3, pp. 286-293. doi:10.1134/s0040579516030027.

12. Gershuni G. Z., Zhukhovitskii E. M. Konvektivnaya Ustoichivost' Neszhimaemoi Zhidkosti [Convective Stability of An Incompressible Fluid]. Moscow, Nauka, 1972 (In Russian).

Received 13/III/2016;

received in revised form 25/V/2016;

accepted 09/IX/2016.

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2016, 2016. T. 20, № 3. С. 567-577 ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1483

УДК 532.51

ПЛОСКАЯ КОНВЕКЦИЯ ВЯЗКОЙ НЕСЖИМАЕМОЙ ЖИДКОСТИ ПРИ ЗАДАННОЙ ТЕПЛООТДАЧЕ НА СВОБОДНОЙ ГРАНИЦЕ

С. С. Власова1, Е. Ю. Просвиряков2

1 Казанский национальный исследовательский технический университет им. А. Н. Туполева, Россия, 420111, Казань, ул. Карла Маркса, 10.

2 Институт машиноведения УрО РАН,

Россия, 620049, Екатеринбург, ул. Комсомольская, 34.

Аннотация

Найдено точное стационарное решение краевой задачи, описывающее конвективное движение вязкой несжимаемой жидкости в плоском слое при квадратичном нагреве свободной поверхности в приближение Сток-са. Линеаризация уравнений Обербека—Буссинеска позволяет описать движение жидкости в точках экстремумов давления и температуры. Выведено условие, при котором наблюдается противотечение (два встречных потока) в жидкости. При наличии застойной точки в жидкости наблюдается шесть незамкнутых вихрей.

Ключевые слова: точное решение, закон Ньютона-Рихмана, тепловая конвекция, уравнения Обербека-Буссинеска, противотечение.

Декларация о финансовых и других взаимоотношениях. Работа выполнена при поддержке фонда содействия развитию малых форм предприятий в научно-технической сфере (программа УМНИК), договор № 8389 ГУ2/2015. Все авторы принимали участие в разработке концепции статьи и в написании рукописи. Авторы несут полную ответственность за предоставление окончательной рукописи в печать. Окончательная версия рукописи была одобрена всеми авторами. Авторы не получали гонорар за статью.

ОКСГО8

Светлана Сергеевна Власова: http://orcid.org/0000-0002-1458-8760 Евгений Юрьевич Просвиряков: http://orcid.org/0000-0002-2349-7801

Поступила в редакцию 13/111/2016; в окончательном варианте — 25/У/2016; принята в печать — 09/1Х/2016.

© 2016 Самарский государственный технический университет. Образец для цитирования

Vlasova S. S., Prosviryakov E. Yu. Two-dimensional convection of an incompressible viscous fluid with the heat exchange on the free border // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. T. 20, № 3. С. 567-577. doi: 10.14498/vsgtu1483. Сведения об авторах

Светлана Сергеевна Власова ([email protected]; автор, ведущий переписку), аспирант, каф. машиноведения и инженерной графики.

Евгений Юрьевич Просвиряков (к.ф.-м.н.; [email protected]), ведущий научный сотрудник, лаб. прикладной механики.