Two-dimensional convection of an incompressible viscous fluid with the heat exchange on the free border

Full Text

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: V l a s o v a S. S., P r o s v i r y a k o v 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 (sveta040790@yandex.ru; Corresponding Author), Postgraduate Student, Dept. of Machinery Engineering & Engineering Graphics. Evgeny Yu. Prosviryakov (Cand. Phys. & Math. Sci.; evgen_pros@mail.ru), Leading Researcher, Laboratory of Applied Mechanics. 567 V l a s o v a S. S., P r o s v i r y a k o v E. Yu. 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: Vx ∂Vx ∂Vx ∂P ∂ 2 Vx ∂ 2 Vx + Vz =- +ν + , ∂x ∂z ∂x ∂x2 ∂z 2 ∂Vz ∂Vz ∂P Vx + Vz =- + ν∆Vz + gβT, ∂x ∂z ∂z ∂Vx ∂Vz + = 0, ∂x ∂z ∂T ∂T ∂2T ∂2T . Vx + Vz =χ + ∂x ∂z ∂x2 ∂z 2 (1) 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 ρ; T is a deviation from the average temperature; ν, χ, β are the dissipative coefficients of the kinematic viscosity, heat diffusivity and thermal Figure 1. Schematic view of the model problem 568 Two-dimensional convection of an incompressible viscous fluid. . . expansion of fluid, respectively; g is the acceleration of gravity; ∆= ∂2 ∂2 + ∂x2 ∂z 2 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 T = T0 (z) + T11 (z) , P = P0 (z) + P11 (z) . 2 2 (2) 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. β 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 , T11 , P0 : du dw d2 u = u2 + w + P11 , = -u, ∂z 2 ∂z ∂z dP0 d2 w dw = ν 2 + gβT0 - w , ∂z ∂z ∂z dP11 dT0 d2 T0 = gβT11 , w = χ T11 + , dz dz dz 2 dT11 d2 T11 2uT11 + w =χ . dz dz 2 ν (3) 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: Θ is for the temperature, Vz and Vx are for 3 4 and gβΘh the velocities gβΘh respectively. Thus, the accounting of geometrical νl νl2 anisotropy of the task results in the formation of a flow and availability of two scales for the velocities, since gβΘh4 gβΘh3 = δ ; νl2 νl the pressure for an incompressible fluid divided by the constant density is gβΘh; δ = h/l. We write down the system (3) in the dimensionless form: dP11 dw dT11 1 d2 T11 = T11 , = -u, Gr δ 2 2uT11 + w = dz dz dz Pr dz 2 2T 2 dT 1 d dw dP 0 0 0 2 2d w Gr δ 2 w = δ2 + , Gr δ w = - + δ + T0 , dz Pr dz 2 dz dz dz 2 du d2 u Gr δ 2 u2 + δ = -P11 + δ 2 2 . dz dz (4) 569 V l a s o v a S. S., P r o s v i r y a k o v E. Yu. 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 gβΘh5 Gr = ν 2 l2 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 δ 2 1. In this case, the received exact solutions are fair with any Grashof number they satisfy the assessed value Gr 1/δ 2 where Gr ∈ (0; 1/δ 2 ). 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: d2 T11 dP11 d2 u dw = 0, = T , = P11 , = -u, 11 2 2 dz dz dz dz 2 d2 T0 dP0 2 2d w = -δ T , = T + δ . 11 0 dz 2 dz dz 2 (5) 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: T11 = C1 z + C2 , P11 = C1 z2 + C2 z + C3 , 2! z4 z3 z2 + C2 + C3 + C4 z + C5 , 4! 3! 2 z5 z4 z3 z2 (6) w = -C1 - C2 - C3 - C4 - C5 z + C6 , 5! 4! 3! 2! z2 z3 T0 = -C1 δ 2 - C2 δ 2 + C7 z + C8 , 3! 2! 4 z3 z2 2z P0 = -2C1 δ - 2C2 δ 2 + C7 - δ 2 C3 + C8 - δ 2 C4 z + C9 . 4! 3! 2 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 ) 2 with dimensionless variables, it is written as: u = C1 1 T = (1 - x2 ). 2 The boundary conditions for defining the constants of integration of dimensionless solutions (5) can be written as: z = -1 : w = u = 0, 570 T11 = -1, 1 T0 = ; 2 (7) Two-dimensional convection of an incompressible viscous fluid. . . z = 0 : w = 0, P11 = 0, P0 = 1, dT0 = - Bi T0 , dz dT11 = - Bi T11 . (8) 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 = C6 = 0, 1 Bi , C2 = - , C4 = C1 = Bi + 1 Bi + 1 δ 2 Bi(Bi + 3) - 3 Bi(Bi + 1) C7 = , 3!(Bi + 1)2 C9 = 1. 2(4 Bi +15) 3 Bi +10 , C5 = , 5!(Bi + 1) 5!(Bi + 1) 3 (Bi + 1) - δ 2 (Bi + 3) C8 = . 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: Vx = x Vz = - T =- P =- z3 2(4 Bi +15)z 3 Bi +10 Bi z 4 - + + , 4!(Bi +1) 3!(Bi +1) 5!(Bi +1) 5!(Bi +1) z4 (4 Bi +15) z 2 (3 Bi +10) z Bi z 5 + - - , 5!(Bi +1) 4!(Bi +1) 5!(Bi +1) 5!(Bi +1) Bi δ 2 z 3 δ2z2 δ 2 Bi (Bi +3) - 3 Bi (Bi +1) + + z+ 3!(Bi +1) 2!(Bi +1) 3!(Bi +1)2 3(Bi +1) - δ 2 (Bi +3) x2 (Bi z - 1) + + , 3!(Bi +1)2 2!(Bi +1) 2 δ 2 Bi (Bi +3) - 3 Bi (Bi +1) z 2 2 Bi δ 2 z 4 2δ 2 z 3 + + + 4!(Bi +1) 3!(Bi +1) 4! (Bi +1)2 + 3 (Bi +1) - δ 2 (Bi +3) 2 (4 Bi +15) δ 2 - z+ 5!(Bi +1) 3! (Bi +1)2 Bi z 2 z x2 +1+ - . 2!(Bi +1) 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: u = (z + 1)f (z) = (z + 1) Bi z 3 (Bi + 4)z 2 - + 4!(Bi +1) 4!(Bi +1) 571 V l a s o v a S. S., P r o s v i r y a k o v E. Yu. + (Bi +4)z 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 ∈ (-∞; -10/3) ∪ (-5/2; -1) ∪ (-1; +∞) . When Bi = 0, the assumed cubic polynomial degenerates into a linear function. Two solutions are possible in the one case with Bi ∈ (-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: ψ = -xz(z + 1)2 - Bi z 2 (2 Bi +5)z 3 Bi +10 + - 5!(Bi +1) 5!(Bi +1) 5!(Bi +1) and vorticity: Ωy = x Bi z 3 z2 2(4 Bi +15) - + . 3!(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.0 0.000024 0.000144 0.000048 0.000288 z 0.0 0.000408 0.000336 0.000288 0.000672 0.000576 0.00024 0.000528 0.000528 0.000816 0.00012 0.000768 0.000264 0.000744 0.000456 0.000552 0.00096 0.000168 0.000984 0.000864 0.001056 0.2 0.000792 0.000456 0.000768 0.000936 0.000696 0.000696 0.000984 0.00060.000384 0.001032 0.000888 0.000336 0.001104 0.000576 0.000912 0.000792 0.00096 0.001008 0.000816 0.001008 0.00108 0.00084 0.000648 0.001032 0.000648 0.00048 0.000192 0.000144 0.000432 0.000936 0.000864 0.00084 0.00072 0.000504 0.000672 0.000912 0.000264 0.4 0.000888 0.000552 0.000312 0.000744 0.00072 0.000480.0006 0.000624 0.000408 0.000312 0.000624 0.000432 0.000072 0.000504 0.00036 0.000216 0.00036 0.00024 0.00012 0.000216 0.000384 0.000072 0.6 0.0368 0.0368 0.01760.024 0.0032 0.04 0.032 0.0288 0.0272 0.0192 0.0096 0.0336 0.0112 0.0384 0.0336 0.0224 0.0352 0.0384 0.0288 0.0144 0.0304 0.0256 0.0208 0.0352 0.0144 0.0032 0.0304 0.032 0.008 0.2 0.0208 0.0272 0.024 0.016 0.0224 0.008 0.01760.0128 0.0096 0.0048 0 0.0016 0.0016 0.6 0.008 0 0.008 0.0112 0.0096 0.0064 0.0032 0.0032 0.8 0.0112 0.000024 0.5 0.0 x Figure 2. 0.0096 0.000024 0.0016 0.0048 0.0064 0.0096 0.008 1.0 0.0112 0 1.0 0.0048 0.0048 0.0064 0.8 1.0 0.0192 0.0128 0.0064 0.0016 z 0 572 0.0016 0.0064 0.40.0112 0.0048 0.000192 0.000096 0.000048 0.000024 0.000168 0.000096 0.016 0.0256 0.5 1.0 1.0 0 0.0032 0.5 0.00640.008 0.0048 0.0032 0.0016 0.00960.0112 0.0 0.5 1.0 x Isolines of the current function (left) and the function of vorticity (right) with Bi = -2.7 Two-dimensional convection of an incompressible viscous fluid. . . 0.0 0.00011 0.00077 0.00033 0.0 0.00121 0.000220.00066 zz 0.00176 0.0033 0.00264 0.0033 0.00407 0.00231 0.00286 0.001540.00066 0.20.004070.00352 0.00198 0.00143 0.00363 0.0044 0.00451 0.00418 0.00231 0.00055 0.00319 0.00484 0.00418 0.00275 0.00132 0.00374 0.00462 0.00462 0.00308 0.00209 0.00374 0.00495 0.00429 0.00198 0.00132 0.0044 0.00473 0.00396 0.00341 0.00308 0.00451 0.40.00484 0.00242 0.002420.00352 0.00473 0.00385 0.00429 0.00385 0.00088 0.00011 0.00165 0.00275 0.00143 0.00396 0.00319 0.00253 0.00286 0.00088 0.00363 0.00176 0.00297 0.00297 0.00187 0.00341 0.6 0.0022 0.00099 0.00253 0.00154 0.00187 0.0022 0.00099 0.00044 0.00165 0.00121 0.00033 0.0011 0.8 0.00044 0.0011 0.00077 0.00055 0.00022 1.0 0 1.0 0 0.5 0.0 xx 0.5 1.0 0.1035 0.09 0.0495 0.0045 0.063 0.0765 0.099 0.108 0.0855 0.002640.00209 0.099 0.1125 0.0855 0.0360.0540.072 0.018 0.036 0.018 0.20.09450.0810.06750.054 0.1035 0.09 0.108 0.0765 0.0585 0.0405 0.0945 0.081 0.063 0.045 0.009 0.072 0.0585 0.0450.03150.0225 0.027 0.4 zz 0.0675 0.0495 0.0315 0.009 0.0405 0.027 0.0135 0.0135 0.0225 0.0045 0.6 0 0 0.0045 0.8 0.0225 0.04050.0315 0.0135 0.0225 0.0135 0.018 0.009 0.045 0.036 0.027 0.018 1.00.0495 1.0 0.5 0.036 0.0315 0.045 0.009 0.0045 0.0 xx 0.027 0.04050.0495 0.5 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 + 21)/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]: 5z 4 + 20z 3 - 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 δ. 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 f (z) = z(z + 1)2 - Bi z 2 (2 Bi +5)z 3 Bi +10 + - . 5!(Bi +1) 5!(Bi +1) 5!(Bi +1) 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 ∈ (-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. 573 V l a s o v a S. S., P r o s v i r y a k o v E. Yu. 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 z ∈ [-1; 0] takes negative values with Biot Bi ∈ [-1; 0). The received T11 solution takes up the constant value when Bi = 0: T11 = -1. The additive component is monotonously increasing with Bi ∈ (-∞; -1) ∪ (0; +∞), otherwise it is monotonously decreasing. Function T11 vanishes in the interval of the definition in the point z = 1/Bi with Bi ∈ [-1; 0). Now let’s consider the properties of the background temperature and pressure. The background temperature can be presented in the following form: T0 = Bi C2 δ 2 z3 z2 - C2 δ 2 - Bi C8 z + C8 . 3! 2! The existence of an odd number of solutions is equivalent to an inequality being solved: C2 δ 2 T0 (-1)T0 (0) = C8 C8 (Bi + 1) - (Bi + 3) < 0. 3! Due to the boundary conditions (7) and (8), the temperature can be presented in another form: T0 = 1 (z + 1) + -z 2 Biδ 2 (Bi +1) + zδ 2 (Bi +1)(Bi +3)- 2 3!(Bi + 1)2 - 3 Bi(Bi +1) - δ 2 (Bi +3) = (z + 1) 1 g (z) . = + 2 3!(Bi +1)2 Let’s investigate the function g(z). The function g(z) is a square polynomial function depending on two parameters: δ and Bi. The function g(z) takes a zero value in the interval of the definition when an inequality is being solved: g(-1)g(0) = (3Bi + 3Bi2 + 3δ 2 + Biδ 2 )(3Bi + 3Bi2 + 6δ 2 + 6Biδ 2 + 2Bi2 δ 2 ) < 0. The function g(z) 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 T0 is not more than 0.5, the number of zeros of the function T0 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 δ = 1 (Fig. 4) and when Bi = -0.4 and δ = 0.48 (Fig. 5). It is obvious that the square law change of P11 converts into a linear dependence P11 = -z when Bi = 0. The function P11 reduces to zero within the interval [-1; 0) in the point z = 2/ Bi. Since the square curve turns into a linear dependence, the function P11 monotonously decreases. Let’s consider the background pressure function P0 (z). The function P0 (z) is a quartic polynomial depending on two parameters: δ and Bi. The function 574 Two-dimensional convection of an incompressible viscous fluid. . . 0.0 0.259 0.481 0.592 0.2 0.444 0.222 0.185 z 0.6 0.074 0.074 0.148 0.8 0.222 1.0 0 1.0 0.5 0.37 0.5 1.116 1.0665 0.6 1.1295 0.0 1.0575 1.0 1.0 1.0935 x 1.044 1.062 1.053 0.0 0.5 1.0215 1.0755 1.1205 1.0 1.035 1.107 1.098 0.037 0.185 1.08 1.125 1.1115 1.134 1.0845 0.407 0.444 1.1025 1.089 1.0305 0.8 0.333 0.296 0.259 1.0485 0.037 0.111 1.0035 1.026 1.0395 1.071 1.008 0.4 0.148 0.111 1.0125 1.017 0.2 0.296 0.4 z 0.333 0.37 0.0 0.407 0.555 0.629 0.518 0.5 1.0 x Figure 4. Isolines of the functions of temperature and pressure with Bi = -2 and δ = 1 z 0.0 0.208 0.104 0.26 0.078 0.39 0.26 0.2 0.1560 0.182 0.234 0.338 0.104 0.494 0.078 0.4 0.208 0.416 0.182 0.026 0.13 0.312 0.6 0.13 0.026 0.234 0.442 0.8 0.364 0.052 0.156 1.0 0.052 1.0 0.468 0.286 0.5 0.468 0.494 0.0 0.01.008 0.442 0.26 0.104 0.13 0.26 0.026 0.182 0.338 0.156 0.234 0.0520.104 0.208 0.234 0.364 0.052 0.13 0.156 0.286 0.026 0.39 0.182 0.078 0.078 0.312 0.416 0.5 0.208 0 1.0 1.008 0.987 0.945 0.2 0.924 1.05 1.092 0.966 0.4 1.05 0.819 1.092 z 0.84 0.6 1.134 1.029 0.798 1.071 0.861 1.0 1.155 1.0 x 1.134 0.714 1.071 1.113 0.8 0.903 0.756 0.777 0.735 0.672 0.651 0.567 0.5 0.609 0.63 0.546 0.588 0.693 0.0 0.5 0.882 1.113 1.029 1.155 1.0 x Figure 5. Isolines of the function of temperature pressure with Bi = -0.4 and δ = 0.48 P0 (z) takes one zero value in the interval of the definition z ∈ [-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 δ and Bi in the considered interval. The existence of the extremum points in the interval of the definition z ∈ [-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 δ 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 NavierStokes 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. 575

About the authors

Svetlana S Vlasova

Kazan National Research Technical University named after A. N. Tupolev

Email: sveta040790@yandex.ru
(sveta040790@yandex.ru; Corresponding Author), Postgraduate Student, Dept. of Machinery Engineering & Engineering Graphics 10, Karl Marx st., Kazan, 420111, Russian Federation

Evgeny Yu Prosviryakov

Institute of Engineering Science, Ural Branch of RAS

Email: evgen_pros@mail.ru
(Cand. Phys. & Math. Sci.; evgenpros@mail.ru), Leading Researcher, Laboratory of Applied Mechanics 34, Komsomolskaya st., Ekaterinburg, 620049, Russian Federation

References

Supplementary files