Theoretical investigation of coherent structures in stratified turbulent channel flow

A theoretical model is constructed of turbulent stratified flow in a flat horizontal channel with allowance for coherent structures that arise in it. The ordered part is separated from the turbulence of the flow and to describe the Reynolds-type equations are derived. The remaining part of the turbulence is taken into account parametrically in the form of an effective exchange coefficient. The flow is divided into a core, in which the ordered structures are manifested quite clearly, and wall regions, in which ordered large-scale structures are weakly manifested. To study the coherent structures in the core of the flow, an approach analogous to one already used to model ordered structures in open flows [4] is used. Monin-Obukhov scaling theory is used to describe the turbulence in the wall region.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 32–38, May–June, 1991.     

Numerical simulation of expansion/compression effect on shock induced turbulent flow

The interaction of homogeneous and isotropic turbulence with a shock wave is observed by solving the Reynolds-averaged Navier–Stokes equations with the k–? turbulence model. All turbulent fluctuations are measured at the period of expansion in the turbulent field and during compression by the reflected shock on turbulent field, and it is observed that the longitudinal turbulent velocity fluctuation is enhanced more at the period of expansion due to incident shock wave movement far from the turbulent field. The amplification of the turbulent kinetic energy (TKE) level in the shock/turbulence interaction depends on the shock wave strength and the longitudinal velocity difference across the shock wave. On decreasing the longitudinal velocity difference across the shock, the turbulent kinetic energy (TKE) level is less amplified. The TKE level is amplified by the factor of 1.5–1.8 in the shock/turbulence interaction where the dissipation rate of TKE decreases in all cases of shock/turbulence interaction. After the shock/turbulence interaction, the turbulent dissipative-length scale is amplified slightly and the amplification of the length scales decreases when increasing the shock strength. To cite this article: M.A. Jinnah, K. Takayama, C. R. Mecanique 333 (2005).     

Interaction of a normal shock wave with a compressible turbulent flow

A speckle photographic method, which is sensitive to changes of gradients in fluid density, is applied for analyzing a compressible turbulent air flow with density fluctuations. Spatial correlation coefficients, turbulent length scales, and energy spectra are determined under the assumption of homogeneous isotropic turbulence. The experiments are performed in a shock tube where the flow is passed through a turbulence grid. Measurements are taken before and after the turbulent regime interacts with the normal shock wave reflected from the tube's end wall. Amplification of the turbulence intensity by the shock interaction process is verified quantitatively and is shown to be restricted to the lower wave numbers in the spectrum.A version of this paper was presented at the 11th Symposium on Turbulence, University of Missouri-Rolla, Oct. 17–19, 1988.To Professor Dr.-Ing. Klaus Gersten on the occasion of his 60th birthday     

Wave turbulence and intermittency in directional wave fields

The evolution of surface gravity waves is driven by nonlinear interactions that trigger an energy cascade similarly to the one observed in hydrodynamic turbulence. This process, known as wave turbulence, has been found to display anomalous scaling with deviation from classical turbulent predictions due to the emergence of coherent and intermittent structures on the water surface. In the ocean, waves are spread over a wide range of directions, with a consequent attenuation of the nonlinear properties. A laboratory experiment in a large wave facility is presented to discuss the sensitivity of wave turbulence on the directional properties of model wave spectra. Results show that the occurrence of coherent and intermittent structures become less likely with the broadening of the wave directional spreading. There is no evidence, however, that intermittency completely vanishes.     

Turbulent air flow over a moving curved surface

A numerical model of turbulent air flow over a curved surface is described. The model is based on two-dimensional nonlinear Reynolds equations and continuity equations written in a coordinate system moving with the profile of the curved surface. The Reynolds stresses are represented in the form of the product of the isotropic turbulent viscosity coefficient, which increases linearly with height, and the deformation tensor of the mean velocity field. Flow over a stationary sinusoidal surface and a sinusoidal gravity wave on water is simulated. The structure of the velocity and pressure wave fields is obtained. The differences in flow over stationary and moving surfaces are analyzed.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 20–24, September–October, 1986.     

Averaged equations and wave jumps describing the propagation of stokes waves with slowly varying parameters

The equations describing the stationary envelope of periodic waves on the surface of a liquid of constant or variable depth are investigated. Methods previously used for investigating the propagation of solitons [1–5] are extended to the case of periodic waves. The equations considered are derived from the cubic Schrödinger equation assuming slow variation of the wave parameters. In using these equations it is sometimes necessary to introduce wave jumps. By analogy with the soliton case a wave jump theory in accordance with which the jumps are interpreted as three-wave resonant interactions is considered. The problems of Mach reflection from a vertical wall and the decay of an arbitrary wave jump are solved. In order to provide a basis for the theory solutions describing the interaction of two waves over a horizontal bottom are investigated. The averaging method [6] is used to derive systems of equations describing the propagation of one or two interacting wave's on the surface of a liquid of constant or variable depth. These systems have steady-state solutions and can be written in divergence form.The author wishes to thank A. G. Kulikovskii and A. A. Barmin for useful discussions.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 113–121, September–October, 1989.     

Comparison of the RANS and PDF methods for air-particle flows

Two simulation methods, namely Reynolds-Averaged Navier–Stokes (RANS) equations, and Probability Distribution Function (PDF) are currently widely used for the modeling of multiphase flows. These two approaches are supplemented with appropriate closure equations that take into account all the pertinent forces and interaction effects on the solid particles, such as: particle–turbulence interactions; turbulence modulation; particle–particle interactions; particle–wall interactions; gravitation, drag and lift forces. The two methods have been used in order to simulate the turbulent particulate flow in upward pipes. The flow domain in all cases was a cylindrical pipe and the computations were carried for upward pipe flow. Monodisperse as well as polydisperse mixtures of particles have been considered. In general, the average velocity results obtained from the two methods are in close agreement, because the methods predict well the average velocity distribution of the carrier fluid as well as the solids. Thus, the differences in the average axial velocities predicted by the methods are not substantial. Differences in the turbulence intensity are more significant. A comparison of the numerical results obtained shows the relative importance of retaining the diffusion terms in both the axial and radial directions in the RANS method. Also the comparisons of the results show the relative effect of the lift forces in the distribution of solid particles.     

Interaction of wall turbulence with a compliant surface

The effect of various parameters of the compliant surface on the interaction with wall turbulence is analyzed using the monoharmonic approximation. It is shown that the interaction is resonant in character and that for certain values of the parameters a considerable reduction in turbulent skin friction is possible.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 67–72, July–August, 1990.     

Two-scale semiempirical theory of turbulent boundary layers and jets

To characterize the turbulence of boundary layers in the energy-bearing interval of wave numbers several turbulence scales are sometimes used (for example, [1, 2]). In particular, the universality of the semiempirical model of turbulence [2] can be extended in this way. A turbulence model with one equation (energy balance of the turbulence) has been constructed and used [3–6] and it has been established that the number of problems that can be solved for a universal choice of the values of the empirical coefficients increases appreciably if not one but two turbulent scales are used. In the present paper, it is shown that the introduction of a second scale makes it possible to take into account the interaction of shear layers in flows with two shear layers (for example, a channel or jet), and also to take into account the influence of turbulence of an external flow on a boundary layer. The interaction of shear layers is taken into account in theories containing a transport equation for the turbulent frictional stress _t (for example, [7]), in which the essence of the interaction reduces to diffusion of _t from layer to layer. In the present paper, a predominant volume interaction effect is assumed. It takes the form of a difference between the interaction of large-scale vortices with a shear deformation motion in flows with one and two shear layers, and also in the presence of turbulence in an external flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 8, pp. 17–25, November–December, 1982.     

Energy-balance equation for turbulent pulsations in the theory of the free turbulent boundary layer

Semiempirical expressions are proposed for the coefficient of turbulent viscosity and for the scale of turbulence in the equations for the free turbulent boundary layer in an incompressible fluid, these equations consisting of the equation of continuity, the equations of motion, and the equation for the average energy balance in the turbulent pulsations. The advantage of the expressions over the existing ones is that the two empirical constants in the equations have nearly the same values for circular and plane turbulent streams and also for a turbulent boundary layer at the edge of a semiinfinite homogeneous flow with a stationary fluid. The mean-energy distribution and the mean energy of the turbulent pulsations computed in this paper agree well with the experimental values.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 75–79, November–December, 1970.     

Diffraction of a plane shock wave by a disk that enters a compressible liquid at an inclination

A study is made of the three-dimensional problem of the interaction between a disk that enters water asymmetrically and a shock wave that is moving toward the disk. The water is assumed to be a perfect compressible liquid and the flow adiabatic. The changes in the flow parameters and the state are determined by numerical integration of the equations that describe the flow. A three-dimensional version of the finite-difference scheme of [1] is used in accordance with the method of [2]. The influence of the intensity of the shock wave on the drag coefficient of the disk and the shape of the free surface is investigated.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 185–187, January–February, 1990.     

Korteweg-de vries-type equations in the theory of surface waves in electrically conducting liquids

Equations are obtained which describe the propagation of long waves of small, but finite amplitude in an ideal weakly conducting liquid and on the basis of these equations the influence of MHD interaction effects on the characteristics of the solitary waves is investigated. The wave equations are derived under less rigorous constraints on the external magnetic field and the MHD interaction parameter than in [1–3]. It is shown that the evolution of the free surface is described by the KdV-Burgers or KdV equations with a dissipative perturbation, and that the propagation velocity of the solitary waves depends on the strength of the external magnetic field.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 177–180, November–December, 1989.     

On the Interaction of Turbulent Shear Layers with Harmonic Perturbations

The problem of coherent perturbations in a turbulent shear layer is considered for the purpose of developing a mathematical model based on a triple decomposition that extracts the coherent components of random fluctuations. The governing equations for the mean and the coherent parts of flow are derived, assuming the eddy-viscosity equivalence for the random part of flow, and solved by iterations to provide a coupled solution of the problem as a whole. Calculations agree well with experimental data in the upstream part of the layer where the mean–coherent flow interaction is the most important. In this region, the interaction changes the mean flow velocity distribution in such a manner that the neutral stability curve is shifted upstream relative to its position in the undisturbed layer and the perturbation intensity decreases further downstream. Experiments show that the coherent waves suppress the turbulent Reynolds stress production downstream of this region, but the model fails to predict the layer spreading correctly probably due to an inadequate turbulence closure of the mean flow. For the case of a turbulent mixing layer, we suggest a new closure relation that takes into account this coherent-random interaction.     

Theory of weak turbulence of waves on a fluid surface

The nonlinear interaction of waves in a fluid of finite depth is discussed. Forbidden decay processes in the gravitational portion of the spectrum are eliminated from the Hamiltonian by means of a canonical transformation. This provides an opportunity to obtain a kinetic equation which takes into account scattering of capillary waves by gravitational waves, in addition to decays in the subsystem of gravitational waves. The distribution N_k P^1/2h^1/4k^–4 is obtained for capillary waves in shallow water with constant flow of energy P with respect to the spectrum in the space of the wave numbers k. The interaction of the gravitational and capillary turbulence spectra is discussed. An induced distribution of gravitational waves is found which results from their interaction with capillary waves. It is an increasing function of the wave numbers q in the region bounded by the capillary constant k_o, N_q q^9/4 (q < k_o). The coupling of spectra in the gravitational and capillary regions and the conversion from slightly turbulent distributions to universal distributions are discussed.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 6, pp. 97–106, November–December, 1974.     

Localized Coherent Structures in the Boundary Layer

A Blasius laminar boundary layer and a steady turbulent boundary layer on a flat plate in an incompressible fluid are considered. The spectral characteristics of the Tollmien—Schlichting (TS) and Squire waves are numerically determined in a wide range of Reynolds numbers. Based on the spectral characteristics, relations determining the three–wave resonance of TS waves are studied. It is shown that the three–wave resonance is responsible for the appearance of a continuous low–frequency spectrum in the laminar region of the boundary layer. The spectral characteristics allow one to obtain quantities that enter the equations of dynamics of localized perturbations. By analogy with the laminar boundary layer, the three–wave resonance of TS waves in a turbulent boundary layer is considered.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A Robust High-Resolution Split-Type Compact FD Scheme for Spatial Direct Numerical Simulation of Boundary-Layer Transition

In this paper an efficient split-type Finite-Difference (FD) scheme with high modal resolution – most important for the streamwise convection terms that cause wave transport and interaction – is derived for a mixed Fourier-spectral/FD method that is designed for the spatial direct numerical simulation (DNS) of boundary-layer transition and turbulence. Using a relatively simple but thorough and instructive modal analysis we discuss some principal trouble sources of the related FD discretization. The new scheme is based on a 6th-order compact FD discretization in streamwise and wall-normal direction and the classical 4th-order Runge–Kutta time-integration scheme with symmetrical final corrector step. Exemplary results of a fundamental-(K-) type breakdown simulation of a strongly decelerated Falkner–Skan boundary layer (Hartree parameter _H = – 0.18) using 70 mega grid points in space are presented up to the early turbulent regime (Re_,turb 820). The adverse pressure gradient gives rise to local separation zones during the breakdown stage and intensifies final breakdown by strong amplification of (background) disturbances thus enabling rapid transition at moderate Reynolds number. The appearance and dynamics of small-scale vortical structures in early turbulence basically similar to the large-scale structures at transition can be observed corroborating Kachanov's hypothesis on the importance of the K-regime of breakdown for coherent structures in turbulence.     

Thermo-mechanical modeling of turbulent heat transfer in gas–solid flows including particle collisions

A thermo-mechanical turbulence model is developed and used for predicting heat transfer in a gas–solid flow through a vertical pipe with constant wall heat flux. The new four-way interaction model makes use of the thermal k_θ–τ_θ equations, in addition to the hydrodynamic k–τ transport, and accounts for the particle–particle and particle–wall collisions through a Eulerian/Lagrangian formulation. The simulation results indicate that the level of thermal turbulence intensity and the heat transfer are strongly affected by the particle collisions. Inter-particle collisions attenuate the thermal turbulence intensity near the wall but somewhat amplify the temperature fluctuations in the pipe core region. The hydrodynamic-to-thermal times-scale ratio and the turbulent Prandtl number in the region near the wall increase due to the inter-particle collisions. The results also show that the use of a constant or the single-phase gas turbulent Prandtl number produces error in the thermal eddy diffusivity and thermal turbulent intensity fields. Simulation results also indicate that the inter-particle contact heat conduction during collision has no significant effect in the range of Reynolds number and particle diameter studied.     

Simulation of ordered structures in open turbulent flows

Extensive experimental material [1–4] indicates that ordered (coherent) structures play an important part in determining the nature of the flow, the generation of Reynolds stresses and turbulence energy, and the transport of heat, momentum, and passive admixtures in a turbulent flow. In the present paper, a model is constructed for describing coherent structures in which, given the profile of the mean velocity, one can determine the characteristic sizes, the propagation velocities, and also the frequency and amplitude characteristics of these ordered motions. The model is based on the analogy between the ordered formations and secondary flows in a subsidiary laminar flow whose velocity profile is the same as the turbulent profile of the mean velocity. The influence of small-scale pulsations is described by the introduction of the coefficient of turbulent viscosity. In the framework of the model, numerical calculations are made for two-dimensional turbulent flows in a mixing layer, a jet, and a wake behind a cylinder. The results of the calculations are compared with experimental data.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 45–52, July–August, 1981.     

Calculation of principal characteristics of turbulent streams in a state of structural equilibrium

A review of articles on the study of turbulent streams having transverse displacement, in which a turbulent energy balance equation is used, is contained in [1]. Levin [2] proposed a certain development of Rotta's method [3] making it possible to determine the characteristics of the average flow and the radial distribution of pulsation magnitudes. However, in this article the scale of the turbulence (the quantityl) was given as an empirical function of the coordinates. At the same time it is clear that the distribution of the turbulence scale depends on the conditions of the problem. A special differential equation proposed in [4,5] describing the variation in time and space of the quantityl has the drawback that in deriving this equation it is necessary to invoke additional hypotheses which are difficult to test experimentally. In the present article, along with the velocity of the average flow, the pressure, and the pulsation magnitudes, the scale of the turbulence is considered as an important characteristic of the stream, determined by the reference system which consists of the Reynolds equations, continuity equations, and equations for the component of the Reynolds stress tensor. Rotta's approximate semiempirical relations and an experimental relation for the single-point correlation coefficient between the turbulent pulsations in velocity are used for closure of the system obtained. An approximate calculation is given for the principal average and pulsation characteristics of the flow for the region of the stream where the turbulence is in a state of structural equilibrium [6]. A comparison of the calculated and experimental data is presented.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 95–99, January–February, 1973.