Large-eddy simulation of circular cylinder flow at subcritical Reynolds number: Turbulent wake and sound radiation

The flows past a circular cylinder at Reynolds number 3900 are simulated using large-eddy simulation(LES) and the far-field sound is calculated from the LES results. A low dissipation energy-conserving finite volume scheme is used to discretize the incompressible Navier–Stokes equations. The dynamic global coefficient version of the Vreman's subgrid scale(SGS) model is used to compute the sub-grid stresses. Curle's integral of Lighthill's acoustic analogy is used to extract the sound radiated from the cylinder. The profiles of mean velocity and turbulent fluctuations obtained are consistent with the previous experimental and computational results. The sound radiation at far field exhibits the characteristic of a dipole and directivity. The sound spectra display the-5/3 power law. It is shown that Vreman's SGS model in company with dynamic procedure is suitable for LES of turbulence generated noise.     

Investigation of the microstructure a turbulent jet in a wake

The study of the characteristics of the turbulence in the boundary layer and in free jets is one of the most important problems of the aerodynamics of viscous fluids. The accumulation of information on the pulsation characteristics of jet flows and the establishment of the corresponding governing laws may serve to verify the basic hypotheses of the semiempirical theories of turbulence, and also for the development of more advanced computational methods. In many cases the measurement of the pulsation characteristics of turbulent jets is of practical interest.The studies made up till now [1–5] of the microstructure of turbulent flow in the primary region of submerged axisymmetric jets have made it possible to obtain several interesting results. In particular, in addition to the average velocity profiles, hot-wire anemometric equipment has been used to measure the normal and tangential Reynolds stresses and also the intermittency factor in cross sections of the jet, the distribution of the intensity of the longitudinal and lateral velocity pulsations along the axis, the correlation coefficients and the corresponding integral turbulence scales, etc. These measurements have made it possible to draw several important conclusions on the mechanism of turbulent exchange, on the order of the terms omitted in the equation of motion, and on the semiempirical theories of turbulence [6–9].The common deficiency of the studies mentioned above is that near the boundary of a submerged jet, where the average velocity is practically equal to zero, the intensity of the pulsations is so great that it makes the reliability of the results obtained by means of the hotwire anemometer questionable. In this connection Townsend [6] indicated the advisability of studying the microstructure of a turbulent jet issuing into a low-velocity ambient flow.The present study had as its objective the investigation of the microstructure of the primary region of an axisymmetric jet in a wake flow over quite a broad range of the flow ratio parameter m=u/u₀;here u₀ is the average velocity at the nozzle exit, u is the velocity of the ambient stream. For various values of the parameter m in the primary region of the jet measurements were made of the profiles of the three components of the pulsation velocity and the Reynolds shear stresses, and also the values of the average velocity and two components of the pulsation velocity at a large number of points on the jet axis. The measured profiles of the Reynolds shear stresses were compared with the corresponding profiles calculated on the basis of the boundary layer equations from the experimentally determined average velocity profiles. For two values of the parameter m, in one of the sections of the jet measurements were made of the correlation coefficients of the longitudinal components of the pulsation velocity and the variation across the jet of the integral turbulence scale was determined.The results obtained give an idea of the influence of the parameter m on the characteristics of the turbulent jet in an ambient stream.     

Plate cascade in subsonic unsteady gas flow

Accounting for fluid compressibility creates serious difficulties in solving the problem of oscillations of a grid of thin, slightly curved profiles in a subsonic stream. The problem has been solved in [1–3] for a widely-spaced cascade without stagger whose profiles oscillate in phase opposition. The phenomenon of aerodynamic (acoustic) resonance, which may arise in a grid in the direction transverse to the stream for definite values of the stream velocity and profile oscillation frequency, was discovered in [2]. An approximate solution of the problem in which account is not taken of the effect of the vortex trails on the gas flow has been obtained in [4]. In [5, 6] Meister studied in the exact linear formulation the problem of unsteady gas motion through an unstaggered cascade of semi-infinite plates. In [7] Meister considered a grid of profiles with finite chords, but the problem solution was not completed. The problem of subsonic gas flow through a staggered lattice whose profiles oscillate following a single law with constant phase shift was solved most completely in the studies of Kurzin [8, 9] using the method of integral equations. A method of solving the problem for the case of arbitrary harmonic oscillation laws for the lattice profiles was indicated in [10]. The results of the calculation of the unsteady aerodynamic forces for the particular case of a plate cascade without stagger are presented in [9,11], and the possibility of the occurrence of aerodynamic resonance in the cascade in the directions transverse to and along the stream is indicated.Another method of solving the problem is given in [12], in which the more general problem of unsteady subsonic gas flow through a three-dimensional cascade of plates is solved. In the present study this method is applied to the solution of the problem of oscillations of staggered plate cascades in a two-dimensional subsonic gas flow. The results are presented of an electronic computer calculation of the unsteady aerodynamic characteristics of the cascade profiles, which show the essential influence of fluid compressibility on these characteristics. In particular, a sharp decrease of the aerodynamic damping in the acoustic resonance regimes is obtained.     

Influence of acoustics and the flow regime in the boundary layer on nozzle walls on the mixing layer of an immersed jet

Experimental determinations were made of the width of the mixing layer and the level of turbulent pulsations in the initial section of a subsonic circular immersed jet for different parameters of the boundary layer on the nozzle walls and in the presence of acoustic excitation. It was established that the rate of expansion of the turbulent mixing layer depends on the flow regime in the boundary layer. For laminar initial boundary layer, external acoustic excitation can lead to a decrease in the expansion velocity of the mixing layer and of the intensity of the velocity pulsations on the jet axis within the initial section. If the frequency and amplitude of acoustic excitation at which a decrease in the rate of expansion of the mixing layer and of the pulsation intensity was observed remained unchanged, the influence of the acoustics disappeared when the boundary layer became turbulent. The acoustic vibrations influenced the subsonic jets by generating vortex perturbations when they interact with the edge of the nozzle.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 36–42, November–December, 1982.We are grateful to K. I. Artamonov, now deceased, for support and discussing the results, and O. I. Navoznov and S. F. Agafonov for help in organizing and performing the experiments.     

基于比拟理论的翼型扰流声场数值模拟

从二维模型方程的全离散形式出发,重点分析了差分格式的色散特性和各向异性效应,证实迎风紧致格式比对称格式有更好的色散和各向同性特性,故有利于声场的数值模拟,并采用三阶迎风紧致格式（ＵＣＤ３）和四阶对称紧致格式（ＳＣＤ４）计算了绕ＮＡＣＡ００１２翼型的可压缩非定常流场,并将此流场作为近场声源,运用声学比拟理论对气动声进行模拟。     

Interaction of a turbulent boundary layer with a shock wave at transonic flow velocities

An asymptotic theory of the interaction of a turbulent boundary layer on a plate with a normal shock wave of low intensity has been constructed in various studies [1–4] under the assumption that the averaged velocity of the particles in the boundary layer in front of the interaction region satisfies a logarithmic law. In the present paper a different approach to this problem is proposed based on a power law of the velocity in the undisturbed boundary layer. The obtained results give different estimates for not only the sizes of the characteristic flow regions in the interaction region but also for the shock intensity leading to boundary layer separation.     

Estimate of the correlation between pressure pulsations and the divergence of the velocity in subsonic flows of variable density

In flows with variable density, the turbulence energy equation contains a large number of correlations, about which little is at present known [1]. One of the least studied is the correlation between the pressure and the divergence of the velocity. Usually, this correlation is ignored [2, 3]. The aim of the present paper is to estimate the pulsations of the divergence of the velocity and the correlation with the pressure pulsations in a subsonic turbulent flow with variable density. Three cases are considered: 1) mixing of gases having different densities, 2) diffusion combustion, 3) combustion of a homogeneous mixture. It is assumed that the Mach number is small, the Reynolds number large, and the coefficients of molecular diffusion and thermal diffusivity equal; external forces are absent.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 4–11, May–June, 1979.     

Investigation of downstream and sideline subsonic jet noise using Large Eddy Simulation

The sound fields radiated by Mach number 0.6 and 0.9, circular jets with Reynolds numbers varying from 1.7×10³ to 4×10⁵ are investigated using Large Eddy Simulations. As the Reynolds number decreases, the properties of the sound radiation do not change significantly in the downstream direction, whereas they are modified in the sideline direction. At low Reynolds numbers, for large angles downstream from the jet axis, the acoustic levels are indeed remarkably lower and a large high-frequency part of the sound spectra vanishes. For all Reynolds numbers, the downstream and the sideline sound spectra both appear to scale in frequency with the Strouhal number. However their peak amplitudes vary following two different velocity exponents according to the radiation direction. The present observations suggest the presence of two sound sources: a Reynolds number-dependent source, predominant for large radiation angles, connected to the randomly-developing turbulence, and a deterministic source, radiating downstream, related to a mechanism intrinsic to the jet geometry, which is still to be comprehensively described. This view agrees well with the experimental results displaying two distinguishable components in turbulent mixing noise [1, 2].     

Turbulent boundary layer on a catalytic surface in a nonequilibrium dissociating gas

The majority of the studies which consider the flow of a dissociating gas in a turbulent boundary layer are devoted to the investigation of either frozen or equilibrium flows on a flat plate.The frozen turbulent boundary layer has been studied by Dorrance [1], Kutateladze and Leont'ev [2], and Lapin and Sergeev [3]. A study of the effect of catalytic recombination processes at the plate surface on the heat transfer in a frozen turbulent boundary layer was made by Lapin [4].Kosterin and Koshmarov [5], Ginzburg [6], Dorrance [7], and Lapin [8] have studied the turbulent boundary layer on a plate in equilibrium dissociating gas.The calculation of the heat transfer in a turbulent boundary layer on a catalytic plate surface with nonequilibrium dissociation was made by Kulgein [9]. In this study the nonequilibrium nature of the dissociation process was taken into account only in the laminar sublayer, while the flow in the turbulent core was considered frozen. The solution was found numerically using a computer by means of a laborious iteration process.The present paper reports a method for calculating the turbulent boundary layer on a flat catalytic plate with arbitrary dissociation rate. The method, constructed using the assumptions customary for turbulent boundary layer theory, is a successive approximation method. Good convergence of the method is assured by the fact that the effect of the nonequilibrium nature of the dissociation process on the parameter distribution in the boundary layer and, consequently, on the friction and heat transfer may be allowed for merely by finding corrections, usually relatively small, to the distribution of these parameters in the equilibrium or frozen flows. The basis of the study is the two-layer scheme of the turbulent boundary layer. The Prandtl and Schmidt numbers and also their turbulent analogs are taken equal to unity. As the model of the dissociating gas we use the Lighthill model of the ideal dissociating gas [10], extended by Freeman [11] to nonequilibrium flows.     

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.     

Wall shear stress measurement in a turbulent pipe flow using ultrasound Doppler velocimetry

A turbulent boundary layer of a water flow is investigated by means of pulsed ultrasound Doppler velocimetry. The advantage of this method is the acquisition of complete velocity profiles along the sound propagation line within very short time intervals. The shear stress velocity, used for normalizing the velocity profiles, was determined by fitting the profiles to the universal profiles in a turbulent boundary layer obtained from Prandtl's mixing length theory. A coordinate transformation in the near-wall region is proposed to allocate the velocity data to "true" wall distances. From the experimental values of the wall shear stress velocity, the friction factors for a turbulent pipe flow are calculated and compared to the Blasius law. The overall error in measurement was estimated to NJ.4%.     

Compressibility effects due to turbulent fluctuations

This paper deals with intrinsic effects of compressibility, i.e. with dilatation fluctuations in response to pressure fluctuations. Three different types of turbulent flows are considered in more detail: homogeneous turbulent shear flow, wall-bounded turbulent shear flow and shock/turbulence interaction. A survey of the present knowledge in this field, mainly based on DNS data, is given. Using the linear inviscid perturbation equations a direct link between fluctuations of dilatation and of velocity in the direction of mean shear is presented for homogeneous shear flow. This relation might form the basis for a more universal pressure-dilatation model. It is conjectured that the insignificance of intrinsic compressibility effects in wall-bounded supersonic shear flow is mainly due to the impermeability constraint of the wall. To this end, a linear stability analysis of supersonic channel flow along cooled, but permeable walls has been performed based on Coleman et al.'s [5] mean flow data. It shows an increase in the moduli of eigenfunctions related to compressibility, like pressure, and in moduli of quantities derived from eigenfunctions such as ‘pressure dilatation’ and squared dilatation. Although these results do not prove our hypothesis they provide hints in this direction. Shock/turbulence interaction is viewed as a source of compressibility. Former DNS data of Hannappel and Friedrich [10] for shock/isotropic turbulence interaction showing the effect of compressibility on the amplification of fluctuations are interpreted based on linear perturbation equations.     

Pulsating motions of a solid particle in turbulent flow

The relations governing the transverse pulsations of a spherical particle in a turbulent flow are examined. On the basis of the resulting relations an estimate is made of the effect of several factors on the intensity of the velocity pulsations of the solid component in dispersed flow through a conduit.     

Turbulence intensity and turbulent transport characteristics behind gratings in pipes

The magnitudes of the variance ² of the concentration distribution C(y) of a material being mixed in a stream behind turbulizing gratings in a pipe are measured. The impurity (helium) was inserted in the air stream via a thin pipe mounted directly behind one of the grating orifices. The intensity of the turbulent stream velocity pulsations behind the grating was measured by the method of helium diffusion. Estimates of the magnitudes of the turbulent diffusion coefficient D and the diffusion scale of turbulencel were given on the basis of these measurements. Characteristics of the turbulent stream behind gratings in wind tunnels have been investigated in [1–3], for example. The process of turbulent diffusion behind gratings has been studied to a lesser degree. The flow directly behind the grating is of greatest interest for practice, but in the majority of papers the measurements have been made at a distance x from the grating which exceeds (10–20)M, where M is the distance, or step, between the centers of the grating orifices. Moreover, as a rule, the gratings investigated possessed low hydraulic drag and the intensity of stream turbulence ahead of them was considerably lower than ahead of the technical mixing apparatus.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 36–42, January–February, 1973.     

On the Lighthill relationship and sound generation from isotropic turbulence

In 1952 Lighthill (1952) developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman (1952) applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded-time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley (1994) revisited Proudman's work and applied the Lighthill relationship to evaluate the radiated acoustic power directly from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale-based Reynolds number 2.0×10³R 12.7×10³. Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 and 3.3 and the second-order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68–3.68.This research was supported by the National Aeronautics and Space Administration under NASA Contract No. NAS1-19480 while the first author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665, U.S.A.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Effect of particles suspended in a fluid flow on the decay of isotropic turbulence

Dynamic equations have been obtained for the two-point double correlations of the fluctuation velocities of a fluid and the particles suspended in it at low volume concentrations of the solid phase. In the case of uniform isotropic turbulence these equations can be considerably simplified. The final period of decay of isotropic turbulence has been studied in detail. At this stage in the case of high-inertia particles the inhomogeneous-fluid turbulence is similar to the turbulence of a homogeneous fluid (without particles) in the sense that the presence of the particles affects only the fluctuation energy but leaves unchanged the spatial scales of turbulence and the spatial energy spectrum function. The suspended particles lead to exponential damping of the turbulent pulsations.Little theoretical information is available on the hydrodynamics of a suspension of fine particles in a turbulent liquid or gas. Research has been mainly confined to the behavior of the individual particles in a given turbulence field [1]. The problem of the turbulent motion of the mixture as a whole has been examined by Barenblatt [2], who derived the equations of motion of the mixture, using Kolmogorov's hypothesis to close them. Hinze [3] has also attempted to derive equations for turbulent pulsations of the mixture. However, as Murray showed [4], Hinze' s equations contradict Newton' s third law.The effect of suspended particles on the turbulence of a two-phase flow is governed by the noncorrespondence of the local velocities of the particles and the medium. The forces of resistance to the motion of the particles relative to the fluid lead to additional dissipation of fluctuation energy and decay of turbulence [2]. On the other hand, if the averaged velocities of particles and medium do not correspond, the suspended particles may also have a destabilizing effect [5, 6], causing energy transfer from the averaged to the pulsating motion. Below we shall consider the case where the averaged velocities of the two phases coincide, i.e., we shall deal only with the first of the two above-mentioned effects.The authors thank G.I. Barenblatt for his useful advice.     

Boundary layer in the vicinity of the stagnation point of a body of axial symmetry in a turbulent stream

The flow in the boundary layer in the vicinity of the stagnation point of a flat plate is examined. The outer stream consists of turbulent flow of the jet type, directed normally to the plate. Assumptions concerning the connection between the pulsations in velocity and temperature in the boundary layer and the average parameters chosen on the basis of experimental data made it possible to obtain an isomorphic solution of the boundary layer equations. Equations are obtained for the friction and heat transfer at the wall in the region of gradient flow taking into account the effect of the turbulence of the impinging stream. It is shown that the friction at the wall is insensitive to the turbulence of the impinging stream, while the heat transfer is significantly increased with an increase in the pulsations of the outer flow. These properties are confirmed by the results of experimental studies [1–4].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 83–87, September–October, 1973.     

Noise Investigation of a High Subsonic, Moderate Reynolds Number Jet Using a Compressible Large Eddy Simulation

This study investigates the noise radiated by a subsonic circular jet with a Mach number of 0.9 and a Reynolds number of 65000 computed by a compressible Large Eddy Simulation (LES). First, it demonstrates the feasibility of using LES to predict accurately both the flow field and the sound radiation on a domain including the acoustic field. Mean flow parameters, turbulence intensities, velocity spectra and integral length scales are in very good agreement with experimental data. The noise generated by the jet, provided directly by the simulation, is also consistent with measurements in terms of sound pressure spectra, levels and directivity. The apparent location of the sound sources is at the end of the potential core in accordance with some experimental observations at similar Reynolds numbers and Mach numbers. Second, the noise generation mechanisms are discussed in an attempt to connect the flow field with the acoustic field. This study shows that for the simulated moderate Reynolds number jet, the predominant sound radiation in the downstream direction is associated with the breakdown of the shear layers in the central jet zone. Received 24 January 2002 and accepted 16 July 2002 Published online 3 December 2002 RID="*" ID="*" A preliminary version of some of the results presented here was reported in AIAA Paper 2000–2009 presented at the 6th AIAA/CEAS Aeroacoustics Conference in Lahaina, Hawaii, June 2000. Computing time was supplied by the Institut du Développement et des Ressources en Informatique Scientifique (IDRIS – CNRS). Communicated by T.B. Gatski     

Statistical characteristics of the diffusion of a chemically active additive in a turbulent mixing zone

In the article a numerical solution of the connected system of the equations of turbulent transfer for the fields of the velocity and concentration of a chemically active additive is used to calculate a number of the second moments of the concentration field in a flat mixing zone. The system of transfer equations is derived from the equations for a common function of the distribution of the fields of the pulsations of the velocity and the concentration [1] and is simplified in the approximation of the boundary layer. A closed form of the transfer equations is obtained on the level of three moments, using the hypothesis of four moments [2] and its generalized form for mixed moments of the field of the velocity and the field of a passive scalar. The differential operator of the closed system of the equations of turbulent transfer for the fields of the velocity and the concentration is found by a method of closure not of the parabolic type but of a weakly hyperbolic type [3]. An implicit difference scheme proposed in [4] is used for the numerical solution. The results of the numerical solution are compared with the experimental data of [5].