Thermophysics and Aeromechanics

Gavrilov, A. A.;Rudyak, V. Ya.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040016

Abstract Fully developed turbulent pipe flows of power-law fluids are studied by means of direct numerical simulation. Two series of calculations at generalised Reynolds numbers of approximately 10000 and 20000 were carried out. Five different power law indexes n from 0.4 to 1 were considered. The distributions of components of Reynolds stress tensor, averaged viscosity, viscosity fluctuations, and measures of turbulent anisotropy are presented. The friction coefficient predicted by the simulations is in a good agreement with the correlation obtained from experiment. Flows of power-law fluids exhibit stronger anisotropy of the Reynolds stress tensor compared with the flow of Newtonian fluid. The turbulence anisotropy becomes more significant with the decreasing flow index n. An increase in apparent viscosity away from the wall leads to the damping of the wall-normal velocity pulsations. The suppression of the turbulent energy redistribution between the Reynolds stress tensor components observed in the simulations leads to a strong domination of the axial velocity pulsations. The damping of wall-normal velocity pulsations leads to a reduction of the fluctuating transport of momentum from the core toward the wall, which explains the effect of drag reduction.

Kulik, V. M.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040028

Abstract The algorithm for calculation of dynamic compliance of multilayer coatings was developed. The compliance modulus and phase lag of coating surface motion vs. the current pressure depend on viscoelastic properties of materials, ratio of wavelength to layer thickness λ/H, and ratio of wave velocity to propagation velocity of shear vibrations in the base layer V / C t,2 0 Dynamic compliance of the two-layer coating consisting of a thick base layer and thin durable outer layer was calculated. The elasticity modulus of the outer layer ranged up to eight values of elastic modulus of the inner layer; the density of the outer layer either remained equal to the density of the inner layer or increased proportionally to the elastic modulus. Depending on V / C t,2 0 two scenarios of compliant coating interaction with the turbulent flow were distinguished: resonant and broadband ones. It is shown that the vibration properties of two-layer coatings can be significantly better than the properties of the monolayer coatings. This makes it possible either to increase the coating strength or to work efficiently at lower velocities.

Eghlima, Z.;Mansour, K.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S086986431604003X

Abstract This paper describes a numerical solution of the bow shock shape ahead of some blunt and sharp axisymmetric noses containing sphere, blunt cone, and sharp cone at steady transonic flow in the Mach number range of 1.01 to 1.2. For validating the results, one sphere and three blunt cones are modeled, and their shock standoff distance is compared with other experimental and numerical studies. The flow over other noses with similar geometric parameters is then solved and compared with each other. In this study, the Reynolds-averaged Navier—Stokes equations are solved using the Spalart—Allmaras turbulence model. The purpose of this study is to determine the shock standoff distance for some blunt and sharp noses at low supersonic free flight speed. The shock standoff distance is determined from the Mach number curve on the symmetry line. The present numerical simulations reach down to M8=1.01 a range where it is almost very difficult to set in experimental studies. The shock wave locations were found to agree well with previous numerical and experimental studies. Our results are closer to the experimental results compared to other numerical studies. In addition, the results for shock standoff distances over paraboloids in these speed ranges have not been previously published as far as we know.

Aleksandrov, V. Yu.;Arefyev, K. Yu.;Ilchenko, M. A.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040041

Abstract Intensification of mixing between the gaseous working body ejected through a jet nozzle with ambient medium is an important scientific and technical problem. Effective mixing can increase the total efficiency of power and propulsion apparatuses. The promising approach, although poorly studied, is generation of acoustic self-oscillation inside the jet nozzle: this impact might enhance the decay of a supersonic jet and improve the mixing parameters. The paper presents peculiar properties of acoustic self-excitation in jet nozzle. The paper presents results of experimental study performed for a model injector with a set of plates placed into the flow channel, enabling the excitation of acoustic self-oscillations. The study reveals the regularity of under-expanded supersonic jet decay in submerged space for different flow modes. Experimental data support the efficiency of using the jet nozzle with acoustic self-oscillation in application to the systems of gas fuel supply. Experimental results can be used for designing new power apparatuses for aviation and space industry and for process plants.

Changcharoen, W.;Samruaisin, P.;Eiamsa-ard, P.;Eiamsa-ard, S.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040053

Abstract The influence of co/counter dual-twisted tapes (CoT/CT) on heat transfer rate in a circular tube has been investigated experimentally. In the experiment, the dual-twisted tapes are placed at the entry of the test tube in two arrangements: (1) each of dual twisted tape was twisted in the same direction that can produce co-swirl flow at the entry and (2) each of dual twisted tape was twisted in the opposite direction that can produce counter-swirl flow. Dual tapes were twisted in three different twist ratios (y/w = 3, 4, and 5) for generating different swirl intensities at the entry of the test section while the single twisted tape (ST) was also the test for comparison. The aim at using the dual twisted tapes is to create co/counter-rotating swirl flows having a significant influence on the flow turbulence intensity at the entry section leading to higher heat transfer enhancement. Average Nusselt numbers of CoT/CT are determined and also compared with those obtained from other similar cases, i.e., ST. The experimental results on the heat transfer rates indicated that the tubes with the dual twisted tapes (CoT/CT) are higher than those with the single tape at the entry section (x/D = 0 to 10). The heat transfer rates at longer distance became lower due to high interaction of each swirl. In addition, the mean Nusselt number and friction factor for the swirl generator created by the CT is nearly similar to CoT results.

Shagapov, V. Sh.;Yumagulova, Yu. A.;Shepelkevich, O. A.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040065

Abstract The solution was obtained for a problem of gas hydrate growth in water with dissolved gas. The rate of hydrate formation depends on gas diffusion to the contact with gas hydrate. Three versions of problem configuration were considered: planar, radial, and spherical symmetry. For these cases, the values of the self-similarity coordinate were obtained: this parameter controls the growth of gas hydrate in water with gas dissolved and the level of top temperature in the hydrate zone. Analysis was performed for temperature fields related to heat release during hydrate formation.

Pavlov, P. A.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040077

Abstract Methods for calculating the geometric characteristics of the vapor phase in explosive wall boiling-up processes on a metal wall are analyzed. A monotonic growth of superheat in the liquid above the equilibrium evaporation temperature is specified. We show that the choice of the model for bubble interaction has a profound influence on the geometric characteristics which define the value of the heat flux. Computer simulation was employed to obtain the dependence of dry area on time in two interaction models. We have found that, for a model with instantaneous bubble coalescence, the dry area can be evaluated by the Kolmogorov formula using a correction factor for the most probable triple interaction. An approximation of the distribution length of wetting line over the lifetime of wetting-line segments is obtained. The possibility of using the obtained data for calculation of rapid condensation is analyzed.

Bondareva, N. S.;Sheremet, M. A.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040089

Abstract Numerical simulation of transient melting regimes inside an enclosure in the presence of a local heat source has been carried out. Mathematical model formulated in terms of dimensionless variables such as stream function, vorticity, and temperature has been numerically solved by finite difference method. Effects of the Rayleigh number 4·105 ≤ Ra ≤ 5·107, Stefan number 2.21 ≤ Ste ≤ 5.53, and dimensionless time on velocity and temperature fields as well as on the local Nusselt number along the heat source surface have been analyzed in detail. The transient effects of the considered process at high values of the Rayleigh number have been identified.

Jabbari, F.;Saedodin, S.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040090

Abstract In this study, the rate of conduction-radiation heat transfer between two thick concentric spheres is analytically investigated because of the ever-increasing importance of radiation shield applications. The heat transfer rate, the percentage of reduction in heat transfer, temperature and emissivity of surfaces are calculated when one and two thick radiation shields are placed between two thick spheres. The calculations show that the use of a radiation shield with a lower emissivity is better than two radiation shields with a higher emissivity to reduce the heat transfer rate. In addition, optimal combinations of radiation shields with different materials are proposed.

Shishkin, A. V.;Vostrikov, A. A.

2016 Thermophysics and Aeromechanics

doi: 10.1134/S0869864316040107

Abstract The oxidation of zirconium by supercritical CO2 (20.2 and 14.1 MPa, and 823 K) results in the formation of a layer of the monoclinic ZrO2 and amorphous carbon. In the constant electric field (E = 293 kV/m) when a sample of zirconium is the anode, the formation of amorphous carbon leads to a growth of leakage current and a twelve-fold increase in the rate of zirconium oxidation.