

2022 year, number 4
A.V. Dovgal, M.M. Katasonov, V.V. Kozlov, A.M. Pavlenko
Khristianovich Institute of Theoretical and Applied Mechanics SB RAS, Novosibirsk, Russia
Keywords: boundary layer, hydrodynamic stability, transition to turbulence, localized flow disturbances
Abstract >>
A review of experimental results obtained with the participation of the authors of this article, as well as in other studies on the generation and development of localized perturbations in two and threedimensional boundary layers, is presented. The accompanying wave phenomena and the contribution due to the spatialtemporal deformations of the velocity field of the initial laminar flow to its transition into turbulent state in boundary layers on the surface of aircraft, including small unmanned aerial vehicles (UAVs), are discussed.

V.I. Kornilov
Khristianovich Institute of Theoretical and Applied Mechanics SB RAS, Novosibirsk, Russia
Keywords: axisymmetric body, turbulent boundary layer, distributed blowing, perforated region, friction
Abstract >>
Results of numerical and experimental investigations of the efficiency of distributed air blowing through a perforated region of the surface of an axisymmetric body of revolution in an essentially incompressible flow with the Reynolds number Re_{L} = 4,24·10^{6} are reported. The blowing factor C_{b} is varied in the interval 00.00885. The Reynolds number Re^{**} based on the momentum thickness δ^{**} ahead of the perforated region is 5600. As the streamwise coordinate х increases, stable reduction of local friction is observed up to the distance of 600 δ^{**} from the blowing region; the maximum value of local friction reaches 56.5% directly in the region of blowing with the maximum intensity. It is found that spontaneous blowing, which occurs due to the natural difference between the barometric and static pressures in the wind tunnel test section, can also ensure reduction of the skin friction coefficient whose value in the abovenoted flow region in the main operation regime is 28.5%. Specific features of the flow around the body of revolution in the case of a rough perforated surface is analyzed through numerical simulations, and it is demonstrated that it is necessary and important to take into account this factor accompanying the blowing process.

I.I. Mazhul, Yu.P. Gounko
Khristianovich Institute of Theoretical and Applied Mechanics SB RAS, Novosibirsk, Russia
Keywords: numerical simulation, supersonic flow, longitudinal slots, air bypass, airintake start
Abstract >>
The paper presents the results of numerical simulation of starting conditions for a supersonic air intake using the bypass of part of the captured air flow through longitudinal slots. The airintake duct consists of a tapered inlet section and a constantcrosssection throat built using flat surfaces, and it contains longitudinal bypass slots stretching along the flat bottom surface. Various types of the sectional bypass of the air flow are investigated, and the change of the airintakeduct flow structure depending on the flow bypass mode is demonstrated. Numerical modeling of the threedimensional flow was performed for freestream Mach number M_{∞} = 4 on the basis of Reynoldsaveraged NavierStokes equations and the κω SST model of turbulence.

A.I. Fedorchenko^{1,2}, F. Marsik^{1}, V.I. Terekhov^{3}, V.V. Terekhov^{3}
^{1}Institute of Thermomechanics CAS, Prague, Czechia ^{2}Kutateladze Institute of Thermophysics SB RAS, Novosibirsk ^{3}Kutateladze Institute of Thermophysics SB RAS, Novosibirsk, Russia
Keywords: heat transfer, impinging liquid jet, selfsimilar solution, impulse jet
Abstract >>
The results of an analytical study of unsteady heat transfer in the vicinity of a stagnation point of an axisymmetric impinging liquid jet are presented. A selfsimilar solution is obtained, which allows the analysis of heat transfer behavior as a function of the Reynolds number, distance from the nozzle to the surface, and characteristic time. The behavior of the Nusselt number for the solution in asymptotically limiting cases is shown: the Froude number Fr >> 1 (small distances from the nozzle to the barrier) and Fr << 1 (large distances). It is noted that the effect of the Froude number can be quite significant and should be taken into account when predicting heat transfer. Significant intensification of heat transfer is shown for a pulsed impinging jet at short pulse durations.

P. Eiamsaard^{1}, K. Wongcharee^{2}, K. Kunnarak^{2}, Kumar Manoj^{3}, S. Eiamsaard^{2}, M. Pimsarn^{4}
^{1}Rajabhat Rajanagarindra University, Chachoengsao, Thailand ^{2}Mahanakorn University of Technology, Bangkok, Thailand ^{3}DIT University, Dehradun Uttarakhand, India ^{4}King Mongkut’s Institute of Technology Ladkrabang, Bangkok, Thailand
Keywords: heat exchanger, heat transfer visualization, swirling impinging jet, triple spiralcorrugated nozzle, twistedtape
Abstract >>
Thermal visualization using a thermochromic liquid crystal sheet was performed to investigate and compare heat transfer behavior of a swirling impinging jet emerging from a triple spiralcorrugated nozzle with a twisted tape ( SIJ with TT ), a SIJ issuing through a triple spiralcorrugated nozzle, and a typical/conventional impinging jet ( CIJ ) emerging from a smooth straight circular nozzle. The experimental results showed that the stronger jet recirculation close to the wall due to the swirl flow created by triple spiralcorrugated nozzle and twisted tape inserts ( SIJ with TT ) support in reducing the value of Nusselt number between the stagnation region and surroundings. For the SIJ with TT , the average Nusselt number (Nu) performed by the jets with L/d_{H} shows that the magnitude of the heat transfer coefficient was significantly enhanced with increasing twist ratios. The maximal value of the Nusselt the SIJ with TT was obtained at the lowest twist ratio ( y/W = 2.0). This was up to 5.7 % and 35.5 % higher than those of the SIJ with TT at y/W = 4.0, and the CIJ, respectively.

M.V. Salnikov^{1}, K.S. Kinzin^{1,2}, V.A. Ivaschenko^{1}, I.V. Naumov^{1,2}, R.I. Mullyadzhanov^{1,2}
^{1}Kutateladze Institute of Thermophysics SB RAS, Novosibirsk, Russia ^{2}Novosibirsk State University, Novosibirsk, Russia
Keywords: direct numerical simulation, NavierStokes equations, lattice Boltzmann equations, swirl flow, vortex breakdown
Abstract >>
The paper presents a systematic view for a swirl flow produced within a cylindrical container with fixed walls and a rotating top lid using two theoretical approaches: the lattice Boltzmann equations and NavierStokes equations. The flow modes for this confined liquid flow were tested as functions of two parameters: Reynolds number and the cylinder height to radius ratio. The lattice convergence of both solutions was analyzed. The simulation data exhibit compliance between these two variants and with the available experimental data, including the case of the development of recirculation flow at the cylinder axis (which corresponds to the event of vortex breakdown).

A.A. Kirdyashkin, A.G. Kirdyashkin, V.E. Distanov, I.N. Gladkov
Sobolev Institute of Geology and Mineralogy SB RAS, Novosibirsk,Russia
Keywords: free convective flows, subduction zone, plume conduit, melt, primary chamber, heat power
Abstract >>
A model of a thermochemical plume in the subduction zone is presented here, based on experimental modeling and theoretical analysis, as well as geological and geophysical data. A thermochemical plume originates at the boundary between the upper and lower mantle when a chemical additive is present in the crustal layer of the subducting oceanic lithospheric plate, which lowers the melting point of the crustal layer substance. The structure of thermal gravitational flows in the plume conduit melted out in the crustal layer was established due to experimental modeling. The heat power at the plume base and the heat power transferred by the plume conduit to the surrounding mantle were determined. The depth of location of a primary magma chamber, from which the thermochemical plume originates and rises to the surface of the continent, where the volcano is formed, was determined depending on the crustal layer thickness, the rate of subduction, and the angle of subducting lithospheric plate inclination. The primary magma chamber is formed in a region where the rates of subduction and melting of the crustal layer are equal in magnitude and directed oppositely. A model of the origin of a thermochemical plume on the primary chamber roof is presented, and the conditions of the plume conduit outcrop are determined. Based on theoretical modeling, geological and geophysical data on the depth of the primary chamber origin, the thermal and hydrodynamic conditions for the existence of a thermochemical plume incipient from the primary chamber and responsible for volcano formation on the surface were found.

N.S. Belousova^{1,2}, O.G. Glotov^{1,2}
^{1}Voevodsky Institute of Chemical Kinetics and Combustion SB RAS, Novosibirsk,Russia ^{2}Novosibirsk State Technical University, Novosibirsk,Russia
Keywords: combustion, titanium particles, fragmentation, aerodynamic drag coefficient, particle motion
Abstract >>
The results are reported on the study of burning titanium particles with the diameter of 200  550 µm while free fall in air. Video with the rates of 25 and 500 fps captures the particle trajectories. The mathematical procedure for averaging of particle trajectories was developed. The trajectories can be classified into classes with the particle size variation within each class about 20  30 µm. The evolutions of particle coordinate x(t) and particle velocity v(t) as functions of time t (until the event of particle fragmentation) were presented as generalized functions for each of the sizeclasses. The comparison of empirical curves x(t) and v(t) with the analytical solutions for this problem about a spherical particle motion under gravity and aerodynamic drag forces gives the effective aerodynamic drag coefficient for a burninginair titanium particle in the format C_{d} = A /Re, where Re is the Reynolds number. For the particles within the tested size range, the parameter А is independent of the particle diameter and equals about 61 ± 2 for the particle motion with the Reynolds number varying from 1 to 10. If we take the aerodynamic drag coefficient C_{d} = 61/Re and the air viscosity equal to 5,07·10^{5} Pa·s, then the motion of burning titanium particles in air can be described by the known analytical solution for a problem of spherical particle motion driven by forces of gravity and drag within the accuracy provided by the empirical equations.

V.A. Kuzmin, I.A. Zagrai, N.A. Shmakova
Vyatka State University, Kirov,Russia
Keywords: solidpropellant rocket engine, thermal radiation, heterogeneous combustion products, emission characteristics, radiation energy flux density, emissivity factor
Abstract >>
Computational experiments in the field of thermal radiation emitted by combustion products of rocket engines allow tracing the effect of determining factors on the emission characteristics, which serves to plan, predict, and interpret a physical experiment. The paper considers the influence of the condensed phase (aluminum oxide particles), the gas phase and its individual components (H_{2}O and CO_{2}) on spectral and integral radiation energy flux densities and emissivity factors in various sections of the flow of heterogeneous combustion products of a model solidpropellant rocket engine (SPRE). A comparison of the emission characteristics of incomplete combustion products and combustion products is performed for the chamber and initial section of the plume of the firststage SPRE of Trident II (D5). The calculations of spectral radiation intensity of an initial section of the plume are compared with the results of other authors.

E.Yu. Shadrin^{1}, I.S. Sadkin^{1,2}, E.P. Kopyev^{1}, I.S. Anufriev^{1,2}, V.V. Leshchevich^{3}, S.Yu. Shimchenko^{3}
^{1}Kutateladze Institute of Thermophysics SB RAS, Novosibirsk, Russia ^{2}Novosibirsk State Technical University, Novosibirsk, Russia ^{3}Luikov Heat and Mass Transfer Institute of NAS, Minsk, Belarus
Keywords: superheated steam, spraying burner, liquid fuel, steam, atomization, shadow photography method, disperse composition
Abstract >>
The characteristics of a fuel spray atomized by a steam jet were studied using the method of shadow photography. The studied method of forming a twophase flow for dispersing liquid fuel allows application of a wide range of hydrocarbons and increases the service life of the combustion equipment due to the absence of fuel spraying nozzles. Using a longfocus macroscopic lens, the dispersed composition of spent engine oil was measured at various fuel supply frequencies: 10, 25, and 40 Hz. The dispersed phase velocity was determined using the PIValgorithms; it amounted to 60 m/s for all studied regimes. It is shown that the frequency of liquid fuel supply does not affect the size and velocity of the formed fuel droplets.

B. Boudraa, R. Bessaih
University of Mentouri Brothers, Constantine 1, Algeria
Keywords: mixture model, wavy wall, laminar, forced convection
Abstract >>
In the current study, we conducted a numerical analysis of waterTiO_{2} nanofluidand entropy generation in a wavy channel under constant heat flow( q^{"}) by using a twophase mixture model. The analyses were carried under a laminar forced convection flow condition. Reynolds numbers (Re) are considered in the range of 50≤Re≤600, with the volume fraction of nanoparticles φ are considered in the range of 1 ≤ φ ≤ 5 %. The governing equations are solved by using AnsysFluent software 14.5.The validation of the outcomes has demonstrated a strong consensus between the results and the literature's data. The effect of various Re, φ, and wavy amplitudes (α) on the flow behavior and heat transfer are examined. Moreover,the distribution of the static temperature, streamlines, total entropy generation S_{g,t}, Bejan number (Be) contours have been presented and discussed. Results showed that the heat transfer rate improves when φ , Re, and α increase. The heat transfer is enhanced when using a wavy wall compared to a straight wall. The performance evaluation criterion (PEC) increases with φ and Re. Therefore, it is recommended to use large values for each of φ and Re in the wavy channel from engineering and economics perspectives. Concerning S_{g,t}, S_{g,t,h} Be, we conclude that these parameters decrease when φ and Re increase, except S_{g,v} is increased with φ and Re.

R. Lazarova^{1}, G.E. Georgiev^{1}, A.N. Cherepanov^{2}, V. Dyakova^{1}
^{1}Institute of Metal Science, Equipment and Technologies with Hydro and
Aerodynamics Centre. "Acad. A. Balevki" – Bulgarian Academy of Sciences, Sofia, Bulgaria ^{2}Khristianovich Institute of Theoretical and Applied Mechanics SB RAS, Novosibirsk, Russia
Keywords: electronbeam treatment, nanoparticles, modifier, crystal structure
Abstract >>
The paper presents experimental and theoretical study of crystal nucleation and growth in aluminum and structure development in the aluminum melt after the introduction of a cubicshaped nanoparticle modifier and using the electronbeam surface treatment method. The output of this study is the rate of solid phase nucleation as a function of TiCN nanoparticle size dispersed in molten aluminum. The numerical simulation of crystal structure growth for a sample of AlSi12Cu2NiMg alloy treated by electron beam was performed using the MAGMASOFT computer code.

A.B. Kruglov, V.I. Rachkov, I.G. Merinov, V.S. Kharitonov, L.P. Paredes
National Research Nuclear University “MEPhI”, Moscow,Russia
Keywords: pulse heating method, thermal conductivity coefficient of lead, impurities in lead, thermal resistance of contact of lead melt and steel
Abstract >>
The article presents the results of measuring the coefficient of thermal conductivity of lead in the temperature range of 3501000 °C using the pulse heating method. The methodology of processing experimental data is described. The estimates of the experimental data error are given. The difference in the content of impurities in the lead samples is shown to have an in significant effect on the thermal conductivity coefficient of the lead melt. The deviation of the experimental data on the thermal conductivity of lead from the proposed interpolating dependence does not exceed ±2 %. The obtained data are compared with the known recommended dependences for calculating the thermal conductivity coefficient of lead.

A.Sh. Agazhanov, R.N. Abdullaev, D.A. Samoshkin, Yu.M. Kozlovskii
Kutateladze Institute of Thermophysics SB RAS, Novosibirsk, Russia
Keywords: magnesiumlithium alloy, eutectic composition, thermal conductivity, thermal diffusivity, specific heat capacity, density, linear thermal expansion coefficient, relative elongation
Abstract >>
A comprehensive study of the thermophysical properties of an ultralight eutectic MgLi alloy (with a lithium content of 23 at.%), promising for use in the aerospace industry has been carried out. New reliable experimental data have been obtained for thermal conductivity, thermal diffusivity, specific heat capacity, density, linear thermal expansion coefficient and relative elongation of the alloy in the temperature range of 300  759 ÷781 K. The obtained results have been compared with wellknown literature data on thermal conductivity and thermal expansion for MgLi system alloys. Approximation equations and a table of reference data of the studied properties are presented. The relative density change during melting of eutectic has been determined.

S.D. Sleptsov^{1}, N.A. Savvinova^{2}
^{1}Kutateladze Institute of Thermophysics SB RAS Novosibirsk, Novosibirsk, Russia ^{2}NorthEastern Federal University, Yakutsk, Russia
Keywords: ice, melting, Stefan problem, anisotropic scattering, selectivity, substrate
Abstract >>
Numerical calculation of the influence of various types of opaque substrates on ice melting was carried out using mathematical modeling methods. It is shown that while maintaining a constant low temperature of the left boundary of the heatconducting substrate, ice melting does not depend on its thermophysical properties. The threshold value of thermophysical properties, below which the dependence on the substrate must be taken into account, is found.

Abstract >>
4 August, 2007 is the 75th anniversary of the well known scientist, scientific administrator and head of scientific school of Tymen Division of Khristianovich Institute of Theoretical and Applied Mechanics of the Siberian Branch of the Russian Academy of Sciences, Doctor of Physical and Mathematical Sciences, Professor Amir A. Gubaidullin.

