CC BY
12
Ovil Aviation High Technologies
Vol. 22, No. 05, 2019
ТРАНСПОРТ
05.22.01 - Транспортные и транспортно-технологические системы страны, ее регионов и городов, организация производства на транспорте;
05.22.08 - Управление процессами перевозок;
05.22.13 - Навигация и управление воздушным движением;
05.22.14 - Эксплуатация воздушного транспорта
UDC 621.1.016+532.526
DOI: 10.26467/2079-0619-2019-22-5-8-18
MATHEMATICAL MODELING AND THE STUDY OF EXCHANGE PROCESSES IN DISPERSE BOUNDARY LAYER CONTROL ACTIONS
A.A. BONDARENKO1, V.N. KOVALNOGOV2, R.V. FEDOROV2, A.V. CHUKALIN2
1Ulyanovsk Civil Aviation Institute named after Chief Marshal of Aviation B.P. Bugaev,
Ulyanovsk, Russia 2Ulyanovsk State Technical University, Ulyanovsk, Russia
The work was supported by the scholarship of the President of the Russian Federation for the project
№ SP-2069.2018.1
A significant interest of researchers is attracted to the effective management and forecasting of exchange processes in the boundary layer, which are key for the implementation of effective and reliable equipment. Modeling of exchange processes occurring in a high-speed dispersed boundary layer with external influences is a very difficult task. Mathematical modeling allows us to develop reliable devices and engines for the fields of aircraft, energy, shipbuilding with minimal costs for its creation. Despite the interest of numerous groups of researchers around the scientific projects and a large number of works, the current theory of the boundary layer is imperfect. This may be due to several circumstances: firstly, the theory of single-phase turbulent flows of continuous media is far from being completed, secondly, turbulent flows with dispersed impurities in the form of particles greatly complicate the already intricate flow pattern. Interest in dispersed flows is particularly relevant due to the fact that almost all gas-dynamic flows contain a certain concentration of particles, and their impact can provoke significant changes in the structure of the boundary layer and affect the intensity of exchange processes. The article proposes a two-fluid mathematical model describing the motion of a high-speed dispersed boundary layer on a surface with hemispherical damping cavities. The use of hemispherical damping cavities allows to reduce turbulent exchange in the boundary layer, which makes it possible to control the intensity of metabolic processes. The possibility of a significant reduction of turbulent heat transfer and friction in the dispersed boundary layer is established. The proposed method of impact on the turbulent transport in the boundary layer will improve the equipment and installations, including GTU and GTE used in various industries of our country, such as energy, aircraft, shipbuilding.
Key words: turbulent transport, hemispherical damping cavities, mathematical modeling, boundary layer, heat transfer.
INTRODUCTION
The technical progress is based on the constant equipment efficiency and performance upgrade. The actual trend in efficiency raise is the management of exchange processes in the boundary layer. The management method uses the decrease in intensity of highly turbulent pulses in the boundary layer achieved due to turbulence energy damping in hemispheric damping cavities. The implementation of the method may also result in the boundary layer lamination, the reduction in heat-exchange, the decrease of vibration processes generation, and aerodynamic noise pollution decrease. The theoretical grounds being developed will allow to show the gas-dynamic conditions within dynamically unsteady boundary layer and to optimize the processes for particular technical devices and systems.
Vol. 22, No. 05, 2019
Ovil Aviation High Technologies
THE THEORETICAL GROUNDS AND RESEARCH METODS OF DISPERSED FLOWS
The research of dispersed flows with phase transitions is carried out at the research centers of the USA, the UK, Germany, Denmark, China, Japan, Israel and Russia. The theoretical development is intense for the aspects of heat and mass exchange and fluid and gas mechanics connected with the dispersed two-phase flows.
The developer of mathematical models and methods faces certain challenges. Firstly, numerous physical processes taking place in the dispersed flows have to be taken into consideration. Secondly, physical processes in conditions of high temperatures are not always indisputable, so the modelling methods developed may insufficiently adequate results. Thirdly, the model cannot be comprehensive, it has to have certain application boundaries and limits. Thus, the research supposes thorough analyses and accentuation on the most critical physical aspects for the problem solution.
At present, the active research on the dispersed flows has led to a significant knowledgebase. The paper [2] provides a review of both theoretic and experimental research on the turbulent dispersed flows around bodies of different shapes: sphere, cylinder, plates and cascaded airfoils. It has been proved that the presence of particles in the flow may cause significant erosive wear, gas dynamic spatter, icing, light emission increased heat exchange and overheating. Such phenomena, for instance, appear when aircraft are travelling through dust-filled atmosphere, also in case of powerplants in operation [5, 6, 7]. The basic phenomena in the context of dispersed flows are spotted in the boundary layer. They lead to the increased intensity of metabolic processes, vortex shedding and turbulent wake formation [8] In order to correctly assess the physical processes in the dispersed boundary layer it is necessary to consider the velocity distribution, the temperatures (non-isothermal condition), and the particle concentration depending on the flow field around the surface and the surface type. The papers [2, 9] classify the possible disperse flow field regimes using the Stokes number StKf for the mean motion, which characterizes the velocity relaxation for gasses and particles in the vicinity of the boundary/wall [2]. The papers [10, 11] observe the influence of the particles on the boundary-layer drag and heat-exchange noticing the quality changes in metabolic processes depending on the particle concentration, size and flow angles.
According to the latest tendencies, the two-phase flow calculation contains the modeling of impulse, the mass and heat transfer for each of the phases, and also the modeling of phase interaction within the flow. The current two-phase flow mathematic models are divided into two classes: Euler (two-fluid) and Euler-Lagrange.
The theory and application of two-fluid model calculation is shown in research papers [1, 2, 12].
Two-fluid models have the advantage of using similar equations to describe the homogeneous and non-homogeneous flows thus enabling to solve the system of equations using one numerical approach and, as follows, to increase the timing efficiency of the modelling adequacy processing and check. The criticisms of such models are: the solution error resultant from the disperse phase approximation: the missing data of the motion of separate particles in the non-homogeneous flow; the cumbersome boundary conditions for non-homogeneous flows on the surfaces causing the flow restrictions. The advantage of Euler-Lagrange model is the detailed approach. The model takes into consideration the complete information about the environment, including the particle motion and their interaction in the vector flow with the known temperature field. The principal difficulties are connected with the specific turbulent nature of the motion: the interaction of the dispersed cluster particles with the external phase turbulent vortices and with one another as a result of the particle contact with the boundary/wall and setting on the wall; the influence of turbulent fluctuations on the phase transition speed; the back reaction of the particles on the turbulence [1]. Having analyzed both
Ovil Aviation High Technologies
Vol. 22, No. 05, 2019
approaches, we should acknowledge that the criticisms of the first class of models are the advantages of the second class of models, and vice versa.
The paper by L.I. Zaichik and V.M. Alipchenko studies the statistic theory and continual method of the modelling of hydrodynamics and heat-exchange in two-phase flows. The research thoroughly investigates the basic theoretical problems caused by the particle motion in the turbulent flows, the phenomenon of particle accumulation in the near-wall layer. The authors claim that statistic models can be used for effective accounting of the particle interaction with the small-scale vortices to determine the sub-grid turbulence within the large-scale vortex method.
Considering the challenges mentioned above, it is impossible to come to direct numerical method of solution for the system of differential equations. Making equations for the instant parameters is problematic for modeling of the dispersed flow turbulence. At the same time it is known that outer effects on the dispersed flow typical of gas turbine engines lead to a number of specific features of turbulent transfer, which as a matter of principle cannot be considered neither using a classical model of Prandtl "mixing length" model, nor the widely differential two-parameter "k-e" used turbulence model. To cope with the difficulties mentioned above it is necessary to develop the model which will adequately describe the motion and the heat-exchange of the dispersed flow in the boundary layer considering the outer effects.
To implement the model, it is essential to develop as an analytical kernel the two-fluid model of the boundary layer. The original mathematical model takes into consideration the inner sources of the heat and the amount of motion, and also the turbulent transfer in the boundary layer and the outer factors in action. Moreover, the immediate effect of the particles on the carrying medium will be considered by the boundary conditions and corresponding members, which are conditional to the particle impact.
The intermediate impact will be considered using the eddy diffusivities of momentum and heat.
The two-fluid model [4, 13], describing the high-speed boundary layer with outer impact, will be shown as differential equations:
The differential equation of motion:
(1)
The differential equation of energy:
/ди
дТ
U- + Vox
(2)
The differential equation of continuity:
дт дх ду
(3)
Pressure-volume relation
Vol. 22, No. 05, 2019
Civil Aviation High Technologies
where T - temperature; t - time; u, v - longitudinal velocity component, cross-stream velocity component; x, y -longitudinal and transverse coordinate; p - density; cp - specific heat capacity at constant pressure; X - code of heat conductivity; XT - eddy diffusivity of heat; p - coefficient of dynamic viscosity ; pT - eddy diffusivity of momentum; qv - internal heat source intensity value; sv -eddy momentum flux intensity value; p - flow pressure; R - individual gas constant.
To solve the system of equations it is necessary to find the values of XT, qv, sv.
The parameters which characterize the presence of the disperse phase in the boundary layer are considered through the thermal and aerodynamic particle influence [4]:
qv = a(Tm -T); s =%m -u), (5)
where a and b are the complexes characterizing the particle parameters; Tsm - normal temperature distribution; usm - normal velocity distribution
6as Ps
a = PT (6)
Pe ds
b = 0,75psPcfs|us -u Pe ds
where as,cfs - unit surface conductance and dispersed medium condensed particle resistance coefficient; ps - condensate phase density (the mass of condensed particles per the unit of the medium volume); pB - density of particle matter; us, Ts - particle velocity and temperature; ds - particle diameter.
The complexes a and b are averaged within the boundary layer, as they encounter insignificant changes (at different iterations at the moment of clarification of the medium flow field), relative to qv and sv.
The profiles of temperature Tsm and velocity usm and the particles of the representative fraction in the sections of the boundary layer are determined using the equations [4]
u = u
sm smn
Фu +0 -Фи ) —
u
• T =T
' sm smxi
T
ФT +(1-Фг ))
(9)
where - coefficients varied from 0 to1 and characterizing the particle diameter and cross travel
speed. The more is the cross travel speed and its diameter, and the less is the boundary layer gauge, the bigger is the coefficient, and, accordingly, vice versa. The research conditions were = = = 1.
To find the value of the coefficients a, b, , the particle trajectory and parameters are calculated using the methods listed in [4].The methods suggested allow to estimate and calculate the boundary layer for the case of the high-speed dispersed flow around the surface. The most common disperse phase flows are the CO2 (exhaust gas) flows. The sizes of small solid particles in the fluidized state, formed by the fuel combustion range from 10-8 to 10-6 depending on the fuel type, grade and the combustion process. The exhaust gases are seldom emitted into the
Научный Вестник МГТУ ГА_Том 22, № 05, 2019
Civil Aviation High Technologies Vol. 22, No. 05, 2019
atmosphere, as they contain an enormous amount of heat energy. For different types of engines the exhaust gases are used for heating and inner and outer joints ice accretion prevention (compressor guide vanes, casing), in energy efficient gas turbine engines the exhaust gases are used for heat energy recuperation by passing the exhaust gas heat to the air which is later used for fuel combustion. Papers [15, 16] introduce the ways of heat regeneration and recuperation, and also the schemes of heat regeneration in aviation gas turbine engines. The analyses show that a significant part of the exhaust gas energy (kinetic and heat) is wasted before it reaches the areas of potential use. Thus, there is a task of methods development and research for maintaining the exhaust gases energy so that to increase the efficiency of their use.
For the conditions in question, the methods stated in [17] determine the thermophysical properties of CO2
f = (Г/Г0)0-717; Sl = (T/T0)0939 ; f = (7/70)0-236 (10)
Mo сро л0
One of the effective and perspective ways of the boundary layer influence is the use of semispherical damping cavities. (figure 1) The presence of damping cavities on the surface allows to reduce the turbulent-flow resistance and the heat-exchange of the solid body in flow. The outer impacts on the boundary layer, including the semispherical damping cavities impact may be expressed in terms of eddy diffusivity of momentum and eddy diffusivity of heat determined by means of the modified Prandtl "mixing length" model:
AT = ^Tcp/Pr; ^T = pl2du/dy. (11)
where Pr - Prantdl number; l - the mixing length, calculated as:
I = ^y[1 - exp(-pv*y/26^]. (12)
where x - coefficient showing the intensity of the turbulent transfer in the boundary layer; v* -dynamic speed.
The semispherical damping cavities will have the immediate influence on the coefficient x, as they affect the intensity of the boundary layer transfer. The decreased intensity of the turbulent transfer occurs due to the turbulent vortices suppression in the semispherical damping cavities (figure 1).
To model the impact of semispherical damping cavities on the boundary layer in a broad band of gas-dynamic conditions, it is suggested that the coefficient expressed the intensity of the turbulent transfer momentum.
« ^1-8.4-105A*vf2exp(1-n)
«о 1 + 21 4 dum/dx ' u0(du/dy)
(13)
where A*v - empirical coefficient, which characterizes the impact of semispherical damping cavities depending on the cavity volume; f - perforation relative area; n - the number of perforation holes; u0 -velocity scanning value; um - velocity at the boundary layer ceiling.
Vol. 22, No. 05, 2019
Ovil Aviation High Technologies
Fig. 1. Diagram of the damping cavity (a) and the damping surface (b) with two perforations, hemispherical form: 1 - perforations; 2 - wall; 3 - hemispherical damping cavity
The coefficient was determined using the approximation of the experiment data. The method of cubic/isometric approximation was used to get the coefficient with minimum error [18]. The empiric data on the reduced turbulent friction on the surface under the impact of semispherical damping cavities were used as approximating values [19].
THE RESULTS OF THE RESEARCH OF DISPERSED BOUNDARY LAYER ON THE SURFACE WITH SEMISPHERICAL DAMPING CAVITIES
To study the exchange processes and conditions of semispherical damping cavities impact the modelling has been carried out using the two-fluid modified mathematical model of high-speed dispersed boundary layer.
The experiment conditions: the flat surface in the dispersed flow without the inertial particle fallout = 1; = 1; tys = 0. The surface length given was 0.6 m. The semispherical damping cavities with the volume of 0.575 cm3 , their contact with the flow was via two perforated holes 0=0.8 mm, the relative perforation area on the surface f = 0.0012; u = 40 m/s; T = 640 K; R = 189 J/° kgK; the coefficient of particle aerodynamic impact G = 0.510-6; G = 0.5-10-8 (figure 2).
Fig. 2. Study of the coefficient of friction resistance in the boundary layer of the dispersed flow in the flow of smooth (a) and damping surface (b): ■ - aerodynamic coefficient of action of particles G = 0.5 10-6; ▲ - aerodynamic coefficient of action of particles G = 0.5 10-8; ♦ - pure gas flowing smooth surface
Ovil Aviation High Technologies
Vol. 22, No. 05, 2019
The presence of dispersed phase in the flow significantly increases its resistance, for the conditions in question, the resistance turbulent friction coefficient has increased by 2.7% at G = 0.510-8, and at particle impact G = 0.510-6 it increased by 9.8%. The impact of semispherical damping cavities on the boundary layer allows to reduce the turbulent friction significantly. With the impact of cavities the turbulent friction resistance coefficient is decreased by 3.1% at G = 0.510-6, and in case of the particle impact G = 0.5 10-8 the decrease reaches 17%.
Also, the research on heat loss of the turbulent dispersed flow for the smooth and damping surfaces has been carried out (figure 3).
Fig. 3. The study of heat transfer in the boundary layer of a dispersed flow on a smooth (a) and damping (b) surface: ■ - aerodynamic coefficient of the impact of particles G = 0.5 10-8, with a flow around the damping surface; ▲ - aerodynamic coefficient of the impact of particles G = 0.5 10-8, when flowing around a smooth surface;
♦ - clean gas, when flowing smooth surface
Considering the results, it is worth mentioning that the use of cavities allows to reduce the heat loss on the surface, thus reducing the heat impact on the surface and the flow heat loss. The decrease of heat exchange in the dispersed boundary layer with the aerodynamic particle impact coefficient G = 0.5 10-8 reaches 24,3%, but in case of the particle impact coefficient increase to G = 0.5 10-6 the efficiency of the impact is significantly reduced and reaches 13.6%. Such decrease results from the increase in the boundary layer turbulence due to the aerodynamic impact of the particles.
The damping cavities allow to significantly reduce the resistance and the heat exchange on the surface, which makes their use appropriate in order to increase the device efficiency.
THE RESULTS ANALYSES AND CONCLUSION
The use of semispherical damping cavities leads to significant decrease in the turbulent transfer in the boundary layer, which allows to reduce the energy of turbulent vortices and to reduce the coefficient of turbulent friction Cf/Cfo by 17%. The decrease of heat exchange using cavities for loosely dusted flows reached Nu/Nu0 24.3%, and for dusty flows the given value reached Nu/Nu0 13.6%. The example of the effective application of semispherical damping cavities is the Patent of the Russian Federation №170277 "the combustion chamber of the gas turbine engine with laminaring panels" developed to enhance the aviation gas turbine engine efficiency. The outer case of the combustion chamber implies the laminaring panels with semispherical damping cavities which facilitate the turbulent transfer in the boundary layer. The mathematical model suggested allows to forecast the effectiveness of semispherical damping cavities application on the surfaces subjected to high-speed dispersed flows.
Том 22, № 05, 2019_Научный Вестник МГТУ ГА
Vol. 22, No. 05, 2019 Civil Aviation High Technologies
REFERENCES
1. Varaksin, A.Yu. (2003). Turbulentnyye techeniya gaza s tverdymi chastitsami [The Turbulent gas flows with solid particles]. Moscow: FIZMATLIT, 192 p. (in Russian)
2. Varaksin, A.Yu. (2018). Gas-Solid Flows Past Bodies. High Temperature, vol. 56, no. 2, pp. 275-295. DOI: 10.7868/S0040364418020175
3. Pakhomov, M.A. and Terekhov, V.I. (2016). The effect of drop evaporation on gas turbulence and heat transfer for two-phase flow behind sudden pipe expansion. High Temperature, vol. 54, no. 3, pp. 330-337. DOI: 10.7868/S0040364416030157
4. Kovalnogov, N.N. (1996). Pogranichnyy sloy v potokakh s intensivnymi vozdeystviyami [The boundary layer in streams with intense effects]. Ulyanovsk, 240 p. (in Russian)
5. Mikhatulin, D.S., Polezhayev, Yu.V. and Reviznikov, D.L. (2007). Teploobmen i razrusheniye tel v sverkhzvukovom geterogennom potoke [Heat transfer and destruction of bodies in supersonic heterogeneous flow]. Moscow: Yanus-K, 392 p. (in Russian)
6. Varaksin, A.Yu. (2017). Teoriya i proyektirovaniye gazoturbinnykh i kombinirovannykh ustanovok [Theory and design of gas turbines and combined installation]. Moscow: Izd-vo MGTU im. N.E. Baumana, 680 p. (in Russian)
7. Perelman, R.G. and Pryakhin, V.V. (1986). Eroziya elementov parovykh turbin [Erosion of steam turbine elements]. Moscow: Energoatomizdat, 181 p. (in Russian)
8. Pakhomov, M.A. and Terekhov, V.I. (2016). Structure of the nonisothermal swirling gas-droplet flow behind an abrupt tube expansion. Fluid Dynamics, vol. 51, no. 1, pp. 70-80.
9. Tsirkunov, Yu.M. (2001). Gas-particle flows around bodies - key problems, modeling and numerical analysis. Proc. Fourth International Conference on Multiphase Flow, in E. Michaelides (Ed.), May 27 - June 1, New Orleans, LA, USA. - CD ROM Proc. ICMF'2001, paper no. 607, pp. 1-31.
10. Lashkov, V.A. (2012). Koeffitsiyenty vosstanovleniya skorosti pri udare tverdykh chastits gazovzvesi o poverkhnost tela: dis. ... dokt. fizmat. nauk [Ratios of speed recovery upon impact of solid particles of a gas suspension on the body surface: D. Phys-Math. Sc. Thesis]. St. Petersburg: SPbGU, 379 p. (in Russian)
11. Polezhayev, Yu.V., Repin, I.V. and Mikhatulin, D.S. (1992). Heat-transfer in a heterogeneous supersonic flow. High Temperature, vol. 30, no. 6, pp. 949-954.
12. Terekhov, V.I. and Pakhomov, M.A. (2011). Effect of particles on the flow structure and dispersion of solid impurities in a two-phase axisymmetric jet. Technical Physics. The Russian Journal of Applied Physics, vol. 56, no. 10, pp. 1406-1414.
13. Zaychik, L.I. and Pershukov, V.A. (1996). Problemy modelirovaniya gazodispersnykh turbu-lentnykh techeniy s goreniyem ili fazovymi perekhodami (obzor) [Problems of modeling gas-dispersed turbulent flows with combustion or phase transitions (review)]. Izv. RAN. MZHG, no. 5, pp. 3-19. (in Russian)
14. Zaychik, L.I., Alipchenkov, V.M. (2007). Statisticheskiye modeli dvizheniya chastits v turbulentnoy zhidkosti [Statistical models of particle motion in a turbulent fluid]. Moscow: FIZMATLIT, 312 p. (in Russian)
15. Getman, V.V., Lezhneva, N.V. (2014). Metody utilizatsii teploty ukhodyashchikh gazov ot energeticheskikh ustanovok [Methods of utilization of the heat of flue gases from power plants]. Vestnik kazanskogo tekhnologicheskogo universiteta [Herald of Kazan Technological University], pp. 104-107. (in Russian)
16. Mamayev, B.I. (1991). Puti povysheniya k.p.d. aviatsionnoy gazovoy turbiny [Ways to improve the efficiency of aviation gas turbine]. Gazoturbinnyye i kombinirovannyye ustanovki: tezisy dokl. Vsesoyuznoy konferentsii перевод названия конференции. Moscow: MGTU, pp. 43-44. (in Russian)
Civil Aviation High Technologies
Vol. 22, No. 05, 2019
17. Kovalnogov, N.N. (2010). Nature of gas influence on gas-dynamic temperature stratification efficiency. Vestnik Ulyanovskogo gosudarstvennogo tekhnicheskogo universiteta, no. 2, pp. 67-69. (in Russian)
18. Gortyshov, Yu.F., Dresvyannikov, F.N. and Idiatullin, N.S. (1985). Teoriya i tekhnika teplofizicheskogo eksperimenta [Theory and technique of thermophysical experiment]. Moscow: Energoatomizdat, 360 p. (in Russian)
19. Kovalnogov, V.N., Fedorov, R.V., Chukalin, A.V., Bondarenko, A.A. and Miroshin,
A.N. (2018). Eksperimentalnoye issledovaniye pogranichnogo sloya na perforirovannoy poverkhnosti s dempfiruyushchimi polostyami [Experimental study of the boundary layer on a perforated surface with damping cavities]. Trudy sedmoy Rossiyskoy natsionalnoy konferentsiipo teploobmenu. Moscow: переводMEI, vol. 2, pp. 373-376. (in Russian)
INFORMATION ABOUT THE AUTHORS
Alexander A. Bondarenko, Candidate of Technical Sciences, Dean of the Flight Operation and Air Traffic Control Faculty, Ulyanovsk Institute of Civil Aviation, baa-69@yandex.ru.
Vladislav N. Kovalnogov, Doctor of Technical Sciences, Professor, Head of the Thermal and Fuel Energy Chair of Ulyanovsk State Technical University, kvn@ulstu.ru.
Ruslan V. Fedorov, Candidate of Technical Sciences, Associate Professor of the Thermal and Fuel Energy Chair of Ulyanovsk State Technical University, r.fedorov@ulstu.ru.
Andrey V. Chukalin, Research Associate, Thermal and Fuel Energy Chair of Ulyanovsk State Technical University, chukalin.andrej@mail.ru.
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И ИССЛЕДОВАНИЕ ОБМЕННЫХ ПРОЦЕССОВ В ДИСПЕРСНОМ ПОГРАНИЧНОМ СЛОЕ С УПРАВЛЯЮЩИМИ ВОЗДЕЙСТВИЯМИ
1 2 2 2 А.А. Бондаренко , В.Н. Ковальногов , Р.В. Федоров , А.В. Чукалин
1 Ульяновский институт гражданской авиации имени Главного маршала авиации
Б.П. Бугаева», г. Ульяновск, Россия 2Ульяновский государственный технический университет,
г. Ульяновск, Россия
Работа выполнена при финансовой поддержке Стипендии Президента РФ
по проекту № СП-2069.2018.1
К эффективному управлению и прогнозированию обменных процессов в пограничном слое, являющихся ключевыми для реализации эффективного и надежного оборудования, привлечен значительный интерес исследователей. Моделирование обменных процессов, протекающих в высокоскоростном дисперсном пограничном слое с внешними воздействиями, является весьма сложной задачей. Математическое моделирование позволяет разрабатывать надежные устройства и двигатели для областей авиастроения, энергетики, судостроения с минимальными издержками на их создание. Несмотря на интерес многочисленных групп исследователей по всему миру и множество работ, существующая теория пограничного слоя несовершенна. Это может быть связано с несколькими обстоятельствами: во-первых, разработанная теория однофазных турбулентных течений несовершенна и на сегодняшний день содержит множество эмпирических зависимостей; во-вторых, турбулентные потоки с дисперсными примесями в виде частиц сильно осложняют и без того замысловатую картину течения. Интерес к дисперсным потокам особенно актуален вследствие того, что практически все газодинамические течения содержат некоторую концентрацию частиц, а их воздействие может спровоцировать значительные изменения структуры пограничного слоя и повлиять на интенсивность обменных процессов. В статье предложена двухжидкостная математическая модель, описывающая движение высокоскоростного дисперсного пограничного слоя на поверхности с полусферическими демпфирующими полостями. Применение полусферических
Vol. 22, No. 05, 2019
Civil Aviation High Technologies
демпфирующих полостей позволяет снижать турбулентный обмен в пограничном слое, что дает возможность управления интенсивностью обменных процессов. Установлена возможность существенного снижения турбулентного теплообмена и трения в дисперсном пограничном слое. Предложенный способ воздействия на турбулентный перенос в пограничном слое позволит усовершенствовать оборудование и установки, в том числе газотурбинных установок и газотурбинных двигателей, применяемые в разных областях промышленности нашей страны, таких как авиастроение, энергетика, судостроение.
Ключевые слова: турбулентный перенос, полусферические демпфирующие полости, математическое моделирование, пограничный слой, теплообмен.
СПИСОК ЛИТЕРАТУРЫ
1. Вараксин А.Ю. Турбулентные течения газа с твердыми частицами. М.: ФИЗМАТЛИТ, 2003. 192 с.
2. Вараксин А.Ю. Обтекание тел дисперсных газовыми потоками // Теплофизика высоких температур. 2018. Т. 56, № 2. С. 282-305. DOI: 10.7868/S0040364418020175
3. Пахомов М.А., Терехов В.И. Влияние испарения капель на турбулентность газа и теплообмен при течении двухфазного потока за внезапным расширением трубы // Теплофизика высоких температур. 2016. Т. 54, № 3. С. 352-359. DOI: 10.7868/S0040364416030157
4. Ковальногов Н.Н. Пограничный слой в потоках с интенсивными воздействиями. Ульяновск: УлГТУ, 1996. 245 с.
5. Михатулин Д.С., Полежаев Ю.В., Ревизников Д.Л. Теплообмен и разрушение тел в сверхзвуковом гетерогенном потоке. М.: ЯНУС-К, 2007. 392 с.
6. Вараксин А.Ю., Арбеков А.Н., Иванов В.Л. Теория и проектирование газотурбинных и комбинированных установок / Под ред. А.Ю. Вараксина. М.: Изд-во МГТУ им. Н.Э. Баумана, 2017. 680 с.
7. Перельман Р.Г., Пряхин В.В. Эрозия элементов паровых турбин. М.: Энергоатом-издат, 1986. 184 с.
8. Пахомов М.А., Терехов В.И. Структура течения в неизотермическом закрученном газокапельном потоке за внезапным расширением трубы // Изв. РАН. МЖГ. 2016. № 1. С. 69-78.
9. Tsirkunov Yu.M. Gas-Particle Flows around Bodies - Key Problems, Modeling and Numerical Analysis // Proc. Fourth Int. Conf. on Multiphase Flow (ICMF'01). New Orleans. USA. 2001. Paper no. 607. Pp. 1-31.
10. Лашков В.А. Коэффициенты восстановления скорости при ударе твердых частиц газовзвеси о поверхность тела: дисс. ... докт. физ.-мат. наук. СПб.: СПбГУ, 2012. 379 с.
11. Полежаев Ю.В., Репин И.В., Михатулин Д.С. Теплообмен в сверхзвуковом гетерогенном потоке // Теплофизика высоких температур. 1992. Т. 30, № 6. С. 1147-1153.
12. Терехов В.И., Пахомов М.А. Влияние частиц на структуру течения и дисперсию твердой примеси в двухфазной осесимметричной струе // Журнал техническая физика. 2011. Т. 81, № 10. С. 27-35.
13. Зайчик Л.И., Першуков В.А. Проблемы моделирования газодисперсных турбулентных течений с горением или фазовыми переходами (обзор) // Изв. РАН. МЖГ. 1996. № 5. С. 3-19.
14. Зайчик Л.И., Алипченков В.М. Статистические модели движения частиц в турбулентной жидкости. М.: ФИЗМАТЛИТ, 2007. 312 с.
15. Гетман В.В., Лежнева Н.В. Методы утилизации теплоты уходящих газов от энергетических установок // Вестник Казанского технологического университета. 2013. Т. 16, № 12. С. 104-107.
СМ! Aviation High Technologies
Vol. 22, No. 05, 2019
16. Мамаев Б.И. Пути повышения к.п.д. авиационной газовой турбины // Газотурбинные и комбинированные установки: тезисы докл. Всесоюзной конференции. М.: МГТУ, 1991. С. 43-44.
17. Ковальногов Н.Н. Влияние природы газа на эффективность газодинамической температурной стратификации // Вестник Ульяновского государственного технического университета. 2010. № 2 (50). С. 67-69.
18. Гортышов Ю.Ф., Дресвянников Н.С. Теория и техника теплофизического эксперимента / Под ред. В.К. Щукина. М.: Энергоатомиздат, 1985. 360 с.
19. Ковальногов В.Н., Федоров Р.В., Чукалин А.В. Экспериментальное исследование пограничного слоя на перфорированной поверхности с демпфирующими полостями / В.Н. Ковальногов, Р.В. Федоров, А.В. Чукалин, А.А. Бондаренко, А.Н. Мирошин // Труды седьмой Российской национальной конференции по теплообмену. М.: Издательский дом МЭИ, 2018. Т. 2. С. 373-376.
СВЕДЕНИЯ ОБ АВТОРАХ
Бондаренко Александр Аркадьевич, кандидат технических наук, декан факультета летной эксплуатации и управления воздушным движением Ульяновского института гражданской авиации имени Главного маршала авиации Б.П. Бугаева, baa-69@yandex.ru.
Ковальногов Владислав Николаевич, доктор технических наук, профессор, заведующий кафедрой тепловой и топливной энергетики Ульяновского государственного технического университета, kvn@ulstu.ru.
Федров Руслан Владимирович, кандидат технических наук, доцент кафедры тепловой и топливной энергетики Ульяновского государственного технического университета, r.fedorov@ulstu.ru.
Чукалин Андрей Валентинович, научный сотрудник кафедры тепловой и топливной энергетики Ульяновского государственного технического университета,
сhukalin.andrej@mail.ru.
Поступила в редакцию 24.05.2019 Received 24.05.2019
Принята в печать 24.09.2019 Accepted for publication 24.09.2019
