CC BY
31
UDC 517.95: 519.63
Lamtyugova S. N.1, Sidorov M. V.2, Sytnykova I. V.3
1Ph.D in Physics and Mathematics, Assistant of the Department of Advanced Mathematics of O. M. Beketov National University of
Urban Economy in Kharkiv, Kharkiv, Ukraine 2Ph.D in Physics and Mathematics, Associate Professor, Associate Professor of the Department of Applied Mathematics of Kharkiv
National University of Radio Electronics, Kharkiv, Ukraine 3Ph.D in Pedagogics, Associate Professor of the Department of Advanced Mathematics of O. M. Beketov National University of
Urban Economy in Kharkiv, Kharkiv, Ukraine
METHOD OF NUMERICAL ANALYSIS OF THE PROBLEM OF STATIONARY FLOW PAST BODIES OF REVOLUTION BY VISCOUS
FLUID
Context. The nonlinear stationary problem of flow past a body of revolution by a viscous incompressible fluid is examined in this article.
Objective. The purpose of this work is to develop a new method of numerical analysis of stationary problems of flow around bodies of revolution by viscous incompressible fluid.
Method. The mathematical model of the process under consideration is a nonlinear boundary value problem for the stream function obtained by the transition from the system of Navier-Stokes equations to one nonlinear equation of the fourth order. A special feature of the formulation the task of the flow past body is that the boundary value problem is considered in an infinite region and both boundary conditions on the boundary of the streamlined body and the condition at infinity are imposed for the stream function. Using the structural method (the ^-functions method), the task solution structure, that exactly satisfies all the boundary conditions of the task, and also guarantees the necessary behavior of the stream function at infinity, is constructed. Two approaches are proposed to approximate the uncertain components of the structure. The first approach is based on the use of the successive approximations method, which makes it possible to reduce the solution of the initial nonlinear task to the solution of a sequence of linear boundary value problems. These linear tasks are solved by the Bubnov-Galerkin method at each step of the iteration process. The second approach for approximating the uncertain components of the structure is based on the usage of the nonlinear Galerkin method and it is proposed to use it in the case of divergence of successive approximations. In this case, the solution of the initial nonlinear task reduces to solving a system of nonlinear algebraic equations.
Results. A computational experiment was carried out for the task of flow past a sphere, an ellipsoid of rotation and two articulated ellipsoids for various Reynolds numbers.
Conclusions. The conducted experiments have confirmed the efficiency of the proposed method of numerical analysis of stationary problems of flow around bodies of revolution by viscous incompressible fluid. The prospects for further research may consist in using the method developed for the implementation of semi-discrete and projection methods for solving non-stationary problems.
Keywords: flow past bodies, stationary flow, viscous fluid, stream function, ^-functions method, successive approximations method, Galerkin method.
NOMENCLATURE
dQ - boundary of a streamlined body;
E - Stokes operator in a spherical coordinate system;
F - class of functions which have generalized derivatives up to the second order inclusive and are quadratically summable together with the derivatives with respect to Qj; on the boundary dQ they meet to the homogeneous boundary conditions of the corresponding task;
Jn (Z) - Gegenbauer functions of the first kind;
n - outer normal to dQ;
r, 6 , 9 - variables of the spherical coordinate system; Re - Reynolds number;
u - solution of the auxiliary task with homogeneous boundary conditions;
4°), , un), - - sequence of approximations to the
generalized solution of the task for u , which is considered in the finite domain;
u* - generalized solution of the task for u , which is considered in the finite domain;
© Lamtyugova S. N., Sidorov M. V., Sytnykova I. V., 2018
DOI 10.15588/1607-3274-2018-1-6
Ux - unperturbed fluid velocity at infinity;
v = Re-1 - viscosity factor;
{9k (r, 6)} - a complete system of particular solutions
of the equation e 2U = 0 with respect to the exterior of a sphere of finite radius;
{t j (r, 6)} - a complete system of particular solutions of the equation e2u = 0 relative to the domain {ra(r,6) < M};
(r, 6)} - sequence of functions that is complete concerning whole plane;
®1, ® 2 - undefined components of the solution structure;
y = y(r, 6) - stream function;
ra - a sufficiently smooth function, built with the help of a constructive apparatus of the ^-functions theory;
ra = 0 - normalized equation of dQ ;
Q - flow domain;
Qn - bounded domain {0 < ra(r, 6) < Mn } ;
Q1 - any finite part of Q .
INTRODUCTION
The phenomena, which are observed in the atmosphere and the hydrosphere, the problems of hydrodynamics, aerodynamics, heat power engineering, chemical kinetics, and biomedicine, can be often researched within the framework of an incompressible viscous fluid model. The tasks, which represent the practical interest, as a rule, are described by the nonlinear Navier-Stokes equations. In the case, when the task possesses symmetry properties and can be reduced to a two-dimensional one, it is convenient to enter the stream function instead of components of the fluid velocity. Methods for solving the external tasks (problems of flow past bodies) for equations with respect to the stream function have not been developed sufficiently, which is due to the high order and nonlinearity of the equations, as well as the limitlessness of the domain in which the equations are considered. The structural R-functions method, which has been used in this work, allows to construct the so-called structures of the solution of boundary value problems are bundles of functions which exactly satisfy all the boundary conditions of the task. The geometry of the area is taken into account exactly. Approximating further the undefined components of the solution structure by some projection method, we obtain an approximate solution in an analytical form, which facilitates its further usage for finding various flow characteristics. Thus, the development of new, as well as the improvement of existing methods of mathematical modeling of stationary problems of flow past bodies by viscous incompressible fluid using the R-function method, is a scientific issue of current interest.
The object of the investigation is the stationary hydrodynamic process of flow past bodies by a viscous incompressible fluid, described by a nonlinear equation with respect to the stream function.
The subject of the study is a mathematical model of the stationary problem of flow past a body by viscous incompressible fluid and the method of its numerical analysis.
The purpose of the work is to develop a new method of mathematical simulation of stationary flow past bodies of revolution by viscous incompressible fluid on the basis of the R-function method.
1 PROBLEM STATEMENT
We consider the nonlinear steady-state problem of flow past body of rotation by viscous incompressible fluid in a spherical coordinate system [1, 2]:
Z72 1
vE y = -
' sin 0
5y 5Ey 5y 5Ey^ 50 dr Hr 50 J +
r 2 sin 0
f 2ctg0^-2 40Ï Eymn ,
dr r 50
I =o,
ISQ ' 3n
= 0,
dQ
lim y • r 2 =1 Ux sin2 0 .
r 2
(1) (2) (3)
52 y sin 0 5 f 1 5y where Ey =-— + -
dr
2
2 50 I sin 0 50
2 REVIEW OF THE LITERATURE
In general, the methods of finite differences and finite elements are used for numerical analysis of hydrodynamic problems. The main advantage of these methods is the relative simplicity of implementation, however, when moving to a new domain (especially complex geometry), it is necessary to generate again and adjust the calculated grid under the task. Furthermore, the demolition of the conditions on infinity on a certain contour, which is located at a finite distance from the streamlined body, gives the additional errors in the approximate solution. Using the R functions method [3, 4] of V. L. Rvachev, the Academician of Ukrainian National Academy of Sciences, it is possible to take precisely into account the geometric and analytical information included in the statement of the problem.
The R-functions method in hydrodynamics problems was used in [5-10], but the tasks for calculating steady flows of an ideal fluid, stationary and nonstationary flows of a viscous fluid in limited domains or viscous liquids in the presence of helical symmetry were considered. The tasks of flow past bodies of rotation by a stationary flow of a viscous fluid using the R-function method were solved in [11-16]. The task of calculating the external slow flow past bodies by a viscous incompressible fluid (Stokes linearization) is solved in [11-13]. Usage of R-functions method, successive approximations and Galerkin-Petrov method to calculate axisymmetric stationary flows of a viscous incompressible fluid was proposed in works [14, 15]. The task of mass transfer of a body of rotation with a uniform translational flow was considered in [16].
In this work the R-functions method, successive approximations and Galerkin method are proposed to apply for mathematical modeling of nonlinear stationary task of flow past body of revolution by a viscous incompressible fluid.
3 MATERIALS AND METHODS
For solving the task (1)-(3), we introduce the function
[9]
fM( x) =
Mx
1 - exp-, 0 < x < M ;
x - M
1, x > M (M = const > 0), which meets the following conditions:
Isq
1) raM > 0 in Q ; 2) = 0;
3)
5ra
M
5n
= -1; 4) raM = 1, if raM > M.
dQ
In the task (1)-(3) we will make a replacement y = U0 + u , where U0 is solution of the linear task (Stokes approximation):
E 2u0 = 0 in Q,
(4)
r
r
^ lôQ
= 0
du 0
5n
=0
dQ
lim r 2uo =1 Ux sin2 6.
r^+x 2
(5)
(6)
Method for solving the task (4)-(6) was developed in the works [11-13].
For function u we obtain a boundary value problem with homogeneous boundary conditions
Z72 1
vE u = -
r 2 sin 6
d(uo + u) dE(uo + u)
56
d(uo + u) dE (uo + u)
dr
56
r2 sin6
5r
2ctg6
d(uo + u) dr
2 d(uo + u) '7 56
E(uo + u) in Q .
I a dM
uL~> = o, —
l5Q 5n
=o
dQ
_2
lim r u = o .
r^+x
(7)
(8)
(9)
We use the method of successive approximations to solve task (7)-(9). If the initial approximation is given u^^
and k approximation un) is built, then the new (k +1)
approximation u^k+1) will be found as a solution of the linear problem
vE Vk+1 =
r 2 sin 6
5(uo + u{k)) dE(uo + u^))
(k h
56 dr
5(uo + u{k)) dE(uo + u{k)) 1
5r
r2 sin6
2ctg6
56
d(uo + ujk )) dr
2 5(uo + u^)) "
56
(k+1)
E(uo + u^k)) in Qn du^k+1)
dQ
= o.
5n
= o.
(io)
(11)
dQ
continued by zero outside Qn .
A theorem on the convergence of successive approximations (10)—(11) at small Reynolds numbers for each
n in the energy norm to a unique generalized solution u* e F of the task (7)-(9), which is examined in a bounded domain Qn , was formulated and proved in work [15].
According to the Bubnov-Galerkin method, at each iteration in the domain Qn, we will seek the approximate solution of task (10)—(11) in the form
(k+1) = n.N
N
= 1' j=1
(k+1)
'n, j Yj
* j
(12)
Numbers cn^j1"1, j = 1'. .' N, are solution of the system
of linear algebraic equations
N
X c«k+1)v([* j, * ] + (K* j, (■ )) = ( f(k+1), * ), i = 1, N
= ( f (k+1)
j=1 where
k = _aî ^ctg6A+1a
^ r 2 56 r dr
2ctg6 5 + 2 _^1Îa_ 2ctg6 5 2 5
,2 56 r dr
,2 56 r dr
(f(k+1),*i) = J E(uo + unk))Êk^-Q dr 56
Qn
_ J E(uo + u?)) d(uo + un))^drd6 + Q o n 56 dr
n
h J 2ctg6 d(uo + ^)) E(uo + u^))*idrd6_ J Pi"
_ j 2 d(uo + )) E(uo + unk))*;drd6.
Q
r 56
We construct a system of coordinate functions {(■}.
For this we use a complete system of particular solutions of 2
the equation E u = o [17] and the structure method of ^-functions [3, 4].
In work [13] it was proved that the function of the form
u = a>M °1 + roM (1 _roM )® 2.
(13)
is the structure of the solution of the boundary value task (7)-(9), i.e. for any choice of sufficiently smooth functions
®1 and ®2 ( ®1 • r_2 ^ o at r ^ +x ) it exactly satisfies the boundary conditions (8) and the condition at infinity (9).
Approximations of the indeterminate components ®1 and ® 2 of structure (13) will be sought in the form
®1 1 = Za k ^k, ® 2
k=1
m2 =
= ZP j=1
y t j,
1
Q
n
1
r
where
(r,6)} = {r1-J (cos6), k = 2,3,
r3~kJk(cos6), k = 4,5,...},
{t (r, 6)} = {rJ2 (cos 6), J3 (cos 6),
rJJj (cos6), rJ+2Jj (cos6), j = 2,3,...}.
Then, the complete sequence of functions concerning whole plane has the form:
{^ (r, 6)} = {raM (r, 6)9k (r, 6),
raM (r,6)(1 (r,6))tj (r,6)}.
In case of divergence of successive approximations we will use nonlinear Galerkin method. An approximate solution M«,W of the task (7)-(9) which is considered in the finite domain Q n is sought in the form (12). As a result, we obtain a system of nonlinear equations, each of which is a quadratic function with respect to cn, j. As the initial approximation, we choose a set cn, j , which corresponding to the solution of the linear task (Stokes approximation), or, for large Reynolds numbers, to the solution which obtained for smaller Reynolds numbers.
4 EXPERIMENTS
A computational experiment was carried out for the tasks of flow past a sphere, ellipsoids of rotation with different semi-axes, and two articulated ellipsoids. The double integrals in the systems for the determination cn, j were
calculated approximately by the Gauss' formula with 50 nodes for each variable. For the method of successive
approximations, calculations were ceased at e = 10"6. It was experimentally established that for Re > 10 the successive approximations diverge. In this case, the nonlinear Galerkin method was used to obtain the approximate solution.
Calculations for the tasks of flow past a sphere
x2 + y2 + z2 = 1 at Ux = 1, M = 10, m1 = 10, m2 = 14,
2 2 X2 + y2 + z 2 = ,
Re = 10; 20; 25 , ellipsoid of revolution ^2 _ 1
at Ux = 1, M = 10, m1 = 10, m2 = 14, Re = 10; 30; 50 , two articulated ellipsoids, limited by surfaces
(X -1)2 + y 2 + z 2 = 1, (X +1)2 + y2 + z 2 = 1, at U „ = 1,
2
22
M = 10, m1 = 6, m2 = 10, Re = 5; 10; 20, were performed. 5 RESULTS
Fig. 1 shows the streamlined contours of the obtained approximate solution for the task of flow past a sphere. The detailed pictures of streamlined contours and vector velocity field behind the sphere are shown in Fig. 2.
Fig. 3 shows the streamlined contours of the obtained approximate solution for the task of flow past an ellipsoid of rotation. The detailed pictures of streamlined contours and vector velocity field behind the ellipsoid are shown in Fig. 4.
Fig. 5 shows the streamlined contours of the obtained approximate solution for the task of flow past two articulated ellipsoids. The detailed pictures of streamlined contours and vector velocity field behind ellipsoids and in the hollow between them are shown in Fig. 6.
Figure 2 - Detailed pictures of streamlined contours and vector velocity field behind the sphere:
a - Re = 10. b- Re = 20. c- Re = 25
Figure 3 - Streamlined contours of obtained approximate solution for flow past an ellipsoid of rotation: a - Re = 10, b - Re = 30, c - Re = 50
Figure 4 - Detailed pictures of streamlined contours and vector velocity field behind the ellipsoid: a - Re = 10, b - Re = 30, c - Re = 50
Figure 5 - Streamlined contours of obtained approximate solution for flow past two articulated ellipsoids:
a - Re = 5 , b - Re = 10 , c - Re = 20
Figure 6 - Detailed pictures of streamlined contours and vector velocity field behind ellipsoids and in the hollow between them:
a - Re = 5, b - Re = 10, c - Re = 20
6 DISCUSSION
For small Reynolds numbers, the flow is symmetrical, without the formation of a detachment zone in the aft domain of the body. With an increase in the Reynolds number the flow pattern changes: the secondary vortices appear behind the body and their size and intensity increase.
In the task of flow past a sphere, it is found that the secondary vortices behind the body appear with an increase in the Reynolds number up to « 20 - 25. For the task of flow past an ellipsoid of rotation with the ratio of the semi-a
axes of an ellipsoid — > 2, the secondary vortices behind
the body do not appear independently of the Reynolds number, but appear with a decrease in the ratio of the semi-axes. For two articulated ellipsoids, a computational experiment showed that as the Reynolds number increases up to «10, vortices appear behind the body and in the hollow at the junction of the ellipsoids.
The obtained results correlate with the known results of physical experiments [18] and the results obtained by other authors [19, 20], which indicates the effectiveness of the developed numerical method.
CONCLUSIONS
A new numerical method for the calculating the external flows of viscous incompressible fluid is proposed. The proposed method is based on the joint application of the R-functions method, successive approximations and Galerkin method. The algorithm of the method does not change when the geometry of the domain is changed, and the structure of the solution accurately takes into account both the boundary conditions on the boundary of the streamlined body and the condition at infinity. It is experimentally established that the iterative process begins to diverge at Reynolds numbers Re > 10. In this case, to approximate the undetermined components of the structure of the solution of the nonlinear stationary problem of flow past a body of revolution by a viscous incompressible fluid in a spherical coordinate system the nonlinear Galerkin method was applied. For various Reynolds numbers, the stationary problem of flow past a body of rotation in a spherical
coordinate system for a sphere, ellipsoid, and two articulated ellipsoids has been solved numerically. For each body Reynolds numbers, at which secondary vortices arise behind the body, are defined.
ACKNOWLEDGMENTS
The work was carried out at the department of Applied Mathematics at Kharkiv National University of Radio Electronics and department of Advanced Mathematics at
0. M. Beketov National University of Urban Economy in Kharkiv within the framework of collaborative scientific research conducted by the departments.
REFERENCES
1. Лойцянский Л. Г. Механика жидкости и газа / Л. Г. Лойцянс-кий. - М. : Дрофа, 2003. - 840 с.
2. Химическая гидродинамика : справочное пособие / [А. М. Кутепов, А. Д. Полянин, З. Д. Запрянов, А. В. Вязьмин, Д. А. Казенин]. - М. : Квантум, 1996. - 336 с.
3. Рвачев В. Л. Теория ^-функций и некоторые ее приложения / В. Л. Рвачев. - К.: Наук. думка, 1982. - 552 с.
4. Shapiro V. Semi-analytic geometry with ^-functions / V. Shapiro // Acta Numerica. - 2007. - V. 16. - P. 239-303.
5. Колосова С. В. Применение проекционных методов и метода ^-функций к решению краевых задач в бесконечных областях: дисс. ... к. ф.-м. н.: 01.01.07 - Вычислительная математика / С. В. Колосова. - Харьков: ХИРЭ, 1972. - 85 с.
6. Колосова С. В. Применение метода ^-функций к расчету плоских течений вязкой жидкости / С. В. Колосова, М. В. Сидоров // Вюник Харювського нащонального ушверситету. Сер. Прикл. матем. i мех. - 2003. - № 602. - С. 61-67.
7. Суворова И. Г. Компьютерное моделирование осесимметрич-ных течений жидкости в каналах сложной формы / И. Г. Суворова // Вестник национального технического университета «ХПИ». - 2004. - № 31. - С. 141-148.
8. Тевяшев А. Д. Об одном подходе к математическому моделированию плоских стационарных течений вязкой несжимаемой жидкости в конечных односвязных областях / А. Д. Тевя-шев, Н. В. Гибкина, М. В. Сидоров // Радиоэлектроника и информатика. - 2007. - № 2 (37). - С. 50-57.
9. Максименко-Шейко К. В. Математическое моделирование теплообмена при движении жидкости по каналам с винтовым типом симметрии методом ^-функций / К. В. Максименко-Шейко // Доповда Нащонально!' академп наук Украши. - 2005. -№ 9. - С. 41-46.
10. Artyukh A. Mathematical modeling and numerical analysis of nonstationary plane-parallel flows of viscous incompressible fluid by ^-functions and Galerkin method / A. Artyukh, M. Sidorov // Econtechmod. - 2014. - Vol. 3, No 3. - P. 3-11.
11. Ламтюгова С. Н. Математическое моделирование линеаризованных задач обтекания в сферической и цилиндрической системах координат / С. Н. Ламтюгова // Вюник Заж^зького нацюнального ушверситету. Сер. Фiзико-математичнi науки. -2012. - № 1. - С. 112-122.
12. Ламтюгова С. М. Застосування методу Л-функцш до розра-хунку зовшшшх повшьних течш в'язкм рщини / С. М. Ламтюгова, М. В. Сидоров // Вiдбiр i обробка iнформацiï. - 2012. -№ 36 (112). - С. 56-62.
13. Lamtyugova S. N. Numerical analysis of the external slow flows of a viscous fluid using the ^-function method / S. N. Lamtyugova, M. V. Sidorov // Journal of Engineering Mathematics. - 2015. -Vol. 91, No. 1. - P. 59-79. D0I:10.1007/s10665-014-9746-x
14. Колосова С. В. Применение итерационных методов к решению внешних задач гидродинамики / С. В. Колосова, С. Н. Ламтюгова, М. В. Сидоров // Радиоэлектроника и информатика. - 2012. - № 3. - С. 13-17.
15. Ламтюгова С. Н. Применение итерационных методов к расчету обтекания тел стационарным потоком вязкой жидкости / С. Н. Ламтюгова // Радиоэлектроника и информатика. - 2015. -№ 2. - С. 49-56.
16. Колосова С. В. Об одном методе численного анализа вязких течений, усложненных массообменом (задача обтекания) / С. В. Колосова, С. Н. Ламтюгова, М. В. Сидоров // Радиоэлектроника и информатика. - 2014. - № 1 (64). - С. 25-30.
17. Полянин А. Д. Справочник по линейным уравнениям математической физики / А. Д. Полянин. - М. : Физматлит, 2001. -576 с.
18. Ван-Дайк М. Альбом течений жидкости и газа / М. Ван-Дайк. -М. : Мир, 1986. - 184 с.
19.Batchelor G. K. An introduction to fluid dynamics / G. K. Batchelor. - Cambridge : Cambridge University Press, 1967. -615 p.
20. Taamneh Y. CFD Simulations of Drag and Separation Flow Around Ellipsoids / Y. Taamneh // Jordan Journal of Mechanical and Industrial Engineering. - 2011. - Vol. 5, No. 2. - P. 129-132.
Article was submitted 16.07.2017.
After revision 05.09.2017.
Ламтюгова С. М.1, Сидоров М. В.2, Ситникова Ю. В.3
'Канд. фiз.-мат. наук, асистент кафедри вищо! математики Харгавського нацюнального ушверситету мюького господарства iменi О. М. Бекетова, Харгав, Украша
2Канд. фiз.-мат. наук, доцент, доцент кафедри прикладно! математики Харгавського нацюнального ушверситету радюелектрошки, Харгав, Украша
3Канд. пед. наук, доцент кафедри вищо! математики Харювського нацюнального ушверситету мюького господарства iменi О. М. Бекетова, Харгав, Украша
МЕТОД ЧИСЕЛЬНОГО АНАЛ1ЗУ ЗАДАЧ1 ОБТ1КАННЯ Т1Л ОБЕРТАННЯ СТАЦ1ОНАРНИМ ПОТОКОМ В'ЯЗКО1 Р1ДИНИ
Актуальтсть. У статт розглядаеться нелшшна стацюнарна задача обтжання тша обертання в'язкою нестисливою рщиною.
Мета роботи - розробка нового методу чисельного аналiзу стацюнарних задач обтжання тш обертання в'язкою нестисливою рщиною.
Метод. Математичною моделлю розглядуваного процесу служить нелшшна крайова задача для функцп течп, отримана переходом вщ системи рiвнянь Нав'е-Стокса до одного нелшшного рiвняння четвертого порядку. Особливютю постановки задачi обтжання е те, що крайова задача розглядаеться в нескшченнш обласп i для функцп течп ставляться як крайовi умови на межi тша, що обтжаеться, так i умова на нескшченносп. За допомогою структурного методу (методу Л-функцш) будуеться структура розв'язку задачу яка точно задовольняе вс крайовi умови задачу в тому чи^ гарантуе потрiбну поведшку функцп течп на нескшченносп. Для апроксимацп невизначених компонент структури пропонуеться два пщходи. Перший пщхщ заснований на використанш методу послщовних набли-жень, який дозволяе звести розв'язування вихщно! нелшшно! задачi до розв'язування послщовносп лшшних крайових задач. На кожному крощ ггерацшного процесу щ лшшш задачi розв'язуються методом Бубнова-Гальоркша. Другий пiдхiд до апроксимацп невизначених компонент структури заснований на застосуванш нелшшного методу Гальоркша i його пропонуеться використовувати в разi розбiжностi послiдовних наближень. У цьому випадку розв'язування вихщно! нелiнiйноi задачi зводиться до розв'язування системи нелшшних алгебра!чних рiвнянь.
Результати. Обчислювальний експеримент проведений для задачi обтiкання сфери, елшсоща обертання i двох з'еднаних елшсодав для рiзних чисел Рейнольдса.
Висновки. Проведенi експерименти пiдтвердили працездатнiсть запропонованого методу чисельного аналiзу стацюнарних задач обтжання тiл обертання в'язкою нестисливою рщиною. Перспективи подальших дослщжень можуть полягати у використаннi розробле-ного методу при реалiзацii полудюкретних i проекцiйних методiв розв'язання нестацюнарних задач.
Ключовi слова: задача обтжання, стацiонарна течiя, в'язка рiдина, функщя течii, метод Л-функцш, метод послщовних наближень, метод Гальоркша.
Ламтюгова С. Н.1, Сидоров М. В.2, Ситникова Ю. В.3
'Канд. физ.-мат. наук, ассистент кафедры высшей математики Харьковского национального университета городского хозяйства имени А. Н. Бекетова, Харьков, Украина
2Канд. физ.-мат. наук, доцент, доцент кафедры прикладной математики Харьковского национального университета радиоэлектроники, Харьков, Украина
3Канд. пед. наук, доцент кафедры высшей математики Харьковского национального университета городского хозяйства имени А. Н. Бекетова, Харьков, Украина
МЕТОД ЧИСЛЕННОГО АНАЛИЗА ЗАДАЧИ ОБТЕКАНИЯ ТЕЛ ВРАЩЕНИЯ СТАЦИОНАРНЫМ ПОТОКОМ ВЯЗКОЙ ЖИДКОСТИ
Актуальность. В статье рассматривается нелинейная стационарная задача обтекания тела вращения вязкой несжимаемой жидкостью.
Цель работы - разработка нового метода численного анализа стационарных задач обтекания тел вращения вязкой несжимаемой жидкостью.
Метод. Математической моделью рассматриваемого процесса служит нелинейная краевая задача для функции тока, полученная переходом от системы уравнений Навье-Стокса к одному нелинейному уравнению четвертого порядка. Особенностью постановки
задачи обтекания является то, что краевая задача рассматривается в бесконечной области и для функции тока ставятся как краевые условия на границе обтекаемого тела, так и условие на бесконечности. С помощью структурного метода (метода ^-функций) строится структура решения задачи, которая точно удовлетворяет всем краевым условиям задачи, в том числе гарантирует нужное поведение функции тока на бесконечности. Для аппроксимации неопределенных компонент структуры предлагается два подхода. Первый подход основан на использовании метода последовательных приближений, который позволяет свести решение исходной нелинейной задачи к решению последовательности линейных краевых задач. На каждом шаге итерационного процесса эти линейные задачи решаются методом Бубнова-Галеркина. Второй подход к аппроксимации неопределенных компонент структуры основан на применении нелинейного метода Галеркина и его предлагается использовать в случае расходимости последовательных приближений. В этом случае решение исходной нелинейной задачи сводится к решению системы нелинейных алгебраических уравнений.
Результаты. Вычислительный эксперимент проведен для задачи обтекания сферы, эллипсоида вращения и двух сочлененных эллипсоидов для различных чисел Рейнольдса.
Выводы. Проведенные эксперименты подтвердили работоспособность предложенного метода численного анализа стационарных задач обтекания тел вращения вязкой несжимаемой жидкостью. Перспективы дальнейших исследований могут заключаться в использовании разработанного метода при реализации полудискретных и проекционных методов решения нестационарных задач.
Ключевые слова: задача обтекания, стационарное течение, вязкая жидкость, функция тока, метод ^-функций, метод последовательных приближений, метод Галеркина.
REFERENCES
1. Lojcjanskij L. G. Mehanika zhidkosti i gaza. Moscow, Drofa, 2003, 840 p.
2. Kutepov A. M., Poljanin A. D., Zaprjanov Z. D., Vjaz'min A. V., Kazenin D. A. Himicheskaja gidrodinamika: spravochnoe posobie. Moscow, Kvantum, 1996, 336 p.
3. Rvachev V. L. Teorija R-funkcij i nekotorye ee prilozhenija. Kiev, Nauk. dumka, 1982, 552 p.
4. Shapiro V. Semi-analytic geometry with R-functions, Acta Numérica, 2007, V. 16, pp. 239-303.
5. Kolosova S. V. Primenenie proekcionnyh metodov i metoda R-funkcij k resheniju kraevyh zadach v beskonechnyh oblastjah: diss. ... k. f.-m. n.: 01.01.07. Vychislitel'naja matematika. Har'kov, HIRJe, 1972, 85 p.
6. Kolosova S. V., Sidorov M. V. Primenenie metoda R-funkcij k raschetu ploskih techenij vjazkoj zhidkosti, Visnyk Harkivs'kogo nacional'nogo universytetu. Ser. Prykl. matem. i meh., 2003, No. 602, pp. 61-67.
7. Suvorova I. G. Komp'juternoe modelirovanie osesimmetrichnyh techenij zhidkosti v kanalah slozhnoj formy, Vestnik nacional 'nogo tehnicheskogo universiteta «HPI», 2004, No. 31, pp. 141-148.
8. Tevjashev A. D., Gibkina N. V., Sidorov M. V. Ob odnom podhode k matematicheskomu modelirovaniju ploskih stacionarnyh techenij vjazkoj neszhimaemoj zhidkosti v konechnyh odnosvjaznyh oblastjah, Radiojelektronika i informatika, 2007, No. 2 (37), pp. 50-57.
9. Maksimenko-Shejko K. V. Matematicheskoe modelirovanie teploobmena pri dvizhenii zhidkosti po kanalam s vintovym tipom simmetrii metodom R-funkcij, Dopovidi Nacional'noi' akademii'nauk Ukrai'ny, 2005, No. 9, pp. 41-46.
10. Artyukh A., Sidorov M. Mathematical modeling and numerical analysis of nonstationary plane-parallel flows of viscous
incompressible fluid by R-functions and Galerkin method, Econtechmod, 2014, Vol. 3, No 3, pp. 3-11.
11 .Lamtyugova S. N. Matematicheskoe modelirovanie linearizovannyh zadach obtekanija v sfericheskoj i cilindricheskoj sistemah koordinat, Visnyk Zaporiz'kogo nacional'nogo universytetu. Ser. Fizyko-matematychni nauky, 2012, No. 1, pp. 112-122.
12. Lamtjugova S. M., Sidorov M. V. Zastosuvannja metodu R-funkcij do rozrahunku zovnishnih povil'nih techij v'jazkoi ridini, Vidbir i obrobka informacii, 2012, No. 36 (112), pp. 56-62.
13. Lamtyugova S. N., Sidorov M. V. Numerical analysis of the external slow flows of a viscous fluid using the R-function method, Journal of Engineering Mathematics, 2015, Vol. 91, No. 1, pp. 59-79. DOI: 10.1007/s10665-014-9746-x
14. Kolosova S. V., Lamtyugova S. N., Sidorov M. V. Primenenie iteracionnyh metodov k resheniju vneshnih zadach gidrodinamiki, Radiojelektronika i informatika, 2012, No. 3, pp. 13-17.
15. Lamtyugova S. N. Primenenie iteracionnyh metodov k raschetu obtekanija tel stacionarnym potokom vjazkoj zhidkosti, Radiojelektronika i informatika, 2015, No. 2, pp. 49-56.
16. Kolosova S. V., Lamtyugova S. N., Sidorov M. V. Ob odnom metode chislennogo analiza vjazkih techenij, uslozhnennyh massoobmenom (zadacha obtekanija), Radiojelektronika i informatika, 2014, No. 1 (64), pp. 25-30.
17.Poljanin A. D. Spravochnik po linejnym uravnenijam matematicheskoj fiziki. Moscow, Fizmatlit, 2001, 576 p.
18. Van-Dajk M. Al'bom techenij zhidkosti i gaza. Moscow, Mir, 1986, 184 p.
19. Batchelor G. K An introduction to fluid dynamics, Cambridge: Cambridge University Press, 1967, 615 p.
20. Taamneh Y. CFD Simulations of Drag and Separation Flow Around Ellipsoids, Jordan Journal of Mechanical and Industrial Engineering, 2011, Vol. 5, No. 2, pp. 129-132.
