Wavy flow of a liquid film in the presence of a cocurrent turbulent gas flow

Wavy downflow of viscous liquid films in the presence of a cocurrent turbulent gas flow is analyzed theoretically. The parameters of two-dimensional steady-state traveling waves are calculated for wide ranges of liquid Reynolds number and gas flow velocity. The hydrodynamic characteristics of the liquid flow are computed using the full Navier-Stokes equations. The wavy interface is regarded as a small perturbation, and the equations for the gas are linearized in the vicinity of the main turbulent flow. Various optimal film flow regimes are obtained for the calculated nonlinear waves branching from the plane-parallel flow. It is shown that for high velocities of the cocurrent gas flow, the calculated wave characteristics correspond to those of ripple waves observed in experiments.     

Interaction of a turbulent stream of gas with a liquid film

We consider the turbulent motion of a gas in contact with a liquid film next to a wall. We assume that the stream of gas excites in the liquid a complex system of motions which are analogous in principle to the motions in the near-wall zone of a homogeneous turbulent stream with transverse shear. As a result of these motions, the gas stream has considerable turbulence even at the gas-film boundary. On this assumption, we calculate the relation between the pressure drop and the average gas velocity and find that it is in satisfactory qualitative and quantitative agreement with experimental results. As our scale of turbulence at the boundary, we took a linear variation as a function of the film thickness, which enabled us to describe the available experimental results satisfactorily, making use of two empirical constants.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 67–74, March–April, 1976.In conclusion, the authors take this opportunity to express their grateful recollection of conversations with the late Professor P. A. Semenov, who drew their attention to this problem, and also to thank G. G. Chernyi and G. A. Lyubimov for their comments and their interest in the work.     

Primary instabilities of liquid film flow sheared by turbulent gas stream

Evolution of excited waves on a viscous liquid film has been investigated experimentally for the annular gas–liquid flow in a vertical tube. For the first time the dispersion relations are obtained experimentally for linear waves on liquid film surface in the presence of turbulent gas flow. Both cocurrent and countercurrent flow regimes are investigated. As an example of comparison with theory, the experimental data are compared to the results of calculations based on the Benjamin quasi-laminar model for turbulent gas flow. The calculation results are found to be in good agreement with experiments for moderate values of film Reynolds number.     

Characteristics of solitary waves on a falling liquid film sheared by a turbulent counter-current gas flow

We report an experimental investigation of a falling water film sheared by a turbulent counter-current air flow in an inclined rectangular channel. Film thickness and wave velocity measurements associated with visual observation are conducted to study the influence of the air flow on controlled traveling waves consisting of a large wave hump preceded by capillary ripples. First, we focus on the variation of the shape, amplitude and velocity of the waves as the gas velocity is gradually increased. We demonstrate that the amplitude of the main hump grows substantially even for moderate gas velocities, whereas modification of the wave celerity becomes significant above a specific gas velocity around 4 m/s, associated with an alteration of the capillary region. The influence of the gas flow on 3D secondary instabilities of the solitary waves detected in a previous study Kofman et al. (2014), namely rugged or scallop waves, is also investigated. We show that the capillary mode is damped while the inertial mode is enhanced by the interfacial shear. Next, the gas velocity is increased until the onset of upstream-moving patterns referred to as flooding in our experiments. At moderate inclination angles (typically < 7 ^○), flooding occurs for a gas velocity around 8 m/s and is initiated at the scallop wave crests by a backward wave-breaking phenomenon preceded by the onset of ripples on the flat residual film separating two waves. At high inclination angle, a rapid development of solitons is observed as the air velocity is increased preventing the waves to turn back. Finally, at high liquid Reynolds number, sudden and intermittent events are triggered consisting of very large amplitude waves that go back upwards very fast. These “slugs” either extend over the whole width of the channel or are very localized and can thus potentially evolve towards atomization.     

Wave flow of a liquid film together with a flow of gas

The wave flow of a thin layer of viscous liquid in conjunction with a flow of gas was considered in a linear formulation earlier [1, 2]. In this paper the problem of the wave flow of a liquid film together with a gas flow is solved in a nonlinear setting. On this basis relationships are derived for calculating the parameters of the film and the hydrodynamic quantities.Ivanovo. Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 12–18, January–February, 1972.     

Evaporation of polydispersed droplets in a highly turbulent channel flow

A model experiment for the study of evaporating turbulent two-phase flows is presented here. The study focuses on a situation where pre-atomized and dispersed droplets vaporize and mix in a heated turbulent flow. The test bench consists in a channel flow with characteristics of homogeneous and isotropic turbulence where fluctuations levels reach very high values (25% in the established zone). An ultrasonic atomizer allows the injection of a mist of small droplets of acetone in the carrier flow. The large range diameters ensure that every kind of droplet behavior with regards to turbulence is possible. Instantaneous concentration fields of the vaporized phase are extracted from fluorescent images (PLIF) of the two phase flow. The evolution of the mixing of the acetone vapor is analyzed for two different liquid mass loadings. Despite the high turbulence levels, concentration fluctuations remain significant, indicating that air and acetone vapor are not fully mixed far from the injector.     

Unsteady turbulent flow of a gas suspension in a channel under conditions of injection and forced pressure oscillations

A turbulent flow of a suspension of solid particles in a gas is considered. The suspension is located in a channel with permeable walls (the pressure at the left end face of the channel follows a sinusoidal law). The flow considered here reflects the principal features of the flow in the combustion chamber of a solid-propellant rocket motor. The unsteady flow of the gas suspension is described by using the Eulerian-Lagrangian approach. A stochastic variant of the discrete-trajectory approach is used for modeling the particle motion. The influence of the condensed phase on the turbulence characteristics and acoustic oscillations of the parameters of the working medium in the channel in the case of injection is discussed. The calculated results are compared with data obtained in a physical experiment.     

Rarefied gas flow past a sphere with injection

Hypersonic rarefied gas flow over the windward face of a sphere is considered in the presence of distributed injection from the surface of the body. A similar problem was previously solved in [1–3] within the framework of continuum mechanics and in [4] on the basis of model kinetic equations. In the present study the calculations were carried out using the Monte Carlo method of direct statistical modeling [5, 6]. The injected gas was the same as the free-stream gas. A simple monatomic gas model with a rigid sphere interaction potential was employed. The reflection of the molecules from the surface of the body was assumed to be diffuse with total energy accommodation. The calculation procedure using weighting factors is described in [7]. The influence of injection on the mechanical and thermal effect of the gas flow on the body is investigated for various degrees of rarefaction of the medium and injection rates.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 175–179, July–August, 1990.     

Linear stability in the flow of a liquid film concurrent with a gas flow

The two-phase flow of liquid films are often encountered in practice, but the number of theoretical papers devoted to this problem is limited. The problem of the linear stability of a viscous liquid film subjected to a gas flow has been formulated in [1] and, in somewhat different form, in [2]. The linear stability of plane-parallel motion in films has been studied analytically in [1–8] for some limiting cases. The range of validity of the analytic approaches remains an open question. Therefore, an exact numerical analysis of flow stability over a fairly broad range is required. In the present paper a separate solution of the problem for the gas and the liquid is shown to be possible. The Orr-Sommerfeld equation has been integrated numerically, and the results are compared to the results of analytic calculations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 143–146, January–February, 1976.The author is grateful to É. É. Markovich for directing the work and to V. Ya. Shkadov for his interest in the work and many useful comments.     

Dynamics of a liquid film sheared by a co-flowing turbulent gas

Consider the dynamics of a thin laminar liquid film flowing over an inclined wall in the presence of a co-flowing turbulent gas. The solution to the full two-phase flow problem poses substantial technical difficulties. However, by making appropriate assumptions, the solution process can be simplified and can provide valuable insights. The assumptions allow us to solve the gas and liquid problems independently. Solving for the gas flow reduces to finding perturbations to pressure and tangential stresses at the interface, influencing the liquid problem through the boundary conditions. We analyze the effect of gas flow on the liquid problem by developing an integral-boundary-layer model, which is valid up to moderate liquid Reynolds numbers. We seek solitary-wave solutions of this model under the influence of gas flow via a pseudo-arclength continuation method. Our computations demonstrate that as a general trend, the wave speed increases with increasing the gas shear and the liquid flow rate. Further insight into the problem is provided via time-dependent computations of the integral-boundary-layer model.     

Evaporation of bi-component droplets in a heated, highly turbulent flow

This work aims to understand the phenomena that occur in a combustion chamber where multi-component fuel droplets are injected. Many evaporation models exist but the influence of turbulence on spray vaporization is not yet well understood. This study gives a useful database to improve these models. The objective of the work is to measure the dispersion and the evaporation of bi-component (octane/3-pentanone) droplets and the resulting vapor mixing in a well-known, heated, highly turbulent channel flow. The carrier flow shows high turbulence levels, flat profiles for the mean velocity and the velocity fluctuations. The injected droplets have a large variety of behaviors due to the large polydispersion and to the turbulence. The evolution of 3-pentanone liquid concentration, mass flux, and droplet clusters are described. Mean concentration, fluctuations of concentration, and mixing of the vapor phase are characterized.     

Nonlinear waves in a liquid film entrained by a turbulent gas stream

In the long-wavelength approximation and on the basis of a simplified system of equations analogous to the one considered by Shkadov and Nabil' [1, 2], an investigation is made into waves of finite amplitude in thin films of a viscous liquid on the walls of a channel in the presence of a turbulent gas stream. A bibliography on the linear stability of such plane-parallel flows can be found in [3–5]. The nonlinear stability is considered in [6]. A stationary periodic solution is sought in the form of a Fourier expansion whose coefficients are found near the upper curve of neutral stability by Newton's method and near the lower branch of the stability curve by the method of Petviashvili and Tsvelodub [7, 8].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No, 2, pp. 37–42, March–April, 1981.I thank V. Ya. Shkadov for supervising the work and all the participants of G. I. Petrov's seminar for a helpful discussion.     

Interfacial instability in turbulent flow over a liquid film in a channel

We revisit the stability of a deformable interface that separates a fully-developed turbulent gas flow from a thin layer of laminar liquid. Although this problem has received considerable attention previously, a model that requires no fitting parameters and that uses a base-state profile that has been validated against experiments is, as yet, unavailable. Furthermore, the significance of wave-induced perturbations in turbulent stresses remains unclear. To address these outstanding issues, we investigate this problem and introduce a turbulent base-state velocity that requires specification of a flow rate or a pressure drop only; no adjustable parameters are necessary. This base state is validated extensively against available experimental data as well as the results of direct numerical simulations. In addition, the effect of perturbations in the turbulent stress distributions is investigated, and demonstrated to be small for cases wherein the liquid layer is thin. The detailed modelling of the liquid layer also elicits two unstable modes, ‘interfacial’ and ‘internal’, with the former being the more dominant of the two. We show that it is possible for interfacial roughness to reduce the growth rate of the interfacial mode in relation to that of the internal one, promoting the latter, to the status of most dangerous mode. Additionally, we introduce an approximate measure to distinguish between ‘slow’ and ‘fast’ waves, the latter being the case for ‘critical-layer’-induced instabilities; we demonstrate that for the parameter ranges studied, the large majority of the waves are ‘slow’. Finally, comparisons of our linear stability predictions are made with experimental data in terms of critical parameters for onset of wave-formation, wave speeds and wavelengths; these yield agreement within the bounds of experimental error.     

Interaction of a gas-droplet turbulent jet with a cocurrent high-velocity high-temperature gas flow

A mathematical model and a method for calculating a gas-droplet turbulent jet with allowance for velocity nonequilibrium and virtual mass of the condensed phase during turbulent fluctuations and also heat and mass transfer within the three-temperature scheme are developed. Methodical calculations are performed. The results of these calculations are in reasonable agreement with available experimental data. The structure of the gas-droplet jet in a cocurrent high-velocity high-temperature gas flow is studied by numerical methods. The ratio of intensities of heat and mass transfer between the phases and turbulent diffusion transfers of substances is found to be different at the initial, transitional, and basic segments of the jet. This difference is responsible for the nonmonotonic axial distribution of vapor density and the lines of the halved mass flow of the condensed phase. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 85–94, May–June, 2008.     

Heat transfer in a wall jet propagating in a highly turbulent cocurrent flow

Heat transfer in a jet propagating in a cocurrent flow has been studied over wide ranges of the injection ratio (m=U_s/U₀<1 and m>1) and flow turbulence (Tu₀=0.2–25%). It is shown experimentally that for m<1, a 1% increase in turbulence leads to a 1% increase in heat transfer, and the wall adiabatic temperature and the relative heat-transfer function should be taken into account in heat-transfer calculations. For m>1, the flow turbulence does not affect the heat transfer and the heat production can be calculated using the laws typical of jet flows. Kutateladze Institute of Thermal Physics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 3, pp. 119–125, May–June, 1998.     

DNS of turbulent heat transfer in a channel flow with a high spatial resolution

Direct numerical simulations of turbulent heat transfer in a channel flow are performed to investigate the effects of Reynolds and Prandtl numbers on higher-order turbulence statistics such as a turbulent Prandtl number and the budget for the dissipation rate of the temperature variance. The Reynolds numbers based on the friction velocity and the channel half width are 180 and 395, and the molecular Prandtl numbers Pr’s 0.71–10.0. Careful attention is paid to ensure accuracy of the higher-order statistics through the use of a high spatial resolution comparable to Batchelor length scale. The wall-asymptotic value of the turbulent Prandtl number is mostly independent of Reynolds number for the current range of Pr’s. The budget for the dissipation rate of the temperature variance has been computed, and the negligible effect of a Reynolds number on the sum of all source and sink terms in near-wall region in the current computational range is found. This result is quite similar to the one in the budget for the dissipation rate of turbulent energy. In addition, a priori test for existing models is also performed to assess the Pr dependence on the individual terms and their summations in the budget.     

The concentration field in a turbulent channel flow with polymer injection at the wall

Instantaneous concentration profiles have been measured in turbulent water channel flows at 5 axial locations immediately downstream of a line, wall injection of a dyed 700 ppm polymer solution and for comparison, dyed water. Concentration was deduced from a line of fluoresced radiation that was stimulated by a laser beam directed through the dyed injectant and normal to the channel wall. Both statistical and time-resolved results show how the turbulent mixing is modified and damped when the injectant is a polymer solution. A version of this paper was presented at the 11th Symposium on Turbulence, University of Missouri-Rolla, Oct. 17–19, 1988